PMCCPMCCPMCC

Search tips
Search criteria 

Advanced


 
Formats
Article sections
Authors
Related links
Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
 
PLoS One. 2011; 6(11): e27380.
Published online 2011 November 8. doi:  10.1371/journal.pone.0027380
PMCID: PMC3210806
Copyright Neiman et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Identifying Temporal Codes in Spontaneously Active Sensory Neurons
Alexander B. Neiman,1,2 David F. Russell,1,3 and Michael H. Rowe1,3*
1Neuroscience Program, Ohio University, Athens, Ohio, United States of America
2Department of Physics and Astronomy, Ohio University, Athens, Ohio, United States of America
3Department of Biological Sciences, Ohio University, Athens, Ohio, United States of America
Michal Zochowski, Editor
University of Michigan, United States of America
* E-mail: rowe/at/ohio.edu
Conceived and designed the experiments: MHR DFR ABN. Performed the experiments: MHR DFR. Analyzed the data: MHR DFR ABN. Contributed reagents/materials/analysis tools: MHR ABN DFR. Wrote the paper: ABN MHR DFR.
Received July 29, 2011; Accepted October 15, 2011.
Abstract
The manner in which information is encoded in neural signals is a major issue in Neuroscience. A common distinction is between rate codes, where information in neural responses is encoded as the number of spikes within a specified time frame (encoding window), and temporal codes, where the position of spikes within the encoding window carries some or all of the information about the stimulus. One test for the existence of a temporal code in neural responses is to add artificial time jitter to each spike in the response, and then assess whether or not information in the response has been degraded. If so, temporal encoding might be inferred, on the assumption that the jitter is small enough to alter the position, but not the number, of spikes within the encoding window. Here, the effects of artificial jitter on various spike train and information metrics were derived analytically, and this theory was validated using data from afferent neurons of the turtle vestibular and paddlefish electrosensory systems, and from model neurons. We demonstrate that the jitter procedure will degrade information content even when coding is known to be entirely by rate. For this and additional reasons, we conclude that the jitter procedure by itself is not sufficient to establish the presence of a temporal code.
A fundamental question in sensory neuroscience is how information is encoded in spike trains. The question often takes the form of distinguishing between rate codes, in which information is encoded in terms of the number of spikes within an encoding window, and temporal codes, in which the position of spikes within an encoding window carries information beyond that available from the number of spikes in the window [1]. Temporal codes are usually associated with nonlinear relations between the Fourier components of a stimulus and a neuronal response [1], [2], i.e. correlations between a particular frequency component of a stimulus and higher-frequency components of the response. These nonlinear relations provide information about the stimulus beyond that provided by linear correlations within the frequency band of the stimulus. In contrast, rate coding can be nonlinear, but it is characterized by a lack of correlation between Fourier components of the stimulus and higher-frequency components of the response, or by the fact that such nonlinear correlations, when present, do not provide any additional information about the stimulus. The pioneering work of Adrian [3] provided clear evidence that cutaneous sensory afferents use firing rate to encode stimulus intensity (a concise history of this work and related issues is in [4]). More recent work on a number of sensory systems has provided equally compelling evidence that precise spike timing can carry information beyond that available from measures of firing rate (e.g., [5]–[17] among many others).
An additional consideration is that primary afferent neurons in a variety of sensory systems exhibit an ongoing background discharge. Examples include vestibular afferents [18], [19], and electroreceptor afferents in several aquatic species [20]–[22]. Such background firing can arise from a variety of mechanisms including intrinsic oscillators, intrinsic noise, or random synaptic events. The resulting discharges span the spectrum from highly periodic to completely random spike sequences. Several studies have attempted to relate the properties of this background discharge to the stimulus encoding properties of afferents, by stimulating a system with time-varying Gaussian noise, and assessing information transmission based on various information metrics calculated from their responses (reviewed in [4], [10], [23]).
To assess the relative importance of firing rate versus precise spike timing in stimulus encoding, a computational procedure is often used in which the time of each spike is “jittered” by the addition of a variable time offset, chosen randomly from a zero-mean distribution [6], [20], [24]–[26]. The jittering produces a surrogate data set for which information metrics can be computed and compared to the same metrics computed from the original data. If the addition of jitter significantly decreases the information transmission and/or encoding efficiency of the afferent, as happens, for example, for some vestibular afferents [24], then the existence of a temporal encoding scheme is inferred.
However, the distinction between a rate code and a timing code can be problematic for a number of reasons. First, as discussed by Theunissen and Miller [1], the use of spike timing to encode transient or high frequency components of a stimulus can be consistent with a rate coding scheme, e.g. [6], [27]. Nor does the use of a temporal encoding scheme require high spike timing precision. Even in the case of a highly periodic spontaneously firing neuron, which like all self-sustained oscillators is inherently nonlinear, the response magnitude at different points in the neuron's cycle (its phase response curve) can be closely related to its linear response function [28], [29]. Weak stimuli can be linearly encoded in the instantaneous firing rate of a periodically firing neuron, and this encoding can be accounted for within the framework of linear response theory [28], [29]. Thus, the intrinsic timing precision of a periodically firing neuron is not necessarily indicative of a temporal code as understood in the current neuroscience literature.
Second, the linear stimulus reconstruction technique [1], [23], [30] that is typically used in conjunction with the jitter procedure treats a neuron as a linear “black box” whose transfer function is tuned to minimize the mean square error of stimulus estimation. This technique essentially assumes a rate code, since the stimulus is estimated by convolving a spike train with the response function of the optimal linear filter. Adding external noise in the form of jitter is equivalent to a distortion of the transfer function of the optimal filter. Thus, conclusions about the existence of a nonlinear time code drawn solely from application of a linear stimulus reconstruction technique may be questionable.
Third, the rationale for jitter analysis is based on the assumption that the standard deviation (SD) of the jitter distribution is small relative to the duration of an “encoding window”, so that the number of spikes within the window is unaffected, and only their temporal position within the window is altered. Thus, the SD of the jitter is normally chosen to be much smaller than the characteristic time scale of the stimulus on the assumption that this will be less than the duration of the encoding window. However, since the duration of the encoding window itself is never determined, this assumption cannot be validated, and so the results of artificial jittering should be interpreted with caution.
Here, we develop an analytical framework that provides a detailed, quantitative assessment of the effects of artificial jitter on spike train metrics commonly used to analyze sensory encoding: coefficient of variation, serial correlations, power spectral density, transfer functions, and coherence functions. This theoretical analysis allows us to specify precisely the relationships between these metrics as calculated for original and jittered spike trains. Using this framework, we show that jitter alters the higher order statistics of spike trains by introducing spurious serial correlations among interspike intervals. This can alter encoding properties. More importantly, we show that for weak stimuli and linear responses, jitter merely increases the noise in the background discharge. This occurs independently of any applied lower-frequency stimulus, and with minimal effects on stimulus-response gain. The additional noise from jitter results in suppression of the stimulus-response coherence (or the linear reconstruction kernel), and consequently of the mutual information rate, as estimated with the linear reconstruction technique. We illustrate these theoretical results by applying them to a model neuron with gamma-distributed interspike intervals, to a phase model of a periodically firing neuron, and to experimental data from vestibular and electroreceptor afferents. Although we focus on spontaneously active neurons, the theory, results, and conclusions we develop have broad applicability to analyses of sensory encoding.
Experimental
Turtle
The activity of vestibular posterior canal afferents was recorded in in vitro preparations of red-eared turtles, Trachemys (Pseudemys) scripta elegans, of 10–13 cm carapace length, as in [31]. Turtles were sacrificed by decapitation, and further dissection was done in a bath of oxygenated turtle Ringer's solution. After removal of the dorsal cranium, the brainstem was transected at the meso-thalamic junction, and the rostral portions discarded. A small hole was drilled in the bone overlying the posterior canal, approximately 2–3 mm from the posterior ampulla, to allow placement of a mechanical probe on the posterior semicircular duct. Stimuli consisted of indentations of the posterior semicircular duct using this probe. The head was then placed in a humidified recording chamber that was continuously infused with mixture of 95% O2/5% CO2. All procedures were approved by the Ohio University Institutional Animal Care and Use Committee (IACUC) (protocol number L01-35). Afferent spikes were recorded with glass micropipettes filled with 2 M NaCl and having electrical impedances of 50–100 MΩ. The electrodes were inserted into the posterior division of the VIIIth nerve along the antero-dorsal margin [32]. Signals from the electrodes were amplified, digitized at 10 kHz, and stored for offline analysis using Spike2™ software (Cambridge Electronic Design, CED).
Paddlefish
The spontaneous firing of electroreceptor afferents of paddlefish (Polyodon spathula) was recorded in in vivo preparations, as in [21]. This passive electrosensory system has cutaneous electroreceptors of the cathodally excited “Lorenzinian” ampullary type, like sharks and rays. Paddlefish are named for their “rostrum”, a flattened sensory appendage shaped like a canoe paddle, covered with electroreceptors, projecting forward from the head. A fish was held in an all-plastic chamber, maintained by a stream of oxygenated 22°C water into the mouth. A special advantage of paddlefish is that water flow or turbulence around electroreceptors on the rostrum, which might disturb spontaneous afferent firing, could be stilled by partitioning the chamber bath using a slab of agarose across the base of the rostrum. The cranium was opened dorsally to expose the sensory ganglion of the anterior lateral line cranial nerve, on one or both sides. A tungsten microelectrode was advanced into this ganglion to record single-unit spikes from an afferent's cell body. If used, electrical stimuli were applied from a 2 mm dipole electrode, or between two large chlorided Ag plates at the chamber ends, connected to a low-noise linear constant-current electrical isolator, commanded by a CED computer interface replaying arbitrary waveforms such as band-limited Gaussian noise. Data acquisition was similar to that for turtle. Data were from experiments at University of Missouri-St. Louis, under an IACUC-approved animal use protocol (W01-13) there.
The duration of Gaussian stimuli was 300–500 s for the turtle posterior canal afferents and 180 s for the paddlefish electroreceptor afferents.
Data analysis
Data analyses were performed offline using custom software programmed in MATLAB. The same analyses were used for both experimental recordings and numerical simulations. Definitional equations for the analyses are included for clarity.
Spontaneous discharge statistics
Three metrics were used to characterize spontaneous discharge: the coefficient of variation of the interspike interval distribution, serial correlation coefficients in ISI sequences, and the power spectral density of a spike train. Given a sequence of spike times t1,t2,...,tK, the corresponding sequence of interspike intervals (ISIs) is An external file that holds a picture, illustration, etc. Object name is pone.0027380.e001.jpg. The variability of an ISI distribution was characterized using the coefficient of variation, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e002.jpg, where An external file that holds a picture, illustration, etc. Object name is pone.0027380.e003.jpg is the mean ISI, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e004.jpg is the SD of the ISI distribution, and An external file that holds a picture, illustration, etc. Object name is pone.0027380.e005.jpg denotes averaging over k intervals. The mean firing rate is An external file that holds a picture, illustration, etc. Object name is pone.0027380.e006.jpg.
Serial correlation coefficients (SCCs) are derived from the normalized ISI autocorrelation function, and estimate the average degree to which ISIs are correlated or anticorrelated with other ISIs in the sequence. SCC values were calculated as An external file that holds a picture, illustration, etc. Object name is pone.0027380.e007.jpg, where m denotes the number of intervening intervals (lags), and ranged from 0 to 100 [33]. SCCs can range from −1 (perfect anticorrelation) to +1 (perfect correlation), while a value of 0 signifies no correlation. Spike generation is referred to as a renewal process if all SCCs in the spike sequence for m≥1 are 0, i.e., the ISIs are statistically independent.
Power spectral density (PSD), An external file that holds a picture, illustration, etc. Object name is pone.0027380.e008.jpg, is a measure of the distribution of a signal's energy in the frequency domain. It is particularly useful for identifying periodicities in a signal, expressed by peaks at particular frequencies. For purposes of PSD calculations, each neuronal spike train was represented as a sequence of Dirac delta functions centered at spike times An external file that holds a picture, illustration, etc. Object name is pone.0027380.e009.jpg from which the mean firing rate An external file that holds a picture, illustration, etc. Object name is pone.0027380.e010.jpg has been subtracted, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e011.jpg. The PSD of a spike train, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e012.jpg, has units of (spikes/s)2/Hz or simply Hz. The delta functions were approximated by rectangular pulses of height An external file that holds a picture, illustration, etc. Object name is pone.0027380.e013.jpg, where An external file that holds a picture, illustration, etc. Object name is pone.0027380.e014.jpg is the sampling interval, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e015.jpgs. The PSD was then estimated using the Welch method (function pwelch in the MATLAB Statistical Toolbox) with a 2.048 s Hamming window.
Information measures
We used two approaches to assess information encoding in neural responses to external stimuli. The first was a conventional linear reconstruction technique that estimates the lower bound of the mutual information rate, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e016.jpg [23], [34]. In this approach, a Gaussian stimulus, s(t), is applied to a neuron, and an estimate of the stimulus, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e017.jpg, is obtained from the neural response by convolving the output spike train with an optimal linear filter that minimizes the SD of the noise in the reconstruction, calculated as An external file that holds a picture, illustration, etc. Object name is pone.0027380.e018.jpg. The characteristics of the optimal filter are specified by its transfer function, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e019.jpg, where An external file that holds a picture, illustration, etc. Object name is pone.0027380.e020.jpg is the cross-spectral density of the stimulus and response, and An external file that holds a picture, illustration, etc. Object name is pone.0027380.e021.jpg is the PSD of the stimulated spike train (response). The lower bound of the mutual information rate is estimated from the SR coherence function [35] as:
A mathematical equation, expression, or formula. Object name is pone.0027380.e022.jpg
(1)
where fc is the stimulus cutoff frequency. SR coherence is a normalized measure of stimulus-response cross-correlation at different frequencies, defined as:
A mathematical equation, expression, or formula. Object name is pone.0027380.e023.jpg
(2)
where An external file that holds a picture, illustration, etc. Object name is pone.0027380.e024.jpg is the PSD of the stimulus. Calculations of SR coherence were done using the MATLAB function mscohere, with windowing as for PSDs.
The quality of the reconstruction is quantified by the coding fraction, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e025.jpg, defined [30] as:
A mathematical equation, expression, or formula. Object name is pone.0027380.e026.jpg
(3)
where A is the SD of the stimulus, and An external file that holds a picture, illustration, etc. Object name is pone.0027380.e027.jpg is the SD of the reconstruction noise, which can be calculated from the SR coherence [30] as: An external file that holds a picture, illustration, etc. Object name is pone.0027380.e028.jpg. The coding fraction ranges from 0 (encoding on a chance level) to 1 (perfect encoding).
The second information metric we used was response-response (RR) coherence, which provides an estimate of the upper limit of the mutual information rate [26], [36], [37]. In this method, a neuron is stimulated by a sequence of identical segments of a Gaussian noise stimulus. Each stimulus segment results in a response An external file that holds a picture, illustration, etc. Object name is pone.0027380.e029.jpg. The average coherence between responses is:
A mathematical equation, expression, or formula. Object name is pone.0027380.e030.jpg
(4)
where An external file that holds a picture, illustration, etc. Object name is pone.0027380.e031.jpg is the cross-spectral density of m-th and n-th responses, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e032.jpg is the PSD of m-th response, and angled brackets indicate averaging over the ensemble of responses. The square root of RR coherence serves as the upper bound of SR coherence, so that the following inequality holds: An external file that holds a picture, illustration, etc. Object name is pone.0027380.e033.jpg.
Numerical jitter analysis
A jittered response was obtained by adding independent zero-mean Gaussian random time offsets An external file that holds a picture, illustration, etc. Object name is pone.0027380.e034.jpg to each spike time. After computing the SR coherence for the jittered spike train, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e035.jpg, we obtained the lower bound of the mutual information rate, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e036.jpg, and the coding fraction, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e037.jpg, and compared them to An external file that holds a picture, illustration, etc. Object name is pone.0027380.e038.jpg and An external file that holds a picture, illustration, etc. Object name is pone.0027380.e039.jpgof the original spike train. This was repeated for various values of jitter SD, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e040.jpg, its magnitude. The tilde symbol denotes measures calculated from the jittered spike train.
Theory for jittered spike trains
The derivations of equations used in this section to express the exact relationships between statistical metrics of original and jittered spike trains are given in Appendix S1. In the analysis that follows, we assumed that both the stimulus and response are stationary stochastic processes, and we used zero-mean Gaussian-distributed jitter with values that followed the real-valued characteristic function:
A mathematical equation, expression, or formula. Object name is pone.0027380.e041.jpg
. (5)
The coefficient of variation of a jittered spike train is:
A mathematical equation, expression, or formula. Object name is pone.0027380.e042.jpg
The SD of the jitter can conveniently be expressed in units of the SD of the original ISIs, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e043.jpg, where An external file that holds a picture, illustration, etc. Object name is pone.0027380.e044.jpg is a dimensionless scaling parameter. Using this substitution, the CV of a jittered spike train becomes:
A mathematical equation, expression, or formula. Object name is pone.0027380.e045.jpg
. (6)
The serial correlation coefficients of the jittered ISIs, expressed in terms of the SCCs of the original ISI sequence, are given by:
A mathematical equation, expression, or formula. Object name is pone.0027380.e046.jpg
(7)
where An external file that holds a picture, illustration, etc. Object name is pone.0027380.e047.jpg is the Kronecker delta function: An external file that holds a picture, illustration, etc. Object name is pone.0027380.e048.jpg for An external file that holds a picture, illustration, etc. Object name is pone.0027380.e049.jpg, and 0 otherwise. The power spectral density function of a jittered spike train is:
A mathematical equation, expression, or formula. Object name is pone.0027380.e050.jpg
(8)
where An external file that holds a picture, illustration, etc. Object name is pone.0027380.e051.jpg is the PSD of the original spike train. For normally distributed jitter with the characteristic function given by Eq.(5), the PSD of a jittered spike train will be:
A mathematical equation, expression, or formula. Object name is pone.0027380.e052.jpg
(9)
The cross-spectral density of a stimulus and a jittered spike train is:
A mathematical equation, expression, or formula. Object name is pone.0027380.e053.jpg
(10)
The SR coherence of a jittered spike train is:
A mathematical equation, expression, or formula. Object name is pone.0027380.e054.jpg
(11)
and the RR coherence of a jittered spike train is:
A mathematical equation, expression, or formula. Object name is pone.0027380.e055.jpg
(12)
Equations 6–12 are exact and allow an investigator to calculate the metrics of a jittered spike train while bypassing the actual numerical procedure of jittering. For example, Eq. 11 with Eq. 1 allow analytical computation of the lower bound estimate of the mutual information rate for jittered spike trains, based solely on measures of coherence, firing rate and PSD of the original spike train.
Model neurons with gamma-distributed ISIs
In a class of models, we simulated a spike train as a renewal process, i.e., where all ISI durations are independent, with an ISI probability density function (PDF) given by the gamma distribution, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e056.jpg, where the parameter L, called the order of the gamma distribution, sets its shape, and θ is a scaling parameter. The CV, mean ISI, and the ISI variance are, respectively, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e057.jpg, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e058.jpg, and An external file that holds a picture, illustration, etc. Object name is pone.0027380.e059.jpg. The PSD of the gamma spike train can be calculated exactly as:
A mathematical equation, expression, or formula. Object name is pone.0027380.e060.jpg
(13)
For large values of L (i.e., small values of CV), the PSD shows a sharp peak at a frequency corresponding to the mean firing rate, and smaller peaks at higher harmonics. Thus, for large L, the gamma neuron serves as a simple model of a regularly firing neuron. The special case where An external file that holds a picture, illustration, etc. Object name is pone.0027380.e061.jpg corresponds to a Poisson (random) spike train with a uniform power spectrum An external file that holds a picture, illustration, etc. Object name is pone.0027380.e062.jpg.
A stimulus, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e063.jpg, was introduced as Gaussian modulation of the firing rate, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e064.jpg, where An external file that holds a picture, illustration, etc. Object name is pone.0027380.e065.jpg is the spontaneous firing rate, set to 50 or 100 spk/s. The Gaussian stimulus was band-limited to below a cutoff frequency An external file that holds a picture, illustration, etc. Object name is pone.0027380.e066.jpgHz, and had a flat PSD, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e067.jpg, where A is the stimulus SD. Using this stimulus for a total duration of 600 s, we numerically generated sequences of ISIs for a rate-modulated gamma distribution.
Phase model for neuronal oscillator
The dynamics of a periodically firing neuron responding to a stimulus, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e068.jpg, can be represented in terms of its phase [28], [29], [38], [39] as:
A mathematical equation, expression, or formula. Object name is pone.0027380.e069.jpg
(14)
where An external file that holds a picture, illustration, etc. Object name is pone.0027380.e070.jpg is the phase, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e071.jpg is the spontaneous firing rate, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e072.jpg is a 2 π-periodic function known as the phase resetting curve (PRC), An external file that holds a picture, illustration, etc. Object name is pone.0027380.e073.jpgis a white Gaussian noise, and D is the noise intensity. Such a model generates a spike when the phase variable crosses the value of An external file that holds a picture, illustration, etc. Object name is pone.0027380.e074.jpg with a positive slope. We used the so-called type-I PRC where An external file that holds a picture, illustration, etc. Object name is pone.0027380.e075.jpg, with spontaneous rate An external file that holds a picture, illustration, etc. Object name is pone.0027380.e076.jpgspk/s. The noise intensity was An external file that holds a picture, illustration, etc. Object name is pone.0027380.e077.jpg, so that in the absence of any stimulus the model generated a spike train with CV = 0.045. The stimulus s(t) was identical to the one used for the gamma neuron, i.e. Gaussian, band-limited to below 20 Hz, with SD = A. The equation for ϕ was solved numerically using the explicit Euler method with a time step of 0.01 ms.
We developed an analytical framework (Methods, Appendix S1) to investigate the effects of jitter on spike train metrics commonly used to analyze sensory encoding and higher order statistics. We applied this analytical framework to data from model neurons and also experimental data from two types of sensory afferent neurons, to test the efficacy of jitter in distinguishing rate coding from temporal coding.
Examples of jitter influences on stimulus encoding
To set the stage, Fig. 1 illustrates the effects of jitter on experimental data from two representative examples of afferent neurons stimulated by relatively weak external Gaussian noise. The first example is a turtle posterior canal afferent (PCA) stimulated by mechanical indentation of the posterior semicircular duct. The stimulus had a SD of 6.7 µm and was band limited with an upper cutoff frequency An external file that holds a picture, illustration, etc. Object name is pone.0027380.e078.jpg Hz (see Methods). The second example is a paddlefish ampullary electroreceptor afferent (EA) stimulated by a spatially uniform electric field with SD = 0.70 µV/cm and cutoff frequency An external file that holds a picture, illustration, etc. Object name is pone.0027380.e079.jpg Hz. The PCA, with a background firing rate of 21.6 spk/s and CV = 0.21, had a higher intrinsic noise level than the EA which had a background firing rate of 49.0 spk/s and CV = 0.13. These differences are representative of the two afferent populations (see legend of Fig. 2).
Figure 1
Figure 1
Effects of external spike time jitter on stimulus-response (SR) coherence.
Figure 2
Figure 2
Reduction of mutual information and coding fraction due to spike time jitter.
SR coherence functions for these two afferents illustrate the effects of jitter in the frequency domain. For the PCA, jitter with SD≥10 ms resulted in a significant decrease of SR coherence, but smaller (e.g. 3 ms) jitter had minimal effects (Fig. 1a). In contrast, the EA showed a high sensitivity to jitter: SDs as small as 3–5 ms suppressed coherence dramatically (Fig. 1b). Note the excellent correspondence between the results obtained by numerically jittering the afferent spike times (shading) and the results obtained by analytic calculations using Eq. 11 (dotted and dashed lines).
Jitter analysis applied to samples of 10 PCAs and 10 EAs confirmed that, for a given jitter SD, PCAs were less sensitive to jitter than EAs. For the original responses, the lower bound of the mutual information rate (ILB) was An external file that holds a picture, illustration, etc. Object name is pone.0027380.e087.jpg bit/s for the PCA sample, and An external file that holds a picture, illustration, etc. Object name is pone.0027380.e088.jpg bit/s for the EA sample, while the coding fraction (An external file that holds a picture, illustration, etc. Object name is pone.0027380.e089.jpg) was An external file that holds a picture, illustration, etc. Object name is pone.0027380.e090.jpg for the PCAs, and An external file that holds a picture, illustration, etc. Object name is pone.0027380.e091.jpg for the EAs (mean ± SD). Small jitter (SD = 5 ms) had little effect on PCAs, producing a An external file that holds a picture, illustration, etc. Object name is pone.0027380.e092.jpg% reduction of ILB values, and a An external file that holds a picture, illustration, etc. Object name is pone.0027380.e093.jpg% reduction of An external file that holds a picture, illustration, etc. Object name is pone.0027380.e094.jpg values (circles, Fig. 2a,d). However, for the EAs (triangles, Fig. 2a,d), jitter of identical SD had a much larger effect, reducing ILB by An external file that holds a picture, illustration, etc. Object name is pone.0027380.e095.jpg% and An external file that holds a picture, illustration, etc. Object name is pone.0027380.e096.jpg by An external file that holds a picture, illustration, etc. Object name is pone.0027380.e097.jpg%. Because the SD of this jitter was much shorter than the time scale of the stimulus (An external file that holds a picture, illustration, etc. Object name is pone.0027380.e098.jpg = 100 ms for PCAs, 50 ms for EAs), this result could indicate that EAs employ a temporal code to adequately represent the stimulus. In contrast, precise spike timing seems to be less important for PCAs [31].
The percent reduction for both the lower bound of the mutual information rate and the coding fraction increased sigmoidally with increasing jitter SD (Fig. 2). When jitter SD is expressed in units of time, most effects are seen to occur over a 10-fold size range, although the boundaries of the effective range differs between the two types of afferents (Fig. 2a,d). To determine if jitter effects were simply related to the cutoff frequency of the stimulus, the reductions of ILB and γ were also plotted (Fig. 2b, e) after normalizing the jitter SD to the time scale of the stimulus, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e099.jpg. The fact that the curves for the PCA and EA (in Fig. 2b, e) do not superimpose indicates that jitter effects are not simply related to the proximity of the jitter SD to An external file that holds a picture, illustration, etc. Object name is pone.0027380.e100.jpg. To assess any relationship between jitter size and the intrinsic variability of a neuron, we also normalized jitter SD to the SD of the afferent ISI distributions. This type of normalization resulted in nearly identical information and coding loss curves for the two afferents (Fig. 2c,f). This indicates that the intrinsic variability of a neuron is the principal determinant of its sensitivity to external jitter. For both afferent samples, decreases in the mutual information rate and the coding fraction exceeded 50% at the point where the jitter SD equaled the SD of the afferent ISI distributions. Since both ILB and An external file that holds a picture, illustration, etc. Object name is pone.0027380.e101.jpg showed a similar strong dependence on jitter SD, we consider only ILB in the following sections.
Influences of jitter on serial correlations among ISIs
Adding jitter to a spike train indeed increases the variability of its ISI distribution, which is reflected in the increase of CV with jitter SD according to Eq. 6. For renewal processes of similar mean firing rate, the effect of jitter is stronger for less variable spike trains, i.e., those with smaller CV values (more uniform ISIs).
Jitter also alters the serial correlation coefficients among ISIs in two different ways. First, the striking result of Eq. 7 is that jitter introduces negative serial correlations into ISI sequences that were originally generated by a renewal process. For example, with jitter, the SCC of adjacent ISIs, i.e., where lag, m = 1, becomesAn external file that holds a picture, illustration, etc. Object name is pone.0027380.e102.jpg, where An external file that holds a picture, illustration, etc. Object name is pone.0027380.e103.jpg is the SCC of the adjacent ISIs in original spike train, and An external file that holds a picture, illustration, etc. Object name is pone.0027380.e104.jpg, where An external file that holds a picture, illustration, etc. Object name is pone.0027380.e105.jpg is the jitter SD and An external file that holds a picture, illustration, etc. Object name is pone.0027380.e106.jpg is the SD of the original ISI distribution. If the original spike train is a renewal process, where An external file that holds a picture, illustration, etc. Object name is pone.0027380.e107.jpg, then the first SCC of the jittered spike train, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e108.jpg, becomes negative and approaches −0.5 for large values of ε An external file that holds a picture, illustration, etc. Object name is pone.0027380.e109.jpg. Thus, jitter converts a renewal process to a non-renewal process in which adjacent ISIs are negatively correlated, i.e., short ISIs will tend to be followed by long ISIs and vice versa. This introduction of anti-correlated sequential ISIs can be understood qualitatively as a consequence of the jitter values that are added to spike times being drawn from a distribution with a mean. For any pair of spikes, a large absolute value of jitter added to the first spike time is more likely to be followed by a jitter value closer to the mean being added to the second spike time, due to the phenomenon of regression (reversion) to the mean. Thus, the interval between the spike pair is modified by a pair of numbers that themselves tend to be anti-correlated.
Second, a quite different effect of jitter is observed for spike trains that are originally generated by a non-renewal process. In this case, the non-zero SCC values of the original spike train are suppressed by jitter: An external file that holds a picture, illustration, etc. Object name is pone.0027380.e110.jpg. Since An external file that holds a picture, illustration, etc. Object name is pone.0027380.e111.jpg, this scaling parameter in the denominator means that larger jitter brings the SCC values An external file that holds a picture, illustration, etc. Object name is pone.0027380.e112.jpg closer to zero. This effect is also stronger for less-variable spike trains.
Fig. 3 illustrates these effects of jitter on SCCs of spontaneous spike trains of a typical turtle PCA and a typical paddlefish EA. These data are complementary in that turtle PCAs show renewal statistics [40] whereas paddlefish EAs exhibit extended-range serial correlations due to interaction of two distinct types of embedded oscillators [21], [41]. Jitter introduced a negative SCC at the first lag (m = 1) for the renewal afferent (PCA, Fig. 3a), and the value of this SCC became more negative for larger jitter SD, as expected. By contrast, for the non-renewal afferent (EA, Fig. 3b), larger jitter resulted in greater suppression of SCCs for m>1, as expected.
Figure 3
Figure 3
Influence of spike time jitter on serial correlation coefficients (SCCs).
Influences of jitter on spike train power spectra
The power spectrum of a regularly firing neuron has a main peak at a fundamental frequency corresponding to the mean firing rate, and broader peaks at harmonics of the fundamental, as seen in Fig. 4 (labels F and H) for a model neuron with gamma-distributed ISIs, and in Fig. 5 for spontaneous PCA and EA spike trains. As the CV decreases, the peak at the fundamental frequency becomes narrower and higher (Fig. 4, red line in a1 vs. b1, for CVs of 0.05 vs. 0.18), other factors being equal. The discharges of the gamma neuron model (Fig. 4) and the PCA (Fig. 5a1) are both renewal processes, i.e. lacking any serial ISI correlations. In contrast, the EA's spontaneous discharge is non-renewal due to the interaction of multiple internal oscillators [21], and so the PSD of the EA (Fig. 5a2) shows several additional peaks besides at the mean firing rate (asterisk at 44.3 Hz), including one due to the epithelial oscillations (dot at 26 Hz), and peaks at second-order combinations of these fundamentals (44.3±26 Hz) [21].
Figure 4
Figure 4
Effect of spike time jitter on the power spectral density (PSD) of spontaneous neural discharge of two different model gamma neurons calculated from Eq.13.
Figure 5
Figure 5
Effect of spike time jitter on PSD of spontaneous discharge of one PCA (a1) and one EA (a2).
As Figs. 4 and ​and55 show, the PSDs of original and jittered spike trains converge for extremely low frequencies (An external file that holds a picture, illustration, etc. Object name is pone.0027380.e123.jpg). They also converge for very high frequencies (An external file that holds a picture, illustration, etc. Object name is pone.0027380.e124.jpg) because the PSD, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e125.jpg, approaches An external file that holds a picture, illustration, etc. Object name is pone.0027380.e126.jpg for original as well as jittered spike trains (Eq.11). At intermediate frequencies, jitter increases power, except around the mean firing rate and its higher harmonics, where power is suppressed by the jitter. In other words, external jitter both enhances intrinsic variability and suppresses intrinsic periodicity of discharges.
A non-obvious outcome of our analysis is that, for a regularly firing neuron, even small millisecond-level jitter results in a dramatic increase in PSD power at frequencies below the mean firing rate, if An external file that holds a picture, illustration, etc. Object name is pone.0027380.e127.jpg. This increase of power is important for our purposes because the major fraction of stimulus power often lies in this low frequency range. To further quantify this, we considered the power ratio of the PSDs of the jittered and original spike trains, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e128.jpg, where An external file that holds a picture, illustration, etc. Object name is pone.0027380.e129.jpg indicates a power gain produced by jitter. The power enhancement at low frequencies is more dramatic for afferents with more regular discharge, as seen in comparing An external file that holds a picture, illustration, etc. Object name is pone.0027380.e130.jpg for the gamma neuron with CV = 0.05 (Fig. 4a2) or CV = 0.18 (Fig. 4b2).
Small (2–5 ms) jitter caused analogous changes in the PSD of an EA (Fig. 5a2), including a more than 10-fold power gain at low frequencies (Fig. 5b2). Similar but less pronounced effects were observed for the PCA for small jitter (2–5 ms) (Fig. 5a1 and 5b1). This EA and PCA had matching CV = 0.18. Larger jitter (10 ms) completely abolished this EA's PSD peaks.
Our analytical results provide a clear explanation of this power gain at low frequency, as indicated by the excellent correspondence between numerical results from direct jittering procedures and analytical results from Eq. 9 (gray shading vs. superimposed lines in Figs. 5a1 and 5a2). For low frequencies and small jitter, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e131.jpg, the Taylor expansion of Eq. 9 with terms up to An external file that holds a picture, illustration, etc. Object name is pone.0027380.e132.jpg gives a PSD for the jittered spike train of: An external file that holds a picture, illustration, etc. Object name is pone.0027380.e133.jpg. The first term in this expansion indicates an increase in power at low frequencies for the jittered spike train, proportional to the square of the frequency. This is visible in Figs. 4a1 and 4b1 in the steep initial slopes of jittered curves, just above zero frequency. A similar Taylor expansion of the power ratio gives:
A mathematical equation, expression, or formula. Object name is pone.0027380.e134.jpg
(15)
which shows that at low frequencies the power ratio also scales as the square of the frequency.
The magnitude of the power gain is determined jointly by the jitter SD and by the variability of the original spike train, i.e. the term inside brackets in Eq.15. At zero frequency, the PSD of a spike train is determined by the CV, the mean firing rate, and the sum of the SCCs: An external file that holds a picture, illustration, etc. Object name is pone.0027380.e135.jpg [33]. For a regular neuron with a small CV, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e136.jpg at low frequencies (An external file that holds a picture, illustration, etc. Object name is pone.0027380.e137.jpg). Thus, from Eq. 15, the increased power gain due to jitter is more dramatic both for more-regular neurons with smaller values of CV, and for neurons with higher values of the mean firing rate.
If present, negative SCCs further suppress low-frequency variability, resulting in reduced spectral power at low frequencies, compared to a renewal process of equal CV and mean firing rate [42], [43]. Consequently the presence of negative SCCs will cause the effects of jitter to be stronger even if the CV is relatively large, e.g. for the EA in Figs. 3 and ​and55.
Comparison of jitter effects on SR coherence and transfer functions
Besides reshaping the PSD of spontaneous spike trains (Figs. 4 and ​and5),5), the addition of jitter also decreases a response metric to external stimulation, the SR coherence. This is easy to see from the definition of the SR coherence function, Eq. 2, where the spike train PSD is in the denominator, such that the jitter-induced increase of power within the frequency band containing the stimulus leads to a decrease of SR coherence.
For low frequencies and Gaussian jitter, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e138.jpg, the SR coherence function of a jittered spike train can be expanded to a Taylor series,
A mathematical equation, expression, or formula. Object name is pone.0027380.e139.jpg
This shows that the reduction of SR coherence by jitter is lessened by neuron-specific variability, which is proportional to An external file that holds a picture, illustration, etc. Object name is pone.0027380.e140.jpg as discussed above. For a given value of jitter SD, the reduction of SR coherence and the mutual information rate will be greater for a regular spike train, because the magnitude of spike train power An external file that holds a picture, illustration, etc. Object name is pone.0027380.e141.jpg will be smaller (An external file that holds a picture, illustration, etc. Object name is pone.0027380.e142.jpg) within the low frequency band of a stimulus (An external file that holds a picture, illustration, etc. Object name is pone.0027380.e143.jpg). This explains the difference between the magnitude of the jitter effects at low frequencies observed for the PCA and EA in Fig. 1.
On the other hand, another metric of responses to stimulation, the transfer function An external file that holds a picture, illustration, etc. Object name is pone.0027380.e144.jpg, is less affected by jitter. An external file that holds a picture, illustration, etc. Object name is pone.0027380.e145.jpg is a ratio expressing the response magnitude of a linear system relative to the power of a stimulus, at different frequencies. It is not normalized to the spike train PSD, and so is less affected by jitter. The expansion of An external file that holds a picture, illustration, etc. Object name is pone.0027380.e146.jpg to a Taylor series for An external file that holds a picture, illustration, etc. Object name is pone.0027380.e147.jpg, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e148.jpg, shows that the effect of jitter on the transfer function does not depend on the variability of the original spike train at all, and is small for An external file that holds a picture, illustration, etc. Object name is pone.0027380.e149.jpg. This has been observed experimentally [24].
Considerations from Linear Response Theory
Further insight into the effects of jitter came from using linear response theory [44] to approximate the SR coherence function [28], [43], [45]. In this approach, for weak stimuli, the PSD of a stimulated spike train is approximated as the sum of the PSD of the spontaneous discharge An external file that holds a picture, illustration, etc. Object name is pone.0027380.e150.jpg and the PSD of the stimulus weighted with the square of the transfer function An external file that holds a picture, illustration, etc. Object name is pone.0027380.e151.jpg: An external file that holds a picture, illustration, etc. Object name is pone.0027380.e152.jpg. Consequently the SR coherence becomes:
A mathematical equation, expression, or formula. Object name is pone.0027380.e153.jpg
At low frequencies and small jitter magnitudes, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e154.jpg, only the term for the PSD of the spontaneous discharge is strongly affected by the jitter. Thus, for weak stimuli, the suppression of SR coherence at low frequencies by jitter is explained entirely by jitter's effect on the PSD of spontaneous discharge, without any reference to the stimulus.
Model neuron with gamma-distributed ISIs
To demonstrate explicitly that information carried by rate-modulation of a spike train is sensitive to small external jitter, we constructed spike trains from ISI sequences generated from a gamma distribution, with the spike rate modulated by a slow Gaussian stimulus (Methods). To mimic the situation of the mammalian vestibular afferents studied in Sadeghi et al. [24], we constructed spike trains with a spontaneous rate of An external file that holds a picture, illustration, etc. Object name is pone.0027380.e155.jpg spk/s and An external file that holds a picture, illustration, etc. Object name is pone.0027380.e156.jpg. The firing rate was modulated by Gaussian noise with a cutoff frequency of An external file that holds a picture, illustration, etc. Object name is pone.0027380.e157.jpg Hz (Methods). Fig. 6a shows that jitter with SD as small as 1 ms significantly suppressed the SR coherence and, consequently, the mutual information rate. This demonstrates clearly that reduced stimulus encoding resulting from the addition of small external spike time jitter can be observed in the absence of any temporal code.
Figure 6
Figure 6
Effects of small-magnitude spike time jitter on stimulus–response (SR) coherence for a model gamma neuron.
How does stimulus magnitude affect the sensitivity of a spike train to external jitter? For a weak Gaussian noise stimulus alone, the low-frequency power in the spike train will increase with stimulus amplitude and the CV of the spike train will increase quadratically with stimulus SD (Fig. 7a). The effect of jitter alone can be represented as an increase in ISI variability (CV). According to our analysis (above, Fig. 2c,f), the sensitivity of a spike train to external jitter decreases as the variability of the original spike train increases. Thus, our analysis predicts that the effects of jitter of a given magnitude will become smaller as stimulus SD increases. This prediction was borne out for a model gamma neuron (Fig. 6b), and for an EA (Fig. 7b): for a fixed value of jitter SD, the percentage of information loss due to the jitter decreased for larger values of stimulus SD.
Figure 7
Figure 7
Effect of stimulus SD on coefficient of variation (CV) and mutual information.
Effect of jitter on nonlinear responses
Finally, we wish to determine if the effects of jitter on linear responses can be dissociated from the effects of jitter on non-linear responses, where precise spike timing carries information in addition to that carried by firing rate. As the amplitude of a stimulus grows, the response of a neuron becomes progressively nonlinear, and a linear encoding model is no longer optimal [36]. One approach to revealing such nonlinear responses in neurons is to repeatedly present an identical segment of a noise waveform, so-called “frozen noise”. If the stimulus is strong enough, individual spikes become time locked to particular stimulus features resulting in stereotypical neural responses to repeated stimulus presentations [46]–[49]. We studied the effects of jitter on stimulus-induced synchronization using a phase neuron model (Eq. 14). Fig. 8 shows raster plots of this phase model's spike times in response to repeated presentations of 600 s segments of weak or strong “frozen noise” stimuli, as well as the effects of adding 2 ms jitter. For the weak stimulus alone (8a, upper block), spike times varied across stimulus trials because of intrinsic noise in the system, and the jitter has no apparent effect (8a, lower block). However, for the stronger stimulus (8b, upper block), the phase model's firing was tightly locked to the stimulus so that spike times were synchronized into well-defined temporal patterns, reproduced reliably across the ensemble of stimulus trials. For this stronger stimulus, the 2 ms jitter clearly degraded the spike time synchronization (8b, lower block).
Figure 8
Figure 8
Effects of spike time jitter on linear (a) and nonlinear (b) responses from a phase neuron model (Eq.14).
We quantified the effects of jitter on this cross-trial synchronization by comparing SR and RR coherences for the original and jittered spike trains (Eqs. 2, 4, 11, 12). For a weak stimulus, SR and RR coherence were essentially identical and can hardly be distinguished in Fig. 9a (black vs. blue solid lines), indicating that a linear stimulus encoding model is appropriate. Small (2 ms) spike time jitter resulted in identical and significant reduction of both coherence functions (Fig. 9a, black vs. blue dotted lines), that was nearly complete by the stimulus cutoff at An external file that holds a picture, illustration, etc. Object name is pone.0027380.e160.jpg Hz.
Figure 9
Figure 9
Stimulus-response (SR) and response-response (RR) coherence for the linear and nonlinear regimes shown in Fig. 8.
For a stronger but otherwise identical noise stimulus (Fig. 9b), the response was clearly nonlinear, indicated by RR coherence (blue solid line) being larger than SR coherence (black solid line) across the whole frequency range. The large values of RR coherence at high frequencies up to 500 Hz, well beyond the stimulus cutoff at An external file that holds a picture, illustration, etc. Object name is pone.0027380.e161.jpg Hz, clearly reflected the stimulus-induced spike synchronization and the resulting small trial-to-trial spike time variability seen in Fig. 8b (upper block). For this nonlinear response, small jitter (SD = 2 ms) affected SR and RR coherence in different ways. There was insignificant suppression of SR coherence (Fig. 9b, black dotted line) by such small jitter across the entire frequency range, due to the large intra-trial variability of ISIs, that is, their large CV imposed by the stimulus. On the other hand, RR coherence (blue dotted line) was strongly suppressed, but only in the high-frequency band, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e162.jpg, indicating that small jitter values contaminated only the trial-to-trial spike synchronization.
Spike generation is an inherently noisy process, due to various internal sources of noise in neurons [46], [50]–[52]. Consequently the temporal precision of neuronal responses is always somewhat degraded by intrinsic spike time jitter. Indeed, estimating and removing this inherent jitter from neural responses has been shown to significantly improve stimulus reconstructions and estimates of neural transfer functions computed from spiking responses [9]. A reverse approach is often used to assess the degree to which the application of artificial jitter (noise) to the timing of spikes in a spike train degrades stimulus encoding [6], [20], [24]–[26]. We have developed an analytical framework that allows quantitative assessment of the effects of artificial spike timing jitter on both the spontaneous and stimulus-driven discharges of sensory neurons. This framework allows an efficient and analytical assessment of the effects of jitter on various spike train metrics, without requiring the actual numerical jittering procedure. In particular, our theory applies to information metrics estimated using the linear stimulus reconstruction technique, as in numerous neuroscience studies.
We have applied this analytical framework to experimental data from turtle vestibular and paddlefish electroreceptor afferents, as well as to model neurons. Our main results regarding the effects of added spike time noise on sensory encoding are: 1. Jitter increases the variability of spontaneous discharges, as measured by the CV. However, jitter also drastically reshapes the correlation metrics of a spike train, e.g. serial correlation coefficients, and the power spectrum. 2. The relative size of jitter is a critical parameter, as jitter reduces stimulus-response coherence in proportion to the ratio of the jitter SD to the intrinsic variability of the neuron's discharge. 3. Jitter reduces both the mutual information rate and the coding fraction of neuronal responses, even in cases where information is linearly encoded. For example, a gamma model neuron with the rate modulated by the stimulus (that is, with a strict rate-coding scheme explicitly imposed) is sensitive to jitter in the same manner as sensory afferents. 4. For non-linear responses, where spikes are synchronized to stimulus events, small amounts of jitter that have minimal effect on SR coherence can indeed significantly reduce the cross-trial (RR) coherence of repeated responses, and non-linear encoding. Based on these results, we conclude that the degradation of sensory encoding resulting from added spike time noise, as estimated with the linear reconstruction technique, does not provide unequivocal evidence for a temporal code.
Distortions from jitter
We showed that jitter alters the correlation structure of spike sequences. Jitter introduces negative serial ISI correlations into renewal processes, due to the zero-mean nature of the jitter size distribution. On the other hand, jitter degrades serial ISI correlations that are already present. Such distortions by jitter have not been described previously, but since the distinction between renewal and nonrenewal processes is considered fundamental by computational neuroscientists (reviewed in [53], [54]), these distortions may complicate the interpretation of jitter effects. Although negative serial correlations have been shown to reduce background noise in neural discharges [42], [43], [55]–[58], which can enhance linear encoding, the full impact of ISI correlations on CNS processing of sensory information and animal behavior is not known.
What determines a neuron's susceptibility to small external jitter, and what is a meaningful criterion for “small”?
The rationale behind the use of jitter to demonstrate temporal encoding is that the jitter SD, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e163.jpg, is assumed to be small relative to the duration of the neuron's encoding window. Thus, jitter alters the position of spikes within the window (a temporal code), but not the number (a rate code). In the absence of any direct information about the length of the encoding window, the jitter amplitude is typically chosen to be smaller than the characteristic time scale of the stimulus. For a Gaussian noise stimulus that is band-limited to a cutoff frequency An external file that holds a picture, illustration, etc. Object name is pone.0027380.e164.jpg, this criterion for small jitter is then: An external file that holds a picture, illustration, etc. Object name is pone.0027380.e165.jpg, or An external file that holds a picture, illustration, etc. Object name is pone.0027380.e166.jpg. For example, jitter with an SD of 3 ms would be considered to be small for a stimulus with a cutoff frequency of 20 Hz, since An external file that holds a picture, illustration, etc. Object name is pone.0027380.e167.jpg.
In our results, jitter of this small magnitude resulted in a significant reduction of information measures for electroreceptor afferents, but not for vestibular afferents (Fig. 2a,d). This appears consistent with the degree of jitter-induced suppression of information metrics being determined by the intrinsic variability of a spontaneously active neuron, such that neurons with a more regular discharge appear more susceptible to small jitter.
The suppressive effect of jitter on stimulus encoding becomes obvious when the jitter SD is approximately equal to the standard deviation of the neuron's ISI distribution, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e168.jpg, that is, when An external file that holds a picture, illustration, etc. Object name is pone.0027380.e169.jpg (Fig. 2c,f). This result establishes a completely different criterion for what constitutes small jitter. For the example of the electroreceptor afferents with An external file that holds a picture, illustration, etc. Object name is pone.0027380.e170.jpg3–5 ms, jitter with An external file that holds a picture, illustration, etc. Object name is pone.0027380.e171.jpg ms is not small by the new criterion, and indeed significantly reshapes the power spectrum and the SR coherence function of the afferent. On the other hand, for the vestibular afferents, a 3 ms jitter is small relative to its An external file that holds a picture, illustration, etc. Object name is pone.0027380.e172.jpg value (10–40 ms), and consequently it does not significantly affect the afferent's response.
Normalizing the magnitude of the jitter to the value of An external file that holds a picture, illustration, etc. Object name is pone.0027380.e173.jpg, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e174.jpg, provides a universally applicable means of scaling the jitter magnitude for purposes of evaluating its effects on information theoretic metrics, as shown in Fig. 2c,f. This suggests that the effect of jitter is essentially independent of spontaneous discharge regularity per se, if normalized to it, and that the appropriate criterion for considering jitter to be small is unrelated to the time scale of the stimulus, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e175.jpg. Thus, if scaled appropriately, artificial spike time noise has consistent effects on sensory encoding no matter whether a neuron's spontaneous firing is noisy or highly periodic, high or low frequency, renewal or non-renewal.
Why does small-amplitude jitter affect the encoding of low-frequency stimuli?
We examined this question in detail. With the linear reconstruction technique, both the lower bound estimate of mutual information rate and the coding fraction are expressed in terms of SR coherence (Eqs. 1–3). Thus, the suppression of information encoding due to artificial jitter can be understood entirely in terms of the relationship between SR coherence and jitter SD. Our analysis shows that for low frequencies and small jitter amplitudes, i.e. when the product of the frequency and the jitter SD is small (An external file that holds a picture, illustration, etc. Object name is pone.0027380.e176.jpg), the suppression of SR coherence is due to an increase in the power spectral density of the jittered response, while the cross-spectrum is relatively unaffected. That is, artificial jitter significantly enhances power at frequencies lower than the mean firing rate of a regularly firing neuron, and much lower than the inverse of the jitter SD. For progressively weaker stimuli, this jitter-induced power gain is increasingly stimulus independent. Since the response PSD increases without any corresponding increase in the cross-spectrum, the SR coherence is reduced for the jittered responses, and consequently the lower bound estimate of the mutual information rate, and the coding fraction, are reduced also.
The same argument can be made for the optimal reconstruction filter, An external file that holds a picture, illustration, etc. Object name is pone.0027380.e177.jpg, which is calculated as the ratio of the cross-spectrum to the response PSD (Methods), in the frequency domain. Here also, jittered responses exhibit enhanced power in the background noise, which reduces the magnitude of the reconstruction filter, and leads to a reduction of the mutual information rate and the coding fraction. Thus, the sensitivity of an encoding process to small-amplitude jitter, as estimated with the linear reconstruction technique, can be explained completely by jitter-induced transformations of the response spectral characteristics.
The example of Fig. 9a clearly shows that for a weak stimulus and a neural response in which all information is linearly encoded, dramatic suppression of the lower bound mutual information rate by 2 ms jitter can occur. Our example with a rate-modulated gamma neuron (Fig. 6a) further demonstrates that small jitter has essentially the same effect on responses that are explicitly rate encoded as it does on the afferent neurons used in this study. Thus, the sensitivity to small external jitter of stimulus encoding, as estimated by linear reconstruction, cannot per se be taken as evidence for a temporal encoding scheme.
Jitter effect on nonlinear encoding
The effects of jitter on linear and non-linear encoding were quite distinct. For non-linear responses manifested as stimulus-induced synchronization of neuronal firing (as Fig. 8b), estimates of the lower bound of the mutual information rate were essentially unaffected by jitter (Fig. 9b). Instead, jitter dramatically reduced the cross-trial (RR) coherence between responses to a repeated stimulus, but only at frequencies above the stimulus band, corresponding to the time scale of the trial-to-trial spike synchronization induced by the stimulus. Thus, jitter does disrupt coding schemes based on spike timing, when present.
A full assessment of the effects of jitter requires that it also be applied to measures of non-linear encoding, such as direct estimates of the mutual information rate [7], [8], [59], or RR coherence and the upper bound of the mutual information rate associated with it [23], [26], [37].
The variability of neurons with sparse spontaneous activity, e.g. thalamic neurons [8] or whisker primary afferents [25], [27], as estimated from a single presentation of a noise stimulus, is mainly determined by the stimulus itself. Although our theory also applies to such neurons, the use of cross-trial (RR) variability to assess intrinsic variability may be more appropriate than the CV metric used here for afferents with robust spontaneous activity. Indeed, for neurons with sparse spontaneous activity, lower bound information estimates may significantly underestimate the true mutual information rate [10]. Measures of nonlinear encoding must be implemented instead [7], [8], [26], [37], [59].
Supporting Information
Appendix S1
Derivation of relationships between original and jittered spike trains.
(DOC)
Click here for additional data file.(151K, doc)
Acknowledgments
The authors are grateful to Ellengene Peterson for helpful comments on the manuscript. A.B.N. is grateful to Benjamin Lindner for valuable discussions.
Footnotes
Competing Interests: The authors have declared that no competing interests exist.
Funding: This work was supported by National Institutes of Health (NIH) DC 04922 and DC 05063. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
References
1. Theunissen F, Miller JP. Temporal encoding in nervous systems: a rigorous definition. J Comput Neurosci. 1995;2:149–162. [PubMed]
2. Dayan P, Abbott LF. Cambridge, Mass: Massachusetts Institute of Technology Press; 2001. Theoretical Neuroscience: Computational and Mathematical Modeling of Neural Systems.
3. Adrian ED, Zotterman Y. The impulses produced by sensory nerve-endings: Part II. The response of a Single End-Organ. J Physiol. 1926;61:151–171. [PubMed]
4. Rieke F, Warland DK, de Ruyter van Steveninck R, Bialek W. Cambridge, Mass: MIT Press; 1997. Spikes: Exploring the Neural Code.
5. Montemurro MA, Panzeri S, Maravall M, Alenda A, Bale MR, et al. Role of precise spike timing in coding of dynamic vibrissa stimuli in somatosensory thalamus. J Neurophysiol. 2007;98:1871–1882. [PubMed]
6. Rokem A, Watzl S, Gollisch T, Stemmler M, Herz AV, et al. Spike-timing precision underlies the coding efficiency of auditory receptor neurons. J Neurophysiol. 2006;95:2541–2552. [PubMed]
7. Nemenman I, Lewen GD, Bialek W, de Ruyter van Steveninck RR. Neural coding of natural stimuli: information at sub-millisecond resolution. PLoS Comput Biol. 2008;4:e1000025. [PMC free article] [PubMed]
8. Foffani G, Morales-Botello ML, Aguilar J. Spike timing, spike count, and temporal information for the discrimination of tactile stimuli in the rat ventrobasal complex. J Neurosci. 2009;29:5964–5973. [PubMed]
9. Aldworth ZN, Miller JP, Gedeon T, Cummins GI, Dimitrov AG. Dejittered spike-conditioned stimulus waveforms yield improved estimates of neuronal feature selectivity and spike-timing precision of sensory interneurons. J Neurosci. 2005;25:5323–5332. [PubMed]
10. Aldworth ZN, Dimitrov AG, Cummins GI, Gedeon T, Miller JP. Temporal encoding in a nervous system. PLoS Comput Biol. 2011;7:e1002041. [PMC free article] [PubMed]
11. Butts DA, Weng C, Jin J, Alonso JM, Paninski L. Temporal precision in the visual pathway through the interplay of excitation and stimulus-driven suppression. J Neurosci. 2011;31:11313–11327. [PMC free article] [PubMed]
12. de Ruyter van Steveninck RR, Lewen GD, Strong SP, Koberle R, Bialek W. Reproducibility and variability in neural spike trains. Science. 1997;275:1805–1808. [PubMed]
13. Gollisch T, Meister M. Rapid neural coding in the retina with relative spike latencies. Science. 2008;319:1108–1111. [PubMed]
14. Huetz C, Philibert B, Edeline JM. A spike-timing code for discriminating conspecific vocalizations in the thalamocortical system of anesthetized and awake guinea pigs. J Neurosci. 2009;29:334–350. [PubMed]
15. Koppl C. Phase locking to high frequencies in the auditory nerve and cochlear nucleus magnocellularis of the barn owl, Tyto alba. J Neurosci. 1997;17:3312–3321. [PubMed]
16. Rathbun DL, Warland DK, Usrey WM. Spike timing and information transmission at retinogeniculate synapses. J Neurosci. 2010;30:13558–13566. [PMC free article] [PubMed]
17. Rieke F, Bodnar DA, Bialek W. Naturalistic stimuli increase the rate and efficiency of information transmission by primary auditory afferents. Proc Biol Sci. 1995;262:259–265. [PubMed]
18. Goldberg JM. The vestibular end organs: morphological and physiological diversity of afferents. Curr Opin Neurobiol. 1991;1:229–235. [PubMed]
19. Goldberg JM. Afferent diversity and the organization of central vestibular pathways. Exp Brain Res. 2000;130:277–297. [PubMed]
20. Engelmann J, Gertz S, Goulet J, Schuh A, von der Emde G. Coding of stimuli by ampullary afferents in Gnathonemus petersii. J Neurophysiol. 2010;104:1955–1968. [PubMed]
21. Neiman AB, Russell DF. Two distinct types of noisy oscillators in electroreceptors of paddlefish. J Neurophysiol. 2004;92:492–509. [PubMed]
22. Tricas TC, New JG. Sensitivity and response dynamics of elasmobranch electrosensory primary afferent neurons to near threshold fields. J Comp Physiol A. 1998;182:89–101. [PubMed]
23. Borst A, Theunissen FE. Information theory and neural coding. Nat Neurosci. 1999;2:947–957. [PubMed]
24. Sadeghi SG, Chacron MJ, Taylor MC, Cullen KE. Neural variability, detection thresholds, and information transmission in the vestibular system. J Neurosci. 2007;27:771–781. [PubMed]
25. Jones LM, Depireux DA, Simons DJ, Keller A. Robust temporal coding in the trigeminal system. Science. 2004;304:1986–1989. [PMC free article] [PubMed]
26. Massot C, Chacron MJ, Cullen KE. Information transmission and detection thresholds in the vestibular nuclei: single neurons vs. population encoding. J Neurophysiol. 2011;105:1798–1814. [PubMed]
27. Jones LM, Lee S, Trageser JC, Simons DJ, Keller A. Precise temporal responses in whisker trigeminal neurons. J Neurophysiol. 2004;92:665–668. [PMC free article] [PubMed]
28. Schleimer JH, Stemmler M. Coding of information in limit cycle oscillators. Phys Rev Lett. 2009;103:248105. [PubMed]
29. Ermentrout GB, Galan RF, Urban NN. Relating neural dynamics to neural coding. Phys Rev Lett. 2007;99:248103. [PMC free article] [PubMed]
30. Gabbiani F. Coding of time-varying signals in spike trains of linear and half-wave rectifying neurons. Network-Computation in Neural Systems. 1996;7:61–85.
31. Rowe MH, Neiman AB. Information analysis of posterior canal afferents in the turtle, Trachemys scripta elegans. Brain Res. 2011 doi: 10.1016/j.physletb.2003.10.071. [PMC free article] [PubMed]
32. Brichta AM, Goldberg JM. Responses to efferent activation and excitatory response-intensity relations of turtle posterior-crista afferents. J Neurophysiol. 2000;83:1224–1242. [PubMed]
33. Cox DR, Lewis PAW. London: Methuen; 1966. The Statistical Analysis of Series of Events.
34. Bialek W, Rieke F, de Ruyter van Steveninck RR, Warland D. Reading a neural code. Science. 1991;252:1854–1857. [PubMed]
35. Theunissen F, Roddey JC, Stufflebeam S, Clague H, Miller JP. Information theoretic analysis of dynamical encoding by four identified primary sensory interneurons in the cricket cercal system. J Neurophysiol. 1996;75:1345–1364. [PubMed]
36. Theunissen F, Miller JP. Temporal encoding in nervous systems: a rigorous definition. J Comput Neurosci. 1995;2:149–162. [PubMed]
37. Chacron MJ. Nonlinear information processing in a model sensory system. J Neurophysiol. 2006;95:2933–2946. [PubMed]
38. Yoshimura K, Arai K. Phase reduction of stochastic limit cycle oscillators. Phys Rev Lett. 2008;101:154101. [PubMed]
39. Teramae JN, Nakao H, Ermentrout GB. Stochastic phase reduction for a general class of noisy limit cycle oscillators. Phys Rev Lett. 2009;102:194102. [PubMed]
40. Aldworth ZN, Dimitrov AG, Cummins GI, Gedeon T, Miller JP. Temporal encoding in a nervous system. PLoS Comput Biol. 2011;7:e1002041. [PMC free article] [PubMed]
41. Neiman AB, Russell DF. Models of stochastic biperiodic oscillations and extended serial correlations in electroreceptors of paddlefish. Phys Rev E Stat Nonlin Soft Matter Phys. 2005;71:061915. [PubMed]
42. Chacron MJ, Longtin A, Maler L. Negative interspike interval correlations increase the neuronal capacity for encoding time-dependent stimuli. J Neurosci. 2001;21:5328–5343. [PubMed]
43. Chacron MJ, Lindner B, Longtin A. Noise shaping by interval correlations increases information transfer. Phys Rev Lett. 2004;92:080601. [PubMed]
44. Hanggi P, Thomas H. Stochastic processes: Time evolution, symmetries and linear response. Physics Reports. 1982;88:207–319. doi: 10.1016/0370-1573(82)90045-X.
45. Chacron MJ, Maler L, Bastian J. Electroreceptor neuron dynamics shape information transmission. Nat Neurosci. 2005;8:673–678. [PubMed]
46. Mainen ZF, Sejnowski TJ. Reliability of spike timing in neocortical neurons. Science. 1995;268:1503–1506. [PubMed]
47. Teramae JN, Tanaka D. Robustness of the noise-induced phase synchronization in a general class of limit cycle oscillators. Phys Rev Lett. 2004;93:204103. [PubMed]
48. Ermentrout GB, Galan RF, Urban NN. Reliability, synchrony and noise. Trends Neurosci. 2008;31:428–434. [PMC free article] [PubMed]
49. Goldobin DS, Pikovsky A. Synchronization and desynchronization of self-sustained oscillators by common noise. Phys Rev E Stat Nonlin Soft Matter Phys. 2005;71:045201. [PubMed]
50. Bryant HL, Jr, Marcos AR, Segundo JP. Correlations of neuronal spike discharges produced by monosynaptic connections and by common inputs. J Neurophysiol. 1973;36:205–225. [PubMed]
51. Berry MJ, Warland DK, Meister M. The structure and precision of retinal spike trains. Proc Natl Acad Sci U S A. 1997;94:5411–5416. [PubMed]
52. Marsalek P, Koch C, Maunsell J. On the relationship between synaptic input and spike output jitter in individual neurons. Proc Natl Acad Sci U S A. 1997;94:735–740. [PubMed]
53. Farkhooi F, Strube-Bloss MF, Nawrot MP. Serial correlation in neural spike trains: experimental evidence, stochastic modeling, and single neuron variability. Phys Rev E Stat Nonlin Soft Matter Phys. 2009;79:021905. [PubMed]
54. vila-Akerberg O, Chacron MJ. Nonrenewal spike train statistics: causes and functional consequences on neural coding. Exp Brain Res. 2011;210:353–371. [PubMed]
55. Nguyen H, Neiman AB. Spontaneous dynamics and response properties of a Hodgkin-Huxley-type neuron model driven by harmonic synaptic noise. Eur Phys J Spec Top. 2010;187:179–187. [PMC free article] [PubMed]
56. Engel TA, Helbig B, Russell DF, Schimansky-Geier L, Neiman AB. Coherent stochastic oscillations enhance signal detection in spiking neurons. Phys Rev E Stat Nonlin Soft Matter Phys. 2009;80:021919. [PubMed]
57. Ratnam R, Nelson ME. Nonrenewal statistics of electrosensory afferent spike trains: implications for the detection of weak sensory signals. J Neurosci. 2000;20:6672–6683. [PubMed]
58. Chacron MJ, Lindner B, Longtin A. Noise shaping by interval correlations increases information transfer. Phys Rev Lett. 2004;92:080601. [PubMed]
59. Strong SP, de Ruyter van Steveninck RR, Bialek W, Koberle R. On the application of information theory to neural spike trains. Pac Symp Biocomput. 1998;621-632 [PubMed]
Articles from PLoS ONE are provided here courtesy of
Public Library of Science