Original Research ARTICLE
Biophysical basis of the sound analog membrane potential that underlies coincidence detection in the barn owl
- 1Department of Biology, University of Maryland, College Park, MD, USA
- 2Systems Biology, Osaka Bioscience Institute, Suita, Japan
- 3Division of Biology, California Institute of Technology, Pasadena, CA, USA
Interaural time difference (ITD), or the difference in timing of a sound wave arriving at the two ears, is a fundamental cue for sound localization. A wide variety of animals have specialized neural circuits dedicated to the computation of ITDs. In the avian auditory brainstem, ITDs are encoded as the spike rates in the coincidence detector neurons of the nucleus laminaris (NL). NL neurons compare the binaural phase-locked inputs from the axons of ipsi- and contralateral nucleus magnocellularis (NM) neurons. Intracellular recordings from the barn owl's NL in vivo showed that tonal stimuli induce oscillations in the membrane potential. Since this oscillatory potential resembled the stimulus sound waveform, it was named the sound analog potential (Funabiki et al., 2011). Previous modeling studies suggested that a convergence of phase-locked spikes from NM leads to an oscillatory membrane potential in NL, but how presynaptic, synaptic, and postsynaptic factors affect the formation of the sound analog potential remains to be investigated. In the accompanying paper, we derive analytical relations between these parameters and the signal and noise components of the oscillation. In this paper, we focus on the effects of the number of presynaptic NM fibers, the mean firing rate of these fibers, their average degree of phase-locking, and the synaptic time scale. Theoretical analyses and numerical simulations show that, provided the total synaptic input is kept constant, changes in the number and spike rate of NM fibers alter the ITD-independent noise whereas the degree of phase-locking is linearly converted to the ITD-dependent signal component of the sound analog potential. The synaptic time constant affects the signal more prominently than the noise, making faster synaptic input more suitable for effective ITD computation.
The ability to tell the direction of the sound source, or sound localization, is a fundamental auditory function in many animal species. Among various species examined (see Klump, 2000; Heffner and Heffner, 2003; for reviews), the barn owl, which can locate its prey in the total darkness purely on the basis of acoustic cues (Payne, 1971; Konishi, 1973), shows great sound localization acuity, with a minimum discriminable angle of a few degrees (Knudsen et al., 1979; Bala et al., 2003). The auditory system of the barn owl computes the interaural time difference (ITD) to determine the azimuthal location of the sound source (Konishi, 1993). In birds, the ITD, or the time difference of sound arrival between two ears, is computed in a specialized neural circuit: axons from the cochlear nucleus magnocellularis (NM) form delay lines, and neurons of the nucleus laminaris (NL) detects coincident inputs from ipsi- and contralateral NM axons (Jeffress, 1948; Carr and Konishi, 1990). Physiological studies showed that barn owls' NL neurons vary their discharge rates with changes in ITDs of less than 10 μs (Carr and Konishi, 1990; Peña et al., 1996). A wide variety of highly-specialized cellular, synaptic, and network mechanisms underlie this temporal acuity (see Grothe et al., 2010; Ashida and Carr, 2011, for recent reviews).
In our previous reports, we investigated the cellular properties of owl NL neurons using in vivo intracellular recordings (Funabiki et al., 2011) and modeling (Ashida et al., 2007; Funabiki et al., 2011). When tonal stimuli were presented, the membrane potential of the NL neuron oscillated at the same frequency as the stimulus tone and was thus named the “sound analog potential” (Funabiki et al., 2011). The amplitude of this sound-induced oscillation changes periodically with ITD, showing an almost linear relationship with the spike rate of the NL neuron. Previous modeling results demonstrated that convergence of phase-locked excitatory synaptic inputs from NM fibers gives rise to oscillatory membrane potentials (Kempter et al., 1998; Slee et al., 2010; see Figure 1A for an example) and that the sound analog potential observed in NL can be quantitatively reproduced if biologically plausible parameters are chosen (Ashida et al., 2007; Funabiki et al., 2011). How presynaptic, synaptic, and postsynaptic properties, however, affect the formation of the sound analog potential remains to be elucidated.
Figure 1. Example traces of the model membrane potential with default parameters. The default parameters used for this example are summarized in Table 1. (A) Simulated membrane potential oscillating at 4 kHz. (B) AC component of the simulated membrane potential. Amplitude = 1.25 mV (see Materials and Methods for definition). Average peak-to-peak height = 2.50 mV. (C) Noise component of the simulated membrane potential. Amplitude (measured by the time-averaged standard deviation) is 1.03 mV in this example.
In the accompanying paper, we focused on the theoretical relationship between phase-locked inputs and membrane potential oscillations (Ashida et al., 2013). The key factors which can affect the formation of the sound analog potentials include (1) the frequency of the tonal stimulus, (2) the mean firing rate of NM neurons, (3) the number of converging NM fibers per NL neuron, (4) the average degree of phase-locking of these fibers, (5) the time scale of unitary synaptic input, and (6) the membrane properties of the NL neuron. Among these factors, the membrane properties have been extensively studied. A number of studies have shown that ITD coding in the avian NL and in the mammalian medial superior olive (MSO) is dynamically controlled by the low threshold potassium (KLVA) channels (Manis and Marx, 1991; Reyes et al., 1994; Svirskis et al., 2002; Rothman and Manis, 2003; Day et al., 2008; Gai et al., 2009; Jercog et al., 2010; Mathews et al., 2010), hyperpolarization-activated cation channels (Yamada et al., 2005; Khurana et al., 2011), and by optimized fast sodium channels (Kuba et al., 2006; Ashida et al., 2007; Scott et al., 2010). The accompanying paper (Ashida et al., 2013) examines the first key factor, the effect of tonal frequency on the formation of the sound analog potential. In the present paper, we focus on the remaining four factors (2–5 above) and systematically examine how the sound analog potential and ITD coding of the NL neuron depends on these presynaptic and synaptic factors.
Phase-locked inputs and oscillatory potentials
Our model consists of NM fibers and a single-compartment NL cell body. The phase-locked spiking activity of each NM fiber is modeled as an inhomogeneous Poisson process with a periodic intensity function oscillating at the stimulus sound frequency. Spikes of converging NM fibers are summed and then modified by the synaptic and membrane filters, producing an oscillation in the membrane potential of the NL neuron, namely, the sound analog potential (Figure 1A). The mean synaptic input is referred to as the “DC” component, while the main oscillation component at the stimulus frequency is called the “AC” or “signal” component (Figure 1B) because the NL neuron changes its spike rate almost linearly with the amplitude of the AC component (Funabiki et al., 2011). All the other frequency components, including higher harmonics, are regarded as “noise” (Figure 1C) because they do not encode ITDs (Ashida et al., 2007; Slee et al., 2010). Analytical expressions that relate model parameters to the DC, AC, and noise components are summarized in Table 1. The accompanying paper (Ashida et al., 2013) shows that with default parameters (Table 1) the linear approximation gives good predictions for the DC, AC, and noise of simulated synaptic inputs and oscillatory membrane potentials. In the following sections, we examine how sound analog potentials are controlled by the number of presynaptic NM fibers, the mean firing rate of these fibers, their average degree of phase-locking measured by vector strength (VS), and the synaptic time scale measured by the half peak width of the unitary synaptic input. Note that, in our simulations, we assumed that ipsi- and contralateral NM inputs arrived perfectly in-phase, resulting in the maximum oscillation amplitude. The ITD-dependence of the sound analog potential will be examined in the section titled Implications for ITD coding.
Mean Spike Rate of NM
The effect of the average NM spike rate λ0 is fairly simple. As the linear approximation theory indicates (Equations 1–5 in Table 1), traces of simulated membrane potentials become less noisy as the mean spike rate of presynaptic NM fibers increases (Figure 2A). This observation is confirmed by the power spectral density curves (Figure 2B). Increases in the NM spike rate reduce noise but retain AC and higher harmonics (Figures 2C,D), provided that the average input level (DC) is kept constant. Note that the power spectrum density is discontinuous at frequencies f = nfs (fs: signal frequency; n = 0, ±1, ±2,…), corresponding to the DC, AC, and higher harmonics (see Ashida et al., 2013, for further analytical formulations). The noise amplitudes of the synaptic input and the membrane potential decrease linearly with the square root of the spike rate of NM fibers (Equations 3 and 5). The typical spike rate of an NM neuron is over 400 spikes/s (Peña et al., 1996), which is much greater than that of an auditory nerve fiber (Köppl and Yates, 1999). The high spiking rate of NM neurons thus contributes to noise reduction in NL neurons.
Figure 2. Dependence of synaptic input in NL on the firing rate of presynaptic NM fibers. (A) Simulated traces of the model membrane potential. The number above each trace shows the output spike rate of NM fibers. The traces become less noisy as the NM rate increases. (B) Power spectral densities of the four traces shown in (A). Low frequency noise components decrease with increasing NM rates, while peaks at the input frequency and higher harmonics remain unchanged. (C) Dependence of the DC, AC, and noise amplitudes of the simulated synaptic input on the mean spike rate of NM fibers. (D) Dependence of the AC and noise amplitudes of the simulated membrane potential on the mean spike rate of NM fibers. Solid lines in (C) and (D) are obtained from analytical calculations (Equations 1–5). Vertical broken gray lines in (C) and (D) show the default parameter (500 Hz) used in our simulations.
Number of Converging NM Fibers
The dependence of the AC and noise components on the number M of presynaptic fibers is similar to the dependence on the mean spike rate λ0 (Equations 2–5 in Table 1). Traces of simulated membrane potential become less noisy as the number M of presynaptic NM fibers increases (Figure 3A). Power spectral density curves support this observation (Figure 3B). The noise amplitude decreases linearly with the square root of M (Figures 3C,D), as predicted by the theoretical calculations (Equations 3 and 5). The convergence of large numbers of NM fibers results in stable sinusoidal inputs (Figure 3A, bottom). By contrast, when the number of presynaptic fibers is small, the overall potential waveform is distinct from a pure sinusoid and each unitary synaptic event becomes discernible (Figure 3A, top) even though the calculated AC component itself is the same (Figures 3C,D). An NL cell receives a few hundred NM afferents (Carr and Boudreau, 1993), whereas an NM neuron receives 1–4 auditory nerve inputs via large endbulb synapses (Carr and Boudreau, 1991). Thus, the membrane potential waveforms should be different between NM and NL, potentially reflecting the difference in their computational roles in ITD coding.
Figure 3. Dependence of the synaptic input in NL on the number of presynaptic NM fibers. (A) Simulated traces of the model membrane potential. The number above each trace shows the numbers of NM fibers. The traces become less noisy as the number of NM fibers increases. (B) Power spectral densities of the four traces shown in (A). Low frequency noise components decrease with increasing numbers of NM fibers, while peaks at the input frequency and higher harmonics remain unchanged. Although the trace with 3 NM fibers (A, top) looks considerably different from the other three traces, its AC component (B, top) has the same amplitude as the other three. (C) Dependence of the DC, AC, and noise amplitudes of the simulated synaptic input on the number of NM fibers. (D) Dependence of the AC and noise amplitudes of the simulated membrane potential on the number of NM fibers. Solid lines in (C) and (D) are obtained from analytical calculations (Equations 1–5). Vertical broken gray lines in (C) and (D) show the default parameter (300 fibers) used in our simulations.
Degree of Phase-Locking
where N is the total number of spikes in the sequence, and tj is the timing of the j-th spike. Note that the VS takes a value between 0 and 1. VS = 1 means that all the spikes occurred at a certain phase of the reference frequency f (i.e., perfect phase-locking) and a VS = 0 implies that the spike sequence has no phase preference.
Increases in the VS of the NM inputs lead to a gain in the AC component without altering the noise (Figures 4A,B). Theoretical calculations (Equations 2 and 4) and simulations (Figures 4C,D) indicate a linear relationship between AC and VS. This can be explained as follows: VS can be calculated as the absolute value of the Fourier component of the spike sequence at the stimulus frequency, normalized by the total number of spikes, regarding each spike as a delta function (Ashida et al., 2010). Within the regime where the synaptic and the membrane processes act as linear filters, the Fourier components of the resulting synaptic conductance and the membrane potential at the stimulus frequency are still linear with the VS. Consequently, the AC component, which is linearly related to the Fourier component, is linear to the VS. The prominent phase-locking property of NM fibers up to 8 kHz (Sullivan and Konishi, 1984; Köppl, 1997) is thus a fundamental component of high frequency ITD coding of the barn owl.
Figure 4. Dependence of the synaptic input in NL on the degree of phase-locking of presynaptic NM fibers. (A) Simulated traces of the model membrane potential. The number above each trace shows the vector strength of NM fibers. The traces show larger oscillations (higher AC amplitudes) as VS increases. (B) Power spectral densities of the four traces shown in (A). Peaks at the input frequency (4 kHz) and higher harmonics increase with increasing VS, while the other frequency components remain unchanged. (C) Dependence of the DC, AC, and noise amplitudes of the simulated synaptic input on VS. (D) Dependence of the AC and noise amplitudes of the simulated membrane potential on VS. Solid lines in (C) and (D) are obtained from analytical calculation without higher harmonics included. The dotted black line in (C) is obtained from analytical calculations with the second harmonic included. Vertical broken gray lines in (C) and (D) show the default parameter (VS = 0.6) used in our simulations.
Higher harmonics, which are regarded as part of the noise, also increase with VS (Figure 4B). Hence the noise amplitude in the input conductance is no longer constant for VS > 0.7 (Figure 4C). The noise component in the membrane potential, however, is almost independent of VS (Figure 4D), because the low-pass properties of the membrane effectively filter out higher harmonics of 8 kHz and above (see accompanying paper: Ashida et al., 2013, for how higher harmonics change with VS). It should be noted that even when all the input spikes are perfectly phase-locked (i.e., VS = 1), noise does not disappear (Figures 4C,D). This is due to the cycle-to-cycle variability. With our default parameters (M = 300 fibers spiking at λ0 = 500 Hz locked to fs = 4 kHz tonal stimulus), for example, the probability P that each fiber spikes at a certain stimulus cycle can be calculated as P = λ0/fs = 500/4000 = 0.125. Then the average number of inputs counted in a single stimulus cycle is MP = 37.5, with a standard deviation being . Therefore, even with perfect phase-locking, cycle-to-cycle variability of spike count would be on the order of 5.73/37.5 = 15%.
Synaptic Time Constant
Changes in the synaptic time scale (measured by the half peak width W of the unitary synaptic input) affect both the AC and noise responses (Figures 5A,B). Changes in W shift the filtering property of the synaptic input (Figure 6A). Slowing down the synaptic process reduces high frequency components. To conserve the 4 kHz signal component, the half peak width W should be equal to or smaller than the order of 0.1 ms (Figures 5C,D, 6B; see Funabiki et al., 2011, for related discussion). Decrease in W (i.e., speeding up the synaptic process) generally results in an increase in the AC and noise components of the input conductance (Figure 5C). As Equation 3 indicates, the noise component of the synaptic conductance blows up to infinity with W approaching to zero (Figure 5C). In the membrane potential (Figure 5D), however, the increase in noise with decreasing W is much slower because the increasing noise consists mostly of high frequency components (Figure 6A), which are filtered out by the low-pass effects of the membrane. As a result, the AC is more sensitive than the noise to the change in the width W of the unitary synaptic input (Figure 5D).
Figure 5. Dependence of synaptic input in NL on the synaptic time scale. (A) Simulated traces of the model membrane potential. The number above each trace shows the half peak widths W of the unitary synaptic input modeled by an alpha function (see Table 1). The traces show larger oscillations (higher AC amplitudes) as W decreases. (B) Power spectral densities of the four traces shown in (A). Higher frequency components decrease with increasing W. (C) Dependence of the DC, AC, and noise amplitudes of the simulated synaptic input on the synaptic time scale W. (D) Dependence of the AC and noise amplitudes of the simulated membrane potential on the synaptic time scale W. Solid lines in (C) and (D) are obtained from analytical calculations (Equations 1–5). Vertical broken gray lines in (C) and (D) show the default parameter (W = 0.1) used in our simulations.
Figure 6. Frequency properties of the synaptic filter. (A) Fourier transform |Fα(f)| of the synaptic filter α(t) (see Table 1). Normalized curves with W = 0.05, 0.10, 0.20, and 0.40 are shown. The synaptic filter becomes more likely to reduce high frequency components as the synaptic time scale W becomes smaller. (B) Synaptic filter at 4 kHz showing non-linear dependence on W. The vertical broken gray line in (B) shows the default parameter (W = 0.1) used in our simulations. (C) Comparison of synaptic filters at different sound stimulus frequencies fs (0.5, 1, 2, 4, and 8 kHz). Filter strength exceeds 0.5 (broken gray line), if and only if the inequality W < k/(2πfs) is satisfied. Note that this critical W-value is dependent on the frequency.
The half peak width W of the unitary synaptic input is linear to the time constant τ of the alpha function (Table 1). The AC component of the input conductance AG = (2rDG)/(1 + (2πfsτ)2) (see accompanying paper, Ashida et al., 2013, for detailed derivations). It is therefore necessary that W < k/(2πfs) in order to retain an effective amount of AC input (Figure 6C). This is roughly equivalent to W < 0.4/fs, with fs being the stimulus sound frequency. If W is larger than this criterion, the AC component is effectively filtered out by the synaptic process (e.g., Figure 5A, W = 0.4 ms). To reproduce the sound analog potential observed in vivo (AC = 1–2 mV, Funabiki et al., 2011), the synaptic time scale should be at least several times faster than that observed in chick high frequency NL cells in vitro (W = 0.2–0.3 ms, frequency range = 2.5–3.3 kHz, Kuba et al., 2005).
Implications for ITD Coding
In the preceding sections, we investigated the dependence of the sound analog membrane potential on the parameters of the phase-locked synaptic inputs. In this section, we examine how these parameters affect ITD coding. In the avian auditory brainstem, ITDs are compensated by the NM-axonal delay lines and computed by the coincidence detector neurons in NL (Carr and Konishi, 1990; Köppl and Carr, 2008). The phase difference δ between the bilateral synaptic inputs from NM to NL reflects the ITD (Figure 7A; see Funabiki et al., 2011, for more discussion). The AC component (4 kHz) of the total synaptic input changes periodically with the phase difference δ (Figures 7A,B), while the noise component is independent of δ (Figure 7B). The second harmonic (8 kHz), whose period is half of the main AC signal, is less than 0.1 mV and shows different dependence on the phase δ (Figure 7B).
Figure 7. Schematic representation of the binaural synaptic input to the NL neuron. (A) Summation of the synaptic input from ipsi- and contralateral NM fibers. The oscillation amplitude of the total synaptic input in NL is maximal when the two inputs arrive perfectly in-phase, while it becomes smaller when the two inputs are out of phase. For clarity, onset effects, higher harmonics, and noise components are not included in this schematic figure. (B) The oscillation amplitude of the total NM inputs. The amplitude of the AC component changes periodically with the phase difference δ of the two inputs according to the Equation: Amp(δ) = |A cos(δ/2)|, with A being the maximum AC amplitude (Ashida et al., 2007). The second harmonic also changes periodically with δ, but its period is half of that of the main signal (AC).
In the previous sections, and in our accompanying paper (Ashida et al., 2013), we used a single compartment NL model without a spike generator to focus on the fundamental properties of the sound analog membrane potential. In this section, however, we use a two-compartment NL model (Ashida et al., 2007; Funabiki et al., 2011) in which the cell body (soma) receives synaptic input and the axonal node generates spikes (Figure 8A), to study how sound analog potentials are converted to the output spike rate of the NL neuron. With the default input parameters, the spike rate of the model neuron is modulated periodically with the phase δ of the ipsi- and contralateral inputs, indicating a linear conversion of the AC signals into spike rate (Figure 8B; see Funabiki et al., 2011, for detailed discussion on the linear conversion). Generally, the model neuron shows the highest spike rate when the two inputs arrive perfectly in-phase (δ = 0 in Figure 8B). We hereafter refer to this maximum spike rate as the “in-phase rate.” When the inputs arrive in perfect anti-phase (δ = ±π in Figure 8B), the spike rate becomes lowest, which we call the “out-of-phase rate.” Note that the out-of-phase rate is the spike rate driven primarily by noise, while the in-phase rate is driven by both AC and noise components.
Figure 8. Reponses of the two-compartment NL model to 4 kHz input. (A) Two-compartment NL neuron model, with a soma and a node interconnected by an axonal conductance. The somatic compartment has the same amount of leak and low-voltage-activated potassium (KLVA) conductances as the single compartment model. In addition to the leak and KLVA, the nodal compartment has Na and high-voltage-activated potassium (KHVA) conductances to generate spikes. (B) Spike rate of the model neuron plotted against the phase difference δ between ipsi- and contralateral model NM inputs. We define the spike rate with δ = 0 as the “in-phase rate,” and δ = ±π as the “out-of-phase rate.” In-phase and out-of-phase rates are plotted against the mean spike rate of NM fibers (C), the total number of NM fibers (D), the degree of phase-locking of the NM fibers (E), and the half peak width W of the unitary synaptic input conductance (F). Vertical broken gray lines in (C–F) show the default parameters (λ0 = 500 Hz, M = 300 fibers, r = 0.6, W = 0.1 ms) used in our simulations.
Varying the mean spike rate (λ0) of NM, the number (M) of NM fibers, the degree (r) of phase-locking of these fibers, and the synaptic time scale (W), we calculated the in-phase and out-of-phase rates of the model neuron. The out-of-phase rate increases with both decreasing NM spike rate (Figure 8C) and decreasing numbers of NM inputs (Figure 8D), because of the increasing noise (Figures 2D, 3D). The in-phase-rate is higher than the out-of-phase rate by a few hundred Hz (Figures 8C,D), provided that the number of NM fibers exceeds 100 and their spike rates are over a few hundred Hz. Changes in VS do not alter the out-of-phase rate (Figure 8E) because the noise component without higher harmonics is independent of VS (Figure 4D). The out-of-phase curve is flat even near VS = 1 (Figure 8E), indicating a small contribution of the second and higher harmonics in our 4 kHz NL model (Figure 7B). The in-phase rate increases almost linearly with VS (Figure 8E), in agreement with the linear conversion of the AC component to spike rate found in NL neurons in vivo (Funabiki et al., 2011). Since both the AC and noise components increase with decreasing synaptic time scale W (Figure 5D), both the in-phase and out-of-phase rates also increase with decreasing W (Figure 8F). The increase in the in-phase-rate is, however, more prominent, because W is more likely to affect the AC amplitude than the noise amplitude (Figure 5D). Therefore, the modulation depth, which is the difference between the in-phase and out-of-phase rates, becomes higher for smaller W (i.e., faster synaptic inputs). In our simulation, W needed to be equal or smaller than 0.1 milliseconds for the 4 kHz NL model neuron to obtain modulation depths of a few hundred Hz, consistent with measured spike rates in barn owl's NL (Peña et al., 1996; Funabiki et al., 2011).
Phase-locking, or the generation of action potentials at a certain phase of the reference signal, contributes to fine temporal information coding in the auditory system (Oertel, 1999; Ashida et al., 2010; Brette, 2012). In the owl's auditory brainstem, phase-locked synaptic input sequences from the NM axons are filtered by synaptic and membrane processes, leading to oscillatory membrane potentials in NL (Gerstner et al., 1996; Kempter et al., 1998; Ashida et al., 2007; Funabiki et al., 2011). The phase differences between inputs from ipsi- and contralateral NM, which are the consequence of ITD, lead to the periodic changes in the AC component of the binaural input (Figure 7). The NL neuron acts as a linear converter of AC signals into spike rate (Funabiki et al., 2011), and thus changes spike rate periodically with ITD or with the input phase difference (Figure 8B). Similar oscillatory membrane potentials were found in auditory hair cells (Russell and Sellick, 1978; Dallos, 1985). Unlike auditory hair cells, however, NL neurons are located a few synapses away from the cochlea. The sound analog potential in NL is therefore affected by various presynaptic, synaptic, and postsynaptic factors.
In this paper, we focused particularly on the sound analog membrane potential induced by pure tones. Since the NL neuron is narrowly tuned to frequency (Peña et al., 2001), the input to each NL neuron is expected to be dominated by the AC component at or near its characteristic frequency even if the stimulus sound is broadband. Nevertheless, for more complex stimuli than simple tones, multiple frequency components might affect ITD coding in more complicated fashion. For example, different frequency components show different dependence on the synaptic time constant (Figure 6C). Thus, the optimal synaptic time scale for pure tones could differ from that for broadband stimuli. More modeling and physiological investigation will be necessary to examine how complex stimuli is presented and processed in owls' NL.
In our series of numerical simulations, we fixed the average input level (DC) in order to focus on the effects of the AC and noise components on the ITD coding performance of the model NL neuron. This assumption is based on the following observations: (1) the DC component is irrelevant to ITD coding (Funabiki et al., 2011); (2) the NL neuron increases its spike rate even with a small amount of external current injection (Funabiki et al., 2011); and (3) the DC level should be carefully chosen so that the model neuron shows good sensitivity to AC signals (Ashida et al., 2007). For these reasons, we presumed that the DC input of the NL neuron in vivo should be optimized to efficiently compute ITD-dependent AC signals. In our simulations we determined the DC level so that the difference between the in-phase and out-of-phase spike rates exceeded 180 spikes/s (Figure 8B), as was observed in vivo (Peña et al., 1996). Large increases or decreases in the constant DC should lead to over- or under-excitability of the neuron, resulting in the degradation of the overall ITD coding performance (see also Funabiki et al., 2011, for related discussion on the effects of DC).
Changes in both the number (Figure 2) and spike rate (Figure 3) of presynaptic NM fibers alter the noise component but not the AC signal. In the ITD coding pathway of the barn owl, an NM neuron receives 1–4 auditory nerve terminals (Carr and Boudreau, 1991), while an NL neuron is estimated to receive a few hundred NM afferents (Carr and Boudreau, 1993) and a neuron in the anterior part of the dorsal lateral lemniscus (LLDa) is estimated to receive 2–10 NL inputs (Fischer and Konishi, 2008). The primary role of the NM neuron is to convey precise temporal information from auditory nerves to NL (Sullivan and Konishi, 1984; Köppl, 1997), whereas coincidence detector neurons in NL compare bilateral NM inputs using the sound analog potentials, whose amplitudes systematically change with ITD (Funabiki et al., 2011). The LLDa provides finer ITD tuning than NL by reducing noise (Fischer and Konishi, 2008). The optimal convergence number of input fibers for 4 kHz ITD computation in NL is predicted to be on the order of a few hundred (Figure 8D), whereas the optimal number for NM and LLDa, which have different computational roles, should be smaller than NL.
In contrast to the owl's NL neurons, neurons in the gerbil's MSO, the mammalian counterpart of NL, have recently been reported to receive less than 10 excitatory synaptic inputs (Couchman et al., 2010). Thus, the membrane potential of the MSO neuron may not show the clear sinusoidal waveforms that characterize responses in the owl's NL. Whereas the NL neuron receives only slow GABAergic inhibition (Burger et al., 2011; Coleman et al., 2011), the gerbil's MSO neurons receive fast glycinergic inhibition (Magnusson et al., 2005), which plays an essential role in the ITD coding in the MSO (Brand et al., 2002; Pecka et al., 2008). The difference in computational strategies in MSO and NL may include differences in the number of inputs in these neurons.
The degree of phase-locking, as measured by VS, has a linear relationship to the resulting AC signal (Figure 4), but does not affect the noise, which originates from the cycle-to-cycle variability of input spike counts. The in-phase rate, therefore, monotonically increases with VS, while the out-of-phase rate remains almost independent of VS (Figure 8E). Thus, the prominent phase-locking ability observed in owl's NM (Köppl, 1997) leads directly to the high computational efficiency in NL. In the accompanying paper (Ashida et al., 2013), we pointed out that perfect phase-locking may not always be beneficial to ITD coding because of the detrimental effects of higher harmonics. In our 4 kHz model neuron, however, the amplitude of the second harmonic is less than 10% of the signal amplitude (Figure 7B) because of the low-pass property of the membrane (compare Figures 4C,D). Thus, the negative contribution higher harmonics should be limited to NL neurons with best frequencies of 2 kHz or below.
Short synaptic time constants are one of the most important features of high-frequency coincidence detection. Auditory coincidence detector neurons in the chicken NL (Kuba et al., 2005) and gerbil MSO (Couchman et al., 2010) receive extremely fast synaptic inputs with time constants of far below a millisecond. The synaptic time scale, measured by the half peak width of unitary synaptic input (W), affects both AC and noise components. Although these components of the input conductance increase with decreasing synaptic time constant (Figure 5C), the effect of the increasing noise on the membrane potential is smaller than that of the AC component (Figure 5D), because the membrane filter cuts out high frequency noise that appears with decreasing W (Figure 6A). In consequence, the modulation between in-phase and out-of phase rates becomes higher for smaller W (Figure 8F). In order to obtain a sufficient AC signal, W needs to be smaller than 0.4/fs (ms), with fs being the signal frequency (in kHz). This implies that the higher the cell's signal frequency, the faster the synaptic input process should be (Slee et al., 2010). Previous studies of NL (Grau-Serrat et al., 2003) and mammalian MSO (Mathews et al., 2010; Khurana et al., 2011) identified various mechanisms underlying the submillisecond accuracy of synaptic inputs. Jercog et al. (2010) showed that synaptic inputs from the ipsilateral side are faster than those from the contralateral side, contributing to the fine ITD tuning of MSO neurons. Further investigation is necessary to determine how fast the actual owl's NL synapses are in vivo.
Membrane time constants of about 0.1 ms or less are essential for high frequency signal processing in NL (Gerstner et al., 1996; Kempter et al., 1998). Impedance analyses (Gutfreund et al., 1995; Ashida et al., 2007) showed that a neuron that handles high frequency signals needs to be both leaky (i.e., having a low input resistance) and electrotonically compact (i.e., having a small capacitance). Short dendrites (Carr and Konishi, 1990; Carr and Boudreau, 1993; Agmon-Snir et al., 1998; Kuokkanen et al., 2010) and large KLVA conductances (Kuba et al., 2005) are consistent with these requirements. It should be noted that even with membrane properties that allow high frequency signals up to several kHz, the membrane still acts as a low-pass filter (see accompanying paper: Ashida et al., 2013) and efficiently eliminates higher harmonics (Figures 4C,D, 5C,D, 7B). Recent simulation results suggest that the kinetic properties of the KLVA, such as increased activation with depolarization, may improve ITD computation in the gerbil's MSO neuron (Svirskis et al., 2002; Day et al., 2008; Gai et al., 2009; Jercog et al., 2010; Mathews et al., 2010). Dynamic interplay between KLVA and hyperpolarization-activated cation currents has been suggested to stabilize synaptic inputs in the MSO (Khurana et al., 2011). How these conductances dynamically affect the computation of high frequency AC signals in owl's NL remains to be elucidated.
In the owl's NL, ITD tuning has been found up to about 7 kHz (Carr and Konishi, 1990). In our preliminary simulations, however, it was hard to reproduce ITD tuning at 7 kHz (data not shown), primarily due to the decrease in AC amplitude with increasing frequency (see also Figure 4 of our accompanying paper: Ashida et al., 2013). Therefore, additional mechanisms and finer tuning should be incorporated in the actual ITD coding of the barn owl. Given that VS decreases with frequency (Köppl, 1997), high-frequency NL neurons would require more inputs to reduce noise (Figure 3D), as well as faster synaptic inputs (Figure 5D) and smaller membrane time constants (Kempter et al., 1998) to preserve AC signals. In our modeling, we focused on the NM-NL circuit, where only excitatory synaptic inputs were considered. In the avian brainstem, however, both NM and NL neurons receive GABAergic inhibitory inputs from the superior olivary nucleus (Carr et al., 1989; Burger et al., 2005; Coleman et al., 2011; see Burger et al., 2011, for a recent review), which play an important role in refining ITD tuning in NL (Fujita and Konishi, 1991; Funabiki et al., 1998; Yang et al., 1999; Nishino et al., 2008). Inhibitory inputs might contribute in stabilizing the formation of sound analog potential.
In addition, neuronal activity in the Mauthner cell cap (Korn and Faber, 1975), hippocampus (Radman et al., 2007), and in the cortex (Anastassiou et al., 2011) can be affected by inhomogeneous distributions of the extracellular field potential. This effect, called the ephaptic coupling, provides a strong basis of synchronous spike activities in these brain areas. In the owl's NL, the neurophonic, or the extracellular field potential, is strongly correlated with tonal stimuli whose amplitudes are in the millivolt range (see Kuokkanen et al., 2010; and references therein). If ephaptic coupling occurs between the neurophonic and NL neurons, it could be an additional mechanism of mediating the sound analog potential.
Materials and Methods
The detailed modeling procedure and its analytical consideration are provided in our previous (Ashida et al., 2007; Funabiki et al., 2011) and accompanying (Ashida et al., 2013) papers. Here we summarize the fundamental equations of our model NM-NL system.
Modeling Phase-Locked Synaptic Input to NL
In order to model the phase-locked NM activity, we used an inhomogeneous Poisson process with a time-dependent intensity function λ(t) = 2πλ0 pκ(2πfst), where λ0 is the average intensity, pκ is the von-Mises distribution function (parameterized by the concentration parameter κ) and fs is the stimulus sound frequency. The model NL neuron receives inputs from M fibers that are phase-locked to the stimulus tone with a VS of r = I1(κ)/I0(κ), with In being the Modified Bessel function of order n. In order to model the unitary synaptic conductance in NL, we used an alpha-function α(t) = (Ht/τ) exp(1 − t/τ), with H being the peak height and τ being the time constant. Note that the half peak width W of the alpha function is linearly related to τ as W = 2.446τ All the synaptic conductance input is linearly summed and delivered to the soma of the model NL neuron. The model equations and parameters used in this paper are summarized in Table 1.
In the present paper, we changed the average spike rate of NM fibers (λ0), the number of converging fibers to the model NL neuron (M), the average degree of phase-locking (r) of the NM fibers, and the synaptic time scale (W). In numerical simulations where λ0, M, or W was changed, the peak height H of the unitary synaptic input was re-scaled to conserve the total input conductance. For example, when the mean spike rate of NM fibers was set at twice (M = 1000 Hz) the default rate (M = 500 Hz), the peak height H was reduced to the half (0.65 nS) of the default value (1.3 nS).
Single Compartment NL Neuron Model
In order to examine the parameter dependence of the AC and noise components of the synaptic input and the membrane potential of the model NL neuron, we used a conductance-based single compartment model. The dynamics of the membrane potential of the soma VS is described by a first-order differential Equation:
where CS is the membrane capacitance, ILS = gLS(EL − VS) is the leak current, IKLVAS = ḡKLVAS d(VS, t)(EK − VS) is the KLVA current, IaxonS is the axonal current, and Isyn = gsyn(Esyn − VS) is the model synaptic input. For the single compartment model, the axonal current was fixed to zero. The activation variable d(V, t) of the KLVA conductance obeys the first-order differential Equation:
Membrane properties, including the maximum KLVA conductance ḡKLVA, were fixed in our simulations.
ITD Dependence and Two-Compartment NL Neuron Model
A two-compartment NL neuron model was used to examine the parameter dependence of output spike rate (see Ashida et al., 2007; for detailed explanations for the two-compartment model). The model consists of the soma and the node connected by the axonal resistance (Figure 8A). The model equations and parameters are the same as in our previous study (Funabiki et al., 2011). The somatic variables and parameters are the same as those used for the single compartment model. Sodium and high-voltage-activated potassium (KHVA) conductances were introduced in the nodal compartment to generate spikes. The membrane potential of the node VN is described as:
where ILN = gLN(EL − VN) is the leak current, IKLVAN = ḡKLVAN d(VN, t)(EK − VN) and IKHVAN = ḡKLVAN n(VN, t)(EK − VN) are, respectively, the KLVA and KHVA currents, INaN = ḡNaN m(VN, t)h(VN, t)(ENa − VN) is the fast sodium current, and IaxonS = −IaxonN = gaxon(VN − VS) is the axonal current. The KHVA conductance has only the activation variable n(V, t), while the sodium conductance has both activation m(V, t) and inactivation h(V, t) variables. These variables obey the equation , where x stands for n, m, or h. The model equations and parameters used in this paper are summarized in Table 2.
In constructing synaptic inputs, half the NM fibers were assumed to be from the ipsilateral side, and the remaining half of the NM fibers from the contralateral side. NM fibers from each side were assumed to be phase-locked to generate the oscillatory conductance, and the NM inputs from the two sides were summed with a phase difference δ (see Figure 7A). The spike rate of the model neuron with δ = 0 was called the “in-phase rate,” while that with δ = π was referred to as the “out-of-phase rate” (see Figure 8A).
Analytical results obtained in the accompanying paper (Ashida et al., 2013) are summarized in Table 1. Equations 1–5 describe the parameter dependence of the average (DG), the signal part (AG), and the noise part (NG) of the synaptic conductance input, as well as the signal (AV) and noise (NV) components of the membrane potential.
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.
We thank Masakazu Konishi and Kousuke Abe for valuable comments and discussion. This work was supported by NIH DC00436 to Catherine E. Carr, by NIH P30 DC04664 to the University of Maryland Center for the Comparative and Evolutionary Biology of Hearing, by fellowships from the Alexander von Humboldt Foundation and the Hanse-Wissenschaftskolleg to Go Ashida and Catherine E. Carr, and by a postdoctoral fellowship from JSPS and Grant-in-Aid for Scientific Research (B) to Kazuo Funabiki. Go Ashida is also supported by the Cluster of Excellence “Hearing4all” at the University of Oldenburg.
ITD, interaural time difference; NM, nucleus magnocellularis; NL, nucleus laminaris; VS, vector strength; EPSG, excitatory post synaptic conductance; KHVA, high voltage activated potassium; KLVA, low voltage activated potassium; LLDa, anterior part of the dorsal lateral lemniscus; MSO, medial superior olive.
Ashida, G., Abe, K., Funabiki, K., and Konishi, M. (2007). Passive soma facilitates submillisecond coincidence detection in the owl's auditory system. J. Neurophysiol. 97, 2267–2282. doi: 10.1152/jn.00399.2006
Ashida, G., Funabiki, K., and Carr, C. E. (2013). Theoretical foundations of the sound analogue membrane potential that underlies coincidence detection in the barn owl. Front. Comput. Neurosci. 7: 151. doi: 10.3389/fncom.2013.00151
Ashida, G., Wagner, H., and Carr, C. E. (2010). “Processing of phase-locked spikes and periodic signals,” in Analysis of Parallel Spike Trains, eds S. Grün and S. Rotter (New York, NY: Springer), 59–74.
Burger, R. M., Cramer, K. S., Pfeiffer, J. D., and Rubel, E. W. (2005). Avian superior olivary nucleus provides divergent inhibitory input to parallel auditory pathways. J. Comp. Neurol. 481, 6–18. doi: 10.1002/cne.20334
Burger, R. M., Fukui, I., Ohmori, H., and Rubel, E. W. (2011). Inhibition in the balance: binaurally coupled inhibitory feedback in sound localization circuitry. J. Neurophysiol. 106, 4–14. doi: 10.1152/jn.00205.2011
Carr, C. E., and Boudreau, R. E. (1993). Organization of the nucleus magnocellularis and the nucleus laminaris in the barn owl: encoding and measuring interaural time differences. J. Comp. Neurol. 334, 337–355. doi: 10.1002/cne.903340302
Coleman, W. L., Fischl, M. J., Weimann, S. R., and Burger, R. M. (2011). GABAergic and glycinergic inhibition modulate monaural auditory response properties in the avian superior olivary nucleus. J. Neurophysiol. 105, 2405–2420. doi: 10.1152/jn.01088.2010
Couchman, K., Grothe, B., and Felmy, F. (2010). Medial superior olivary neurons receive surprisingly few excitatory and inhibitory inputs with balanced strength and short-term dynamics. J. Neurosci. 30, 17111–17121. doi: 10.1523/JNEUROSCI.1760-10.2010
Day, M. L., Doiron, B., and Rinzel, J. (2008). Subthreshold K+ channel dynamics interact with stimulus spectrum to influence temporal coding in an auditory brain stem model. J. Neurophysiol. 99, 534–544. doi: 10.1152/jn.00326.2007
Funabiki, K., Ashida, G., and Konishi, M. (2011). Computation of interaural time difference in the owl's coincidence detector neurons. J. Neurosci. 31, 15245–15256. doi: 10.1523/JNEUROSCI.2127-11.2011
Funabiki, K., Koyano, K., and Ohmori, H. (1998). The role of GABAergic inputs for coincidence detection in the neurons of nucleus laminaris of the chick. J. Physiol. (Lond.) 508, 851–869. doi: 10.1111/j.1469-7793.1998.851bp.x
Gai, Y., Doiron, B., Kotak, V., and Rinzel, J. (2009). Noise-gated encoding of slow inputs by auditory brain stem neurons with a low-threshold K+ current. J. Neurophysiol. 102, 3447–3460. doi: 10.1152/jn.00538.2009
Goldberg, J. M., and Brown, P. B. (1969). Response of binaural neurons of dog superior olivary complex to dichotic tonal stimuli: some physiological mechanisms of sound localization. J. Neurophysiol. 32, 613–636.
Jercog, P. E., Svirskis, G., Kotak, V. C., Sanes, D. H., and Rinzel, J. (2010). Asymmetric excitatory synaptic dynamics underlie interaural time difference processing in the auditory system. PLoS Biol. 8: e1000406. doi: 10.1371/journal.pbio.1000406
Kempter, R., Gerstner, W., van Hemmen, J. L., and Wagner, H. (1998). Extracting oscillations: neuronal coincidence detection with noisy periodic spike input. Neural Comput. 10, 1987–2017. doi: 10.1162/089976698300016945
Khurana, S., Remme, M. W. H., Rinzel, J., and Golding, N. L. (2011). Dynamic interaction of Ih and IK−LVA during trains of synaptic potential in principal neurons of the medial superior olive. J. Neurosci. 31, 8936–8947. doi: 10.1523/JNEUROSCI.1079-11.2011
Klump, G. M. (2000). “Sound localization in birds,” in Comparative Hearing: Birds and Reptiles, eds R. J. Dooling, R. R. Fay, and A. N. Popper (New York, NY: Springer), 249–307. doi: 10.1007/978-1-4612-1182-2_6
Koyano, K., Funabiki, K., and Ohmori, H. (1996). Voltage-gated ionic currents and their roles in timing coding in auditory neurons of the nucleus magnocellularis of the chick. Neurosci. Res. 26, 29–45.
Kuba, H., Yamada, R., Fukui, I., and Ohmori, H. (2005). Tonotopic specialization of auditory coincidence detection in nucleus laminaris of the chick. J. Neurosci. 25, 1924–1934. doi: 10.1523/JNEUROSCI.4428-04.2005
Kuokkanen, P. T., Wagner, H., Ashida, G., Carr, C. E., and Kempter, R. (2010). Theoretical and experimental signal-to-noise ratio of the neurophonic potential in nucleus laminaris of the barn owl (Tyto alba). J. Neurophysiol. 104, 2274–2290. doi: 10.1152/jn.00395.2010
Magnusson, A. K., Kapfer, C., Grothe, B., and Koch, U. (2005) Maturation of glycinergic inhibition in the gerbil medial superior olive after hearing onset. J. Physiol. 568, 497–512. doi: 10.1113/jphysiol.2005.094763
Mathews, P. J., Jercog, P. E., Rinzel, J., Scott, L. L., and Golding, N. L. (2010). Control of submillisecond synaptic timing in binaural coincidence detectors by KV1 channels. Nat. Neurosci. 13, 601–611. doi: 10.1038/nn.2530
Nishino, E., Yamada, R., Kuba, H., Furuta, T., Kaneko, T., and Ohmori, H. (2008). Sound-intensity-dependent compensation for the small interaural time difference cue for sound source localization. J. Neurosci. 28, 7153–7164. doi: 10.1523/JNEUROSCI.4398-07.2008
Pecka, M., Brand, A., Behrend, O., and Grothe, B. (2008) Interaural time difference processing in the mammalian medial superior olive: the role of glycinergic inhibition. J. Neurosci. 28, 6914–6925. doi: 10.1523/JNEUROSCI.1660-08.2008
Radman, T., Su, Y., An, J. H., Parra, L. C., and Bikson, M. (2007) Spike timing amplifies the effect of electric fields on neurons: implications for endogenous field effects. J. Neurosci. 27, 3030–3036. doi: 10.1523/JNEUROSCI.0095-07.2007
Rothman, J. S., and Manis, P. B. (2003). The roles potassium currents play in regulating the electrical activity of ventral cochlear nucleus neurons. J. Neurophysiol. 89, 3097–3113. doi: 10.1152/jn.00127.2002
Scott, L. L., Mathews, P. J., and Golding, N. L. (2010) Perisomatic voltage-gated sodium channels actively maintain linear synaptic integration in principal neurons of the medial superior olive. J. Neurosci. 30, 2039–2050. doi: 10.1523/JNEUROSCI.2385-09.2010
Svirskis, G., Kotak, V., Sanes, D. H., and Rinzel, J. (2002). Enhancement of signal-to-noise ratio and phase locking for small inputs by a low-threshold outward current in auditory neurons. J. Neurosci. 22, 11019–11025.
Yamada, R., Kuba, H., Ishii, T. M., and Ohmori, H. (2005). Hyperpolarization-activated cyclic nucleotide-gated cation channels regulate auditory coincidence detection in nucleus laminaris of the chick. J. Neurosci. 25, 8867–8877. doi: 10.1523/JNEUROSCI.2541-05.2005
Keywords: phase-locking, sound localization, auditory brainstem, periodic signals, oscillation, owl
Citation: Ashida G, Funabiki K and Carr CE (2013) Biophysical basis of the sound analog membrane potential that underlies coincidence detection in the barn owl. Front. Comput. Neurosci. 7:102. doi: 10.3389/fncom.2013.00102
Received: 04 December 2012; Accepted: 07 July 2013;
Published online: 08 November 2013.
Edited by:Misha Tsodyks, Weizmann Institute of Science, Israel
Reviewed by:Maoz Shamir, Boston University, USA
G. Björn Christianson, University College London, UK
Copyright © 2013 Ashida, Funabiki and Carr. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Go Ashida, Cluster of Excellence “Hearing4all,” AG Computational Neuroscience, Fakultät VI, Carl von Ossietzky Universität Oldenburg, Carl-von-Ossietzky-Str., 9-11, D-26129 Oldenburg, Germany e-mail: email@example.com