Original Research ARTICLE
Frequency discrimination and stimulus deviance in the inferior colliculus and cochlear nucleus
- 1Auditory Neurophysiology Laboratory, Institute of Neuroscience of Castilla y León, University of Salamanca, Salamanca, Spain
- 2Department of Neurobiology, Institute of Life Sciences, The Interdisciplinary Center for Neural Computation and the Edmond and Lily Safra Center for Brain Sciences, The Hebrew University of Jerusalem, Jerusalem, Israel
- 3Department of Cell Biology and Pathology, Faculty of Medicine, University of Salamanca, Salamanca, Spain
Auditory neurons that exhibit stimulus-specific adaptation (SSA) decrease their response to common tones while retaining responsiveness to rare ones. We recorded single-unit responses from the inferior colliculus (IC) where SSA is known to occur and we explored for the first time SSA in the cochlear nucleus (CN) of rats. We assessed an important functional outcome of SSA, the extent to which frequency discriminability depends on sensory context. For this purpose, pure tones were presented in an oddball sequence as standard (high probability of occurrence) or deviant (low probability of occurrence) stimuli. To study frequency discriminability under different probability contexts, we varied the probability of occurrence and the frequency separation between tones. The neuronal sensitivity was estimated in terms of spike-count probability using signal detection theory. We reproduced the finding that many neurons in the IC exhibited SSA, but we did not observe significant SSA in our CN sample. We concluded that strong SSA is not a ubiquitous phenomenon in the CN. As predicted, frequency discriminability was enhanced in IC when stimuli were presented in an oddball context, and this enhancement was correlated with the degree of SSA shown by the neurons. In contrast, frequency discrimination by CN neurons was independent of stimulus context. Our results demonstrated that SSA is not widespread along the entire auditory pathway, and suggest that SSA increases frequency discriminability of single neurons beyond that expected from their tuning curves.
Auditory neurons displaying stimulus-specific adaptation (SSA) decrease their response to high probability stimuli (standards) while maintaining responsiveness to rare ones (deviants, Ulanovsky et al., 2003). SSA is correlated with behavioral habituation (Netser et al., 2011; Gutfreund, 2012) and it has been proposed to underlie sensory memory for stimulation history (Ulanovsky et al., 2004). Neurons showing SSA have been found in the mammalian auditory pathway from the inferior colliculus (IC) up to the cortex (Ulanovsky et al., 2003; Pérez-González et al., 2005; Anderson et al., 2009; Malmierca et al., 2009; von der Behrens et al., 2009; Yu et al., 2009; Antunes et al., 2010; Lumani and Zhang, 2010; Reches et al., 2010; Taaseh et al., 2011; Zhao et al., 2011; Ayala and Malmierca, 2012; Duque et al., 2012; Anderson and Malmierca, 2013) as well as in bird midbrain and forebrain (Reches and Gutfreund, 2008, 2009). Originally, SSA was suggested to emerge in the auditory cortex and to be transmitted downstream to subcortical nuclei through the corticofugal pathway (Nelken and Ulanovsky, 2007), as subcortical SSA is mostly confined to the non-lemniscal regions (Malmierca et al., 2009; Antunes et al., 2010; Duque et al., 2012), the main target of the corticofugal projections (Malmierca and Ryugo, 2011). However, it has been recently shown that cortical deactivation does not affect SSA neither in the non-lemniscal auditory thalamus (Antunes and Malmierca, 2011) nor in the IC (Anderson and Malmierca, 2013), while SSA in lemniscal regions is minimal (Malmierca et al., 2009; Antunes et al., 2010; Bäuerle et al., 2011). Thus, SSA may be computed independently in the non-lemniscal pathway and in primary auditory cortex. Thus far, the existence of SSA has not been explored in auditory nuclei below the IC, where the lemniscal and non-lemniscal divisions first emerge.
Frequency discrimination has been widely explored in psychoacoustics (Nelson and Kiester, 1978; Sinnott et al., 1985; Syka et al., 1996; Talwar and Gerstein, 1998, 1999; Shofner, 2000; Witte and Kipke, 2005; Walker et al., 2009), but few studies tested the detection of frequency deviants by single neurons (Ulanovsky et al., 2003; Malmierca et al., 2009; von der Behrens et al., 2009). SSA has already been shown to result in a change in frequency discrimination performance by single neurons (Ulanovsky et al., 2003; Malmierca et al., 2009) but this relationship has not been thoroughly explored.
The main goal of our study is to compare the relationships between frequency discrimination and SSA in two neuronal populations; one at the IC that it is already known to exhibit SSA and the other at a lower auditory structure, the cochlear nucleus (CN) where SSA has not been explored thus far. For this purpose we assessed whether the probabilistic context affects frequency discrimination as judged by signal detection theory (Green and Swets, 1966) based on distributions of spike counts, and to what extent changes in frequency discriminability reflect the degree of SSA in these two stations. We show that SSA and the enhancement in neurometric frequency discrimination in the IC are strongly correlated and that both depend on the frequency separation and deviant probability in similar ways. Our results also demonstrated that SSA and context-dependent neuronal sensitivity are not present in CN supporting the hypothesis that SSA first emerge in non-lemniscal IC.
Materials and Methods
Experiments were performed on 71 adult female rats (Rattus norvergicus, Rj: Long–Evans) with body weights between 160 and 270 g. All experimental procedures were carried out at the University of Salamanca with the approval of, and using methods conforming to the standards of, the University of Salamanca Animal Care Committee. Anesthesia was induced (1.5 g/kg, i.p., 20% solution) and maintained (0.5 g/kg, i.p. given as needed) with urethane. Urethane was chosen as an anesthetic because of its effects on multiple aspects of neural activity, including inhibition and spontaneous firing, are known to be less than those of barbiturates and other anesthetic drugs (Hara and Harris, 2002). The respiration was maintained artificially (SAR-830/P Ventilator) monitoring the end-tidal CO2 level (CapStar-100). For this purpose, the trachea was cannulated and atropine sulfate (0.05 mg/kg, s.c.) was administered to reduce bronchial secretions. Body temperature was maintained at 38 ± 1°C by means of a heating blanket. Details of surgical procedures have been described previously (Hernández et al., 2005; Pérez-González et al., 2005; Malmierca et al., 2009; Antunes et al., 2010). The animal was placed inside a sound-attenuated room in a stereotaxic frame in which the ear bars were replaced by a hollow speculum that accommodated a sound delivery system.
Acoustic Stimuli and Electrophysiological Recording
Extracellular single unit responses were recorded from neurons in the IC and CN in two separate sets of experiments. For the IC recordings, a craniotomy was performed to expose the cerebral cortex overlying the IC and a tungsten electrode (1 – 2 MΩ) (Merrill and Ainsworth, 1972) was lowered through the cortex by means of a piezoelectric microdrive (Burleigh 6000 ULN). Neuron identification in the IC was based on stereotaxic coordinates, physiological criteria of tonotopicity, and response properties (Rees et al., 1997; LeBeau et al., 2001; Malmierca et al., 2003; Hernández et al., 2005; Pérez-González et al., 2005, 2006). An electrode dorsoventral penetration (with an angle of 20° from the frontal plane) through the central nucleus of the IC is identified by the stepwise progression from low to high frequencies (Malmierca et al., 2008) and by the constant presence of tonically firing units (Rees et al., 1997). Typical responses of the neurons in the cortices of the IC (i.e., dorsal, lateral, and rostral) are characterized by longer response latencies, predominantly on-phasic firing patterns and weaker tonic firing than those from the central nucleus. Cortical IC neurons commonly display broadly tuned, W-shaped, or other complex tuning curves (Lumani and Zhang, 2010; Geis et al., 2011; Duque et al., 2012) and a clear topographic organization of the frequencies along the dorsal cortex is not present (Malmierca et al., 2008; Lumani and Zhang, 2010). For the recording of CN neurons, part of the cerebellum was carefully aspirated to visually localize the dorsal cochlear nucleus (DCN). Glass micropipettes filled with 2 M NaCl (15 – 25 MΩ) or tungsten electrodes (1 – 2 MΩ) were advanced into the DCN. For some IC experiments and most of the CN recordings, an electrolytic lesion (10 – 15 μA for 10 – 15 s) was applied for subsequent histological verification of the recording site. Brains were fixed using a mixture of 1% paraformaldehyde and 1% glutaraldehyde diluted in 0.4 M phosphate buffer saline (0.5% NaNO3 in PBS). After fixation, tissue was cryoprotected in 30% sucrose and sectioned in the coronal or sagital plane at a thickness of 40 μm on a freezing microtome. Slices were Nissl stained with 0.1% cresyl violet to facilitate identification of cytoarchitectural boundaries. The CN units were assigned to one of the two main divisions (dorsal or ventral) of the nucleus using as reference the standard sections from a rat brain atlas (Paxinos and Watson, 2007).
Acoustic stimuli were delivered through a sealed acoustic system (Rees, 1990; Rees et al., 1997) using two electrostatic loudspeakers (TDT-EC1) driven by two TDT-ED1 modules. Search stimuli were pure tones or noise bursts monaurally delivered under computer control using TDT System 2 hardware (Tucker-Davis Technologies) and custom software (Faure et al., 2003; Pérez-González et al., 2005, 2006; Malmierca et al., 2008). The output of the system at each ear was calibrated in situ using a ¼ inch condenser microphone (Brüel and Kjær 4136, Nærum, Denmark) and a DI-2200 spectrum analyzer (Diagnostic Instruments Ltd., Livingston, Scotland, UK). The maximum output of the TDT system was flat from 0.3 to 5 kHz (~100 ± 7 dB SPL) and from 5 to 40 kHz (90 ± 5 dB SPL). The highest frequency produced by this system was limited to 40 kHz. The second and third harmonic components in the signal were 40 dB or more below the level of the fundamental at the highest output level (Hernández et al., 2005; Malmierca et al., 2009).
Action potentials were recorded with a BIOAMP amplifier (TDT), the 10× output of which was further amplified and bandpass-filtered (TDT PC1; fc: 0.5 – 3 kHz) before passing through a spike discriminator (TDT SD1). Spike times were logged at one microsecond resolution on a computer by feeding the output of the spike discriminator into an event timer (TDT ET1) synchronized to a timing generator (TDT TG6). Stimulus generation and on-line data visualization were controlled with custom software. Spike times were displayed as dot rasters sorted by the acoustic parameter varied during testing.
Once a neuron was isolated, the monoaural frequency response area (FRA), i.e., the combination of frequencies and intensities capable of evoking a response, was obtained by an automated procedure with 5 stimulus repetitions at each frequency (from 0.5 to 40 kHz, in 20 – 30 logarithmic steps) and intensity step (steps of 10 dB) presented randomly at a repetition rate of 4 Hz. The stimuli used to generate the tuning curves were pure tones with duration of 75 ms. The neuronal response to the combination of frequencies and intensities was plotted using MATLAB software (Mathworks, Inc.) and the best frequency (BF) and threshold for each neuron were identified.
Stimulus Presentation Paradigms
For all neurons, stimuli were presented in an oddball paradigm similar to that used to record mismatch negativity responses in human studies (Näätänen, 1992), and more recently in the cat auditory cortex (Ulanovsky et al., 2003, 2004), rat IC (Malmierca et al., 2009; Pérez-González et al., 2012) and auditory thalamus (Antunes et al., 2010; Antunes and Malmierca, 2011). Briefly, we presented two stimuli consisting of pure tones at two different frequencies (f1 and f2), that elicited a similar firing rate and response pattern at the same level of 10 – 40 dB SPL above threshold. Both frequencies were within the excitatory response area previously determined for the neuron. A train of 400 stimulus presentations containing both frequencies was delivered in three different sequences (Figure 1). The repetition rate of the train of stimuli for the IC neurons was 4 Hz, as it has been previously demonstrated to be suitable to elicit SSA in IC neurons of the rat (Malmierca et al., 2009). In the CN recordings, we explored repetitions rates of 4, 8, 12, and 20 Hz. Due to the different repetition rates used, the duration of the pure tones was 75 ms for the IC recordings (Hernández et al., 2005) and 25 ms for the CN recordings (in a few recordings at 4 and 8 Hz, tones lasted 75 ms as well), including a 5 ms rise/fall ramp for both cases.
Figure 1. The oddball stimulation paradigm. Two frequencies (f1, f2) were presented pseudo-randomly with different probabilities of occurrence. In Sequence 1 (S1), f1, and f2 occurred with equal probability (p50%), which served as a control condition. This condition is useful to see the neuron's tuning to the frequencies chosen. For the oddball condition, the probability of the frequencies was modified such that one frequency (f1, circles) was the standard tone, occurring with high probability, and the other (f2, squares) was the deviant tone, with low probability of occurrence (Sequence 2, S2). The probabilities of f1 and f2 were reversed in Sequence 3 (S3) in order to have each frequency presented as deviant and standard. We tested two probabilities for the deviant tone (pDev), 30 and 10%, so the corresponding probabilities for the standard (pStd) were 70 and 90%, respectively.
As shown in Figure 1, in Sequence 1 (S1) both frequencies were presented with the same probability of occurrence (equiprobable condition; p(f1) = p(f2) = 50%). In Sequence 2 (S2), one frequency (f1) was presented as the standard (i.e., high probability within the sequence: 90 or 70%); interspersed randomly among the standards were the deviant stimuli (i.e., low probability: 10 or 30%, respectively) at the second frequency (f2). After obtaining one data set, the relative probabilities of the two stimuli were reversed, with f2 as the standard and f1 as the deviant (S3). Sequences 2 and 3 constitute what we refer to as an oddball condition. The responses to the standard and deviant stimuli were normalized to spikes per stimulus, to account for the different number of presentations in each condition, because of the different probabilities. We tested several frequency separations between f1 and f2, expressed as frequency contrast Δf = (f2 − f1)/(f2 × f1)1/2 (Ulanovsky et al., 2003). As the frequency pairs were chosen to evoke similar firing rates in responses to both tones, Δf ranged from 0.02 to 3. The Δf values were grouped into three intervals: Δf ≤ 0.07, 0.07< Δf ≤ 0.2 and Δf > 0.2 (≤ 0.101, 0.101< Δf ≤ 0.288 and Δf > 0.288 octaves, respectively), in order to approximate to the values used in other studies, i.e., Δf = 0.04, 0.10, and 0.37 (Ulanovsky et al., 2003, 2004; Malmierca et al., 2009). The same paradigm was repeated changing the probability of the deviant tone (pDev = 10%, 30%) or the Δf. For the CN experiments, we only tested pDev = 10% and Δf = 0.1. The presentation of sequences at different deviant probabilities and at different repetition rates was randomized.
We measured the sharpness of the FRA of IC neurons calculating the bandwidth (BW) and Q-values at 10 and 40 dB SPL above the threshold as in our previous work (Hernández et al., 2005; Malmierca et al., 2009). The BW at n dB expresses the difference in kHz between the lower (FL) and upper (FU) frequencies of the FRA (BWn = FU − FL). The Q-value is calculated as the characteristic frequency (CF) divided by the BW at n dB above threshold (Qn = CF/BW).
The amount of SSA was quantified by two indices that have been used in previous studies (Ulanovsky et al., 2003, 2004; Malmierca et al., 2009; Antunes et al., 2010; Antunes and Malmierca, 2011; Pérez-González et al., 2012). The first index was the Frequency-Specific SSA Index (SI) defined as: SI(fi) = [d(fi) − s(fi)]/[d(fi) + s(fi)], where i = 1 or 2 and d(fi) and s(fi) are responses (as normalized spike counts) to frequency fi when it was deviant or standard, respectively. The second one was the Common-SSA Index (CSI) defined as CSI = [d(f1) + d(f2) − s(f1) − s(f2)]/[d(f1) + d(f2) + s(f1) + s(f2)], where d(f) and s(f) are responses to each frequency f1 or f2 when they were the deviant (d) or standard (s) stimulus, respectively. These indices reflect the extent to which the neuron responds more strongly to the frequencies when they are deviant compared to when they are standard. The possible SI and CSI values range from −1 to +1, being positive if the response to the deviant stimulus is greater and negative if the response to the standard stimulus is greater.
To estimate the neuronal sensitivity we performed a receiver operating characteristic (ROC) analysis (Tanner and Swets, 1954; Cohn et al., 1975; Fawcett, 2006; for a review of the use of ROC in psychometric and neurometric data analysis, see Stüttgen et al., 2011). This analysis has been previously used to measure the ability of CN units to signal changes in intensity (Shofner and Dye, 1989) and the sensitivity of IC units to interaural-time differences and binaural correlation (Skottun et al., 2001; Shackleton et al., 2003, 2005; Gordon et al., 2008). It is assumed that when different stimuli elicit different firing rates the response of a neuron provides the basis for discriminating between them. However, there is also a substantial variability in the responses to each stimulus, so the distributions of firing rates to similar stimuli overlap, and thus discrimination based upon firing rate will only be correct on a proportion of trials. The ROC analysis allows us to calculate the performance of the best possible discriminator between the two frequencies which is based on spike counts only. This discriminator is a function of the two probability distributions of spike counts in response to the two stimuli.
The ROC plots the probability of correct detection of f2 against the probability of “false alarm” detection of f2 when f1 occurred. Since detection is assumed to be based on spike counts only, trials have to be classified to one or the other frequency based solely on the evoked spike count. Thus, any discriminator between the two frequencies consists, in practice, of a list of spike counts that are assigned to frequency f1, with all other spike counts assigned to frequency f2 (we do not need to consider so-called “randomized rules” here, because we are only interested in the integral of the ROC, see below). In many studies, ROCs are calculated by a threshold on spike counts: all spike counts below the threshold are assigned to one frequency, and those above the threshold to the other. However, the lemma of Neyman and Pearson (Maris, 2012) requires spike counts to be assigned to frequencies based on their likelihoods, the ratio p(n|f2)/p(n|f1). For an optimal decision rule, a threshold is selected, and all spike counts whose likelihood is larger than that threshold are assigned to f2 (with the others assigned to f1). The probabilities of correct decision and false alarm for this decision rule can then be calculated in a straightforward manner. The ROC is obtained by calculating these probabilities while varying the threshold.
Then, we calculated the area under the ROC curve (AUC) as an estimate of the neural discriminability of frequency. The AUC corresponds to the probability of correct stimulus detection expected from an ideal observer in a two-alternative forced-choice psychophysical task (Green and Swets, 1966; Fawcett, 2006). Thus, sensitivity measured as AUC varies between 0.5 and 1, where 0.5 occurs when the spike count distributions for frequencies f1 and f2 are identical, and 1 indicates complete separation of the distributions. To compensate for sampling bias, we corrected each AUC value by performing 10,000 permutations of the original spike count distributions, assigned randomly to either f1 or f2, calculated the corresponding AUCs, and subtracted their mean value from the original AUC. Due to this correction some of the AUC values we report are smaller than 0.5. We also used the permutations test to estimate the probability of the AUC being significantly larger than 0.5. This way, we obtained one AUC value for the equiprobable condition (S1) and two AUC values for the oddball conditions (S2, S3). We used the mean AUC of S2 and S3 for the analyses instead of the maximum value as in previous works (Ulanovsky et al., 2003; Malmierca et al., 2009), in order to avoid an upward bias.
The CSI values were tested against zero by bootstrapping (1000 samples) in order to estimate a 95% confidence interval. Typically, CSI values smaller than 0.1 were not statistically different from zero (85% of all cases with CSI < 0.1 and 15% of the cases with CSI > 0.1). Thus, CSI values within the range of −0.1 to 0.1 were considered be due to random fluctuations in spike counts. This procedure provided a CSI cutoff comparable to other values previously set with different criteria (e.g., CSI = 0.18 for auditory thalamus of the rat; Antunes et al., 2010). It may be somewhat smaller than the cutoff in thalamus because of the lower variability in the responses of IC neurons (e.g., Chechik et al., 2006).
To investigate how frequency sensitivity is affected by the stimulation context we recorded the response of 224 well isolated single units in the IC and 51 units in the CN using an oddball paradigm. The frequency contrast (Δf ≤ 0.07, 0.07< Δf ≤ 0.2, Δf > 0.2) and probability of the deviant tone (pDev = 30% or 10%) were varied in IC recordings, and the repetition rate (4, 8, 12, and 20 Hz) in the CN. Additionally, an equiprobable context (p(f1) = p(f2) = 50%) was tested as control condition in both sets of experiments.
Neurons in the IC Show Different Degrees of SSA and Stimulus Discriminability
As might be expected from our previous studies (Pérez-González et al., 2005; Malmierca et al., 2009), neurons in the IC exhibited different degrees of SSA. Figure 2 shows the distribution of the CSI under different stimulus conditions in the current sample. The distributions of CSI are skewed toward positive values, and their medians are significantly different from zero (Signed Rank Test; p < 0.05) regardless of the condition tested (Figure 2). Positive CSI values reflect a stronger response to the deviant tone than to the standard one. The effects of frequency separation and deviant probability were tested using a Two-Way ANOVA on Δf × probability. There was a main effect of Δf (F(2, 489) = 18, p = 0) and of probability condition [F(1, 489) = 39, p = 0]. The interaction just failed to reach significance [F(2, 489) = 2.5, p = 0.08]. Post-hoc comparisons showed that the most positive CSI values were observed for deviant probability of 10% at the two highest frequency contrast intervals; 0.07< Δf ≤ 0.2 and Δf > 0.2. For the 10% probability condition, the CSIs increased significantly with increased frequency separation: CSI10%/Δf > 0.2 > CSI10%/0.07 < Δf ≤ 0.2 > CSI10%/Δf ≤ 0.07. On the other hand, the post-hoc comparisons did not show a significant difference between the average CSIs in the 30% condition and different frequency separations. There was also a significant difference due to changes in deviant probability for the two highest frequency separation intervals: CSI10%/0.07 < Δf ≤ 0.2 > CSI30%/0.07 < Δf ≤ 0.2; CSI10%/Δf > 0.2 > CSI30%/Δf > 0.2. This trend was emphasized by the higher percentage of neurons with CSI values larger than 0.1 when deviant probability was 10% compared to 30% (percentages indicated in Figure 2). From the six groups, only seven neurons (3.1%) showed CSI ≤ −0.1.
Figure 2. Distribution of stimulus-specific adaptation indices of IC neurons. Histograms of the common SSA index (CSI) displayed according to the frequency separation intervals (columns: Δf ≤ 0.07, 0.07 < Δf ≤ 0.2, Δf > 0.2) and the probability of the deviant tone (rows: pDev = 30 and 10%). The CSI was calculated from the responses recorded in S2 and S3. A CSI = 0 indicates an equal neuronal response when the tones were presented as deviant or standard, while positive and negative values represent higher firing rates when the tones were deviant or standard, respectively. For each stimulus condition the CSI values were tested against zero (solid line). The numbers next to the dashed line indicate the value of the median and the statistical significance (Signed Rank Test; *p < 0.05, ***p < 0.001). The distributions moved toward positive values when Δf was larger and pDev smaller. The percentages indicate the amount of neurons with CSI > 0.1.
Examples of individual IC neurons exhibiting different CSI values are shown in Figures 3–5. The deviant probability for the three cases was 10% and the frequencies tested (f1, f2) in these examples were separated by 0.144 octaves (Δf = 0.1) around its BF at an intensity of 10 – 50 dB above threshold.
Figure 3. Example of a non-adapting neuron in the IC. (A) Narrow FRA in color code for response magnitude. The tested frequencies (f1: 8.7 kHz, f2: 9.6 kHz, white crosses) were chosen around the BF (8.8 kHz) (arrowhead), at 45 dB SPL. (B) PSTH of the accumulated response to all the frequencies (0.5 – 40 kHz) and intensities (0 – 80 dB SPL) presented (1 ms bins). (C) Rate-level function at BF. (D–G) The responses of the neuron for each pair of stimuli for each of the three sequences (S1, S2, S3) are shown as dot raster plots (D), PSTH (3 ms bins) (E), spike probability distributions (F), and ROC curves (G). In the dot raster each dot represents the occurrence of a spike. The black bar under the PSTH and dot raster plots indicates the duration of the stimulus (75 ms). The probability of each frequency for each sequence is indicated on the upper left of the (E) panels. In the ROC curves (G) the dashed line corresponds to random guessing or no stimulus discrimination (AUC = 0.5), indicating complete overlap of the spike probability distributions. The red line represents the ROC curve calculated using the recorded data, the curves plotted in gray were obtained with the permutation method of the original spike count distributions, and the black line is represents the mean ROC curve of permutations. A total of 10,000 permutations were calculated, but for visual clarity only 100 curves are displayed. For each ROC curve, the area under the ROC curve (AUC) is shown corresponding to the original AUC value minus the mean AUC from permutations, as well as, the significance value for AUC > 0.5 (Permutation test; *p < 0.05). The repetition rate was 4 Hz and the frequency separation was 0.141 octaves. This neuron did not show SSA (CSI = 0.04, Bootstrapping; p > 0.05), displaying a very similar response to f1 and f2 across the three sequences regardless the probability of each tone.
Figure 3 illustrates a neuron with a CSI not significantly different from zero (CSI = 0.04; p > 0.05). This neuron had a narrow FRA with a narrow bandwidth both at 10 and 40 dB SPL above its threshold (BF = 8.8 kHz, Q10 = 5.62 and Q40 = 1.22) (Figure 3A). It had an onset firing pattern (Figure 3B) and showed a mixed/complex rate-level function (Figure 3C). The responses elicited in the equiprobable condition (S1) and oddball condition (S2 and S3) are shown as dot rasters (Figure 3D) as well as the peristimulus time histograms (PSTH) (Figure 3E). Figure 3F displays the corresponding spike-count distributions and the ROC curve is shown in Figure 3G. This neuron displayed a very robust and reliable response across the 400 stimulus presentations. In consequence, its spike count distributions were very different from Poisson distributions: while the average spike count is about 1, the probability of having zero spike counts is much smaller than that of either frequency evoking a single spike (for a Poisson distribution, these two probabilities should be approximately equal when the mean spike count is close to 1). The spike-count distributions for f1 and f2 were very similar, overlapping almost completely (Figure 3F), although the average spike count was slightly larger for f2 than for f1. The large overlap between these distributions resulted in AUC values very close to 0.5, but the very low variability resulted in an AUC that was significantly larger than 0.5 in the equiprobable condition. When f2 was the deviant, this difference was maintained, but when f1 was the deviant, the average spike count in response to f2 decreased slightly, enough to render the AUC not significantly different from 0.5 (AUC(S1) = 0.55, p = 0; AUC(S2) = 0.55, p = 0.01; AUC(S3) = 0.49, p = 0.65). Thus, the frequency discrimination capability of this neuron was poor in an equiprobable context and did not improve in an oddball stimulation context, consistent with its low CSI.
By contrast, neurons with high CSI values fired significantly differently in response to deviant and standard tones in the oddball condition. The neuron illustrated in Figure 4 had a CSI = 0.88 (p < 0.05). It was tuned to a wide range of frequencies (Figure 4A) reflected by its low Q-values (BF = 10 kHz, Q10 = 0.74 and Q40 = 0.27). This neuron also had an onset firing pattern, although it showed a large variability of first spike latency (FSL) (Figure 4B) and had a non-monotonic rate-level function (Figure 4C). During the equiprobable presentation of the tones (S1), this neuron adapted its response to both frequencies, and had a very low probability to respond at all (P≥1sps = 0.005). In the oddball condition, responses to the standard tone remained extremely sparse, but deviant trials did evoke a few spikes with higher probability. Thus, the overlap between the spike-count distributions was reduced substantially (probability of firing ≥ 1 sps in response to the deviant/standard was 0.4/0.017 and 0.23/0.014 for S2 and S3, respectively). As a result, the AUCs for the oddball conditions were higher than for the equiprobable condition [AUC(S1) = 0.5, p = 0.5; AUC(S2) = 0.7, p = 0; AUC(S3) = 0.6, p = 0].
Figure 4. Example of an adapting neuron in the IC. (A) FRA of a broadly tuned neuron with a BF of 10 kHz (arrowhead). The frequencies tested are indicated by the white crosses around the BF (f1: 9.5 kHz, f2: 9.6 kHz), at 0 dB SPL. (B–G) Same format as in Figure 3. This neuron showed strong SSA (CSI = 0.88, Bootstrapping; p < 0.05) reducing its firing to the high probability tone in S2 and S3 while still responding to the low probability one across most stimulus presentations. This differential firing is reflected in an AUC larger than 0.5 (Permutation test; *p < 0.05) in oddball sequences (S2 and S3).
The examples shown in Figures 3 and 4 are extreme cases, and neurons in the IC showed a continuous distribution of SSA as depicted in Figure 2. For example, Figure 5 illustrates a partially-adapting neuron (CSI = 0.5; p < 0.05) tuned to high frequencies (Figure 5A) and with a non-monotonic rate-level function (Figure 5C). The bandwidth of the FRA increased between 10 and 40 dB above threshold, respectively (BF = 29.9 kHz, Q10 = 7.81 and Q40 = 2.25). This neuron showed a poor, although significant discrimination capability at the equiprobable condition [AUC(S1) = 0.57, p = 0] which improved in the oddball condition [AUC(S2) = 0.61, p = 0; AUC(S3) = 0.76, p = 0].
Figure 5. Example of a partially-adapting neuron in the IC. (A) FRA from a neuron with a BF of 29.9 kHz (arrowhead). The frequencies tested were f1: 27.1 kHz and f2: 30 kHz, at 70 dB SPL. (B–G) Same format as in Figure 3. This neuron displayed a significant level of SSA (CSI = 0.5, Bootstrapping; p < 0.05), responding to both tones across the 400 stimulus trials. Although, the neuron displayed significant discriminability in the equiprobable condition (S1, AUC = 0.57) (Permutation test; *p < 0.05), this was improved under the oddball sequences (S2, AUC = 0.61; S3, AUC = 0.76).
Frequency Discriminability Depends on Stimulus Context in the IC
IC neurons were able to discriminate very similar frequencies even when both tones had the same probability of occurrence (p(f1) = p(f2) = 50%). The tested frequencies were selected online to evoke similar response magnitudes. Nevertheless, the noise in the estimation of response rates resulted in some imbalance between the responses to the two frequencies, leading to significant discriminability between them. The discriminability elicited under the equiprobable condition (AUC50%) across the three Δf intervals significantly differed from a mere random discrimination (AUC = 0.5, Signed Rank Test; p < 0.001) (Figure 6A). Furthermore, more than half of the neurons from each frequency separation group had AUC50% significantly larger than 0.5 (p < 0.05) (Figure 6A, indicated in percentage). In a substantial number of cases, AUC50% exceeded 0.71 (24.1, 23.2, and 41.6% for the three Δf groups), the generally accepted definition of a threshold (Green and Swets, 1966). Neurons with AUCs above this threshold for the smallest frequency contrast interval (Δf ≤ 0.07) had narrower bandwidths (Q10 = 6.23 ± 5.43) that the rest of neurons (Q10 = 3.53 ± 5.47) (Signed Rank Test; p < 0.05).
Figure 6. Neurometric performance under equiprobable and oddball conditions of IC neurons. (A) Distributions of the AUC values for the equiprobable condition (AUC50%) indicating the median (dashed line) significantly differs from 0.5 (Signed Rank Test; ***p < 0.001). The percentage of neurons whose AUC50% was higher than 0.5 is indicated for each panel (Permutation test; p < 0.05). (B) Scatter plots showing the neurometric performance for frequency discrimination expressed as percentage correct under the oddball condition (rows: pDev = 30, 10%) versus the equiprobable condition (pf1 = pf2 = 50%), for each frequency contrast interval (columns: Δf ≤ 0.07, 0.07 < Δf ≤ 0.2, Δf > 0.2). Separately are represented the neurons with CSI ≤ 0.1 (gray circles) and CSI > 0.1 (dark crosses). For the oddball condition, the percentage correct corresponds to the mean AUC value obtained from S2 and S3. The number of neurons above and below the diagonal line (equal performance in both conditions) is indicated by the inset on the bottom right of each panel. (C) Sensitivity curves of individual neurons expressed as percentage of change elicited when pDev = 10% and 30% regarding the pf1 = pf2 = 50% condition.
To address the central question of this paper, Figure 6B compares the percent correct (as estimated by AUC) in the oddball and equiprobable condition for each neuron. For the oddball condition, we used the mean discriminability (AUCoddball) from the values elicited in the two oddball sequences since there was not significant difference in the AUC values elicited in S2 and S3 (Rank Sum Test; p > 0.05). Neurons whose discriminability was unaffected in the oddball condition fell along the diagonal line. Neurons under the diagonal line showed a better discriminability in the equiprobable condition. By contrast, neurons that improved their discriminability in the oddball paradigm were located above the diagonal. Neurons with CSI > 0.1 are marked by crosses, the others are marked by circles. When pDev = 10%, there was a larger proportion of neurons with CSI > 0.1 than neurons with CSI ≤ 0.1 that showed improved discriminability in the oddball condition (χ2 = 58.6, df = 1, p < 0.001), but these proportions did not depend on frequency separation (χ 2 = 5.4, df = 2, p = 0.07). For this probability condition, the AUCs of neurons with CSI ≤ 0.1 were slightly, although significantly, smaller in the oddball than in the equiprobable condition (46/96, neurons above and below the bisecting line, respectively, for all frequency separation classes together). This effect was due presumably to the poorer sampling of the spike count histograms for the deviant stimuli in the oddball condition. On the other hand, AUCoddball increased substantially for neurons with CSI > 0.1 (161/59 neurons above and below the bissecting line, respectively). The increase resulted in many neurons whose frequency discrimination was below threshold in the equiprobable condition (AUC50% < 0.71) and that exceed threshold in the oddball conditions (AUCoddball > 0.71). Within this subset of neurons, there are cases in which the neurometric performance reached values close to 100% correct in the oddball condition. Such cases were much more common at the largest frequency contrasts (0.07< Δf ≤ 0.2 and Δf > 0.2). For pDev = 30%, the discriminability did not change consistently relative to the equiprobable condition, and proportions of neurons with slight increase or decrease in discriminability were as common in the different frequency difference classes (χ2 = 5.4, df = 2, p = 0.07) and among CSI classes (χ2 = 3.9, df = 1, p = 0.05).
In order to verify whether the same trend was observed at the level of single neurons, we obtained the individual “sensitivity curves” for the neurons that were tested under all probabilities conditions (50, 30, and 10%) and for the same frequency pairs (Figure 6C). The discriminability increment was expressed as the percentage of change in AUCoddball relative to the discriminability displayed under the equiprobable condition (AUC50%). These sensitivity curves revealed a considerable diversity in the neuronal performance. Both neuron identity and stimulus probability had a significant effect on the discrimination capability for the intermediate Δf interval [Two-Way ANOVA on stimulus probability × neuron, significant main effect of stimulus probability: F(2, 128) = 7.7, p < 0.001], but for the smallest and largest Δf the main effect of stimulus probability was not significant.
Since some neurons under the equiprobable condition showed significant discriminability values that exceeded a mere random response (Figure 6A), we took into account this neuron-specific tuning. We calculated the discriminability enhancement index (DEI) as the difference between the discriminability elicited in the oddball condition and that elicited in the equiprobable one (DEI = AUCoddball − AUC50%). DEI ranges from −0.5 to 0.5, with positive values indicating an improvement in discriminating two stimuli under an oddball context. The comparison of the mean population values of DEI across all stimulus combinations (Two-Way ANOVA, stimulus probability × Δf) demonstrated that it was affected by the frequency separation [F(2, 489) = 5.72, p < 0.01] but not by stimulus probability [F(1, 489) = 3.71, p = 0.055], with no interaction between those factors [F(2, 489) = 2.1, p = 0.12] (Figure 7).
Figure 7. Stimulus discriminability enhancement of IC neurons across different stimulus conditions. Box plots of the discriminability enhancement under the oddball condition (DEI) showing the mean (dashed line) and the median (solid line) values, as well as, the 5th and 95th outliers. All the mean values were positive (except for the 10%/Δf ≤ 0.07 condition), reflecting a better stimulus discrimination when one of the frequencies is presented as a deviant tone, that is, with low probability of ocurrence (30 or 10%). The DEIs were only affected by the frequency separation [F(2, 489) = 5.72, p < 0.01] (Two-Way ANOVA, deviant probability × frequency separation).
IC Neurons with High SSA Showed a Greater Discriminability Enhancement Under Oddball Conditions
Finally, we analyzed the relationship between the two metrics used to quantify the neuronal responses in order to explore whether or not the change in stimulus discrimination can be predicted by their SSA index. This analysis demonstrated a strong positive correlation between the degree of adaptation (CSI) and the enhancement in the frequency discriminability (DEI) shown by neurons under the condition with the lowest deviant probability, that is, when pDev = 10% (Spearman's rho; p < 0.001) (Figure 8). The great majority of neurons with CSI < 0.1 had discrimination indices clustered around the origin (gray circles). By contrast, most neurons with CSI > 0.1 (crosses) had a positive DEI, indicating that adapting neurons had better frequency discrimination for oddball sequences, and furthermore, there was a tendency for larger CSI values to be associated with larger DEI values.
Figure 8. The stimulus discriminability reflects the degree of stimulus-specific adaptation exhibited by IC neurons. Correlation of the stimulus discriminability enhancement (DEI) and the SSA index (CSI). In gray circles are represented the neurons with CSI ≤ 0.1 and in dark crosses those with CSI > 0.1. The linear correlation was reflected by the Spearman's correlation coefficient (rs), whose strength varied according to the frequency contrast (columns) and probability of the deviant tone (rows). *p < 0.05, ***p < 0.001.
Relation Between the Width of Frequency Tuning and the SSA or Discriminability Exhibited by IC Neurons
Previous reports demonstrated a differential expression of SSA through the lemniscal and non-lemniscal subdivisions of the IC (Pérez-González et al., 2005; Malmierca et al., 2009; Ayala and Malmierca, 2012; Duque et al., 2012) and medial geniculate body (MGB) (Antunes et al., 2010) of the rat. Neurons in the cortical regions of the IC exhibit broader FRAs than the ones from the central nucleus and the broader the response area is, the higher the SSA levels are (Duque et al., 2012). In order to test whether or not this relationship is found in our neuronal sample, we analyzed the width of response areas as a function of the level of SSA.
Figure 9A displays the bandwidths at 10 and 40 dB SPL above threshold (BW10, 40, respectively) for the lowest deviant probability (pDev = 10%) as a function of the CSI. The group of CSI ≤ 0.1 included all neurons that were considered to lack SSA. The other CSI cutoffs were selected to have approximately equal-size groups. It is interesting to note that there were neurons with very broad bandwidth already at 10 dB above threshold. We performed an analysis of covariance of BW, with level above threshold (10 or 40 dB SPL) and frequency separation as qualitative factors and CSI as a continuous factor. We found a highly significant effect of CSI [F(1, 688) = 37.5, p = 0]. The slope of the dependence of BW on CSI indicated that BW increased on average by 6.6 kHz as CSI increased from zero to one. The main effect of frequency separation was not significant, [F(2, 688) = 0.54, p = 0.6], while the level above threshold had, as expected, a significant effect on BW [F(1, 688) = 77.4, p = 0]. There was a significant interaction between the CSI slope and level above threshold [F(1, 688) = 5.8, p = 0.01], and post-hoc comparison indicated that CSI slopes at 10 and 40 dB above threshold were significantly different (p < 0.05).
Figure 9. Bandwidth of frequency response areas and SSA level of IC neurons. (A) Box plots of the bandwidth (BW) values grouped into CSI ranges for the three Δf intervals with the mean (dashed line) and median values (solid line) indicated, as well as, the 5th and 95th outliers. There is an increment in the bandwidths as neurons have higher CSI [F(1, 688) = 37.5, p = 0] and as the level increased [F(1, 688) = 77.4, p = 0] (Analysis of covariance of BW with level, Δf and CSI as factors). (B) Box plots of BW values grouped into DEI ranges. The BW increases for neurons that displayed higher discriminability improvement under the oddball condition [F(1, 688) = 26.7, p = 0]. Same format as panel (A). (Analysis of covariance of BW with level, Δf and DEI as factors).
As expected from the positive correlation between DEI and CSI (Figure 8), a significant effect of DEI on BW was also found [F(1, 688) = 26.7, p = 0] (Figure 9B). In consequence, a greater neuronal discriminability in the oddball condition is associated with a wider frequency integration range. DEI also had significant interaction with level above threshold [F(1, 688) = 3.9, p = 0.048] (analysis of covariance of BW with level, frequency separation, and now with DEI as a continuous factor).
In selected cases, we made electrolytic lesions in the IC and determined that we recorded neurons from central nucleus (n = 9) and from cortical regions (n = 16). Within this very limited sample, the central nucleus neurons had a CSI of 0.11 ± 0.21 and a DEI of −0.001 ± 0.12 (median ± SD). For the cortical neurons, the CSI and DEI were of 0.34 ± 0.3 and of 0.03 ± 0.13, respectively. However, this number of histological localizations was insufficient to guarantee a reliable study to correlate SSA and discriminability degree across the different subdivisions of the IC. Figure 10A showed an example of a typical lesion located in the lateral cortex of the IC (Loftus et al., 2008; Malmierca et al., 2011).
Figure 10. Histological identification of recording sites of IC and CN neurons. (A) Example of recording sites marked with an electrolytic lesion (arrowheads) in the lateral cortex of the IC at 1.9 mm lateral, according to Paxinos and Watson (2007). Two different tracts are indicated by Tr #1 and Tr #2. (B,C) Recording sites (arrowheads) and tracts (Tr) located in DCN (at 11.52 mm from bregma) and VCN (at 11.04 mm from bregma), respectively. The slices were Nissl stained and cut at 40 μm in a sagital (A) and coronal plane (B,C). Scale bars of 500 μm. D, dorsal; C, caudal; M, medial.
CN Neurons do not Exhibit SSA and Their Frequency Discriminability is not Sensitive to a Probability Context
Since SSA is present in the IC, we wanted to explore whether SSA is already ubiquitously expressed earlier. We recorded 51 CN neurons to test whether SSA is exhibited by single-units and if so, whether adaptation strength correlates with neuronal sensitivity as shown for the IC neurons.
A total of 44 neurons out of 51 were localized and assigned to the ventral cochlear nucleus (VCN) (n = 10) or DCN (n = 34). The histological reconstruction for the remaining 7 neurons was not possible. Figure 10B shows the electrolytic lesion in a Nissl-stained section, illustrating the recording site of the neuron displayed in Figure 11 and located in the DCN. Another example of recording site, in the VCN, is shown in Figure 10C. The recorded neurons had a wide variety of firing patterns and rate-level functions, as have been described in detail before (Stabler et al., 1996). More than the half of neurons in the DCN (21/34) displayed non-monotonic rate-level functions (5/10 in the VCN). Our sample of DCN neurons included chopper (n = 13), primary like (n = 10), pause/build (n = 6), and onset (n = 5) firing patterns. In the VCN, all firing patterns except the choppers were present (primary like, n = 5; pause/build, n = 2; onset, n = 3). Figure 11 shows the response of a DCN neuron with a typical V-shaped FRA with a low-frequency tail. The evoked activity was robust across the 400 trials of the equiprobable (S1) and deviant sequences (S2, S3). The neuron showed significant discriminability under the equiprobable condition (AUC(S1) = 0.6, p < 0.05) which did not greatly improve under the oddball sequences (AUC(S2) = 0.62, AUC(S3) = 0.53, p < 0.05). This neuron failed to show SSA at a repetition rate of 4 Hz, as well as at faster stimuli presentation rates of 8 and 20 Hz (CSI = 0, p > 0.05).
Figure 11. Example of a CN neuron. The format for all panels is the same as in Figures 3–5. (A) The BF was 23.4 kHz (indicated by the arrowhead) and the tested frequencies (f1 = 23.4 kHz and f2 = 25.9 kHz, white crosses) differ by 0.144 octaves. The PSTH of the accumulated response to all the frequencies (0.5 – 40 KHz) and intensities (0 – 80 dB SPL) presented (1 ms bins), as well as, the rate-level function at BF are shown in (B) and (C), respectively. The neuron exhibited a low CSI = 0.1 (Bootstrapping; p > 0.05) with a robust response across trials (D,E) and its frequency discriminability was not sensitive to the oddball condition (F,G). The AUC values in all conditions were slightly higher than 0.5 (Permutation test; *p < 0.05).
SSA was not present in the neuronal population recorded in CN (Figure 12). We used faster repetition rates than in the IC since SSA seems to increase monotonically with stimulation rate (Malmierca et al., 2009; Antunes et al., 2010; Zhao et al., 2011; Patel et al., 2012). Regardless of the extreme repetitions rates, the strength of the neuronal response was equal for deviants and for standards stimuli (Signed Rank Test; p > 0.05) (Figure 12A) resulting in SI values clustered around zero (Figure 12B). There were no differences between the CSIs elicited by VCN and DCN neurons for any repetition rate tested (Rank Sum Test; p > 0.05). FSL is also affected by probability condition in IC, being shorter to the deviant stimulus regardless of the frequency tested (f1 or f2) (Malmierca et al., 2009; Pérez-González and Malmierca, 2012; Pérez-González et al., 2012). For CN neurons, the vast majority of FSL to deviant and to standard was almost equal and no differences in the median FSL between them was observed (Signed Rank Test; p > 0.05) (Figure 12C). The median of the FSLs was 9.62 ± 4.9 ms (range: 3.4 – 29.5 ms) and 9.82 ± 4.7 ms (range: 3.9 – 28.2 ms) for deviant and standard tone, respectively. These latencies are clearly shorter than the latencies of IC neurons (FSL to deviant: 26.1 ± 13.2 ms; range: 7.5 – 72 ms, FSL to standard: 29.6 ± 13.2 ms, range: 7.3 – 74.5 ms; from Malmierca et al., 2009). Although some neurons showed significant CSI > 0.1 (0.11 – 0.28) at 4 (n = 3), 8 (n = 4), 12 (n = 2), and 20 Hz (n = 5) (most of them from the DCN, n = 5), the average SSA indices were not significantly different from zero (Signed Rank Test; p > 0.05) nor they were sensitive to the rate of stimulation (Kruskal–Wallis Test; p > 0.05) (Figure 12D). Finally, CSI was not affected by increasing the frequency separation from Δf = 0.1 to Δf ≥ 0.2 (0.2 − 0.37) (Signed Rank Test; p > 0.05) (Figure 12E).
Figure 12. CN neurons do not exhibit SSA. (A) Spike count to deviant (dev) and to standard (std) stimuli elicited at different repetition rates represented in a color code. This color code is the same for (B–E) panels. (4 Hz, n = 47; 8 Hz, n = 29; 12 Hz, n = 9; 20 Hz, n = 26) (B) Scattergraph of the Frequency-Specific SSA indices (SI) for f1 and f2 presented at different repetitions rates. This index reflects the normalized spikes counts elicited when each frequency was the deviant tone regarding the response evoked when it was the standard one. (C) Scattergraph of the median first spike latency (FSL) for f1 and f2 when they were the deviant (FSLdev) or the standard (FSLstd) stimulus. In the right column, are displayed the box plots of the median FSL for the population of neurons to deviant (red) and to standard (blue) stimulus across repetition rates. The n for this panel is the double of the number of neurons tested, since two frequencies were tested as deviant and as standard stimulus for each neuron. (D,E) Box plots of the Common–SSA index (CSI) when increasing the rate of stimulation and scattergraph of CSI when varying the frequency separation factor (Δf), respectively.
In parallel with the lack of SSA, frequency discrimination was not affected by changes in tone probability in this neuronal population. The estimated correct detection in the oddball condition remained very similar to that elicited in the equiprobable one for most of the CN neurons (Figure 13A), and no improvement in frequency discriminability was elicited at the population level for any repetition rate group (Signed Rank Test; p > 0.05) (Figure 13B). Thus, the DEI was essentially zero and insensitive to increments in the repetition rate (Kruskal–Wallis Test; p > 0.05) (Figure 13C). As expected, it was not correlated with the SSA index (Spearman's correlation) (Figure 13D).
Figure 13. Stimulus sensitivity in CN neurons. (A) Neurometric boxes displaying the percentage of correct identifications in the equiprobable and oddball condition (pDev = 10%) for different repetitions rates (4, 8, 12, and 20 Hz). (B) Box plots indicating that no improvement was elicited under the oddball condition regarding the equiprobable one (Signed Rank Test; p > 0.05). (C) Box plots that indicate that the discrimination enhancement indexes (DEI) remained at zero and did not change across repetition rate (Kruskal–Wallis Test; p > 0.05). (D) No positive correlation between the CSI and DEI was found (Spearman's correlation coefficient, rs; *p < 0.05).
Our study demonstrates that sensitivity to frequency in IC neurons but not in CN neurons depends on probability context. Changes in frequency discriminability in IC neurons reflected the level of SSA they exhibit. Both the CSI and DEI values increased with frequency separation and DEI tended to be positive (Figures 7 and 8). The lack of effect of probability context in CN was related to the lack of SSA in the neuronal sample we recorded from (Figure 12).
Strong SSA is Exhibited by IC Neurons but not by CN Neurons
The strength of SSA reported here is similar to that reported previously by Malmierca et al. (2009) for IC neurons. This is not surprising, since we used a similar experimental preparation including animal model, parameters, and paradigm of stimulation (presentation rate: 4 Hz; tone duration: 75 ms; random presentation of tones). Other studies also have examined SSA in the IC, although as these studies have used different stimulation paradigms (e.g., Pérez-González et al., 2005; Lumani and Zhang, 2010), different metrics to quantify SSA (Pérez-González et al., 2005) or different stimulus repetition rates (Zhao et al., 2011), a quantitative comparison is difficult. According to the sample of histological verifications of the recording sites (Figure 10A) and taking into account the distribution of CSI (Figure 2), we recorded neurons from the central nucleus, as well as from the cortical non-lemniscal regions of IC. SSA varies as a continuum throughout the entire IC and it is strong and widespread in the non-lemniscal regions of the IC (Malmierca et al., 2009; Duque et al., 2012) and MGB (Antunes et al., 2010), being low or almost absent in the lemniscal subdivisions, the central nucleus of the IC and ventral MGB. Also, the neurons in the cortex of the IC exhibit broader response areas than those from the central nucleus (Malmierca et al., 2008, 2009; Geis et al., 2011; Duque et al., 2012). In agreement with these results, we show here that neurons with wider bandwidths (values as high as 30–40 kHz) showed the strongest SSA (Figure 9). Thus, the convergence of ascending, narrowly tuned frequency inputs with different frequency selectivity could be a major mechanism underlying SSA. In support to this idea, Taaseh et al. (2011) and Mill et al. (2011) showed that individual inputs showing simple fatigue could result in SSA. Beyond this mechanism, SSA could be further refined through the local inhibitory circuits and descending inputs from higher auditory centers. In this respect, a modulatory role of postsynaptic GABAA receptors in shaping SSA in the IC has already been demonstrated (Pérez-González and Malmierca, 2012; Pérez-González et al., 2012).
Considering that (1) the IC is the locus of convergence for most inputs originating at lower auditory brainstem nuclei and the locus where the lemniscal and non-lemniscal pathways appears (Malmierca et al., 2003; Lee and Sherman, 2010, 2011), and that (2) our results demonstrated the lack of widespread SSA at the CN (Figure 12), it is tempting to suggest that cells exhibiting strong SSA in the subcortical pathways first emerge in the non-lemniscal IC. Two possible confounds currently limit this hypothesis. First, a decrease in the responsiveness and changes in the response variability of auditory cortical neurons caused by the anesthesia (Kisley and Gerstein, 1999; Harris et al., 2011) could also result in the absence of strong SSA in CN. This would be the case, for example, if SSA in the CN were dependent on descending projections for its generation. However, this possibility seems unlikely since previous studies have demonstrated that SSA at the IC (Anderson and Malmierca, 2013) and MGB (Antunes and Malmierca, 2011) persist even if the corticofugal pathway is reversibly deactivated. Second, the CN has multiple distinct physiological response types which are well-correlated with anatomical and cellular characteristics. While neurons were recorded in both VCN and DCN, currently there is no detailed classification of CN neurons in the anesthetized rat. Therefore, we cannot rule out that we recorded from all response types in this study. Indeed, because across-frequency integration seems to be important in SSA, CN neuronal types that show frequency convergence and that project to the IC, e.g., some multipolar cells or small cells from the cap area (Winter and Palmer, 1995; Jiang et al., 1996; Palmer et al., 1996; Malmierca et al., 2002; Cant and Benson, 2003) might potentially be capable of showing high levels of SSA as well. To rule out this possibility a detailed morphological and physiological study is necessary in the future. Finally, the presence of SSA in brainstem nuclei between the CN and IC also remains to be tested.
Neuronal Sensitivity of IC
We show here that the vast majority of IC neurons discriminate between nearby tones around BF even when they occur with equal probability (Figures 6A,B; upper left quadrants). While many of the AUCs were significantly larger than 0.5, they also tended to be smaller than 0.71, the standard definition of a psychophysical threshold. Note that other pairs of frequencies within the FRA with the same frequency difference could give rise to larger AUCs. Thus, our results for the equiprobable case should be seen as a lower bound on the frequency discrimination capabilities of IC neurons. Even with the biased selection of frequencies to test, a small population reached AUC > 0.71 for frequency separations as small as Δf = 0.07, very close to the psychophysical thresholds of rats (e.g., Talwar and Gerstein, 1998, 1999). Interestingly, these neurons also showed narrower bandwidths than the rest. The narrow bandwidth could result in large changes in firing rates for nearby frequencies, leading to the high AUCs (Gordon et al., 2008). Such narrowly-tuned inputs could also account for the hyperacuity of MGB and cortical neurons in oddball conditions, as previously reported (Ulanovsky et al., 2003; von der Behrens et al., 2009).
Frequency discrimination in IC depended on context, being larger when the stimuli had decreased probability (Figures 6B,C and 7). Robust increases in discrimination in the oddball conditions occurred for the lower deviant probability (pDev = 10%) and the larger frequency separations (Δf > 0.07). Note that we used here the mean AUC across the two oddball sequences, rather than the maximal one as used previously (Ulanovsky et al., 2003; Malmierca et al., 2009). The mean AUC is a more conservative estimate of frequency discriminability, and its use may explain why we did not observe extreme discrimination performance as reported previously in IC (Malmierca et al., 2009 their Figure 7, neurons in the upper left corner). Either way, these results emphasize the influence of context on sensory processing as early as in the IC as has been demonstrated before for the processing of interaural phase (Spitzer and Semple, 1991, 1993, 1998; McAlpine et al., 2000), level differences (Sanes et al., 1998), monaural frequency transitions (Malone and Semple, 2001) and simulated motion (Wilson and O'Neill, 1998).
Frequency Discriminability Enhancement Reflects the Degree of SSA
We found a strong correlation between the discriminability enhancement and the degree of SSA, but only for the condition with the lowest deviant probability (pDev = 10%) and larger frequency separation (Δf > 0.07) (Figure 8). These are also the conditions that had higher CSI. This positive correlation is expected from the design of the experiment. The two frequencies were selected to evoke equivalent responses in the equiprobable condition, and in the oddball condition they were expected to evoke different responses. In consequence, we expected a substantial overlap between the spike count distributions in the equiprobable condition, but a decreasing overlap in the oddball condition. Indeed, DEI depended on deviant probability and frequency separation very similarly to CSI (Figure 7). Finally, the absence of SSA and null enhancement in deviant detectability by CN neurons reinforced the notion that deviant discriminability is a functional consequence of SSA (Figures 12 and 13).
Nevertheless, we also found that neurons with CSI ≤ 0.1 showed a significant decrease in AUC in the oddball condition (Figure 6B). This decrease was due to larger corrections for the AUCs obtained under the oddball condition than to the AUCs under the equiprobable one. This trend should be seen as a negative bias in the estimation of the AUC under oddball conditions. Given that as a rule AUC increased with decreasing deviant probability, our conclusions should be considered as conservative.
Neuronal responses in auditory cortex are plastic at many different time scales (Condon and Weinberger, 1991; Kilgard and Merzenich, 1998, 2002; Fritz et al., 2003; Ulanovsky et al., 2003; Froemke et al., 2007). Here we demonstrate that neurons in IC show some sort of short-term plasticity under similar conditions to neurons in MGB and A1. As previously suggested (Antunes et al., 2010), the non-lemniscal regions of the IC could transmit SSA to the non-lemniscal MGB neurons, which in turn would project to the superficial layers of AC (Cetas et al., 1999; Huang and Winer, 2000; Anderson et al., 2009). Neurons in the medial division of the MGB have large-diameter axons that are known to terminate primarily in layer I of the auditory cortex in both primary and secondary cortical fields. For example, in the somatosensory cortex Cauller and Connors (1994) observed strong excitatory effects on pyramidal cells present in layers II, III and V to be mediated by long horizontal axons located in layer I. Further experiments are required in order to check this possibility. Thus, at all levels of the auditory pathways, context-dependence of the responses could serve for adjusting the neural code to match the statistics of the input signal to produce an efficient representation of auditory scene. Similarly, the changes in responses as a function of tone probability could serve in the processes of auditory scene analysis. Indeed, auditory stream segregation is also sensitive to frequency separation and presentation rate (e.g., Fishman et al., 2004; Fishman and Steinschneider, 2010). Moreover, there is evidence suggesting the involvement of pre-attentive neural process in auditory stream segregation (Winkler et al., 2003). Thus, SSA in IC may increase the saliency of low-probability signals, helping to segregate them by reducing the ambiguity of the neuronal representations for downstream read-out mechanisms.
Interestingly, our results suggest that the initial locus for the computation of SSA is not at the very first stations of the auditory pathway, e.g., the CN. Thus, the picture of the auditory system that emerges here reinforces the idea that the initial coding of sounds is purely based on their short-term physical characteristics, and sensitivity to longer contexts that is required for higher-order processing, efficient coding, and auditory scene analysis appears only later.
Manuel S. Malmierca and David Pérez-González designed research; Yaneri A. Ayala performed research; Yaneri A. Ayala, Daniel Duque, David Pérez-González and Israel Nelken analyzed data; Yaneri A. Ayala, David Pérez-González, Israel Nelken and Manuel S. Malmierca wrote the paper.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Financial support was provided by the Spanish MEC (BFU2009-07286), EU (EUI2009-04083, in the framework of the ERA-NET NEURON Network of European Funding for Neuroscience Research) to Manuel S. Malmierca. Israel Nelken was supported by a grant from the Israeli Ministry of Ministry of Health in the framework of the ERA-NET Network. Daniel Duque held a fellowship from the Spanish MINECO (BES-2010-035649). David Pérez-González held a postdoctoral fellowship from the Botín Foundation. We thank Drs. Alan Palmer and Trevor M. Shackleton for their constructive comments on a previous version.
Bäuerle, P., von der Behrens, W., Kössl, M., and Gaese, B. H. (2011). Stimulus-specific adaptation in the gerbil primary auditory thalamus is the result of a fast frequency-specific habituation and is regulated by the corticofugal system. J. Neurosci. 31, 9708–97022.
Cetas, J. S., de Venecia, R. K., and McMullen, N. T. (1999). Thalamocortical afferents of Lorente de Nó: medial geniculate axons that project to primary auditory cortex have collateral branches to layer I. Brain Res. 830, 203–208.
Cohn, T. E., Green, D. G., and Tanner, W. P. Jr. (1975). Receiver operating characteristic analysis. Application to the study of quantum fluctuation effects in optic nerve of Rana pipiens. J. Gen. Physiol. 66, 583–616.
Duque, D., Pérez-González, D., Ayala, A. Y., Palmer, A. R., and Malmierca, M. S. (2012). Topographic distribution, frequency and intensity dependence of stimulus specific adaptation in the inferior colliculus of the rat. J. Neurosci. 32, 17762–17774.
Faure, P. A., Fremouw, T., Casseday, J. H., and Covey, E. (2003). Temporal masking reveals properties of sound-evoked inhibition in duration-tuned neurons of the inferior colliculus. J. Neurosci. 23, 3052–3065.
Fishman, Y. I., Arezzo, J. C., and Steinschneider, M. (2004). Auditory stream segregation in monkey auditory cortex: effects of frequency separation, presentation rate, and tone duration. J. Acoust. Soc. Am. 116, 1656–1670.
Gordon, N., Shackleton, T. M., Palmer, A. R., and Nelken, I. (2008). Responses of neurons in the inferior colliculus to binaural disparities: insights from the use of Fisher information and mutual information. J. Neurosci. Methods 169, 391–404.
Hernández, O., Espinosa, N., Pérez-González, D., and Malmierca, M. S. (2005). The inferior colliculus of the rat: a quantitative analysis of monaural frequency response areas. Neuroscience 132, 203–217.
LeBeau, F. E., Malmierca, M. S., and Rees, A. (2001). Iontophoresis in vivo demonstrates a key role for GABA(A) and glycinergic inhibition in shaping frequency response areas in the inferior colliculus of guinea pig. J. Neurosci. 21, 7303–7312.
Loftus, W. C., Malmierca, M. S., Bishop, D. C., and Oliver, D. L. (2008). The cytoarchitecture of the inferior colliculus revisited: a common organization of the lateral cortex in rat and cat. Neuroscience 154, 196–205.
Malmierca, M. S., Blackstad, T. W., and Osen, K. K. (2011). Computer-assisted 3-D reconstructions of Golgi-impregnated neurons in the cortical regions of the inferior colliculus of rat. Hear. Res. 274, 13–26.
Malmierca, M. S., Hernández, O., Falconi, A., Lopez-Poveda, E. A., Merchan, M., and Rees, A. (2003). The commissure of the inferior colliculus shapes frequency response areas in rat: an in vivo study using reversible blockade with microinjection of kynurenic acid. Exp. Brain Res. 153, 522–529.
Malmierca, M. S., Izquierdo, M. A., Cristaudo, S., Hernández, O., Pérez-González, D., Covey, E., et al. (2008). A discontinuous tonotopic organization in the inferior colliculus of the rat. J. Neurosci. 28, 4767–4776.
Malmierca, M. S., and Ryugo, D. K. (2011). “Descending connections of auditory cortex to the midbrain and brainstem,” in The Auditory Cortex, eds J. A. Winer and C. E. Schreiner (New York, NY: Springer), 189–208.
McAlpine, D., Jiang, D., Shackleton, T. M., and Palmer, A. R. (2000). Responses of neurons in the inferior colliculus to dynamic interaural phase cues: evidence for a mechanism of binaural adaptation. J. Neurophysiol. 83, 1356–1365.
Mill, R., Coath, M., Wennekers, T., and Denham, S. L. (2011). A neurocomputational model of stimulus-specific adaptation to oddball and Markov sequences. PLoS Comput. Biol. 7:e1002117. doi: 10.1371/journal.pcbi.1002117
Patel, C. R., Redhead, C., Cervi, A. L., and Zhang, H. (2012). Neural sensitivity to novel sounds in the rat's dorsal cortex of the inferior colliculus as revealed by evoked local field potentials. Hear. Res. 286, 41–54.
Pérez-González, D., Hernández, O., Covey, E., and Malmierca, M. S. (2012). GABA(A)-mediated inhibition modulates stimulus-specific adaptation in the inferior colliculus. PLoS ONE 7:e34297. doi: 10.1371/journal.pone.0034297
Pérez-González, D., Malmierca, M. S., Moore, J. M., Hernández, O., and Covey, E. (2006). Duration selective neurons in the inferior colliculus of the rat: topographic distribution and relation of duration sensitivity to other response properties. J. Neurophysiol. 95, 823–836.
Shackleton, T. M., Skottun, B. C., Arnott, R. H., and Palmer, A. R. (2003). Interaural time difference discrimination thresholds for single neurons in the inferior colliculus of guinea pigs. J. Neurosci. 23, 716–724.
Shofner, W. P., and Dye, R. H. Jr. (1989). Statistical and receiver operating characteristic analysis of empirical spike-count distributions: quantifying the ability of cochlear nucleus units to signal intensity changes. J. Acoust. Soc. Am. 86, 2172–2184.
Skottun, B. C., Shackleton, T. M., Arnott, R. H., and Palmer, A. R. (2001). The ability of inferior colliculus neurons to signal differences in interaural delay. Proc. Natl. Acad. Sci. U.S.A. 98, 14050–14054.
Spitzer, M. W., and Semple, M. N. (1993). Responses of inferior colliculus neurons to time-varying interaural phase disparity: effects of shifting the locus of virtual motion. J. Neurophysiol. 69, 1245–1263.
Spitzer, M. W., and Semple, M. N. (1998). Transformation of binaural response properties in the ascending auditory pathway: influence of time-varying interaural phase disparity. J. Neurophysiol. 80, 3062–3076.
Stabler, S. E., Palmer, A. R., and Winter, I. M. (1996). Temporal and mean rate discharge patterns of single units in the dorsal cochlear nucleus of the anesthetized guinea pig. J. Neurophysiol. 76, 1667–1688.
von der Behrens, W., Bäuerle, P., Kossl, M., and Gaese, B. H. (2009). Correlating stimulus-specific adaptation of cortical neurons and local field potentials in the awake rat. J. Neurosci. 29, 13837–13849.
Keywords: SSA, deviant sensitivity, change detection, mismatch negativity, non-lemniscal pathway, ROC analysis
Citation: Ayala YA, Pérez-González D, Duque D, Nelken I and Malmierca MS (2013) Frequency discrimination and stimulus deviance in the inferior colliculus and cochlear nucleus. Front. Neural Circuits 6:119. doi: 10.3389/fncir.2012.00119
Received: 10 September 2012; Accepted: 19 December 2012;
Published online: 14 January 2013.
Edited by:Eric D. Young, Johns Hopkins University, USA
Reviewed by:Sarah L. Pallas, Georgia State University, USA
Edward L. Bartlett, Purdue University, USA
Copyright © 2013 Ayala, Pérez-González, Duque, Nelken and Malmierca. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in other forums, provided the original authors and source are credited and subject to any copyright notices concerning any third-party graphics etc.
*Correspondence: Manuel S. Malmierca, Laboratory for the Neurobiology of Hearing, Auditory Neurophysiology Unit (Lab 1), Institute of Neuroscience of Castilla y León, University of Salamanca, C/ Pintor Fernando Gallego 1, 37007 Salamanca, Spain. e-mail: email@example.com