A biophysical modelling platform of the cochlear nucleus and other auditory circuits: From channels to networks

A Computational Model of Auditory Chirp-Velocity Sensitivity and Amplitude-Modulation Tuning in Inferior Colliculus Neurons

We demonstrate a model of chirp-velocity sensitivity in the inferior colliculus (IC) that retains the tuning to amplitude modulation (AM) that was established in earlier models. The mechanism of velocity sensitivity is sequence detection by octopus cells of the posteroventral cochlear nucleus, which have been proposed in physiological studies to respond preferentially to the order of arrival of cross-frequency inputs of different amplitudes. Model architecture is based on coincidence detection of a combination of excitatory and inhibitory inputs. Chirp-sensitivity of the IC output is largely controlled by the strength and timing of the chirp-sensitive octopus-cell inhibitory input. AM tuning is controlled by inhibition and excitation that are tuned to the same frequency. We present several example neurons that demonstrate the feasibility of the model in simulating realistic chirp-sensitivity and AM tuning for a wide range of characteristic frequencies. Additionally, we explore the systematic impact of varying parameters on model responses. The proposed model can be used to assess the contribution of IC chirp-velocity sensitivity to responses to complex sounds, such as speech.

An Anatomical and Physiological Basis for Coincidence Detection Across Time Scales in the Auditory System

Coincidence detection is a common neural computation that identifies co-occurring stimuli by integration of inputs. In the auditory system, octopus cells act as coincidence detectors for complex sounds that include both synchronous and sequenced combinations of frequencies. Octopus cells must detect coincidence on both the millisecond and submillisecond time scale, unlike the average neuron, which integrates inputs over time on the order of tens of milliseconds. Here, we show that octopus cell computations in the cell body are shaped by inhibition in the dendrites, which adjusts the strength and timing of incoming signals to achieve submillisecond acuity. This mechanism is crucial for the fundamental process of integrating the synchronized frequencies of natural auditory signals over time.

High-resolution volumetric imaging constrains compartmental models to explore synaptic integration and temporal processing by cochlear nucleus globular bushy cells

Globular bushy cells (GBCs) of the cochlear nucleus play central roles in the temporal processing of sound. Despite investigation over many decades, fundamental questions remain about their dendrite structure, afferent innervation, and integration of synaptic inputs. Here, we use volume electron microscopy (EM) of the mouse cochlear nucleus to construct synaptic maps that precisely specify convergence ratios and synaptic weights for auditory nerve innervation and accurate surface areas of all postsynaptic compartments. Detailed biophysically based compartmental models can help develop hypotheses regarding how GBCs integrate inputs to yield their recorded responses to sound. We established a pipeline to export a precise reconstruction of auditory nerve axons and their endbulb terminals together with high-resolution dendrite, soma, and axon reconstructions into biophysically detailed compartmental models that could be activated by a standard cochlear transduction model. With these constraints, the models predict auditory nerve input profiles whereby all endbulbs onto a GBC are subthreshold (coincidence detection mode), or one or two inputs are suprathreshold (mixed mode). The models also predict the relative importance of dendrite geometry, soma size, and axon initial segment length in setting action potential threshold and generating heterogeneity in sound-evoked responses, and thereby propose mechanisms by which GBCs may homeostatically adjust their excitability. Volume EM also reveals new dendritic structures and dendrites that lack innervation. This framework defines a pathway from subcellular morphology to synaptic connectivity, and facilitates investigation into the roles of specific cellular features in sound encoding. We also clarify the need for new experimental measurements to provide missing cellular parameters, and predict responses to sound for further in vivo studies, thereby serving as a template for investigation of other neuron classes.

Mammalian octopus cells are direction selective to frequency sweeps by excitatory synaptic sequence detection

Octopus cells are remarkable projection neurons of the mammalian cochlear nucleus, with extremely fast membranes and wide-frequency tuning. They are considered prime examples of coincidence detectors but are poorly characterized in vivo. We discover that octopus cells are selective to frequency sweep direction, a feature that is absent in their auditory nerve inputs. In vivo intracellular recordings reveal that direction selectivity does not derive from across-frequency coincidence detection but hinges on the amplitudes and activation sequence of auditory nerve inputs tuned to clusters of hot spot frequencies. A simple biophysical octopus cell model excited with real nerve spike trains recreates direction selectivity through interaction of intrinsic membrane conductances with the activation sequence of clustered excitatory inputs. We conclude that octopus cells are sequence detectors, sensitive to temporal patterns across cochlear frequency channels. The detection of sequences rather than coincidences is a much simpler but powerful operation to extract temporal information.

Dynamic Heterogeneity Shapes Patterns of Spiral Ganglion Activity

Neural response properties that typify primary sensory afferents are critical to fully appreciate because they establish and, ultimately represent, the fundamental coding design used for higher-level processing. Studies illuminating the center-surround receptive fields of retinal ganglion cells, for example, were ground-breaking because they determined the foundation of visual form detection. For the auditory system, a basic organizing principle of the spiral ganglion afferents is their extensive electrophysiological heterogeneity establishing diverse intrinsic firing properties in neurons throughout the spiral ganglion. Moreover, these neurons display an impressively large array of neurotransmitter receptor types that are responsive to efferent feedback. Thus, electrophysiological diversity and its neuromodulation are a fundamental encoding mechanism contributed by the primary afferents in the auditory system. To place these features into context, we evaluated the effects of hyperpolarization and cAMP on threshold level as indicators of overall afferent responsiveness in CBA/CaJ mice of either sex. Hyperpolarization modified threshold gradients such that distinct voltage protocols could shift the relationship between sensitivity and stimulus input to reshape resolution. This resulted in an "accordion effect" that appeared to stretch, compress, or maintain responsivity across the gradient of afferent thresholds. cAMP targeted threshold and kinetic shifts to rapidly adapting neurons, thus revealing multiple cochleotopic properties that could potentially be independently regulated. These examples of dynamic heterogeneity in primary auditory afferents not only have the capacity to shift the range, sensitivity, and resolution, but to do so in a coordinated manner that appears to orchestrate changes with a seemingly unlimited repertoire.SIGNIFICANCE STATEMENT How do we discriminate the more nuanced qualities of the sound around us? Beyond the basics of pitch and loudness, aspects, such as pattern, distance, velocity, and location, are all attributes that must be used to encode acoustic sensations effectively. While higher-level processing is required for perception, it would not be unexpected if the primary auditory afferents optimized receptor input to expedite neural encoding. The findings reported herein are consistent with this design. Neuromodulation compressed, expanded, shifted, or realigned intrinsic electrophysiological heterogeneity to alter neuronal responses selectively and dynamically. This suggests that diverse spiral ganglion phenotypes provide a rich substrate to support an almost limitless array of coding strategies within the first neural element of the auditory pathway.

In-vitro Recordings of Neural Magnetic Activity From the Auditory Brainstem Using Color Centers in Diamond: A Simulation Study

Magnetometry based on nitrogen-vacancy (NV) centers in diamond is a novel technique capable of measuring magnetic fields with high sensitivity and high spatial resolution. With the further advancements of these sensors, they may open up novel approaches for the 2D imaging of neural signals in vitro. In the present study, we investigate the feasibility of NV-based imaging by numerically simulating the magnetic signal from the auditory pathway of a rodent brainstem slice (ventral cochlear nucleus, VCN, to the medial trapezoid body, MNTB) as stimulated by both electric and optic stimulation. The resulting signal from these two stimulation methods are evaluated and compared. A realistic pathway model was created based on published data of the neural morphologies and channel dynamics of the globular bushy cells in the VCN and their axonal projections to the principal cells in the MNTB. The pathway dynamics in response to optic and electric stimulation and the emitted magnetic fields were estimated using the cable equation. For simulating the optic stimulation, the light distribution in brain tissue was numerically estimated and used to model the optogenetic neural excitation based on a four state channelrhodopsin-2 (ChR2) model. The corresponding heating was also estimated, using the bio-heat equation and was found to be low (<2°C) even at excessively strong optic signals. A peak magnetic field strength of ∼0.5 and ∼0.1 nT was calculated from the auditory brainstem pathway after electrical and optical stimulation, respectively. By increasing the stimulating light intensity four-fold (far exceeding commonly used intensities) the peak magnetic signal strength only increased to 0.2 nT. Thus, while optogenetic stimulation would be favorable to avoid artefacts in the recordings, electric stimulation achieves higher peak fields. The present simulation study predicts that high-resolution magnetic imaging of the action potentials traveling along the auditory brainstem pathway will only be possible for next generation NV sensors. However, the existing sensors already have sufficient sensitivity to support the magnetic sensing of cumulated neural signals sampled from larger parts of the pathway, which might be a promising intermediate step toward further maturing this novel technology.

Assessing temporal responsiveness of primary stimulated neurons in auditory brainstem and cochlear implant users

The reasons why clinical outcomes with auditory brainstem implants (ABIs) are generally poorer than with cochlear implants (CIs) are still somewhat elusive. Prior work has focused on differences in processing of spectral information due to possibly poorer tonotopic representation and higher channel interaction with ABIs than with CIs. In contrast, this study examines the hypothesis that a potential contributing reason for poor speech perception in ABI users may be the relative lack of temporal responsiveness of the primary neurons that are stimulated by the ABI. The cochlear nucleus, the site of ABI stimulation, consists of different neuron types, most of which have much more complex responses than the auditory nerve neurons stimulated by a CI. Temporal responsiveness of primary stimulated neurons was assessed in a group of ABI and CI users by measuring recovery of electrically evoked compound action potentials (ECAPs) from single-pulse forward masking. Slower ECAP recovery tended to be associated with poorer hearing outcomes in both groups. ABI subjects with the longest recovery time had no speech understanding or even no hearing sensation with their ABI device; speech perception for the one CI outlier with long ECAP recovery time was well below average. To the extent that ECAP recovery measures reveal temporal properties of the primary neurons that receive direct stimulation form neural prosthesis devices, they may provide a physiological underpinning for clinical outcomes of auditory implants. ECAP recovery measures may be used to determine which portions of the cochlear nucleus to stimulate, and possibly allow us to enhance the stimulation paradigms.

Periodicity Pitch Perception

This study presents a computational model to reproduce the biological dynamics of "listening to music." A biologically plausible model of periodicity pitch detection is proposed and simulated. Periodicity pitch is computed across a range of the auditory spectrum. Periodicity pitch is detected from subsets of activated auditory nerve fibers (ANFs). These activate connected model octopus cells, which trigger model neurons detecting onsets and offsets; thence model interval-tuned neurons are innervated at the right interval times; and finally, a set of common interval-detecting neurons indicate pitch. Octopus cells rhythmically spike with the pitch periodicity of the sound. Batteries of interval-tuned neurons stopwatch-like measure the inter-spike intervals of the octopus cells by coding interval durations as first spike latencies (FSLs). The FSL-triggered spikes synchronously coincide through a monolayer spiking neural network at the corresponding receiver pitch neurons.

Neuronal population model of globular bushy cells covering unit-to-unit variability

Computations of acoustic information along the central auditory pathways start in the cochlear nucleus. Bushy cells in the anteroventral cochlear nucleus, which innervate monaural and binaural stations in the superior olivary complex, process and transfer temporal cues relevant for sound localization. These cells are categorized into two groups: spherical and globular bushy cells (SBCs/GBCs). Spontaneous rates of GBCs innervated by multiple auditory nerve (AN) fibers are generally lower than those of SBCs that receive a small number of large AN synapses. In response to low-frequency tonal stimulation, both types of bushy cells show improved phase-locking and entrainment compared to AN fibers. When driven by high-frequency tones, GBCs show primary-like-with-notch or onset-L peristimulus time histograms and relatively irregular spiking. However, previous in vivo physiological studies of bushy cells also found considerable unit-to-unit variability in these response patterns. Here we present a population of models that can simulate the observed variation in GBCs. We used a simple coincidence detection model with an adaptive threshold and systematically varied its six parameters. Out of 567000 parameter combinations tested, 7520 primary-like-with-notch models and 4094 onset-L models were selected that satisfied a set of physiological criteria for a GBC unit. Analyses of the model parameters and output measures revealed that the parameters of the accepted model population are weakly correlated with each other to retain major GBC properties, and that the output spiking patterns of the model are affected by a combination of multiple parameters. Simulations of frequency-dependent temporal properties of the model GBCs showed a reasonable fit to empirical data, supporting the validity of our population modeling. The computational simplicity and efficiency of the model structure makes our approach suitable for future large-scale simulations of binaural information processing that may involve thousands of GBC units.

In the auditory system, specialized neuronal circuits process various types of acoustic information. A group of neurons, called globular bushy cells (GBCs), faithfully transfer timing information of acoustic signals to their downstream neurons responsible for the perception of sound location. Previous physiological studies found representative activity patterns of GBCs, but with substantial individual variations among them. In this study, we present a population of models, instead of creating one best model, to account for the observed variations of GBCs. We varied all six parameters of a simple auditory neuron model and selected the combinations of parameters that led to acceptable activity patterns of GBCs. In total, we tested more than half a million combinations and accepted ~11600 GBC models. Temporal spiking patterns of real GBCs depend on the sound frequency, and our model population was able to replicate this trend. The model used here is computationally efficient and can thus serve as a building block for future large-scale simulations of auditory information processing.

Intrinsic physiological properties underlie auditory response diversity in the avian cochlear nucleus

Sensory systems exploit parallel processing of stimulus features to enable rapid, simultaneous extraction of information. Mechanisms that facilitate this differential extraction of stimulus features can be intrinsic or synaptic in origin. A subdivision of the avian cochlear nucleus, nucleus angularis (NA), extracts sound intensity information from the auditory nerve and contains neurons that exhibit diverse responses to sound and current injection. NA neurons project to multiple regions ascending the auditory brain stem including the superior olivary nucleus, lateral lemniscus, and avian inferior colliculus, with functional implications for inhibitory gain control and sound localization. Here we investigated whether the diversity of auditory response patterns in NA can be accounted for by variation in intrinsic physiological features. Modeled sound-evoked auditory nerve input was applied to NA neurons with dynamic clamp during in vitro whole cell recording at room temperature. Temporal responses to auditory nerve input depended on variation in intrinsic properties, and the low-threshold K⁺ current was implicated as a major contributor to temporal response diversity and neuronal input-output functions. An auditory nerve model of acoustic amplitude modulation produced synchrony coding of modulation frequency that depended on the intrinsic physiology of the individual neuron. In Primary-Like neurons, varying low-threshold K⁺ conductance with dynamic clamp altered temporal modulation tuning bidirectionally. Taken together, these data suggest that intrinsic physiological properties play a key role in shaping auditory response diversity to both simple and more naturalistic auditory stimuli in the avian cochlear nucleus. NEW & NOTEWORTHY This article addresses the question of how the nervous system extracts different information in sounds. Neurons in the cochlear nucleus show diverse responses to acoustic stimuli that may allow for parallel processing of acoustic features. The present studies suggest that diversity in intrinsic physiological features of individual neurons, including levels of a low voltage-activated K⁺ current, play a major role in regulating the diversity of auditory responses.

Brain-like emergent auditory learning: A developmental method

Compared with machine audition, the human auditory system can recognize speech accurately and quickly. This paper proposes a new developmental network (DN) that simulates the human auditory system and constructs an artificial auditory model for speech recognition. The new model simulates each key element of the human auditory pathway as a deep network; in particular, an additional layer in the network is considered to simulate the function of the superior colliculus in the thalamus for speech context integration. The mel-frequency cepstral coefficient (MFCC) is used to extract the features of the speech signal as the sensory input of the DN. The emergent feature of DN model provides an explanation of how such internal neurons represent the short speech context when they are not supervised by the external world. The experimental results show that the recognition rates of English words and phrases can be improved significantly compared to those reported in the existing literature. The proposed DN model provides a new method to solve difficult problems, such as universal speech recognition, in traditional machine audition systems. Meanwhile, the same learning principle can potentially be used in or adapted to other computational contexts and applications.

Dual Coding of Frequency Modulation in the Ventral Cochlear Nucleus

Frequency modulation (FM) is a common acoustic feature of natural sounds and is known to play a role in robust sound source recognition. Auditory neurons show precise stimulus-synchronized discharge patterns that may be used for the representation of low-rate FM. However, it remains unclear whether this representation is based on synchronization to slow temporal envelope (ENV) cues resulting from cochlear filtering or phase locking to faster temporal fine structure (TFS) cues. To investigate the plausibility of those encoding schemes, single units of the ventral cochlear nucleus of guinea pigs of either sex were recorded in response to sine FM tones centered at the unit's best frequency (BF). The results show that, in contrast to high-BF units, for modulation depths within the receptive field, low-BF units (<4 kHz) demonstrate good phase locking to TFS. For modulation depths extending beyond the receptive field, the discharge patterns follow the ENV and fluctuate at the modulation rate. The receptive field proved to be a good predictor of the ENV responses for most primary-like and chopper units. The current in vivo data also reveal a high level of diversity in responses across unit types. TFS cues are mainly conveyed by low-frequency and primary-like units and ENV cues by chopper and onset units. The diversity of responses exhibited by cochlear nucleus neurons provides a neural basis for a dual-coding scheme of FM in the brainstem based on both ENV and TFS cues.SIGNIFICANCE STATEMENT Natural sounds, including speech, convey informative temporal modulations in frequency. Understanding how the auditory system represents those frequency modulations (FM) has important implications as robust sound source recognition depends crucially on the reception of low-rate FM cues. Here, we recorded 115 single-unit responses from the ventral cochlear nucleus in response to FM and provide the first physiological evidence of a dual-coding mechanism of FM via synchronization to temporal envelope cues and phase locking to temporal fine structure cues. We also demonstrate a diversity of neural responses with different coding specializations. These results support the dual-coding scheme proposed by psychophysicists to account for FM sensitivity in humans and provide new insights on how this might be implemented in the early stages of the auditory pathway.