Signatures of Somatic Inhibition and Dendritic Excitation in Auditory Brainstem Field Potentials

Joshua H. Goldwyn; Myles McLaughlin; Eric Verschooten; Philip X. Joris; John Rinzel

doi:10.1523/JNEUROSCI.0600-17.2017

Abstract

Extracellular voltage recordings (V_e; field potentials) provide an accessible view of in vivo neural activity, but proper interpretation of field potentials is a long-standing challenge. Computational modeling can aid in identifying neural generators of field potentials. In the auditory brainstem of cats, spatial patterns of sound-evoked V_e can resemble, strikingly, V_e generated by current dipoles. Previously, we developed a biophysically-based model of a binaural brainstem nucleus, the medial superior olive (MSO), that accounts qualitatively for observed dipole-like V_e patterns in sustained responses to monaural tones with frequencies >∼1000 Hz (Goldwyn et al., 2014). We have observed, however, that V_e patterns in cats of both sexes appear more monopole-like for lower-frequency tones. Here, we enhance our theory to accurately reproduce dipole and non-dipole features of V_e responses to monaural tones with frequencies ranging from 600 to 1800 Hz. By applying our model to data, we estimate time courses of paired input currents to MSO neurons. We interpret these inputs as dendrite-targeting excitation and soma-targeting inhibition (the latter contributes non-dipole-like features to V_e responses). Aspects of inferred inputs are consistent with synaptic inputs to MSO neurons including the tendencies of inhibitory inputs to attenuate in response to high-frequency tones and to precede excitatory inputs. Importantly, our updated theory can be tested experimentally by blocking synaptic inputs. MSO neurons perform a critical role in sound localization and binaural hearing. By solving an inverse problem to uncover synaptic inputs from V_e patterns we provide a new perspective on MSO physiology.

SIGNIFICANCE STATEMENT Extracellular voltages (field potentials) are a common measure of brain activity. Ideally, one could infer from these data the activity of neurons and synapses that generate field potentials, but this “inverse problem” is not easily solved. We study brainstem field potentials in the region of the medial superior olive (MSO); a critical center in the auditory pathway. These field potentials exhibit distinctive spatial and temporal patterns in response to pure tone sounds. We use mathematical modeling in combination with physiological and anatomical knowledge of MSO neurons to plausibly explain how dendrite-targeting excitation and soma-targeting inhibition generate these field potentials. Inferring putative synaptic currents from field potentials advances our ability to study neural processing of sound in the MSO.

Introduction

Extracellular voltage (V_e) recordings, also known as field potentials, can provide a view of neural activity across multiple locations in a brain area. The biophysics of synaptic and neural generators of field potentials are well understood (for review, see Buzsáki et al., 2012), but interpreting field potentials remains a challenge. Field potentials are shaped in complicated ways by volume conduction (Kajikawa and Schroeder, 2015), morphology of neurons (Pettersen and Einevoll, 2008), their spatial arrangement (Klee and Rall, 1977), and the statistics of synaptic inputs (Lindén et al., 2011), among other factors.

Inferring underlying neural activity from field potentials has been accomplished in studies that skillfully wove together mathematical modeling with knowledge of neurophysiology and anatomy (Rall and Shepherd, 1968: olfactory bulb; Nicholson and Llinas, 1971: cerebellum; Kuokkanen et al., 2010: nucleus laminaris; and Fernández-Ruiz et al., 2013: dentate gyrus). We follow this approach to analyze V_e recorded in the auditory brainstem of cats in response to monaural pure tones. We adopt a standard terminology and refer to these field potentials as the auditory neurophonic (Weinberger et al., 1970).

A motivation for studying the auditory neurophonic has been the difficulty of accessing individual medial superior olive (MSO) neurons with in vivo preparations. Traditional extracellular techniques in MSO have yielded small samples (Goldberg and Brown, 1969; Guinan et al., 1972a; Yin and Chan, 1990; Day and Semple, 2011) and are contaminated by the neurophonic due to its large amplitude and the small size of MSO action potentials (Scott et al., 2007). Examples of single-unit patch-clamp recordings (Franken et al., 2015) and juxtacellular recordings (van der Heijden et al., 2013; Plauška et al., 2016) have recently been reported, but the technical challenges of these methods remain. Another approach has been to record from the axons of MSO neurons at locations isolated from the neurophonic (Bremen and Joris, 2013), but with this approach there is less certainty regarding the location and identity of recorded neurons. Given its large size and ease of recording, the neurophonic provides a useful view of the MSO with possible applicability across species (including humans).

Early investigators of the neurophonic understood that postsynaptic currents in MSO neurons likely are the dominant generators of these field potentials (Galambos et al., 1959; Biedenbach and Freeman, 1964; Tsuchitani and Boudreau, 1964). A second motivation to study the neurophonic is, therefore, to observe synaptic inputs to MSO because they may be difficult to access with other techniques. MSO neurons have a simple and stereotyped bipolar morphology (Stotler, 1953; Rautenberg et al., 2009), the dendritic arbors of neighboring neurons extend approximately in parallel, and synaptic excitation primarily targets these dendrites (Cant and Casseday, 1986; Couchman et al., 2012). For sounds played to one ear, there is a flow of excitatory current into dendrites on one side of the MSO (current sink) and a corresponding outflow of current through the soma and opposite dendrite (current source). These facts suggest that MSO neurons act (qualitatively) as current dipoles. Our previous work has supported this view (Mc Laughlin et al., 2010; Goldwyn et al., 2014).

This dipole theory is a helpful but incomplete account of V_e responses. In particular, we observed that sustained V_e responses to low-frequency tones are more coherent in space (more “monopole-like”) than the “dipole-like” responses to higher-frequency tones (∼1000 Hz and above). We argue that the “dipole theory” is incomplete because it does not consider inhibitory inputs that are known to converge with excitatory inputs on MSO neurons. The role of inhibition in MSO processing is a matter of debate. Some have surmised that inhibition precedes excitation and thereby delays the effect of excitatory inputs to MSO (Brand et al., 2002; Pecka et al., 2008; Myoga et al., 2014), but juxtacellular and intracellular recordings do not support this view (van der Heijden et al., 2013; Franken et al., 2015). If inhibition could be studied via field potentials, such recordings may enable the formulation and evaluation of alternative hypotheses.

We demonstrate that it may be possible to infer essential features of synaptic inputs to MSO from brainstem field potentials. We construct a mathematical model of an MSO neuron that receives both dendrite-targeting excitation and soma-targeting inhibition and use field potentials to determine the time courses of these paired inputs. Simulations accurately reproduce the diverse and frequency-dependent spatiotemporal V_e patterns observed in experiments. Furthermore, aspects of the inferred synaptic currents are consistent with properties of excitatory and inhibitory inputs to MSO. Our model-based analysis of the neurophonic provides a route to “see” excitatory and inhibitory inputs to MSO in field potentials. Importantly, we use the model to show that soma-targeting inhibition contributes a monopole-like component to neurophonic responses that is necessary to accurately model onset responses and sustained portions of responses to low-frequency tones. This prediction can be tested in experiments. In particular, our theory predicts that blocking inhibition pharmacologically (Brand et al., 2002; Pecka et al., 2008; Jercog et al., 2010; Roberts et al., 2013; Myoga et al., 2014; Franken et al., 2015) would remove monopole-like features from neurophonic responses.

Materials and Methods

Neurophonic recordings.

We analyzed extracellular voltage recordings from adult cats of both sexes anesthetized with sodium pentobarbital at doses titrated to achieve an areflexive state. These experiments were first described by Mc Laughlin et al. (2010). We summarize our methods here and refer readers to that publication for further details.

We recorded extracellular voltage signals using a planar array of five electrodes (quartz-platinum, 2–4 MΩ). We advanced the electrode array through the auditory brainstem in steps of 50 or 100 μm using a TREC microdrive. We selected an angle of approach (15° to 30° mediolateral relative to the midsagittal plane) with the intention that the electrode track would run parallel to the orientation of dendrites in the MSO (the short axis of the MSO). All signals were filtered (10 Hz to 10 kHz) by the TREC headstages. The 10 Hz lower cutoff eliminates the “DC” component of the signal, a point we discuss in more detail below when comparing onset and ongoing portions of neurophonic responses.

We presented monaural tone bursts over a range of stimulus frequencies (100–2500 Hz for contralateral inputs). We analyzed a subset of these data, focusing on the range of frequencies over which neurophonic responses were largest and most reliably evoked (600–1800 Hz, for contralateral stimulation). Tone bursts to the ipsilateral side were, in most cases, presented at a slightly higher frequency (i.e., 610 Hz instead of 600 Hz, 710 Hz instead of 700 Hz, etc.). For ease of presentation, we ignore this small frequency difference when discussing results in text and figures. The stimulus level was 70 dB SPL for all recording sessions analyzed in this study.

In two recording sessions, for two frequencies (700 and 1300 Hz), we repeated the same stimulus 100 times and averaged the responses. In these two sessions the tone duration was 50 ms. In all other recording sessions, each sound stimulus was presented once for a duration of 2000 ms.

In our exploration of these data, we identified a subset of recordings that exhibited qualitatively similar spatiotemporal patterns. Our goal for this study was to provide an account of these characteristic response patterns. As such, in this paper we present data from five recording sessions (4 animals) and exclude data from nine recording sessions (5 animals). In each recording session, we obtained five simultaneous recordings (using the 5 electrode array). We present data from one electrode per recording session and do not make comparisons across electrodes in this study. Data from these or similar recording sessions were presented previously by Mc Laughlin et al. (2010) and Goldwyn et al. (2014). In Goldwyn et al. (2014) we limited ourselves to analyzing responses to tones with frequency 1000 Hz or higher because these responses exhibited dipole-like features most strongly.

In some cases, we obtained histological information that allowed us to identify the location of the MSO along the electrode track (Table 1; Mc Laughlin et al. (2010)).

View this table:

Table 1.

Spatial parameter values used to model neurophonic responses

Spatiotemporal V_e patterns.

The basic unit of analysis in this paper is the spatiotemporal pattern of V_e. To construct these patterns, we aggregated the (non-simultaneous) V_e responses that were obtained at multiple locations in the brainstem by advancing the electrode microdrive in steps of 50 or 100 μm through the brainstem. A prominent feature of neurophonic responses to pure tone stimuli is a temporal oscillation at the tone frequency. One period of oscillation (the inverse of the tone frequency; 1.67 ms for 600 Hz tone, 1.43 ms for 700 Hz tone, etc.) sets a natural time scale for analysis. The spatiotemporal V_e patterns that we constructed and analyzed have spatial extents of ∼2000 μm (smallest is 1500 μm, largest is 2500 μm) and temporal extents of one period of oscillation (time in milliseconds depends on the frequency of the pure tone stimulus).

We distinguished between “onset” and “ongoing” portions of neurophonic responses. The ongoing portion, which began ∼20 ms poststimulus onset, consisted of sustained oscillations. The oscillation frequency matched the frequency of the acoustic tone burst and the amplitude of oscillations was stable over the duration of the response. The onset portion of the V_e signal comprised the first ∼20 ms of the response. It also exhibited oscillations, but with amplitudes that grew (after ∼5 ms latency) and attenuated before settling into the sustained ongoing response.

The transition from onset to sustained response is likely due, in part, to physiological processes such as adaptation in synaptic transmission and neural firing. In addition, the filtering properties of our recording equipment affected the recorded signal in significant ways. Specifically, the TREC headstages were not DC coupled (10 Hz lower cutoff frequency) and thus removed “steady-state” and low-frequency components (<10 Hz) of neurophonic responses. High-pass filtering impedes interpretation of field potential data because the time-average of ongoing portions of V_e recordings has 0 mean (no “baseline”; Herreras, 2016). We make the observation that during a brief portion of the response near signal onset, there is transient evidence of the signal baseline, even after high-pass filtering; Figure 1 shows a demonstration. In some analyses, we use this onset portion of V_e responses to overcome ambiguities associated with interpreting the ongoing portion of responses.

Figure 1.

Demonstration of the action of a 10 Hz high-pass filter. A, A signal composed of a constant baseline and 250 Hz sine wave that have an instantaneous onset at t = 10 ms (signal is for demonstration purposes only, it is not actual data). B, Same signal after processing with a 10 Hz high-pass filter. Filter is implemented by solving, in discrete time, a differential equation to describe an RC circuit: RC(ds/dt − ds̃/dt) = s̃, where s is the original signal and s̃ is the filtered signal. The constant RC is equal to 1/2πf_c, where f_c = 10 Hz is the cutoff frequency (Bartiromo and De Vincenzi, 2016).

In most recording sessions, we obtained a single V_e recording 2000 ms long at multiple spatial locations within the brainstem. There is variability (“noise”) in the recording, but we isolated a meaningful signal from the ongoing portion of the response by averaging with respect to the period of oscillation. Specifically, we associated recording time with a corresponding phase value (relative to some arbitrary starting phase), partitioned each period into small phase “bins” and averaged V_e values that shared similar phase values. For ease of analysis, we used linear interpolation to resolve N = 99 distinct phase values so that all spatiotemporal patterns were the same size, regardless of stimulus frequency. We refer to this process as “cycle-averaging.” We illustrate the process of creating a cycle-averaged spatiotemporal V_e pattern in Figure 2. Onset responses were not stationary over time, so we did not cycle-average onset responses. Instead, to study onset responses, we selected a single cycle of the response that contained the time at which V_e reached its maximum amplitude (computed over all spatial locations).

Figure 2.

Construction of spatiotemporal V_e patterns from neurophonic responses. A, Portion of ongoing V_e response recorded at one location in the brainstem. Vertical gray lines mark cycles with respect to 1200 Hz tone (period T = 0.83 ms). B, Responses are aligned with respect to the period of oscillation and averaged at each phase value to create the cycle-averaged V_e response (black line). C, Cycle-averaged V_e responses constructed at all recording locations are aggregated to create a spatiotemporal V_e pattern. D, Spatiotemporal V_e patterns accurately reconstructed from three Fourier modes (see Materials and Methods).

In two recording sessions, we presented 100 repetitions of a single tone frequency (700 or 1300 Hz, 50 ms duration). In these cases, we could average over the repeated trials to enhance the signal-to-noise ratio in our data. After averaging responses over trials, we created spatiotemporal V_e patterns from the onset and ongoing portions of these responses, as described above.

Approximation of V_e spatiotemporal patterns in the Fourier domain.

To simplify our model-fitting procedure (described below), we found it helpful to represent cycle-averaged V_e responses in the Fourier domain (temporal frequency). Specifically, we computed Fourier coefficients from spatiotemporal V_e patterns by evaluating the discrete Fourier transform (Kido, 2015): The notation V_e^CA (x, k) emphasizes that we used cycle-averaged data in this calculation, not the “raw” V_e recordings (except in analyses of onset responses, for which we used one cycle chosen based on the peak amplitude of V_e). In this equation, x identifies the spatial location (depth in the brainstem) of the recording electrode, k is the discrete time point within a cycle, and N is the period in discrete time units (N = 99 in our analyses). The index j identifies the discrete frequency associated with each Fourier coefficient (f_discrete = jf_tone, where f_discrete is the discrete frequency and f_tone is the stimulus frequency). The Fourier coefficient a₀ represents the 0 frequency “DC” component of the V_e response, the Fourier coefficient a₁ and its complex conjugate a_N₋₁ represent V_e responses at the fundamental frequency (same as the stimulus frequency). Complex conjugate pairs with higher index values, such as a₂ and a_N₋₂, are Fourier coefficients for higher harmonics (multiples of the stimulus frequency).

After computing the discrete Fourier coefficients, we reconstructed cycle-averaged V_e^CA responses using the inverse discrete Fourier transform (Kido, 2015): The 0 frequency coefficient a₀ was only required for analysis of onset responses. The 10 Hz lower cutoff frequency in the recording equipment removes this component from ongoing portions of V_e responses as discussed above (the TREC headstages were not DC-coupled).

Importantly, we found that we could set a_j = 0 for all but a few Fourier modes, and still accurately reconstruct V_e^CA using the above equation. This observation was key to simplifying our model-fitting procedure (described below). Specifically, we found that three pairs of Fourier coefficients {a_j, a_N−j} for j = 1, 2, and 3 sufficed to accurately reconstruct ongoing spatiotemporal V_e^CA patterns. For onset patterns, we used these three pairs and, in addition, retained the constant term a₀. An example V_e^CA pattern and its Fourier-based reconstruction (using 3 pairs of Fourier coefficients) are shown in Figure 2 (bottom row).

We measured error in the Fourier approximation by treating cycle-average patterns as matrices (space-time) of V_e values and using the Frobenius norm to define relative error as follows: where, for an arbitrary matrix M with elements m_ij, the Frobenius norm is

Relative errors of this approximation for varying numbers of Fourier coefficients are shown in Figure 3(top row). We chose to use three Fourier modes (the 3 pairs {a_j, a_N−j} for j = 1, 2, and 3), plus the 0 frequency coefficient for onset responses, in our analyses and modeling because this approximation was satisfactorily accurate across recordings and stimulus frequencies. Figure 3, bottom row, reports relative errors across all recording sessions and stimulus frequencies for the three-mode Fourier reconstruction. The error was greater for responses to lower-frequency tones when compared with higher-frequency tones. Also, onset responses were not as accurately approximated as ongoing responses. This is expected because onset V_e patterns were computed from a single cycle of response whereas ongoing V_e patterns were computed from cycle-averaged data.

Figure 3.

Relative error in Fourier domain reconstruction of V_e patterns. Relative error decreases with the number of Fourier modes used in approximation of V_e patterns for onset responses to contralateral tones (A1), ongoing responses to contralateral tones (B1), onset responses to ipsilateral tones (C1), and ongoing responses to ipsilateral tones (D1). Error bars in top row are mean and standard deviation across five recording sessions and 12 stimulus frequencies (600–1800 Hz, increments of 100 Hz). In all subsequent analyses, we use three Fourier modes (and additionally for onset responses, the 0 frequency Fourier mode). Detailed results for relative error in Fourier domain reconstruction of V_e patterns using three modes (and a 0 frequency mode for onset responses) for onset responses to contralateral tones (A2), ongoing responses to contralateral tones (B2), onset responses to ipsilateral tones (C2), and ongoing responses to ipsilateral tones (D2). Solid line in bottom row is mean over all recording sessions and colored dots are results for each recording session (key at right).

Mathematical model.

Our goal was to model extracellular voltage patterns generated by (simulated) neural activity. Moreover, we sought a model that could be fit to spatiotemporal V_e patterns recorded in experiments. As in our previous studies of the neurophonic (Mc Laughlin et al., 2010; Goldwyn et al., 2014), we assumed that extracellular voltage in the auditory brainstem varied in a direction parallel to the orientation of dendrites in the MSO and was relatively constant in other directions. In addition, we assumed the spatiotemporal pattern of membrane current in a neuron I_m generates V_e responses according to Poisson's equation (electrostatic approximation; Mitzdorf, 1985). Thus, our model for V_e included a single spatial dimension and Poisson's equation reduced to: where r_e is electrical resistivity in the extracellular medium (assumed to be isotropic and homogeneous). The terminals of the model neuron are denoted x_C₁ and x_C₂. We abbreviate the right-hand side of the equation as J(x, t). The function J(x, t) is proportional to membrane current between the terminals of the cable, and 0 beyond the terminals of the cable. We defined boundaries in the extracellular domain to be two points, x_G₁ and x_G₂, that are distant from the terminals of the cable: x_G1 ≪ x_C1 and x_G2 ≫ x_C2. We imposed the boundary conditions V_e(x_G1, t) = V_e(x_G2, t) = 0. In other words x_G₁ and x_G₂ represent the location of electric ground.

The starting point for our model of membrane current I_m was the classical description of a neuron as a cable with a passive (leaky) membrane (Rall, 1977). The membrane potential as a function of position along the neuron x and time t, denoted V_m(x, t), satisfies the linear partial differential equation (cable equation): where c_m is membrane capacitance per unit length, r_m is membrane resistance, r_i is intracellular (axial) resistance per unit length, and E_lk is the leak reversal potential. The value of E_lk has no effect on our results, so we set E_lk = 0 mV for simplicity (or equivalently, we think of V_m in the above equation as deviation from a resting potential. The term I_in represents synaptic input currents and will be discussed in more detail below. Nonlinearities such as spike-generating sodium current and voltage-gated currents that are known to be present in MSO neurons (e.g., low-threshold potassium) are not included in our model. Spikes in MSO neurons are small (Scott et al., 2007) and typically difficult to identify in extracellular recordings, so we did not expect that they contribute appreciably to neurophonic responses. In our previous work, we observed that the contribution of low-threshold potassium current to V_e responses was small relative to synaptic currents (Goldwyn et al., 2014). We would not expect, therefore, that our results would change substantially if this current were included.

Following standard practice we introduced the membrane time constant parameter τ_m = r_mc_m and the space constant parameter λ = and reformulated Equation 6 as The spatial domain of the cable is a finite interval [x_C₁, x_C₂] and we imposed “sealed end” (no flux) boundaries at the end points = 0.

The cable equation expresses a conservation of current relationship. In particular, membrane current in this model is the sum of capacitive, leak, and input currents and can be written (after rearranging terms): Our approach to simulating V_e can now be stated briefly. For given input current I_in and parameter values (τ_m, λ, and others), we computed V_m by solving Equation 7. We then computed I_m using Equation 8. Last, with I_m determined, we specified the right-hand side of Equation 5 [denoted as J(x, t)] and found V_e(x, t) using a Green's function (Tuckwell, 1988): This procedure establishes a straightforward and analytical connection between MSO neuron activity (cable equation for V_m) and the auditory neurophonic (Poisson's equation for V_e). The critical “missing link” in this sequence is the presumed knowledge of the input current I_in. Our experimental measurements were extracellular (we recorded V_e). We did not have access to the synaptic currents driving neural activity in the MSO. The challenge, therefore, was to infer I_in to produce simulated V_e patterns that accurately reproduced in vivo data. We describe our method for determining I_in in the Model fitting procedure section, below.

Equation 9 establishes how the transmembrane current I_m(x, t) of one MSO neuron (modeled as a passive cable) generates an extracellular voltage response V_e(x, t). The auditory neurophonic, like most field potentials, represents the combined activity of many neurons that are “near” the recording electrode. We replaced I_m(x, t) in Equation 5, therefore, with membrane current summed over many neurons: To simplify this population-level description of membrane currents, we considered an idealized group of MSO neurons. Neurons in this local subpopulation were modeled as identical cables, receiving identical inputs I_in(x, t), and oriented in parallel to one another along a common spatial dimension (x-axis). In other words, I_m^(j) ≡ I_m and V_m^(j) followed Equation 7 for all neurons in the subpopulation. Thus, the aggregate membrane current summed over a population of n neurons is as follows: To allow for some variation in the spatial position of each model neuron, we assumed the center of each neuron Δx_j was drawn (independently) from a Gaussian distribution with mean x_c (the center of the MSO) and variance σ².

In the limit of a large population of neurons, this sum could be replaced with the convolution integral There are two additional parameters in this model of aggregate membrane current. The amount of “spatial jitter” in the population of MSO neurons is σ and the number of neurons presumed to contribute to recorded field potentials is n. As we explain below, however, n and the extracellular resistance r_e did not need to be specified when fitting simulated V_e responses to data.

Model fitting procedure.

As discussed above, spatiotemporal V_e patterns are the basic unit of data that we considered in this study. Our goal, then, was to find parameter values and input currents so that solutions of Equation 9 (simulated V_e) matched recordings of the auditory neurophonic. Two simplifications made this model fitting procedure tractable. First, because spatiotemporal V_e patterns were accurately approximated in the Fourier domain, it sufficed to fit a small number of Fourier coefficients. Second, we stipulated that I_in(x, t) was the sum of two point currents localized to two distinct spots on the cable: The symbol δ(x − x_i), where = 1 or 2, is the Dirac delta function. It is 0 everywhere except if x = x_i, has area equal to 1, and thus represents the assumption that synaptic inputs are point sources (localized and not spread out along the cable). Unless otherwise stated, the location of the first input (x₁) is the center of the cable, and the second input targets an off-center position (x₂, 75 μm from the cable center for onset responses and 150 μm from the cable center for ongoing responses).

In this form, I_in(x, t) has a natural interpretation based on known properties of MSO neurons. Excitation evoked by monaural sound stimuli are segregated (ventromedial side of MSO for contralateral sounds, the dorsolateral side of MSO for ipsilateral inputs). Inhibition to MSO neurons is known to primarily target the soma (Kapfer et al., 2002) and excitation is known to primarily target dendrites (Couchman et al., 2012). The input currents can, therefore, be tentatively identified as soma-targeting inhibition (I₁) and dendrite-targeting excitation (I₂). Figure 4 illustrates the arrangement of synaptic inputs to MSO neurons and the idealized configuration of input currents in the model. Under this interpretation, fitting the model can be viewed as a test of a “working hypothesis” that the combined contributions of inhibition and excitation to MSO neurons are the dominant generators of the spatiotemporal V_e patterns observed in vivo.

Figure 4.

Illustration of synaptic inputs. A, MSO neurons receive dendrite-targeting excitation that arrives via spherical bushy cells (SBC) and soma-targeting inhibition that arrives via the MNTB and LNTB. Excitatory inputs are segregated onto opposite dendritic arbors based on the side of stimulation. The side of stimulation also determines the source of inhibition (MNTB for contralateral stimulation; LNTB for ipsilateral stimulation). B, Idealized model of an MSO neuron. The neuron model consists of a leaky cable [linear V_m dynamics; see Eq. 7 and two input currents (I₁ and I₂)]. One input arrives at the center of the cable and the other arrives away from the center. The side of the second input I₂ changes with side of sound stimulus. A plausible interpretation is that I₁ represents soma-targeting inhibition and I₂ represents dendrite-targeting excitation.

Fitting I₁(t) and I₂(t) was tractable when we worked in the Fourier (temporal frequency) domain. Because V_e responses were accurately approximated by three frequency components (and a constant term for onset responses), we could rewrite Equation 13 in terms of Fourier coefficients and sinusoidal temporal dynamics as follows: where f is the frequency of the sound stimulus and the Fourier coefficients α_j and β_j are complex-valued functions. The zero frequency “DC” components α₀ and β₀ were non-zero only when modeling onset responses, for reasons discussed above (10 Hz lower frequency cutoff in V_e recordings). With I_in expressed as this sum of complex exponentials, we could solve the preceding equations to determine V_e. This formulation admits complex-valued solutions, with coefficients that include both magnitude and phase information to completely describe the waveform of V_e (which is known to be real-valued). To make comparisons to data, therefore, we extracted the real-valued portion of the calculated V_e.

This representation of I_in substantially simplified the problem of fitting simulated V_e responses to data because I_in(x, t) is completely determined by a small number of Fourier coefficients (α_j and β_j). Additionally, all equations in the model are linear in time. As a result, orthogonality of the complex exponential function ensures that Fourier coefficient pairs can be fit “individually”, i.e., values of α_j and β_j do not depend on values of α_k and β_k for j ≠ k.

When possible, we chose the values of other parameters in the model to reflect the known physiological and anatomical properties of MSO neurons. The physical length of the cable model neuron is 430 μm (Stotler, 1953). We imagine this cable represents two dendrites of length 200 μm each and a soma of length 30 μm. The space constant of the cable is therefore set to be λ = 200 μm, or approximately the length of one dendrite (Mathews et al., 2010). The membrane time constant is τ = 0.3 ms and reflects the exceptionally fast dynamics of MSO neurons (Golding and Oertel, 2012). These parameter values were fixed for all model fits (regardless of MSO, tone frequency, and contralateral or ipsilateral side of tone presentation).

We selected the remaining parameters by hand to obtain satisfactory fits to the data. They included the location of the cable center (Eq. 14, x₁), the “spatial jitter” of the population (Eq. 12, σ), and the locations of electric ground (Eq. 9, x_G₁ and x_G₂). We selected different parameter values for each recording session (i.e., each MSO), but the parameter values did not change with tone frequency. We found it necessary to use different values for onset and ongoing portions of responses, and for contralaterally- and ipsilaterally-presented tones. These parameter values are reported in Table 1. We consider possible reasons for these differences in parameter values in the Discussion.

When reporting results, we often speak of input currents obtained using the model. To be precise, our model fitting procedure estimates the quantity nI_in(x, t). The parameters for membrane resistance (Eq. 6, r_m), electrical resistivity of the extracellular domain (Eq. 5, r_e) and the number of neurons in the neuron population (Eq. 12, n) enter as multiplicative factors. All equations in the model are linear, however, so these parameters have no effect on solutions other than rescaling the amplitude of I_in. We find it convenient, therefore, to refer to this quantity as “input current” and trust this will not cause confusion. We typically report normalized values of this quantity, so there is no need to specify the values of these scaling factors.

With all parameters selected, we used a global optimization algorithm in MATLAB (lsqnonlin) to find Fourier coefficients α_j and β_j. We sought to minimize the error between simulated and recorded spatial-temporal patterns of V_e. Specifically, for each frequency component used in the Fourier representation of V_e, we solved the minimization problems (j = 0, 1, 2, 3 for onset data, j = 1, 2, 3 for ongoing data): where a_j^(data) are Fourier coefficients defined as in Equation 1 from spatiotemporal V_e data and a_j^(model) are Fourier coefficients extracted from simulated V_e responses. The index i on the spatial variable indicates the finite set of locations (electrode positions) at which V_e recordings were obtained. Once α_j and β_j were determined, we reconstructed the input currents I₁(t) and I₂(t) in the time domain by extracting the real part of the complex-valued functions described by Equation 14.

Although we fit models by minimizing error in the Fourier domain, we use the relative error measure introduced above (Equation 3) when reporting the quality of model fits. In particular, we treated cycle-averaged neurophonic responses as matrices V_e^CA(data) and computed the relative error of the simulated responses V_e^CA(model) patterns as follows: We chose a soma-targeting position for I₁ and dendrite-targeting input for I₂ based on known patterns of synaptic inputs to MSO. We do not claim to be identifying unique and globally optimal solutions to this minimization problem. As we discuss below, similar fits can be obtained for different locations of input currents.

All computations were performed on a laptop computer using the MATLAB scientific computing software package (RRID:SCR_001622).

Results

“Dipole” and “non-dipole” features in V_e responses

The standard conceptual model for neurophonic responses, what we have termed the dipole theory, posits that membrane currents in MSO neurons are the dominant generators of auditory brainstem V_e. This account caricatures MSO neurons as current dipoles because dendrite-targeting excitation generates inward (“sink”) current and the return of this current to the extracellular domain through the soma and opposite dendrite generates outward (“source”) current.

The V_e responses shown in Figure 5 exhibit features that we classify as dipole-like, i.e., consistent with the dipole theory. These data (which are, in fact, the average of 100 repeated trials) were recorded in response to a 1300 Hz pure tone presented to the ear contralateral to the side of the brain in which the data were collected. The three V_e time courses in Figure 5A were recorded at locations separated by increments of 300 μm. From histology we estimated that the green curve was recorded near the MSO center and the blue and red curves were recorded medial and lateral to the MSO center, respectively (negative depth indicates more ventromedial location, positive depth indicates more dorsolateral). These V_e time courses oscillate at 1300 Hz, which matches the frequency of the stimulating tone. Such sustained sound-evoked oscillations are the signature of the auditory neurophonic. The essential dipole-like feature of these responses is the anti-phase relationship between the temporal oscillations in V_e signals recorded medial to the MSO center (blue curve, negative depths) and lateral to the MSO (red curve, positive depths). This anti-phase relationship is present in both onset and ongoing responses (Fig. 5A1,A2).

Figure 5.

Dipole-like neurophonic response to 1300 Hz contralateral tone. A, Extracellular voltage (V_e) recordings at three locations. Portions of onset responses are shown in A1 and portions of ongoing responses are shown in A2. On opposite sides of the MSO, V_e oscillated roughly out of phase. Position relative to MSO estimated from histology (negative depth indicates more ventromedial location; positive depth indicates more dorsolateral). B, Spatiotemporal V_e patterns constructed from data. Onset pattern (B1) was created from one cycle, indicated by the shaded region in A1. Ongoing pattern (B2) was created from cycle-averaged data. Horizontal black lines in B1 and B2 mark the center of the MSO, estimated from histology. C, Simulated V_e patterns obtained from fitting model to onset responses (C1) and ongoing responses (C2). Black bars indicate location of model MSO neuron. Color scale is same as corresponding responses in B. D, Input currents obtained in fitting model to neurophonic responses. For onset responses (D1): soma-targeting I₁ was a source (positive-valued, outward current) and dendrite-targeting I₂ was a sink (negative-valued, inward current). For ongoing responses (D2): I₁ and I₂ have zero mean and thus cannot be identified as sinks or sources.

We visualize the spatiotemporal dynamics of V_e responses using two-dimensional color maps. These “patterns” depict one cycle of the onset V_e response (Fig. 5B1) and the cycle-averaged ongoing response (Fig. 5B2). One can visually identify these V_e patterns as dipole-like by observing that temporal oscillations in V_e at depths lateral to the MSO (larger depths on the y-axis) appear to be anti-phase to oscillations at more medial positions (smaller depth values). Indeed, if one were to draw a cross-section at a fixed moment in time, one would observe a sharp transition (and a null in V_e) around a depth of 1000 μm. This position corresponds approximately (based on histological observations) to the location of the MSO center, which we mark with a black line in Figure 5B1 and B2.

We will discuss modeling results in more detail below, but here we provide initial evidence that we accurately reproduced both onset (Fig. 5C1) and ongoing (Fig. 5C2) spatiotemporal patterns of V_e responses to monaural tones using our mathematical model. The currents that drove activity in the model neuron (and were obtained from fitting the model to V_e data) are shown in Figure 5D1 and D2. For this dipole-like pattern, the amplitude of the dendrite-targeting current I₂ was larger than the amplitude of the soma-targeting current I₁.

A primary motivation for the present work was that we observed diverse and frequency-dependent V_e responses to monaural tones that are not always “dipole-like.” In particular, sustained responses to lower-frequency tones tended to exhibit monopole-like features as illustrated with an example in Figure 6. These are V_e responses to a 700 Hz contralaterally presented tone (average of responses in 100 trials). The time courses in Figure 6, A1 and A2, are from three locations in and around the MSO. The prominent temporal oscillation at the tone frequency (700 Hz in this case) is again apparent.

Figure 6.

Monopole-like features in neurophonic response to 700 Hz contralateral tone. Format is same as Figure 5, but data were obtained in a different recording session and using a lower frequency tone stimulus. A, Extracellular voltage (V_e) recordings at three locations. Portions of onset responses are shown in A1 and portions of ongoing responses are shown in A2. On opposite sides of the MSO, V_e oscillated roughly in an in-phase relation. B, spatiotemporal V_e patterns constructed from data. Onset pattern (B1) was created from one cycle, indicated by the shaded region in A1. Ongoing pattern (B2) was created from cycle-averaged data. Horizontal black lines in B1 and B2 mark the center of the MSO, estimated from histology. C, Simulated V_e patterns obtained from fitting model to onset responses (C1) and ongoing responses (C2). Black bars indicated location of model MSO neuron. Color scale is same as corresponding responses in B. D, Input currents obtained in fitting model to neurophonic responses. For onset responses (D1): soma-targeting I₁ was a source (positive-valued, outward current) and dendrite-targeting I₂ was a sink (negative-valued, inward current).

The spatial distribution of these V_e responses differed from the previous dipole-like example. Extracellular voltages obtained at positions medial (blue) and lateral (red) to the MSO did not exhibit an anti-phase relationship. Rather, V_e responses at all three locations oscillated nearly in phase with one another. The coherent oscillation was present in both the onset and ongoing portions of the response. Consequently, the dominant feature of the spatiotemporal V_e pattern of these responses was a temporal oscillation that was coherent across recording depths (see especially the ongoing pattern in Fig. 6B2). The largest amplitude of the spatially coherent temporal oscillation was at a recording depth near 1100 μm. This position was near the center of the MSO, as estimated from histology.

In Figure 6C, we show that simulated V_e patterns accurately reproduced the V_e patterns observed in onset and ongoing portions of this dataset. Reproducing these monopole-like V_e patterns required strong soma-targeting current I₁ (Fig. 6D). Soma-targeting current produced spatially-coherent V_e responses in the model with input current at the center of the cable offset by return currents distributed symmetrically across the cable. The dendrite-targeting current I₂ was also present and, as we will show, the combination of a soma-targeting input (that creates monopole-like V_e patterns) and a dendrite-targeting input (that creates dipole-like V_e patterns) sufficed to accurately simulate non-dipole-like neurophonic responses.

The preceding examples illustrated our parsimonious and biophysically-based update to the dipole theory. In the remainder of this work, we will demonstrate that the model flexibly reproduced diverse and frequency-dependent neurophonic responses. In addition, we will argue that it is plausible to view dendrite-targeting current in the model as a signature of excitatory inputs to MSO neurons and soma-targeting current as a signature of inhibitory inputs to MSO.

Simulated V_e patterns reproduce diverse and frequency-dependent neurophonic responses

We identified five MSOs (in 4 cats) in which neurophonic responses exhibited qualitatively similar spatiotemporal V_e patterns. These patterns varied with stimulus frequency and with side of sound presentation. Contralateral responses at three frequencies are in Figure 7 and ipsilateral responses are in Figure 8. In all cases we found that the model could be fit to data so that simulated V_e patterns accurately reproduced the diverse and frequency-dependent neurophonic responses observed in vivo. Spatiotemporal V_e patterns obtained from data and simulations are interleaved in alternating columns in Figures 7 and 8. All V_e patterns shown in these figures are ongoing portions of responses.

Figure 7.

Spatiotemporal V_e patterns created from ongoing portions of responses to contralateral monaural tones. Experimental data and corresponding simulation results are interleaved column-wise. Responses are shown for three frequencies (600, 1200, and 1800 Hz; columns) from five recording sessions (rows). Temporal extent of each V_e pattern is one cycle of the response. Color map is rescaled independently for each MSO, as indicated in key at right. Responses to 600 Hz exhibit a monopole-like spatially-coherent oscillation. Responses to 1200 Hz exhibit a dipole-like pattern: V_e oscillations at different depths are anti-phasic. Responses to 1800 Hz exhibit a wave-like non-dipole pattern.

Figure 8.

Spatiotemporal V_e patterns created from responses to ongoing portions of ipsilateral monaural tones. Format is same as Figure 7: responses are shown for three frequencies (600, 1200, and 1800 Hz; columns) from five recording sessions (rows) with experimental data and simulation results interleaved column-wise. Temporal extent of each V_e pattern is one cycle of the response. Color map is rescaled independently for each MSO, as indicated in key at right. Some features of these responses are similar to responses to contralateral tones; responses to 600 Hz tones tend to appear monopole-like and responses to 1200 Hz tones tend to appear dipole-like. Some responses appear more complicated and idiosyncratic; for example, see MSO #1 at 600 Hz.

The dominant feature of low-frequency responses (600 Hz, left columns) was a monopole-like temporal oscillation that is coherent along the spatial (vertical) dimension. Responses to 1200 Hz tones (central columns) exhibit dipole-like patterns, apparent in the (nearly) anti-phase relation between temporal oscillations at small recordings depths relative to larger depths. Extracellular voltage responses to high-frequency tones (1800 Hz, right columns), exhibited patterns that were dipole-like in some ways, but also had the appearance of “traveling waves”.

Model fits to contralateral responses were accurate across all recording sessions and stimulus frequencies (Fig. 9A,B). Ipsilateral responses exhibited, on occasion, more complex patterns (see MSO 1 response to 600 Hz tone in Fig. 8). As a result, model fits to data were more accurate for contralateral responses than ipsilateral responses (Fig. 9C,D). The error we are reporting in Figure 9 measured differences between model fits and cycle-averaged data, but model parameters were determined from comparisons between simulations and the Fourier domain approximation of the data (see Materials and Methods). The accuracy of the model was limited, therefore, by the accuracy of the Fourier approximation. In other words, the relative errors reported in Figure 3 are a lower bound for relative error in the model.

Figure 9.

Relative error in fitting simulated spatiotemporal V_e patterns to data. A–D, Colored dots indicate results for each recording session (color code given in key at right) and solid black line is the error averaged over the five MSOs. The model used in these cases included two input currents (one targeting center of neuron; one targeting off-center). Dotted line shows mean relative error when the model included a single input targeting the center of the model neuron (soma only configuration). Dashed line shows mean error when a single input targeting a position away from the center of the model neuron was used (dendrite only configuration).

We distinguished dipole-like from “non-dipole-like” V_e features primarily on the basis of visually inspecting spatiotemporal patterns like those shown in Figures 7 and 8. To make this classification more objective, we fit versions of the model that included input current at only a single location. We reasoned that if a spatiotemporal V_e pattern could be accurately approximated by a model that included an input current targeting one off-center (“dendrite”) location, then the V_e response could reasonably be labeled dipole-like. Alternatively, if a spatiotemporal V_e pattern could be accurately approximated by a model that included one input current targeting the center (“soma”) of the model neuron, then the V_e response could be termed monopole-like.

By this measure, ongoing V_e responses transitioned from monopole-like for low-frequency tones to dipole-like for higher-frequency tones (Fig. 9B,D). Specifically, the “soma input only” model (dotted line) yielded more accurate fits to V_e data at low frequencies relative to the “dendrite input only model” (dashed line). This relationship reversed at frequencies >∼1000 Hz. In all cases, the “standard” model with two inputs yielded more accurate fits to data than models with one input only. This is not surprising, but it emphasizes our point that the dipole theory does not provide a comprehensive and quantitative account of neurophonic responses. Instead, these results support the “working hypothesis” that V_e patterns are shaped by the combination of dendrite-targeting excitation and soma-targeting inhibition.

We also used the model to accurately reproduce onset portions of neurophonic response to monaural tones across the same range of frequencies. Onset V_e patterns are not displayed, but see Figure 9, A and C, for relative errors and Figures 5 and 6 (left columns) for examples. Onset responses tended to require two inputs to be accurately modeled (errors in fits using dendrite-only or soma-only models were much higher than errors in fits using both inputs for all frequencies studied).

In sum, model fitting of V_e responses required a combination of dendrite-targeting and soma-targeting currents to accurately reproduce the diverse and frequency-dependent neurophonic responses observed in vivo. The contribution of soma-targeting current was most apparent for responses to lower-frequency tones, which appeared monopole-like. The contribution of dendrite-targeting current was most apparent for ongoing responses to higher-frequency tones, which appeared dipole-like.

Onset responses identify I₁ and I₂ as source and sink current, respectively

As illustrated by the preceding results, we can use our idealized modeling approach to accurately reproduce extracellular voltage recordings. The key step in the model-fitting procedure was the determination of two input currents: I₁(t) and I₂(t) in Equation 13. By construction of the model, there was a natural interpretation of these input currents: I₁ as soma-targeting inhibition and I₂ as dendrite-targeting excitation to MSO neurons. To support this interpretation, we examined these inputs currents further.

To reasonably identify I₁ as inhibitory current and I₂ as excitatory current, we expected that I₁ should represent an outward current (source) and I₂ should represent an inward current (sink) in the model. In other words, following standard convention, we expected I₁ to take positive values and I₂ to take negative values. The input currents obtained from model fits of ongoing responses had no baseline (DC) component (as discussed in Materials and Methods). They oscillated with a time-average of zero and, consequently, could not be classified as sink or source current. We focused, therefore, on I₁ and I₂ obtained from fits to onset portions of responses. There was a baseline component in these currents that could possibly offset I₁ and I₂ to be positive- or negative-valued.

The ranges of I₁ and I₂ values obtained in all model fits to onset responses are displayed in Figure 10. For the parameter sets we used, input currents fit to onset responses were segregated (for the most part) into negative-valued I₁ and positive-valued I₂. In other words, the soma-targeting current (I₁) was a source (outward current) and the dendrite-targeting input (I₂) was a sink (inward current). Because inhibition is primarily soma-targeting in MSO and excitation is primarily dendrite-targeting, we could plausibly conclude that I₁ reflects inhibitory inputs to MSO and I₂ reflects excitatory inputs to MSO.

Figure 10.

Range of I₁ and I₂ values obtained from model fits to spatiotemporal V_e patterns. Results for contralateral (left column) and ipsilateral (right column) onset responses in five recording sessions (rows). Scale of I₁ and I₂ is undetermined (see Materials and Methods), so we simply divide y-axis into positive values (source current by convention) and negative values (sink current). Most I₂ values are negative (sink, inward current) and most I₁ values are positive (source, outward current).

Characteristics of putative inhibition and putative excitation

Putative inhibition attenuates with tone frequency and putative excitation exhibits a “best frequency”

We have argued that the inputs to the model can plausibly be viewed as inhibition (I₁, targets center of model neuron) and excitation (I₂, targets away from center of model neuron). To assess this interpretation, we measured the amplitude of these putative currents obtained from model fits to each V_e response (5 MSOs, 12 frequencies in increments of 100 Hz from 600 to 1800 Hz, ongoing responses). We defined amplitude as peak-to-trough difference in the waveforms I₁(t) and I₂(t). To facilitate comparisons across different MSOs, we normalized the amplitude measure in each MSO and for each side of sound input (contralateral or ipsilateral ear) to the maximum amplitude value obtained for I₁ over all frequencies.

For the parameter sets we used to fit V_e, we found consistent trends across recordings sessions. The amplitude of I₁ decreased as the frequency of the sound stimulus increased, as shown in Figure 11A and B. There were a few exceptions at the very lowest frequencies for contralateral responses, but overall the amplitude profile of I₁ was “low-pass”. In contrast, the amplitude profiles of I₂ were distinctly non-monotonic. We defined a “best frequency” for each MSO and each side of stimulation as the maximum value of the amplitude profiles. Note this best frequency refers only to the response at which amplitude of putative excitatory inputs to MSO neurons was largest, it is not a measure of firing rate (spiking) in MSO. The median best frequency for contralateral and ipsilateral responses was 1300 Hz.

Figure 11.

Peak-to-trough amplitude of input currents obtained from fitting model to spatiotemporal V_e patterns. Amplitudes are normalized for each MSO to the maximum amplitude of I₁ for all frequencies. I₁ amplitudes for contralateral responses (A) and ipsilateral responses (B). I₂ amplitudes for contralateral responses (C) and ipsilateral responses (D). Color code for recording sessions is given in B. Solid line is the mean over recording sessions.

We cannot directly relate the amplitude-frequency profiles in Figure 11A–D to known synaptic currents because we only measured extracellular voltage in these experiments. Nonetheless, these amplitude-frequency profiles are consistent with known properties of synaptic inputs to MSO. For instance, Couchman et al. (2010) found the average decay time constant for evoked EPSCs in vitro to be 270 μs and for IPSCs to be 1.76 ms (gerbil MSO). Some have suggested that IPSCs could have fast (submillisecond time scale) kinetics (Brand et al., 2002), but in vitro evidence consistently identify postsynaptic inhibitory currents and potentials with millisecond-scale dynamics; (Grothe and Sanes, 1994; Jercog et al., 2010; Fischl et al., 2012; Roberts et al., 2013). Direct measurements of synaptic currents in vivo have not been reported, but EPSPs in gerbil are fast relative to IPSPs (Franken et al., 2015). Because dendrite-targeting excitation in MSO is fast and temporally precise (Joris et al., 1994a), it is plausible that EPSCs can contribute to the neurophonic on a cycle-by-cycle basis in response to relatively high tone frequencies. Soma-targeting inhibition is slower, thus we expect the contribution of inhibition to the neurophonic to have low-pass properties. IPSPs temporally summate at frequencies >200–300 Hz (Grothe and Sanes, 1994; Roberts et al., 2013), although some cycle-by-cycle oscillatory component is visible at 800 Hz in the study by Roberts et al. (2013) (in vitro recordings in gerbil).

Many properties of synaptic transmission shape the dependence of current amplitude on stimulus frequency (stochastic vesicle release, temporal jitter, adaptation, etc.). The tonotopic organization of MSO and its inputs likely contributes to the non-monotonic profile of I₂. The best frequency of 1300 Hz may reflect the frequency-tuning of excitatory inputs to MSO neurons in the region of the brainstem from which we recorded. In additional recordings (data not shown), we sampled neurophonic responses from several locations in the dorsoventral plane of the MSO and found that the frequency at which V_e amplitudes were maximal varied in a manner consistent with the presumed tonotopic axis of MSO.

Putative inhibition precedes putative excitation

Last, we examined the temporal dynamics of I₁ and I₂ waveforms for ongoing responses. Recent in vitro work (in gerbils) has shown that inhibitory inputs to MSO neurons can precede excitatory events by several hundred microseconds (Roberts et al., 2013), as proposed by Brand et al. (2002). We were interested to know, therefore, if a similar temporal ordering of putative inhibition (I₁) and putative excitation (I₂) emerged from fitting spatiotemporal neurophonic responses.

In response to low-frequency tones, we observed that I₂ time courses appeared to include brief events reminiscent of EPSCs. The timing of these events coincided with the appearance of negative-going regions of V_e in recordings and simulations. The correspondence between negative-going I₂ and V_e is expected because we associate a transient dip in I₂ with a dendrite-targeting current sink. Figure 12A–C provides an example of a response to a low-frequency tone in which this temporal ordering of excitation and inhibition is apparent. Arrows in Figure 12, A and B, mark the location of the brief negative-going events in the V_e patterns.

Figure 12.

Peak of putative inhibition precedes peak of putative excitation. A, Spatiotemporal V_e pattern for responses in MSO 2 to ongoing portion of contralateral 600 Hz tone. B, Simulated V_e pattern obtained from fitting model to data in A. Arrows in A and B mark faint region of negative-going V_e that we identify as a possible “sink event.” C, Model inputs that create the simulated V_e pattern in B. Triangles mark maximum of I₁ and minimum of I₂. Peak of I₁ (putative inhibition) precedes peak of I₂ (putative excitation) by Δt = 135 μs in this example. D, E, Peak of putative inhibition precedes peak excitation for model fits to ongoing responses for all tone frequencies in four of five recording sessions (exception is MSO #4, which shows more irregular timing relationships, data not shown).

Because inward currents are negative-valued by convention, the minimum value of I₂ represents the maximum value of the putative excitatory current. We refer to this as “peak excitation”. Similarly, we defined “peak inhibition” as the maximum value of I₁ (outward currents are positive, by convention). We then measured the difference Δt between the times at which peak inhibition and peak excitation occurred; using the convention that positive Δt indicates that peak inhibition precedes peak excitation in the model. In the example shown in Figure 12C, peak inhibition preceded peak excitation by Δt = 135 μs.

We found the temporal ordering of putative inhibition and putative excitation to be consistent across tone frequency and side of sound stimulus for four of the five MSOs we analyzed. Results for tones presented to the contralateral ear are shown in Figure 12D and results for tones presented to the ipsilateral ear are shown in Figure 12E. MSO number 4 is excluded from this figure because timing results were not consistent across frequency.

Although Roberts et al. (2013) used a different experimental preparation and methodology, these findings concur (circumstantially) with their work and provide evidence that inhibition precedes excitation in MSO neurons. Roberts et al. (2013) used a thick slice preparation to stimulate afferent inputs to gerbil MSO neurons and measured onset times of synaptic events, which were defined in their study as the 20% rise time of PSPs. They observed that the onset of inhibitory PSPs preceded the onset of excitatory PSPs by 380 μs, on average, for contralateral stimuli and 320 μs for ipsilateral stimuli.

Consideration of alternate input configurations

We have found that a model with two input locations accurately replicated diverse and frequency-dependent neurophonic responses observed in vivo in response to monaural pure tones. The model is minimal in the sense that models with a single input could not accurately fit recorded V_e patterns (recall the high relative errors for the “soma only” and “dendrite only” models in Fig. 9). It is also biophysically based in the sense that we can plausibly associate the two inputs with dendrite-targeting excitation and soma-targeting inhibition in MSO neurons.

We have not, however, ruled out other possible input configurations. In addition to the “standard” model (dendrite-targeting and soma-targeting inputs), we constructed two alternative models. In the “same dendrite” input configuration, two inputs targeted the same dendrite (75 and 150 μm from the center of the model neuron). In the “opposite dendrite” input configuration, two inputs targeted positions on the model neuron ±75 μm relative to its center for onset responses (±150 μm for ongoing responses). Schematics of these model configurations are in the left column of Figure 13.

Figure 13.

Consideration of alternative input configurations. Spatiotemporal V_e patterns from data (A) and models (B–D) are constructed from onset responses to contralateral 700 Hz tone. Input configurations are illustrated in schematics at left. The standard model (B) includes a soma-targeting input (I₁) and a dendrite-targeting input (I₂). The alternative models included two inputs that targeted the same dendrite (C) or two inputs that targeted opposite dendrites (D). Spatiotemporal V_e patterns are accurately reproduced by all models (center column). Input currents obtained from model fits are shown in right column. Currents with positive values are identified as sources (outward current) and those with negative values currents are sinks (inward current).

We found that these models replicated V_e responses with the same accuracy as our standard model. The similarity of modeling fitting results for all three input configurations highlights the difficulty of unambiguously inferring synaptic currents and neural dynamics from extracellular voltage data. It also, in our view, reinforces the importance of using available anatomical and physiological knowledge to construct and interpret a model of V_e responses.

Consider the onset response to 700 Hz contralateral tone first depicted in Figure 6. The spatiotemporal V_e pattern is reproduced in central column Figure 13 along with simulated V_e patterns produced by three different input configurations. The simulated V_e patterns are nearly identical.

The input currents obtained from fitting each model are shown in the right column of Figure 13. We labeled positive values as source current and negative values as sink currents, following standard convention. Sink current is an inward current that depolarizes the model neuron so it is natural, as we suggested previously, to associate sinks in the model with synaptic excitation to MSO neurons. Similarly, we associate sources (outward currents) with synaptic inhibition. Following this interpretation, we characterized the same dendrite model as having excitation and inhibition targeting the same dendrite (Fig. 13C). We characterized the opposite dendrite model as having inhibitory inputs to both dendrites with similar time courses (Fig. 13D). Although we cannot state with certainty that these configurations are unrealistic, we view the standard model (dendrite-targeting excitation and soma-targeting inhibition) as most consistent with known properties of MSO neurons and their synaptic inputs.

Discussion

A minimal model of the MSO reproduces diverse and frequency-dependent neurophonic responses

Acoustic stimulation evokes prominent and sustained V_e responses in the auditory brainstem (Tsuchitani and Boudreau, 1964; Mc Laughlin et al., 2010). These field potentials, known as the auditory neurophonic, have characteristic spatial profiles that can resemble patterns of V_e created by current dipoles (Galambos et al., 1959; Biedenbach and Freeman, 1964; Tsuchitani and Boudreau, 1964; Guinan et al., 1972b; Mc Laughlin et al., 2010; Goldwyn et al., 2014).

We constructed an idealized, but physiologically plausible, model of the MSO to simulate these V_e responses. There were two input currents that drove activity in the model. One input targeted the center of the model neuron and the other targeted an off-center location. Based on known morphology of MSO neurons and their synaptic inputs, we associated these inputs with soma-targeting inhibition and dendrite-targeting excitation, respectively. We determined time courses of these putative synaptic currents to accurately reproduce onset and ongoing portions of neurophonic responses to monaural pure tones presented over a range of frequencies. This is an advance over our previous modeling work which, for the most part, made qualitative comparisons between simulated and recorded V_e responses and focused on dipole-like responses to tones with frequencies ≥1000 Hz (Goldwyn et al., 2014).

A testable working hypothesis: inhibition explains non-dipole-like neurophonic responses

Due to the dipole-like spatial patterning of V_e responses and bipolar morphology of MSO neurons, early investigators hypothesized that dendrite-targeting excitation is the dominant generator of the neurophonic (Galambos et al., 1959; Biedenbach and Freeman, 1964). We drew attention here to the fact that auditory neurophonic responses can exhibit “non-dipole” features, particularly in response to low-frequency tones (<1000 Hz, approximately; Fig. 6). Simulations using dendrite-targeting excitation alone did not produce non-dipole V_e patterns, so we adopted the working hypothesis (informed by the physiology and anatomy of MSO neurons) that inhibition is soma-targeting (Kapfer et al., 2002) and that inhibition combines with dendrite-targeting excitation to shape neurophonic responses.

When we fit this parsimonious model to field potential data, we obtained time courses of soma-targeting and dendrite-targeting currents. Several observations supported our interpretation of these inputs as putative inhibition and excitation, respectively. First, for onset responses, dendrite-targeting inputs were negative (Fig. 10). These are sinks in the model (i.e., inward currents), and can plausibly be associated with excitation. Soma-targeting inputs, in contrast, were positive in most cases. These are sources (i.e., outward current), and can be associated with inhibition. Second, putative inhibition tended to attenuate with increasing tone frequency, whereas the putative excitatory current exhibited a maximum amplitude in response to ∼1300 Hz tones (Fig. 11). These results accord with our expectations because inhibitory currents undergo temporal summation in response to high-frequency tones due to their relatively slow (millisecond-scale) kinetics (Grothe and Sanes, 1994; Franken et al., 2015), whereas the amplitude profile of excitation may be shaped by the tonotopic organization of MSO and its inputs. Indeed, in unpublished recordings (to be the subject of a future report) in which we oriented the five electrode array in a plane that sampled different dorsoventral levels, we saw evidence of tonotopic tuning in neurophonic responses. Third, we found that inhibition preceded excitation (Fig. 12). A similar temporal ordering of MSO inhibition and excitation has been observed in vitro (Grothe, 1994; Roberts et al., 2013). Importantly, our observation that inhibition may underlie monopole-like features in the neurophonic could be tested in future experiments. We would expect that blocking inhibition (for instance, using pharmacological methods as done by Brand et al., 2002; Pecka et al., 2008; Franken et al., 2015) would transform V_e responses to appear more dipole-like.

Evaluation of modeling approach

Our plausible and quantitative model demonstrates how excitatory and IPSCs, as distributed across bipolar MSO neurons, can generate the neurophonic. That being said, we have not ruled out alternative hypotheses. For example, we know that different subpopulations of MSO neurons are activated depending on tone frequency (tonotopy of the MSO). It is possible, therefore, that frequency-dependent changes in V_e responses reflect the spatial distribution of MSO activity and/or the particular ways our electrodes sampled these activity patterns throughout the brainstem (i.e., our choice of insertion angle). Another possibility is that nearby brain regions contribute to auditory neurophonic responses. The lateral superior olive is one candidate (Biedenbach and Freeman, 1964; Clark and Dunlop, 1968). Others have suggested that the dipole-like spatial profile of the neurophonic reflects delay lines as postulated by the Jeffress model of sound localization (Jeffress, 1948), see Bojanowski et al. (1989) for cat, and related work in the auditory brainstem of birds (Sullivan and Konishi, 1986; Köppl and Carr, 2008; Carr et al., 2015). In preliminary analyses of neurophonic recordings at different rostrocaudal positions in the MSO, we have not observed evidence of delay lines in V_e responses along that axis.

There are also cellular and biophysical details we excluded from our model, including morphological complexity of dendrites and axons, heterogeneity of MSO response dynamics (Baumann et al., 2013; Remme et al., 2014), and voltage-gated currents, such as low-threshold potassium current (Svirskis et al., 2003) and spike-generating sodium current (Scott et al., 2010). These details shape V_e responses (Reimann et al., 2013; Ness et al., 2016 for contribution of voltage-gated currents to field potentials) and they underlie more accurate models of MSO dynamics, but we have not observed that they change the qualitative aspects of neurophonic responses that we explored in this study. To support this statement, we performed simulations of a more biophysically detailed model that includes voltage-gated low threshold potassium current and some morphological structure (soma region has larger diameter than dendrite region). Details can be found in the studies by Mathews et al. (2010) and Goldwyn et al. (2014). Simulations of this model, with a sinusoidal excitatory input current targeting one dendrite, are shown in Figure 14, left column. The V_e response in Figure 14B1 exhibits a dipole-like spatial pattern. The non-synaptic membrane currents evoked in response to the dendrite-targeting input are composed primarily of source currents near the site of input and in the soma region (Fig. 14C1). We repeated this simulation, but using a homogeneous and passive version of the model. We fixed the conductance variable associated with the low-threshold potassium current to its value for V_m equal to the resting membrane potential, and we reduced the diameter of the soma so that it matched the diameter of the dendrite regions. With these changes, we converted the biophysical model to a passive cable model similar to what we considered in this study. Nonsynaptic membrane currents are distributed differently for this model compared with the biophysical model (Fig. 14C2). In particular, the homogeneous model lacks the prominent source currents in the soma region of the cell. Nevertheless, simulated V_e responses are qualitatively similar and exhibit the characteristic dipole-like pattern Figure 14B2. We suggest the following interpretation of our work: the time courses of the putative synaptic currents that we obtained with the passive cable model may not be quantitatively-precise reflections of excitatory and inhibitory currents in MSO neurons, but they do reveal essential properties of the spatial arrangement and relative timing of sinks and sources that generate the neurophonic.

Figure 14.

Comparison of simulated responses to dendrite-targeting excitation for models with and without biophysical details. Left column, Responses to a model with some biophysical detail (soma region larger diameter than dendrite region, and includes voltage-gated low threshold potassium current throughout the cell). Right column, Responses to model simplified to behave as a passive cable (soma region same diameter as dendrites, and conductance of potassium current is not dynamic). Responses are to a 1000 Hz excitatory sinusoidal input targeting one dendrite (“lower” dendrite in the orientation of these figures). A1, A2, Membrane potential response in time (x-axis) and along the spatial extent of the neuron model (y-axis). B1, B2, Extracellular voltage response. Notice vertical scale (spatial domain) is extended beyond the ends of the model neuron. Location of the model neuron is shown by black vertical bar. C1, C2, Nonsynaptic membrane currents, defined as the sum all transmembrane currents (capacitive, ionic, leak), but excluding input current. Black horizontal line in C indicates location of excitatory input. The spatiotemporal distribution of membrane currents (source currents) differ in the two models, but dipole-like V_e patterns remains qualitatively similar.

In fitting our model to data, we chose parameter values, when possible, based on known properties of MSO neurons. Some parameters, however, were chosen based on numerical experimentation and took different values for contralateral and ipsilateral responses, and for onset and ongoing responses, for the same recording session (Table 1). There are several known differences between contralateral and ipsilateral inputs that may explain the need for different parameter sets. Excitatory inputs to MSO are derived from the same cell types regardless of the ear of stimulation (spherical busy cells; Cant and Casseday, 1986), but the projections of axons from the ipsilateral side are more complex and irregular than those from the contralateral side (Karino et al., 2011). Inhibition evoked by ipsilateral stimuli [which arrives via the lateral nucleus of the trapezoid body (LNTB)] may be less temporally precise than inhibition evoked by contralateral stimuli [which arrives via the medial nucleus of the trapezoid body (MNTB)]. MNTB neurons in cat can show enhanced phase-locking (relative to spike-timing of auditory nerve fibers) as they relay spikes to MSO (Smith et al., 1998). LNTB is a more diverse nucleus than the MNTB (Spirou et al., 1998), but recent evidence in gerbil shows that LNTB neurons can also show enhanced phase-locking, although this was not the case for all LNTB neurons studied (Franken et al., 2015). Although the matter is certainly not settled, overall it seems unlikely that the LNTB provides the degree of temporally-precise inhibition to MSO that is conveyed by MNTB. The likelihood that contralateral inputs to MSO are more orderly arranged and more precisely timed may also explain why, in most cases, simulation results were more accurate for contralateral responses than for ipsilateral responses (Fig. 9). The need for different parameter sets to model ipsilateral and contralateral responses may also reflect possible mismatches in the frequency tuning of ipsilateral and contralateral inputs to MSO neurons (Shamma et al., 1989); see for instance Day and Semple (2011) and Benichoux et al. (2015) for consideration of this “stereausis” phenomenon. Differences in onset and ongoing responses may reflect a surge of well timed excitation at stimulus onset provided by spherical bushy cells (Cant and Casseday, 1986). Spherical bushy responses to tones are marked by strong onset responses followed by adaptation, both in terms of spike probability and spike-time precision (Smith et al., 1993; Joris et al., 1994a,b).

We presented data from five recording sessions (of 14 total) that we identified as sharing similar spatiotemporal V_e response patterns. We excluded the other datasets for a number of reasons. In two MSOs, the electrode track seemed to miss the MSO so that spatial patterns of V_e responses were only partially recorded. In a third case, the data could have been modeled with our method, but responses were not collected for tone frequencies >1500 Hz. Among the six remaining datasets, most were excluded because V_e responses exhibited idiosyncratic changes with tone frequency or side of stimulation. For example, in three cases, the spatial patterns of V_e responses to low-frequency tones differed from V_e responses to high-frequency tones suggesting (possibly) the existence of multiple populations of neurons that generate distinct neurophonic responses in a frequency-dependent manner. Additional work is required to account for idiosyncratic V_e responses that do not match the “typical” responses studied here (Figs. 7, 8).

Field potentials as a route for further study of MSO

Field potentials (V_e) are an important source of data for surveying neural activity in vivo. These measures resist simple interpretation because there is no direct way to identify the neural generators of V_e. Our approach, inspired by pioneering work of Rall and Shepherd (1968) and others, illustrates that we can make plausible connections between synaptic inputs to MSO and sound-evoked brainstem field potentials using a physiologically-informed mathematical model.

We highlighted how soma-targeting inhibition in MSO can create non-dipole-features of auditory neurophonic responses. The role of inhibition in MSO processing of binaural inputs is a topic of considerable interest (for reviews, see Grothe, 2003; Joris and Yin, 2007). In single-unit in vivo preparations, IPSCs cannot be directly recorded (but IPSPs can be detected; Franken et al., 2015). Recordings of V_e may provide a useful alternative perspective. For example, our model-based analysis of V_e data suggest that inhibition may be particularly prominent in low-frequency sustained responses (Fig. 9) and may precede excitation (Fig. 12). In vitro recordings in principal cells of the MSO show that inhibitory events precede excitatory events (Roberts et al., 2013). There is indirect evidence that inhibition acts to delay the effect of incoming excitatory inputs and thereby shift the tuning of MSO neurons to interaural time differences in vivo (Brand et al., 2002; Myoga et al., 2014; but see van der Heijden et al., 2013; Franken et al., 2015).

The interaural time difference at which an MSO neuron fires maximally (“best delay”) has implications for how sound source location is represented across the population of MSO neurons. There is debate (Goodman et al., 2013; Harper et al., 2014) whether the distribution of best delays covers the range of interaural delays imposed by an animal's head size, as envisioned by Jeffress (1948), or even reflects acoustical regularities in the animal's environment (Benichoux et al., 2015), or rather whether MSO neurons are grouped into subpopulations clustered around a distinct set of best delays (Harper and McAlpine, 2004). It has been proposed that inhibitory inputs that precede excitatory inputs alter the tuning of best delays, so understanding the relative dynamics of inhibitory and excitatory inputs to MSO is of great importance to current theories of spatial hearing. Our work shows that model-assisted interpretation of V_e responses may aid in determining the timing of inhibitory inputs.

Further understanding of the neurophonic will also clarify to what extent, if any, the neurophonic can operate as a mechanism for non-synaptic coupling between nearby MSO neurons. In simulations, endogenous V_e generated by MSO neurons can modulate spike timing and thresholds for spike generation (Goldwyn and Rinzel, 2016). Accurate models of V_e generation refine our understanding of how endogenous field potentials influence neural activity.

We studied brainstem field potentials that were recorded intracranially (an invasive procedure). There are common electrophysiological measures of auditory brainstem activity that can be obtained non-invasively. The auditory brainstem response is an important tool for diagnosing auditory function in normal and impaired human hearing. Moreover, phase-locked neural potentials can be measured near the cochlea and on the scalp and have been proposed as tools to investigate human hearing (Snyder and Schreiner, 1984; Kuwada et al., 1986; Shaheen et al., 2015; Verschooten et al., 2015). Presumably, these signals and the auditory neurophonic share neural generators (Caird et al., 1985; Sontheimer et al., 1985). Our observations regarding possible connections between neurophonic responses and MSO neural activity may, therefore, aid interpretation of these non-invasive diagnostic measures of neural activity in the human auditory system.

Footnotes

This work was supported by Grants from KU Leuven BOF (OT-14-118) and Research Foundation Flanders (G0A1113N and G091214N) to P.X.J.
The authors declare no competing financial interests.
Correspondence should be addressed to Dr. Joshua H. Goldwyn, Department of Mathematics and Statistics, Swarthmore College, 500 College Avenue, Swarthmore, PA 19081. jhgoldwyn{at}gmail.com

References

↵
2. Bartiromo R,
3. De Vincenzi M
(2016) Electrical measurements in the laboratory practice: undergraduate lecture notes in physics. Switzerland: Springer.
↵
2. Baumann VJ,
3. Lehnert S,
4. Leibold C,
5. Koch U
(2013) Tonotopic organization of the hyperpolarization-activated current (Ih) in the mammalian medial superior olive. Front Neural Circuits 7:117. doi:10.3389/fncir.2013.00117 pmid:23874271
OpenUrl CrossRef PubMed
↵
2. Benichoux V,
3. Fontaine B,
4. Franken TP,
5. Karino S,
6. Joris PX,
7. Brette R
(2015) Neural tuning matches frequency-dependent time differences between the ears. eLife 4:e06072. doi:10.7554/eLife.06072 pmid:25915620
OpenUrl Abstract/FREE Full Text
↵
2. Biedenbach MA,
3. Freeman WJ
(1964) Click-evoked potential map from the superior olivary nucleus. Am J Physiol 206:1408–1414. pmid:14193563
OpenUrl PubMed
↵
2. Bojanowski T,
3. Hu K,
4. Schwarz DW
(1989) Analogue signal representation in the medial superior olive of the cat. J Otolaryngol 18:3–9. pmid:2921783
OpenUrl PubMed
↵
2. Brand A,
3. Behrend O,
4. Marquardt T,
5. McAlpine D,
6. Grothe B
(2002) Precise inhibition is essential for microsecond interaural time difference coding. Nature 417:543–547. doi:10.1038/417543a pmid:12037566
OpenUrl CrossRef PubMed
↵
2. Bremen P,
3. Joris PX
(2013) Axonal recordings from medial superior olive neurons obtained from the lateral lemniscus of the chinchilla (Chinchilla laniger). J Neurosci 33:17506–17518. doi:10.1523/JNEUROSCI.1518-13.2013 pmid:24174683
OpenUrl Abstract/FREE Full Text
↵
2. Buzsáki G,
3. Anastassiou CA,
4. Koch C
(2012) The origin of extracellular fields and currents: EEG, ECoG, LFP and spikes. Nat Rev Neurosci 13:407–420. doi:10.1038/nrn3241 pmid:22595786
OpenUrl CrossRef PubMed
↵
2. Caird D,
3. Sontheimer D,
4. Klinke R
(1985) Intra-and extracranially recorded auditory evoked potentials in the cat: I. Source location and binaural interaction. Electroencephalogr Clin Neurophysiol 61:50–60. doi:10.1016/0013-4694(85)91072-7 pmid:2408863
OpenUrl CrossRef PubMed
↵
2. Cant NB,
3. Casseday JH
(1986) Projections from the anteroventral cochlear nucleus to the lateral and medial superior olivary nuclei. J Comp Neurol 247:457–476. doi:10.1002/cne.902470406 pmid:3722446
OpenUrl CrossRef PubMed
↵
2. Carr CE,
3. Shah S,
4. McColgan T,
5. Ashida G,
6. Kuokkanen PT,
7. Brill S,
8. Kempter R,
9. Wagner H
(2015) Maps of interaural delay in the owl's nucleus laminaris. J Neurophysiol 114:1862–1873. doi:10.1152/jn.00644.2015 pmid:26224776
OpenUrl Abstract/FREE Full Text
↵
2. Clark GM,
3. Dunlop CW
(1968) Field potentials in the cat medial superior olivary nucleus. Exp Neurol 20:31–42. doi:10.1016/0014-4886(68)90122-2 pmid:5637114
OpenUrl CrossRef PubMed
↵
2. Couchman K,
3. Grothe B,
4. 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 pmid:21159981
OpenUrl Abstract/FREE Full Text
↵
2. Couchman K,
3. Grothe B,
4. Felmy F
(2012) Functional localization of neurotransmitter receptors and synaptic inputs to mature neurons of the medial superior olive. J Neurophysiol 107:1186–1198. doi:10.1152/jn.00586.2011 pmid:22131383
OpenUrl Abstract/FREE Full Text
↵
2. Day ML,
3. Semple MN
(2011) Frequency-dependent interaural delays in the medial superior olive: implications for interaural cochlear delays. J Neurophysiol 106:1985–1999. doi:10.1152/jn.00131.2011 pmid:21775710
OpenUrl Abstract/FREE Full Text
↵
2. Fernández-Ruiz A,
3. Muñoz S,
4. Sancho M,
5. Makarova J,
6. Makarov VA,
7. Herreras O
(2013) Cytoarchitectonic and dynamic origins of giant positive local field potentials in the dentate gyrus. J Neurosci 33:15518–15532. doi:10.1523/JNEUROSCI.0338-13.2013 pmid:24068819
OpenUrl Abstract/FREE Full Text
↵
2. Fischl MJ,
3. Combs TD,
4. Klug A,
5. Grothe B,
6. Burger RM
(2012) Modulation of synaptic input by GABA_B receptors improves coincidence detection for computation of sound location. J Physiol 590:3047–3066. doi:10.1113/jphysiol.2011.226233 pmid:22473782
OpenUrl CrossRef PubMed
↵
2. Franken TP,
3. Roberts MT,
4. Wei L,
5. Golding NL,
6. Joris PX
(2015) In vivo coincidence detection in mammalian sound localization generates phase delays. Nat Neurosci 18:444–452. doi:10.1038/nn.3948 pmid:25664914
OpenUrl CrossRef PubMed
↵
2. Galambos R,
3. Schwartzkopff J,
4. Rupert A
(1959) Microelectrode study of superior olivary nuclei. Am J Physiol 197:527–536. pmid:13826000
OpenUrl Abstract/FREE Full Text
↵
2. Goldberg JM,
3. Brown PB
(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. pmid:5810617
OpenUrl FREE Full Text
↵
2. Golding NL,
3. Oertel D
(2012) Synaptic integration in dendrites: exceptional need for speed. J Physiol 590:5563–5569. doi:10.1113/jphysiol.2012.229328 pmid:22930273
OpenUrl CrossRef PubMed
↵
2. Goldwyn JH,
3. Rinzel J
(2016) Neuronal coupling by endogenous electric fields: cable theory and applications to coincidence detector neurons in the auditory brain stem. J Neurophysiol 115:2033–2051. doi:10.1152/jn.00780.2015 pmid:26823512
OpenUrl Abstract/FREE Full Text
↵
2. Goldwyn JH,
3. Mc Laughlin M,
4. Verschooten E,
5. Joris PX,
6. Rinzel J
(2014) A model of the medial superior olive explains spatiotemporal features of local field potentials. J Neurosci 34:11705–11722. doi:10.1523/JNEUROSCI.0175-14.2014 pmid:25164666
OpenUrl Abstract/FREE Full Text
↵
2. Goodman DF,
3. Benichoux V,
4. Brette R
(2013) Decoding neural responses to temporal cues for sound localization. eLife 2:e01312. doi:10.7554/eLife.01312 pmid:24302571
OpenUrl Abstract/FREE Full Text
↵
2. Grothe B
(1994) Interaction of excitation and inhibition in processing of pure tone and amplitude-modulated stimuli in the medial superior olive of the mustached bat. J Neurophysiol 71:706–721. pmid:8176433
OpenUrl Abstract/FREE Full Text
↵
2. Grothe B
(2003) New roles for synaptic inhibition in sound localization. Nat Rev Neurosci 4:540–550. doi:10.1038/nrn1136 pmid:12838329
OpenUrl CrossRef PubMed
↵
2. Grothe B,
3. Sanes DH
(1994) Synaptic inhibition influences the temporal coding properties of medial superior olivary neurons: an in vitro study. J Neurosci 14:1701–1709. pmid:8126564
OpenUrl Abstract
↵
2. Guinan JJ,
3. Guinan SS,
4. Norris BE
(1972a) Single auditory units in the superior olivary complex: I. Responses to sounds and classifications based on physiological properties. Int J Neurosci 4:101–120. doi:10.3109/00207457209147165
OpenUrl CrossRef
↵
2. Guinan JJ,
3. Norris BE,
4. Guinan SS
(1972b) Single auditory units in the superior olivary complex: II. Locations of unit categories and tonotopic organization. Int J Neurosci 4:147–166. doi:10.3109/00207457209164756
OpenUrl CrossRef
↵
2. Harper NS,
3. McAlpine D
(2004) Optimal neural population coding of an auditory spatial cue. Nature 430:682–686. doi:10.1038/nature02768 pmid:15295602
OpenUrl CrossRef PubMed
↵
2. Harper NS,
3. Scott BH,
4. Semple MN,
5. McAlpine D
(2014) The neural code for auditory space depends on sound frequency and head size in an optimal manner. PLoS One 9:e108154. doi:10.1371/journal.pone.0108154 pmid:25372405
OpenUrl CrossRef PubMed
↵
2. Herreras O
(2016) Local field potentials: myths and misunderstandings. Front Neural Circuits 10:101. doi:10.3389/fncir.2016.00101 pmid:28018180
OpenUrl CrossRef PubMed
↵
2. Jeffress LA
(1948) A place theory of sound localization. J Comp Physiol Psychol 41:35–39. doi:10.1037/h0061495 pmid:18904764
OpenUrl CrossRef PubMed
↵
2. Jercog PE,
3. Svirskis G,
4. Kotak VC,
5. Sanes DH,
6. 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 pmid:20613857
OpenUrl CrossRef PubMed
↵
2. Joris P,
3. Yin TC
(2007) A matter of time: internal delays in binaural processing. Trends Neurosci 30:70–78. doi:10.1016/j.tins.2006.12.004 pmid:17188761
OpenUrl CrossRef PubMed
↵
2. Joris PX,
3. Carney LH,
4. Smith PH,
5. Yin TC
(1994a) Enhancement of neural synchronization in the anteroventral cochlear nucleus: I. Responses to tones at the characteristic frequency. J Neurophysiol 71:1022–1036. pmid:8201399
OpenUrl Abstract/FREE Full Text
↵
2. Joris PX,
3. Smith PH,
4. Yin TC
(1994b) Enhancement of neural synchronization in the anteroventral cochlear nucleus: II. Responses in the tuning curve tail. J Neurophysiol 71:1037–1051. pmid:8201400
OpenUrl Abstract/FREE Full Text
↵
2. Kajikawa Y,
3. Schroeder CE
(2015) Generation of field potentials and modulation of their dynamics through volume integration of cortical activity. J Neurophysiol 113:339–351. doi:10.1152/jn.00914.2013 pmid:25274348
OpenUrl Abstract/FREE Full Text
↵
2. Kapfer C,
3. Seidl AH,
4. Schweizer H,
5. Grothe B
(2002) Experience-dependent refinement of inhibitory inputs to auditory coincidence-detector neurons. Nat Neurosci 5:247–253. doi:10.1038/nn810 pmid:11850629
OpenUrl CrossRef PubMed
↵
2. Karino S,
3. Smith PH,
4. Yin TC,
5. Joris PX
(2011) Axonal branching patterns as sources of delay in the mammalian auditory brainstem: a re-examination. J Neurosci 31:3016–3031. doi:10.1523/JNEUROSCI.5175-10.2011 pmid:21414923
OpenUrl Abstract/FREE Full Text
↵
2. Kido K
(2015) Digital Fourier analysis: fundamentals. New York: Springer.
↵
2. Klee M,
3. Rall W
(1977) Computed potentials of cortically arranged populations of neurons. J Neurophysiol 40:647–666. pmid:874533
OpenUrl FREE Full Text
↵
2. Köppl C,
3. Carr CE
(2008) Maps of interaural time difference in the chicken's brainstem nucleus laminaris. Biol Cybern 98:541–559. doi:10.1007/s00422-008-0220-6 pmid:18491165
OpenUrl CrossRef PubMed
↵
2. Kuokkanen PT,
3. Wagner H,
4. Ashida G,
5. Carr CE,
6. Kempter R
(2010) On the origin of the extracellular field potential in the nucleus laminaris of the barn owl (Tyto alba). J Neurophysiol 104:2274–2290. doi:10.1152/jn.00395.2010 pmid:20685926
OpenUrl Abstract/FREE Full Text
↵
2. Kuwada S,
3. Batra R,
4. Maher VL
(1986) Scalp potentials of normal and hearing-impaired subjects in response to sinusoidally amplitude-modulated tones. Hear Res 21:179–192. doi:10.1016/0378-5955(86)90038-9 pmid:3700256
OpenUrl CrossRef PubMed
↵
2. Lindén H,
3. Tetzlaff T,
4. Potjans TC,
5. Pettersen KH,
6. Grn S,
7. Diesmann M,
8. Einevoll GT
(2011) Modeling the spatial reach of the LFP. Neuron 72:859–872. doi:10.1016/j.neuron.2011.11.006 pmid:22153380
OpenUrl CrossRef PubMed
↵
2. Mathews PJ,
3. Jercog PE,
4. Rinzel J,
5. Scott LL,
6. Golding NL
(2010) Control of submillisecond synaptic timing in binaural coincidence detectors by K_v1 channels. Nat Neurosci 13:601–609. doi:10.1038/nn.2530 pmid:20364143
OpenUrl CrossRef PubMed
↵
2. Mc Laughlin M,
3. Verschooten E,
4. Joris PX
(2010) Oscillatory dipoles as a source of phase shifts in field potentials in the mammalian auditory brainstem. J Neurosci 30:13472–13487. doi:10.1523/JNEUROSCI.0294-10.2010 pmid:20926673
OpenUrl Abstract/FREE Full Text
↵
2. Mitzdorf U
(1985) Current source-density method and application in cat cerebral cortex: investigation of evoked potentials and EEG phenomena. Physiol Rev 65:37–100. pmid:3880898
OpenUrl FREE Full Text
↵
2. Myoga MH,
3. Lehnert S,
4. Leibold C,
5. Felmy F,
6. Grothe B
(2014) Glycinergic inhibition tunes coincidence detection in the auditory brainstem. Nat Commun 5:3790. doi:10.1038/ncomms4790 pmid:24804642
OpenUrl CrossRef PubMed
↵
2. Ness TV,
3. Remme MW,
4. Einevoll GT
(2016) Active subthreshold dendritic conductances shape the local field potential. J Physiol 594:3809–3825. doi:10.1113/JP272022 pmid:27079755
OpenUrl CrossRef PubMed
↵
2. Nicholson C,
3. Llinas R
(1971) Field potentials in the alligator cerebellum and theory of their relationship to purkinje cell dendritic spikes. J Neurophysiol 34:509–531. pmid:4329777
OpenUrl FREE Full Text
↵
2. Pecka M,
3. Brand A,
4. Behrend O,
5. Grothe B
(2008) Interaural time difference processing in the mammalian medial superior olive: the role of glycinergic inhibition. J Neurosci 28:6914–6925. doi:10.1523/JNEUROSCI.1660-08.2008 pmid:18596166
OpenUrl Abstract/FREE Full Text
↵
2. Pettersen KH,
3. Einevoll GT
(2008) Amplitude variability and extracellular low-pass filtering of neuronal spikes. Biophys J 94:784–802. doi:10.1529/biophysj.107.111179 pmid:17921225
OpenUrl CrossRef PubMed
↵
2. Plauška A,
3. Borst JG,
4. van der Heijden M
(2016) Predicting binaural responses from monaural responses in the gerbil medial superior olive. J Neurophysiol 115:2950–2963. doi:10.1152/jn.01146.2015 pmid:27009164
OpenUrl Abstract/FREE Full Text
↵
2. Rall W
(1977) Core conductor theory and cable properties of neurons. In: Handbook of physiology: cellular biology of neurons, pp 39–97. Bethesda, MD: American Physiological Society.
↵
2. Rall W,
3. Shepherd GM
(1968) Theoretical reconstruction of field potentials and dendrodendritic synaptic interactions in olfactory bulb. J Neurophysiol 31:884–915. pmid:5710539
OpenUrl FREE Full Text
↵
2. Rautenberg PL,
3. Grothe B,
4. Felmy F
(2009) Quantification of the three-dimensional morphology of coincidence detector neurons in the medial superior olive of gerbils during late postnatal development. J Comp Neurol 517:385–396. doi:10.1002/cne.22166 pmid:19760600
OpenUrl CrossRef PubMed
↵
2. Reimann MW,
3. Anastassiou CA,
4. Perin R,
5. Hill SL,
6. Markram H,
7. Koch C
(2013) A biophysically detailed model of neocortical local field potentials predicts the critical role of active membrane currents. Neuron 79:375–390. doi:10.1016/j.neuron.2013.05.023 pmid:23889937
OpenUrl CrossRef PubMed
↵
2. Remme MW,
3. Donato R,
4. Mikiel-Hunter J,
5. Ballestero JA,
6. Foster S,
7. Rinzel J,
8. McAlpine D
(2014) Subthreshold resonance properties contribute to the efficient coding of auditory spatial cues. Proc Natl Acad Sci U S A 111:E2339–E2348. doi:10.1073/pnas.1316216111 pmid:24843153
OpenUrl Abstract/FREE Full Text
↵
2. Roberts MT,
3. Seeman SC,
4. Golding NL
(2013) A mechanistic understanding of the role of feedforward inhibition in the mammalian sound localization circuitry. Neuron 78:923–935. doi:10.1016/j.neuron.2013.04.022 pmid:23764291
OpenUrl CrossRef PubMed
↵
2. Scott LL,
3. Hage TA,
4. Golding NL
(2007) Weak action potential backpropagation is associated with high-frequency axonal firing capability in principal neurons of the gerbil medial superior olive. J Physiol 583:647–661. doi:10.1113/jphysiol.2007.136366 pmid:17627992
OpenUrl CrossRef PubMed
↵
2. Scott LL,
3. Mathews PJ,
4. Golding NL
(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 pmid:20147532
OpenUrl Abstract/FREE Full Text
↵
2. Shaheen LA,
3. Valero MD,
4. Liberman MC
(2015) Towards a diagnosis of cochlear neuropathy with envelope following responses. J Assoc Res Otolaryngol 16:727–745. doi:10.1007/s10162-015-0539-3 pmid:26323349
OpenUrl CrossRef PubMed
↵
2. Shamma SA,
3. Shen NM,
4. Gopalaswamy P
(1989) Stereausis: binaural processing without neural delays. J Acoust Soc Am 86:989–1006. doi:10.1121/1.398734 pmid:2794252
OpenUrl CrossRef PubMed
↵
2. Smith PH,
3. Joris PX,
4. Yin TC
(1993) Projections of physiologically characterized spherical bushy cell axons from the cochlear nucleus of the cat: evidence for delay lines to the medial superior olive. J Comp Neurol 331:245–260. doi:10.1002/cne.903310208 pmid:8509501
OpenUrl CrossRef PubMed
↵
2. Smith PH,
3. Joris PX,
4. Yin TC
(1998) Anatomy and physiology of principal cells of the medial nucleus of the trapezoid body (MNTB) of the cat. J Neurophysiol 79:3127–3142. pmid:9636113
OpenUrl Abstract/FREE Full Text
↵
2. Snyder RL,
3. Schreiner CE
(1984) The auditory neurophonic: basic properties. Hear Res 15:261–280. doi:10.1016/0378-5955(84)90033-9 pmid:6501114
OpenUrl CrossRef PubMed
↵
2. Sontheimer D,
3. Caird D,
4. Klinke R
(1985) Intra- and extracranially recorded auditory evoked potentials in the cat: II. Effects of interaural time and intensity differences. Electroencephalogr Clin Neurophysiol 61:539–547. doi:10.1016/0013-4694(85)90973-3 pmid:2415329
OpenUrl CrossRef PubMed
↵
2. Spirou GA,
3. Rowland KC,
4. Berrebi AS
(1998) Ultrastructure of neurons and large synaptic terminals in the lateral nucleus of the trapezoid body of the cat. J Comp Neurol 398:257–272. doi:10.1002/(SICI)1096-9861(19980824)398:2%3C257::AID-CNE7%3E3.0.CO%3B2-%23 pmid:9700570
OpenUrl CrossRef PubMed
↵
2. Stotler WA
(1953) An experimental study of the cells and connections of the superior olivary complex of the cat. J Comp Neurol 98:401–431. doi:10.1002/cne.900980303 pmid:13069629
OpenUrl CrossRef PubMed
↵
2. Sullivan WE,
3. Konishi M
(1986) Neural map of interaural phase difference in the owl's brainstem. Proc Natl Acad Sci U S A 83:8400–8404. doi:10.1073/pnas.83.21.8400 pmid:3022292
OpenUrl Abstract/FREE Full Text
↵
2. Svirskis G,
3. Dodla R,
4. Rinzel J
(2003) Subthreshold outward currents enhance temporal integration in auditory neurons. Biol Cybern 89:333–340. doi:10.1007/s00422-003-0438-2 pmid:14669013
OpenUrl CrossRef PubMed
↵
2. Tsuchitani C,
3. Boudreau JC
(1964) Wave activity in the superior olivary complex of the cat. J Neurophysiol 27:814–827. pmid:14205006
OpenUrl FREE Full Text
↵
2. Tuckwell HC
(1988) Introduction to theoretical neurobiology: volume 1 linear cable theory and dendritic structure. Cambridge, UK: Cambridge UP.
↵
2. van der Heijden M,
3. Lorteije JA,
4. Plauška A,
5. Roberts MT,
6. Golding NL,
7. Borst JG
(2013) Directional hearing by linear summation of binaural inputs at the medial superior olive. Neuron 78:936–948. doi:10.1016/j.neuron.2013.04.028 pmid:23764292
OpenUrl CrossRef PubMed
↵
2. Verschooten E,
3. Robles L,
4. Joris PX
(2015) Assessment of the limits of neural phase-locking using mass potentials. J Neurosci 35:2255–2268. doi:10.1523/JNEUROSCI.2979-14.2015 pmid:25653380
OpenUrl Abstract/FREE Full Text
↵
2. Weinberger NM,
3. Kitzes LM,
4. Goodman DA
(1970) Some characteristics of the “auditory neurophonic”. Experientia 26:46–48. doi:10.1007/BF01900383 pmid:5412293
OpenUrl CrossRef PubMed
↵
2. Yin TC,
3. Chan JC
(1990) Interaural time sensitivity in medial superior olive of cat. J Neurophysiol 64:465–488. pmid:2213127
OpenUrl Abstract/FREE Full Text

In this issue

View Full Page PDF

Citation Tools

Respond to this article

Request Permissions

Keywords

Cited By...

Research Articles

Show more Research Articles

Systems/Circuits

Show more Systems/Circuits

[1] ↵

Bartiromo R,
De Vincenzi M
(2016) Electrical measurements in the laboratory practice: undergraduate lecture notes in physics. Switzerland: Springer.

[3] Bartiromo R,

[4] De Vincenzi M

[5] ↵

Baumann VJ,
Lehnert S,
Leibold C,
Koch U
(2013) Tonotopic organization of the hyperpolarization-activated current (Ih) in the mammalian medial superior olive. Front Neural Circuits 7:117. doi:10.3389/fncir.2013.00117 pmid:23874271
OpenUrl CrossRef PubMed

[7] Baumann VJ,

[8] Lehnert S,

[9] Leibold C,

[10] Koch U

[11] ↵

Benichoux V,
Fontaine B,
Franken TP,
Karino S,
Joris PX,
Brette R
(2015) Neural tuning matches frequency-dependent time differences between the ears. eLife 4:e06072. doi:10.7554/eLife.06072 pmid:25915620
OpenUrl Abstract/FREE Full Text

[13] Benichoux V,

[14] Fontaine B,

[15] Franken TP,

[16] Karino S,

[17] Joris PX,

[18] Brette R

[19] ↵

Biedenbach MA,
Freeman WJ
(1964) Click-evoked potential map from the superior olivary nucleus. Am J Physiol 206:1408–1414. pmid:14193563
OpenUrl PubMed

[21] Biedenbach MA,

[22] Freeman WJ

[23] ↵

Bojanowski T,
Hu K,
Schwarz DW
(1989) Analogue signal representation in the medial superior olive of the cat. J Otolaryngol 18:3–9. pmid:2921783
OpenUrl PubMed

[25] Bojanowski T,

[26] Hu K,

[27] Schwarz DW

[28] ↵

Brand A,
Behrend O,
Marquardt T,
McAlpine D,
Grothe B
(2002) Precise inhibition is essential for microsecond interaural time difference coding. Nature 417:543–547. doi:10.1038/417543a pmid:12037566
OpenUrl CrossRef PubMed

[30] Brand A,

[31] Behrend O,

[32] Marquardt T,

[33] McAlpine D,

[34] Grothe B

[35] ↵

Bremen P,
Joris PX
(2013) Axonal recordings from medial superior olive neurons obtained from the lateral lemniscus of the chinchilla (Chinchilla laniger). J Neurosci 33:17506–17518. doi:10.1523/JNEUROSCI.1518-13.2013 pmid:24174683
OpenUrl Abstract/FREE Full Text

[37] Bremen P,

[38] Joris PX

[39] ↵

Buzsáki G,
Anastassiou CA,
Koch C
(2012) The origin of extracellular fields and currents: EEG, ECoG, LFP and spikes. Nat Rev Neurosci 13:407–420. doi:10.1038/nrn3241 pmid:22595786
OpenUrl CrossRef PubMed

[41] Buzsáki G,

[42] Anastassiou CA,

[43] Koch C

[44] ↵

Caird D,
Sontheimer D,
Klinke R
(1985) Intra-and extracranially recorded auditory evoked potentials in the cat: I. Source location and binaural interaction. Electroencephalogr Clin Neurophysiol 61:50–60. doi:10.1016/0013-4694(85)91072-7 pmid:2408863
OpenUrl CrossRef PubMed

[46] Caird D,

[47] Sontheimer D,

[48] Klinke R

[49] ↵

Cant NB,
Casseday JH
(1986) Projections from the anteroventral cochlear nucleus to the lateral and medial superior olivary nuclei. J Comp Neurol 247:457–476. doi:10.1002/cne.902470406 pmid:3722446
OpenUrl CrossRef PubMed

[51] Cant NB,

[52] Casseday JH

[53] ↵

Carr CE,
Shah S,
McColgan T,
Ashida G,
Kuokkanen PT,
Brill S,
Kempter R,
Wagner H
(2015) Maps of interaural delay in the owl's nucleus laminaris. J Neurophysiol 114:1862–1873. doi:10.1152/jn.00644.2015 pmid:26224776
OpenUrl Abstract/FREE Full Text

[55] Carr CE,

[56] Shah S,

[57] McColgan T,

[58] Ashida G,

[59] Kuokkanen PT,

[60] Brill S,

[61] Kempter R,

[62] Wagner H

[63] ↵

Clark GM,
Dunlop CW
(1968) Field potentials in the cat medial superior olivary nucleus. Exp Neurol 20:31–42. doi:10.1016/0014-4886(68)90122-2 pmid:5637114
OpenUrl CrossRef PubMed

[65] Clark GM,

[66] Dunlop CW

[67] ↵

Couchman K,
Grothe B,
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 pmid:21159981
OpenUrl Abstract/FREE Full Text

[69] Couchman K,

[70] Grothe B,

[71] Felmy F

[72] ↵

Couchman K,
Grothe B,
Felmy F
(2012) Functional localization of neurotransmitter receptors and synaptic inputs to mature neurons of the medial superior olive. J Neurophysiol 107:1186–1198. doi:10.1152/jn.00586.2011 pmid:22131383
OpenUrl Abstract/FREE Full Text

[74] Couchman K,

[75] Grothe B,

[76] Felmy F

[77] ↵

Day ML,
Semple MN
(2011) Frequency-dependent interaural delays in the medial superior olive: implications for interaural cochlear delays. J Neurophysiol 106:1985–1999. doi:10.1152/jn.00131.2011 pmid:21775710
OpenUrl Abstract/FREE Full Text

[79] Day ML,

[80] Semple MN

[81] ↵

Fernández-Ruiz A,
Muñoz S,
Sancho M,
Makarova J,
Makarov VA,
Herreras O
(2013) Cytoarchitectonic and dynamic origins of giant positive local field potentials in the dentate gyrus. J Neurosci 33:15518–15532. doi:10.1523/JNEUROSCI.0338-13.2013 pmid:24068819
OpenUrl Abstract/FREE Full Text

[83] Fernández-Ruiz A,

[84] Muñoz S,

[85] Sancho M,

[86] Makarova J,

[87] Makarov VA,

[88] Herreras O

[89] ↵

Fischl MJ,
Combs TD,
Klug A,
Grothe B,
Burger RM
(2012) Modulation of synaptic input by GABA_B receptors improves coincidence detection for computation of sound location. J Physiol 590:3047–3066. doi:10.1113/jphysiol.2011.226233 pmid:22473782
OpenUrl CrossRef PubMed

[91] Fischl MJ,

[92] Combs TD,

[93] Klug A,

[94] Grothe B,

[95] Burger RM

[96] ↵

Franken TP,
Roberts MT,
Wei L,
Golding NL,
Joris PX
(2015) In vivo coincidence detection in mammalian sound localization generates phase delays. Nat Neurosci 18:444–452. doi:10.1038/nn.3948 pmid:25664914
OpenUrl CrossRef PubMed

[98] Franken TP,

[99] Roberts MT,

[100] Wei L,

[101] Golding NL,

[102] Joris PX

[103] ↵

Galambos R,
Schwartzkopff J,
Rupert A
(1959) Microelectrode study of superior olivary nuclei. Am J Physiol 197:527–536. pmid:13826000
OpenUrl Abstract/FREE Full Text

[105] Galambos R,

[106] Schwartzkopff J,

[107] Rupert A

[108] ↵

Goldberg JM,
Brown PB
(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. pmid:5810617
OpenUrl FREE Full Text

[110] Goldberg JM,

[111] Brown PB

[112] ↵

Golding NL,
Oertel D
(2012) Synaptic integration in dendrites: exceptional need for speed. J Physiol 590:5563–5569. doi:10.1113/jphysiol.2012.229328 pmid:22930273
OpenUrl CrossRef PubMed

[114] Golding NL,

[115] Oertel D

[116] ↵

Goldwyn JH,
Rinzel J
(2016) Neuronal coupling by endogenous electric fields: cable theory and applications to coincidence detector neurons in the auditory brain stem. J Neurophysiol 115:2033–2051. doi:10.1152/jn.00780.2015 pmid:26823512
OpenUrl Abstract/FREE Full Text

[118] Goldwyn JH,

[119] Rinzel J

[120] ↵

Goldwyn JH,
Mc Laughlin M,
Verschooten E,
Joris PX,
Rinzel J
(2014) A model of the medial superior olive explains spatiotemporal features of local field potentials. J Neurosci 34:11705–11722. doi:10.1523/JNEUROSCI.0175-14.2014 pmid:25164666
OpenUrl Abstract/FREE Full Text

[122] Goldwyn JH,

[123] Mc Laughlin M,

[124] Verschooten E,

[125] Joris PX,

[126] Rinzel J

[127] ↵

Goodman DF,
Benichoux V,
Brette R
(2013) Decoding neural responses to temporal cues for sound localization. eLife 2:e01312. doi:10.7554/eLife.01312 pmid:24302571
OpenUrl Abstract/FREE Full Text

[129] Goodman DF,

[130] Benichoux V,

[131] Brette R

[132] ↵

Grothe B
(1994) Interaction of excitation and inhibition in processing of pure tone and amplitude-modulated stimuli in the medial superior olive of the mustached bat. J Neurophysiol 71:706–721. pmid:8176433
OpenUrl Abstract/FREE Full Text

[134] Grothe B

[135] ↵

Grothe B
(2003) New roles for synaptic inhibition in sound localization. Nat Rev Neurosci 4:540–550. doi:10.1038/nrn1136 pmid:12838329
OpenUrl CrossRef PubMed

[137] Grothe B

[138] ↵

Grothe B,
Sanes DH
(1994) Synaptic inhibition influences the temporal coding properties of medial superior olivary neurons: an in vitro study. J Neurosci 14:1701–1709. pmid:8126564
OpenUrl Abstract

[140] Grothe B,

[141] Sanes DH

[142] ↵

Guinan JJ,
Guinan SS,
Norris BE
(1972a) Single auditory units in the superior olivary complex: I. Responses to sounds and classifications based on physiological properties. Int J Neurosci 4:101–120. doi:10.3109/00207457209147165
OpenUrl CrossRef

[144] Guinan JJ,

[145] Guinan SS,

[146] Norris BE

[147] ↵

Guinan JJ,
Norris BE,
Guinan SS
(1972b) Single auditory units in the superior olivary complex: II. Locations of unit categories and tonotopic organization. Int J Neurosci 4:147–166. doi:10.3109/00207457209164756
OpenUrl CrossRef

[149] Guinan JJ,

[150] Norris BE,

[151] Guinan SS

[152] ↵

Harper NS,
McAlpine D
(2004) Optimal neural population coding of an auditory spatial cue. Nature 430:682–686. doi:10.1038/nature02768 pmid:15295602
OpenUrl CrossRef PubMed

[154] Harper NS,

[155] McAlpine D

[156] ↵

Harper NS,
Scott BH,
Semple MN,
McAlpine D
(2014) The neural code for auditory space depends on sound frequency and head size in an optimal manner. PLoS One 9:e108154. doi:10.1371/journal.pone.0108154 pmid:25372405
OpenUrl CrossRef PubMed

[158] Harper NS,

[159] Scott BH,

[160] Semple MN,

[161] McAlpine D

[162] ↵

Herreras O
(2016) Local field potentials: myths and misunderstandings. Front Neural Circuits 10:101. doi:10.3389/fncir.2016.00101 pmid:28018180
OpenUrl CrossRef PubMed

[164] Herreras O

[165] ↵

Jeffress LA
(1948) A place theory of sound localization. J Comp Physiol Psychol 41:35–39. doi:10.1037/h0061495 pmid:18904764
OpenUrl CrossRef PubMed

[167] Jeffress LA

[168] ↵

Jercog PE,
Svirskis G,
Kotak VC,
Sanes DH,
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 pmid:20613857
OpenUrl CrossRef PubMed

[170] Jercog PE,

[171] Svirskis G,

[172] Kotak VC,

[173] Sanes DH,

[174] Rinzel J

[175] ↵

Joris P,
Yin TC
(2007) A matter of time: internal delays in binaural processing. Trends Neurosci 30:70–78. doi:10.1016/j.tins.2006.12.004 pmid:17188761
OpenUrl CrossRef PubMed

[177] Joris P,

[178] Yin TC

[179] ↵

Joris PX,
Carney LH,
Smith PH,
Yin TC
(1994a) Enhancement of neural synchronization in the anteroventral cochlear nucleus: I. Responses to tones at the characteristic frequency. J Neurophysiol 71:1022–1036. pmid:8201399
OpenUrl Abstract/FREE Full Text

[181] Joris PX,

[182] Carney LH,

[183] Smith PH,

[184] Yin TC

[185] ↵

Joris PX,
Smith PH,
Yin TC
(1994b) Enhancement of neural synchronization in the anteroventral cochlear nucleus: II. Responses in the tuning curve tail. J Neurophysiol 71:1037–1051. pmid:8201400
OpenUrl Abstract/FREE Full Text

[187] Joris PX,

[188] Smith PH,

[189] Yin TC

[190] ↵

Kajikawa Y,
Schroeder CE
(2015) Generation of field potentials and modulation of their dynamics through volume integration of cortical activity. J Neurophysiol 113:339–351. doi:10.1152/jn.00914.2013 pmid:25274348
OpenUrl Abstract/FREE Full Text

[192] Kajikawa Y,

[193] Schroeder CE

[194] ↵

Kapfer C,
Seidl AH,
Schweizer H,
Grothe B
(2002) Experience-dependent refinement of inhibitory inputs to auditory coincidence-detector neurons. Nat Neurosci 5:247–253. doi:10.1038/nn810 pmid:11850629
OpenUrl CrossRef PubMed

[196] Kapfer C,

[197] Seidl AH,

[198] Schweizer H,

[199] Grothe B

[200] ↵

Karino S,
Smith PH,
Yin TC,
Joris PX
(2011) Axonal branching patterns as sources of delay in the mammalian auditory brainstem: a re-examination. J Neurosci 31:3016–3031. doi:10.1523/JNEUROSCI.5175-10.2011 pmid:21414923
OpenUrl Abstract/FREE Full Text

[202] Karino S,

[203] Smith PH,

[204] Yin TC,

[205] Joris PX

[206] ↵

Kido K
(2015) Digital Fourier analysis: fundamentals. New York: Springer.

[208] Kido K

[209] ↵

Klee M,
Rall W
(1977) Computed potentials of cortically arranged populations of neurons. J Neurophysiol 40:647–666. pmid:874533
OpenUrl FREE Full Text

[211] Klee M,

[212] Rall W

[213] ↵

Köppl C,
Carr CE
(2008) Maps of interaural time difference in the chicken's brainstem nucleus laminaris. Biol Cybern 98:541–559. doi:10.1007/s00422-008-0220-6 pmid:18491165
OpenUrl CrossRef PubMed

[215] Köppl C,

[216] Carr CE

[217] ↵

Kuokkanen PT,
Wagner H,
Ashida G,
Carr CE,
Kempter R
(2010) On the origin of the extracellular field potential in the nucleus laminaris of the barn owl (Tyto alba). J Neurophysiol 104:2274–2290. doi:10.1152/jn.00395.2010 pmid:20685926
OpenUrl Abstract/FREE Full Text

[219] Kuokkanen PT,

[220] Wagner H,

[221] Ashida G,

[222] Carr CE,

[223] Kempter R

[224] ↵

Kuwada S,
Batra R,
Maher VL
(1986) Scalp potentials of normal and hearing-impaired subjects in response to sinusoidally amplitude-modulated tones. Hear Res 21:179–192. doi:10.1016/0378-5955(86)90038-9 pmid:3700256
OpenUrl CrossRef PubMed

[226] Kuwada S,

[227] Batra R,

[228] Maher VL

[229] ↵

Lindén H,
Tetzlaff T,
Potjans TC,
Pettersen KH,
Grn S,
Diesmann M,
Einevoll GT
(2011) Modeling the spatial reach of the LFP. Neuron 72:859–872. doi:10.1016/j.neuron.2011.11.006 pmid:22153380
OpenUrl CrossRef PubMed

[231] Lindén H,

[232] Tetzlaff T,

[233] Potjans TC,

[234] Pettersen KH,

[235] Grn S,

[236] Diesmann M,

[237] Einevoll GT

[238] ↵

Mathews PJ,
Jercog PE,
Rinzel J,
Scott LL,
Golding NL
(2010) Control of submillisecond synaptic timing in binaural coincidence detectors by K_v1 channels. Nat Neurosci 13:601–609. doi:10.1038/nn.2530 pmid:20364143
OpenUrl CrossRef PubMed

[240] Mathews PJ,

[241] Jercog PE,

[242] Rinzel J,

[243] Scott LL,

[244] Golding NL

[245] ↵

Mc Laughlin M,
Verschooten E,
Joris PX
(2010) Oscillatory dipoles as a source of phase shifts in field potentials in the mammalian auditory brainstem. J Neurosci 30:13472–13487. doi:10.1523/JNEUROSCI.0294-10.2010 pmid:20926673
OpenUrl Abstract/FREE Full Text

[247] Mc Laughlin M,

[248] Verschooten E,

[249] Joris PX

[250] ↵

Mitzdorf U
(1985) Current source-density method and application in cat cerebral cortex: investigation of evoked potentials and EEG phenomena. Physiol Rev 65:37–100. pmid:3880898
OpenUrl FREE Full Text

[252] Mitzdorf U

[253] ↵

Myoga MH,
Lehnert S,
Leibold C,
Felmy F,
Grothe B
(2014) Glycinergic inhibition tunes coincidence detection in the auditory brainstem. Nat Commun 5:3790. doi:10.1038/ncomms4790 pmid:24804642
OpenUrl CrossRef PubMed

[255] Myoga MH,

[256] Lehnert S,

[257] Leibold C,

[258] Felmy F,

[259] Grothe B

[260] ↵

Ness TV,
Remme MW,
Einevoll GT
(2016) Active subthreshold dendritic conductances shape the local field potential. J Physiol 594:3809–3825. doi:10.1113/JP272022 pmid:27079755
OpenUrl CrossRef PubMed

[262] Ness TV,

[263] Remme MW,

[264] Einevoll GT

[265] ↵

Nicholson C,
Llinas R
(1971) Field potentials in the alligator cerebellum and theory of their relationship to purkinje cell dendritic spikes. J Neurophysiol 34:509–531. pmid:4329777
OpenUrl FREE Full Text

[267] Nicholson C,

[268] Llinas R

[269] ↵

Pecka M,
Brand A,
Behrend O,
Grothe B
(2008) Interaural time difference processing in the mammalian medial superior olive: the role of glycinergic inhibition. J Neurosci 28:6914–6925. doi:10.1523/JNEUROSCI.1660-08.2008 pmid:18596166
OpenUrl Abstract/FREE Full Text

[271] Pecka M,

[272] Brand A,

[273] Behrend O,

[274] Grothe B

[275] ↵

Pettersen KH,
Einevoll GT
(2008) Amplitude variability and extracellular low-pass filtering of neuronal spikes. Biophys J 94:784–802. doi:10.1529/biophysj.107.111179 pmid:17921225
OpenUrl CrossRef PubMed

[277] Pettersen KH,

[278] Einevoll GT

[279] ↵

Plauška A,
Borst JG,
van der Heijden M
(2016) Predicting binaural responses from monaural responses in the gerbil medial superior olive. J Neurophysiol 115:2950–2963. doi:10.1152/jn.01146.2015 pmid:27009164
OpenUrl Abstract/FREE Full Text

[281] Plauška A,

[282] Borst JG,

[283] van der Heijden M

[284] ↵

Rall W
(1977) Core conductor theory and cable properties of neurons. In: Handbook of physiology: cellular biology of neurons, pp 39–97. Bethesda, MD: American Physiological Society.

[286] Rall W

[287] ↵

Rall W,
Shepherd GM
(1968) Theoretical reconstruction of field potentials and dendrodendritic synaptic interactions in olfactory bulb. J Neurophysiol 31:884–915. pmid:5710539
OpenUrl FREE Full Text

[289] Rall W,

[290] Shepherd GM

[291] ↵

Rautenberg PL,
Grothe B,
Felmy F
(2009) Quantification of the three-dimensional morphology of coincidence detector neurons in the medial superior olive of gerbils during late postnatal development. J Comp Neurol 517:385–396. doi:10.1002/cne.22166 pmid:19760600
OpenUrl CrossRef PubMed

[293] Rautenberg PL,

[294] Grothe B,

[295] Felmy F

[296] ↵

Reimann MW,
Anastassiou CA,
Perin R,
Hill SL,
Markram H,
Koch C
(2013) A biophysically detailed model of neocortical local field potentials predicts the critical role of active membrane currents. Neuron 79:375–390. doi:10.1016/j.neuron.2013.05.023 pmid:23889937
OpenUrl CrossRef PubMed

[298] Reimann MW,

[299] Anastassiou CA,

[300] Perin R,

[301] Hill SL,

[302] Markram H,

[303] Koch C

[304] ↵

Remme MW,
Donato R,
Mikiel-Hunter J,
Ballestero JA,
Foster S,
Rinzel J,
McAlpine D
(2014) Subthreshold resonance properties contribute to the efficient coding of auditory spatial cues. Proc Natl Acad Sci U S A 111:E2339–E2348. doi:10.1073/pnas.1316216111 pmid:24843153
OpenUrl Abstract/FREE Full Text

[306] Remme MW,

[307] Donato R,

[308] Mikiel-Hunter J,

[309] Ballestero JA,

[310] Foster S,

[311] Rinzel J,

[312] McAlpine D

[313] ↵

Roberts MT,
Seeman SC,
Golding NL
(2013) A mechanistic understanding of the role of feedforward inhibition in the mammalian sound localization circuitry. Neuron 78:923–935. doi:10.1016/j.neuron.2013.04.022 pmid:23764291
OpenUrl CrossRef PubMed

[315] Roberts MT,

[316] Seeman SC,

[317] Golding NL

[318] ↵

Scott LL,
Hage TA,
Golding NL
(2007) Weak action potential backpropagation is associated with high-frequency axonal firing capability in principal neurons of the gerbil medial superior olive. J Physiol 583:647–661. doi:10.1113/jphysiol.2007.136366 pmid:17627992
OpenUrl CrossRef PubMed

[320] Scott LL,

[321] Hage TA,

[322] Golding NL

[323] ↵

Scott LL,
Mathews PJ,
Golding NL
(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 pmid:20147532
OpenUrl Abstract/FREE Full Text

[325] Scott LL,

[326] Mathews PJ,

[327] Golding NL

[328] ↵

Shaheen LA,
Valero MD,
Liberman MC
(2015) Towards a diagnosis of cochlear neuropathy with envelope following responses. J Assoc Res Otolaryngol 16:727–745. doi:10.1007/s10162-015-0539-3 pmid:26323349
OpenUrl CrossRef PubMed

[330] Shaheen LA,

[331] Valero MD,

[332] Liberman MC

[333] ↵

Shamma SA,
Shen NM,
Gopalaswamy P
(1989) Stereausis: binaural processing without neural delays. J Acoust Soc Am 86:989–1006. doi:10.1121/1.398734 pmid:2794252
OpenUrl CrossRef PubMed

[335] Shamma SA,

[336] Shen NM,

[337] Gopalaswamy P

[338] ↵

Smith PH,
Joris PX,
Yin TC
(1993) Projections of physiologically characterized spherical bushy cell axons from the cochlear nucleus of the cat: evidence for delay lines to the medial superior olive. J Comp Neurol 331:245–260. doi:10.1002/cne.903310208 pmid:8509501
OpenUrl CrossRef PubMed

[340] Smith PH,

[341] Joris PX,

[342] Yin TC

[343] ↵

Smith PH,
Joris PX,
Yin TC
(1998) Anatomy and physiology of principal cells of the medial nucleus of the trapezoid body (MNTB) of the cat. J Neurophysiol 79:3127–3142. pmid:9636113
OpenUrl Abstract/FREE Full Text

[345] Smith PH,

[346] Joris PX,

[347] Yin TC

[348] ↵

Snyder RL,
Schreiner CE
(1984) The auditory neurophonic: basic properties. Hear Res 15:261–280. doi:10.1016/0378-5955(84)90033-9 pmid:6501114
OpenUrl CrossRef PubMed

[350] Snyder RL,

[351] Schreiner CE

[352] ↵

Sontheimer D,
Caird D,
Klinke R
(1985) Intra- and extracranially recorded auditory evoked potentials in the cat: II. Effects of interaural time and intensity differences. Electroencephalogr Clin Neurophysiol 61:539–547. doi:10.1016/0013-4694(85)90973-3 pmid:2415329
OpenUrl CrossRef PubMed

[354] Sontheimer D,

[355] Caird D,

[356] Klinke R

[357] ↵

Spirou GA,
Rowland KC,
Berrebi AS
(1998) Ultrastructure of neurons and large synaptic terminals in the lateral nucleus of the trapezoid body of the cat. J Comp Neurol 398:257–272. doi:10.1002/(SICI)1096-9861(19980824)398:2%3C257::AID-CNE7%3E3.0.CO%3B2-%23 pmid:9700570
OpenUrl CrossRef PubMed

[359] Spirou GA,

[360] Rowland KC,

[361] Berrebi AS

[362] ↵

Stotler WA
(1953) An experimental study of the cells and connections of the superior olivary complex of the cat. J Comp Neurol 98:401–431. doi:10.1002/cne.900980303 pmid:13069629
OpenUrl CrossRef PubMed

[364] Stotler WA

[365] ↵

Sullivan WE,
Konishi M
(1986) Neural map of interaural phase difference in the owl's brainstem. Proc Natl Acad Sci U S A 83:8400–8404. doi:10.1073/pnas.83.21.8400 pmid:3022292
OpenUrl Abstract/FREE Full Text

[367] Sullivan WE,

[368] Konishi M

[369] ↵

Svirskis G,
Dodla R,
Rinzel J
(2003) Subthreshold outward currents enhance temporal integration in auditory neurons. Biol Cybern 89:333–340. doi:10.1007/s00422-003-0438-2 pmid:14669013
OpenUrl CrossRef PubMed

[371] Svirskis G,

[372] Dodla R,

[373] Rinzel J

[374] ↵

Tsuchitani C,
Boudreau JC
(1964) Wave activity in the superior olivary complex of the cat. J Neurophysiol 27:814–827. pmid:14205006
OpenUrl FREE Full Text

[376] Tsuchitani C,

[377] Boudreau JC

[378] ↵

Tuckwell HC
(1988) Introduction to theoretical neurobiology: volume 1 linear cable theory and dendritic structure. Cambridge, UK: Cambridge UP.

[380] Tuckwell HC

[381] ↵

van der Heijden M,
Lorteije JA,
Plauška A,
Roberts MT,
Golding NL,
Borst JG
(2013) Directional hearing by linear summation of binaural inputs at the medial superior olive. Neuron 78:936–948. doi:10.1016/j.neuron.2013.04.028 pmid:23764292
OpenUrl CrossRef PubMed

[383] van der Heijden M,

[384] Lorteije JA,

[385] Plauška A,

[386] Roberts MT,

[387] Golding NL,

[388] Borst JG

[389] ↵

Verschooten E,
Robles L,
Joris PX
(2015) Assessment of the limits of neural phase-locking using mass potentials. J Neurosci 35:2255–2268. doi:10.1523/JNEUROSCI.2979-14.2015 pmid:25653380
OpenUrl Abstract/FREE Full Text

[391] Verschooten E,

[392] Robles L,

[393] Joris PX

[394] ↵

Weinberger NM,
Kitzes LM,
Goodman DA
(1970) Some characteristics of the “auditory neurophonic”. Experientia 26:46–48. doi:10.1007/BF01900383 pmid:5412293
OpenUrl CrossRef PubMed

[396] Weinberger NM,

[397] Kitzes LM,

[398] Goodman DA

[399] ↵

Yin TC,
Chan JC
(1990) Interaural time sensitivity in medial superior olive of cat. J Neurophysiol 64:465–488. pmid:2213127
OpenUrl Abstract/FREE Full Text

[401] Yin TC,

[402] Chan JC

Main menu

User menu

Search

Signatures of Somatic Inhibition and Dendritic Excitation in Auditory Brainstem Field Potentials

Abstract

Introduction

Materials and Methods

Neurophonic recordings.

Spatiotemporal V_e patterns.

Approximation of V_e spatiotemporal patterns in the Fourier domain.

Mathematical model.

Model fitting procedure.

Results

“Dipole” and “non-dipole” features in V_e responses

Simulated V_e patterns reproduce diverse and frequency-dependent neurophonic responses

Onset responses identify I₁ and I₂ as source and sink current, respectively

Characteristics of putative inhibition and putative excitation

Putative inhibition attenuates with tone frequency and putative excitation exhibits a “best frequency”

Putative inhibition precedes putative excitation

Consideration of alternate input configurations

Discussion

A minimal model of the MSO reproduces diverse and frequency-dependent neurophonic responses

A testable working hypothesis: inhibition explains non-dipole-like neurophonic responses

Evaluation of modeling approach

Field potentials as a route for further study of MSO

Footnotes

References

In this issue

Citation Manager Formats

Keywords

Responses to this article

Jump to comment:

Related Articles

Cited By...

More in this TOC Section

Research Articles

Systems/Circuits

Main menu

User menu

Search

Signatures of Somatic Inhibition and Dendritic Excitation in Auditory Brainstem Field Potentials

Abstract

Introduction

Materials and Methods

Neurophonic recordings.

Spatiotemporal Ve patterns.

Approximation of Ve spatiotemporal patterns in the Fourier domain.

Mathematical model.

Model fitting procedure.

Results

“Dipole” and “non-dipole” features in Ve responses

Simulated Ve patterns reproduce diverse and frequency-dependent neurophonic responses

Onset responses identify I1 and I2 as source and sink current, respectively

Characteristics of putative inhibition and putative excitation

Putative inhibition attenuates with tone frequency and putative excitation exhibits a “best frequency”

Putative inhibition precedes putative excitation

Consideration of alternate input configurations

Discussion

A minimal model of the MSO reproduces diverse and frequency-dependent neurophonic responses

A testable working hypothesis: inhibition explains non-dipole-like neurophonic responses

Evaluation of modeling approach

Field potentials as a route for further study of MSO

Footnotes

References

In this issue

Citation Manager Formats

Jump to section

Keywords

Responses to this article

Jump to comment:

Related Articles

Cited By...

More in this TOC Section

Research Articles

Systems/Circuits

Spatiotemporal V_e patterns.

Approximation of V_e spatiotemporal patterns in the Fourier domain.

“Dipole” and “non-dipole” features in V_e responses

Simulated V_e patterns reproduce diverse and frequency-dependent neurophonic responses

Onset responses identify I₁ and I₂ as source and sink current, respectively