Abstract
Hidden hearing loss (HHL) is a deficit in auditory perception and speech intelligibility that occurs despite normal audiometric thresholds and results from noise exposure, aging, or myelin defects. While mechanisms causing perceptual deficits in HHL patients are still unknown, results from animal models indicate a role for peripheral auditory neuropathies in HHL. In humans, sound localization is particularly important for comprehending speech, especially in noisy environments, and its disruption may contribute to HHL. In this study, we hypothesized that neuropathies of cochlear spiral ganglion neurons (SGNs) that are observed in animal models of HHL disrupt the activity of neurons in the medial superior olive (MSO), a nucleus in the brainstem responsible for locating low-frequency sound in the horizontal plane using binaural temporal cues, leading to sound localization deficits. To test our hypothesis, we constructed a network model of the auditory processing system that simulates peripheral responses to sound stimuli and propagation of responses via SGNs to cochlear nuclei and MSO populations. To simulate peripheral auditory neuropathies, we used a previously developed biophysical SGN model with myelin defects at SGN heminodes (myelinopathy) and with loss of inner hair cell-SGN synapses (synaptopathy). Model results indicate that myelinopathy and synaptopathy in SGNs give rise to decreased interaural time difference (ITD) sensitivity of MSO cells, suggesting a possible mechanism for perceptual deficits in HHL patients. This model may be useful to understand downstream impacts of SGN-mediated disruptions on auditory processing and to eventually discover possible treatments for various mechanisms of HHL.
Introduction
The ability to determine the location of the source of a sound is critical for all animals. They can more easily find prey, escape from predators and survive other dangers in nature thanks to their sound localization skills. Humans, as well, benefit from this ability to assess their safety and to distinguish speech when competing sounds are present. Unlike visual and somatosensory systems, the auditory system does not map the spatial origin of stimuli onto the sensory epithelium. Instead, the brain uses temporal, spectral, and intensity cues to determine the location of the source of sounds in three-dimensional space (). Locating a sound in the horizontal plane requires precise temporal and intensity information coming from both ears. Integration of binaural information in humans takes place in the superior olivary complex (SOC) located in the brainstem, specifically in the medial superior olive (MSO), where azimuthal sound localization occurs. MSO cells receive binaural excitatory inputs from spherical bushy cells (SBCs) and binaural inhibitory inputs driven by globular bushy cells (GBCs), which act via relay nuclei, the medial and lateral nuclei of the trapezoid body (MNTB and LNTB) (). SBCs and GBCs are located in the cochlear nucleus, the first relay point for signals from the periphery to the central auditory system. Multiple spiral ganglion neurons (SGNs) project from the cochlea to SBCs and GBCs and anatomical studies show that SBCs typically receive input from 2 to 4 SGNs while GBCs receive input from 9 to 69 SGNs (; ; ; ). Integration of input from multiple SGNs, along with specializations in synaptic and intrinsic physiology, enable SBCs and GBCs to respond to sound with more precise phase-locking than the SGN fibers, therefore transmitting precise timing information to MSO cells (). A reduced representation of this circuitry is shown in Figure 1.
FIGURE 1
Humans can resolve interaural time differences (ITDs), the difference in the arrival time of a sound to each ear, as short as 10 μs, and can locate sound sources as precisely as a few degrees of azimuth (
The precise timing of binaural signals is essential to detect the horizontal direction of the sound source. Therefore, the disruption of signaling along peripheral auditory circuits would significantly impair sound localization ability in humans. Many behavioral and electrophysiological studies suggest that humans with normal audiometric thresholds can have problems with encoding and processing binaural cues, giving rise to speech intelligibility and ITD sensitivity deficits. These deficits occur as a result of noise exposure (
Results
In this study, we aimed to bridge the gap between observed peripheral auditory deficits in animal models of HHL and mechanisms underlying perceptual deficits in patients with HHL. For that purpose, we modeled a mammalian SOC circuit that includes cochlear sound processing and SGN, SBC, GBC, and MSO cell populations (Figure 1), focusing on the crucial role that MSO cells play in sound localization. We simulated circuit responses to binaural sound stimuli at different azimuthal locations under multiple peripheral auditory deficit conditions that have been shown to cause HHL in animal models, i.e., myelinopathy at the SGN axons (
We prescribed SGN myelinopathy levels (see also
Experimental results show that the synaptopathy driven by noise-exposure selectively targets IHC-HT SGN synapses (
Effect of Peripheral Neuropathies on Spike Activity and Dynamics in the Superior Olivary Complex Circuit
First, we explored changes in firing rates of different cell populations as a function of varying myelinopathy and synaptopathy levels (Figure 2). The spike rates of SGN fibers relative to the control (0% Lu variation or 0% synaptopathy) were significantly decreased with increasing degrees of myelinopathy (Figure 2A) and both HT and random synaptopathies (Figures 2B,C). As predicted, bigger drops in activity were observed in random (Figure 2C) as compared to HT synaptopathy (Figure 2B).
FIGURE 2

Spike activity and dynamics in all cell populations was disrupted with higher levels of peripheral auditory deficiencies. Relative spike rate (A–C) and relative VS (D–F) measurements of different cell populations (SGN, SBC, GBC, and MSO) in response to 200 Hz, 50 dB sound stimulus for different peripheral neuropathy conditions: myelinopathy (A,D), synaptopathy at HT (B,E) and random synaptopathy (C,F). In myelinopathy, 0% variation of Lu represents a circuit with a homogeneous population of SGNs with 10 μm long Lu, and 100% variation of Lu represents a circuit with a heterogeneous SGN population with 10 μm ≤ Lu ≤ 20 μm. In both synaptopathy scenarios (HT and random), 0% synaptopathy indicates all IHC-SGN synapses are intact. Synaptopathy level of 100% means all IHC-HT SGN synapses are deficient in HT synaptopathy, whereas the same number of synapses are randomly removed in random synaptopathy (1/3rd of all synapses). In all conditions, increasing the degree of the deficiency (myelinopathy or synaptopathy) decreased relative spike rates and relative vector strength (VS, see “Materials and Methods” section) of all cell types, yet with different slopes for different scenarios. However, in all scenarios, decreases were more pronounced for MSO cells.
This outcome is in agreement with both experimental studies on animal models (
Furthermore, comparing population firing rates in both synapse loss conditions (HT and random synaptopathies) suggests that random synaptopathy has a larger impact on the activities of cochlear nucleus cells, whereas HT synaptopathy barely decreases spike rates of SBCs and GBCs (Figures 2B,C). This difference arises from low activity levels of HT SGN fibers at 50 dB (
Next, we investigated the effect of the SGN neuropathies on phase locking of neuronal firing to the sound wave. The relative vector strength (VS, see “Materials and Methods”) of SGN axons was approximately 0.82 in our control case in response to a 200 Hz, 50 dB sound pulse, while SBCs had a relative VS of 0.93 (Figures 2D–F), indicating increased synchrony compared to the SGN input. The GBC and MSO cell responses were also highly phase-locked to the sound wave, with relative VS approaching to 1.0 (Figures 2D–F). These results agree with experimental observations (
In the myelinopathy scenario, the decrease in MSO activity arises from two myelinopathy outcomes at the SGN level: lower and increasingly asynchronous SGN activity. To separately demonstrate the significance of synchronous SGN input on MSO activity, we simulated the activity of all cell types with artificially randomized SGN input. The raster plot of the putative control SGN population in response to 200 Hz, 50 dB sound stimulus is shown in Figure 3A. We then artificially randomized SGN spike times for all myelinopathy levels in order to disrupt synchrony. We jittered each SGN spike time by an amount δ and analyzed two levels of randomization: −1.25 ms ≤ δ ≤ 1.25 ms (i.e., low jitter) and −2.5 ms ≤ δ ≤ 2.5 ms (i.e., high jitter, see Figure 3B for the raster plot of a randomized putative control). Even though SGN spike rates were unchanged (Figure 3C), the activities of downstream cells [SBCs (Figure 3D), GBCs (Figure 3E) and MSOs (Figure 3F)] decreased with increasing jitter, i.e., higher values of δ, for all Lu variations. Moreover, MSO cell activity was most affected, as SGN inputs with low jitter (dotted-blue line in Figure 3F) and high jitter (dashed-red line in Figure 3F) resulted in no spiking activity in the MSO population for any Lu variation. Relative vector strengths in all populations decreased with low jitter compared to the control case for all Lu variations (Figure 3H) and phase-locking to sound was essentially completely abolished in all populations with high jitter regardless of Lu variation (Figure 3I). The profound effect on MSO spiking provides evidence that not only the level of SGN input but also the degree of its synchronization determines MSO activity.
FIGURE 3

Synchronized synaptic input from SGN cells was necessary to maintain activity rates in downstream cell populations. (A) Raster plot showing the activity of a normal SGN population in response to 200 Hz, 50 dB sound stimulus. (B) Raster plot where the spike times shown in panel (A) are randomized by adding a jitter δ varying between –2.5 and 2.5 ms (high jitter case) to the SGN spike times in panel (A). (C) SGN, (D) SBC, (E) GBC, and (F) MSO spike rates for low (–1.25 ms ≤ δ ≤ 1.25 ms) and high jitter and Lu variation levels. Here, 0% variation of Lu represents a circuit with a homogeneous population of SGNs with 10 μm long Lu and 100% variation of Lu represents a circuit with a heterogeneous SGN population with 10 μm ≤ Lu ≤ 20 μm. (G–I) Relative VS of all cell populations at all myelinopathy levels with (G) normal (same figure as Figure 2D), (H) low jitter and (I) high jitter SGN inputs. Increasing randomization of SGN spikes decreased spike rates and relative VS of all cell types at all myelinopathy levels. Note that there is no cyan line in panels (H,I) due to no MSO cell firing in low jitter and high jitter cases (dotted-blue and dashed-red lines in panel (F).
Effects of Peripheral Neuropathies on Interaural Time Difference Coding in the Medial Superior Olive
Since the horizontal locations of sound sources are encoded by the difference in firing rates between ipsi- and contralateral MSO populations (
FIGURE 4

Disruption in the peripheral auditory system decreased MSO spike rates for all ITDs, resulting in smaller differences between left and right MSO activity and hence lower ITD sensitivity. (A–C) Spike rates of left (solid lines) and right (dashed lines) MSO cells as a function of ITD for various levels of (A) myelinopathy at SGN cells, (B) synaptopathy at IHC-HT SGN synapses and (C) random synaptopathy in IHC-SGN synapse population. (D–F) The difference between left and right MSO firing rates (MSOleft–MSOright) at various levels of (D) myelinopathy, (E) synaptopathy at HT synapses and (F) random synaptopathy. In myelinopathy, 0% variation of Lu represents a circuit with a homogeneous population of SGNs with 10 μm long Lu, and 100% variation of Lu represents a circuit with a heterogeneous SGN population with 10 μm ≤ Lu ≤ 20 μm. In both synaptopathy scenarios (HT and random), 0% synaptopathy indicates all IHC-SGN synapses are intact. In HT-Synaptopathy a level of 100% means that all IHC-HT SGN synapses are deficient, whereas the same number of synapses are randomly removed in random synaptopathy (1/3rd of all synapses).
Resonance of Medial Superior Olive Cells Compensates for Effects of Peripheral Neuropathies
The asymmetric bell-shaped curves in Figures 4A–C are hallmarks of MSO cell activity resulting from their ability to detect coincident subthreshold presynaptic signals with a high temporal precision (
FIGURE 5

The effects of different peripheral auditory disruption conditions on the activity of MSO cells were less pronounced at the sound frequencies where MSO cells show resonant responses. (A) Spike probability (color bar) of individual MSO cells with varying synaptic conductances (Y-axis) in response to different sound frequencies (X-axis). Results provide evidence of resonant behavior for sound frequencies of ∼ 300 Hz. (B) The amplitudes in panels (C–K) are the distance between the peak and the trough of ΔR, which is the difference between the activities of left and right MSO cells for different ITDs. (C–K) The amplitudes of ΔR in response to 50 dB (C–E), 65 dB (F–H), and 80 dB (I–K) sound stimuli at frequencies varying from 200 to 1,200 Hz in peripheral auditory disruption scenarios of various levels. In myelinopathy (C,F,I), 0% variation of Lu represents a circuit with a homogeneous population of SGNs with 10 μm long Lu and 100% variation of Lu represents a circuit with a heterogeneous SGN population with 10 μm ≤ Lu ≤ 20 μm. In both synaptopathy scenarios [HT (D,G,J) and random (E,H,K)], 0% synaptopathy indicates all IHC-SGN synapses are intact. Synaptopathy level of 100% means all IHC-HT SGN synapses are deficient in HT synaptopathy, whereas the same number of synapses are randomly removed in random synaptopathy (1/3rd of all synapses).
Next, to determine whether these deficits have differential outcomes on ITD sensitivities of MSOs at resonant compared to non-resonant sound frequencies, we simulated the SOC circuit model with peripheral auditory deficits in response to varying frequencies of sound stimuli. We quantified MSO cells’ ITD sensitivities by the amplitude of the difference between left and right MSO firing rates, ΔR, as a function of ITDs (Figure 5B) and their activity measured by their spike rates (Figure 6). The amplitudes (Figures 5C–K) and the spike rates (Figure 6) of our putative control (0% Lu variation or 0% synaptopathy) were much higher in response to sound around the resonant frequency (∼300 Hz). However, the sound frequency that resulted in this peak mean activity/ITD sensitivity decreased with increasing sound level from ∼400 Hz at 50 dB to ∼200 Hz at 80 dB. This outcome may stem from higher neurotransmitter release rates from IHCs in response to higher sound levels. Myelinopathy and both synaptopathy scenarios had similar effects on ITD sensitivities (Figures 5C–K) and MSO activities (Figure 6): both measures decreased for all peripheral auditory deficit scenarios, showing graded decreases with increased neuropathy levels for sound stimuli at all frequencies. However, the decreases in ITD sensitivity and MSO spike rate were not as severe for sound stimuli at the resonant frequencies compared to other frequencies.
FIGURE 6

Resonance properties of MSO cells resulted in higher MSO spike rates at resonant frequencies for all peripheral auditory disruption conditions. (A–I) MSO spike rates in response to 50 dB (A–C), 65 dB (D–F), and 80 dB (G–I) sound stimuli of frequencies varying from 200 to 1,200 Hz in case of myelinopathy (A,D,G), synaptopathy at HT (B,E,H) and random synaptopathy (C,F,I). Note that the resonant frequency (frequency with the highest MSO activity) decreased with increasing sound levels.
Effects of Globular Bushy Cells-Mediated Inhibition on Interaural Time Difference Coding in the Medial Superior Olive
Our SOC model includes GBC-mediated inhibition to MSO cells, however, competing theories exist as to the role of this inhibition for ITD coding in the MSO (
First, we simulated MSO activity of our putative control case in response to various sound stimuli and compared the properties of ΔR functions of MSO cells with and without GBC-mediated inhibition. To quantify differences in responses, we used two measures: best ITD difference and amplitude difference (Supplementary Figure 1A). We defined best ITD difference as the difference in MSO cells’ best ITDs between the models with and without inhibition [(best ITD)without inhibition − (best ITD)with inhibition], where best ITD is the ITD value at which MSO cells exhibit the highest activity, i.e., ΔR has the highest value. Results showed that eliminating inhibition does not shift the best ITD significantly (Supplementary Figure 1B). Furthermore, differences of ΔR amplitudes (see Figure 5B for the definition of amplitude) between the two models [(Amplitude)without inhibition − (Amplitude)with inhibition] demonstrated that GBC-mediated inhibition only slightly modulated MSO activity, with bigger changes in response to sound stimuli closer to the resonant frequencies, ∼300 Hz (Supplementary Figure 1C). Finally, we simulated myelin defects in the no-inhibition model and observed the same effects of myelinopathy as in the model with GBC-mediated inhibition, specifically MSO activity was decreased with higher levels of Lu variation and this decrease was less-pronounced for sound stimuli with near-resonant frequencies (Supplementary Figure 2).
Discussion
In this study, we built a computational model of mammalian brainstem auditory circuits to understand the impact of peripheral neural deficits, such as SGN myelin defects or loss of IHC-SGN synapses, on binaural auditory processing. Specifically, we explored how the activity of SGNs, cochlear nucleus cells (SBCs and GBCs) and MSO cells are affected by varying degrees of demyelination and synaptopathy of auditory nerve and IHC-SGN synapses, respectively. Motivated by experimental results in animal models, we modeled the degree of demyelination by increasing the range of Lu, the length of the initial unmyelinated segment in SGN axons (
Since the difference in MSO firing rates between the left and right MSO is critical for the detection of the horizontal angle of sound sources, lower activity levels in both MSOs, as we observed in all peripheral auditory deficit scenarios we studied, decreases this difference (Figure 4), presumably causing binaural processing and sound localization deficits. These results are in line with behavioral human studies, which provided evidence that myelin defects affecting the auditory nerve or auditory brainstem generate problems in locating sound (
As MSO cells are known to have resonance properties due to their phasic firing patterns (
In addition, a counterintuitive effect arises in response to higher frequency sounds at higher levels (80 dB SPL, 800 and 1,000 Hz sound stimuli) where the amplitude of ΔR seemed to increase for increasing levels of neuropathies (Figures 5I–K), even though the MSO spike rates did not show the same effect (Figures 6G–I). This counterintuitive outcome may stem from the fact that MSO cells were overly stimulated at higher sound levels, resulting in smaller differences between left and right MSO spike rates as a function of varying ITDs. However, introducing peripheral neuropathies resulted in a decrease in the amount of input MSO cells received, causing a more significant difference between both MSO spike rates for varying ITDs, even though the absolute rate of MSO activity decreased for all ITD levels.
The compensation for peripheral auditory deficits by MSO cells’ resonance properties would only be significant if these deficits occurred at the SGN fibers having CFs corresponding to resonant frequencies of MSOs. There is no study to our knowledge that investigates the CFs of SGN fibers with myelin defects. However, several animal studies provided evidence that hidden hearing loss occurs as a result of synaptopathy at IHC’s with high CF [>10 kHz in mice (
Changes in the phases of firing, relative to the sound wave, of cells upstream of MSO may also contribute to auditory processing deficits. Model results supported this effect as shown by the significant degradation of MSO activity when SGN spike times were jittered (Figure 3). While those results analyzed variations in the phase-locking to the sound wave, as measured by the VS, the possibility remains that the actual phase angle of spikes relative to the sound wave may be disrupted by SGN myelinopathy and synaptopathy. To investigate this, we further computed average phase angles of spikes, relative to sound waves, of all cell types in our SOC circuit model in response to 50 dB SPL sound stimuli with frequencies of 200 Hz (Supplementary Figure 3), 400 Hz (Supplementary Figure 4), and 600 Hz (Supplementary Figure 5) for both myelinopathy (Supplementary Figures 3A–5A) and HT synaptopathy (Supplementary Figures 3B–5B) scenarios. In the SGN, SBC, and GBC populations, we didn’t observe any significant changes in average phase angles for either neuropathy scenario in response to any frequency sound stimuli. Interestingly, the phase of firing of MSO cells increased significantly for higher myelinopathy levels in response to 400 Hz 50 dB sound stimuli, the frequency closest to the resonant frequency (Supplementary Figure 4A). This result may arise from the resonance properties of MSO cells at 400 Hz, where the absolute decrease in spike rates is the highest for higher levels of myelinopathy (Figure 6A), significantly affecting the average firing phases of MSO cells. In myelinopathy, MSO cells only fire in the first few sound cycles (Supplementary Figure 4M) and therefore lose the subsequent spikes that fire at earlier phases relative to sound due to the gradual phase shift over time that is apparent in the putative control (Supplementary Figure 4L). The loss of only low-phase spikes results in an increase in the average phase angle of firing in MSO cells. Further investigation on circuits downstream of MSO cells (e.g., inferior colliculus, medial geniculate body, etc.) would be useful to better understand the effect of myelinopathy-induced phase angle increases on auditory perceptual deficits. This analysis showing minimal effects of peripheral neuropathy on spike phase angle, relative to the sound wave, suggests that the primary mechanism underlying binaural deficits obtained in the model is the loss of excitatory signaling from SGN fibers, that results in lower coincident signals downstream of SGNs, thus decreasing activity levels in cochlear nucleus and MSO cells.
In our SOC circuit, we included binaural inhibitory inputs that MSO cells receive from GBCs via relays through the medial and lateral nuclei of the trapezoid body (MNTB and LNTB). However, the function of these inhibitory signals has remained controversial over the years. In particular, some studies suggest that they play a role in shifting the best ITD of MSO cells toward more contralateral ITD values, contributing to the ITD sensitivity of MSO cells (
In conclusion, our model results showed that mechanisms underlying HHL in animal models, such as myelin defects at SGN fibers or synapse loss at IHC-SGN synapses, may have significant effects on downstream auditory cell responses, such as MSO cell activity, which plays a crucial role in sound localization. Results indicate that loss of SGN spiking activity as well as potentially asynchronous spike timing contribute to a degradation of ITD sensitivity and coding in the MSO. Model results predict that the primary mechanism inducing binaural auditory processing deficits in myelinopathy and synaptopathy conditions is a decrease in SGN firing activity.
These results may provide a reasonable explanation for human auditory deficits where audiometric thresholds are normal but encoding and processing of binaural cues are degraded. In addition, this study may give us insight into possible treatments for HHL scenarios to overcome these binaural processing deficits. To better elucidate the mechanisms of perceptual deficits resulting from peripheral auditory neuropathies, effects on LSO activity patterns should also be investigated, as LSO cells are responsible for localizing high frequency sounds (
Materials and Methods
Peripheral Auditory System Model
The activity of the peripheral auditory system is simulated by a guinea pig model described in
Spiral Ganglion Neurons Fiber Model
A compartmental model for peripheral axons of SGN fibers is modeled as described in
Each IHC-SGN synapse is connected to one SGN fiber. Each release event determined by the peripheral auditory system model triggers a post-synaptic response at the corresponding SGN fiber that is modeled as an external current pulse (
We include three different types of IHC-SGN synapses based on the response to varying sound stimuli levels: low-threshold (LT), medium-threshold (MT), and high-threshold (HT) synapses. We denote the SGNs based on the type of synapses they are connected to (e.g., an SGN fiber connected to a LT synapse is called a LT SGN). In each side, we model 100 low-threshold (LT), 100 medium-threshold (MT), and 100 high-threshold (HT) SGNs.
Cochlear Nucleus Network Structure
Spiral ganglion neurons that are activated by sound stimuli relay this signal to the cochlear nucleus. SBCs and GBCs in the cochlear nucleus receive excitatory inputs from multiple ipsilateral SGNs [2–4 to SBCs and 9–69 to GBCs (
In our reduced cochlear nucleus circuit model, we modeled SBC, GBC, and MSO cell populations, with each population containing 300 neurons (Figure 1). SBC neurons send direct excitatory input to MSO cells and, for simplicity, we assumed GBCs send inhibitory signals directly to MSO cells, as in
Node Dynamics of Spherical Bushy Cells, Globular Bushy Cells, and Medial Superior Olive
For neuron dynamics, we implemented previously developed models used in a cochlear nucleus circuit model of gerbils (
The membrane potentials of SBCs and GBCs are modeled as:
where Cm is the membrane capacitance, gl is the leak conductance, gNa is the Na+ conductance, gKHT and gKLT are high and low threshold K+ conductances, respectively, and gh is the conductance of the hyperpolarization-activated cation current, or H current. Erest stands for the resting membrane potential and Ex represents the Nernst potentials of each ionic current x (for x = Na+, K+, and H) (Table 1). The variables m, h, n, p, w, z, and r are the voltage dependent conductance gating variables expressed as:
and
where,
TABLE 1
| Parameters | SBC/GBC | MSO | |
| Cm (pF) | Membrane capacitance | 12 | 70 |
| Erest (mV) | Resting membrane potential | −65 | −55.8 |
| gl (nS) | Leak conductance | 37 | 13 |
| ENa (mV) | Nernst potential of Na+ | 50 | 56.2 |
| gNa (nS) | Na+ conductance | 4592.8 | 3900 |
| EK (mV) | Nernst potential of K+ | −77 | −90 |
| gKHT (nS) | High threshold K+ conductance | 35.1 | N/A |
| gKLT (nS) | Low threshold K+ conductance | 367.4 | 650 |
| Eh (mV) | Nernst potential of H current | −43 | −35 |
| gh (nS) | H current conductance | 36.7 | 520 |
| A (nS) | Synaptic strength | 13/4.76 | Isyn,e: 54.37/Isyn,i: 5.5 |
Parameters for spherical bushy cells (SBC), globular bushy cells (GBC), and medial superior olive (MSO) cells.
Isyn is the excitatory synaptic current generated by SGN activity expressed as
where N is the number of presynaptic spikes, tsi is the time of the presynaptic SGN spike i and Eex = 0 mV is the reversal potential for excitatory current. A is the synaptic strength and equals 4.76 nS for GBCs, as in
The current balance equation for MSO cells is expressed as:
with parameters as in Equation 1 (Table 1). The variables m, h, w, z, and r are the voltage dependent conductance gating variables that are governed by Equation 3 with steady state activation and time constant functions given by:
Isyn,e and Isyn,i are excitatory and inhibitory synaptic currents received from the SBC and GBC cells, respectively. Isyn,e is described as
where N is the number of presynaptic SBC spikes, tsi is the time of the i-th presynaptic SBC spike, Eex = 0 mV is the reversal potential for excitatory current and A = 54.37 nS is the synaptic strength. The tdelay is the time required for the signal from SBCs to reach MSO cells, which is 1.5 ms for ipsilateral input and 1.6 ms for contralateral input (
Isyn,i is expressed as:
where N is the number of presynaptic GBC spikes, tsi is the time of the i-th presynaptic GBC spike, Eex = −70 mV is the reversal potential for excitatory current and A = 5.5 nS is the synaptic strength. The tdelay is the time required for the signal from inhibitory GBCs to reach MSO cells, which is 1.5 ms for ipsilateral input and 1.0 ms for contralateral input (
Relative Vector Strength Measurement
Vector strength (VS) is a measure to determine the degree of phase-locking of spiking in a neuron population to a sound wave. To calculate VS for each neural population, the phase angle θi relative to the sound wave of each spike i fired by cells in the population is first measured as:
where n is the number of spikes of a cell population, ti is the time of each spike i, and t1 and t2 are the peaks of the sound wave just before and just after ti, respectively. Vector strength (
This measure varies between 0 and 1, where 1 means all spikes are phase-locked to the sound wave at the same angle and 0 means no phase-locking to sound.
A measure of relative VS is derived from VS in order to assess the degree of phase-locking relative to the putative control. It is defined as:
where N is the total number of spikes in the putative control. As in VS, relative VS also varies between 0 and 1. Here, a relative VS = 1 means perfect phase-locking to sound with the same number of spikes relative to the putative control, whereas a relative VS = 0 happens in case of low phase-locking and/or low number of spikes relative to the control (Figure 7).
FIGURE 7

Examples of vector strength (VS) and relative VS measurements for various groups of neurons in response to sound stimulus. In the control case, the spike times of all neurons are phase-locked to the sound wave, i.e., they always fire at the same phase of the sound wave and VS approaches 1. In Case 1, phase-locking is low, resulting in a VS ∼ 0. As the number of spikes is the same as the control in Case 1, VS and relative VS are the same. However, the group of neurons in Case 2 does not fire nearly as much as the control case, leading to low relative VS, even though its VS is high.
Identifying Resonance Properties
To detect resonance properties in the MSO cell model, we simulated MSO cell responses to direct excitatory synaptic input arriving at various frequencies. We used the same node dynamics of MSO cells described before (see Equations 19–31 and Table 1) and only modified excitatory and inhibitory synaptic input (Isyn,e and Isyn,i, respectively, see Equations 30, 31), such that:
and
where
Here, ts is the presynaptic spike time, Eex = 0 mV is the reversal potential for excitatory current, f is the frequency of presynaptic spikes, Vm is the membrane potential and A is the synaptic strength (synaptic conductance). For the resonance study, we varied f and A, and calculated spike probabilities (Figure 5) as follows:
Simulations
We simulated the sound-evoked activity of all cell types in response to a sound stimulus pulse of 200 Hz and 50 dB for 100 ms, unless stated otherwise. Our results show the responses averaged over five sound stimuli.
Conclusion
Some people have difficulty understanding speech in crowded social settings, also known as the “cocktail party problem,” despite having normal hearing thresholds for all sound frequencies. One potential cause of this condition is “hidden hearing loss” (HHL), an auditory disorder resulting from noise exposure, aging or peripheral neuropathy. In this study, we hypothesized that the perceptual deficits caused by HHL arise from sound localization problems due to disrupted inputs to the medial superior olive (MSO), a nucleus in the brainstem that integrates binaural signals to determine the relative timing of sound arrival to both ears and to detect the horizontal angle of the sound source. To explore the impacts HHL has on MSO activity, we simulated MSO circuits that receive signals from both ears affected by two peripheral neuropathies which have been previously shown to cause HHL in animal models: (1) loss of synapses between inner hair cells and auditory nerve fibers, and (2) disruption of auditory-nerve myelin. We provide evidence that both scenarios disrupt the activity of MSO cells correlated with sound localization that may, in turn, result in speech intelligibility deficits. This model may be used to elucidate downstream effects of peripheral neuropathies and to propose possible clinical treatments for HHL.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
MB performed the research and analyzed the data. MB, MR, KG, GC, VB, and MZ designed the research, wrote the manuscript, contributed to the article, and approved the submitted version.
Funding
This work was supported by National Institutes of Health-National Institute on Deafness and Other Communication Disorders: NIDCD RO1DC018500 (GC) and R01DC018284 (MR), R01 DC04084 (KG); National Institute of Biomedical Imaging and Bioengineering: NIBIB R01EB018297 (MZ and VB).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fncir.2022.856926/full#supplementary-material
References
1
BudakM.GroshK.SasmalA.CorfasG.ZochowskiM.BoothV. (2021). Contrasting mechanisms for hidden hearing loss: Synaptopathy vs myelin defects.PLoS Comput. Biol.17:e1008499. 10.1371/journal.pcbi.1008499
2
BernsteinL. R.TrahiotisC. (2016). Behavioral manifestations of audiometrically-defined “slight” or “hidden” hearing loss revealed by measures of binaural detection.J. Acoust. Soc. Am.1403540. 10.1121/1.4966113
3
BharadwajH. M.MasudS.MehraeiG.VerhulstS.Shinn-CunninghamB. G. (2015). Individual differences reveal correlates of hidden hearing deficits.J. Neurosci.352161–2172. 10.1523/JNEUROSCI.3915-14.2015
4
BishopH. I.ZitoK. (2013). The downs and ups of sensory deprivation: evidence for firing rate homeostasis in vivo.Neuron80247–249. 10.1016/j.neuron.2013.10.011
5
ChoiJ. E.SeokJ. M.AhnJ.JiY. S.LeeK. M.HongS. H.et al (2018). Hidden hearing loss in patients with Charcot-Marie-Tooth disease type 1A.Sci. Rep.810335. 10.1038/s41598-018-28501-y
6
EddinsA. C.OzmeralE. J.EddinsD. A. (2018). How aging impacts the encoding of binaural cues and the perception of auditory space.Hear. Res.36979–89. 10.1016/j.heares.2018.05.001
7
EnckeJ.HemmertW. (2018). Extraction of Inter-Aural Time Differences Using a Spiking Neuron Network Model of the Medial Superior Olive.Front. Neurosci.12:140. 10.3389/fnins.2018.00140
8
FurmanA. C.KujawaS. G.LibermanM. C. (2013). Noise-induced cochlear neuropathy is selective for fibers with low spontaneous rates.J. Neurophysiol.110577–586. 10.1152/jn.00164.2013
9
FurstM.LevineR. A. (2015). Hearing disorders in multiple sclerosis.Handb. Clin. Neurol.129649–665. 10.1016/B978-0-444-62630-1.00036-6
10
GaiY.DoironB.RinzelJ. (2010). Slope-based stochastic resonance: how noise enables phasic neurons to encode slow signals.PLoS Comput. Biol.6:e1000825. 10.1371/journal.pcbi.1000825
11
GleichO.SemmlerP.StrutzJ. (2016). Behavioral auditory thresholds and loss of ribbon synapses at inner hair cells in aged gerbils.Exp. Gerontol.8461–70. 10.1016/j.exger.2016.08.011
12
GoldbergJ. M.BrownP. B. (1969). Response of binaural neurons of dog superior olivary complex to dichotic tonal stimuli: some physiological mechanisms of sound localization.J. Neurophysiol.32613–636. 10.1152/jn.1969.32.4.613
13
GrotheB.NeuweilerG. (2000). The function of the medial superior olive in small mammals: temporal receptive fields in auditory analysis.J. Comp. Physiol. A.186413–423. 10.1007/s003590050441
14
GrotheB.ParkT. J. (1998). Sensitivity to interaural time differences in the medial superior olive of a small mammal, the Mexican free-tailed bat.J. Neurosci.186608–6622.
15
GuinanJ. J.GuinnanS. S.NorrisB. E. (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.
16
GuinanJ. J.NorrisB. E.GuinanS. S. (1972b). Single Auditory Units in the Superior Olivary Complex: II: Locations of Unit Categories and Tonotopic Organization.Int. J. Neurosci.4147–166.
17
HarperN. S.McAlpineD. (2004). Optimal neural population coding of an auditory spatial cue.Nature430682–686. 10.1038/nature02768
18
JorisP. X.CarneyL. H.SmithP. H.YinT. C. (1994). Enhancement of neural synchronization in the anteroventral cochlear nucleus. I. Responses to tones at the characteristic frequency.J. Neurophysiol.711022–1036. 10.1152/jn.1994.71.3.1022
19
KandelE. R.SchwartzJ. H.JessellT. M. (2000). Principles of Neural Science.New York: McGraw-Hill Companies.
20
KuenzelT. (2019). Modulatory influences on time-coding neurons in the ventral cochlear nucleus.Hear Res.384:107824. 10.1016/j.heares.2019.107824
21
KuenzelT.BorstJ. G.van der HeijdenM. (2011). Factors controlling the input-output relationship of spherical bushy cells in the gerbil cochlear nucleus.J. Neurosci.314260–4273. 10.1523/JNEUROSCI.5433-10.2011
22
KujawaS. G.LibermanM. C. (2009). Adding insult to injury: cochlear nerve degeneration after “temporary” noise-induced hearing loss.J. Neurosci.2914077–14085. 10.1523/JNEUROSCI.2845-09.2009
23
LibermanM. C. (1991). Central projections of auditory-nerve fibers of differing spontaneous rate. I. Anteroventral cochlear nucleus.J. Comp. Neurol.313240–258. 10.1002/cne.903130205
24
LorteijeJ. A.RusuS. I.KushmerickC.BorstJ. G. (2009). Reliability and precision of the mouse calyx of Held synapse.J. Neurosci.2913770–13784. 10.1523/JNEUROSCI.3285-09.2009
25
Mc LaughlinM.van der HeijdenM.JorisP. X. (2008). How secure is in vivo synaptic transmission at the calyx of Held?J. Neurosci.2810206–10219. 10.1523/JNEUROSCI.2735-08.2008
26
MengX.HuguetG.RinzelJ. (2012). TYPE III EXCITABILITY, SLOPE SENSITIVITY AND COINCIDENCE DETECTION.Discrete Contin. Dyn. Syst. Ser. A322729–2757. 10.3934/dcds.2012.32.2729
27
MiddlebrooksJ. C. (2015). Sound localization.Handb. Clin. Neurol.12999–116. 10.1016/B978-0-444-62630-1.00006-8
28
Mikiel-HunterJ.KotakV.RinzelJ. (2016). High-Frequency Resonance in the Gerbil Medial Superior Olive.PLoS Comput. Biol.12:e1005166. 10.1371/journal.pcbi.1005166
29
MyogaM. H.LehnertS.LeiboldC.FelmyF.GrotheB. (2014). Glycinergic inhibition tunes coincidence detection in the auditory brainstem.Nat. Commun.5:3790.
30
ParkT. J.KlugA.HolinstatM.GrotheB. (2004). Interaural level difference processing in the lateral superior olive and the inferior colliculus.J. Neurophysiol.92289–301.
31
PrendergastG.MillmanR. E.GuestH.MunroK. J.KlukK.DeweyR. S.et al (2017). Effects of noise exposure on young adults with normal audiograms II: Behavioral measures.Hear. Res.35674–86. 10.1016/j.heares.2017.10.007
32
RanceG.RyanM. M.CarewP.CorbenL. A.YiuE.TanJ.et al (2012). Binaural speech processing in individuals with auditory neuropathy.Neuroscience226227–235. 10.1016/j.neuroscience.2012.08.054
33
RemmeM. W.DonatoR.Mikiel-HunterJ.BallesteroJ. A.FosterS.RinzelJ.et al (2014). Subthreshold resonance properties contribute to the efficient coding of auditory spatial cues.Proc. Natl. Acad. Sci. U.S.A.111E2339–E2348. 10.1073/pnas.1316216111
34
RobertsM. T.SeemanS. C.GoldingN. L. A. (2013). mechanistic understanding of the role of feedforward inhibition in the mammalian sound localization circuitry.Neuron78923–935. 10.1016/j.neuron.2013.04.022
35
RyugoD. K.SentoS. (1991). Synaptic connections of the auditory nerve in cats: relationship between endbulbs of held and spherical bushy cells.J. Comp. Neurol.30535–48. 10.1002/cne.903050105
36
ScottL. L.MathewsP. J.GoldingN. L. (2005). Posthearing developmental refinement of temporal processing in principal neurons of the medial superior olive.J. Neurosci.257887–7895. 10.1523/JNEUROSCI.1016-05.2005
37
SentoS.RyugoD. K. (1989). Endbulbs of held and spherical bushy cells in cats: morphological correlates with physiological properties.J. Comp. Neurol.280553–562. 10.1002/cne.902800406
38
SergeyenkoY.LallK.LibermanM. C.KujawaS. G. (2013). Age-related cochlear synaptopathy: an early-onset contributor to auditory functional decline.J. Neurosci.3313686–13694. 10.1523/JNEUROSCI.1783-13.2013
39
SpirouG. A.RagerJ.ManisP. B. (2005). Convergence of auditory-nerve fiber projections onto globular bushy cells.Neuroscience136843–863. 10.1016/j.neuroscience.2005.08.068
40
SteadmanM. (2018). Mark-Steadman/Neural-Representations-Of-Speech: Neural Representations Of Speech (Version v1.1).Geneva: Zenodo.
41
SteadmanM. A.SumnerC. J. (2018). Changes in Neuronal Representations of Consonants in the Ascending Auditory System and Their Role in Speech Recognition.Front. Neurosci.12:671. 10.3389/fnins.2018.00671
42
SumnerC. J.Lopez-PovedaE. A.O’MardL. P.MeddisR. A. (2002). revised model of the inner-hair cell and auditory-nerve complex.J. Acoust. Soc. Am.1112178–2188. 10.1121/1.1453451
43
SvirskisG.KotakV.SanesD. H.RinzelJ. (2002). Enhancement of signal-to-noise ratio and phase locking for small inputs by a low-threshold outward current in auditory neurons.J. Neurosci.2211019–11025.
44
TeichertM.LiebmannL.HübnerC. A.BolzJ. (2017). Homeostatic plasticity and synaptic scaling in the adult mouse auditory cortex.Sci. Rep.7:17423. 10.1038/s41598-017-17711-5
45
TsuchitaniC. (1977). Functional organization of lateral cell groups of cat superior olivary complex.J. Neurophysiol.40296–318. 10.1152/jn.1977.40.2.296
46
TurrigianoG. G. (2008). The self-tuning neuron: synaptic scaling of excitatory synapses.Cell135422–435. 10.1016/j.cell.2008.10.008
47
TurrigianoG. G.LeslieK. R.DesaiN. S.RutherfordL. C.NelsonS. B. (1998). Activity-dependent scaling of quantal amplitude in neocortical neurons.Nature391892–896. 10.1038/36103
48
ValeroM. D.BurtonJ. A.HauserS. N.HackettT. A.RamachandranR.LibermanM. C. (2017). Noise-induced cochlear synaptopathy in rhesus monkeys (Macaca mulatta).Hear. Res.353213–223. 10.1016/j.heares.2017.07.003
49
VianaL. M.O’MalleyJ. T.BurgessB. J.JonesD. D.OliveiraC. A.SantosF.et al (2015). Cochlear neuropathy in human presbycusis: Confocal analysis of hidden hearing loss in post-mortem tissue.Hear Res.32778–88. 10.1016/j.heares.2015.04.014
50
WanG.CorfasG. (2017). Transient auditory nerve demyelination as a new mechanism for hidden hearing loss.Nat. Commun.8:14487. 10.1038/ncomms14487
51
WhittJ. L.PetrusE.LeeH. K. (2014). Experience-dependent homeostatic synaptic plasticity in neocortex.Neuropharmacology7845–54.
52
WinterI. M.RobertsonD.YatesG. K. (1990). Diversity of characteristic frequency rate-intensity functions in guinea pig auditory nerve fibres.Hear Res.45191–202.
53
YangS.WeinerB. D.ZhangL. S.ChoS. J.BaoS. (2011). Homeostatic plasticity drives tinnitus perception in an animal model.Proc. Natl. Acad. Sci. U.S.A.10814974–14979.
54
YinT. C.ChanJ. C. (1990). Interaural time sensitivity in medial superior olive of cat.J. Neurophysiol.64465–488. 10.1152/jn.1990.64.2.465
Summary
Keywords
hidden hearing loss (HHL), synaptopathy, myelin abnormalities, binaural processing deficits, medial superior olive, computational model
Citation
Budak M, Roberts MT, Grosh K, Corfas G, Booth V and Zochowski M (2022) Binaural Processing Deficits Due to Synaptopathy and Myelin Defects. Front. Neural Circuits 16:856926. doi: 10.3389/fncir.2022.856926
Received
17 January 2022
Accepted
23 March 2022
Published
14 April 2022
Volume
16 - 2022
Edited by
Shaowen Bao, University of Arizona, United States
Reviewed by
Randy J. Kulesza, Lake Erie College of Osteopathic Medicine, United States; Adrian Rodriguez-Contreras, City College of New York (CUNY), United States
Updates

Check for updates
Copyright
© 2022 Budak, Roberts, Grosh, Corfas, Booth and Zochowski.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Gabriel Corfas, corfas@med.umich.eduVictoria Booth, vbooth@umich.eduMichal Zochowski, michalz@umich.edu
†These authors have contributed equally to this work
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.