JN Fuel your research with LabChart
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
QUICK SEARCH:   [advanced]


     

J Neurophysiol 94: 1655-1658, 2005. First published April 20, 2005; doi:10.1152/jn.01226.2004
0022-3077/05 $8.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
94/2/1655    most recent
01226.2004v1
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 ISI Web of Science (3)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Wagner, H.
Right arrow Articles by Carr, C. E.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Wagner, H.
Right arrow Articles by Carr, C. E.

REPORT

Microsecond Precision of Phase Delay in the Auditory System of the Barn Owl

Hermann Wagner1,2, Sandra Brill1, Richard Kempter3 and Catherine E. Carr2

1Institute for Biology II, Rheinisch-Westfölische Technische Hochschule Aachen, Aachen, Germany; 2Department of Biology, University of Maryland, College Park, Maryland; and 3Institute for Theoretical Biology, Humboldt University Berlin, and Neuroscience Research Centre, Charité Medical Faculty of Berlin, and Bernstein Center for Computational Neuroscience, Berlin, Germany

Submitted 1 December 2004; accepted in final form 7 April 2005


   ABSTRACT
TOP
 ABSTRACT
INTRODUCTION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
The auditory system encodes time with sub-millisecond accuracy. To shed new light on the basic mechanism underlying this precise temporal neuronal coding, we analyzed the neurophonic potential, a characteristic multiunit response, in the barn owl's nucleus laminaris. We report here that the relative time measure of phase delay is robust against changes in sound level, with a precision sharper than 20 µs. Absolute measures of delay, such as group delay or signal-front delay, had much greater temporal jitter, for example due to their strong dependence on sound level. Our findings support the hypothesis that phase delay underlies the sub-millisecond precision of the representation of interaural time difference needed for sound localization.


   INTRODUCTION
TOP
 ABSTRACT
 INTRODUCTION
GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
The barn owl is well known for its superb sound-localization capabilities (Bala et al. 2003Go; Carr and Konishi 1990Go; Gerstner et al. 1996Go; Keller et al. 1998Go; Kempter et al. 2001Go; Koppl 1997Go; Moiseff and Konishi 1981Go; Pena et al. 1996Go; Sullivan and Konishi 1984Go, 1986Go; Viete et al. 1997Go). The owl uses interaural time difference (ITD) to encode sound azimuth with a behavioral accuracy of <10 µs (Bala et al. 2003Go) and a neuronal sensitivity of 25–100 µs (Bala et al. 2003Go; Moiseff and Konishi 1981Go). The highest overall monaural temporal sensitivity has been measured in the third-order nucleus laminaris (NL) where binaural convergence creates tuning to ITD (Carr and Konishi 1990Go; Reyes et al. 1996Go; Schwarz 1992Go; Sullivan and Konishi 1986Go). The NL exhibits a characteristic frequency-following multiunit response termed the neurophonic (Schwarz 1992Go; Snyder and Schreiner 1984Go; Sullivan and Konishi 1986Go), which we used to study timing. The neurophonic well represents both monaural and binaural temporal sensitivity (Sullivan and Konishi 1986Go).

Earlier measurements of the conduction time from the ear to NL have found values between 2 and 3 ms (Carr and Konishi 1990Go). Sullivan and Konishi (1984)Go already mentioned the importance of phase but did not quantify temporal precision. Koppl (1997)Go quantified temporal precision in the second-order nucleus magnocellularis that provides input to the nucleus laminaris but found that the response delay depended on sound level. The rate of change amounted to ~5 µs/dB around the characteristic frequency. Taking into account that the range of interaural level differences experienced by barn owls is ~20 dB (Keller et al. 1998Go; Viete et al. 1997Go), how is it possible that the ITD can be represented in NL of owls with a neuronal precision much sharper than 100 µs?

The auditory system encodes delay both as conduction time, an absolute time measure (Fitzgerald et al. 2001Go; Goldstein et al. 1971Go; Ruggero 1980Go) and, by phase locking, a relative time cue (Anderson et al. 1971Go; Carr and Konishi 1990Go; Koppl 1997Go; Reyes et al. 1996Go; Sullivan and Konishi 1984Go, 1986Go). Absolute and relative time codes are also known in physics, where a distinction is made between group velocity and phase velocity, thus resulting in group and phase delays (Anderson et al. 1971Go; Fitzgerald et al. 2001Go; Goldstein et al. 1971Go; Koppl 1997Go; Ruggero 1980Go). While group delay describes the latency of the envelope of a band-pass-filtered signal, phase delay refers to the times of occurrence of its peaks and troughs. The high-frequency limit of group delay, the signal-front delay, has also been used as a measure of delay (Fitzgerald et al. 2001Go; Ruggero 1980Go). We analyzed data obtained from the NL of the barn owl to determine which measure of delay was suited for precise and level-independent representation of ITD.

Nine barn owls (Tyto alba pratincola) were used in this study. The procedures conformed to National Institutes of Health guidelines for animal research and were approved by the animal care and use committee of the University of Maryland. In contrast to earlier studies, the analogue waveform of the neurophonic potential in or close to NL was recorded at a sampling period of 20.8 µs with commercial, Epoxylite-coated tungsten electrodes (Frederick Haer, Brunswick, ME) with impedances of 2–8 M{Omega}. Neurophonic recordings had the advantage of being stable for ≥1 h and allowed measurements of local multiunit activity. Specific recording sites were defined by combining stereotaxic techniques, physiological characterization, and histologically verified lesions.

Acoustic stimuli (clicks and noises) were digitally generated by custom-written software ("Xdphys" written in Dr. M. Konishi's lab at the California Institute of Technology, Pasadena, CA) driving a signal-processing system (Tucker Davies Technology, Gainesville, FL). Clicks had a rectangular form of varying intensity [0 dB (corresponding to 65 dB SPL) to 40-dB attenuation] and a duration of two samples (equivalent to 41.6 µs). Only condensation clicks were used. The standard click had 0 dB attenuation.

Neurophonic responses to clicks were recorded in the 3.5- to 7-kHz region of the tonotopically organized NL. The spontaneous activity (10 ms before click presentation) as well as the driven activity (10 ms after click presentation) were stored. Clicks were repeated 128 times (Fig. 1 A). The driven activity contained an oscillatory response (Fig. 1, A and B). Its envelope increased smoothly within ~1 ms and fell off almost symmetrically. The oscillation under the envelope typically exhibited a complex waveform containing several spectral components. Fourier analysis showed that one or two components were <2 kHz (Fig. 1C). Another component was close to the best frequency as obtained from iso-level frequency response curves. Because we wanted to study processes related to frequency tuning, only the high-frequency component was analyzed. Therefore the neurophonic potential was high-pass filtered to reveal the oscillation of the high-frequency component alone (Fig. 1D). Auditory filtering is well described by gammatone functions and their derivatives (Irino and Patterson 2001Go; Tan and Carney 2003Go). Thus the high-pass filtered click-evoked response was fitted with a Gammatone function of order 3 (Fig. 1D).



View larger version (65K):
[in this window]
[in a new window]
  		FIG. 1. Analysis of neurophonic potentials evoked with standard clicks. Data from monaural stimulation. A: time course of neurophonic potentials. The x axis shows the time after recording onset. The click was presented at 10 ms after recording onset. The 128 individual traces are stacked on top of each other. Scale bar: 25 mV. B: mean of the 128 traces shown in A. Scale bar: 2.5 mV. C: normalized amplitude of the Fourier transform (512 points) of the mean response from 9.984 ms (data point 480) to the end of the recording. Note the different spectral components. The dotted line indicates the high-pass filter used for the separation of the components. D: high-pass filtered click-evoked mean response (—) and fit of this response with a Gammatone function of order 3 ( · · · ). The envelope of the Gammatone function is shown by the · - · line. The group delay ({tau}gr) corresponds to the peak of the envelope. The signal-front delay ({tau}fr) and the phase delay ({tau}ph) were determined from the filtered, but unfitted data (see also text). Scale bar: 2.5 mV. E: histogram of the variation of the 3 delay measures obtained from 96 recording sites (176 data sets, 89 left and 87 right stimulation, 128 repetitions each). · · · : phase delay, - - -: signal-front delay, —: group delay. The bin width was 10 µs for the phase delay, 20 µs for the group delay, and 30 µs for the signal-front delay. Note that SDs are plotted, allowing for an analysis with a bin width that is smaller than the sampling interval.

 
In the click-evoked response, typically several peaks and troughs could be distinguished. We considered only local extrema occurring after stimulus presentation and having an amplitude greater than the mean plus 2 SDs of the background noise for at least three consecutive data points. The latency of the first extremum detected in this way was the signal-front delay (Fig. 1D). The group delay was assigned to the maximum of the envelope of the Gammatone function (Fig. 1D). The extremum closest to the group delay, determined in the first of the 128 trials that was within 1 SD of the average group delay, was chosen for the phase delay (Fig. 1D). To estimate the variability or "jitter" of each type of delay, we used the respective extrema in each of the 128 traces (Fig. 1A). The SD of these delays determined over the 128 repetitions was taken as a measure for the variability and was plotted as a data point in a histogram (Fig. 1E). In our sample of 176 data sets, the signal-front delay had the largest jitter (median: 500 µs), whereas the jitter in the phase delay was smallest (median: 10.4 µs). The jitter in the group delay (median: 64 µs) was between these two extremes. Thus the temporal precision necessary for ITD coding (Bala et al. 2003Go; Moiseff and Konishi 1981Go) could be achieved with the phase delay and the group delay but not with the signal-front delay. Therefore we did not consider the signal-front delay further.

In a second experiment, we decreased click amplitudes ≤40 dB. As is typically observed in audition, the group delay increased as the stimulus level decreased (Fig. 2). In the example shown in Fig. 2A, however, peak number 4 at 0 dB attenuation coincided with peak number 3 at 20-dB attenuation and with the barely visible peak number 1 at 40-dB attenuation (vertical line in Fig. 2A). Note that the definition of phase delay allows for a jumping between subsequent peaks (Fitzgerald et al. 2001Go). In 72 data sets obtained from 43 recording sites, phase delay remained essentially constant with a narrow distribution around one sample point (mean: –3 µs, Fig. 2B). In contrast, group delay increased ≤0.6 ms when click level was reduced by 20 dB (Fig. 2B), in agreement with the result in Carr and Konishi (1990)Go.



View larger version (27K):
[in this window]
[in a new window]
  		FIG. 2. Level dependence of click-evoked responses. Data from monaural stimulation. A: when the stimulus level decreased (numbers on left side refer to attenuation in dB), the shape of the waveform remained oscillatory, the amplitude decreased and the response occurred at longer latencies. However, as indicated by the vertical line, the position of the response maxima remained stable. Scale bar: 10 mV. B: level dependence of delays in the pooled data of 72 data sets. Stimulus amplitudes were decreased from 0- to 20-dB attenuation. To determine the changes in phase or group delay induced by changes in click level, the times of arrival of the respective response extrema evoked by the different click intensities were picked from the mean curves and subtracted from each other. For phase delay values, this sometimes included jumps between peaks as explained in the text. Each difference provided a data point for the histograms. —, group delay; · · · , phase delay. Bin widths: 20 µs for the phase delay and 40 µs for the group delay.

 
The low variability and the level invariance of phase delay implied that this delay measure would be the most reliable code for representing the behaviorally relevant interaural time difference (ITD). The "best ITD" at a given recording site may be computed from two phase delays obtained through monaural stimulation (Sullivan and Konishi 1986Go): we subtracted the phase delay for stimulation of the left ear from the phase delay for stimulation of the right ear. This subtraction lead to a "best ITD" that was independent of interaural level difference (ILD) (compare Fig. 3, A for 0 dB ILD with B for 10 dB ILD). On the other hand, ITDs between two group delays depended on the ILD. Even though an ITD computed from the group delay may be zero for 0 dB ILD (Fig. 3A), it changed significantly for 10 dB ILD (Fig. 3B).



View larger version (14K):
[in this window]
[in a new window]
  		FIG. 3. Influence of interaural level difference (ILD) on interaural time difference (ITD). Data from monaural stimulation (A and B) and binaural stimulation (C). A: at equal levels of the clicks presented to the left and right ears (5 dB attenuation: L5, R5), the difference between right and left group delays ({Delta}{tau}gr) and phase delays ({Delta}{tau}ph) were near 0 µs. B: when the level on the left side was decreased by 10 dB (L15), group delay increased, while the phase delay stayed constant. C: in the binaural situation with broadband noise as the stimulus, the positions of the extrema of the ITD curve for an ILD of 0 dB (—) were almost identical to the extrema of the ITD curve measured with an ILD of 10 dB ( · · · ). Events were obtained from an arbitrary setting of a TTL level trigger. Because {Delta}{tau}gr had changed but the ITD curve did not change, the phase delay, but not the group delay is important for representing ITD. Scale bars in A and B: 2.5 mV.

 
Because the delay differences obtained from monaural responses in Fig. 3, A and B, only indirectly represented the binaural situation, we also tested the level dependence of ITD tuning with binaural stimuli in a few (n = 6) cases. In Fig. 3C, we plotted the response at a given recording site as a function of the ITD of binaural noise in contrast to our previous results obtained through neurophonic potentials in response to monaurally presented clicks. The ITD tuning curves for ILDs of 0 and 10 dB were virtually identical in agreement with previous findings that ITD tuning for binaural stimulation is stable under varying conditions of stimulus level (Pena et al. 1996Go; Viete et al. 1997Go). Thus the level tolerance of ITD tuning as shown in Fig. 3C independently demonstrates the importance of phase delay.

Thus phase delay, and not group delay or signal-front delay, appears to underlie ITD tuning. Single-unit recordings from auditory nerve and cochlear nucleus support this conclusion (Koppl 1997Go; Sullivan and Konishi 1984Go). Likewise, auditory nerve fibers of squirrel monkeys show very low phase-delay jitter when stimulated with sinusoids near the characteristic frequency (Anderson et al. 1971Go).

The level independence of phase delay is consistent with a variety of filters proposed for peripheral auditory processing (Irino and Patterson 2001Go; Tan and Carney 2003Go). The remarkable stability of phase delay is consistent with the model of Gerstner et al. (1996)Go and Kempter et al. (2001)Go, who predicted that during development synapses from NM to NL and axonal arbors from NM to NL are selected in such a way that phase delays are similar. These authors also argued that only such a selection allows for using phase delay to code temporal information in the NL and to represent ITD. Note that this conclusion is in line with the existence of a neurophonic potential in adult animals when it is assumed that the neurophonic potential is typically the summed response of an ensemble of magnocellular axons. A coherent summation of responses of different axons is only feasible when we have a coincident arrival of volleys of phase-locked spikes at the borders of NL and a coherent transmission of spikes through the nucleus. In other words, theory predicts that phase delays in different magnocellular axons must be similar.

Timing is important in many neuronal systems. It plays a role in models of learning (spike-timing-dependent plasticity) in feature binding as well as in precise reactions to dynamic stimuli—such as approaching targets. To compare the precision of the different systems, a temporal quality factor is helpful. The coefficient of variation (CV), defined as the quotient of SD and mean, may be a good criterion. In our example, the CV for the phase delay is ~0.01. In systems such as the visual cortex, temporal jitter is much larger (Bair et al. 2002Go; Bisley et al. 2004Go) (CV ~ 0.1), whereas values similar to the CV observed in the owl's NL are found in the electrosensory system (Carr et al. 1986Go) and the auditory system of bats (Covey and Casseday 1991Go) and may be computed from synfire chains (Abeles et al. 1993Go).


   GRANTS
TOP
 ABSTRACT
 INTRODUCTION
 GRANTS
ACKNOWLEDGMENTS
 REFERENCES
 
This research was sponsored by the German Research Foundation (DFG, Wa-606/12, Ke-788/1–3) and by National Institutes of Health Grants DC-000636 to C. E. Carr and by P30 04664.


   ACKNOWLEDGMENTS
TOP
 ABSTRACT
 INTRODUCTION
 GRANTS
 ACKNOWLEDGMENTS
REFERENCES
 
M. Knepper and E. Smith helped with measuring and analyzing click responses. R. Schätte made helpful comments on the manuscript.


   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: H. Wagner, Institut für Biologie II, RWTH Aachen, Kopeinikusstra{beta}e 16, D-52074 Aachen


   REFERENCES
TOP
 ABSTRACT
 INTRODUCTION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Abeles M, Bergman H, Margalit E, and Vaadia E. Spatiotemporal firing pattern in the frontal cortex of behaving monkeys. J Neurophysiol 70: 1629–1638, 1993.[Abstract/Free Full Text]

Anderson DJ, Rose JE, Hind JE, and Brugge JF. Temporal position of discharges in single auditory nerve fibres within the cycle of a sine-wave stimulus: frequency and intensity effects. J Acoust Soc Am 49: 1131–1139, 1971.

Bair W, Cavanaugh JR, Smith MA, and Movshon JA. The timing of response onset and offset in macaque visual neurons. J Neurosci 22: 3189–3205, 2002.[Abstract/Free Full Text]

Bala AD, Spitzer MW, and Takahashi TT. Prediction of auditory spatial acuity from neural images on the owl's auditory space map. Nature 424: 771–774, 2003.[CrossRef][Medline]

Bisley JW, Krishna BS, and Goldberg ME. A rapid and precise on-response in posterior parietal cortex. J Neurosci 24: 1833–1838, 2004.[Abstract/Free Full Text]

Carr CE, Heiligenberg W, and Rose GTJ. A time-comparison circuit in the electric fish midbrain. I. Behavior and physiology. J Neurosci 6: 107–119, 1986.[Abstract]

Carr CE and Konishi M. A circuit for detection of interaural time differences in the brain stem of the barn owl. J Neurosci 10: 3227–3246, 1990.[Abstract]

Covey E and Casseday JH. The monaural nuclei of the lateral lemniscus in an echolocating bat: parallel pathways for analyzing temporal features of sound. J Neurosci 11: 3456–3470, 1991.[Abstract]

Fitzgerald JV, Burkitt AN, Clark GM, and Paolini AG. Delay analysis in the auditory brain stem of the rat: comparison with click latency. Hearing Res 159: 85–100, 2001.[CrossRef][Medline]

Gerstner W, Kempter R, van Hemmen JL, and Wagner H. A neuronal learning rule for submillisecond temporal coding. Nature 383: 76–78, 1996.[CrossRef][Medline]

Goldstein JL, Bear T, and Kiang NYS. A theoretical treatment of latency, group delay, and tuning characteristics for auditory nerve responses to clicks and tones. In: Physiology of the Auditory System, edited by Sachs MB. Baltimore, MD: National Education Consultants, 1971, p. 133–141.

Irino T and Patterson RD. A compressive gammachirp auditory filter for both physiological and psychophysical data. J Acoust Soc Am 109: 2008–2022, 2001.[Medline]

Keller CH, Hartung K, and Takahashi TT. Head-related transfer functions of the barn owl: measurement and neural responses. Hear Res 118: 13–34, 1998.[CrossRef][ISI][Medline]

Kempter R, Leibold C, Wagner H, and van Hemmen JL. Formation of temporal-feature maps by axonal propagation of synaptic learning. Proc Natl Acad Sci USA 98: 4166–4171, 2001.[Abstract/Free Full Text]

Koppl C. Phase locking to high frequencies in the auditory nerve and cochlear nucleus magnocellularis of the barn owl, Tyto alba. J Neurosci 17: 3312–3321, 1997.

Moiseff A and Konishi M. Neuronal and behavioral sensitivity to binaural time differences in the owl. J Neurosci 1: 41–49, 1981.

Pena JL, Viete S, Albeck Y, and Konishi M. Tolerance to sound intensity of binaural coincidence detection in the nucleus laminaris of the owl. J Neurosci 16: 7046–7054, 1996.[Abstract/Free Full Text]

Reyes AD, Rubel EW, and Spain WJ. In vitro analysis of optimal stimuli for phase-locking and time-delayed modulations of firing in avian nucleus laminaris neurons. J Neurosci 16: 993–1007, 1996.[Abstract/Free Full Text]

Ruggero MA. Systematic errors in indirect estimates of basilar membrane travel times. J Acoust Soc Am 67: 707–710, 1980.[CrossRef][ISI][Medline]

Schwarz DWF. Can central neurons reproduce sound waveforms? An analysis of the neurophonic potential in the laminar nucleus of the chicken. J Otolaryng 21: 30–38, 1992.

Snyder RL and Schreiner CE. The auditory neurophonics: basic properties. Hear Res 15: 261–280, 1984.[Medline]

Sullivan WE and Konishi M. Segregation of stimulus phase and intensity coding in the cochlear nucleus of the barn owl. J Neurosci 4: 1787–1799, 1984.[Abstract]

Sullivan WE and Konishi M. Neural map of interaural phase difference in the owl's brain stem. Proc Natl Acad Sci USA 83: 8400–8404, 1986.[Abstract/Free Full Text]

Tan Q and Carney LH. A phenomenological model for the responses of auditory-nerve fibers. II. Nonlinear tuning with a frequency glide. J Acoust Soc Am 114: 2007–2020, 2003.[Medline]

Viete S, Pena JL and Konishi M. Effects of interaural intensity difference on the processing of interaural time difference in the owl's nucleus laminaris. J Neurosci 17: 1815–1824, 1997.[Abstract/Free Full Text]





This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
94/2/1655    most recent
01226.2004v1
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 ISI Web of Science (3)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Wagner, H.
Right arrow Articles by Carr, C. E.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Wagner, H.
Right arrow Articles by Carr, C. E.

HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online
Copyright © 2005 by the The American Physiological Society.