J Neurophysiol 111: 2320 –2331, 2014. First published March 12, 2014; doi:10.1152/jn.00836.2013.

Heterogeneity of intrinsic biophysical properties among cochlear nucleus neurons improves the population coding of temporal information J. Ahn,1 L. J. Kreeger,1 S. T. Lubejko,1 D. A. Butts,1,2 and K. M. MacLeod1,2,3 1 Department of Biology, University of Maryland, College Park, Maryland; 2Neuroscience and Cognitive Science Program, University of Maryland, College Park, Maryland; and 3Center for the Comparative and Evolutionary Biology of Hearing, University of Maryland, College Park, Maryland

Submitted 25 November 2013; accepted in final form 11 March 2014

Ahn J, Kreeger LJ, Lubejko ST, Butts DA, MacLeod KM. Heterogeneity of intrinsic biophysical properties among cochlear nucleus neurons improves the population coding of temporal information. J Neurophysiol 111: 2320 –2331, 2014. First published March 12, 2014; doi:10.1152/jn.00836.2013.—Reliable representation of the spectrotemporal features of an acoustic stimulus is critical for sound recognition. However, if all neurons respond with identical firing to the same stimulus, redundancy in the activity patterns would reduce the information capacity of the population. We thus investigated spike reliability and temporal fluctuation coding in an ensemble of neurons recorded in vitro from the avian auditory brain stem. Sequential patch-clamp recordings were made from neurons of the cochlear nucleus angularis while injecting identical filtered Gaussian white noise currents, simulating synaptic drive. The spiking activity in neurons receiving these identically fluctuating stimuli was highly correlated, measured pairwise across neurons and as a pseudo-population. Two distinct uncorrelated noise stimuli could be discriminated using the temporal patterning, but not firing rate, of the spike trains in the neural ensemble, with best discrimination using information at time scales of 5–20 ms. Despite high cross-correlation values, the spike patterns observed in individual neurons were idiosyncratic, with notable heterogeneity across neurons. To investigate how temporal information is being encoded, we used optimal linear reconstruction to produce an estimate of the original current stimulus from the spike trains. Ensembles of trains sampled across the neural population could be used to predict ⬎50% of the stimulus variation using optimal linear decoding, compared with ⬃20% using the same number of spike trains recorded from single neurons. We conclude that heterogeneity in the intrinsic biophysical properties of cochlear nucleus neurons reduces firing pattern redundancy while enhancing representation of temporal information. cochlear nucleus; optimal linear decoding; reconstruction; temporal coding; spectrotemporal

to distinguish sound features within an everyday cacophony of background noise belies the computational difficulty of representing a complex auditory scene. Temporal envelope modulations in the sound amplitude across the frequency spectrum are fundamental for distinguishing and categorizing human speech, birdsong, and background textural sounds (Joris et al. 2004; Shannon et al. 1995; Singh and Theunissen 2003). Many of the mechanisms involved in the encoding of acoustic temporal phase cues, or “fine structure,” have been elucidated in the avian brain stem neural pathways responsible for calculating interaural time differences involved in sound localization (Carr and Konishi 1990;

OUR SEEMINGLY EFFORTLESS ABILITY

Address for reprint requests and other correspondence: K. M. MacLeod, Dept. of Biology, Univ. of Maryland, College Park, MD 20742 (e-mail: [email protected]). 2320

Overholt et al. 1992; Sullivan and Konishi 1984; Takahashi et al. 1984; for reviews, see Burger and Rubel 2008 and Kuba 2007). In contrast, the spectrotemporal information processing that underlies sound recognition is more poorly understood, but in birds this process begins with cochlear nucleus angularis (NA), which is also the start of the parallel ascending pathway for binaural sound intensity comparisons (Adolphs 1993; Köppl and Carr 2003; Manley et al. 1988; Soares and Carr 2001; Sullivan and Konishi 1984; Takahashi et al. 1984; Warchol and Dallos 1990). Detailed studies of how temporal envelope features are encoded in NA and its ascending targets are key for understanding the encoding of birdsong and other communication and environmental signals (Keller and Takahashi 2000; Narayan et al. 2006; Schneider and Woolley 2010; Steinberg and Peña 2011; Woolley and Casseday 2005; Woolley et al. 2005). Here, we investigated the mechanisms that underlie the temporal reliability of NA neurons during naturalistic input stimulation. Neurons in NA respond with highly reliable patterns when presented with repeated white noise current stimuli in vitro (Bryant and Segundo 1976; Kreeger et al. 2012), and in vivo when presented with repeated frozen acoustic stimulation (Steinberg and Peña 2011), like neurons in other brain areas in vitro (Mainen and Sejnowski 1995; Tripathy et al. 2013) and in vivo (Butts et al. 2010; de Ruyter van Steveninck et al. 1997; Reinagel and Reid 2002). One prediction from these findings is that any two neurons receiving identical, or nearly identical, input should respond with highly correlated spike patterns, assuming these neurons have similar intrinsic properties (Bryant et al. 1973; Ermentrout et al. 2008; Galán et al. 2006; Reyes 2003; Street and Manis 2007). However, if every neuron responded with an identical pattern to a given input, then the information encoded by a population of neurons receiving similar inputs would be highly redundant. A solution to this problem is to introduce some form of neural heterogeneity among the set of synaptic inputs, network connectivity, or intrinsic properties (Friedrich and Laurent 2001; Gupta 2000; Padmanabhan and Urban 2010; Tripathy et al. 2013). Identified neuronal cell types in auditory and other brain areas show an extraordinary diversity of firing properties (Connors and Gutnick 1990; Hodgkin 1948; Markram et al. 2004; Peruzzi et al. 2000), and the cochlear nucleus is no exception (Carr and Soares 2002; Fukui and Ohmori 2003; Golding et al. 1995; Hirsch and Oertel 1988; Oertel et al. 1988, 1990, 2011; Rothman and Manis 2003; Soares et al. 2002; Wu and Oertel 1984). We therefore investigated the relationship between the

0022-3077/14 Copyright © 2014 the American Physiological Society

www.jn.org

TEMPORAL CODING IN THE COCHLEAR NUCLEUS

diversity of the biophysical properties of cochlear nucleus neurons and their encoding of complex stimuli. Using an in vitro brain slice preparation and direct current (DC) injection, we tested how reliably and coherently a pseudo-population (henceforth “population”) of avian cochlear nucleus neurons responded when the neurons were presented with identical fluctuating noise stimuli. We found that the firing correlations across neurons were high, but heterogeneity in the intrinsic properties led to idiosyncratic spike trains in individual neurons. Variations in the spike trains across neurons significantly improved encoding using population information. Thus heterogeneity in the intrinsic biophysical properties of neurons in the cochlear nucleus effectively balances reliability, high population correlation and efficient coding of temporal information in the inputs. MATERIALS AND METHODS

Brain slice preparation. All animal procedures were reviewed and approved by the University of Maryland Institutional Animal Care and Use Committee. Chicken auditory brain stem slices were prepared as described previously (Kreeger et al. 2012). Briefly, chicken embryos incubated 17–18 days were rapidly decapitated, the brain stem was submerged in chilled and oxygenated low sodium artificial cerebral spinal fluid (ACSF) (in mM: 97.5 NaCl, 3 KCl, 2.5 MgCl2, 26 NaHCO3, 2 CaCl2, 1.25 NaH2PO4, 10 dextrose, 3 HEPES, 230 sucrose), and transverse slices (250 ␮m) were made on a vibrating tissue slicer (Leica Microsystems, Wetzler, Germany). Slices were incubated in oxygenated ACSF (in mM: 130 NaCl, 3 KCl, 2 MgCl2, 26 NaHCO3, 2 CaCl2, 1.25 NaH2PO4, 10 dextrose, 3 HEPES) for 30 min at 34°C then held in solution at room temperature. Slices were submerged in a recording chamber and continuously perfused with warmed (29.6 ⫾ 1.4°C), oxygenated ACSF (1–2 ml/min). No synaptic channel blockers were used. Patch-clamp electrophysiology. Whole cell patch-clamp recordings were made from NA cells using infrared/differential interference contrast video microscopy. Initial micropipette resistances were 4 – 8 M⍀ with an internal solution consisting of the following (in mM): 110 potassium-gluconate, 20 KCl, 1 EGTA, 2 MgCl2, 10 HEPES, 2 Na2ATP, 0.3 Na2GTP, 10 phosphocreatine, and 0.1% biocytin. Recordings were made with a Multiclamp 700B patch clamp amplifier (Molecular Devices, Sunnyvale, CA) in current-clamp mode. Stimulation and recordings were controlled by a PC running custom software in IGOR Pro (Wavemetrics, Lake Oswego, OR). Current stimulation and data collection protocols. Flat current steps and noisy current steps were produced using custom routines written in IGOR Pro. The white noise current stimulus S(t) was generated by convolving Gaussian white noise with an exponential function (time constant, 3 ms), which was chosen to simulate the arrival of many small, stochastic and statistically independent synaptic currents, both excitatory and inhibitory (see Kreeger et al. 2012). Briefly, a Gaussian, exponential-filtered (␶ ⫽ 3 ms) noise stimulus was generated using the same “random” seed to generate the same noise token for all experiments (two seeds, one for each noise token). For each neuron in the main data set in this paper (n ⫽ 60 NA neurons, data set in Table 1), we calibrated the stimulus to account for cell-to-cell variation in input resistance and the DC level to fall within the dynamic range of the neuron (Higgs et al. 2006; Kreeger et al. 2012). To determine whether difference in the DC amplitude across neurons influenced the correlation, we collected a data set (n ⫽ 10) using stimuli with identical DC level (300 pA) and observed no differences of the correlation index (CI) using this set vs. the main data set (see methods below). Cell classification. Offline analysis was performed using custom procedures and functions using IGOR Pro, Matlab (Mathworks, Natick, MA) and Microsoft Excel. Classification was determined

2321

Table 1. Shuffled autocorrelogram and cross-correlogram summary statistics

All SAC Population SAC Pairwise CC, all

Peak CI

Delay, ms

Width, ms

n

15.5 ⫾ 8.7 2.63 4.61 ⫾ 2.64

NA NA 2.04 ⫾ 1.54

1.4 ⫾ 0.7 4.55 3.84 ⫾ 1.70

60 neurons 1,660 pairs

Values are means ⫾ SD. CI, correlation index; SAC, shuffled autocorrelation; CC, cross-correlation.

using physiological firing responses to 400- to 500-ms steps and 50to 100-pA step amplitude intervals (Kreeger et al. 2012; Soares et al. 2002). Shuffled autocorrelation and cross-correlation analysis. To quantify the reliability of spiking by a neuron across repeated trials, we calculated a CI using a shuffled autocorrelation (SAC) analysis (Joris et al. 2006; Street and Manis 2007). The SAC analysis is an allintervals histogram across all trials, excluding intervals taken from within the same trials, and divided by the normalizing factor (NF): NF ⫽ N · (N ⫺ 1) · r2⌬␶D where N is the number of trials, r is the mean firing rate, ⌬␶ is the bin width of the correlation function (0.2 ms), and D is duration of the stimulus epoch (4 s). The central peak maximum is reported as the CI, and the width was full-width at half-maximum of the Gaussian fit of the shuffled auto- or cross-correlation (CC) function. To compare firing across neurons or across different noise tokens, we performed CCs that were normalized to firing rate in an analogous method to the SAC, but as a one-sided function to determine lead/lag times, and the delay reported in the text was the peak location along the x-axis. Summary histograms and statistics include results of correlations with central peaks twice the baseline noise amplitude, true of ⬎93% of possible pairwise combinations (1,660 of 1,770). Spike distance metric. The pseudo-population responses generated by two different noise tokens were tested for discriminability using a spike distance metric (van Rossum 2001) in which the “cost” of transforming one train into the other was calculated. To quantify discrimination at different time scales, the time constant of the exponential that was convolved with the spike train was varied (0.01 ms to 1,000 ms). The median Euclidean difference between the two filtered spike trains was used to assign trains to noise 1 or noise 2, and percentage of correctly classified tabulated. The whole procedure was repeated 30 times, and the mean and standard deviation of the classification percentages plotted. Clustering analysis. The neuronal activity was examined for evidence of clustering based on CC data using standard algorithms in MATLAB (Mathworks, Natick, MA). Dendrograms were constructed from the matrix of interneuron “distances,” where the distance was 1/log(CI) of the pairwise CC, and were optimized using the Ward linkage method and divided using the MATLAB default linkage threshold (0.7). Stimulus reconstruction. We used recorded spike trains to infer the stimulus using a standard optimal linear reconstruction approach (Bialek et al. 1991; Borst and Theunissen 1999; Butts et al. 2007; Stanley et al. 1999). To estimate the optimal linear filters, we used the last 32 of 40 repeat trials (discarding the first 8) and 7.2 s of the 8 s of recording over both 4-s stimulus blocks. The remaining 0.8 s was set aside for crossvalidation, which we used to reconstruct the stimulus, given the recorded spike trains over this time. For single neurons, the optimal linear filter h(␶) is calculated from the spike triggered average k(␶), adjusted by the neuron’s response autocorrelation function C(␶), such that k(␶) ⫽ 冱␶= C(␶ ⫺ ␶=) h(␶=). The reverse filter h(␶) can be solved for by expressing k(␶) and h(␶) as vectors: e.g., k ⫽ [k(0) k(1) k(2) . . . ] and h ⫽ [h(0) h(1) h(2) . . . ], where we are using boldface to denote a vector quantity. Expressing the autocorrelation function as a matrix C with corresponding

J Neurophysiol • doi:10.1152/jn.00836.2013 • www.jn.org

2322

TEMPORAL CODING IN THE COCHLEAR NUCLEUS

entries that keep the above equation in register, the optimal linear reconstruction is given by: h ⫽ C⫺1k where C⫺1 is the matrix inverse of C. For each stimulus reconstruction, we used 30 of the 32 spike trains. Because these spike trains are from the same neuron, the optimal linear filter hi corresponding to the ith spike train ri(t) is simply hi ⫽ h/30, and the linear reconstruction is then given by s*(t) ⫽

兺i hi(␶) ri(t ⫺ ␶)

This estimated stimulus is compared with the actual stimulus to calculate quantities such as r2. Because we use 30 repeats to perform each reconstruction (to compare to population reconstructions, below), the values of r2 that we report is the median of 16 different samples of 30 spike trains (from the 32 available). Because 0.8 s is not enough data to calculate reliable measures of the coherence ␥2(f), we use the full 8 s of data for this. Although this quantity is not cross-validated, we expect that it does not represent overfitting, because of the overall agreement of r2 calculated between crossvalidated and non-cross-validated datasets (data not shown). The reconstruction information Irec is then calculated from the coherence using Irec ⫽ 冱f log2 [1 ⫺ ␥2(f)] (Borst and Theunissen 1999). We performed pseudo-population reconstructions using the same approach, although, in this case, the optimal linear filters must take into account response CCs between neurons in addition to their autocorrelations. This involves modifying the vectors k and h and correlation matrix C in Eq. 1 as follows. Given that the spike-triggered average of the ith neuron is written ki, the vector K ⫽ [k1 k2 k3 . . . ]

represents the concatenation of all spike-triggered averages, and likewise with H, which is composed of individual reconstruction filters hi. The overall correlation matrix Q is then block-constructed by placing autocorrelation matrices Cii along the diagonal, and CC matrices Cij between ith and jth neurons at their corresponding locations, such that Q ⫽ [C11 C12 C13 . . . ; C21 C22 C23 . . . ; C31 . . . ]. The individual reconstruction filters hi can then be extracted from H ⫽ Q⫺1 K (which is identical form to h). RESULTS

Reliable firing to temporally modulated stimuli elicits correlated activity across the ensemble of neurons. To determine whether different neurons would respond with temporally correlated spike train patterns when provided identically fluctuating input, we made sequential whole cell patch-clamp recordings from neurons in the avian cochlear NA while directly injecting a current stimulus composed of exponentially-filtered frozen noise superimposed on a DC step. Spike trains from single neurons showed a high degree of reliability across trials with repeated stimulation (Fig. 1A). To quantify this reliability, we computed the CI, taken as the peak value of the SAC (Joris et al. 2006; Street and Manis 2007), an all-interval histogram normalized to the firing rate. Consistent with previous results (Kreeger et al. 2012), the CIs of the SACs were large for all neurons (15.5 ⫾ 8.7, mean ⫾ SD, n ⫽ 60, Fig. 1, B and D), and peak widths were narrow (1.4 ⫾ 0.7 ms).

Fig. 1. Identical noise current stimulation in different neurons in nucleus angularis (NA) results in correlated firing activity. A: spike train response rasters for five different NA neurons to repeated presentations of the 4-s long, filtered Gaussian noise stimulus (top black trace), color coded for each separate neuron, 40 trials each. Bottom: expanded time scale. B: shuffled autocorrelation (SAC) functions corresponding to neurons a– e, top to bottom. Color coding is as in A. Correlations are normalized by firing rate to facilitate comparison across neurons; a value of one indicated uncorrelated spike timing from one trial to the next (see MATERIALS AND METHODS). C: pairwise shuffled cross-correlation (CC) functions. Traces in each plot show the CC between a given neuron (corresponding to the SAC to its left) vs. the other four neurons in this figure (color coded according to second neuron in each pair). D: mean SAC averaged across all 60 neurons to noise 1 shows a large peak (“noise 1”). Responses of single neurons to two uncorrelated noise stimuli show no correlation (flat curve, “noise 1 ⫻ noise 2”). E: mean CC of all pairwise comparisons (n ⫽ 1,770 combinations; CCs with negative delays were flipped to be positive). J Neurophysiol • doi:10.1152/jn.00836.2013 • www.jn.org

TEMPORAL CODING IN THE COCHLEAR NUCLEUS

By treating this group of neurons as a “pseudo-population” (or simply “population”), we can then investigate how a similar population of neurons would have functioned if recorded simultaneously in the intact brain. We compared the sets of spike trains recorded in different neurons to identically fluctuating repeated stimuli. These fluctuations were superimposed on a current pedestal appropriate to the input resistance of each neuron. Close inspection of the spike trains revealed many alignments across cells, but also that each cell responded with its own idiosyncratic response pattern (Fig. 1A, bottom). We quantified the relative spike timing across neurons with a CC function in an analogous manner to the SAC. Pairwise CC analyses showed large, near-zero delay central peaks for a large majority (⬎93%; see MATERIALS AND METHODS) of all possible combinations (Fig. 1C) with CIs ranging from 1.63 to 27.4 (Fig. 2Ai, Table 1). The CIs of the CCs had a median of 3.87, while the median CC width was 3.47 ms (Fig. 2Aii). A consistent lead or lag in the relative spike times resulted in nonzero peak delays, but most peak delays were close to zero: 75% of correlograms had peak delays smaller than 3 ms, with a median peak delay of 1.72 ms (Fig. 2Aiii). When correlograms were averaged across all pairwise combinations, the resulting function had a CI of 3.6, a width of 3.8 ms and a peak delay of 1.0 ms (Fig. 1E). Similar levels of correlation were observed under two conditions: when stimuli were varied across neurons (0.2– 04 nA) to fall in the steepest part of the neuronal dynamic ranges, as shown here, or when DC level was held constant (data not shown). To compare the spike trains across the entire neural ensemble, we assembled a “population raster” by randomly selecting single trains from each neuron’s set (Fig. 2Bi). We then calculated the shuffled correlogram of the population raster. The peak of the “population” correlogram had a CI of 2.63 and width of 4.55 ms (Fig. 2Bii). These data indicate that, based on

intrinsic properties alone, a high degree of cross-neuronal correlation in firing is likely across pairs or groups of neurons receiving identical inputs. We then asked whether the temporal patterning in the neural ensemble was sufficient for discriminating between different input patterns, and, if so, what is the optimal time scale for discrimination? To test this, a second noise stimulus was presented to each neuron with equal fluctuation and DC amplitude but was temporally uncorrelated to the first. This second noise stimulus also evoked reliable and correlated population firing responses (“noise 2” in Fig. 2C), and, as expected, the two sets of spike trains recorded in the same neuron to the different stimuli showed no correlation whatsoever (“noise 1 ⫻ noise 2”, Fig. 1D). To quantify discrimination, we used a spike distance metric and classification procedure on population rasters for each noise stimulus, created as described previously (Fig. 2Ci) (Narayan et al. 2006; Schneider and Woolley 2010; van Rossum 2001). A cost function procedure that converts one spike train to another was used to calculate spike distances between a given train and every other train elicited by each stimulus. Each spike train was classified as “belonging” to noise 1 or noise 2 based on the distance, and the proportion of trains correctly classified was determined. The analysis was repeated while varying the time window of the cost function, which allows the sensitivity of the discrimination for time scale of the firing patterns to be determined: with short time windows the metric behaves like a “coincidence detector,” while with long time windows only differences in mean firing rate are counted. Different stimuli could be classified with ⬎99% accuracy when time windows of 3 to 20 ms were used and fell to chance with time windows greater than 200 ms (Fig. 2Cii). These results indicate that the two ensembles are completely discrim-

12point

Ai

ii

2323

Ci

iii

60

100

150 100

50

50

0

0

0

5

10

15

20

160

Cell #

150

Frequency

Frequency

Frequency

200 200

120 80

noise 1 1 60

noise 2

40 0

CI

5

10

15

0 0

2

Width (ms)

4

6

8

10

1

Delay (ms)

200 ms

ii

Cii 100 % Correct

Cell #

Correlation

B i 60

1 1.6

1.2

Time (ms)

-20

-10

0

10

80

60

40

20

1

10 100 Tau (ms)

1000

Fig. 2. Summary analysis of correlations across neurons in response to identical noise stimuli. A: histograms of peak correlation index (CI; i), width (ii), and peak delay (iii) of the CC functions (n ⫽ 1,660 pairwise comparisons). Bi: population raster composed of spike trains from 60 different NA neurons. Bii: average population correlogram (n ⫽ 100 random iterations). C: responses to two different uncorrelated noise stimuli were discriminable. Ci: population rasters (1 spike train per neuron randomly selected) for noise 1 (red) and noise 2 (blue). Cii: classification of spike trains based on spike distance metric (van Rossum 2001). Values are means ⫾ SD of 30 repeated analyses. Chance classification is 50% correct (gray line). J Neurophysiol • doi:10.1152/jn.00836.2013 • www.jn.org

2324

TEMPORAL CODING IN THE COCHLEAR NUCLEUS

inable when temporal information in the spike trains is considered, but not when firing rate alone is considered. Does cell type define correlated ensembles? One previously documented aspect of heterogeneity among NA neurons is the division of these neurons into cell types based on their in vitro physiological responses to current step stimuli (for detailed descriptions of cell types in NA, see Carr and Soares 2002; Kreeger et al. 2012; Soares et al. 2002). Our laboratory’s previous study showed that the different cell types affected how the input-output functions changed with noise stimuli (Kreeger et al. 2012). To investigate whether these physiological differences influenced the spike patterns in response to the noise stimuli, we assigned 54 NA neurons to one of the five cell types, Tonic I (n ⫽ 9), Tonic II (n ⫽ 20), Tonic III (n ⫽ 10), Damped (n ⫽ 10) and Single-Spiking (n ⫽ 5), while 6 neurons were excluded from analysis as ambiguous. Representative voltage responses to step currents that illustrate each type are shown in Fig. 3A. Were the neurons that were classified as the same cell type more highly correlated than neurons classified as different cell types? To examine this question we looked at the data three ways. First, to visually inspect the data we plotted the CIs from the SAC and pairwise CC analyses in a correlation matrix, with each of the neurons represented in one column and one row and each box representing a pairwise comparison (Fig. 3A). Second, we grouped the CIs of the CCs into cumulative histograms to determine whether the overall distribution of CIs differed between the within-cell-type pairings vs. the across-cell-type pairings (Fig. 3B). Third, we performed a cluster analysis to look for groupings of neurons based on their CC CIs unbiased by cell type (Fig. 3C). The correlation matrix in Fig. 3A grouped the neurons by cell type and then further ordered the rows and columns by the CI of the autocorrelation (Fig. 3A). The CIs of the SACs, represented along the diagonal, were largest for Tonic I and Single spiking neurons and lowest for Tonic III neurons, comporting with previous work (Kreeger et al. 2012). The remainder of the correlation matrix above and below the diagonal contains the CIs of the CCs. To test whether cell type was a factor in the correlations between neurons, we examined the distribution of pairwise correlations between the same cell types (“within-type”, solid

A

lines in Fig. 3B) vs. those between different cell types (“acrosstype”, dashed lines in Fig. 3B). We performed this analysis separately for each cell type (i.e., Tonic I/Tonic I pair distribution vs. Tonic I/non-Tonic I distribution). We found the distributions largely overlapped, suggesting there were no within-type correlation preferences, with the exception of the Tonic I neurons (solid blue line vs. dashed blue lines in Fig. 3B, P ⬍ 0.05, Kolmogorov-Smirnov test). A leftward shift for the within-type distribution was observed for Tonic III neurons, which also had the lowest CI distribution (solid green line in Fig. 3B). The larger correlations amongst Tonic I neurons and the lower correlations amongst Tonic III neurons suggests that the amplitude of the CC is limited by the intrinsic reliabilities of any two individual neurons, rather than due to cell-type-specific similarities in particular firing patterns. While these data suggest that Tonic I neurons, as a group, may fire together as a correlated ensemble, it is not an exclusive ensemble, since Tonic I neurons also correlated strongly with subsets of neurons classified as other cell types. Visually, the correlation matrix in Fig. 3A did not show any obvious patterning that might indicate clustering of the correlations into subpopulations. To test for correlation groupings unbiased by cell type, we performed a clustering analysis using the CI of the pairwise CC (or more precisely, the inverse of the CI) as an estimate of “distance”. The results of this clustering can be pictured with a “dendrogram” plotted in Fig. 3C. The dendrogram is a visual representation of the results of a minimization procedure that grouped neurons based on their distance to one another (see MATERIALS AND METHODS). Each tip, or twig, of the tree represents an individual neuron, and the vertical height of the branch at intersection with the next twig represents the distance between them. Two large “clusters” (indicated by black and gray coloring of the sub-trees) were found. The smaller of the two clusters (containing 11 out of the 54 neurons in the analysis) comprised neurons with the weakest CCs (tallest twigs), including a number of “below criterion” correlation indexes. The larger cluster of neurons (43/54) comprised neurons that had relatively high CCs. We examined the clusters post hoc to determine there was any pattern related to cell type. For visualization, each twig was color-coded by cell type (colored boxes in two staggered rows correspond to

B Tonic I

Tonic II

C Tonic I Tonic II Tonic III Damped Single

Tonic III DampedSingle

1

Tonic I

Tonic II

Tonic III Damped Single

Fraction

0.8 3.5 3 2.5 2 log(CI) 1.5 1 0.5 0

0.6

TonicI TonicII

0.4

TonicIII damped single

0.2 0

101 Cross-correlation peak (log CI)

11 26 52 28 1 7 35 43 9 19 46 25 45 2 15 10 20 18 24 6 33 48 23 42 38 50 32 14 40 29 30 5 31 49 3 12 37 21 44 39 53 27 17 16 41 4 8 47 13 51 36 22 54 34

Fig. 3. A: correlation matrix of the CI for all pairwise comparisons. Each row and column represents a single neuron (n ⫽ 54), grouped by cell type and ordered by the amplitude of the autocorrelations along the diagonal; CCs are symmetrical above and below the diagonal. Hot/cool colors indicate high/low CI. B: cumulative histograms of CI, color-coded by cell type. Solid lines represent the CIs of CCs between the same cell types (within-group); dashed lines represent the CIs of CCs between neurons of one given cell type (color coded) and all other cell types (across-group). C: cluster analysis of the correlation matrix showed two large clusters: one with significant CCs (black lines) and one with poor CCs (gray lines). Each terminal branch corresponds to a neuron and box below, color-coded by cell type as in B. J Neurophysiol • doi:10.1152/jn.00836.2013 • www.jn.org

TEMPORAL CODING IN THE COCHLEAR NUCLEUS

A

B

Optimal linear filter

Spike trains

Stimulus

X Reconstruction estimate

Norm. Stimulus amplitude

the “twigs” in the dendrogram directly above). The large cluster included all five cell types, while the small cluster included four cell types and did not include any Tonic I neurons. Within the larger cluster, there was no evidence of subclustering related to cell type. We will note, however, a small nearest-neighbor effect: neurons that form pairs of “twigs” tended to be of the same cell type (12 of 18 twig pairs) at a rate above chance. Decoding of the temporal fluctuations by individual neurons: optimal linear stimulus reconstruction. The high degree of reliability in the firing of individual neurons suggests they are responding to temporal events in the stimulus, yet each neuron responded somewhat differently. These differences could potentially reduce the redundancy of information carried by different neurons to allow a more complete encoding of the stimulus by the population. To investigate how information about the stimulus is encoded in the firing activity of individual neurons, we used optimal linear reconstruction (Bialek et al. 1991; Borst and Theunissen 1999) to produce an estimate of the original current stimulus from the spike trains (Fig. 4A). For each neuron, the stimulus reconstruction was created using its optimal linear filter (see MATERIALS AND METHODS). These “reverse filters” were generated from a subset of the spike train data and correspond-

ing current stimulus (7.2 of 8 s total of data; Fig. 4B), and then the reconstructed estimate of the stimulus during a reserved 0.8 s of data was generated by convolving each optimal filter with the corresponding spike train data (Fig. 4C, “Recon-1”). Goodness of fit was determined as the CC of the reconstruction and the original stimulus and reported as an “r2” value (coefficient of determination; Fig. 4D). Reconstructions using single neurons accounted for 7–35% of the variation in the stimulus with an average of 22% (SD ⫾ 6%, n ⫽ 60, Fig. 4C). To determine what frequencies the reconstruction estimates were capable of coding, we measured the power spectra and coherence between the noise stimulus and single-neuron reconstructions (Fig. 4, G and H). Coherence was highest at ⬃25 Hz with a peak amplitude of 0.6 (relative to a maximum possible of one) and low beyond 80 Hz (Fig. 4H, blue line). To determine whether different cell types carried different amounts of information about the stimulus, we analyzed coherence and information rate by cell type (Fig. 4I). The accuracy of the reconstruction estimate is strongly influenced by firing rate (Fig. 4D, n ⫽ 60, linear regression, Pearson’s R ⫽ 0.67, P ⬍ 0.001). To adjust for firing rate, we calculated the reconstruction information using the cross-coherence (see MATERIALS AND METHODS) and calculated the information per spike for each cell type (see MATERIALS AND METHODS) (Fig. 4E),

1.2 0.8 0.4 0 −60

−40

−20

0

10

Norm. Stimulus amplitude

Time (ms)

C

Stimulus Recon-1 Recon-30

2

0

−2 0.10

0.15

0.20

0.25

0.30

Time (sec)

0.4 0.3

r2

0.2 0.1

E

4 3 2 1

20

40

60

80

0

100

20

Firing Rate (Hz)

40

60

80

100

H 1.0

Stimulus Recon-1 Recon-30

50

100

150

Frequency (Hz)

200

250

2 1

D S

Cell Type

I

0.8

Coherence

0.6 0.4 0.2

T1 T2 T3 D S

0.6

Recon-1 Recon-30

0.0 0

3

T1 T2 T3

Firing Rate (Hz)

Coherence

Power

16 14 12 10 8 6 4 2 0

4

0

0

0 0

G

5

F

5

Information (bits/spike)

T1 T2 T3 D S

Information (bits/spike)

D

0.4 0.2 0.0

0

50 100 150 200 250 300

Frequency (Hz)

2325

0

50 100 150 200 250 300

Frequency (Hz)

J Neurophysiol • doi:10.1152/jn.00836.2013 • www.jn.org

Fig. 4. Stimulus reconstruction using reverse correlation noise analysis. A: reconstruction procedure: an optimal linear filter was calculated from the spike trains in response to repeated noisy stimulus; the filter was convolved with each spike train and then averaged to generate a reconstructed estimate of the stimulus. B: optimal linear filter for a typical neuron. Spike occurred at 0 ms. C: original stimulus current (gray) and the corresponding reconstruction estimates based on the same neuron as in B (Recon-1, blue) or population of 30 neurons (Recon-30, green). D: reconstruction fits (r2) for single neurons increased with mean firing rate elicited by the stimulus for that neuron. Color coding is by cell type: T1, Tonic I neurons; T2, Tonic II neurons; T3, Tonic III neurons; D, Damped neurons; S, Single-spiking neurons. E: reconstruction information in bits/spike for the same neurons in D, colored by cell type. F: information rates by cell type. Most neurons carried similar information rates, except for lower rate in the Tonic III neurons (P ⬍ 0.05, ANOVA; Tonic III significantly different from Tonic I, P ⬍ 0.05 multiple-comparisons test). Box plot shows the median, interquartile range and extrema. G: power spectra of the original noise stimulus (black) and the reconstruction estimates in C (blue: single neuron; green: 30 neuron population). H: coherence of the reconstruction estimates with the original stimulus (blue: single neuron; green: 30 neuron population). I: coherence plotted by cell type shows similar curves for the Tonic neurons, while single-spiking neurons had better high-frequency coherence and worse lower-frequency coherence.

2326

TEMPORAL CODING IN THE COCHLEAR NUCLEUS

Fig. 5. Stimulus reconstruction and spike timing reliability. Quality of the single-neuron reconstruction (r2) plotted against reliability (A) and timing precision for each neuron (B) is shown. Reliability was the CI, and timing was the width of the SAC for that neuron. The quality of reconstruction was strongly correlated with timing (Pearson’s R ⫽ ⫺0.527), but only weakly correlated with the reliability (R ⫽ 0.122). In contrast, information per spike was correlated with reliability (R ⫽ 0.685; C), but not with timing precision (R ⫽ 0.010; D).

revealing that the information rate differed by cell type (Fig. 4F; P ⫽ 0.0162, one-way ANOVA, degrees of freedom ⫽ 4, F ⫽ 3.37). Tonic III neurons carried the lowest per spike information content (median, 2.06 bits/spike) compared with rates carried by Tonic I and Single-Spiking neurons (2.65 and 2.8 bits/spike, respectively; significantly different from Tonic I under multiple-comparisons test, Bonferroni correction, ␣ ⫽ 0.05). The mean information rate over all neurons was 2.48 ⫾ 0.54 bits/spike (n ⫽ 60). Given that Tonic III neurons had the lowest, and Tonic I neurons the highest, per-spike information rate as well as spike reliability (data not shown; see also Kreeger et al. 2012), we explicitly investigated the relationship between stimulus coding and spike reliability. We compared the reconstruction outcomes with the SAC analysis for each cell. When the fit of the reconstruction (r2) was plotted against reliability (CI of the SAC) for each cell, only a weak correlation was observed (Fig. 5A, all data: Pearson’s R ⫽ 0.122, P ⬎ 0.05; excluding two outliers: Pearson’s R ⫽ 0.349, P ⬍ 0.05), indicating that reliability was a weak predictor of reconstruction fit, while spike precision (width of the SAC) was a better predictor (Fig. 5B, Pearson’s R ⫽ ⫺0.527, P ⬍ 0.01). On the other hand, the per-spike information was highly correlated with reliability (Fig. 5C, Pearson’s R ⫽ 0.685, P ⬍ 0.01), but uncorrelated with precision (Fig. 5D, Pearson’s R ⫽ 0.010, P ⬎ 0.05). These data suggest that reliable spikes are informative spikes, but that overall reconstruction of the stimulus by individual neurons is limited by the firing rate and spike timing precision. Heterogeneity in spike firing properties leads to improved encoding of temporal fluctuations using population information. Despite high reliability in the firing of individual neurons and high interneuronal correlation, there was also a notable heterogeneity in the spike patterns across neurons. To test whether the

heterogeneity would improve stimulus encoding, we generated a reconstruction estimate of the stimulus based on a set of spike trains drawn from the population to compare with that generated from the same number of spike trains sampled from single neurons. We constructed optimal linear filters for each neuron in the context of the population by taking into account the response CCs in addition to its autocorrelation (see MATERIALS AND METHODS), such that each optimal linear filter depends on the other neurons included in the population. One example reconstruction using a population of 30 neurons is shown in Fig. 4 (green traces labeled “Recon-30” in panels C, G, and H). The population estimate showed clear improvement over the estimate based on the same number of spike trains from a single neuron (Fig. 4C) and recovered a broader range of the frequency response (Fig. 4H). The population reconstruction coherence reached a maximum of ⬎90% at ⬃20 Hz, but the largest relative improvement over the single neuron reconstruction occurred between 50 and 150 Hz. To investigate how the reconstruction improved with spike pattern heterogeneity, we created populations by varying the number of different neurons represented in the pool from 1 to 30 while keeping the total number of spike trains constant (Fig. 6A). As the number of unique neurons in the pool increased, the fit of the reconstruction improved significantly

A

n=3

0

n=6

1.0 0

B 0.5

*

*

**

n=15

1.0 0 Time (sec)

n=30

1.0 0

1.0

*

*

*

10

15

30

0.4

r2

0.3

0.2

0.1

0.0

123

56

Number of unique neurons

Fig. 6. Population coding improves the stimulus reconstruction. A: population rasters drawn from different neuron pool number for a total of 30 trials show increasing heterogeneity in the firing patterns. B: reconstruction fit (r2) vs. number of unique neurons in the population group used in the reconstruction estimate (n ⫽ 120 repeated runs per group containing ⱖ 2 unique neurons; n ⫽ 60 for single-neuron group). Increased heterogeneity results in improved reconstruction fit [one-way ANOVA using Welch’s test assuming unequal variances, F(7,313) ⫽ 4.9, P ⬍ 0.001; all groups of n ⱖ 2 neurons were significantly different from single-neuron group, *P ⬍ 0.05, Dunnett multiplecomparisons test]. Box plot represents the median, interquartile range and extrema.

J Neurophysiol • doi:10.1152/jn.00836.2013 • www.jn.org

TEMPORAL CODING IN THE COCHLEAR NUCLEUS

over single-neuron reconstructions, from a mean r2 value of 0.22 ⫾ 0.6 for single neurons to 0.31 ⫾ 0.6 for pairs of neurons, to 0.53 ⫾ 0.03 for a population of 30 neurons (Fig. 6B; P ⬍ 0.05, one-way ANOVA). The information rate also showed a threefold increase using the population spike trains over single-neuron spike trains (data not shown). Using spike trains from pairs of neurons classified as the same cell type improved the reconstruction equally well as using pairs classified as different cell types (data not shown). These data suggest that the heterogeneity in the spiking responses across individual neurons contributes to the encoding of the stimulus by reducing redundancy across the population. Robust population coding of a common signal using nonidentical input fluctuations. Our experiments providing each NA neuron with identically fluctuating inputs allowed us to isolate the effects of intrinsic firing properties on the spiking responses. Several limitations of this experimental approach, however, should be acknowledged. First, only repeated presentations of a stimulus to an individual neuron were truly identical and resulted in extremely high reliability. During sequential recordings, however, all neurons received identical temporal modulations, but the mean current amplitudes were calibrated individually based on each neuron’s input resistance and input-output function (i.e., “f–I curve”). It was possible that these amplitude differences contributed to the spike train variation only indirectly as a result of differences of the intrinsic properties: controlling for absolute stimulus amplitude resulted in slightly worse, although still high, correlations across neurons (data not shown). Second, we were interested as to whether population coding would be improved due to biophysical differences under more naturalistic conditions. During natural auditory processing, any two neurons are unlikely to experience identical synaptic drive, since realistically they receive only partially overlapping sets of inputs (Fig. 7A). To investigate the effect of nonidentical, but correlated, input stimulation on reliability and signal coding, we analyzed the responses of neurons that were stimulated with either repeated presentations of the “signal” alone, or with a set of stimuli generated by adding random current noise to the signal (signal-to-noise ratio, 4:1; see MATERIALS AND METHODS; Fig. 7B). The addition of random noise variation effectively “roved” the amplitude of the stimulus, allowing a test of the hypothesis that biophysical heterogeneity contributes to improved coding, rather than differences in mean level amplitude that were a limitation in the previous experiments. For single neurons, the addition of random noise variation to the stimulus decreased spike reliability by nearly 80% (Fig. 7, C and D). In contrast, for population responses, the addition of random variation decreased the average CC by 42% (Fig. 7, E and F). In the presence of random noise variation, the CI values were similar for single neurons (SACs had a mean CI of 6.52 ⫾ 4.86) and for CCs (mean CI was 4.60 ⫾ 1.39); thus the addition of random noise to the set of responses in single neuron mimicked the level of correlation and spike time jitter found in the population using similar stimuli. If the improvements in population coding were due to biophysical diversity among NA neurons, then we would expect that reconstructions using the “heterogeneous” population (15 different sequentially recorded neurons) would be improved relative to a “homogeneous” population (single neurons presented with signal ⫹ random noise variation, i.e., “冱

2327

stim”; Fig. 7B). If the improvements in coding observed in experiments using identical stimuli were partly or wholly an artifact of uncontrolled variations in amplitude, then we would expect the population effect to be small; that is, the quality of the reconstructions using the heterogeneous and homogeneous (single cell) populations should be similar. When we analyzed the reconstructions in the presence of random noise variation, we found a large population effect: the fit improved using the heterogeneous population by 42% relative to homogeneous (r2 values: heterogeneous, 0.32 ⫾ 0.02 vs. homogeneous/single, 0.22 ⫾ 0.03, P ⬍ 0.001; Fig. 7G). The population effect observed using stimuli with random variation was nearly as large as that using control stimuli with no random variation (⌬r2 for ⫹variation, 0.09; control, 0.11; Fig. 7H). These data support the contention that biophysical heterogeneity amongst neurons contributes to improved encoding in the cochlear nucleus. DISCUSSION

Temporally fluctuating current stimuli elicit highly reliable firing in neurons of the avian cochlear nucleus. We have shown here that the firing was tightly correlated across the population and also found a high degree of heterogeneity between the spike patterns of different neurons. Stimulus-driven spike trains from individual neurons encoded a substantial portion of the temporal modulations in the stimulus. Pooling information across the population improved the coding more than twofold. The improvements in encoding increased with increasing heterogeneity in the spike patterning, due to variations in the spiking responses across different neurons that spanned a greater proportion of the stimulus, reduced redundancy and combined for a more accurate stimulus reconstruction. These experiments using DC injection allowed us to show that the differences in the patterning were largely due to variation in intrinsic firing properties. These results suggest that biophysical heterogeneity in the firing properties of cochlear nucleus neurons within and across cell type is an important potential mechanism for the efficient coding of auditory stimuli. Reliable firing and temporal correlations across neurons efficiently encode salient stimulus components. Highly reliable, stimulus-driven firing suggests that activity faithfully encodes the dynamic fluctuations in the inputs. Our laboratory’s previous study showed that when presented with frozen white-noise current stimuli, NA neurons fired with reliable, stimuluslocked spike patterns (Kreeger et al. 2012), similar to neurons in other brain areas (Mainen and Sejnowski 1995; Street and Manis 2007). In this study we showed that information about the stimulus was efficiently encoded in the firing patterns using optimal linear reconstruction methods (Bialek et al. 1991; Borst and Theunissen 1999; Butts et al. 2007; Stanley et al. 1999). Spiking activity by individual neurons could account for ⬃22% of the stimulus variance, with high coherence with the stimulus in the frequency range of 10 – 80 Hz. Individual neurons had an average information rate of ⬃2.5 bits/spike, comparable to estimates of coding efficiency in other peripheral, auditory or sensory areas (Bialek et al. 1991; Borst and Theunissen 1999; Butts et al. 2010; Rieke et al. 1997; Tripathy et al. 2013). There was a strong positive correlation between a neuron’s reliability and how much information each spike carried: the more reliable the firing, the more informative was each spike.

J Neurophysiol • doi:10.1152/jn.00836.2013 • www.jn.org

TEMPORAL CODING IN THE COCHLEAR NUCLEUS

A

B

random variation

signal

ANF

Σ stim

heterogeneous pop.

{

2328

{

+ NA1

NA2

single neuron

C

D

stimuli

Correlation

0.6

control (signal only)

+variation

35 SAC

0 -20

F

10

Correlation

E control (signal only)

+variation

0 -20

100 ms

+variation 0 20 Delay (ms)

CC

control +variation

0 Delay (ms)

20

0.4

**

(r 2)

pop

** pop

single single

Pop. Effect (Δ r 2)

H

G Reconstruction Fit

control

0

0.2

0.1

0 +Variation

Control

con

+var

Fig. 7. Robust stimulus coding in the presence of interneuronal input variations. A: each NA neuron sums multiple small inputs that may partially overlap with inputs to other NA neurons. B: schematic of experimental protocol using varied input. A standard noisy current (“signal”) was summed with a set of randomly varying noise currents (“random variation”) to generate the set of stimuli presented to the neurons (“冱 stim”). A “heterogeneous” population response was constructed by selecting a single response from each neuron (n ⫽ 15 sequential recordings) presented with a different variation of the 冱 stim (see data shown as red dot raster in E). For comparison, a “homogeneous” population response was constructed from the responses of a single neuron to 15 trials each with a different variation of the 冱 stim (red dot raster in C). “Control” heterogeneous and homogeneous population responses were constructed using the signal alone (not shown in B; blue dot rasters in C and E). C: single-neuron effects of random variation. Top traces: control noise stimulus (blue) overlaid with one example of a noise stimulus with added random variation (red) (signal/noise, 4:1). Rasters show that random variation in the stimulus introduced variability in the timing and likelihood of spiking (red) compared with control (blue). Rows in each condition correspond to different trains from a single neuron. D: average SAC (n ⫽ 15) showed a large decrease in peak amplitude and increase in width with random added noise. E: effect of random variation on population correlations. Stimuli are as in C. Population rasters are in response to presentation of identical signal stimulus (blue). Addition of random variation to the stimulus resulted in some additional variation (red). Rows in each condition correspond to different neurons. F: average cross-correlograms of responses (CC, n ⫽ 105 pairwise comparisons) collected during control stimulation (blue) and during stimulation with added random variation (red). G: optimal linear reconstruction analysis in the presence of random noise variation (⫹variation, red) show a population effect, with larger r2 values for the heterogeneous population (“pop”) than for a homogeneous population of responses from single neurons (“single”) (**P ⬍ 0.001), similar to the population effect seen under control stimulation for these neurons (control, **P ⬍ 0.001) (n ⫽ 15 neurons). H: absolute population effect size (⌬r2 ⫽ mean population r2 ⫺ mean single neuron r2) was similar in the two conditions.

Highly reliable firing has been observed in vivo during recordings of spiking activity in barn owl NA in response to acoustic noise stimuli (Steinberg and Peña 2011). The reliable firing in NA in vivo was related to the extraction of the acoustic envelope, while the cochlear nucleus in the “temporal” pathway, nucleus magnocellularis, had comparatively less reliable firing but was better phase-locked to the fine structure. Neural activity in NA is the likely source, via direct or indirect

ascending projections, of reliable temporal inputs to higher order auditory nuclei that also demonstrate good envelope coding (Keller and Takahashi 2000). Reliable temporal patterning is observed also in the neural activity recorded in vivo in midbrain and forebrain auditory areas in zebra finch and mouse in response to natural communication signals. When temporal information in the neighborhood of 10 ms is used, the neural activity can be used to distinguish between stimuli,

J Neurophysiol • doi:10.1152/jn.00836.2013 • www.jn.org

TEMPORAL CODING IN THE COCHLEAR NUCLEUS

whereas rate information is comparatively uninformative (Holmstrom et al. 2010; Narayan et al. 2006; Schneider and Woolley 2010). These studies suggest that the highly reliable temporal patterning is critical to processing the spectrotemporal envelope information so characteristic of natural communication signals. One would naturally conjecture that, if different reliable neurons were presented with the same temporally modulated stimulus, they should respond similarly, and their patterns of activity would be correlated in time. Our results show that neurons do, in fact, respond with correlated firing to a precision of 3–5 ms. Similar levels of population correlations have been observed among pyramidal neurons of the mammalian dorsal cochlear nucleus (Street and Manis 2007). The reliable firing and correlations occur during the rising phase and peaks of temporal activity. Previous work showed that many NA neurons are particularly sensitive to these rapid fluctuations and fire more when steep fluctuations occur (Kreeger et al. 2012). Similarly, mammalian cochlear nucleus neurons are sensitive to the rate of rise in the voltage and have narrow temporal integration windows, although stellate cells (analogous to the tonic firing neurons in NA) have been described as integrators (McGinley and Oertel 2006). Because NA neurons receive many small inputs from auditory nerve fibers, rapid fluctuations would be expected to arise when there are correlations across the input population. Correlations in the input array are, therefore, propagated to the next level with increased firing and correlations in the cochlear nucleus population output activity. This suggests that cochlear nucleus neurons show synchrony transfer, such that receipt of coactivated inputs results in correlated outputs (Ratté et al. 2013; Reyes 2003). Evidence from in vivo recordings has shown that sound stimulation can drive stimulus-driven synchronous firing (Shivdasani et al. 2009; Woolley et al. 2006). Stimulus correlations in the temporal domain convey salient features in the sounds, such as onsets and amplitude modulations, and, therefore, carry information about the statistical structure of natural auditory scenes. Psychophysical experiments have shown that coincident onsets or rapid coherent amplitude and frequency modulations act as grouping cues to unify a perception of distinct acoustic components as arising from a single source, promoting streaming and segregation (Bregman 1990; McAdams 1989). The physiological basis of these psychoacoustic effects is not yet understood, but temporal coherence in the spiking activity across a wide range of neuronal populations encoding different features could provide a binding cue (Malsburg and Schneider 1986; Shamma et al. 2011). Biophysical heterogeneity reduces information redundancy and improves temporal envelope coding. Highly reliable firing and population correlations indicate efficient coding of salient temporal features. However, if all neurons responded reliably and identically to the same stimulus, resulting in strong temporal correlations, high redundancy would curtail the information transmitted across the population (Sompolinsky et al. 2001; Zohary et al. 1994). Heterogeneity of response or network properties can improve information processing using population codes (Chelaru and Dragoi 2008; Osborne et al. 2008; Shamir and Sompolinsky 2006). Using identical fluctuating stimuli and DC injection in an in vitro preparation, we show that cochlear nucleus neurons demonstrate a remarkable degree of heterogeneity in the firing patterns across different

2329

cells; remarkable since individually, each neuron’s firing activity was exceptionally reliable, precise and fluctuationdriven. Using optimal linear decoding, in this study we showed that pooling the information across populations of auditory neurons improved stimulus coding of the temporal fluctuations, explaining ⬎50% of the stimulus variance with a population of 30 neurons. Furthermore, under more realistic scenarios using inputs that were highly similar but not identical, heterogeneity among the neuronal firing responses still improved the neural coding of a common signal. Studies in the olfactory bulb have shown similar effects that we show here, such that the diversity of the intrinsic properties of mitral cells receiving common inputs improved coding by reducing redundancy and generating heterogeneity in the temporal patterning of activity (Padmanabhan and Urban 2010; Tripathy et al. 2013). Intrinsic biophysical heterogeneity may thus be a common mechanism for improved sensory coding. Heterogeneity within and across physiological cell type. Based on their physiological firing properties in response to step currents, neurons in NA can be divided into five classes: a single-spiking type and several types of repetitive firing neurons (Fukui and Ohmori 2003; Soares et al. 2002). These cell types, however, had relatively little influence on the neuronal firing correlations when stimulated with noisy currents. Instead of grouping into correlated subensembles of similar neurons, the heterogeneity in firing resulted in a correlational structure was distributed and not clustered. The hypothesis that groups of neurons with similar intrinsic properties would selectively form coherently firing ensembles distinct from other neurons with different intrinsic properties (Ermentrout et al. 2008) did not bear out in our data. Pairwise combinations of neurons of the same “cell type” generate reconstruction estimate improvements of equal magnitude as pairwise combinations of neurons from different cell types. The lack of clustering even when ignoring cell type designations suggested that the variation might be more continuous than discrete. These results underscore the idea that, even within defined cell types, there can also lurk substantial variability in the details of the biophysics (Marder and Goaillard 2006). In summary, the reliability of firing and efficient encoding of temporally modulated stimuli by neurons in the avian cochlear nucleus appear to be balanced with the pressure for redundancy reduction achieved through the expression of continuous heterogeneity in the intrinsic biophysical properties within and across cell classes. The intrinsic heterogeneity of membrane mechanisms, as shown here, along with variation in the synaptic input strength, short-term synaptic plasticity and other network heterogeneities could all independently contribute to improving coding at the single neuron and population level. ACKNOWLEDGMENTS Present address of L. J. Kreeger: Neuroscience Graduate Program, University of Texas at Austin, Austin, TX.

GRANTS This work was supported by funding by the National Institute of Deafness and Other Communication Disorders Grant DC10000 (K. M. MacLeod), a grant from the National Organization for Hearing Research (K. M. MacLeod), and a Howard Hughes Undergraduate Research Fellowship (L. J. Kreeger).

J Neurophysiol • doi:10.1152/jn.00836.2013 • www.jn.org

2330

TEMPORAL CODING IN THE COCHLEAR NUCLEUS

DISCLOSURES No conflicts of interest, financial or otherwise, are declared by the author(s).

AUTHOR CONTRIBUTIONS Author contributions: J.A., L.J.K., S.T.L., and K.M.M. performed experiments; J.A., L.J.K., S.T.L., D.A.B., and K.M.M. analyzed data; J.A., D.A.B., and K.M.M. edited and revised manuscript; J.A., L.J.K., S.T.L., D.A.B., and K.M.M. approved final version of manuscript; L.J.K., D.A.B., and K.M.M. conception and design of research; D.A.B. and K.M.M. interpreted results of experiments; D.A.B. and K.M.M. prepared figures; D.A.B. and K.M.M. drafted manuscript.

REFERENCES Adolphs R. Bilateral inhibition generates neuronal responses tuned to interaural level differences in the auditory brainstem of the barn owl. J Neurosci 13: 3647–3668, 1993. Bialek W, Rieke F, de Ruyter van Steveninck RR, Warland D. Reading a neural code. Science 252: 1854 –1857, 1991. Borst A, Theunissen FE. Information theory and neural coding. Nat Neurosci 2: 947–957, 1999. Bregman AS. Auditory Scene Analysis. Cambridge, MA: MIT Press, 1990. Bryant HL, Marcos AR, Segundo JP. Correlations of neuronal spike discharges produced by monosynaptic connections and by common inputs. J Neurophysiol 36: 205–225, 1973. Bryant HL, Segundo JP. Spike initiation by transmembrane current: a white-noise analysis. J Physiol 260: 279 –314, 1976. Burger RM, Rubel EW. Encoding of interaural timing for binaural hearing. In: The Senses: A Comprehensive Reference. Audition, edited by Dallos P, Oertel D. New York: Academic, 2008, vol. 3, p. 613– 630. Butts DA, Desbordes G, Weng C, Jin J, Alonso JM, Stanley GB. The episodic nature of spike trains in the early visual pathway. J Neurophysiol 104: 3371–3387, 2010. Butts DA, Weng C, Jin J, Yeh CI, Lesica NA, Alonso JM, Stanley GB. Temporal precision in the neural code and the timescales of natural vision. Nature 449: 92–95, 2007. Carr CE, Konishi M. A circuit for detection of interaural time differences in the brain stem of the barn owl. J Neurosci 10: 3227–3246, 1990. Carr CE, Soares D. Evolutionary convergence and shared computational principles in the auditory system. Brain Behav Evol 59: 294 –311, 2002. Chelaru MI, Dragoi V. Efficient coding in heterogeneous neuronal populations. Proc Natl Acad Sci U S A 105: 16344 –16349, 2008. Connors BW, Gutnick MJ. Intrinsic firing patterns of diverse neocortical neurons. Trends Neurosci 13: 99 –104, 1990. de Ruyter van Steveninck RR, Lewen GD, Strong SP, Koberle R, Bialek W. Reproducibility and variability in neural spike trains. Science 275: 1805–1808, 1997. Ermentrout GB, Galán RF, Urban NN. Reliability, synchrony and noise. Trends Neurosci 31: 428 – 434, 2008. Friedrich RW, Laurent G. Dynamic optimization of odor representations by slow temporal patterning of mitral cell activity. Science 291: 889 – 894, 2001. Fukui I, Ohmori H. Developmental changes in membrane excitability and morphology of neurons in the nucleus angularis of the chicken. J Physiol 548: 219 –232, 2003. Galán RF, Fourcaud-Trocmé N, Ermentrout GB, Urban NN. Correlationinduced synchronization of oscillations in olfactory bulb neurons. J Neurosci 26: 3646 –3655, 2006. Golding NL, Robertson D, Oertel D. Recordings from slices indicate that octopus cells of the cochlear nucleus detect coincident firing of auditory nerve fibers with temporal precision. J Neurosci 15: 3138 –3153, 1995. Gupta A. Organizing principles for a diversity of GABAergic interneurons and synapses in the neocortex. Science 287: 273–278, 2000. Higgs MH, Slee SJ, Spain WJ. Diversity of gain modulation by noise in neocortical neurons: regulation by the slow afterhyperpolarization conductance. J Neurosci 26: 8787– 8799, 2006. Hirsch JA, Oertel D. Intrinsic properties of neurones in the dorsal cochlear nucleus of mice, in vitro. J Physiol 396: 535–548, 1988. Hodgkin AL. The local electric changes associated with repetitive action in a non-medullated axon. J Physiol 107: 165–181, 1948. Holmstrom LA, Eeuwes LBM, Roberts PD, Portfors CV. Efficient encoding of vocalizations in the auditory midbrain. J Neurosci 30: 802– 819, 2010.

Joris PX, Louage DH, Cardoen L, van der Heijden M. Correlation index: a new metric to quantify temporal coding. Hear Res 216 –217: 19 –30, 2006. Joris PX, Schreiner CE, Rees A. Neural processing of amplitude-modulated sounds. Physiol Rev 84: 541–577, 2004. Keller CH, Takahashi TT. Representation of temporal features of complex sounds by the discharge patterns of neurons in the owl’s inferior colliculus. J Neurophysiol 84: 2638 –2650, 2000. Köppl C, Carr CE. Computational diversity in the cochlear nucleus angularis of the barn owl. J Neurophysiol 89: 2313–2329, 2003. Kreeger LJ, Arshed A, MacLeod KM. Intrinsic firing properties in the avian auditory brain stem allow both integration and encoding of temporally modulated noisy inputs in vitro. J Neurophysiol 108: 2794 –2809, 2012. Kuba H. Cellular and molecular mechanisms of avian auditory coincidence detection. Neurosci Res 59: 370 –376, 2007. Mainen ZF, Sejnowski TJ. Reliability of spike timing in neocortical neurons. Science 268: 1503–1506, 1995. Malsburg C, Schneider W. A neural cocktail-party processor. Biol Cybern 54: 29 – 40, 1986. Manley GA, Köppl C, Konishi M. A neural map of interaural intensity differences in the brain stem of the barn owl. J Neurosci 8: 2665–2676, 1988. Marder E, Goaillard JM. Variability, compensation and homeostasis in neuron and network function. Nat Rev Neurosci 7: 563–574, 2006. Markram H, Toledo-Rodriguez M, Wang Y, Gupta A, Silberberg G, Wu C. Interneurons of the neocortical inhibitory system. Nat Rev Neurosci 5: 793– 807, 2004. McAdams S. Segregation of concurrent sounds. I. Effects of frequency modulation coherence. J Acoust Soc Am 86: 2148 –2159, 1989. McGinley MJ, Oertel D. Rate thresholds determine the precision of temporal integration in principal cells of the ventral cochlear nucleus. Hear Res 216 –217: 52– 63, 2006. Narayan R, Graña G, Sen K. Distinct time scales in cortical discrimination of natural sounds in songbirds. J Neurophysiol 96: 252–258, 2006. Oertel D, Wright S, Cao XJ, Ferragamo M, Bal R. The multiple functions of T stellate/multipolar/chopper cells in the ventral cochlear nucleus. Hear Res 276: 61– 69, 2011. Oertel D, Wu SH, Garb MW, Dizack C. Morphology and physiology of cells in slice preparations of the posteroventral cochlear nucleus of mice. J Comp Neurol 295: 136 –154, 1990. Oertel D, Wu SH, Hirsch JA. Electrical characteristics of cells and neuronal circuitry in the cochlear nuclei studied with intracellular recordings from brain slices. In: Auditory Function, edited by Edelman GM, Gall WE, Cowan WM. New York: Wiley, 1988, p. 313–336. Osborne LC, Palmer SE, Lisberger SG, Bialek W. The neural basis for combinatorial coding in a cortical population response. J Neurosci 28: 13522–13531, 2008. Overholt EM, Rubel EW, Hyson RL. A circuit for coding interaural time differences in the chick brainstem. J Neurosci 12: 1698 –1708, 1992. Padmanabhan K, Urban NN. Intrinsic biophysical diversity decorrelates neuronal firing while increasing information content. Nat Neurosci 13: 1276 –1282, 2010. Peruzzi D, Sivaramakrishnan S, Oliver DL. Identification of cell types in brain slices of the inferior colliculus. Neuroscience 101: 403– 416, 2000. Ratté S, Hong S, De Schutter E, Prescott SA. Impact of neuronal properties on network coding: roles of spike initiation dynamics and robust synchrony transfer. Neuron 78: 758 –772, 2013. Reinagel P, Reid RC. Precise firing events are conserved across neurons. J Neurosci 22: 6837– 6841, 2002. Reyes AD. Synchrony-dependent propagation of firing rate in iteratively constructed networks in vitro. Nat Neurosci 6: 593–599, 2003. Rieke F, Warland D, de Ruyter van Steveninck RR, Bialek W. Spikes: Exploring the Neural Code. Cambridge, MA: MIT Press, 1997. Rothman JS, Manis PB. The roles potassium currents play in regulating the electrical activity of ventral cochlear nucleus neurons. J Neurophysiol 89: 3097–3113, 2003. Schneider DM, Woolley SMN. Discrimination of communication vocalizations by single neurons and groups of neurons in the auditory midbrain. J Neurophysiol 103: 3248 –3265, 2010. Shamir M, Sompolinsky H. Implications of neuronal diversity on population coding. Neural Comput 18: 1951–1986, 2006. Shamma SA, Elhilali M, Micheyl C. Temporal coherence and attention in auditory scene analysis. Trends Neurosci 34: 114 –123, 2011. Shannon RV, Zeng FG, Kamath V, Wygonski J, Ekelid M. Speech recognition with primarily temporal cues. Science 270: 303–304, 1995.

J Neurophysiol • doi:10.1152/jn.00836.2013 • www.jn.org

TEMPORAL CODING IN THE COCHLEAR NUCLEUS Shivdasani MN, Mauger SJ, Rathbone GD, Paolini AG. Neural synchrony in ventral cochlear nucleus neuron populations is not mediated by intrinsic processes but is stimulus induced: implications for auditory brainstem implants. J Neural Eng 6: 065003, 2009. Singh NC, Theunissen FE. Modulation spectra of natural sounds and ethological theories of auditory processing. J Acoust Soc Am 114: 3394 –3411, 2003. Soares D, Carr CE. The cytoarchitecture of the nucleus angularis of the barn owl (Tyto alba). J Comp Neurol 429: 192–205, 2001. Soares D, Chitwood RA, Hyson RL, Carr CE. Intrinsic neuronal properties of the chick nucleus angularis. J Neurophysiol 88: 152–162, 2002. Sompolinsky H, Yoon H, Kang K, Shamir M. Population coding in neuronal systems with correlated noise. Phys Rev E Stat Nonlin Soft Matter Phys 64: 051904, 2001. Stanley GB, Li FF, Dan Y. Reconstruction of natural scenes from ensemble responses in the lateral geniculate nucleus. J Neurosci 19: 8036 – 8042, 1999. Steinberg LJ, Peña JL. Difference in response reliability predicted by spectrotemporal tuning in the cochlear nuclei of barn owls. J Neurosci 31: 3234 –3242, 2011. Street SE, Manis PB. Action potential timing precision in dorsal cochlear nucleus pyramidal cells. J Neurophysiol 97: 4162– 4172, 2007. Sullivan WE, Konishi M. Segregation of stimulus phase and intensity coding in the cochlear nucleus of the barn owl. J Neurosci 4: 1787–1799, 1984.

2331

Takahashi T, Moiseff A, Konishi M. Time and intensity cues are processed independently in the auditory system of the owl. J Neurosci 4: 1781–1786, 1984. Tripathy SJ, Padmanabhan K, Gerkin RC, Urban NN. Intermediate intrinsic diversity enhances neural population coding. Proc Natl Acad Sci U S A 110: 8248 – 8253, 2013. van Rossum MC. A novel spike distance. Neural Comput 13: 751–763, 2001. Warchol ME, Dallos P. Neural coding in the chick cochlear nucleus. J Comp Physiol A 166: 721–734, 1990. Woolley SMN, Casseday JH. Processing of modulated sounds in the zebra finch auditory midbrain: responses to noise, frequency sweeps, and sinusoidal amplitude modulations. J Neurophysiol 94: 1143–1157, 2005. Woolley SMN, Fremouw TE, Hsu A, Theunissen FE. Tuning for spectrotemporal modulations as a mechanism for auditory discrimination of natural sounds. Nat Neurosci 8: 1371–1379, 2005. Woolley SM, Gill PR, Theunissen FE. Stimulus-dependent auditory tuning results in synchronous population coding of vocalizations in the songbird midbrain. J Neurosci 26: 2499 –2512, 2006. Wu S, Oertel D. Intracellular injection with horseradish peroxidase of physiologically characterized stellate and bushy cells in slices of mouse anteroventral cochlear nucleus. J Neurosci 4: 1577–1588, 1984. Zohary E, Shadlen MN, Newsome WT. Correlated neuronal discharge rate and its implications for psychophysical performance. Nature 370: 140 –143, 1994.

J Neurophysiol • doi:10.1152/jn.00836.2013 • www.jn.org