JN AJP: Advances in Physiology Education
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Neurophysiol 91: 632-645, 2004. First published October 8, 2003; doi:10.1152/jn.00779.2003
0022-3077/04 $5.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
91/2/632    most recent
00779.2003v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (8)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Ingham, N. J.
Right arrow Articles by McAlpine, D.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Ingham, N. J.
Right arrow Articles by McAlpine, D.

Spike-Frequency Adaptation in the Inferior Colliculus

Neil J. Ingham and David McAlpine

Department of Physiology, University College London, Gower Street, London, WC1E 6BT, United Kingdom

Submitted 11 August 2003; accepted in final form 29 September 2003


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 ACKNOWLEDGMENTS
 REFERENCES
 
We investigated spike-frequency adaptation of neurons sensitive to interaural phase disparities (IPDs) in the inferior colliculus (IC) of urethane-anesthetized guinea pigs using a stimulus paradigm designed to exclude the influence of adaptation below the level of binaural integration. The IPD-step stimulus consists of a binaural 3,000-ms tone, in which the first 1,000 ms is held at a neuron's least favorable ("worst") IPD, adapting out monaural components, before being stepped rapidly to a neuron's most favorable ("best") IPD for 300 ms. After some variable interval (1–1,000 ms), IPD is again stepped to the best IPD for 300 ms, before being returned to a neuron's worst IPD for the remainder of the stimulus. Exponential decay functions fitted to the response to best-IPD steps revealed an average adaptation time constant of 52.9 ± 26.4 ms. Recovery from adaptation to best IPD steps showed an average time constant of 225.5 ± 210.2 ms. Recovery time constants were not correlated with adaptation time constants. During the recovery period, adaptation to a 2nd best-IPD step followed similar kinetics to adaptation during the 1st best-IPD step. The mean adaptation time constant at stimulus onset (at worst IPD) was 34.8 ± 19.7 ms, similar to the 38.4 ± 22.1 ms recorded to contralateral stimulation alone. Individual time constants after stimulus onset were correlated with each other but not with time constants during the best-IPD step. We conclude that such binaurally derived measures of adaptation reflect processes that occur above the level of exclusively monaural pathways, and subsequent to the site of primary binaural interaction.


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 ACKNOWLEDGMENTS
 REFERENCES
 
Spike-frequency adaptation, in which a neuron's response to a steady-state stimulus is not maintained at its initially high rate of spiking but instead declines over time to a lower, adapted rate is a common feature of many sensory neurons (Hille 1992Go), including those of the auditory system. The classic "primary-like" response of auditory nerve fibers is a welldescribed example of adaptation in the peripheral auditory pathway (Westerman and Smith 1984Go; Yates et al. 1985Go). Adaptation in the visual system appears to function as a gain control mechanism adjusting neural responses to maintain their responsiveness over a wide range of intensities in both the retina (Laughlin 1989Go) and the visual cortex (Ohzawa et al. 1985Go). However, specific roles for adaptation in auditory processing are yet to be established. One possible role for adaptation in the auditory system lies in determining the sensitivity of auditory neurons to the stimulus context. The context in which a sound is heard is likely to influence the way in which the auditory nervous system processes that sound, and a number of recent studies have begun to examine how sensitivity to stimulus context, including stimulus history, emerges in the ascending auditory pathway (Malone and Semple 2001Go; Malone et al. 2002Go; Sanes et al. 1998Go; Spitzer and Semple 1991Go,1998Go; Ulanovsky et al. 2003Go). One of the paradigms used in these studies is interaural phase modulation (IPM). IPM is a binaural stimulus in which the interaural phase disparity (IPD) is modulated at different rates, to different depths, or around different values of interaural phase, to produce the percept of auditory motion. Many neurons in the inferior colliculus (IC), the major auditory nucleus in the midbrain, are sensitive to the context in which IPDs are presented as part of the IPM stimulus. Their responses may be enhanced or suppressed depending on the direction, speed, or center-locus of the virtual motion cues (McAlpine et al. 2000Go; Spitzer and Semple 1991Go, 1995Go, 1998Go). In contrast, responses in the medial superior olive (MSO) appear to reflect the instantaneous interaural disparity of a stimulus (Spitzer and Semple 1995Go; Yin and Chan 1990Go), suggesting that the main source for the creation of sensitivity to dynamic binaural cues in the IC lies within the IC itself (Borisyuk et al. 2002Go).

Subsequent investigations have shown IC neurons to be sensitive to the context in which acoustic cues other than IPDs, such as interaural level differences (Sanes et al. 1998Go) or frequency glides (Malone and Semple 2001Go), are presented, and that the influence of dynamic context on neural responses is even more pronounced in auditory cortex than in the IC (Malone et al. 2002Go). Together, these data suggest a hierarchy of context sensitivity in the ascending auditory nervous system, with MSO neurons showing the least sensitivity (Spitzer and Semple 1995Go, 1998Go), cortical neurons showing the most sensitivity (Malone et al. 2002Go), and IC neurons showing intermediate sensitivity (Malone and Semple 2001Go; McAlpine et al. 2000Go; Spitzer and Semple 1991Go, 1995Go, 1998Go).

Although the phenomenology of stimulus-context dependency in the auditory pathway is increasingly becoming established, the mechanisms that contribute to this sensitivity remain to be determined. In terms of dynamic interaural cues, it was originally suggested that binaurally sensitive inhibition, possibly derived from the dorsal nucleus of the lateral lemniscus, was responsible for the appearance of context sensitivity in the IC (Spitzer and Semple 1995Go, 1998Go). However, McAlpine et al. (2000Go) argued that an adaptation process residing above the level of binaural integration was sufficient to account for the data. In support of this, McAlpine and Palmer (2002Go) demonstrated that blocking GABA-ergic inhibition in IC neurons enhanced differential sensitivity to dynamic cues of IPM, whereas adding tonic inhibition reduced such sensitivity. The explanation posited for this behavior was that increases in the discharge rate that accompanied blockade of inhibition moved neurons closer to an adapting state, whereas decreases in discharge rate that accompanied addition of tonic inhibition moved neurons away from the adapting state. Recent modeling studies explain the responses of IC neurons to a range of static and dynamic interaural timing cues using a combination of ITD-sensitive and ITD-insensitive GABA-ergic inhibition, postinhibitory rebound, as well as an adaptation mechanism intrinsic to the IC in their calculations (Cai et al. 1998aGo,bGo).

As the 1st stage in assessing the specific contribution of IC adaptation mechanisms to the coding of dynamic stimulus cues, we have developed a novel stimulus paradigm that enables adaptation characteristics of neurons at or above the level of binaural integration to be measured free of the influences of adapting monaural components in the pathway. The stimulus consists of a 3,000-ms binaural tone, the first 1,000 ms of which is presented at a neuron's least favorable (worst) IPD. During this time, neural components in the auditory pathway that are insensitive to interaural phase cues will be responding, and adapting to their steady-state level. After 1,000 ms, the interaural phase of the stimulus is rapidly switched to the neuron's most favorable (best) IPD, allowing the time course of neural adaptation to be assessed free from the influences of adapting monaural components. This is possible because neurons sensitive to IPDs are essentially "gated OFF" by the worst IPD epoch, becoming active only when "gated ON" on stepping to the best IPD. IPD is the only stimulus parameter of which we are aware whose influence on neural activity is delayed until such a late stage in the auditory pathway. Stepping stimulus parameters such as sound frequency and level (including interaural level) in this manner would influence the output of one or both cochleae and their target neurons, removing the possibility of examining binaural adaptation free from the influence of changing/adapting monaural components. Our data indicate this to be an important consideration when only the binaural configuration of a stimulus is altered. The adaptation profiles obtained by this method most probably reflect mechanisms intrinsic to the IC neurons themselves (Borisyuk et al. 2002Go), and not the properties of their input neurons in the MSO, given that MSO neurons show little evidence of the differential sensitivity to dynamic cues (Spitzer and Semple 1998Go). Although only interaural phase cues were used to assess adaptation in the IC, the data have more general application because any stimulus parameter that evokes activity in IC will also evoke some form of adapting response to it. These data will thus be useful in constructing neural models that seek to determine the contribution of different auditory centers to analysis of the auditory scene. Elements of these data were previously published in abstract form (Ingham and McAlpine 2001Go).


    METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 ACKNOWLEDGMENTS
 REFERENCES
 
Surgical preparation

All experiments were carried out in accordance with the guidelines of the UK Home Office, under the control of the Animals (Scientific Procedures) Act 1986. Pigmented guinea pigs (Cavia porcella) were anesthetized with urethane (1 g/kg; 25% solution) in preparation for surgical procedures and experimental recordings. Analgesia was induced, and supplemented as required, using 0.1 ml intramuscular injections of Hypnorm (fentanyl citrate/fluanisone). Local anesthesia at surgical sites was achieved using subcutaneous injections of 2% lignocaine. A tracheal cannula was inserted, and body temperature was maintained at 37°C using a thermostatically controlled heating blanket and rectal probe. Animals were mounted in a modified stereotaxic frame situated inside a sound-attenuating booth (IAC, Winchester, UK). Hollow ear speculae allowed the insertion of custommade ear phones and probe-tube microphones to form a sealed pressure-field sound delivery system. Pressure equalization of the middle ear was achieved by sealing high acoustic-impedance cannulae into each bulla. A craniotomy was performed, extending 2–3 mm rostral and caudal of the interaural axis, and 1–4 mm lateral from the midline on the right side. The dura overlying the cortex was removed, allowing microelectrode access through the cortex to the right inferior colliculus, and the cranium sealed with 2% agar.

Single-neuron recordings

Neural recordings were made using parylene-coated tungsten microelectrodes (1–5 M{Omega} impedance), mounted on a piezo-electric stepper motor, and advanced stereotaxically through the inferior colliculus in a dorsoventral fashion. Electrical activity from the microelectrode was filtered (300 Hz –3 kHz) and amplified (variable gain) using an AC differential amplifier and a PC1 spike conditioner (Tucker Davis Technologies, Alachua, FL). Single neurons were isolated using variable frequency and intensity diotic tones. Single spikes were discriminated from background noise using an SD1 spike discriminator (TDT) and time stamped to 1-µs accuracy using an ET1 event timer (TDT), linked to the computer. Single-neuron isolation was confirmed by the consistency of the discriminated spike waveform displayed on a digital oscilloscope.

Stimulus presentation and data analysis

Acoustic stimuli were produced and presented under computer control. Digitally generated dichotic stimuli (TDT, AP2 digital signal processor; at 100- or 48-kHz sampling rate) were converted to analog signals (TDT, DA3-2). The signals were filtered (TDT, FT6; fc = 40 kHz), and attenuated (TDT, PA4), before amplification and delivered to the ears by RadioShack 40-1377 (Fort Worth, TX), or Beyerdynamic DT-48 (Burgess Hill, UK) loudspeaker units fitted with brass tube attachments sealed into the hollow ear speculae supporting the animal. The sound field inside the sealed system was sampled using Knowles FG3452 microphones by a probe tube inserted to within a few millimeters of each tympanic membrane. Probe microphones were previously calibrated against a Bruel and Kjaer Type 4136 1/8-in. microphone, and the sound systems for each ear were flat to within ±5 dB from 50 to 2,000 Hz and matched to within ±5 dB for this range.

Once isolated, a neuron's characteristic frequency (CF) was estimated audiovisually, and confirmed by generating a frequency-versus-level response area, extending 2 octaves above and 4 octaves below the estimated CF, using diotic 50-ms tones covering a range of 60–100 dB, presented at 5 Hz in pseudo-random order. Binaural sensitivity to interaural phase disparity (IPD) cues was assessed using binaural beats of 3-s duration, with 1-Hz difference in the carrier frequency presented to the 2 ears. The initial and final 500-ms periods of the response were omitted from analysis, and the middle 2,000 ms were averaged and plotted on-line with respect to the IPD of the stimulus to generate period histograms. Using the methodology of Goldberg and Brown (1969Go), the vector strength (R) of the response was calculated from the period histogram. The average best phase was calculated as a vector average of the response magnitudes at each point in the cyclic IPD phase histograms. These circular distributions were assessed using the Rayleigh test for the significance of the mean best phase. Where the Rayleigh coefficient < 13.816, the neuron was considered as not significantly modulated by the IPD of the binaural beat, and such units were excluded from further analyses. Best and worst IPD were determined from the IPD period histograms. Because many period histograms were roughly symmetric, the best IPD equated to the mean best phase described above. However, for some neurons, where the IPD-sensitive response produced asymmetric functions, the best IPD was taken to be that evoking the peak response. The worst IPD was generally determined as 0.5 cycles removed from the best phase, or where the phase response dropped to zero at less-favorable IPDs. In other cases, for example, where the plot was asymmetric, or where there was a baseline response, the worst IPD was taken as the IPD evoking the lowest spike activity.

Best and worst IPD values were used to generate IPD-step stimuli designed to determine the time course of adaptation and recovery from adaptation. Digital sound files 3,000 ms in duration ("wav" format, 48-kHz sample rate), were generated using a custom script in Matlab. Two digital sound files were constructed, one for presentation to either ear. Parameters such as carrier frequency, best IPD, worst IPD, and the rate at which the IPD was stepped between these values were under user control. The total phase shift at each ear was minimized for each phase step by advancing the phase by half the required amount (the difference between best and worst IPD in this case) in one ear, and regressing the phase by the same amount in the other ear. This "IPD-step" stimulus is illustrated in Fig. 1, AE, along with a stylized neural response (Fig. 1, F and G). After stimulus onset at worst IPD, neurons usually responded with a train of discharges that decayed exponentially over time from a maximum (Omax) to reach a steady state (Oss) by the end of the first 1,000 ms of the stimulus. At this point, the interaural phase was stepped rapidly from worst IPD to best IPD, and the neuron responded with a new maximal discharge rate (Bmax), which once more decayed (Bd) over the 300-ms best-IPD epoch to reach a new steady-state level (Bss) above that of Oss. At the end of the best-IPD epoch, the IPD was stepped back to the worst IPD and the neuron responded with a reduced discharge rate, equal to that of the initial steady-state level, Oss. After some predetermined time [interstep interval (ISI)], the IPD was again stepped from worst IPD to best IPD for a further 300 ms. The response during this 2nd best-IPD epoch depends on the duration of the interval between the 1st and 2nd IPD steps. For longer intervals, a similar response pattern was produced to that seen during the 1st best-IPD epoch; that is, the discharge rate was initially high (Rmax) and decayed over time to reach a steady state (Rss). For shorter intervals, the adapting region of the response (Rd) was reduced, such that the initial peak of the recovery response (Rmax) was lower (see Fig. 1G, open circles), or not apparent, followed by a reduction in discharge rate back to the same plateau level seen for the 1st best-IPD epoch (Bss). Thus it was not always possible to fit the response to the 2nd step to best IPD with an exponential function.



View larger version (25K):
[in this window]
[in a new window]
 
FIG. 1. A and B: dichotic waveforms for the interaural phase disparity (IPD)–step stimulus at the time of the 1st step to best IPD for a 389-Hz tone. A: between 1,000 ms and 1,006.4 ms the IPD is stepped from worst IPD (0.61 cycles in this example) to best IPD (0.11 cycles). B: between 1293.6 ms and 1300 ms the IPD is stepped back from best IPD to worst IPD. C and D: IPD transitions for A and B, respectively. C: instantaneous phase of the tone in the leading ear was delayed (black line), and the lagging ear waveform was advanced (gray line), in phase terms, to produce the step to best IPD. D: this process was reversed for the step back to worst IPD. IPD was always stepped between configurations by the shortest route (in IPD terms). This process was repeated after some interval to produce the 2nd step to best IPD. E: changes in IPD for an IPD-step stimulus with two 300-ms steps to best IPD separated by 1,000 ms. F: stylized neuronal response to the IPD-step stimulus, illustrating exponential adaptation profiles. See METHODS for further details. G: exponential decay functions illustrating how recovery from adaptation was assessed. Circles indicate the peak discharge rate of the exponential decay fit at each recovery interval. Gray line indicates the exponential rise-to-max function fitted to peak discharges, used to determine the {tau}rec. See METHODS for further details.

 
Poststimulus time histograms (PSTHs), adjusted for the latency of individual neurons, were constructed from the response to 10–50 repeats of the IPD-step stimulus, and plotted in 10-ms bins for off-line analysis. The "onset," "monaural," and "binaural" adaptation components of the PSTHs were extracted and subject to a 3-point smoothing function. A single exponential decay function was fitted to the data, and adaptation time constants for onset adaptation ({tau}on), monaural adaptation ({tau}mon), and binaural adaptation ({tau}bin) were extracted (see Eq. 1)

(1)
where ft is the discharge rate at time t (spikes/s), fss is the discharge rate at plateau, f0 is the discharge rate when t = 0, t is the time (ms) from onset/switch to best IPD (plus unit latency), and {tau} is the time constant of adaptation (ms).

For onset and binaural adaptation profiles, where responses to multiple stimuli with different recovery intervals were recorded, equivalent data bins from each response were averaged before being smoothed and subject to curve fitting (this was not possible for the monaural adaptation because the response to only the longest recovery time was recorded). The curve fits were subject to statistical analysis using ANOVA. Neurons whose ANOVA test results produced values of P > 0.05 were classified as "nonadapting"; the exponential decay function did not produce a significant fit to the physiological data. Curves fitted to the data for the various adaptation profiles were normalized to enable comparison of the curves and rates of decay for individual neurons. The extent to which a neuron's response adapted was assessed by calculating an adaptation ratio, defined as the ratio of the discharge rate at steady state to that at the peak of the adapting profile (e.g., binaural adaptation ratio = Bss/Bmax).

Recovery from adaptation was assessed in 2 ways. Neural responses to the 1st best-IPD epoch were analyzed as an average of the presentations made with each different interval period. Responses evoked during the 2nd epoch at best IPD were fitted with exponential decay functions, as outlined above (Eq. 1) and described by a 2nd adaptation time constant ({tau}bin2). Because responses to the 2nd step at best IPD were based on an average of fewer stimulus presentations compared with the 1st step to best IPD, they were more variable, resulting in fewer neurons producing sufficient data for detailed analysis of the recovery from adaptation. Recovery time constants ({tau}rec) were determined for 37 of the above neurons. Here, peak discharge rates from the function fitted to the 2nd step to best IPD were plotted as a function of the interval between the 2 best-IPD epochs (see Fig. 1G). With longer recovery intervals, the peak discharge rate of the 2nd step to best IPD increased toward that of the 1st step to best IPD. These data were fitted with exponential rise-to-maximum curves (see Eq. 2). The steady-state discharge rate evoked by the 2nd step to best IPD was similar, for all ISIs, to the discharge rate evoked during the 1st step to best IPD

(2)
where ft is the discharge rate at time t (spikes/s), f0 is the baseline discharge rate, t is the recovery interval time (ms), and {tau} is the time constant of recovery {tau}rec (ms).


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 ACKNOWLEDGMENTS
 REFERENCES
 
Adaptation profiles of 120 IPD-sensitive single neurons, isolated from the right IC of 45 guinea pigs, were examined using the IPD-step stimulus. Many of these animals were involved in other ongoing studies of binaural processing, and yielded data using other experimental paradigms.

Single-neuron responses to IPD steps

The IPD-step stimulus was designed to adapt out the influences of monaural components of the ascending auditory pathway up to the level of the IC, so that the time course of adaptation when only the binaural configuration of a sound is changed could be assessed. This method requires neurons to be IPD sensitive. The top row of Fig. 2 shows responses for 4 IPD-sensitive IC neurons to 3 s of a 1-Hz binaural beat stimulus. In each case, the discharge rate (black histogram bars) was modulated with each cycle of the beat (gray lines), indicating neurons to be IPD sensitive. These data, plotted in the 2nd row of Fig. 2 in the form of period histograms, were used to determine a neuron's best IPD [the IPD that evoked the highest discharge rate (filled circle in each period histogram)] and worst IPD [the IPD that evoked the lowest discharge rate (open circle in each period histogram)]. The remaining rows of Fig. 2 show the responses of each neuron to the IPD-step stimulus. Each neuron responded to the onset at worst IPD with an initially high discharge rate that adapted to a much lower rate (to zero in some cases; e.g., Fig. 2, column B) during the subsequent 1,000 ms. At this point, the IPD was stepped rapidly from worst IPD to best IPD, and held at best IPD for 300 ms, before being stepped rapidly back to worst IPD. This initial step to best IPD evoked a 2nd, initially high, response that adapted to a lower steady-state firing rate over the 300-ms duration of the step. This steady-state response was higher than the steady-state firing rate during the first 1,000 ms of the stimulus, even for those neurons for which the firing rate evoked during the first 1,000 ms remained relatively high after stimulus onset (e.g., Fig. 2, column D). After the step back to worst IPD, firing rates fell back to, and sometimes below, firing rates evoked before the step to best IPD.



View larger version (46K):
[in this window]
[in a new window]
 
FIG. 2. AD: responses of 4 IC neurons to binaural beats and to IPD-step stimuli for different intervals between the 2 steps to best IPD. Row 1: PSTHs to binaural beats. IPD is plotted in gray, with respect to left-most ordinate. Row 2: binaural beat responses plotted as period histograms. Open circles indicate worst IPD; filled circles indicate best IPD. Rows 39: responses to IPD steps for interstep intervals (ISIs) of (from top to bottom) 1,000, 300, 100, 30, 10, 3, and 1 ms. Bottom row: responses to monaural (left ear only) presentations of the left channel wav-file used in the IPD-step stimuli.

 
After the step back to worst IPD, a 2nd step to best IPD was initiated after some interval (1,000–1 ms; see Fig. 2, rows 39) for a further 300 ms, before being stepped rapidly back to worst IPD for the remainder of the stimulus. The response to this 2nd IPD step was more variable than the response to the 1st IPD step. Specifically, the initial portion of the response to the 2nd step to best IPD was highly dependent on the interval between the 2 steps. For longer intervals, the response to the 2nd step resembled that to the 1st step, with an initial onset response, followed by decay to a steady state. For shorter intervals, however, the peak of the 2nd best-IPD step response was attenuated, or absent altogether.

Measuring spike-frequency adaptation after a step to best IPD

To assess the time course of adaptation after the step to best IPD, single exponential decay functions (Eq. 1) were fitted to the first 300-ms period of the PSTH over which the interaural phase was stepped from worst to best IPD (adjusted for neural latency, see METHODS). The majority of neurons (111/120, 92.5%) showed an adapting profile to this 1st step (Fig. 3, A and B) that was well described by a single exponential decay function (ANOVA, P < 0.0001). Fitted functions to the raw (Fig. 3A) or normalized (Fig. 3B) data indicate mean and median adaptation time constants for this step to best IPD ({tau}bin) of 61.4 and 46.9 ms, respectively. Time constants were not normally distributed [Kolmogorov–Smirnov (K-S); P < 0.0001] but, rather, were skewed toward lower values of {tau}bin (Fig. 4A). Using boxplot analysis, 4 neurons (dashed lines in Fig. 3B and open bars in Fig. 4A) were classified as outliers. Removing these outliers from the analysis reduced the mean {tau}bin from 61.4 to 52.9 ms and produced a dramatic reduction in the variance of the population data (see Table 1). These outliers were excluded from further analysis. Responses of the remaining 9/120 (7.5%) neurons to the initial step to best IPD were not well described by exponential functions (ANOVA, P > 0.05). These neurons tended to show a "build-up" response to the step to best IPD, or a stepped response that did not adapt (data not shown).



View larger version (43K):
[in this window]
[in a new window]
 
FIG. 3. A and B: exponential decay functions fitted to the 1st step to best IPD; A shows raw fits and B shows normalized fits. C and D: exponential decay functions fitted to the onset of the worst IPD stimulus. E and F: exponential decay functions fitted to the onset of the monaural stimulus. In the normalized functions of B, D, and F, 1.0 indicates the peak response and 0.0 indicates the steady-state response. Dashed lines in each panel indicate functions for which time constants were classified as outliers by boxplot analyses. See RESULTS for further details.

 



View larger version (20K):
[in this window]
[in a new window]
 
FIG. 4. A: distribution of binaural adaptation time constants ({tau}bin). Outlying values, determined by boxplot analysis, are indicated by open bars and are not included in the population mean (±SD) plotted by the circle. B: distribution of adaptation ratios for the 1st step to best IPD. Black circle indicates the mean adaptation ratio (±SD) for the population. Gray circles indicate the mean adaptation ratio (±SD) for the 2nd step to best IPD at each recovery interval (indicated on the right ordinate).

 


View this table:
[in this window]
[in a new window]
 
TABLE 1. Adaptation time constants

 
The degree to which IC neurons adapted during the 1st step to best IPD was assessed by calculating the ratio of the steady-state discharge rate to the peak discharge rate, from the exponential functions fitted to the decaying spike-frequency profiles to produce an adaptation ratio. Adaptation ratios during the 1st step to best IPD were normally distributed (K-S test, P > 0.05), with a mean of 0.42 ± 0.18 (mean ± SD; median = 0.41; see Fig. 4B), indicating that, on average, 58% of the peak discharge rate after a step to best IPD is adapted out. Adaptation ratios ranged from 0 to 0.75, indicating that between 100 and 25% of the peak discharge is adapted out for the population of 111 neurons.

Recovery from adaptation

A neuron's time course of recovery from adaptation was assessed from the response to the 2nd step to best IPD (see columns in Fig. 2). The magnitude of the initial response to this 2nd step was highly dependent on the duration of the interval between the 2 steps, and by fitting exponential rise-to-max functions to the peak discharge rates over the range of ISIs (Fig. 5; and see Fig. 1G), we were able to determine the rate at which neurons recovered from the effects of adaptation to the 1st step. Of the 48 neurons for which sufficient data points were available, 37 showed a significant exponential fit (ANOVA, P < 0.1), and recovery time constants ({tau}rec) derived from these fits (Fig. 6, A and B) indicate that recovery from binaural adaptation was a much slower process than adaptation itself (mean {tau}rec = 225.5 ± 210.6 ms; see distribution of {tau}rec in Fig. 6C). Relative to {tau}bin, {tau}rec was more variable, and there was no correlation between {tau}bin and {tau}rec for individual neurons.



View larger version (27K):
[in this window]
[in a new window]
 
FIG. 5. Composite responses of 2 neurons illustrating the methods used to assess adaptation during recovery and recovery from adaptation. Light gray bars indicate a composite poststimulus time histogram (PSTH) (10-ms binwidth, adjusted for neuronal latency) of responses recorded to stimuli with ISIs ranging from 3 ms to 1,000 ms. Dark gray bars indicate the average of the 7 responses recorded to the 1st step to best IPD, which was fitted with an exponential decay curve (shown) to calculate binaural adaptation ({tau}bin). Subsequent decaying curves indicate exponential decay fits to the 2nd step to best IPD (each producing a 2nd binaural adaptation time constant, {tau}bin2). Open circles indicate the peak of these exponential fits, used to determine the recovery from adaptation. Exponential rise functions were fitted to these points (shown) to determine the time constant for recovery ({tau}rec).

 



View larger version (30K):
[in this window]
[in a new window]
 
FIG. 6. A: time course of recovery from adaptation, obtained by fitting an exponential rise-to-max function to the peak discharge rate of the 2nd step to best IPD for the different ISIs. B: normalized recovery functions. Inset shows the same curves, with the abscissa expanded to the same time scale (300 ms) as for the adaptation curves. C: distribution of recovery time constants. Circle indicates the population mean recovery time constant ({tau}rec; ±SD).

 
Adaptation during recovery from adaptation

Under natural listening conditions it is likely that a neuron is subject to continuously varying states of adaptation and recovery from adaptation, mirroring the time-varying nature of most natural sounds. To assess a neuron's adaptation during different adaptation states, exponential decays were fitted to the response profiles recorded to the 2nd step to best IPD and a 2nd binaural adaptation time constant ({tau}bin2) was calculated (see Fig. 5 and Fig. 1G) in a similar manner to the step to best IPD. Adaptation time constants measured for the 2nd step to best IPD were similar to those measured for the 1st step to best IPD for individual neurons. For recovery intervals for which sufficient data were available, {tau}bin2 was generally in accordance with {tau}bin. Figure 7 compares {tau}bin2 with {tau}bin for ISIs of 1,000 ms (7A), 300 ms (7B), 100 ms (7C), and 30 ms (7D). Linear regression fits to the data for each interval, and for the pooled (Fig. 7E) and average (7F) data indicate a clear relationship between the 2 adaptation time constants. These observations were confirmed by plotting the averaged {tau}bin2 against its recovery interval (Fig. 7G). There was no systematic change in adaptation time constant recorded across the range of recovery intervals tested. Thus whenever IC neurons adapt, they tend to adapt at the same rate irrespective of their current state of adaptation.



View larger version (21K):
[in this window]
[in a new window]
 
FIG. 7. AD: comparison of {tau}bin2 and {tau}bin for ISIs of (A) 1,000 ms, (B) 300 ms, (C) 100 ms, and (D) 30 ms. Solid lines illustrate linear regression. E: pooled data across all ISIs, including those not shown in AD). F: data averaged across all ISIs. Regression line slope was 1.148 and intercept was –7.090 (r2 = 0.633). G: averaged adaptation time constants (± SD), for the 1st step to best IPD and for subsequent steps after different recovery intervals, plotted against the recovery interval (ms).

 
Responses to monaural and diotic presentation: responses at the phase transition

The tones presented to each ear in the IPD-step stimulus contain a rapid phase shift during the transition from worst to best IPD. Even with monaural presentation, there is an audible "blip" when this transition occurs (4 times in the 3-s IPD-step stimulus), a percept that likely arises as a result of the simultaneous phase shift across many peripheral auditory filters—the rapid phase shift is equivalent to a FM. Because we record the activity of single neurons, presumably reflecting the output of very few auditory filters, this "blip" is unlikely to influence greatly, or at all, neural discharge rates. However, to discount any influence of the monaural phase shifts per se on the response of individual neurons to the IPD-step stimulus, we recorded responses to monaural and diotic stimulation that provide for the same rapid phase transitions as the dichotic stimuli, but not for the change in interaural phase cues.

PSTHs for two representative neurons, in response to contralateral monaural stimulation (Fig. 8, A and B) or diotic stimulation in which the phase transition in each ear was in the same direction (Fig. 8, C and D), indicate no influence of the phase transition on evoked discharge rates (cf. with large increases in discharge rate to dichotic transitions in Fig. 2). Analysis of the discharge rate during the 5- to 10-ms period of the phase transition for 23 neurons (monaural or diotic) revealed little or no influence of these phase shifts on evoked discharge rates (Fig. 8, E and F). No systematic increase or decrease in activity could be attributed to the time of the phase transition. This finding is important because it demonstrates that IC neurons are insensitive to the rapid monaural phase shifts that generate the steps between worst and best IPD, and that any changes in discharge rate can be attributed to the gating of activity through the binaural coincidence detectors in the MSO that project to these neurons.



View larger version (32K):
[in this window]
[in a new window]
 
FIG. 8. A and B: responses of 2 neurons to monaural (left ear only) presentation of the IPD-step stimulus. C and D: responses of neurons in A and B to diotic presentation of the left channel stimulus of the IPD-step stimulus. Only response regions surrounding 1,000 ms, the 1st phase shift period, are shown. Period of phase transition is indicated by the black histogram bars. E and F: averaged discharge rates in 10-ms periods beginning at 30, 20, and 10 ms before the phase shift, during the 10 ms (maximum) of the phase shift, and during 10-ms periods beginning at 10, 20, and 30-ms periods after the phase shift are shown for 60 neurons to monaural stimulation (E), and 8 neurons to diotic stimulation (F). Open symbols represent the values for individual neurons. Filled symbols represent the mean (±SD) at each time period. Black symbols indicate the averaged discharge rate during the phase shift. Gray symbols indicate average discharge rates at other periods. No systematic shift in discharge rate during the 10-ms phase transition period can be distinguished in either case.

 
Spike frequency adaptation at stimulus onset

The rationale for the IPD-step stimulus is that it enables adaptation time constants to be measured above the level of binaural integration free from the influence of adaptation mechanism at levels lower than that of binaural integration. This is possible because holding the stimulus at worst IPD effectively gates off the output of binaural neurons in the MSO, evoking little or no activity in the IPD-sensitive IC neurons to which they project. When the stimulus is stepped to best IPD after 1,000-ms components of the binaural pathway below binaural integration in the MSO, and any monaural inputs to the IC, should be fully adapted. Any subsequent adaptation to changes in the binaural configuration therefore reflects adaptation of binaurally sensitive components only. Note that although binaural adaptation could reflect processing in the MSO, which is simply projected to the IC, indirect evidence, such as the relative insensitivity of MSO neurons compared with IC neurons to different directions of motion, for example, argues against this (see DISCUSSION). Here, we compare responses at stimulus onset to those at the step to best IPD.

Most neurons showed adapting profiles to the onset of monaural stimulation, with discharge rates decaying from a peak rate to a lower, steady-state rate. In some cases, monaural steady-state discharge rates were comparable to those at worst IPD (Fig. 2, columns A and B), although, consistent with previous reports (e.g., Kuwada and Yin 1983Go; McAlpine et al. 1996Go; Yin and Kuwada 1983Go), discharge rates evoked by monaural stimulation could be higher (Fig. 2, column C) or lower (Fig. 2, column D) than discharge rates evoked at worst IPD. Similar responses were also observed using diotic stimulus presentations, where the left waveform was presented simultaneously to both ears (data not shown).

To assess adaptation at stimulus onset, responses for the IPD-step stimulus, and responses to monaural presentation of the left-channel wav file, were fitted with single-exponential decay functions (Eq. 1) for the first 300 ms after stimulus onset. Decay functions fitted to the onset at worst IPD ({tau}on) are shown in Fig. 3, C and D (n = 96), and gave mean and median adaptation time constants of 48.1 and 32.0 ms, respectively. Decay functions fitted to the monaural onset ({tau}mon) are shown in Fig. 3, E and F (n = 52) and gave mean and median time constants of 51.7 and 34.5 ms, respectively. Neurons with time constants above 133 ms for the worst IPD at onset condition (n = 7, dashed lines in Fig. 4D) and neurons with time constants above 136 ms (n = 6, dashed lines in Fig. 3F) for the monaural onset condition were classified by boxplot analysis as outliers, and were excluded from subsequent analyses.

There was substantial variability in the peak discharge rate evoked at stimulus onset for both the monaural and onset at worst-IPD conditions. However, the steady-state firing rates to the monaural onset were less variable than the steady-state firing rates to the worst-IPD onset. Visual inspection indicates that neurons adapted more rapidly at stimulus onset, for both the worst-IPD onset and monaural onset conditions, than at the step to best IPD (cf. Fig. 4, B, D, and F), and this is confirmed by the distribution of time constants in Fig. 9, A and B, which also indicate a significant skewing of the distributions for {tau}mon and {tau}on (K-S, P < 0.0001). However, this was less pronounced for adaptation at the step to best IPD (Fig. 9, A and B; open bars), with relatively more neurons showing longer adaptation time constants than for onset at worst IPD (Fig. 9A; closed bars) or monaural onset (Fig. 9B; closed bars).



View larger version (14K):
[in this window]
[in a new window]
 
FIG. 9. A and B: distribution of adaptation time constants for (A) worst IPD at onset ({tau}on) and (B) monaural onset ({tau}mon), respectively. Outlying values are indicated by open bars. Open gray bars overlying these data indicate the distribution of binaural adaptation ({tau}bin) for comparison (replotted from Fig. 4). None of these data sets conforms to a normal distribution [Kolmogorov–Smirnov (K-S) test, P < 0.05)]. CE: correlations, using Kendall's tau-b nonparametric tests, between (C) {tau}bin and {tau}mon; Kendall's tau-b = 0.188 (P = 0.070), (D) {tau}bin and {tau}on; Kendall's tau-b = 0.104 (P = 0.156), and (E) {tau}mon and {tau}on; Kendall's tau-b = 0.429 (P < 0.001**). Time constant values classified as outliers were excluded from these analyses.

 
Adaptation time constants for steps to best IPD were not significantly correlated (Kendall's tau-b nonparametric test for correlation) with either of the onset adaptation time constants (Fig. 9, C and D). Time constants of the adaptation profiles for the monaural onset {tau}mon and the onset at worst IPD {tau}on were significantly correlated (Fig. 9E; P < 0.001). However, there was no appreciable correlation between measurements of the recovery time constant {tau}rec, and {tau}mon or {tau}on for the same neurons (r2 < 0.01).

For the majority of neurons, responses at onset at worst IPD and monaural onset were more strongly adapted than at the step to best IPD. The steady-state rate was just 12 ± 11% (Fig. 10A; filled bars) and 18 ± 13% (Fig. 10B; filled bars), respectively, of the initial discharge rate evoked at stimulus onset, compared with 42 ± 18% for the step to best IPD (Fig. 10, A and B; open bars, reproduced from Fig. 4B). This trend was also evident for individual neurons (Fig. 10, CE). Where comparison of the adaptation ratio for monaural and onset worst IPD was possible (n = 49), adaptation to the onset worst-IPD stimulus was usually more pronounced [i.e., had a lower adaptation ratio (Fig. 10E; 35/49 cases)] than adaptation at monaural onset. The remaining 14 neurons showed equal, or more pronounced, adaptation to the onset of the monaural stimulus (Fig. 10E; data points on or below the line of equivalence). Just 5 of 96 neurons demonstrated more pronounced adaptation to the 1st step to best IPD compared with the onset at worst IPD (Fig. 10D), and just 2 of 51 neurons demonstrated more pronounced adaptation to the 1st step to best IPD compared with the monaural onset (Fig. 10C); that is, the depth of adaptation was normally more pronounced at stimulus onset than at the step to best IPD after 1,000 ms of a worst-IPD stimulus.



View larger version (15K):
[in this window]
[in a new window]
 
FIG. 10. A and B: distribution of adaptation ratios for (A) ARons and (B) ARmon. Open gray bars plot the distribution ARbin for comparison. Filled circles indicate average adaptation ratio. Open gray circles indicate the mean adaptation ratio (±SDs). Adaptation ratios for the step to best IPD and the onset monaural condition were normally distributed (K-S test, P > 0.05). However, adaptation ratios for the onset worst-IPD adaptation were not normally distributed (K-S test, P < 0.05). CE: correlations, for individual neurons, between ARbin and ARmon (C); ARbin and ARons (D); and ARmon and ARons (E). Solid lines indicate equality.

 
One question this raises is whether these differences in the adaptation profile are attributable to the time within the stimulus at which adaptation was assessed, or the IPD. For the majority of neurons, binaural adaptation at onset was measured only for stimuli in which onset coincided with worst IPD. We did not routinely record responses to IPD steps in which stimulus onset coincided with best IPD. However, for 9 neurons for which these data were obtained, adaptation time constants to best IPD at stimulus onset were not correlated with adaptation time constants obtained to the 1st step to best IPD. Further, these neurons adapted to best IPD at onset to a greater extent (16 ± 11%) than they did to the 1st step to best IPD (35 ± 9% for the same 9 neurons; cf. 42 ± 18% for the population). This deeper adaptation to best IPD at stimulus onset was consistent with the adaptation of the same neurons to worst IPD at onset (11 ± 8%; cf. 12 ± 11% for the population) and monaural onset (16 ± 10%; cf. 18 ± 13% for the population). This suggests that the factor determining the more deeply adapting response profile at stimulus onset is related to a neuron's responses to stimulus onset rather than to the IPD at which the onset is presented. Thus adaptation of monaural pathways, including those below the level of binaural integration, will influence binaural adaptation differently depending on when the step in IPD occurs. If the step to best IPD occurs at, or soon after, stimulus onset, the adaptation of monaural components will contaminate the adaptation profile recorded to the step to best IPD. As the time between stimulus onset and the step to best IPD increases, however, the influence of monaural adaptation wanes. Thus a neuron's adaptation to a binaural stimulus is influenced by the adaptation state of monaural components in the auditory pathway.


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 ACKNOWLEDGMENTS
 REFERENCES
 
The main aim of this study was to obtain a quantitative assessment of the temporal properties of spike-frequency adaptation in IPD-sensitive neurons in the inferior colliculus. The IPD-step stimulus allowed adaptation to be measured free from the influence of adaptation in monaural pathways up to the level of the IC; neural responses sensitive to IPD were effectively gated ON and OFF by altering the interaural phase configuration of the stimulus. Although IPD-sensitive IC neurons are not the primary effectors of IPD sensitivity in the central auditory system, they do receive input directly from the MSO neurons responsible for the generation of IPD sensitivity. It is important to note the continuous nature of the acoustic stimulation throughout the entire 3-s duration of the stimulus. There is no "OFF period" when the stimulus switches from worst-IPD to best-IPD configurations. The only stimulus feature to change during these switches is the phase of the tone in either ear; the overall stimulus level remains constant. This means that the hair cells and neurons in the auditory nerve and cochlear nucleus, and inputs to primary binaural neurons, receive constant stimulation. The perceptible "blip" in the signal when the phase step occurs presumably arises because the rapid phase change, similar to a brief frequency glide, and activates multiple auditory filters in the cochlea. However, this appears to have no influence on the responses of single IC neurons (see Fig. 8), and thus is unlikely to have influenced our measurement of adaptation time constants.

Distinguishing central adaptation from peripheral adaptation

No previous study of which we are aware has been able to distinguish between peripheral (in this context, monaural) adaptation and binaural adaptation time constants, and the ability to do so here is an important outcome of the current study. Using an adaptor-probe stimulus, Finlayson and Adam (1997Go) found superior olive neurons to have a binaural adaptation time constant of about 20 ms, considerably shorter than the 51.7 ms observed here, and for adaptation to monaural stimulation to be similar to that for binaural stimulation. Recovery from adaptation was also faster (about 100 ms) than in the current study (mean = 226 ms) and less variable. It is tempting to suggest possible reasons for these differences, including the fact that we have been able to exclude adapting monaural influences using the IPD-step stimulus, whereas this and previous studies of binaural adaptation have not been able to do so. The rapid adaptation in auditory nerve–fiber responses (Westerman and Smith 1987Go; Yates et al. 1985Go), and the rapid recovery from adaptation (Yates et al. 1983Go), suggests that the time course of adaptation in the peripheral nerve fibers might dominate the time course of adaptation in higher centers, unless it is somehow filtered out by neurons at subsequent stages. In the current study, all of the neurons were low-frequency neurons sensitive to IPDs. The binaural brain stem neurons from which they receive their input must, therefore have phase-locked inputs by specialized pathways that retain, even enhance, temporal sensitivity. Nevertheless, the responses of low-frequency auditory nerve fibers also adapt, albeit the rapid component of their adaptation is less pronounced (slower) than for high-frequency fibers (Westerman and Smith 1985Go). If the combined effect of hair cell, auditory nerve, and brain stem adaptation is transmitted to brain stem nuclei, MSO neurons will receive a diminishing input during the initial 1 s of our IPD-switch stimulus, and will likely show adaptation at stimulus onset. This raises the question, aside from adaptation at stimulus onset brought about by adapting inputs to the MSO: Are the ongoing responses of MSO neurons themselves subject to adaptation? Recordings from the MSO in response to IPM stimuli indicate little or no influence of rate, direction, or depth of the phase modulation (Spitzer and Semple 1995Go). This suggests that MSO neurons might not be subject to intrinsic adaptation processes, once adaptation of monaural inputs has been accounted for.

Although a proportion of the adapting response to stimulus onset at worst IPD may have been a consequence of the activation of monaural pathways to the IC, it is likely that this response also reflects activity that "bleeds through" the binaural cross-correlator at stimulus onset. This is evident even in the responses to binaural beats in the current study (e.g., top row of Fig. 2, AD), where there is a large response at stimulus onset, corresponding to zero IPD, that is not evoked during the subsequent cycles of the beat.

Consequences for binaural processing

Spitzer and Semple (1998Go) noted that primary binaural neurons in the superior olivary complex—neurons with phase-locked monaural inputs—were considerably less influenced by the differences in the rate, direction or depth of the IPM stimulus than were neurons in the IC, leading them to conclude that a hierarchy of binaural responses existed, with sensitivity to apparent motion cues emerging between the site of primary binaural integration in the brain stem and primary auditory cortex (Malone and Semple 2001Go; Malone et al. 2002Go). Subsequent investigation positing a mechanism of adaptation-of-excitation as being responsible for the differential sensitivity of IC neurons to acoustic motion cues (Ingham et al. 2001Go; McAlpine and Palmer 2002Go; McAlpine et al. 2000Go). This suggests the existence of a hierarchy of adaptation-of-excitation, with presumably wider significance for auditory processing than simply low-frequency acoustic-motion cues. Several recent studies (Borisyuk et al. 2002Go; Cai et al. 1998aGo,bGo) used arbitrarily determined adaptation time constants in an attempt to explain the differential sensitivity of IC neurons to auditorymotion cues. The measure of binaural adaptation we have obtained could be used in models of IC processing to determine the contribution of adaptation mechanisms in the IC to any time-varying stimulus, whether binaural or monaural.

Adaptation and sensitivity to stimulus context

The original observation by Spitzer and Semple (1991Go) of noncontiguous responses to the dynamic IPD cues of interaural phase modulation suggested that instantaneous IPD was not the only factor that influenced the responses of neurons, but that stimulus history also had an influence on discharge rates. Although McAlpine et al. (2000Go) argued that response history rather than stimulus history was the critical determinant in rendering IC neurons sensitive to different directions or depths of the motion cues in IPM, stimulus history and response history are intimately linked, given that the only means by which a stimulus feature, such as its context, can be encoded in the brain is by the pattern of action potentials generated. In this sense, stimulus history is, necessarily, response history.

What possible role might spike-frequency adaptation perform in auditory processing? A number of studies, most notably in the visual system, have suggested a role for adaptation-of-excitation in scaling neural output to take account of, for example, stimulus variance (Brenner et al. 2000Go; Fairhall et al. 2001Go). It is also well described that complex cells of the visual cortex adapt to the local contrast (Carandini and Ferster 1997Go; Laughlin 1989Go; Ohzawa et al. 1982Go), the effect being to position a neuron's dynamic range of discharge rates over the relevant range of contrasts. Cells in the lateral geniculate and simple cells in the visual cortex are sensitive only to absolute contrast (Ohzawa et al. 1982Go). If mechanisms underlying binaural adaptation recorded in the IC are inherent to the IC and not a reflection of potential adaptation process occurring in the MSO, it suggests the emergence of context sensitivity at the level of the auditory midbrain driven by adaptation. Although the role of this context sensitivity remains to be determined, it appears to influence sensitivity to dynamic time and level cues, as well as frequency modulations, suggesting it has some general function in providing sensitivity to stimulus context, and may render neurons more useful in allowing them to adjust their output to take account of local stimulus conditions.

Is adaptation in the IC additive or multiplicative?

One outstanding question concerns the form of spike-frequency adaptation in the IC, whether it is multiplicative or additive (i.e., whether neurons adapt their response by a constant factor or by a constant number of action potentials). Additive and multiplicative processes are implicated throughout the brain in neural gain control, the scaling of discharge rate to match the properties of the excitatory input to the dynamic range of neuronal firing (e.g., Mitchell and Silver 2003Go). Further, auditory tasks that are accomplished by comparison of activity between populations of neurons, as has been recently suggested for sound localization (McAlpine et al. 2001Go), would require neurons to adapt at the same rate independent of their current adaptation state, a feature of multiplicative adaptation. Although the IPD-step stimulus was not designed to test whether adaptation in the IC is additive or multiplicative, we have preliminary evidence that IC neurons adapt multiplicatively (Ingham et al. 2002Go). In addition, adaptation time constants appear to be independent of the current state of adaptation (see Fig. 7), consistent with a multiplicative process. Finally, the population adaptation profile is better described by a proportional reduction in spike rate than by a subtractive reduction in spike rate. Figure 11A shows the distribution of unadapted discharge rates at onset to the 1st step to best IPD. The adaptation profiles in Fig. 11B lead to the distribution of final, steady-state discharge rates in Fig. 11C. The mean onset and steady-state responses were used to calculate the expected population response if spike-frequency adaptation adapts on average by a constant proportion (Fig. 11D) or by a constant number of spikes (Fig. 11E). The observed steady-state distribution (Fig. 11C) better reflects adaptation by a constant proportion, more consistent with a multiplicative process, than adaptation by a constant number of spikes, more consistent with an additive process. Although this is not conclusive evidence for multiplicative adaptation, it argues that the population of IC neurons always adapts proportionally.



View larger version (19K):
[in this window]
[in a new window]
 
FIG. 11. A: distribution of peak discharge rates derived from the exponential fits to responses recorded to the 1st step to best IPD. B: exponential decay curves illustrating spike frequency adaptation to the initial step to best IPD (replotted from Fig. 3A). Filled circles indicate the mean values (±SD) of peak discharge rate (left) and adapted steady-state discharge rate (right). C: frequency distribution histogram of adapted discharge rate from exponential fits of responses recorded to the 1st step to best IPD. D: distribution of steady-state firing rate predicted by multiplicative gain scaling of peak discharge rates. E: distribution of steady-state firing rate predicted by additive gain offset of peak discharge rates.

 
Cellular mechanisms of adaptation in the IC

The term "adaptation" appears to be used to describe a range of phenomenological observations, and does not provide much insight into any specific mechanisms that might account for neuronal response properties. Potential mechanisms include intrinsic membrane properties involving small-conductance calcium-activated potassium, or SK, channels (Cordoba-Rodriguez et al. 1999Go; Tang et al. 1997Go; Wagner et al. 2001Go; Wikström and El Manira 1998Go). SK channels have been localized within the mammalian IC using in situ hybridization (Stocker and Pedarzani 2000Go). Prolonged depolarization of the cell membrane—such as that seen during high discharge rates—results in the slow activation of an inward potassium current (Xia et al. 1998Go), which builds up over time producing a significant hyperpolarizing effect on a neuron's resting membrane potential. Such a mechanism is supported by in vitro studies indicating that <=25% of intracellularly recorded responses of IC neurons show adapting spike patterns in response to depolarizing current injection (Li et al. 1998Go; Peruzzi et al. 2000Go; Sivaramakrishnan and Oliver 2001Go). An alternative mechanism to account for the adaptation profiles we observe involves neuronal networks, either local or from other midbrain regions, feeding inhibitory input onto IC neurons (Burger and Pollak 2001Go; Chen et al. 1999Go; Kelly and Li 1997Go). However, recordings from a single IC neuron indicate that adaptation is still evident even when GABA-ergic inhibition is blocked using bicuculline methiodide, which did not abolish the neuron's adapting response profile (Fig. 12). This suggests that GABA-ergic inhibition does not underlie the adaptation process observed in the IC, consistent with recent data indicating that blocking GABA-ergic inhibition does not abolish sensitivity to the dynamic interaural phase cues of IPM (McAlpine and Palmer 2002Go). Such sensitivity has been suggested to arise from a process of spike-frequency adaptation (McAlpine et al. 2000Go). Interpretation of this particular result should be treated cautiously, however, because bicuculline methiodide is also known to alter directly SK channel currents (e.g., Tonkovic-Capin et al. 2001Go). A third possibility is the involvement of fast synaptic plasticity, such as that observed in the barrel cortex (Chung et al. 2002Go). However, this fast synaptic plasticity appears to have a slower time course than the adapting response profiles we observe, although this does not necessarily rule it out as a potential mechanism.



View larger version (19K):
[in this window]
[in a new window]
 
FIG. 12. Response of a single IC neuron to an IPD-step stimulus (duration 1,000 ms) before (white), during (black), and after recovery from (gray) iontophoresis of bicuculline methiodide. Adaptation time constants were 74, 93, and 56 ms, respectively.

 
The question remains, how much of the adaptation we observe is intrinsic to IC neurons, and how much is imposed from the level of primary binaural integration in the brain stem nuclei of the medial superior olive? Unlike IC neurons, MSO neurons are largely insensitive to dynamic interaural phase cues (Spitzer and Semple 1995Go), suggesting their responses do not adapt to any great extent. MSO neurons possess certain cellular specializations for precise coincidence detection such as AMPA receptors (reviewed in Trussell 1997Go), ideal for the rapid and secure transmission of high-rate, phase-locked inputs critical to their primary role in binaural coincidence detection. Once binaural integration is complete, however, such precision may no longer be required. Thus any change in sensitivity to dynamic interaural cues may reflect a transition from a purely "cues"-based encoding of interaural disparities in the MSO to a more "context"-based encoding of interaural disparities in the IC. One could conclude that spike-frequency adaptation in the IC is brought about by intrinsic membrane properties of IC neurons themselves, perhaps by specific ion channels, rather that by the action of networks of inhibitory connections or binaurally adapting inputs from the MSO.


    ACKNOWLEDGMENTS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 ACKNOWLEDGMENTS
 REFERENCES
 
We thank Prof. A. R. Palmer and Dr. T. Shackleton of the Medical Research Council (MRC) Institute for Hearing Research, Nottingham, UK for the use of software to drive these experiments.

Present address: N. J. Ingham, Centre for the Neural Basis of Hearing, The Physiological Laboratory, Downing St., Cambridge, CB2 3EG, UK.

GRANTS

This work was supported by an MRC Career Establishment Grant awarded to D. McAlpine.


    FOOTNOTES
 
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

Address for reprint requests and other correspondence: D. McAlpine, Department of Physiology, University College London, Gower Street, London, WC1E 6BT, UK (E-mail: d.mcalpine{at}ucl.ac.uk).


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 ACKNOWLEDGMENTS
 REFERENCES
 
Boehnke SE and McAlpine D. Effects of stimulus level on ITD tuning of low frequency cells in the auditory midbrain. Abstr 24th Midwinter Res Meet Assoc Res Otolaryngol 190, 2001.

Borisyuk A, Semple MN, and Rinzel J. Adaptation and inhibition underlie responses to time-varying interaural phase cues in a model of inferior colliculus neurons. J Neurophysiol 88: 2134–2146, 2002.[Abstract/Free Full Text]

Brenner N, Bialek W, and de Ruyter VS. Adaptive rescaling maximizes information transmission. Neuron 26: 695–702, 2000.[CrossRef][ISI][Medline]

Burger RM and Pollak GD. Reversible inactivation of the dorsal nucleus of the lateral lemniscus reveals its role in the processing of multiple sound sources in the inferior colliculus of bats. J Neurosci 21: 4830–4843, 2001.[Abstract/Free Full Text]

Cai H, Carney LH, and Colburn HS. A model for binaural response p