Skip to main content

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Journal of Neuroscience
  • Log in
  • My Cart
Journal of Neuroscience

Advanced Search

Submit a Manuscript
  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE
PreviousNext
Articles, Systems/Circuits

The Neural Representation of Interaural Time Differences in Gerbils Is Transformed from Midbrain to Cortex

Lucile A.C. Belliveau, Dmitry R. Lyamzin and Nicholas A. Lesica
Journal of Neuroscience 10 December 2014, 34 (50) 16796-16808; DOI: https://doi.org/10.1523/JNEUROSCI.2432-14.2014
Lucile A.C. Belliveau
1Ear Institute, University College London, London WC1E 6BT, United Kingdom and
2Graduate School of Systemic Neurosciences, Ludwig-Maximilians-University Munich, 80802 Munich, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Dmitry R. Lyamzin
1Ear Institute, University College London, London WC1E 6BT, United Kingdom and
2Graduate School of Systemic Neurosciences, Ludwig-Maximilians-University Munich, 80802 Munich, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Nicholas A. Lesica
1Ear Institute, University College London, London WC1E 6BT, United Kingdom and
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

Interaural time differences (ITDs) are the dominant cue for the localization of low-frequency sounds. While much is known about the processing of ITDs in the auditory brainstem and midbrain, there have been relatively few studies of ITD processing in auditory cortex. In this study, we compared the neural representation of ITDs in the inferior colliculus (IC) and primary auditory cortex (A1) of gerbils. Our IC results were largely consistent with previous studies, with most cells responding maximally to ITDs that correspond to the contralateral edge of the physiological range. In A1, however, we found that preferred ITDs were distributed evenly throughout the physiological range without any contralateral bias. This difference in the distribution of preferred ITDs in IC and A1 had a major impact on the coding of ITDs at the population level: while a labeled-line decoder that considered the tuning of individual cells performed well on both IC and A1 responses, a two-channel decoder based on the overall activity in each hemisphere performed poorly on A1 responses relative to either labeled-line decoding of A1 responses or two-channel decoding of IC responses. These results suggest that the neural representation of ITDs in gerbils is transformed from IC to A1 and have important implications for how spatial location may be combined with other acoustic features for the analysis of complex auditory scenes.

  • auditory cortex
  • inferior colliculus
  • interaural time differences
  • population coding
  • spatial hearing

Introduction

The ability to accurately localize sounds is critical for directing behavior, as well as for identifying and segregating individual sources within complex acoustic scenes (Cherry, 1953; Bronkhorst, 2000; Darwin, 2008). The dominant cue for the localization of a low-frequency sound such as speech is the difference in its arrival time at the two ears, referred to as the interaural time difference (ITD; Wightman and Kistler, 1992). ITD sensitivity in the mammalian brain arises in the medial superior olive (MSO) in the auditory brainstem where cells are sensitive to microsecond differences in the arrival time of inputs from the two ears. The spike rates of cells in the MSO and subsequent subcortical processing stages are modulated by ITD, with most cells responding preferentially to sounds with ITDs corresponding to locations in the contralateral hemifield (for review, see Grothe et al., 2010).

While ITD processing in subcortical areas has been extensively studied, there have been relatively few studies of ITD processing in auditory cortex. It is clear that auditory cortex is necessary for ITD processing in both animals and humans, though lesions in either hemisphere cause a contralateral deficit in spatial processing in animals (Jenkins and Masterton, 1982; Jenkins and Merzenich, 1984; Malhotra et al., 2004), while right auditory cortex appears both necessary and sufficient for ITD processing in humans (Yamada et al., 1996; Tanaka et al., 1999). ITD tuning in primary auditory cortex (A1) was first reported several decades ago (Brugge et al., 1969; Brugge and Merzenich, 1973), but the few studies in A1 with large samples that have been performed since have produced inconsistent results: a study in cats reported results similar to those in subcortical areas, with nearly all cells responding preferentially to ITDs corresponding to locations in the contralateral hemifield (Reale and Brugge, 1990), while studies in chinchillas, rabbits, and monkeys reported a weaker contralateral bias with preferred ITDs distributed more evenly across the physiological range (Benson and Teas, 1976; Fitzpatrick et al., 2000; Scott et al., 2009). There have been no direct studies of single-cell ITD sensitivity in human cortex, but recent EEG and MEG studies suggest a strong contralateral bias (Magezi and Krumbholz, 2010; Salminen et al., 2010).

In this study, we characterize the neural representation of ITD in A1 of gerbils, one of the most widely used model species for studies of ITD processing. In gerbils, the vast majority of cells in subcortical structures have preferred ITDs corresponding to locations in the contralateral hemifield (Spitzer and Semple, 1995; Siveke et al., 2006; Pecka et al., 2008; Lesica et al., 2010), consistent with a two-channel representation in which the ITD of a sound is encoded by the difference in the overall activity of the two brain hemispheres (McAlpine et al., 2001). Here we show that the neural representation of ITDs is transformed between inferior colliculus (IC) and A1 such that the preferred ITDs of A1 cells are distributed evenly throughout the physiological range without any contralateral bias. We examine the impact of this transformation on the population coding of ITDs and assess the ability of two-channel and labeled-line codes to account for gerbil behavioral acuity.

Materials and Methods

In vivo recordings.

All procedures were approved under the UK Animals (Scientific Procedures) Act of 1986. Nineteen adult male gerbils (70–90 g, P60–P120) were anesthetized for surgery with an initial injection of a mix of fentanyl, medetomidine, and midazolam or ketamine and xylazine, and the same solution was infused continuously during recording. A small metal rod was mounted on the skull and used to secure the head of the animal in a stereotaxic device in a sound-attenuated chamber. A craniotomy was made over the inferior colliculus or the primary auditory cortex, an incision was made in the dura mater, and a multi-tetrode array (NeuroNexus) was inserted into the brain. The array had four shanks spaced 0.2 mm apart, and each shank had two tetrodes spaced 0.15 mm apart. Recordings were made with a sampling rate of 25 kHz. Only recordings from the central nucleus of the IC and A1 were analyzed. Because the array covered a large area, recording sites in the central nucleus of the IC could be distinguished from those in other areas by comparison of their responses to tones (Aitkin et al., 1975; Syka et al., 2000), and A1 could be distinguished from other fields based on the direction of the tonotopic gradient (Thomas et al., 1993). A1 recordings were made between 1 and 1.5 mm below the cortical surface (most likely layer V; Happel et al. (2010)). We choose to record in layer V because we found the single-unit yield to be higher there than in layer IV in pilot experiments (we did not try other layers). Though it is difficult to say exactly why this would be the case, the fact that, relative to layer IV, layer V cells are large and sparsely packed, and spike with lower rates and less synchronously, may allow for single units to be more easily separated from multi-unit background. In both IC and A1, recordings were targeted to areas with low preferred frequencies.

Spike sorting.

The procedure for the isolation of single-unit spikes consisted of (1) bandpass filtering each channel and the tetrode array between 500 and 5000 Hz; (2) whitening each tetrode, i.e., projecting the signals from the four channels into a space in which they are uncorrelated; (3) identifying potential spikes as snippets with energy (Choi et al., 2006) that exceeded a threshold (with a minimum of 0.7 ms between potential spikes); (4) projecting each of the snippets into the space defined by the first three principal components for each channel; (5) identifying clusters of snippets within this space using KlustaKwik (http://klustakwik.sourceforge.net) and Klusters (Hazan et al., 2006); and (6) quantifying the likelihood that each cluster represented a single unit using isolation distance (Schmitzer-Torbert et al., 2005). Isolation distance assumes that each cluster forms a multidimensional Gaussian cloud in feature space and measures, in terms of the SD of the original cluster, the increase in the size of the cluster required to double the number of snippets within it. The number of snippets in the “noise” cluster (multi-unit activity) for each tetrode was always at least as large as the number of spikes in any single-unit cluster. Only clusters with an isolation distance >20 were classified as single units and included in our analysis.

Sound delivery.

Sounds were generated with a 48 kHz sampling rate, attenuated, and delivered to speakers. Speakers (Etymotic ER2) coupled to tubes were inserted into both ear canals for sound presentation along with microphones for calibration. The frequency response of these speakers measured at the entrance of the ear canal was flat (±5 dB) between 0.2 and 5 kHz. At each recording site, a sequence of tones with different frequencies and intensities with 5 ms cosine on and off ramps were presented to characterize frequency tuning. Speech and broadband noise were then presented at 60 dB SPL with nine different ITDs spanning the physiological range for gerbils (±160 μs in 40 μs steps) to characterize ITD tuning (with positive values of ITD denoting sounds leading at the ear contralateral to the recording site). These sounds were 500 ms in duration and were presented 32 or 64 times each in random order with a 500 ms pause between sounds and 2 ms cosine on and off ramps. Two different tokens of speech were used. Token 1 was presented to all cells in IC and A1 (n = 188 and 906, respectively). Token 2 was presented to all cells recorded in the left A1 under fentanyl, medetomidine, and midazolam (n = 517). Broadband noise was presented to a subset of cells recorded in the left A1 (n = 492) and all cells recorded in the right A1 (n = 100) and all IC cells.

Decoding ITD from spike rates.

To decode responses based on spike rate alone, we used maximum likelihood decoding. The probability that a spike rate r was evoked by an ITD s is given by Bayes' rule as p(s|r) = (p(r|s)p(s))/p(r). Because all ITDs were presented with equal probability, p(s|r) ∝ p(r|s). Thus, the ITD that is most likely to have caused a given response is simply arg maxs p(r|s). We assumed that the distribution of spike rates evoked by a given ITD was Gaussian (with truncation at zero if necessary). This assumption improved performance in cross-validated testing. We did not place any constraints on the mean spike rates at each ITD, i.e., the shape of the ITD tuning curve. The significance of ITD tuning was assessed by decoding responses after randomizing the pairing of responses and ITDs (Monte Carlo resampling). ITD tuning was considered significant if decoding performance was >4 SDs above the mean performance for 100 different sets of shuffled responses.

To decode population responses with a labeled-line decoder, the joint probability of a set of spike rates from N cells, p(r1,r2,…,rN|s), was computed as the product of the probabilities of the spike rate of each cell, ∏n=1Np(rn|s). Conditional independence was assumed because noise correlations were extremely weak (see Results). To decode population responses with a two-channel decoder, populations were split into two groups, and the responses for the second group were flipped with respect to ITD (as the majority of responses were recorded in the left A1 and our sample of responses from the right A1 was relatively small, we used only the responses from the left A1 for population decoding). The joint probability of the set of total spike rates in each hemisphere, p(rleft,rright|s), was then computed as the product of the probabilities of the total spike rate in each hemisphere, p(rleft|s) p(rright|s). For each population size tested, 100 different random subpopulations were drawn, with replacement, from the full sample of cells (bootstrap resampling). For the comparison with behavior, decoding performance was measured for responses to some ITDs that were not actually presented during the experiment. To measure decoding performance for these ITDs, we simulated responses (spike counts) by sampling from a Gaussian distribution with mean and variance taken as a linear interpolation between the values for the closest ITDs that were actually presented (i.e., it was assumed that the mean and variance of the spike rates for each cell were smooth functions of ITD). This approach was chosen to obtain a robust estimate of decoding performance for smaller ΔITDs than were actually presented. An alternative approach would have been to extrapolate to smaller ΔITDs from the performance at ΔITDs that were actually presented, but this would require an assumption about the shape of the performance curve for small ΔITDs, which is difficult to predict directly from our data.

Decoding ITD from spike times.

To decode responses based on spike timing, we used the metric introduced by Victor and Purpura (1996), which measures the distance between two spike trains as the overall cost of the set of operations required to transform one spike train into the other, with possible operations including the insertion of a spike, the deletion of a spike, and the time shift of a spike (Goldberg et al., 2009). By changing the cost of time shifting a spike relative to deleting the spike at one time and inserting it at another, the metric can be used to evaluate the distance between spike trains at different timescales. Decoding using this metric was performed as follows. (1) A single spike train was removed from the full set of all spike trains. (2) The distance between the removed spike train and each of the remaining spike trains in the set was computed across a range of timescales spaced logarithmically between 1 ms and 1 s. (3) For each timescale, the removed spike train was assigned to the sound for which its average distance to the remaining spike trains evoked by that sound was smallest. This process was repeated for all spike trains in the set to obtain a percentage correct for each timescale, and the overall percentage correct was taken as the maximum value across timescales. Decoding based on spike timing was considered significantly better than decoding based on spike rate alone if decoding performance based on spike timing was >4 SDs above the mean performance for decoding based on spike rate computed via bootstrap resampling.

Results

We made multi-tetrode recordings (Fig. 1A) from populations of single units in the IC and A1 of anesthetized gerbils. Our methods for IC recordings have been described in detail previously (Garcia-Lazaro et al., 2013). For A1 recordings (Fig. 1B), we aligned the shanks of the tetrode array along the rostrocaudal axis (approximately parallel to the tonotopic gradient in A1) and recorded from depths between 1 and 1.5 mm below the cortical surface (most likely layer V; Happel et al. (2010)). The direction of the tonotopic gradient in A1 was evident in the multi-unit activity across tetrodes, as illustrated by the frequency response areas (FRAs) for an example recording site shown in Figure 1C. We used a semi-automated clustering algorithm (see Materials and Methods) to isolate single units based on the first three principal components of their spike waveforms across each of the four channels of a tetrode. Clusters corresponding to single units (colors) and multi-unit noise (gray) are shown for an example tetrode in Figure 1D (note that is a 2D projection of a 12D space). We quantified the quality of each cluster based on its isolation distance (Schmitzer-Torbert et al., 2005) and set a threshold value of 20 for a cluster to be classified as a single unit (this value corresponded to detection of ∼90% of spikes from a target neuron with a false alarm rate of ∼1% in paired intracellular and tetrode recordings in hippocampus). The spike waveforms for two example single units with relatively low (24.5) and high (69.9) isolation distances are shown in Figure 1E. The median isolation distance across our sample of single units in A1 was 32.5 (Fig. 1F).

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

Multi-tetrode recordings in gerbil A1. A, A schematic diagram of the electrode arrangement on the multi-tetrode array. Thirty-two electrodes were grouped into eight tetrodes. B, A schematic diagram of the gerbil auditory cortex illustrating the alignment of the tetrode array with respect to the tonotopic gradient in A1 (modified from Thomas et al. (1993)). C, The FRAs for the multi-unit activity on each tetrode from a typical recording site. The first column shows the FRAs for the tetrodes on the most rostral shank, while the last column shows the FRAs for the tetrodes on the most caudal shank. Multi-unit activity was summed across the four electrodes on each tetrode. For each tetrode, the center frequency (CF) estimated from the FRA is indicated. D, Spikes from single units were identified by projecting spike waveforms into principal component space (12 dimensions corresponding to 3 principal components for each electrode). An example 2D projection that illustrates the isolation of different single-unit clusters is shown, along with the isolation distance of each cluster. Single-unit clusters are shown in color, and undifferentiated multi-unit noise is shown in gray. E, Spike waveforms for two single units (overlaid on a sample of multi-unit noise waveforms). F, A histogram of the isolation distances for all of the single units in our A1 sample.

We targeted our recordings to areas with low preferred frequencies. The distributions of best frequencies (BFs) for our samples of IC and A1 cells are shown in Figure 2A. As our main goal was to compare the neural representation of ITDs with existing measures of gerbil behavioral acuity in the localization of a single broadband low-frequency sound source (Lesica et al., 2010), we restricted our analysis to responses to broadband sounds (speech and noise) with ITDs spanning only the physiological range (±160 μs in 40 μs steps; with positive values of ITD denoting sounds leading at the ear contralateral to the recording site). Note that we have chosen to use a range of ITDs that is slightly larger than that measured for gerbils by Maki and Furukawa (2005), as their measurements were made for frequencies above 1.5 kHz and the physiological range of ITDs tends to increase for lower frequencies (Rébillat et al., 2014).

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Responses to speech at different ITDs in gerbil IC and A1. A, The distribution of BFs in our samples of IC and A1 cells. B, The responses of example cells with significant ITD tuning in IC and A1. Each column shows the FRA for one cell, along with raster plots for the responses to speech at five different ITDs spanning the physiological range and the tuning curve showing the mean spike rate in response to speech as a function of ITD. The black line and gray bands on the tuning curve plots indicate the mean ± 1 SD.

The responses to speech with different ITDs from example cells with significant ITD tuning are shown in Figure 2B. Each column shows the FRA for one cell, along with raster plots for the responses to speech at five different ITDs spanning the physiological range and the tuning curve showing the mean spike rate as a function of ITD. To assess the strength of each cell's ITD tuning, we used a decoder to measure the accuracy with which the spike rate on a single trial could be used to infer which of nine possible ITDs evoked it. We considered ITD tuning to be significant if decoding accuracy was more than 4 SDs above the mean performance for shuffled responses.

The fraction of cells with significant ITD tuning for speech was higher in IC than in A1 (IC: 117/188 cells, 62%; A1: 239/517, 46%). This difference was not due to differences in the distribution of BFs in the two populations; the fraction of cells with significant tuning in our entire A1 sample was the same as that of random subsamples of A1 cells with BFs matched to those of our sample of IC cells (45 ± 3%). The fraction of cells with significant ITD tuning for broadband noise was also higher in IC than in A1 (IC: 134/188 cells, 71%; A1: 133/203, 65%). For all subsequent analyses in this study, only cells with significant ITD tuning were included.

Best ITDs in A1 are distributed evenly across the physiological range

To compare the representation of ITDs in IC and A1, we began by measuring the spike rate tuning curve for each cell's response to speech at different ITDs. All of the example IC cells shown in Figure 2B responded most strongly to the ITD corresponding to the contralateral edge of the physiological range (+160 μs), while each example A1 cell had a different preferred ITD. These examples were representative of the IC and A1 in general; each row in the images in Figure 3A shows the ITD tuning curve for one cell (with cells sorted by best ITD), and the histograms in Figure 3B show the distributions of best ITDs across all cells.

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Best ITDs in A1 are distributed evenly across the physiological range. A, The ITD tuning curves for speech for all significantly tuned cells in our samples of IC and A1 cells. Each row shows the ITD tuning curve for one cell. All tuning curves were normalized to have the same maximum and minimum for plotting. Cells were sorted by best ITD for plotting. B, The histograms of the best ITDs for speech for all significantly tuned cells in our samples. C, The same A1 tuning curves shown in A with cells sorted either by BF or decoding performance. Decoding performance was measured as the percentage of single trial responses that were assigned to the correct ITD by a spike rate decoder (the chance level was 1/9).

The majority of IC cells in our sample had a best ITD of +160 μs, while the best ITDs in A1 were evenly distributed across the physiological range. To quantify the degree to which the distribution of best ITDs in each area was biased toward ipsilateral or contralateral values, we measured the percentage of cells with best ITDs in the contralateral hemifield (cells with best ITD = 0 were ignored). The distribution of best ITDs in the IC was strongly biased toward the contralateral side, with 83% of cells having best ITDs in the contralateral hemifield. In contrast, the distribution of best ITDs in A1 was unbiased, with only 53% of cells having best ITDs in the contralateral hemifield.

The best ITDs of A1 cells were unrelated to their BFs, as illustrated in Figure 3C, left, which shows the same ITD tuning curves for all A1 cells as in Figure 3A, but with cells sorted by BF. For any particular BF, there were cells with a range of different best ITDs, and, across the entire population, best ITD and BF were uncorrelated (r = 0.03, p = 0.61). Figure 3C, right, also shows the same tuning curves, but with the cells sorted by ITD decoding performance. For any particular level of performance, there were cells with a range of different best ITDs, though there was a weak, but significant, correlation between best ITD and decoding performance across the entire population (r = 0.18, p = 0.003), indicating that ITD tuning was slightly stronger for cells with best ITDs corresponding to contralateral locations than for cells with best ITDs corresponding to ipsilateral locations.

ITD tuning is consistent across different sounds

We next investigated whether ITD tuning was consistent across sounds with different spectrotemporal properties by comparing ITD tuning curves for speech and broadband noise in both IC and A1, as well as for two different speech tokens in A1. Figure 4A shows the raster plots for the responses of an example A1 cell to the different sounds at five different ITDs, along with the tuning curves showing the mean spike rate as a function of ITD. The ITD tuning for this example cell was consistent across all three sounds, with the strongest responses evoked by ITDs near 0, corresponding to locations near the midline.

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

ITD tuning is consistent across different sounds. A, The responses of an example cell from A1 with significant ITD tuning for two different segments of speech and broadband noise. Each column shows the raster plots for the responses to one sound at five different ITDs spanning the physiological range and the tuning curve showing the mean spike rate as a function of ITD. The black line and gray bands on the tuning curve plots indicate the mean ± 1 SD. B, The histograms of the correlation coefficients between each cell's ITD tuning curves for speech segment 1 and broadband noise in the IC and A1, and the two speech segments in A1. Only cells with significant ITD tuning for both of the sounds being compared were included. The median value across each sample of cells is noted on each histogram.

To quantify the similarity of ITD tuning across sounds for each cell, we measured the correlation coefficient between ITD tuning curves. As shown in Figure 4B, ITD tuning curves were highly similar across sounds for nearly all cells in both IC and A1; the median correlation between ITD tuning curves for speech and noise was 0.97 in IC and 0.87 in A1, and the median correlation between ITD tuning curves for two different segments of speech was 0.89 in A1. This suggests that the transformation of the representation of ITDs from IC to A1 is a general phenomenon that will be evident for any complex sound.

Spike timing carries relatively little information about ITDs

Studies of IC and A1 responses have shown that for the coding of spectral notches and interaural level differences in high-frequency sounds, spike timing contains substantial information beyond that in spike rate alone (Furukawa and Middlebrooks, 2002; Chase and Young, 2006), but the role of spike timing in coding ITDs in low-frequency sounds is not yet clear. To assess the role of spike timing in coding ITDs, we compared the performance of a decoder based on spike rate alone with that of a decoder that used a distance metric to consider the full spike train at the optimal timescale for each cell (Victor and Purpura, 1996).

Figure 5A shows the raster plots for the responses of an example A1 cell to speech at three different ITDs. For this cell, the timing of some spiking events varied with ITD (see arrows), and considering spike timing in addition to spike rate resulted in a 50% improvement in decoding performance (21% correct for the timing decoder and 14% correct for the rate decoder, for nine possible ITDs). This cell was, however, not typical of either IC or A1; as shown in Figure 5B, the improvement in decoding in both IC and A1 that resulted from considering spike timing in addition to spike rate was relatively small for both speech and noise. In IC, the improvement in the performance of the timing decoder over the rate decoder was significant for 55% of cells for speech and 60% of cells for noise, but the median improvement for those cells with significant improvement was only 11% for speech and 10% for noise. The improvement in A1 was higher than in IC for speech (50% of cells significant, median improvement 18%), and similar for noise (57% significant, median improvement 12%). These results suggest that spike timing is unlikely to play a major role in the coding of ITDs in either IC or A1.

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5.

Spike timing carries relatively little information about ITDs. A, The responses of an example cell from A1 with significant information about ITD in spike timing. The raster plots show the responses to speech at three different ITDs. Spiking events that are unique to a particular ITD are marked with arrows. B, Scatter plots showing the percentage improvement in decoding ITD that resulted from considering spike timing in addition to spike rate versus the percentage correct for spike rate alone for responses to speech and noise in IC and A1. The values for cells for which the improvement was statistically significant are shown as black circles. The percentage of cells for which the improvement was significant and the median improvement across those cells is noted on each histogram.

ITD tuning in A1 is qualitatively similar under different anesthesia

All of the responses described above were recorded under a mix of fentanyl, medetomidine, and midazolam (FMM). As the responses of neurons in gerbil A1 are known to vary with brain state (Ter-Mikaelian et al., 2007), we also made recordings under a mix of ketamine and xylazine (KX) to determine whether our observations of ITD tuning in A1 were dependent on our choice of anesthesia. In general, ITD tuning in A1 was much weaker under KX: only 39 of 289 cells (13%) had significant ITD tuning for speech (compared with 46% under FMM), and only 37 of 289 cells (12%) had significant ITD tuning for broadband noise (compared with 65% under FMM). However, as shown in Figure 6, the qualitative nature of ITD tuning under KX was similar to that under FMM: for those cells with significant ITD tuning for speech, best ITDs were distributed across the physiological range with no bias toward ITDs corresponding to locations in the contralateral hemifield (Figs. 6A,B), and spike timing carried relatively little information about ITD (Fig. 6C). All further analyses of ITD tuning in A1 described below were performed only on responses recorded under FMM.

Figure 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6.

ITD tuning in A1 is qualitatively similar under different anesthesia. A, The ITD tuning curves for speech for all significantly tuned cells in our sample of A1 cells recorded under ketamine and xylazine, plotted as in Figure 3A. B, The histogram of the best ITDs for speech for all significantly tuned cells in our sample of A1 cells recorded under ketamine and xylazine. C, Scatter plots showing the percentage improvement in decoding ITD that resulted from considering spike timing in addition to spike rate versus the percentage correct for spike rate alone for responses to speech and noise in A1 recorded under ketamine and xylazine, plotted as in Figure 5B.

ITD tuning in left and right A1 are similar

All of the responses described above were recorded from the left A1. To verify that ITD tuning in A1 was similar in both brain hemispheres, we made additional recordings from the right A1. ITD tuning in the right A1 was somewhat weaker than that in the left A1: 32 of 100 cells (32%) had significant ITD tuning for speech (compared with 46% in the left A1), and 46 of 100 cells (46%) had significant ITD tuning for broadband noise (compared with 65% in the left A1). As in the left A1, the best ITDs for cells in the right A1 were distributed across the physiological range with only a weak bias toward ITDs corresponding to locations in the contralateral hemifield (Figs. 7A,B), and spike timing carried relatively little information about ITD (Fig. 7C).

Figure 7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 7.

ITD tuning in left and right A1 are similar. A, The ITD tuning curves for speech for all significantly tuned cells in our sample of cells recorded in the right A1, plotted as in Figure 3A. B, The histogram of the best ITDs for speech for all significantly tuned cells in our sample of cells recorded in the right A1. C, Scatter plots showing the percentage improvement in decoding ITD that resulted from considering spike timing in addition to spike rate versus the percentage correct for spike rate alone for responses to speech and noise in the right A1, plotted as in Figure 5B.

Two-channel decoding of population responses in A1 results in a loss of information

The difference in the distributions of best ITDs in IC and A1 suggest a fundamental difference in the coding of ITDs at the population level. We considered two different population codes for ITD: a “two-channel” code (McAlpine et al., 2001; Lüling et al., 2011; Day and Delgutte, 2013) that considers only the total spike rate in each brain hemisphere [also known as a summed code (Lesica et al., 2010) or hemispheric code (Goodman et al., 2013)], and a “labeled-line” code that considers the tuning of individual cells [also known as a distributed code (Lesica et al., 2010) or pattern code (Day and Delgutte, 2013; Goodman et al., 2013].

We have shown previously that because the ITD tuning curves of most cells in gerbil IC are similar, a two-channel decoder performs almost as well as a labeled-line decoder at inferring the ITD of the sound that evoked a particular single trial population response (Lesica et al., 2010). However, for a population with more heterogeneous tuning curves, considering only the total spike rate in each hemisphere can impair decoding performance (Day and Delgutte, 2013; Goodman et al., 2013). We compared the performance of labeled-line and two-channel decoders on IC and A1 responses to speech and noise for populations of increasing size. Rather than constrain decoding to a particular computation (e.g., the difference in total spike rate between the two hemispheres), we used a maximum likelihood approach to infer which of nine possible ITDs evoked each single trial population response based on the joint distributions of spike rates in each hemisphere (two-channel) or in individual cells (labeled line; Miller and Recanzone, 2009; Day and Delgutte, 2013). Because the noise correlations between pairs of simultaneously recorded cells in both IC and A1 were extremely weak (IC: 0.006 ± 0.018, n = 9112; A1: 0.004 ± 0.012, n = 3606), we assumed that the spike rates of individual cells were conditionally independent (Garcia-Lazaro et al., 2013).

As shown in Figure 8, A and B, the labeled-line decoder performed well on both IC and A1 responses to speech and noise, with no loss of information between the two areas (IC: median performance of 94% correct for speech for the largest populations, 93% for noise; A1: 93% for speech, 91% for noise). In contrast, the performance of the two-channel decoder on A1 responses was much worse than its performance on IC responses (IC: 74% for speech, 73% for noise; A1: 49% for speech, 54% for noise). Thus, for a labeled-line code based on the spike rates of individual cells, the information about ITDs present in IC is preserved in A1, but for a two-channel population code based on the total spike rate in each hemisphere, there is a substantial loss of information about ITDs between IC and A1.

Figure 8.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 8.

Two-channel and labeled-line decoding of ITD from population responses. A, The performance of labeled-line and two-channel decoders on IC and A1 responses to speech for populations of increasing size. Performance was measured as the percentage of single trial responses that were assigned to the correct ITD by the decoder. The chance level (1/9) is indicated. The black line and colored bands indicate the mean ± 2 SDs of the performance for 100 different random subpopulations of each size drawn from the full sample of cells. B, The performance of labeled-line and two-channel decoders on IC and A1 responses to noise, plotted as in A.

It should be noted that although the performance of the two-channel decoder on A1 responses was relatively poor, it was still well above chance; although the ITD tuning curve peaks for A1 cells are evenly distributed throughout the physiological range, there is still a significant monotonic modulation of the total population spike rate with ITD, though this modulation is much weaker than that in the IC (Fig. 9A). It should also be noted that the labeled-line decoder significantly outperformed the two-channel decoder not only on A1 responses, but also on IC responses. This result is consistent with recent studies, suggesting that the heterogeneity of tuning curves in IC can carry significant information about ITD (Day and Delgutte, 2013; Goodman et al., 2013).

Figure 9.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 9.

Comparing decoding of population responses and behavior. A, The tuning curve showing the total mean spike rate of our populations of cells in IC and A1 in response to speech and noise as a function of ITD. The black lines and colored bands on the tuning curve plots indicate the mean ± 1 SD. B, A comparison of the performance of labeled-line and two-channel decoders to gerbil behavior. The black dots show the actual performance of gerbils in lateralizing low-frequency noise bursts as a function of the difference in ITD approximated from the angle of separation between two speakers centered on the midline (median performance for 7 gerbils from Lesica et al., 2010). The colored lines show the performance of labeled-line and two-channel decoders on IC and A1 responses to pairs of noise bursts centered on the midline. The dots on the colored lines correspond to ΔITDs that were actually tested experimentally, the remainder of the values was obtained after interpolating ITD tuning curves as described in Materials and Methods. Performance for ΔITD = X μs was assessed by decoding responses to noise with ITD = ±(X/2) μs. C, The colored lines show the performance of labeled-line and two-channel decoders on IC and A1 responses to pairs of noise bursts centered on different ITDs. Performance for center ITD = X μs was assessed by decoding responses to noise with ITD = X ± 10 μs for the IC and ITD = X ± 20 μs for A1. The black line and colored bands indicate the mean ± 2 SDs of the performance for 100 different bootstrap samples of cells from the full populations.

Both two-channel and labeled-line decoding of population responses are sufficient to explain behavior

Our decoding results demonstrate that a labeled-line code carries substantially more information about ITDs than a two-channel code in gerbil A1. However, there is no guarantee that the code that is most informative about a particular sound feature is the one that underlies its perception, especially in the cortex where, presumably, the same neural circuitry is used to analyze many different features (Brette, 2010). One approach to rule in or rule out different candidate codes for a particular feature is to determine whether they are sufficient to account for behavioral performance (Jacobs et al., 2009).

There have only been a few behavioral studies of sound localization in gerbils (Heffner and Heffner, 1988; Maier and Klump, 2006; Maier et al., 2008; Lesica et al., 2010; Carney et al., 2011; Lingner et al., 2012). Figure 9B shows the accuracy with which gerbils lateralized low-frequency noise bursts as a function of the difference in ITD approximated from the angle of separation between two speakers across the midline (median performance for seven gerbils from Lesica et al. (2010)). We used the labeled-line and two-channel decoders described above to simulate the same behavioral task and infer which of two possible ITDs centered on the midline evoked each single trial population response to noise. Surprisingly, we found that both labeled-line and two-channel decoding were sufficient to reproduce this performance, even for A1 responses (Fig. 9B).

Though ability of gerbils to use ITDs to localize sounds has only been tested for pairs of sounds centered on the midline, it is known that behavioral acuity tends to decrease for sounds centered on more lateral locations in many mammals, for example, humans (Mossop and Culling, 1998) and rabbits (Ebert et al., 2008). We examined the ability of the labeled-line and two-channel decoders to infer which of two possible ITDs evoked each single trial population response as a function of the ITD on which the sounds were centered (different ΔITDs were used for IC and A1 so that the performance of the labeled-line decoder for pairs of sounds centered on the midline was ∼90% correct for both brain areas). As shown in Figure 9C, while decoder performance on IC responses decreased for more lateral sounds as expected, decoder performance on A1 responses was relatively consistent across the physiological range.

Discussion

We have shown that the neural representation of ITDs in gerbils is transformed from IC to A1. In the IC, we found that most cells responded maximally to ITDs corresponding to the contralateral edge of the physiological range, consistent with previous studies of ITD processing in different subcortical stages in gerbils (Spitzer and Semple, 1995; Siveke et al., 2006; Pecka et al., 2008; Lesica et al., 2010). In contrast, the preferred ITDs of A1 cells were distributed evenly throughout the physiological range, with an equal number of cells preferring ITDs corresponding to ipsilateral and contralateral locations. This transformation in the distribution of preferred ITDs resulted in a loss of information between IC and A1 when using a two-channel decoder that considered only the total spike rate in each brain hemisphere, but not when using a labeled-line decoder that considered the tuning of individual cells. However, despite this loss of information, the two-channel decoder was still sufficient to reproduce gerbil behavioral performance.

Our analysis has revealed several aspects of the neural representation of ITDs in A1 that appear inconsistent with existing behavioral data. First, decoding of ITD from A1 activity results in performance that is far better than that observed behaviorally. Behavioral performance may be expected to be worse than decoder performance for activity from the early stages of peripheral processing, but differences of this degree in cortex are more surprising. Second, both left and right A1 appear to have a complete representation of azimuthal space (i.e., best ITDs span the full physiological range). Thus, it is unclear why a lesion to either the left or right A1 would cause a deficit in the localization of only contralateral sounds, as is the case in several mammals (Jenkins and Masterton, 1982; Jenkins and Merzenich, 1984; Kavanagh and Kelly, 1987; Malhotra et al., 2004). Third, while decoder performance on IC responses decreased for sounds with ITDs corresponding to more lateral locations, consistent with behavioral observations in several mammals (Mossop and Culling, 1998; Ebert et al., 2008), decoder performance on A1 responses was similar for ITDs corresponding to medial and lateral locations.

One possible explanation for the apparent mismatch between the neural representation of ITDs in A1 and the existing behavioral data in gerbils is that A1 is not actually required for or involved in the localization of single sound sources in a quiet background. While A1 seems to play a role in localization in most mammals that have been tested (King and Middlebrooks, 2011), it does not appear necessary for sound localization in rats (Kelly and Kavanagh, 1986). It may also be that if localization was tested under more difficult (e.g., reverberant) conditions, a better match between A1 activity and behavioral performance would become apparent. Another possible explanation for the disconnect between the neural representation of ITDs in A1 and behavior arises when one considers that the role of cortex is presumably to combine information about different acoustic features for the analysis of complex auditory scenes. While a given population of subcortical cells can be specialized for the processing of a particular feature, cortical populations may need to process multiple stimulus features simultaneously. Thus, the representation of ITDs in A1 (and/or the manner in which information from A1 is decoded in higher cortical areas, which may differ from the decoders we tested) may not be specialized for sound localization per se, but rather for the general processing of complex scenes, allowing, for example, sound sources from different locations to be processed by different subpopulations of cells (Middlebrooks and Bremen, 2013) and facilitating the allocation of attentional resources to enhance or suppress activity related to a given source (Lee and Middlebrooks, 2011).

How does the transformation of ITD tuning between IC and A1 in gerbils compare with that in other species?

The change in the distribution of ITD tuning curve peaks from strongly contralaterally biased in IC to unbiased in A1 makes gerbils unique among species for which ITD tuning in midbrain and cortex has been systematically studied (for review, see Vonderschen and Wagner, 2014). There is a clear transformation of the representation of ITDs in the auditory pathway of barn owls, but in the opposite direction, with narrowly tuned midbrain inputs converging to form broadly tuned channels in the forebrain (Vonderschen and Wagner, 2009, 2012). The differences in ITD tuning between IC and A1 in other mammals are not as clear as those in gerbils, and results differ across species. Studies in cats have reported a relatively strong contralateral bias in both IC and A1 (Reale and Brugge, 1990; Yin and Chan, 1990), while in rabbits, best ITDs span the physiological range with a relatively weak contralateral bias throughout the entire auditory pathway (Fitzpatrick et al., 2000), though tuning curves get sharper in more central stations (Fitzpatrick et al., 1997). Our results are perhaps most similar to those from chinchillas and primates; in chinchillas, there appears to be a strong contralateral bias in the IC (Bremen and Joris, 2013), but only a relatively weak contralateral bias in A1 (Benson and Teas, 1976). While there have been no systematic studies of ITD tuning in the primate midbrain, studies of spatial sensitivity in the IC suggest a strong contralateral bias (Groh et al., 2001, 2003; Zwiers et al., 2004), while ITD tuning in A1 exhibits a relatively weak contralateral bias (Scott et al., 2009).

Our results may also have implications for the study of ITD processing in humans. EEG and MEG studies in humans investigated the representation of ITDs based on measurements of the change in overall cortical activity in response to a change in ITD (Magezi and Krumbholz, 2010; Salminen et al., 2010). In a labeled-line representation, a change in ITD in either direction should cause an increase in overall activity (as the sensory drive is directed toward an unadapted neuronal subpopulation), while in a two-channel representation, overall activity within a given hemisphere should increase with change in ITD in one direction, and decrease with a change in ITD in the other direction. Both studies found that the sign of the change in overall activity depended on the direction of the change in ITD, and, thus, argued for a two-channel representation. However, our A1 data demonstrate that a relatively coarse two-channel representation can coexist with a much more sensitive labeled-line representation and suggest that more detailed studies may be required to determine the true nature of the cortical representation of ITDs in humans.

What neural mechanisms underlie the transformation between IC and A1?

The neural circuitry that facilitates the transformation of the neural representation of ITDs between IC and A1 is not yet clear. It is possible to transform a distribution of best ITDs with a strong contralateral bias into an unbiased one either through the addition of inputs with opposing preferences or through the subtraction of inputs with similar preferences but different tuning curves (Groh et al., 2003). In principle, either of these possibilities could be implemented between IC and A1, even within a single brain hemisphere, using either the small subpopulation of cells in each IC with best ITDs corresponding to locations in the ipsilateral hemifield, or the heterogeneity of tuning curves in the majority of cells with best ITDs corresponding to locations in the contralateral hemifield.

If the transformation between IC and A1 does involve integration across the two brain hemispheres, it is likely through callosal connections (Budinger et al., 2000), as the projections from the IC to the auditory thalamus and from the thalamus to A1 are predominantly ipsilateral (Winer and Schreiner, 2005, 2011). The possibility that callosal connections play a role in shaping ITD tuning in A1 could explain why a unilateral cortical lesion results in a behavioral deficit for only contralateral locations (Jenkins and Masterton, 1982; Jenkins and Merzenich, 1984; Malhotra et al., 2004); without callosal inputs, the residual sensitivity in the remaining A1 may be for ITDs corresponding to ipsilateral locations only. Efforts to identify the neural circuitry that underlie the transformation from IC to A1 should begin by determining the highest stage at which the distribution of best ITDs still has a strong contralateral bias. While we cannot be certain of the layer in which our A1 recordings were made, it was most likely layer V; it is possible that the distribution of best ITDs in layer IV still has a strong contralateral bias (indeed, recent studies have suggested that the responses of cells in layer IV of A1 are simply amplified versions of their thalamic inputs; Li et al., 2013), and that the transformation takes place between layer IV and layer II/III, or between layer II/III and layer V.

How does the ITD tuning in gerbil IC observed in this study compare with that observed previously?

The results of our population decoding analysis of IC responses differ somewhat from those of a similar analysis that we performed in a previous study (Lesica et al., 2010). In the previous study, the performance of the labeled-line and two-channel decoders was nearly identical, whereas in the current study, the performance of the labeled-line decoder was substantially better than that of the two-channel decoder. The difference in the performance of the two-channel decoder in the two studies is consistent with the differences in the distributions of best ITDs in the two populations of cells that were studied. In the previous study, we found that nearly all IC cells had best ITDs corresponding to the contralateral edge of the physiological range and, because of this homogeneity, very little information was lost when ignoring the tuning of individual cells and decoding only the total activity in the population. In the current study, however, 17% of the cells in our IC sample had best ITDs corresponding to locations in the ipsilateral hemifield, and this heterogeneity affected the performance of the two-channel decoder.

We believe that the difference in the distributions of best ITDs in the two studies is due to a difference in the fraction of the IC that was sampled during our recordings. In the previous study, we used a bundle of concentrically arranged electrodes that spanned a relatively small area and targeted the recordings to the rostromedial quadrant of the IC where the dominant input is provided by the MSO (Cant and Benson, 2006), while in the current study, we used a much larger electrode array and sampled a larger fraction of the IC. Thus, the results of the present study are likely a more accurate reflection of the processing of ITDs in the IC as a whole, and are consistent with other recent studies suggesting that differences in the ITD tuning of IC cells carry significant information (Day and Delgutte, 2013; Goodman et al., 2013).

Footnotes

  • This work was supported by the Wellcome Trust and the Munich Center for Neurosciences. We thank D. McAlpine, M. Pecka, and J. Bizley for helpful discussions.

  • The authors declare no competing financial interests.

  • This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.

  • Correspondence should be addressed to Nicholas Lesica, Ear Institute, University College London, Gower Street, London WC1E 6BT, UK. n.lesica{at}ucl.ac.uk

This article is freely available online through the J Neurosci Author Open Choice option.

References

  1. ↵
    1. Aitkin LM,
    2. Webster WR,
    3. Veale JL,
    4. Crosby DC
    (1975) Inferior colliculus. I. Comparison of response properties of neurons in central, pericentral, and external nuclei of adult cat. J Neurophysiol 38:1196–1207, pmid:1177012.
    OpenUrlAbstract/FREE Full Text
  2. ↵
    1. Benson DA,
    2. Teas DC
    (1976) Single unit study of binaural interaction in the auditory cortex of the chinchilla. Brain Res 103:313–338, doi:10.1016/0006-8993(76)90801-5, pmid:1252921.
    OpenUrlCrossRefPubMed
  3. ↵
    1. Bremen P,
    2. Joris PX
    (2013) Axonal recordings from medial superior olive neurons obtained from the lateral lemniscus of the chinchilla (Chinchilla laniger). J Neurosci 33:17506–17518, doi:10.1523/JNEUROSCI.1518-13.2013, pmid:24174683.
    OpenUrlAbstract/FREE Full Text
  4. ↵
    1. Brette R
    (2010) On the interpretation of sensitivity analyses of neural responses. J Acoust Soc Am 128:2965–2972, doi:10.1121/1.3488311, pmid:21110592.
    OpenUrlCrossRefPubMed
  5. ↵
    1. Bronkhorst A
    (2000) The cocktail party phenomenon: a review of research on speech intelligibility in multiple-talker conditions. Acust Acta Acust 86:117–128.
    OpenUrl
  6. ↵
    1. Brugge JF,
    2. Merzenich MM
    (1973) Responses of neurons in auditory cortex of the macaque monkey to monaural and binaural stimulation. J Neurophysiol 36:1138–1158, pmid:4761724.
    OpenUrlFREE Full Text
  7. ↵
    1. Brugge JF,
    2. Dubrovsky NA,
    3. Aitkin LM,
    4. Anderson DJ
    (1969) Sensitivity of single neurons in auditory cortex of cat to binaural tonal stimulation; effects of varying interaural time and intensity. J Neurophysiol 32:1005–1024, pmid:5347704.
    OpenUrlFREE Full Text
  8. ↵
    1. Budinger E,
    2. Heil P,
    3. Scheich H
    (2000) Functional organization of auditory cortex in the Mongolian gerbil (Meriones unguiculatus). IV. Connections with anatomically characterized subcortical structures. Eur J Neurosci 12:2452–2474, doi:10.1046/j.1460-9568.2000.00143.x, pmid:10947822.
    OpenUrlCrossRefPubMed
  9. ↵
    1. Cant NB,
    2. Benson CG
    (2006) Organization of the inferior colliculus of the gerbil (Meriones unguiculatus): differences in distribution of projections from the cochlear nuclei and the superior olivary complex. J Comp Neurol 495:511–528, doi:10.1002/cne.20888, pmid:16498677.
    OpenUrlCrossRefPubMed
  10. ↵
    1. Carney LH,
    2. Sarkar S,
    3. Abrams KS,
    4. Idrobo F
    (2011) Sound-localization ability of the Mongolian gerbil (Meriones unguiculatus) in a task with a simplified response map. Hear Res 275:89–95, doi:10.1016/j.heares.2010.12.006, pmid:21147208.
    OpenUrlCrossRefPubMed
  11. ↵
    1. Chase SM,
    2. Young ED
    (2006) Spike-timing codes enhance the representation of multiple simultaneous sound-localization cues in the inferior colliculus. J Neurosci 26:3889–3898, doi:10.1523/JNEUROSCI.4986-05.2006, pmid:16611804.
    OpenUrlAbstract/FREE Full Text
  12. ↵
    1. Cherry EC
    (1953) Some experiments on the recognition of speech, with one and with two ears. J Acoust Soc Am 5:975–979, doi:10.1121/1.1907229.
    OpenUrlCrossRef
  13. ↵
    1. Choi JH,
    2. Jung HK,
    3. Kim T
    (2006) A new action potential detector using the MTEO and its effects on spike sorting systems at low signal-to-noise ratios. IEEE Trans Biomed Eng 53:738–746, doi:10.1109/TBME.2006.870239, pmid:16602581.
    OpenUrlCrossRefPubMed
  14. ↵
    1. Darwin CJ
    (2008) Listening to speech in the presence of other sounds. Philos Trans R Soc Lond B Biol Sci 363:1011–1021, doi:10.1098/rstb.2007.2156, pmid:17827106.
    OpenUrlAbstract/FREE Full Text
  15. ↵
    1. Day ML,
    2. Delgutte B
    (2013) Decoding sound source location and separation using neural population activity patterns. J Neurosci 33:15837–15847, doi:10.1523/JNEUROSCI.2034-13.2013, pmid:24089491.
    OpenUrlAbstract/FREE Full Text
  16. ↵
    1. Ebert CS Jr.,
    2. Blanks DA,
    3. Patel MR,
    4. Coffey CS,
    5. Marshall AF,
    6. Fitzpatrick DC
    (2008) Behavioral sensitivity to interaural time differences in the rabbit. Hear Res 235:134–142, doi:10.1016/j.heares.2007.11.003, pmid:18093767.
    OpenUrlCrossRefPubMed
  17. ↵
    1. Fitzpatrick DC,
    2. Batra R,
    3. Stanford TR,
    4. Kuwada S
    (1997) A neuronal population code for sound localization. Nature 388:871–874, doi:10.1038/42246, pmid:9278047.
    OpenUrlCrossRefPubMed
  18. ↵
    1. Fitzpatrick DC,
    2. Kuwada S,
    3. Batra R
    (2000) Neural sensitivity to interaural time differences: beyond the Jeffress model. J Neurosci 20:1605–1615, pmid:10662850.
    OpenUrlAbstract/FREE Full Text
  19. ↵
    1. Furukawa S,
    2. Middlebrooks JC
    (2002) Cortical representation of auditory space: information-bearing features of spike patterns. J Neurophysiol 87:1749–1762, pmid:11929896.
    OpenUrlAbstract/FREE Full Text
  20. ↵
    1. Garcia-Lazaro JA,
    2. Belliveau LA,
    3. Lesica NA
    (2013) Independent population coding of speech with sub-millisecond precision. J Neurosci 33:19362–19372, doi:10.1523/JNEUROSCI.3711-13.2013, pmid:24305831.
    OpenUrlAbstract/FREE Full Text
  21. ↵
    1. Goldberg DH,
    2. Victor JD,
    3. Gardner EP,
    4. Gardner D
    (2009) Spike train analysis toolkit: enabling wider application of information-theoretic techniques to neurophysiology. Neuroinformatics 7:165–178, doi:10.1007/s12021-009-9049-y, pmid:19475519.
    OpenUrlCrossRefPubMed
  22. ↵
    1. Goodman DF,
    2. Benichoux V,
    3. Brette R
    (2013) Decoding neural responses to temporal cues for sound localization. Elife 2:e01312, doi:10.7554/eLife.01312, pmid:24302571.
    OpenUrlAbstract/FREE Full Text
  23. ↵
    1. Groh JM,
    2. Trause AS,
    3. Underhill AM,
    4. Clark KR,
    5. Inati S
    (2001) Eye position influences auditory responses in primate inferior colliculus. Neuron 29:509–518, doi:10.1016/S0896-6273(01)00222-7, pmid:11239439.
    OpenUrlCrossRefPubMed
  24. ↵
    1. Groh JM,
    2. Kelly KA,
    3. Underhill AM
    (2003) A monotonic code for sound azimuth in primate inferior colliculus. J Cogn Neurosci 15:1217–1231, doi:10.1162/089892903322598166, pmid:14709238.
    OpenUrlCrossRefPubMed
  25. ↵
    1. Grothe B,
    2. Pecka M,
    3. McAlpine D
    (2010) Mechanisms of sound localization in mammals. Physiol Rev 90:983–1012, doi:10.1152/physrev.00026.2009, pmid:20664077.
    OpenUrlAbstract/FREE Full Text
  26. ↵
    1. Happel MF,
    2. Jeschke M,
    3. Ohl FW
    (2010) Spectral integration in primary auditory cortex attributable to temporally precise convergence of thalamocortical and intracortical input. J Neurosci 30:11114–11127, doi:10.1523/JNEUROSCI.0689-10.2010, pmid:20720119.
    OpenUrlAbstract/FREE Full Text
  27. ↵
    1. Hazan L,
    2. Zugaro M,
    3. Buzsáki G
    (2006) Klusters, NeuroScope, NDManager: a free software suite for neurophysiological data processing and visualization. J Neurosci Methods 155:207–216, doi:10.1016/j.jneumeth.2006.01.017, pmid:16580733.
    OpenUrlCrossRefPubMed
  28. ↵
    1. Heffner RS,
    2. Heffner HE
    (1988) Sound localization and use of binaural cues by the gerbil (Meriones unguiculatus). Behav Neurosci 102:422–428, doi:10.1037/0735-7044.102.3.422, pmid:3395452.
    OpenUrlCrossRefPubMed
  29. ↵
    1. Jacobs AL,
    2. Fridman G,
    3. Douglas RM,
    4. Alam NM,
    5. Latham PE,
    6. Prusky GT,
    7. Nirenberg S
    (2009) Ruling out and ruling in neural codes. Proc Natl Acad Sci U S A 106:5936–5941, doi:10.1073/pnas.0900573106, pmid:19297621.
    OpenUrlAbstract/FREE Full Text
  30. ↵
    1. Jenkins WM,
    2. Masterton RB
    (1982) Sound localization: effects of unilateral lesions in central auditory system. J Neurophysiol 47:987–1016, pmid:7108581.
    OpenUrlFREE Full Text
  31. ↵
    1. Jenkins WM,
    2. Merzenich MM
    (1984) Role of cat primary auditory cortex for sound-localization behavior. J Neurophysiol 52:819–847, pmid:6512590.
    OpenUrlAbstract/FREE Full Text
  32. ↵
    1. Kavanagh GL,
    2. Kelly JB
    (1987) Contribution of auditory cortex to sound localization by the ferret (Mustela putorius). J Neurophysiol 57:1746–1766, pmid:3598629.
    OpenUrlAbstract/FREE Full Text
  33. ↵
    1. Kelly JB,
    2. Kavanagh GL
    (1986) Effects of auditory cortical lesions on pure-tone sound localization by the albino rat. Behav Neurosci 100:569–575, doi:10.1037/0735-7044.100.4.569, pmid:3741608.
    OpenUrlCrossRefPubMed
  34. ↵
    1. King AJ,
    2. Middlebrooks JC
    (2011) in The auditory cortex, Cortical representation of auditory space, eds Winer JA, Schreiner CE (Springer, New York), pp 329–341.
  35. ↵
    1. Lee CC,
    2. Middlebrooks JC
    (2011) Auditory cortex spatial sensitivity sharpens during task performance. Nat Neurosci 14:108–114, doi:10.1038/nn.2713, pmid:21151120.
    OpenUrlCrossRefPubMed
  36. ↵
    1. Lesica NA,
    2. Lingner A,
    3. Grothe B
    (2010) Population coding of interaural time differences in gerbils and barn owls. J Neurosci 30:11696–11702, doi:10.1523/JNEUROSCI.0846-10.2010, pmid:20810890.
    OpenUrlAbstract/FREE Full Text
  37. ↵
    1. Li LY,
    2. Li YT,
    3. Zhou M,
    4. Tao HW,
    5. Zhang LI
    (2013) Intracortical multiplication of thalamocortical signals in mouse auditory cortex. Nat Neurosci 16:1179–1181, doi:10.1038/nn.3493, pmid:23933752.
    OpenUrlCrossRefPubMed
  38. ↵
    1. Lingner A,
    2. Wiegrebe L,
    3. Grothe B
    (2012) Sound localization in noise by gerbils and humans. J Assoc Res Otolaryngol 13:237–248, doi:10.1007/s10162-011-0301-4, pmid:22245918.
    OpenUrlCrossRefPubMed
  39. ↵
    1. Lüling H,
    2. Siveke I,
    3. Grothe B,
    4. Leibold C
    (2011) Frequency-invariant representation of interaural time differences in mammals. PLoS Comput Biol 7:e1002013, doi:10.1371/journal.pcbi.1002013, pmid:21445227.
    OpenUrlCrossRefPubMed
  40. ↵
    1. Magezi DA,
    2. Krumbholz K
    (2010) Evidence for opponent-channel coding of interaural time differences in human auditory cortex. J Neurophysiol 104:1997–2007, doi:10.1152/jn.00424.2009, pmid:20702739.
    OpenUrlAbstract/FREE Full Text
  41. ↵
    1. Maier JK,
    2. Klump GM
    (2006) Resolution in azimuth sound localization in the Mongolian gerbil (Meriones unguiculatus). J Acoust Soc Am 119:1029–1036, doi:10.1121/1.2159429, pmid:16521765.
    OpenUrlCrossRefPubMed
  42. ↵
    1. Maier JK,
    2. Kindermann T,
    3. Grothe B,
    4. Klump GM
    (2008) Effects of omni-directional noise-exposure during hearing onset and age on auditory spatial resolution in the Mongolian gerbil (Meriones unguiculatus)–a behavioral approach. Brain Res 1220:47–57, doi:10.1016/j.brainres.2008.01.083, pmid:18343357.
    OpenUrlCrossRefPubMed
  43. ↵
    1. Maki K,
    2. Furukawa S
    (2005) Acoustical cues for sound localization by the Mongolian gerbil, Meriones unguiculatus. J Acoust Soc Am 118:872–886, doi:10.1121/1.1944647, pmid:16158644.
    OpenUrlCrossRefPubMed
  44. ↵
    1. Malhotra S,
    2. Hall AJ,
    3. Lomber SG
    (2004) Cortical control of sound localization in the cat: unilateral cooling deactivation of 19 cerebral areas. J Neurophysiol 92:1625–1643, doi:10.1152/jn.01205.2003, pmid:15331649.
    OpenUrlAbstract/FREE Full Text
  45. ↵
    1. McAlpine D,
    2. Jiang D,
    3. Palmer AR
    (2001) A neural code for low-frequency sound localization in mammals. Nat Neurosci 4:396–401, doi:10.1038/86049, pmid:11276230.
    OpenUrlCrossRefPubMed
  46. ↵
    1. Middlebrooks JC,
    2. Bremen P
    (2013) Spatial stream segregation by auditory cortical neurons. J Neurosci 33:10986–11001, doi:10.1523/JNEUROSCI.1065-13.2013, pmid:23825404.
    OpenUrlAbstract/FREE Full Text
  47. ↵
    1. Miller LM,
    2. Recanzone GH
    (2009) Populations of auditory cortical neurons can accurately encode acoustic space across stimulus intensity. Proc Natl Acad Sci U S A 106:5931–5935, doi:10.1073/pnas.0901023106, pmid:19321750.
    OpenUrlAbstract/FREE Full Text
  48. ↵
    1. Mossop JE,
    2. Culling JF
    (1998) Lateralization of large interaural delays. J Acoust Soc Am 104:1574–1579, doi:10.1121/1.424369, pmid:9745740.
    OpenUrlCrossRefPubMed
  49. ↵
    1. Pecka M,
    2. Brand A,
    3. Behrend O,
    4. Grothe B
    (2008) Interaural time difference processing in the mammalian medial superior olive: the role of glycinergic inhibition. J Neurosci 28:6914–6925, doi:10.1523/JNEUROSCI.1660-08.2008, pmid:18596166.
    OpenUrlAbstract/FREE Full Text
  50. ↵
    1. Reale RA,
    2. Brugge JF
    (1990) Auditory cortical neurons are sensitive to static and continuously changing interaural phase cues. J Neurophysiol 64:1247–1260, pmid:2258745.
    OpenUrlAbstract/FREE Full Text
  51. ↵
    1. Rébillat M,
    2. Benichoux V,
    3. Otani M,
    4. Keriven R,
    5. Brette R
    (2014) Estimation of the low-frequency components of the head-related transfer functions of animals from photographs. J Acoust Soc Am 135:2534–2544, doi:10.1121/1.4869087, pmid:24815237.
    OpenUrlCrossRefPubMed
  52. ↵
    1. Salminen NH,
    2. Tiitinen H,
    3. Yrttiaho S,
    4. May PJ
    (2010) The neural code for interaural time difference in human auditory cortex. J Acoust Soc Am 127:EL60–EL65, doi:10.1121/1.3290744, pmid:20136180.
    OpenUrlCrossRefPubMed
  53. ↵
    1. Schmitzer-Torbert N,
    2. Jackson J,
    3. Henze D,
    4. Harris K,
    5. Redish AD
    (2005) Quantitative measures of cluster quality for use in extracellular recordings. Neuroscience 131:1–11, doi:10.1016/j.neuroscience.2004.09.066, pmid:15680687.
    OpenUrlCrossRefPubMed
  54. ↵
    1. Scott BH,
    2. Malone BJ,
    3. Semple MN
    (2009) Representation of dynamic interaural phase difference in auditory cortex of awake rhesus macaques. J Neurophysiol 101:1781–1799, doi:10.1152/jn.00678.2007, pmid:19164111.
    OpenUrlAbstract/FREE Full Text
  55. ↵
    1. Siveke I,
    2. Pecka M,
    3. Seidl AH,
    4. Baudoux S,
    5. Grothe B
    (2006) Binaural response properties of low-frequency neurons in the gerbil dorsal nucleus of the lateral lemniscus. J Neurophysiol 96:1425–1440, doi:10.1152/jn.00713.2005, pmid:16571733.
    OpenUrlAbstract/FREE Full Text
  56. ↵
    1. Spitzer MW,
    2. Semple MN
    (1995) Neurons sensitive to interaural phase disparity in gerbil superior olive: diverse monaural and temporal response properties. J Neurophysiol 73:1668–1690, pmid:7643174.
    OpenUrlAbstract/FREE Full Text
  57. ↵
    1. Syka J,
    2. Popelár J,
    3. Kvasnák E,
    4. Astl J
    (2000) Response properties of neurons in the central nucleus and external and dorsal cortices of the inferior colliculus in guinea pig. Exp Brain Res 133:254–266, doi:10.1007/s002210000426, pmid:10968227.
    OpenUrlCrossRefPubMed
  58. ↵
    1. Tanaka H,
    2. Hachisuka K,
    3. Ogata H
    (1999) Sound lateralisation in patients with left or right cerebral hemispheric lesions: relation with unilateral visuospatial neglect. J Neurol Neurosurg Psychiatry 67:481–486, doi:10.1136/jnnp.67.4.481, pmid:10486395.
    OpenUrlAbstract/FREE Full Text
  59. ↵
    1. Ter-Mikaelian M,
    2. Sanes DH,
    3. Semple MN
    (2007) Transformation of temporal properties between auditory midbrain and cortex in the awake Mongolian gerbil. J Neurosci 27:6091–6102, doi:10.1523/JNEUROSCI.4848-06.2007, pmid:17553982.
    OpenUrlAbstract/FREE Full Text
  60. ↵
    1. Thomas H,
    2. Tillein J,
    3. Heil P,
    4. Scheich H
    (1993) Functional organization of auditory cortex in the Mongolian gerbil (Meriones unguiculatus). I. Electrophysiological mapping of frequency representation and distinction of fields. Eur J Neurosci 5:882–897, doi:10.1111/j.1460-9568.1993.tb00940.x, pmid:8281300.
    OpenUrlCrossRefPubMed
  61. ↵
    1. Victor JD,
    2. Purpura KP
    (1996) Nature and precision of temporal coding in visual cortex: a metric-space analysis. J Neurophysiol 76:1310–1326, pmid:8871238.
    OpenUrlAbstract/FREE Full Text
  62. ↵
    1. Vonderschen K,
    2. Wagner H
    (2009) Tuning to interaural time difference and frequency differs between the auditory arcopallium and the external nucleus of the inferior colliculus. J Neurophysiol 101:2348–2361, pmid:19261709.
    OpenUrlAbstract/FREE Full Text
  63. ↵
    1. Vonderschen K,
    2. Wagner H
    (2012) Transformation from a pure time delay to a mixed time and phase delay representation in the auditory forebrain pathway. J Neurosci 32:5911–5923, doi:10.1523/JNEUROSCI.5429-11.2012, pmid:22539852.
    OpenUrlAbstract/FREE Full Text
  64. ↵
    1. Vonderschen K,
    2. Wagner H
    (2014) Detecting interaural time differences and remodeling their representation. Trends Neurosci 37:289–300, doi:10.1016/j.tins.2014.03.002, pmid:24726910.
    OpenUrlCrossRefPubMed
  65. ↵
    1. Wightman FL,
    2. Kistler DJ
    (1992) The dominant role of low-frequency interaural time differences in sound localization. J Acoust Soc Am 91:1648–1661, doi:10.1121/1.402445, pmid:1564201.
    OpenUrlCrossRefPubMed
  66. ↵
    1. Winer JA,
    2. Schreiner CE
    (2005) The inferior colliculus (Springer, New York).
  67. ↵
    1. Winer JA,
    2. Schreiner CE
    (2011) The auditory cortex (Springer, New York).
  68. ↵
    1. Yamada K,
    2. Kaga K,
    3. Uno A,
    4. Shindo M
    (1996) Sound lateralization in patients with lesions including the auditory cortex: comparison of interaural time difference (ITD) discrimination and interaural intensity difference (IID) discrimination. Hear Res 101:173–180, doi:10.1016/S0378-5955(96)00144-X, pmid:8951442.
    OpenUrlCrossRefPubMed
  69. ↵
    1. Yin TC,
    2. Chan JC
    (1990) Interaural time sensitivity in medial superior olive of cat. J Neurophysiol 64:465–488, pmid:2213127.
    OpenUrlAbstract/FREE Full Text
  70. ↵
    1. Zwiers MP,
    2. Versnel H,
    3. Van Opstal AJ
    (2004) Involvement of monkey inferior colliculus in spatial hearing. J Neurosci 24:4145–4156, doi:10.1523/JNEUROSCI.0199-04.2004, pmid:15115809.
    OpenUrlAbstract/FREE Full Text
Back to top

In this issue

The Journal of Neuroscience: 34 (50)
Journal of Neuroscience
Vol. 34, Issue 50
10 Dec 2014
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Advertising (PDF)
  • Ed Board (PDF)
Email

Thank you for sharing this Journal of Neuroscience article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
The Neural Representation of Interaural Time Differences in Gerbils Is Transformed from Midbrain to Cortex
(Your Name) has forwarded a page to you from Journal of Neuroscience
(Your Name) thought you would be interested in this article in Journal of Neuroscience.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
View Full Page PDF
Citation Tools
The Neural Representation of Interaural Time Differences in Gerbils Is Transformed from Midbrain to Cortex
Lucile A.C. Belliveau, Dmitry R. Lyamzin, Nicholas A. Lesica
Journal of Neuroscience 10 December 2014, 34 (50) 16796-16808; DOI: 10.1523/JNEUROSCI.2432-14.2014

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Request Permissions
Share
The Neural Representation of Interaural Time Differences in Gerbils Is Transformed from Midbrain to Cortex
Lucile A.C. Belliveau, Dmitry R. Lyamzin, Nicholas A. Lesica
Journal of Neuroscience 10 December 2014, 34 (50) 16796-16808; DOI: 10.1523/JNEUROSCI.2432-14.2014
del.icio.us logo Digg logo Reddit logo Twitter logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Keywords

  • auditory cortex
  • inferior colliculus
  • interaural time differences
  • population coding
  • spatial hearing

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.

Related Articles

Cited By...

More in this TOC Section

Articles

  • Choice Behavior Guided by Learned, But Not Innate, Taste Aversion Recruits the Orbitofrontal Cortex
  • Maturation of Spontaneous Firing Properties after Hearing Onset in Rat Auditory Nerve Fibers: Spontaneous Rates, Refractoriness, and Interfiber Correlations
  • Insulin Treatment Prevents Neuroinflammation and Neuronal Injury with Restored Neurobehavioral Function in Models of HIV/AIDS Neurodegeneration
Show more Articles

Systems/Circuits

  • Behavioral state-dependent modulation of prefrontal cortex activity by respiration
  • Basolateral amygdala astrocytes are engaged by the acquisition and expression of a contextual fear memory
  • Inference of Electrical Stimulation Sensitivity from Recorded Activity of Primate Retinal Ganglion Cells
Show more Systems/Circuits
  • Home
  • Alerts
  • Visit Society for Neuroscience on Facebook
  • Follow Society for Neuroscience on Twitter
  • Follow Society for Neuroscience on LinkedIn
  • Visit Society for Neuroscience on Youtube
  • Follow our RSS feeds

Content

  • Early Release
  • Current Issue
  • Issue Archive
  • Collections

Information

  • For Authors
  • For Advertisers
  • For the Media
  • For Subscribers

About

  • About the Journal
  • Editorial Board
  • Privacy Policy
  • Contact
(JNeurosci logo)
(SfN logo)

Copyright © 2023 by the Society for Neuroscience.
JNeurosci Online ISSN: 1529-2401

The ideas and opinions expressed in JNeurosci do not necessarily reflect those of SfN or the JNeurosci Editorial Board. Publication of an advertisement or other product mention in JNeurosci should not be construed as an endorsement of the manufacturer’s claims. SfN does not assume any responsibility for any injury and/or damage to persons or property arising from or related to any use of any material contained in JNeurosci.