The roles potassium currents play in regulating the electrical activity of ventral cochlear nucleus neurons

Molecular and functional profiling of cell diversity and identity in the lateral superior olive, an auditory brainstem center with ascending and descending projections

The lateral superior olive (LSO), a prominent integration center in the auditory brainstem, contains a remarkably heterogeneous population of neurons. Ascending neurons, predominantly principal neurons (pLSOs), process interaural level differences for sound localization. Descending neurons (lateral olivocochlear neurons, LOCs) provide feedback into the cochlea and are thought to protect against acoustic overload. The molecular determinants of the neuronal diversity in the LSO are largely unknown. Here, we used patch-seq analysis in mice at postnatal days P10-12 to classify developing LSO neurons according to their functional and molecular profiles. Across the entire sample (n = 86 neurons), genes involved in ATP synthesis were particularly highly expressed, confirming the energy expenditure of auditory neurons. Two clusters were identified, pLSOs and LOCs. They were distinguished by 353 differentially expressed genes (DEGs), most of which were novel for the LSO. Electrophysiological analysis confirmed the transcriptomic clustering. We focused on genes affecting neuronal input-output properties and validated some of them by immunohistochemistry, electrophysiology, and pharmacology. These genes encode proteins such as osteopontin, Kv11.3, and Kvβ3 (pLSO-specific), calcitonin-gene-related peptide (LOC-specific), or Kv7.2 and Kv7.3 (no DEGs). We identified 12 "Super DEGs" and 12 genes showing "Cluster similarity." Collectively, we provide fundamental and comprehensive insights into the molecular composition of individual ascending and descending neurons in the juvenile auditory brainstem and how this may relate to their specific functions, including developmental aspects.

Dynamics of interaction between IH and IKLT currents to mediate double resonances of medial superior olive neurons related to sound localization

Neurons in the medial superior olive (MSO) exhibit high frequency responses such as subthreshold resonance, which is helpful to sensitively detect a small difference in the arrival time of sounds between two ears for precise sound localization. Recently, except for the high frequency depolarization resonance mediated by a low threshold potassium (IKLT) current, a low frequency hyperpolarization resonance mediated by a hyperpolarization-activated cation (IH) current is observed in experiments on the MSO neurons, forming double resonances. The complex dynamics underlying double resonances are studied in an MSO neuron model in the present paper. Firstly, double resonances similar to the experimental observations are simulated as the resting membrane potential is between half-activation voltages of IH and IKLT currents, and stimulation current (IZAP) with large amplitude and exponentially increasing frequency is applied. Secondly, multiple effective factors to modulate double resonances are obtained. Especially, the decrease of time constant of IKLT current and increase of conductance of IH and IKLT currents can enhance the depolarization resonance frequency for precise sound localization. Last, different frequency responses of slow IH and fast IKLT currents in formation of the resonances are acquired. A middle phase difference between IZAP and IKLT currents appears at a high frequency, and the interaction between the positive part of IZAP and the negative IKLT current forms the depolarization resonance. Interaction between the negative part of IZAP and positive IH current with a middle phase difference results in hyperpolarization resonance at a low frequency. Furthermore, the phase difference between IZAP and resonance current can well explain the increase of depolarization resonance frequency modulated by the increase of conductance of IH or IKLT currents. The results present the dynamical and biophysical mechanisms for the double resonances mediated by two currents in the MSO neurons, which is helpful to enhance the depolarization resonance frequency for precise sound localization.

Effects of transient, persistent, and resurgent sodium currents on excitability and spike regularity in vestibular ganglion neurons

Vestibular afferent neurons occur as two populations, regular and irregular, that provide distinct information about head motions. Differences in spike timing regularity are correlated with the different sensory responses important for vestibular processing. Relative to irregular afferents, regular afferents have more sustained firing patterns in response to depolarizing current steps, are more excitable, and have different complements of ion channels. Models of vestibular regularity and excitability emphasize the influence of increased expression of low-voltage-activated potassium currents in irregular neurons. We investigated the potential impact of different modes of voltage-gated sodium (NaV) current (transient, persistent, and resurgent) in cell bodies from vestibular ganglion neurons (VGNs), dissociated and cultured overnight. We hypothesized that regular VGNs would show the greatest impact of persistent (non-inactivating) NaV currents and of resurgent NaV currents, which flow when NaV channels are blocked and then unblocked. Whole-cell patch clamp experiments showed that much of the NaV current modes is carried by NaV1.6 channels. With simulations, we detected little substantial effect in any model VGN of persistent or resurgent modes on regularity of spike timing driven by postsynaptic current trains. For simulated irregular neurons, we also saw little effect on spike rate or firing pattern. For simulated regular VGNs, adding resurgent current changed the detailed timing of spikes during a current step, while the small persistent conductance (less than10% of transient NaV conductance density) strongly depolarized resting potential, altered spike waveform, and increased spike rate. These results suggest that persistent and resurgent NaV current can have a greater effect on the regular VGNs than on irregular VGNs, where low-voltage-activated K conductances dominate.

Comprehensive analysis of cellular specializations that initiate parallel auditory processing pathways in mice

The cochlear nuclear complex (CN) is the starting point for all central auditory processing and comprises a suite of neuronal cell types that are highly specialized for neural coding of acoustic signals. To examine how their striking functional specializations are determined at the molecular level, we performed single-nucleus RNA sequencing of the mouse CN to molecularly define all constituent cell types and related them to morphologically- and electrophysiologically-defined neurons using Patch-seq. We reveal an expanded set of molecular cell types encompassing all previously described major types and discover new subtypes both in terms of topographic and cell-physiologic properties. Our results define a complete cell-type taxonomy in CN that reconciles anatomical position, morphological, physiological, and molecular criteria. This high-resolution account of cellular heterogeneity and specializations from the molecular to the circuit level illustrates molecular underpinnings of functional specializations and enables genetic dissection of auditory processing and hearing disorders with unprecedented specificity.

Rapid, Activity-Dependent Intrinsic Plasticity in the Developing Zebra Finch Auditory Cortex

The acoustic environment an animal experiences early in life shapes the structure and function of its auditory system. This process of experience-dependent development is thought to be primarily orchestrated by potentiation and depression of synapses, but plasticity of intrinsic voltage dynamics may also contribute. Here, we show that in juvenile male and female zebra finches, neurons in a cortical-level auditory area, the caudal mesopallium (CM), can rapidly change their firing dynamics. This plasticity was only observed in birds that were reared in a complex acoustic and social environment, which also caused increased expression of the low-threshold potassium channel Kv1.1 in the plasma membrane and endoplasmic reticulum (ER). Intrinsic plasticity depended on activity, was reversed by blocking low-threshold potassium currents, and was prevented by blocking intracellular calcium signaling. Taken together, these results suggest that Kv1.1 is rapidly mobilized to the plasma membrane by activity-dependent elevation of intracellular calcium. This produces a shift in the excitability and temporal integration of CM neurons that may be permissive for auditory learning in complex acoustic environments during a crucial period for the development of vocal perception and production.SIGNIFICANCE STATEMENT Neurons can change not only the strength of their connections to other neurons, but also how they integrate synaptic currents to produce patterns of action potentials. In contrast to synaptic plasticity, the mechanisms and functional roles of intrinisic plasticity remain poorly understood. We found that neurons in the zebra finch auditory cortex can rapidly shift their spiking dynamics within a few minutes in response to intracellular stimulation. This plasticity involves increased conductance of a low-threshold potassium current associated with the Kv1.1 channel, but it only occurs in birds reared in a rich acoustic environment. Thus, auditory experience regulates a mechanism of neural plasticity that allows neurons to rapidly adapt their firing dynamics to stimulation.

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.

Principal neuron diversity in the murine lateral superior olive supports multiple sound localization strategies and segregation of information in higher processing centers

Principal neurons (PNs) of the lateral superior olive nucleus (LSO) in the brainstem of mammals compare information between the two ears and enable sound localization on the horizontal plane. The classical view of the LSO is that it extracts ongoing interaural level differences (ILDs). Although it has been known for some time that LSO PNs have intrinsic relative timing sensitivity, recent reports further challenge conventional thinking, suggesting the major function of the LSO is detection of interaural time differences (ITDs). LSO PNs include inhibitory (glycinergic) and excitatory (glutamatergic) neurons which differ in their projection patterns to higher processing centers. Despite these distinctions, intrinsic property differences between LSO PN types have not been explored. The intrinsic cellular properties of LSO PNs are fundamental to how they process and encode information, and ILD/ITD extraction places disparate demands on neuronal properties. Here we examine the ex vivo electrophysiology and cell morphology of inhibitory and excitatory LSO PNs in mice. Although overlapping, properties of inhibitory LSO PNs favor time coding functions while those of excitatory LSO PNs favor integrative level coding. Inhibitory and excitatory LSO PNs exhibit different activation thresholds, potentially providing further means to segregate information in higher processing centers. Near activation threshold, which may be physiologically similar to the sensitive transition point in sound source location for LSO, all LSO PNs exhibit single-spike onset responses that can provide optimal time encoding ability. As stimulus intensity increases, LSO PN firing patterns diverge into onset-burst cells, which can continue to encode timing effectively regardless of stimulus duration, and multi-spiking cells, which can provide robust individually integrable level information. This bimodal response pattern may produce a multi-functional LSO which can encode timing with maximum sensitivity and respond effectively to a wide range of sound durations and relative levels.

Dendritic Degeneration and Altered Synaptic Innervation of a Central Auditory Neuron During Age-related Hearing Loss

Cellular morphology and synaptic configuration are key determinants of neuronal function and are often modified under pathological conditions. In the first nucleus of the central auditory system, the cochlear nucleus (CN), principal bushy neurons specialize in processing temporal information of sound critical for hearing. These neurons alter their physiological properties during aging that contribute to age-related hearing loss (ARHL). The structural basis of such changes remains unclear, especially age-related modifications in their dendritic morphology and the innervating auditory nerve (AN) synapses. Using young (2-5 months) and aged (28-33 months) CBA/CaJ mice of either sex, we filled individual bushy neurons with fluorescent dye in acute brain slices to characterize their dendritic morphology, followed by immunostaining against vesicular glutamate transporter 1 (VGluT1) and calretinin (CR) to identify innervating AN synapses. We found that dendritic morphology of aged bushy neurons had significantly reduced complexity, suggesting age-dependent dendritic degeneration, especially in neurons with predominantly non-CR-expressing synapses on the soma. These dendrites were innervated by AN bouton synapses, which were predominantly non-CR-expressing in young mice but had increased proportion of CR-expressing synapses in old mice. While somatic AN synapses degenerated substantially with age, as quantified by VGluT1-labeled puncta volume, no significant difference was observed in the total volume of dendritic synapses between young and old mice. Consequently, synaptic density on dendrites was significantly higher in old mice. The findings suggest that dendritic degeneration and altered synaptic innervation in bushy neurons during aging may underlie their changed physiological activity and contribute to the development of ARHL.

Muscarinic Acetylcholine Receptors Modulate HCN Channel Properties in Vestibular Ganglion Neurons

Efferent modulation of vestibular afferent excitability is linked to muscarinic signaling cascades that close low-voltage-gated potassium channels (i.e., KCNQ). Here, we show that muscarinic signaling cascades also depolarize the activation range of hyperpolarization-activated cyclic-nucleotide gated (HCN) channels. We compared the voltage activation range and kinetics of HCN channels and induced firing patterns before and after administering the muscarinic acetylcholine receptor (mAChR) agonist oxotremorine-M (Oxo-M) in dissociated vestibular ganglion neurons (VGNs) from rats of either sex using perforated whole-cell patch-clamp methods. Oxo-M depolarized HCN channels' half-activation voltage (V 1/2) and sped up the rate of activation near resting potential twofold. HCN channels in large-diameter and/or transient firing VGN (putative cell bodies of irregular firing neuron from central epithelial zones) had relatively depolarized V 1/2 in control solution and were less sensitive to mAChR activation than those found in small-diameter VGN with sustained firing patterns (putatively belonging to regular firing afferents). The impact of mAChR on HCN channels is not a direct consequence of closing KCNQ channels since pretreating the cells with Linopirdine, a KCNQ channel blocker, did not prevent HCN channel depolarization by Oxo-M. Efferent signaling promoted ion channel configurations that were favorable to highly regular spiking in some VGN, but not others. This is consistent with previous observations that low-voltage gated potassium currents in VGN are conducted by mAChR agonist-sensitive and -insensitive channels. Connecting efferent signaling to HCN channels is significant because of the channel's impact on spike-timing regularity and nonchemical transmission between Type I hair cells and vestibular afferents.SIGNIFICANCE STATEMENT Vestibular afferents express a diverse complement of ion channels. In vitro studies identified low-voltage activated potassium channels and hyperpolarization-activated cyclic-nucleotide gated (HCN) channels as crucial for shaping the timing and sensitivity of afferent responses. Moreover, a network of acetylcholine-releasing efferent neurons controls afferent excitability by closing a subgroup of low-voltage activated potassium channels on the afferent neuron. This work shows that these efferent signaling cascades also enhance the activation of HCN channels by depolarizing their voltage activation range. The size of this effect varies depending on the endogenous properties of the HCN channel and on cell type (as determined by discharge patterns and cell size). Simultaneously controlling two ion-channel groups gives the vestibular efferent system exquisite control over afferent neuron activity.

Low-voltage Activating K+ Channels in Cochlear Afferent Nerve Fiber Dendrites

Cochlear afferent nerve fibers (ANF) are the first neurons in the ascending auditory pathway. We investigated the low-voltage activating K⁺ channels expressed in ANF dendrites using isolated rat cochlear segments. Whole cell patch clamp recordings were made from the dendritic terminals of ANFs. Outward currents activating at membrane potentials as low as -64 mV were observed in all dendrites studied. These currents were inhibited by 4-aminopyridine (4-AP), a blocker known to preferentially inhibit low-voltage activating K⁺ currents (I_KL) in CNS auditory neurons and spiral ganglion neurons. When the dendritic I_KL was blocked by 4-AP, the EPSP decay time was significantly prolonged, suggesting that dendritic I_KL speeds up the decay of EPSPs and likely modulates action potentials of ANFs. To reveal molecular subtype of dendritic I_KL, α-dendrotoxin (α-DTX), a selective inhibitor for K_v1.1, K_v1.2, and K_v1.6 containing channels, was tested. α-DTX inhibited 23±9% of dendritic I_KL. To identify the α-DTXsensitive and α-DTX-insensitive components of I_KL, immunofluorescence labeling was performed. Strong K_v1.1- and K_v1.2-immunoreactivity was found at unmyelinated dendritic segments, nodes of Ranvier, and cell bodies of most ANFs. A small fraction of ANF dendrites showed K_v7.2- immunoreactivity. These data suggest that dendritic I_KL is conducted through K_v1.1and K_v1.2 channels, with a minor contribution from K_v7.2 and other as yet unidentified channels.

Distinct neuronal types contribute to hybrid temporal encoding strategies in primate auditory cortex

Studies of the encoding of sensory stimuli by the brain often consider recorded neurons as a pool of identical units. Here, we report divergence in stimulus-encoding properties between subpopulations of cortical neurons that are classified based on spike timing and waveform features. Neurons in auditory cortex of the awake marmoset (Callithrix jacchus) encode temporal information with either stimulus-synchronized or nonsynchronized responses. When we classified single-unit recordings using either a criteria-based or an unsupervised classification method into regular-spiking, fast-spiking, and bursting units, a subset of intrinsically bursting neurons formed the most highly synchronized group, with strong phase-locking to sinusoidal amplitude modulation (SAM) that extended well above 20 Hz. In contrast with other unit types, these bursting neurons fired primarily on the rising phase of SAM or the onset of unmodulated stimuli, and preferred rapid stimulus onset rates. Such differentiating behavior has been previously reported in bursting neuron models and may reflect specializations for detection of acoustic edges. These units responded to natural stimuli (vocalizations) with brief and precise spiking at particular time points that could be decoded with high temporal stringency. Regular-spiking units better reflected the shape of slow modulations and responded more selectively to vocalizations with overall firing rate increases. Population decoding using time-binned neural activity found that decoding behavior differed substantially between regular-spiking and bursting units. A relatively small pool of bursting units was sufficient to identify the stimulus with high accuracy in a manner that relied on the temporal pattern of responses. These unit type differences may contribute to parallel and complementary neural codes.

Neurons in auditory cortex show highly diverse responses to sounds. This study suggests that neuronal type inferred from baseline firing properties accounts for much of this diversity, with a subpopulation of bursting units being specialized for precise temporal encoding.

Predicting Spike Features of Hodgkin-Huxley-Type Neurons With Simple Artificial Neural Network

Hodgkin-Huxley (HH)-type model is the most famous computational model for simulating neural activity. It shows the highest accuracy in capturing neuronal spikes, and its model parameters have definite physiological meanings. However, HH-type models are computationally expensive. To address this problem, a previous study proposed a spike prediction module (SPM) to predict whether a spike will take place 1 ms later based on three voltage values with intervals of 1 ms. Although SPM does well, it fails to evaluate the informative features of the spike. In this study, the feature prediction module (FPM) based on simple artificial neural network (ANN) was proposed to predict spike features including maximum voltage, minimum voltage, and dropping interval. Nine different HH-type models were adopted whose firing patterns cover most of the firing behaviors observed in the brain. Voltage and spike feature samples under constant external input current were collected for training and testing. Experiment results illustrated that the combination of SPM and FPM can accurately predict the spiking part of different HH-type models and can generalize to unseen types of input current. The combination of SPM and FPM may offer a possible way to simulate the action potentials of biological neurons with high accuracy and efficiency.

Sensitivity to Envelope Interaural Time Differences: Modeling Auditory Modulation Filtering

For amplitude-modulated sound, the envelope interaural time difference (ITD_ENV) is a potential cue for sound-source location. ITD_ENV is encoded in the lateral superior olive (LSO) of the auditory brainstem, by excitatory-inhibitory (EI) neurons receiving ipsilateral excitation and contralateral inhibition. Between human listeners, sensitivity to ITD_ENV varies considerably, but ultimately decreases with increasing stimulus carrier frequency, and decreases more strongly with increasing modulation rate. Mechanisms underlying the variation in behavioral sensitivity remain unclear. Here, with increasing carrier frequency (4-10 kHz), as we phenomenologically model the associated decrease in ITD_ENV sensitivity using arbitrarily fewer neurons consistent across populations, we computationally model the variable sensitivity across human listeners and modulation rates (32-800 Hz) as the decreasing range of membrane frequency responses in LSO neurons. Transposed tones stimulate a bilateral auditory-periphery model, driving model EI neurons where electrical membrane impedance filters the frequency content of inputs driven by amplitude-modulated sound, evoking modulation filtering. Calculated from Fisher information in spike-rate functions of ITD_ENV, for model EI neuronal populations distinctly reflecting the LSO range in membrane frequency responses, just-noticeable differences in ITD_ENV collectively reproduce the largest variation in ITD_ENV sensitivity across human listeners. These slow to fast model populations each generally match the best human ITD_ENV sensitivity at a progressively higher modulation rate, by membrane-filtering and spike-generation properties producing realistically less than Poisson variance. Non-resonant model EI neurons are also sensitive to interaural intensity differences. With peripheral filters centered between carrier frequency and modulation sideband, fast resonant model EI neurons extend ITD_ENV sensitivity above 500-Hz modulation.

Kv1 channels regulate variations in spike patterning and temporal reliability in the avian cochlear nucleus angularis

Diverse physiological phenotypes in a neuronal population can broaden the range of computational capabilities within a brain region. The avian cochlear nucleus angularis (NA) contains a heterogeneous population of neurons whose variation in intrinsic properties results in electrophysiological phenotypes with a range of sensitivities to temporally modulated input. The low-threshold potassium conductance (G_KLT) is a key feature of neurons involved in fine temporal structure coding for sound localization, but a role for these channels in intensity or spectrotemporal coding has not been established. To determine whether G_KLT affects the phenotypical variation and temporal properties of NA neurons, we applied dendrotoxin-I (DTX), a potent antagonist of Kv1-type potassium channels, to chick brain stem slices in vitro during whole cell patch-clamp recordings. We found a cell-type specific subset of NA neurons that was sensitive to DTX: single-spiking NA neurons were most profoundly affected, as well as a subset of tonic-firing neurons. Both tonic I (phasic onset bursting) and tonic II (delayed firing) neurons showed DTX sensitivity in their firing rate and phenotypical firing pattern. Tonic III neurons were unaffected. Spike time reliability and fluctuation sensitivity measured in DTX-sensitive NA neurons was also reduced with DTX. Finally, DTX reduced spike threshold adaptation in these neurons, suggesting that G_KLT contributes to the temporal properties that allow coding of rapid changes in the inputs to NA neurons. These results suggest that variation in Kv1 channel expression may be a key factor in functional diversity in the avian cochlear nucleus.NEW & NOTEWORTHY The dendrotoxin-sensitive voltage-gated potassium conductance typically associated with neuronal coincidence detection in the timing pathway for sound localization is demonstrated to affect spiking patterns and temporal input sensitivity in the intensity pathway in the avian auditory brain stem. The Kv1-family channels appear to be present in a subset of cochlear nucleus angularis neurons, regulate spike threshold dynamics underlying high-pass membrane filtering, and contribute to intrinsic firing diversity.

Dendritic synapse geometry optimizes binaural computation in a sound localization circuit

Clustering of synapses allows neurons to overcome attenuation of electrical signals at dendrites. However, we show in avian binaural coincidence detectors computing interaural time difference for sound localization that clustering of synapses rather promotes the dendritic attenuation but augments the intensity tolerance of the binaural computations. Using glutamate uncaging, we found in the neurons that synapses were clustered at distal dendritic branches. Modeling revealed that this strengthened sublinear integration within a dendritic tree but enabled the integration of signals from different trees when inputs grow stronger, preventing monoaural output and maintaining the dynamic range of binaural computation. The extent of this clustering differed according to dendritic length and frequency tuning of neurons, being most prominent for long dendrites and low-frequency tuning. This ensures binaural spatial hearing for wide intensity and frequency ranges, highlighting the importance of coupling of synapse geometry with dendritic morphology and input frequency in sensory signal processing.

Similarities in the Biophysical Properties of Spiral-Ganglion and Vestibular-Ganglion Neurons in Neonatal Rats

The membranes of auditory and vestibular afferent neurons each contain diverse groups of ion channels that lead to heterogeneity in their intrinsic biophysical properties. Pioneering work in both auditory- and vestibular-ganglion physiology have individually examined this remarkable diversity, but there are few direct comparisons between the two ganglia. Here the firing patterns recorded by whole-cell patch-clamping in neonatal vestibular- and spiral ganglion neurons are compared. Indicative of an overall heterogeneity in ion channel composition, both ganglia exhibit qualitatively similar firing patterns ranging from sustained-spiking to transient-spiking in response to current injection. The range of resting potentials, voltage thresholds, current thresholds, input-resistances, and first-spike latencies are similarly broad in both ganglion groups. The covariance between several biophysical properties (e.g., resting potential to voltage threshold and their dependence on postnatal age) was similar between the two ganglia. Cell sizes were on average larger and more variable in VGN than in SGN. One sub-group of VGN stood out as having extra-large somata with transient-firing patterns, very low-input resistance, fast first-spike latencies, and required large current amplitudes to induce spiking. Despite these differences, the input resistance per unit area of the large-bodied transient neurons was like that of smaller-bodied transient-firing neurons in both VGN and SGN, thus appearing to be size-scaled versions of other transient-firing neurons. Our analysis reveals that although auditory and vestibular afferents serve very different functions in distinct sensory modalities, their biophysical properties are more closely related by firing pattern and cell size than by sensory modality.

Spike-rate adaptation in a computational model of human-shaped spiral ganglion neurons

OBJECTIVES

The purpose of this study is to develop a biophysical model of human spiral ganglion neurons (SGNs) that includes voltage-gated hyperpolarization-activated cation (HCN) channels and low-threshold potassium voltage-gated, delayed-rectifier low-threshold potassium (KLT) channels, providing for a more complete simulation of spike-rate adaptation, a feature of most spiking neurons in which spiking activity is reduced in response to sustained stimulation.

METHODS

Our model incorporates features of spike-rate adaptation reported from in vivo studies, whilst also displaying similar behaviour to existing models of human SGNs, including the dependence of electrically evoked thresholds on the polarity of electrical pulses.

RESULTS

Hypothesizing that the mode of stimulation intracellular or extracellular predicts features of spike-rate adaptation similar to in vivo studies, including the influence of stimulus intensity and pulse-rate, we find that the mode of stimulation alters features of spike-rate adaptation. In particular, the reduction in spiking over time with sustained input was generally greater for extracellular, compared to intracellular, stimulation, when simulating a multi-compartment SGN with human morphological features. In contrast, time-constants of spike-rate adaption reported for in vivo data did not fit our predicted responses, highlighting the need for a more complete physiological understanding of the factors contributing to spike-rate adaptation in electrically stimulated human SGNs.

CONCLUSION

Our model extends previous computational models of SGNs with human morphology with ionic channels accounting for features of spike-rate adaptation.

SIGNIFICANCE

The significance of this work resides in the ability to improve the modeling of cochlear implant (CI) stimulation and its effects on neural responses. This will help develop novel, and perhaps personalised, stimulation strategies to reduce variability in CI user outcomes.

A phenomenological spiking model for octopus cells in the posterior-ventral cochlear nucleus

Octopus cells in the posteroventral cochlear nucleus exhibit characteristic onset responses to broad band transients but are little investigated in response to more complex sound stimuli. In this paper, we propose a phenomenological, but biophysically motivated, modeling approach that allows to simulate responses of large populations of octopus cells to arbitrary sound pressure waves. The model depends on only few parameters and reproduces basic physiological characteristics like onset firing and phase locking to amplitude modulations. Simulated responses to speech stimuli suggest that octopus cells are particularly sensitive to high-frequency transients in natural sounds and their sustained firing to phonemes provides a population code for sound level.

Effects of Several Therapeutic Agents on Mammalian Vestibular Function: Meclizine, Diazepam, and JNJ7777120

Management of vestibular dysfunction may include treatment with medications that are thought to act to suppress vestibular function and reduce or eliminate abnormal sensitivity to head motions. The extent to which vestibular medications act centrally or peripherally is still debated. In this study, two commonly prescribed medications, meclizine and diazepam, and a candidate for future clinical use, JNJ7777120, were evaluated for their effects on short latency compound action potentials generated by the peripheral vestibular system and corresponding central neural relays (i.e., vestibular sensory-evoked potentials, VsEPs). The effects of the selected drugs developed slowly over the course of two hours in the mouse. Findings indicate that meclizine (600 mg/kg) and diazepam (> 60 mg/kg) can act on peripheral elements of the vestibular maculae whereas diazepam also acts most effectively on central gravity receptor circuits to exert its suppressive effects. The novel pharmacological agent JNJ7777120 (160 mg/kg) acts in the vestibular periphery to enhance macular responses to transient stimuli (VsEPs) while, hypothetically, suppressing macular responses to sustained or slowly changing stimuli.

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.

Auditory Brainstem Models: Adapting Cochlear Nuclei Improve Spatial Encoding by the Medial Superior Olive in Reverberation

Listeners typically perceive a sound as originating from the direction of its source, even as direct sound is followed milliseconds later by reflected sound from multiple different directions. Early-arriving sound is emphasised in the ascending auditory pathway, including the medial superior olive (MSO) where binaural neurons encode the interaural-time-difference (ITD) cue for spatial location. Perceptually, weighting of ITD conveyed during rising sound energy is stronger at 600 Hz than at 200 Hz, consistent with the minimum stimulus rate for binaural adaptation, and with the longer reverberation times at 600 Hz, compared with 200 Hz, in many natural outdoor environments. Here, we computationally explore the combined efficacy of adaptation prior to the binaural encoding of ITD cues, and excitatory binaural coincidence detection within MSO neurons, in emphasising ITDs conveyed in early-arriving sound. With excitatory inputs from adapting, nonlinear model spherical bushy cells (SBCs) of the bilateral cochlear nuclei, a nonlinear model MSO neuron with low-threshold potassium channels reproduces the rate-dependent emphasis of rising vs. peak sound energy in ITD encoding; adaptation is equally effective in the model MSO. Maintaining adaptation in model SBCs, and adjusting membrane speed in model MSO neurons, 'left' and 'right' populations of computationally efficient, linear model SBCs and MSO neurons reproduce this stronger weighting of ITD conveyed during rising sound energy at 600 Hz compared to 200 Hz. This hemispheric population model demonstrates a link between strong weighting of spatial information during rising sound energy, and correct unambiguous lateralisation of a speech source in reverberation.

Differential contributions of voltage-gated potassium channel subunits in enhancing temporal coding in the bushy cells of the ventral cochlear nucleus

Temporal coding precision of bushy cells in the ventral cochlear nucleus (VCN), critical for sound localization and communication, depends on the generation of rapid and temporally precise action potentials (APs). Voltage-gated potassium (Kv) channels are critically involved in this. The bushy cells in rat VCN express Kv1.1, 1.2, 1.3, 1.6, 3.1, 4.2, and 4.3 subunits. The Kv1.1 subunit contributes to the generation of a temporally precise single AP. However, the understanding of the functions of other Kv subunits expressed in the bushy cells is limited. Here, we investigated the functional diversity of Kv subunits concerning their contributions to temporal coding. We characterized the electrophysiological properties of the Kv channels with different subunits using whole cell patch-clamp recording and pharmacological methods. The neuronal firing pattern changed from single to multiple APs only when the Kv1.1 subunit was blocked. The Kv subunits, including the Kv1.1, 1.2, 1.6, or 3.1, were involved in enhancing temporal coding by lowering membrane excitability, shortening AP latencies, reducing jitter, and regulating AP kinetics. Meanwhile, all the Kv subunits contributed to rapid repolarization and sharpening peaks by narrowing half-width and accelerating fall rate, and the Kv1.1 subunit also affected the depolarization of AP. The Kv1.1, 1.2, and 1.6 subunits endowed bushy cells with a rapid time constant and a low input resistance of membrane for enhancing spike timing precision. The present results indicate that the Kv channels differentially affect intrinsic membrane properties to optimize the generation of rapid and reliable APs for temporal coding.NEW & NOTEWORTHY This study investigates the roles of Kv channels in effecting precision using electrophysiological and pharmacological methods in bushy cells. Different Kv channels have varying electrophysiological characteristics, which contribute to the interplay between changes in the membrane properties and regulation of neuronal excitability which then improve temporal coding. We conclude that the Kv channels are specialized to promote the precise and rapid coding of acoustic input by optimizing the generation of reliable APs.

Nonlinear effects of intrinsic dynamics on temporal encoding in a model of avian auditory cortex

Neurons exhibit diverse intrinsic dynamics, which govern how they integrate synaptic inputs to produce spikes. Intrinsic dynamics are often plastic during development and learning, but the effects of these changes on stimulus encoding properties are not well known. To examine this relationship, we simulated auditory responses to zebra finch song using a linear-dynamical cascade model, which combines a linear spectrotemporal receptive field with a dynamical, conductance-based neuron model, then used generalized linear models to estimate encoding properties from the resulting spike trains. We focused on the effects of a low-threshold potassium current (K_LT) that is present in a subset of cells in the zebra finch caudal mesopallium and is affected by early auditory experience. We found that K_LT affects both spike adaptation and the temporal filtering properties of the receptive field. The direction of the effects depended on the temporal modulation tuning of the linear (input) stage of the cascade model, indicating a strongly nonlinear relationship. These results suggest that small changes in intrinsic dynamics in tandem with differences in synaptic connectivity can have dramatic effects on the tuning of auditory neurons.

Experience-dependent developmental plasticity involves changes not only to synaptic connections, but to voltage-gated currents as well. Using biophysical models, it is straightforward to predict the effects of this intrinsic plasticity on the firing patterns of individual neurons, but it remains difficult to understand the consequences for sensory coding. We investigated this in the context of the zebra finch auditory cortex, where early exposure to a complex acoustic environment causes increased expression of a low-threshold potassium current. We simulated responses to song using a detailed biophysical model and then characterized encoding properties using generalized linear models. This analysis revealed that this potassium current has strong, nonlinear effects on how the model encodes the song’s temporal structure, and that the sign of these effects depend on the temporal tuning of the synaptic inputs. This nonlinearity gives intrinsic plasticity broad scope as a mechanism for developmental learning in the auditory system.

Modulation of Neuronal Potassium Channels During Auditory Processing

The extraction and localization of an auditory stimulus of interest from among multiple other sounds, as in the ‘cocktail-party’ situation, requires neurons in auditory brainstem nuclei to encode the timing, frequency, and intensity of sounds with high fidelity, and to compare inputs coming from the two cochleae. Accurate localization of sounds requires certain neurons to fire at high rates with high temporal accuracy, a process that depends heavily on their intrinsic electrical properties. Studies have shown that the membrane properties of auditory brainstem neurons, particularly their potassium currents, are not fixed but are modulated in response to changes in the auditory environment. Here, we review work focusing on how such modulation of potassium channels is critical to shaping the firing pattern and accuracy of these neurons. We describe how insights into the role of specific channels have come from human gene mutations that impair localization of sounds in space. We also review how short-term and long-term modulation of these channels maximizes the extraction of auditory information, and how errors in the regulation of these channels contribute to deficits in decoding complex auditory information.

Predicting Response of Spontaneously Firing Afferents to Prosthetic Pulsatile Stimulation

Pulsatile electrical stimulation is used in neural prostheses such as the vestibular prosthesis. In a healthy vestibular system, head motion is encoded by changes in the firing rates of afferents around their spontaneous baseline rate. For people suffering from bilateral vestibular disorder (BVD), head motion no longer modulates firing rate. Vestibular prostheses use a gyroscope to detect head motion and stimulate neurons directly in a way that mimics natural modulation. Proper restoration of vestibular function relies on the ability of stimulation to evoke the same firing patterns as the healthy system. For this reason, it is necessary to understand what firing rates are produced for different stimulation parameters. Two stimulation parameters commonly controlled in pulsatile neuromodulation are pulse rate and pulse amplitude. Previous neural recording experiments in the vestibular nerve contradict widely held assumptions about the relationship between pulse rates and evoked spike activity, and the relationship between pulse amplitude and neural activity has not been explored. Here we use a well-established computational model of the vestibular afferent to simulate responses to different pulse rates and amplitudes. We confirm that our simulated neural results agree with the existing experimental data. Finally, we developed the "Action Potential Collision" (APC) equation that defines induced firing as a function of spontaneous firing rate, pulse rate, and pulse amplitude. We show that this relationship can successfully predict simulated vestibular activity by accounting for interactions between pulses and spontaneous firing.

How to pick a peak: Pitch and peak shifting in temporal models of pitch perception

The standard autocorrelation model of pitch perception posits that the pitch of a stimulus can be predicted from the first major peak of a summary autocorrelation function (SACF) after the zero-delay peak. Models based on this theory are capable of predicting a wide range of pitch phenomena. There are, however, a number of cases where the approach fails. Two examples are noise edge pitch (NEP) and the pitch induced by the mistuning of a single component of an otherwise harmonic stimulus. Hartmann, Cariani, and Colburn [(2019). J. Acoust. Soc. Am. 145, 1993-2008] recently proposed the use of multiple SACF peaks in the estimation process. This enables prediction of the NEP but suppresses the shift associated with a mistuned harmonic. A functional model is developed that can predict both of these pitch phenomena. The multiple-peak framework is extended with a non-standard peak-selection method that associates a delay time to a given peak in a manner that takes into account the entire shape of the bump surrounding the peak. This effectively shifts the peak location slightly for non-harmonic stimuli. A possible physiological mechanism that could induce such peak shifting is discussed, and the model is tested against existing psychophysical data.

Modeling the Effect of Temperature on Membrane Response of Light Stimulation in Optogenetically-Targeted Neurons

Optogenetics is revolutionizing Neuroscience, but an often neglected effect of light stimulation of the brain is the generation of heat. In extreme cases, light-generated heat kills neurons, but mild temperature changes alter neuronal function. To date, most in vivo experiments rely on light stimulation of neural tissue using fiber-coupled lasers of various wavelengths. Brain tissue is irradiated with high light power that can be deleterious to neuronal function. Furthermore, absorbed light generates heat that can lead to permanent tissue damage and affect neuronal excitability. Thus, light alone can generate effects in neuronal function that are unrelated to the genuine "optogenetic effect." In this work, we perform a theoretical analysis to investigate the effects of heat transfer in rodent brain tissue for standard optogenetic protocols. More precisely, we first use the Kubelka-Munk model for light propagation in brain tissue to observe the absorption phenomenon. Then, we model the optothermal effect considering the common laser wavelengths (473 and 593 nm) used in optogenetic experiments approaching the time/space numerical solution of Pennes' bio-heat equation with the Finite Element Method. Finally, we then modeled channelrhodopsin-2 in a single and spontaneous-firing neuron to explore the effect of heat in light stimulated neurons. We found that, at commonly used light intensities, laser radiation considerably increases the temperature in the surrounding tissue. This effect alters action potential size and shape and causes an increase in spontaneous firing frequency in a neuron model. However, the shortening of activation time constants generated by heat in the single firing neuron model produces action potential failures in response to light stimulation. We also found changes in the power spectrum density and a reduction in the time required for synchronization in an interneuron network model of gamma oscillations. Our findings indicate that light stimulation with intensities used in optogenetic experiments may affect neuronal function not only by direct excitation of light sensitive ion channels and/or pumps but also by generating heat. This approach serves as a guide to design optogenetic experiments that minimize the role of tissue heating in the experimental outcome.

Nitric oxide regulates the firing rate of neuronal subtypes in the guinea pig ventral cochlear nucleus

The gaseous free radical, nitric oxide (NO) acts as a ubiquitous neuromodulator, contributing to synaptic plasticity in a complex way that can involve either long term potentiation or depression. It is produced by neuronal nitric oxide synthase (nNOS) which is presynaptically expressed and also located postsynaptically in the membrane and cytoplasm of a subpopulation of each major neuronal type in the ventral cochlear nucleus (VCN). We have used iontophoresis in vivo to study the effect of the NOS inhibitor L-NAME (L-NG-Nitroarginine methyl ester) and the NO donors SIN-1 (3-Morpholinosydnonimine hydrochloride) and SNOG (S-Nitrosoglutathione) on VCN units under urethane anaesthesia. Collectively, both donors produced increases and decreases in driven and spontaneous firing rates of some neurones. Inhibition of endogenous NO production with L-NAME evoked a consistent increase in driven firing rates in 18% of units without much effect on spontaneous rate. This reduction of gain produced by endogenous NO was mirrored when studying the effect of L-NAME on NMDA(N-Methyl-D-aspartic acid)-evoked excitation, with 30% of units showing enhanced NMDA-evoked excitation during L-NAME application (reduced NO levels). Approximately 25% of neurones contain nNOS and the NO produced can modulate the firing rate of the main principal cells: medium stellates (choppers), large stellates (onset responses) and bushy cells (primary-like responses). The main endogenous role of NO seems to be to partly suppress driven firing rates associated with NMDA channel activity but there is scope for it to increase neural gain if there were a pathological increase in its production following hearing loss.

Excitatory-Inhibitory Synaptic Coupling in Avian Nucleus Magnocellularis

The activity of neurons is determined by the balance between their excitatory and inhibitory synaptic inputs. Neurons in the avian nucleus magnocellularis (NM) integrate monosynaptic excitatory and polysynaptic inhibitory inputs from the auditory nerve, and transmit phase-locked output to higher auditory centers. The excitatory input is graded tonotopically, such that neurons tuned to higher frequency receive fewer, but larger, axon terminals. However, it remains unknown how the balance between excitatory and inhibitory inputs is determined in NM. We here examined synaptic and spike responses of NM neurons during stimulation of the auditory nerve in thick brain slices of chicken of both sexes, and found that the excitatory-inhibitory balance varied according to tonotopic region, ensuring reliable spike output across frequencies. Auditory nerve stimulation elicited IPSCs in NM neurons regardless of tonotopic region, but the dependence of IPSCs on intensity varied in a systematic way. In neurons tuned to low frequency, IPSCs appeared and increased in parallel with EPSCs with elevation of intensity, which expanded dynamic range by preventing saturation of spike generation. On the other hand, in neurons tuned to higher frequency, IPSCs were smaller than EPSCs and had higher thresholds for activation, thus facilitating high-fidelity transmission. Computer simulation confirmed that these differences in inhibitory input were optimally matched to the patterns of excitatory input, and enabled appropriate level of neuronal output for wide intensity and frequency ranges of sound in the auditory system.SIGNIFICANCE STATEMENT Neurons in nucleus magnocellularis encode timing information of sound across wide intensity ranges by integrating excitatory and inhibitory synaptic inputs from the auditory nerve, but underlying synaptic mechanisms of this integration are not fully understood. We here show that the excitatory-inhibitory relationship was expressed differentially at each tonotopic region; the relationship was linear in neurons tuned to low-frequency, expanding dynamic range by preventing saturation of spike generation; by contrast inhibitory input remained much smaller than excitatory input in neurons tuned to higher frequency, thus ensuring high-fidelity transmission. The tonotopic regulation of excitatory and inhibitory input optimized the output across frequencies and intensities, playing a fundamental role in the timing coding pathway in the auditory system.

Experience- and Sex-Dependent Intrinsic Plasticity in the Zebra Finch Auditory Cortex during Song Memorization

For vocal communicators like humans and songbirds, survival and reproduction depend on highly developed auditory processing systems that can detect and differentiate nuanced differences in vocalizations, even amid noisy environments. Early auditory experience is critical to the development of these systems. In zebra finches and other songbirds, there is a sensitive period when young birds memorize a song that will serve as a model for their own vocal production. In addition to learning a specific tutor's song, the auditory system may also undergo critical developmental processes that support auditory perception of vocalizations more generally. Here, we investigate changes in intrinsic spiking dynamics among neurons in the caudal mesopallium, a cortical-level auditory area implicated in discriminating and learning species-specific vocalizations. A subset of neurons in this area only fire transiently at the onset of current injections (i.e., phasic firing), a dynamical property that can enhance the reliability and selectivity of neural responses to complex acoustic stimuli. At the beginning of the sensitive period, just after zebra finches have fledged from the nest, there is an increase in the proportion of caudal mesopallium neurons with phasic excitability, and in the proportion of neurons expressing Kv1.1, a low-threshold channel that facilitates phasic firing. This plasticity requires exposure to a complex, noisy environment and is greater in males, the only sex that sings in this species. This shift to more phasic dynamics is therefore an experience-dependent adaptation that could facilitate auditory processing in noisy, acoustically complex conditions during a key stage of vocal development.SIGNIFICANCE STATEMENT Auditory experience early in life shapes how humans and songbirds perceive the vocal communication sounds produced by their species. However, the changes that occur in the brain as this learning takes place are poorly understood. In this study, we show that in young zebra finches that are just beginning to learn the structure of their species' song, neurons in a key cortical area adapt their intrinsic firing patterns in response to the acoustic environment. In the complex, cocktail-party-like environment of a colony, more neurons adopt transient firing dynamics, which can facilitate neural coding of songs amid such challenging conditions.

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.

Classification of neurons in the adult mouse cochlear nucleus: Linear discriminant analysis

The cochlear nucleus (CN) transforms the spike trains of spiral ganglion cells into a set of sensory representations that are essential for auditory discriminations and perception. These transformations require the coordinated activity of different classes of neurons that are embryologically derived from distinct sets of precursors. Decades of investigation have shown that the neurons of the CN are differentiated by their morphology, neurotransmitter receptors, ion channel expression and intrinsic excitability. In the present study we have used linear discriminant analysis (LDA) to perform an unbiased analysis of measures of the responses of CN neurons to current injections to objectively categorize cells on the basis of both morphology and physiology. Recordings were made from cells in brain slices from CBA/CaJ mice and a transgenic mouse line, NF107, crossed against the Ai32 line. For each cell, responses to current injections were analyzed for spike rate, spike shape, input resistance, resting membrane potential, membrane time constant, hyperpolarization-activated sag and time constant. Cells were filled with dye for morphological classification, and visually classified according to published accounts. The different morphological classes of cells were separated with the LDA. Ventral cochlear nucleus (VCN) bushy cells, planar multipolar (T-stellate) cells, and radiate multipolar (D-stellate) cells were in separate clusters and separate from all of the neurons from the dorsal cochlear nucleus (DCN). Within the DCN, the pyramidal cells and tuberculoventral cells were largely separated from a distinct cluster of cartwheel cells. principal axes, whereas VCN cells were in 3 clouds approximately orthogonal to this plane. VCN neurons from the two mouse strains overlapped but were slightly separated, indicating either a strain dependence or differences in slice preparation methods. We conclude that cochlear nucleus neurons can be objectively distinguished based on their intrinsic electrical properties, but such distinctions are still best aided by morphological identification.

Coding strategies in the otolith system differ for translational head motion vs. static orientation relative to gravity

The detection of gravito-inertial forces by the otolith system is essential for our sense of balance and accurate perception. To date, however, how this system encodes the self-motion stimuli that are experienced during everyday activities remains unknown. Here, we addressed this fundamental question directly by recording from single otolith afferents in monkeys during naturalistic translational self-motion and changes in static head orientation. Otolith afferents with higher intrinsic variability transmitted more information overall about translational self-motion than their regular counterparts, owing to stronger nonlinearities that enabled precise spike timing including phase locking. By contrast, more regular afferents better discriminated between different static head orientations relative to gravity. Using computational methods, we further demonstrated that coupled increases in intrinsic variability and sensitivity accounted for the observed functional differences between afferent classes. Together, our results indicate that irregular and regular otolith afferents use different strategies to encode naturalistic self-motion and static head orientation relative to gravity.

Cost of auditory sharpness: Model-Based estimate of energy use by auditory brainstem "octopus" neurons

Octopus cells (OCs) of the mammalian auditory brainstem precisely encode timing of fast transient sounds and tone onsets. Sharp temporal fidelity of OCs relies on low resting membrane resistance, which suggests high energy expenditure on maintaining ion gradients across plasma membrane. We provide a model-based estimate of energy consumption in resting and spiking OCs. Our results predict that a resting OC consumes up to 2.6 × 10⁹ ATP molecules (ATPs) per second which remarkably exceeds energy consumption of other CNS neurons. Glucose usage by all OCs in the rat is nevertheless low due to their low number. Major part of the OCs energy use results from the ion mechanisms providing for the low membrane resistance: hyperpolarization-activated mixed cation conductance and low-voltage activated K⁺-conductance. Spatially ordered synapses-a feature of the OCs allowing them to compensate for asynchrony of the synaptic input-brings only a 12% energy saving to OCs excitability cost. Only 13% of total OC energy used for an AP generation (1.5 × 10⁷ ATPs) is associated with the AP generation in the axon initial segment, 64%-with synaptic currents processing and 23%-with keeping resting potential.

Selective attenuation of Ether-a-go-go related K+ currents by endogenous acetylcholine reduces spike-frequency adaptation and network correlation

Most neurons do not simply convert inputs into firing rates. Instead, moment-to-moment firing rates reflect interactions between synaptic inputs and intrinsic currents. Few studies investigated how intrinsic currents function together to modulate output discharges and which of the currents attenuated by synthetic cholinergic ligands are actually modulated by endogenous acetylcholine (ACh). In this study we optogenetically stimulated cholinergic fibers in rat neocortex and find that ACh enhances excitability by reducing Ether-à-go-go Related Gene (ERG) K+ current. We find ERG mediates the late phase of spike-frequency adaptation in pyramidal cells and is recruited later than both SK and M currents. Attenuation of ERG during coincident depolarization and ACh release leads to reduced late phase spike-frequency adaptation and persistent firing. In neuronal ensembles, attenuating ERG enhanced signal-to-noise ratios and reduced signal correlation, suggesting that these two hallmarks of cholinergic function in vivo may result from modulation of intrinsic properties.

Functional Development of Principal Neurons in the Anteroventral Cochlear Nucleus Extends Beyond Hearing Onset

Sound information is transduced into graded receptor potential by cochlear hair cells and encoded as discrete action potentials of auditory nerve fibers. In the cochlear nucleus, auditory nerve fibers convey this information through morphologically distinct synaptic terminals onto bushy cells (BCs) and stellate cells (SCs) for processing of different sound features. With expanding use of transgenic mouse models, it is increasingly important to understand the in vivo functional development of these neurons in mice. We characterized the maturation of spontaneous and acoustically evoked activity in BCs and SCs by acquiring single-unit juxtacellular recordings between hearing onset (P12) and young adulthood (P30) of anesthetized CBA/J mice. In both cell types, hearing sensitivity and characteristic frequency (CF) range are mostly adult-like by P14, consistent with rapid maturation of the auditory periphery. In BCs, however, some physiological features like maximal firing rate, dynamic range, temporal response properties, recovery from post-stimulus depression, first spike latency (FSL) and encoding of sinusoid amplitude modulation undergo further maturation up to P18. In SCs, the development of excitatory responses is even more prolonged, indicated by a gradual increase in spontaneous and maximum firing rates up to P30. In the same cell type, broadly tuned acoustically evoked inhibition is immediately effective at hearing onset, covering the low- and high-frequency flanks of the excitatory response area. Together, these data suggest that maturation of auditory processing in the parallel ascending BC and SC streams engages distinct mechanisms at the first central synapses that may differently depend on the early auditory experience.

Soma-axon coupling configurations that enhance neuronal coincidence detection

Coincidence detector neurons transmit timing information by responding preferentially to concurrent synaptic inputs. Principal cells of the medial superior olive (MSO) in the mammalian auditory brainstem are superb coincidence detectors. They encode sound source location with high temporal precision, distinguishing submillisecond timing differences among inputs. We investigate computationally how dynamic coupling between the input region (soma and dendrite) and the spike-generating output region (axon and axon initial segment) can enhance coincidence detection in MSO neurons. To do this, we formulate a two-compartment neuron model and characterize extensively coincidence detection sensitivity throughout a parameter space of coupling configurations. We focus on the interaction between coupling configuration and two currents that provide dynamic, voltage-gated, negative feedback in subthreshold voltage range: sodium current with rapid inactivation and low-threshold potassium current, I_KLT. These currents reduce synaptic summation and can prevent spike generation unless inputs arrive with near simultaneity. We show that strong soma-to-axon coupling promotes the negative feedback effects of sodium inactivation and is, therefore, advantageous for coincidence detection. Furthermore, the feedforward combination of strong soma-to-axon coupling and weak axon-to-soma coupling enables spikes to be generated efficiently (few sodium channels needed) and with rapid recovery that enhances high-frequency coincidence detection. These observations detail the functional benefit of the strongly feedforward configuration that has been observed in physiological studies of MSO neurons. We find that I_KLT further enhances coincidence detection sensitivity, but with effects that depend on coupling configuration. For instance, in models with weak soma-to-axon and weak axon-to-soma coupling, I_KLT in the axon enhances coincidence detection more effectively than I_KLT in the soma. By using a minimal model of soma-to-axon coupling, we connect structure, dynamics, and computation. Although we consider the particular case of MSO coincidence detectors, our method for creating and exploring a parameter space of two-compartment models can be applied to other neurons.

Brain cells (neurons) are spatially extended structures. The locations at which neurons receive inputs and generate outputs are often distinct. We formulate and study a minimal mathematical model that describes the dynamical coupling between the input and output regions of a neuron. We construct our model to reflect known properties of neurons in the auditory brainstem that play an important role in our ability to locate sound sources. These neurons are known as coincidence detectors because they are most likely to respond when they receive simultaneous inputs. We use simulations to explore coincidence detection sensitivity throughout the parameter space of input-output coupling and to identify the coupling configurations that are best for neural coincidence detection. We find that strong forward coupling (from input region to output region), enhances coincidence detection sensitivity in our model and that low-threshold potassium current further improves coincidence detection. Our study is significant in that we detail how cell structure affects neuronal dynamics and, consequently, the ability of neurons to perform as temporally-precise coincidence detectors.

Enhanced Activation of HCN Channels Reduces Excitability and Spike-Timing Regularity in Maturing Vestibular Afferent Neurons

Vestibular ganglion neurons (VGNs) transmit information along parallel neuronal pathways whose signature distinction is variability in spike-timing; some fire at regular intervals while others fire at irregular intervals. The mechanisms driving timing differences are not fully understood but two opposing (but not mutually exclusive) hypotheses have emerged. In the first, regular-spiking is inversely correlated to the density of low-voltage-gated potassium currents (I _KL). In the second, regular spiking is directly correlated to the density of hyperpolarization-activated cyclic nucleotide-sensitive currents (I _H). Supporting the idea that variations in ion channel composition shape spike-timing, VGNs from the first postnatal week respond to synaptic-noise-like current injections with irregular-firing patterns if they have I _KL and with more regular firing patterns if they do not. However, in vitro firing patterns are not as regular as those in vivo Here we considered whether highly-regular spiking requires I _H currents and whether this dependence emerges later in development after channel expression matures. We recorded from rat VGN somata of either sex aged postnatal day (P)9-P21. Counter to expectation, in vitro firing patterns were less diverse, more transient-spiking, and more irregular at older ages than at younger ages. Resting potentials hyperpolarized and resting conductance increased, consistent with developmental upregulation of I _KL Activation of I _H (by increasing intracellular cAMP) increased spike rates but not spike-timing regularity. In a model, we found that activating I _H counter-intuitively suppressed regularity by recruiting I _KL Developmental upregulation in I _KL appears to overwhelm I _H These results counter previous hypotheses about how I _H shapes vestibular afferent responses.SIGNIFICANCE STATEMENT Vestibular sensory information is conveyed on parallel neuronal pathways with irregularly-firing neurons encoding information using a temporal code and regularly-firing neurons using a rate code. This is a striking example of spike-timing statistics influencing information coding. Previous studies from immature vestibular ganglion neurons (VGNs) identified hyperpolarization-activated mixed cationic currents (I _H) as driving highly-regular spiking and proposed that this influence grows with the current during maturation. We found that I _H becomes less influential, likely because maturing VGNs also acquire low-voltage-gated potassium currents (I _KL), whose inhibitory influence opposes I _H Because efferent activity can partly close I _KL, VGN firing patterns may become more receptive to extrinsic control. Spike-timing regularity likely relies on dynamic ion channel properties and complementary specializations in synaptic connectivity.

A low-threshold potassium current enhances sparseness and reliability in a model of avian auditory cortex

Birdsong is a complex vocal communication signal, and like humans, birds need to discriminate between similar sequences of sound with different meanings. The caudal mesopallium (CM) is a cortical-level auditory area implicated in song discrimination. CM neurons respond sparsely to conspecific song and are tolerant of production variability. Intracellular recordings in CM have identified a diversity of intrinsic membrane dynamics, which could contribute to the emergence of these higher-order functional properties. We investigated this hypothesis using a novel linear-dynamical cascade model that incorporated detailed biophysical dynamics to simulate auditory responses to birdsong. Neuron models that included a low-threshold potassium current present in a subset of CM neurons showed increased selectivity and coding efficiency relative to models without this current. These results demonstrate the impact of intrinsic dynamics on sensory coding and the importance of including the biophysical characteristics of neural populations in simulation studies.

Maintaining a stable mental representation of an object is an important task for sensory systems, requiring both recognizing the features required for identification and ignoring incidental changes in its presentation. The prevailing explanation for these processes emphasizes precise sets of connections between neurons that capture only the essential features of an object. However, the intrinsic dynamics of the neurons themselves, which determine how these inputs are transformed into spiking outputs, may also contribute to the neural computations underlying object recognition. To understand how intrinsic dynamics contribute to sensory coding, we constructed a computational model capable of simulating a neural response to an auditory stimulus using a detailed description of different intrinsic dynamics in a higher-order avian auditory area. The results of our simulation showed that intrinsic dynamics can have a profound effect on processes underlying object recognition. These findings challenge the view that patterns of connectivity alone account for the emergence of stable object representations and encourage greater consideration of the functional implications of the diversity of neurons in the brain.

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.

Role of Oscillations in Auditory Temporal Processing: A General Model for Temporal Processing of Sensory Information in the Brain?

We review the role of oscillations in the brain and in the auditory system showing that the ability of humans to distinguish changes in pitch can be explained as a precise analysis of temporal information in auditory signals by neural oscillations. The connections between auditory brain stem chopper neurons construct neural oscillators, which discharge spikes at various constant intervals that are integer multiples of 0.4 ms, contributing to the temporal processing of auditory cochlear output. This is subsequently spatially mapped in the inferior colliculus. Electrophysiological measurements of auditory chopper neurons in different species show oscillations with periods which are integer multiples of 0.4 ms. The constant intervals of 0.4 ms can be attributed to the smallest synaptic delay between interconnected simulated chopper neurons. We also note the patterns of similarities between microcircuits in the brain stem and other parts of the brain (e.g., the pallidum, reticular formation, locus coeruleus, oculomotor nuclei, limbic system, amygdala, hippocampus, basal ganglia and substantia nigra), dedicated to the processing of temporal information. Similarities in microcircuits across the brain reflect the importance of one of the key mechanisms in the information processing in the brain, namely the temporal coupling of different neural events via coincidence detection.

Biological plausibility and stochasticity in scalable VO2 active memristor neurons

Neuromorphic networks of artificial neurons and synapses can solve computationally hard problems with energy efficiencies unattainable for von Neumann architectures. For image processing, silicon neuromorphic processors outperform graphic processing units in energy efficiency by a large margin, but deliver much lower chip-scale throughput. The performance-efficiency dilemma for silicon processors may not be overcome by Moore's law scaling of silicon transistors. Scalable and biomimetic active memristor neurons and passive memristor synapses form a self-sufficient basis for a transistorless neural network. However, previous demonstrations of memristor neurons only showed simple integrate-and-fire behaviors and did not reveal the rich dynamics and computational complexity of biological neurons. Here we report that neurons built with nanoscale vanadium dioxide active memristors possess all three classes of excitability and most of the known biological neuronal dynamics, and are intrinsically stochastic. With the favorable size and power scaling, there is a path toward an all-memristor neuromorphic cortical computer.

Diverse Intrinsic Properties Shape Functional Phenotype of Low-Frequency Neurons in the Auditory Brainstem

In the auditory system, tonotopy is the spatial arrangement of where sounds of different frequencies are processed. Defined by the organization of neurons and their inputs, tonotopy emphasizes distinctions in neuronal structure and function across topographic gradients and is a common feature shared among vertebrates. In this study we characterized action potential firing patterns and ion channel properties from neurons located in the extremely low-frequency region of the chicken nucleus magnocellularis (NM), an auditory brainstem structure. We found that NM neurons responsible for encoding the lowest sound frequencies (termed NMc neurons) have enhanced excitability and fired bursts of action potentials to sinusoidal inputs ≤10 Hz; a distinct firing pattern compared to higher-frequency neurons. This response property was due to lower amounts of voltage dependent potassium (KV) conductances, unique combination of KV subunits and specialized sodium (NaV) channel properties. Particularly, NMc neurons had significantly lower KV1 and KV3 currents, but higher KV2 current. NMc neurons also showed larger and faster transient NaV current (INaT) with different voltage dependence of inactivation from higher-frequency neurons. In contrast, significantly smaller resurgent sodium current (INaR) was present in NMc with kinetics and voltage dependence that differed from higher-frequency neurons. Immunohistochemistry showed expression of NaV1.6 channel subtypes across the tonotopic axis. However, various immunoreactive patterns were observed between regions, likely underlying some tonotopic differences in INaT and INaR. Finally, using pharmacology and computational modeling, we concluded that KV3, KV2 channels and INaR work synergistically to regulate burst firing in NMc.

Functional roles of Kv1-mediated currents in genetically identified subtypes of pyramidal neurons in layer 5 of mouse somatosensory cortex

We used voltage-clamp recordings from somatic outside-out macropatches to determine the amplitude and biophysical properties of putative Kv1-mediated currents in layer 5 pyramidal neurons (PNs) from mice expressing EGFP under the control of promoters for etv1 or glt. We then used whole cell current-clamp recordings and Kv1-specific peptide blockers to test the hypothesis that Kv1 channels differentially regulate action potential (AP) voltage threshold, repolarization rate, and width as well as rheobase and repetitive firing in these two PN types. We found that Kv1-mediated currents make up a similar percentage of whole cell K⁺ current in both cell types, and only minor biophysical differences were observed between PN types or between currents sensitive to different Kv1 blockers. Putative Kv1 currents contributed to AP voltage threshold in both PN types, but AP width and rate of repolarization were only affected in etv1 PNs. Kv1 currents regulate rheobase, delay to the first AP, and firing rate similarly in both cell types, but the frequency-current slope was much more sensitive to Kv1 block in etv1 PNs. In both cell types, Kv1 block shifted the current required to elicit an onset doublet of action potentials to lower currents. Spike frequency adaptation was also affected differently by Kv1 block in the two PN types. Thus, despite similar expression levels and minimal differences in biophysical properties, Kv1 channels differentially regulate APs and repetitive firing in etv1 and glt PNs. This may reflect differences in subcellular localization of channel subtypes or differences in the other K⁺ channels expressed. NEW & NOTEWORTHY In two types of genetically identified layer 5 pyramidal neurons, α-dendrotoxin blocked approximately all of the putative Kv1 current (on average). We used outside-out macropatches and whole cell recordings at 33°C to show that despite similar expression levels and minimal differences in biophysical properties, Kv1 channels differentially regulate action potentials and repetitive firing in etv1 and glt pyramidal neurons. This may reflect differences in subcellular localization of channel subtypes or differences in the other K⁺ channels expressed.

Extraction of Inter-Aural Time Differences Using a Spiking Neuron Network Model of the Medial Superior Olive

The mammalian auditory system is able to extract temporal and spectral features from sound signals at the two ears. One important cue for localization of low-frequency sound sources in the horizontal plane are inter-aural time differences (ITDs) which are first analyzed in the medial superior olive (MSO) in the brainstem. Neural recordings of ITD tuning curves at various stages along the auditory pathway suggest that ITDs in the mammalian brainstem are not represented in form of a Jeffress-type place code. An alternative is the hemispheric opponent-channel code, according to which ITDs are encoded as the difference in the responses of the MSO nuclei in the two hemispheres. In this study, we present a physiologically-plausible, spiking neuron network model of the mammalian MSO circuit and apply two different methods of extracting ITDs from arbitrary sound signals. The network model is driven by a functional model of the auditory periphery and physiological models of the cochlear nucleus and the MSO. Using a linear opponent-channel decoder, we show that the network is able to detect changes in ITD with a precision down to 10 μs and that the sensitivity of the decoder depends on the slope of the ITD-rate functions. A second approach uses an artificial neuronal network to predict ITDs directly from the spiking output of the MSO and ANF model. Using this predictor, we show that the MSO-network is able to reliably encode static and time-dependent ITDs over a large frequency range, also for complex signals like speech.

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

Models of the auditory brainstem have been an invaluable tool for testing hypotheses about auditory information processing and for highlighting the most important gaps in the experimental literature. Due to the complexity of the auditory brainstem, and indeed most brain circuits, the dynamic behavior of the system may be difficult to predict without a detailed, biologically realistic computational model. Despite the sensitivity of models to their exact construction and parameters, most prior models of the cochlear nucleus have incorporated only a small subset of the known biological properties. This confounds the interpretation of modelling results and also limits the potential future uses of these models, which require a large effort to develop. To address these issues, we have developed a general purpose, biophysically detailed model of the cochlear nucleus for use both in testing hypotheses about cochlear nucleus function and also as an input to models of downstream auditory nuclei. The model implements conductance-based Hodgkin-Huxley representations of cells using a Python-based interface to the NEURON simulator. Our model incorporates most of the quantitatively characterized intrinsic cell properties, synaptic properties, and connectivity available in the literature, and also aims to reproduce the known response properties of the canonical cochlear nucleus cell types. Although we currently lack the empirical data to completely constrain this model, our intent is for the model to continue to incorporate new experimental results as they become available.