Abstract
GABA and glycine are the major inhibitory transmitters that attune neuronal activity in the CNS of mammals. The respective transmitters are mostly spatially separated, that is, synaptic inhibition in the forebrain areas is mediated by GABA, whereas glycine is predominantly used in the brainstem. Accordingly, inhibition in auditory brainstem circuits is largely mediated by glycine, but there are few auditory synapses using both transmitters in maturity. Little is known about physiological advantages of such a two-transmitter inhibitory mechanism. We explored the benefit of engaging both glycine and GABA with inhibition at the endbulb of Held-spherical bushy cell synapse in the auditory brainstem of juvenile Mongolian gerbils. This model synapse enables selective in vivo activation of excitatory and inhibitory neuronal inputs through systemic sound stimulation and precise analysis of the input (endbulb of Held) output (spherical bushy cell) function. The combination of in vivo and slice electrophysiology revealed that the dynamic AP inhibition in spherical bushy cells closely matches the inhibitory conductance profile determined by the glycine-R and GABAA-R. The slow and potent glycinergic component dominates the inhibitory conductance, thereby primarily accounting for its high-pass filter properties. GABAergic transmission enhances the inhibitory strength and shapes its duration in an activity-dependent manner, thus increasing the inhibitory potency to suppress the excitation through the endbulb of Held. Finally, in silico modeling provides a strong link between in vivo and slice data by simulating the interactions between the endbulb- and the synergistic glycine-GABA-conductances during in vivo-like spontaneous and sound evoked activities.
- AP inhibition
- auditory signal processing
- GABAA/glycine receptors
- in vivo electrophysiology
- modeling
- slice electrophysiology
Introduction
Inhibitory synapses in the adult mammalian CNS provide a functional balance of neuronal activity through release of GABA or glycine. GABA mediates inhibition mainly in the cortex and cerebellum, whereas fast glycinergic transmission is predominantly restricted to the brainstem and spinal cord (Fritschy et al., 1994; Legendre, 2001; Huang et al., 2007; Lehmann et al., 2012). Similar patterns of respective neurotransmitters were also described for the central auditory system (Prieto et al., 1994; Kurt et al., 2006; Caspary et al., 2008; Friauf et al., 2011; Sanes and Kotak, 2011).
However, despite this general spatial differentiation, a number of special central synapses engage glycine and GABA signaling in terms of synaptic vesicle content, corelease, and activation of postsynaptic receptors (Triller et al., 1987; Burger et al., 1991; Bohlhalter et al., 1994; Todd et al., 1996; Dumba et al., 1998; Jonas et al., 1998). In the auditory system, GABA is associated with development when it evokes a transient membrane depolarization and ensuing calcium signaling (Sanes and Friauf, 2000; Kullmann and Kandler, 2001; Kullmann et al., 2002; Milenkovic et al., 2007; Witte et al., 2014). Inhibitory terminals in the lateral superior olive change the transmitter phenotype from GABA to glycine after hearing onset (Kotak et al., 1998; Nabekura et al., 2004; Kim and Kandler, 2010). Yet, persistent corelease of both transmitters was demonstrated at some brainstem synapses (Awatramani et al., 2005; Apostolides and Trussell, 2013). Spherical bushy cells (SBCs) in the anteroventral cochlear nucleus (AVCN) receive two types of acoustically triggered inputs: Excitation through the endbulbs of Held, and inhibitory inputs from AVCN D-stellate cells (Smith and Rhode, 1989; Campagnola and Manis, 2014) and dorsal cochlear nucleus tuberculoventral cells (Wickesberg and Oertel, 1990; Saint Marie et al., 1991). Slice experiments revealed a conspicuous inhibitory role of glycinergic signaling in SBCs (Harty and Manis, 1996; Lim et al., 1999, 2000; Xie and Manis, 2013), whereas in vivo recordings demonstrated a contribution of both glycine and GABA to attunement of firing properties (Caspary et al., 1994; Gai and Carney, 2008; Dehmel et al., 2010). Although synaptic inputs to SBCs have been extensively characterized in slice preparations from mice and rats, the function of GABAA-R remained elusive (but for GABAB-R, see Lim et al., 2000; Chanda and Xu-Friedman, 2010). Moreover, a wealth of in vivo studies have been conducted in the AVCN of Mongolian gerbils, but matching characterization of inhibitory conductances is still lacking.
Here, we explored the interaction of glycine-and GABAergic inhibition at the endbulb of Held synapse using selective activation of inhibitory neuronal inputs through systemic sound stimulation. By performing juxtacellular recordings in vivo and whole-cell measurements in acute slices, we identified the mechanisms of SBC's coinhibition in juvenile Mongolian gerbils. The efficacy of glycine-GABA transmission was assessed over a range of physiologically relevant conductance interactions by implementing an in silico model. Our results provide evidence for synergistic GABA enhancement of predominantly glycinergic inhibition. The slow kinetics of glycine-GABA inhibitory currents shows a frequency-dependent buildup, thereby providing efficient dynamic control of AP firing.
Materials and Methods
The experimental procedures were approved by the Saxonian district Government Leipzig (T 93/11, T 84/12, T 67/13 and TVV 06/09) and conducted according to the European Communities Council Directive (86/609/EEC).
In vivo experiments.
Recordings were performed in 9 Mongolian gerbils (Meriones unguiculatus) of either sex aged postnatal days 25–33 (P25–P33). Data were acquired on juvenile animals to circumvent the prolonged postnatal development of inhibitory transmission in auditory brainstem nuclei (Awatramani et al., 2005; Luján et al., 2008; Friauf et al., 2011; Milenkovic and Rübsamen, 2011). Recordings were performed in a sound-attenuated chamber (Type 400, Industrial Acoustic) while the animal was stabilized in a custom-made stereotaxic apparatus positioned on a vibration-isolated table.
Surgical preparation.
Animals were anesthetized with an initial intraperitoneal injection of a mixture of ketamine hydrochloride (0.18 mg/g body weight; Ketamin, Ratiopharm) and xylazine hydrochloride (7 μg/g body weight; Rompun, Bayer). Anesthesia was maintained by additional subcutaneous application of one-third of the initial dose, approximately every 60–90 min during recording sessions. Gerbils were fixed in the stereotaxic device by a metal bolt glued to the skull exposed around the bregma point. The recording electrode was inserted through a hole (diameter 0.5 mm) drilled 1.8–2 mm lateral to the midline and 2 mm caudal to the lambdoid suture. The reference electrode was placed in the superficial cerebellum through a second hole (diameter 0.5 mm) located on the midline. The AVCN was approached dorsally by tilting the animal at 12°–18° to the midsagittal plane. Body temperature was kept between 36.5°C and 37.5°C with a temperature-controlled heating pad.
Acoustic stimulation.
Auditory stimuli were digitally generated using custom-written MATLAB software (MathWorks). The stimuli were transferred to a D/A converter (RP2.1 real-time processor, 97.7 kHz sampling rate, 24 bit, Tucker-Davis Technologies) and delivered through a custom-made earphone (acoustic transducer: DT 770 pro, Beyer Dynamics) fitted with a plastic tube (length 35 mm, diameter 5 mm) positioned into the outer ear canal at a distance of ∼4 mm to the eardrum.
The unit's frequency response area was obtained by a pseudorandom presentation of pure tone pulses (100 ms duration, 5 ms rise-fall time, 200 ms interstimulus interval) derived from a matrix of predefined frequency/intensity pairs (20 frequencies on a log scale, 10 intensity levels on linear scale, 3–5 repetitions). Based on the unit's frequency response area, its characteristic frequency (CF, frequency to which neuron is most sensitive), the excitatory and inhibitory response threshold and the area of the high-frequency inhibitory sideband were assessed. For in-depth analyses, neuronal discharge activities were also recorded during repetitive stimulation (200–300 times) of single pure tones (1) at the unit's CF and (2) within the inhibitory sideband (100 ms duration, 5 ms rise-fall time, 300 ms interstimulus interval). Sound pressure level was set 20 dB above the respective thresholds.
Single-unit recording.
Considering the tonotopic organization of the AVCN (Kopp-Scheinpflug et al., 2002; Dehmel et al., 2010), the multiunit activity at the rostral pole was initially assessed with low impedance glass micropipettes (GB150F-10, Science Products, 1–5 mΩ filled with 3M KCl) to narrow down the target area (frequencies <5 kHz). Then, juxtacellular single-unit recordings (electrodes 7–10 mΩ) were performed on spherical bushy cells, which were identified according to the characteristic waveform (Pfeiffer, 1966; Winter and Palmer, 1990; Englitz et al., 2009; Typlt et al., 2010), and the primary-like peristimulus time histogram pattern (Blackburn and Sachs, 1989).
Recorded voltage signals were amplified (Neuroprobe 1600, A-M Systems), and digitized at a sampling rate of 97.7 kHz (24 bit, RP2.1, Tucker-Davis Technologies). Signals were bandpass filtered between 50 Hz and 5 kHz using a zero-phase forward and reverse digital IIR filter and stored for offline analysis using custom-written MATLAB software (version 8.1, MathWorks). Recordings were selected for further analysis according to three criteria: (1) SD (AP height)/mean (AP height) <20% (10.4 ± 3.9%, n = 16), (2) uniform waveforms, and (3) signal-to-noise ratio >12 (16.1 ± 5.1, n = 16).
Spontaneous and stimulus evoked AP firing in SBCs is characterized by the typical voltage signature composed of two or three components (P-A-B; see Fig. 1A), indicating the discharge of the presynaptic endbulb (P), the postsynaptic EPSP (A), and the postsynaptic AP (B) (Englitz et al., 2009; Typlt et al., 2010). The components were distinguished by means of hierarchical clustering after dimension reduction using principal component analysis. For the waveforms containing all three components, the inflection point at the rising flank of the voltage signal was determined as the EPSP position by calculating the local minimum of the first derivative preceding the AP peak. For signals with an EPSP separated from the AP, the first local maximum preceding the AP was considered the EPSP peak. Time between the EPSP and the peak of the AP indicated the EPSP-to-AP transition time. Average rising slope of the EPSP was calculated between 20% and 80% of the EPSP amplitude.
For histological verification of recording site, Fluorogold was iontophoretically injected at the end of recording session (4 μA for 7 min). The animal was perfused 4–6 h thereafter with 0.9% NaCl solution followed by 5% PFA. After overnight postfixation, the brain was cut on a vibratome (Microm HM 650) and the tissue sections (100 μm thick) were visualized under the confocal laser scanning microscope (TCS PS5, Leica).
Slice preparation.
Coronal slices (170 μm) containing rostral AVCN were cut from P22–P33 gerbils of either sex (Milenkovic et al., 2007; Dietz et al., 2012). The rationale was to perform the experiments on a rather developed system matching the acquired in vivo data. The brainstem was sliced with a vibratome, in low-calcium ACSF solution containing the following (in mm): 125 NaCl, 2.5 KCl, 0.1 CaCl2, 3 MgCl2, 1.25 NaH2PO4, 25 NaHCO3, 25 glucose, 2 sodium pyruvate, 3 myoinositol, 0.5 ascorbic acid, continuously bubbled with 5% CO2 and 95% O2, pH 7.4. Slices were incubated in the standard recording solution (ACSF same as for slicing, except CaCl2 and MgCl2 were changed to 2 and 1 mm, respectively) for 30 min at 37°C and stored at room temperature until recording. Experiments were made at nearly physiological temperature (33 ± 0.5°C).
Perforated-patch and whole-cell recordings.
Gramicidin perforated-patch recordings on large SBCs were acquired as previously described (Witte et al., 2014). Patch pipettes were pulled with Narishige PC-10 vertical puller from filamented borosilicate glass capillaries (Science Products) to have resistances of 5–6 mΩ when filled with a K+-gluconate-based internal solution. The potassium-based solution was used to avoid the potential effect of Cs+ ions on the transport activity of KCC2 (Kakazu et al., 1999). Pipettes were tip-filled as follows (mm): 97.5 potassium-gluconate, 32.5 KCl, 1 MgCl2, 10 HEPES, 5 EGTA, 10 HEPES, pH 7.34 with KOH (290 mOsm adjusted with sucrose). The remainder of the pipette was back-filled with the same solution including gramicidin (30 μg/ml gramicidin A, Sigma) and 25 μm ATTO 488. The latter was used to confirm that the perforated patch was not ruptured. After the experiment, the rupture of the perforated-patch yielded depolarizing IPSPs (due to 34.5 mm Cl−), thus confirming that the gramicidin ionophores were impermeable for Cl− and that recordings were not affected by the [Clpip−]. The progress of perforation was evaluated by monitoring the steady-state current responses to a −5 mV voltage command. The series resistance typically reached a steady level (mean Rs = 18.3 ± 0.7 mΩ, bridge = 21.9 ± 0.7 mΩ, n = 5) within 30–40 min after the giga-seal formation and the experiments were started thereafter.
The recordings of synaptically evoked IPSCs (eIPSCs) on SBCs were conducted as described previously (Dietz et al., 2012). The pipettes had resistances of 4–5 mΩ when filled with (mm) as follows: 140 CsMeSO3, 20 TEA-Cl, 3.3 MgCl2, 10 HEPES, 0.1 EGTA, 5 QX-314-Cl, 5 phosphocreatine, 2 ATP disodium salt, 0.3 GTP disodium salt, and 0.2% biocytin (pH 7.3 with CsOH). For current-clamp recordings of eIPSPs the pipette solution contained the following (mm): 127 potassium-gluconate, 3 KCl, 1 MgCl2, 10 HEPES, 0.2 EGTA, 5 phosphocreatine, 2 ATP disodium salt, 0.3 GTP disodium salt, pH adjusted to 7.3 with KOH. IPSCs and IPSPs were evoked by electrical stimulation of afferent fibers through a bipolar theta glass electrode (Sutter instruments, tip Ø 5 μm) filled with bath solution and placed 30–60 μm adjacent to the cell. The stimulus intensity was slowly increased until set at the lowest value that reliably evoked IPSCs of stable amplitudes within a train. Pulse stimuli (100 μs) were generated by a stimulator (Master 8) and delivered via an isolated stimulus unit (AMPI Iso-flex) to evoke either single events or train-responses at different frequencies. Current-clamp responses were recorded at −62 to −66 mV, approximately corresponding to the resting membrane potential of bushy cells (McGinley and Oertel, 2006; Price and Trussell, 2006; Milenkovic et al., 2007), whereas voltage-clamp measurements were done from Vhold = −71 mV. Pharmacological inhibition of glutamate (50 μm AP-5, 10 μm NBQX) and GABAB receptors (3 μm CGP55845) was performed in all experiments. Offline correction of voltages was done for junction potentials of 11 mV (VC) and 14 mV (CC).
Glycine and GABA (500 μm) were prepared in HCO3−-containing bath solution and pressure applied over the soma of recorded neuron using a Picospritzer (General Valve). As in all other experiments, extracellular solution was supplemented with the GABAB antagonist CGP 55845. The constant stimulation conditions were assured by controlling the pipette diameter, application pressure and duration, and distance from the cell (3 μm, 5 psi, 5 ms, 10 μm, respectively). The perfusion was turned off just before each puff application to avoid unequal dilution of the agonist. Significance of the responses was determined by using the z test, i.e., level of acceptance was set at z > 3.3 (Vm) or z < −3.3 (Im), which corresponds to p < 0.001 [z = (A-BL)/SDBL, with A being the maximal amplitude of the response, BL the mean of the baseline (2 s before stimulation), and SDBL the SD of the baseline]. Pressure ejection of ACSF under the same condition evoked no response (mean ACSF response = −0.34 ± 0.2 pA, z = −0.21 ± 0.13, p > 0.73, n = 4), whereas glycine and GABA evoked significant membrane currents in the same neurons.
The recordings were acquired using a Multiclamp 700B amplifier (Molecular Devices). In voltage-clamp recordings, series resistance was compensated by 50% to remaining Rs of 4–7 mΩ. Bridge balance and pipette capacitance neutralization were adjusted throughout the experiment in current-clamp recordings. Recorded signals were digitized at 20/50 kHz and filtered with a 6 kHz Bessel low-pass filter. Data were examined with pClamp 10 software (Molecular Devices) followed by the analysis with custom-written MATLAB routines.
Data analysis.
Mean peak amplitudes, 10–90% rise times, and decay time constants were analyzed from averaged traces (>7 repetitions) with custom-written MATLAB routines. IPSC decay phase was fitted with biexponential function (between 95% and 5% of the peak amplitude). The weighted τ decay for biexponential fitting was calculated as follows: τ = (Afast × τfast + Aslow × τslow)/(Afast + Aslow), where Afast and Aslow are amplitudes at t = 0 and τfast and τslow are the fast and slow time constants, respectively. Calculation of IPSC conductance was based on the Cl−-driving force of 20 mV determined by the Vhold−EIPSC = −71 to −50.8 mV (n = 6). EIPSC was controlled at the beginning of each experiment.
The average depolarization slope was determined from the first derivative, between 20% and 80% of the interval from the time point of the current injection to the AP initiation point (Depolsucc), or alternatively to the maximum amplitude of subthreshold depolarization (Depolfail). The AP initiation was computed from the maximum of the second derivative, representing the time point of the fastest membrane potential increase, estimated as the AP onset.
Post hoc labeling of biocytin-filled neurons with Cy2-conjugated streptavidin was used to morphologically confirm recordings from SBCs (Milenkovic et al., 2009; Dietz et al., 2012). Images were generated with a confocal laser scanning microscope (TCS PS5, Leica).
Statistics.
Shapiro–Wilk test (Shapiro and Wilk, 1965) for platykurtic samples, or Shapiro–Francia test (Shapiro and Francia, 1972) for leptokurtic samples were used to determine whether the null hypothesis of Gaussian distribution is a reasonable assumption regarding the data distribution. Assumption of homoscedasticity was tested using Levene's test (Levene, 1960), and datasets were compared with the appropriate t test or ANOVA with post hoc pairwise comparisons (Holm–Sidak test) (Sigma Plot 11, Systat Software). Two-way ANOVA was used to compare the data in Figure 5Ci. Control values before superfusion of strychnine and control values before superfusion of SR95531 were compared for all stimulation frequencies. As there were no differences for a given frequency before drug superfusion, the data were pooled and presented as single control histogram bar for each frequency. Thereafter, the main effect of drug perfusion and possible interaction with input frequency were tested with two-way ANOVA followed by pairwise comparisons (Holm–Sidak test).
Repeated-measures (RM) ANOVA was applied to test for effects of drugs (i.e., before-drug, drug, and after-drug conditions). Mann–Whitney U test was applied in one case where the data showed skewed non-Gaussian distribution (see Fig. 2A). Data are reported regarding the distribution as mean ± SEM or median (1 quartile, 3 quartile), unless otherwise noted.
Modeling.
All simulations were performed with NEURON (Hines and Carnevale, 1997), imported as the module pyNEURON (implemented by Uri Cohen) in Python 2.7.3 under Win7–64bit. The biophysical properties of the SBC model were as previously published (Kuenzel et al., 2011) and in accordance with the model by Rothman et al. (1993). An implementation of the sodium conductance (Rothman et al., 1993) was kindly provided by Marek Rudnicki, TU Munich. Other ion-channel models were identical to those described previously (Rothman and Manis, 2003). The model parameters are represented in detail in Table 1. The present model also included (1) a somatic compartment, (2) an axon hillock/first segment compartment containing all voltage-activated sodium conductances, and (3) a stretch of axon.
Two conductance point processes were attached to the somatic compartment to act as the excitatory endbulb terminal (Erev = 0 mV) and as the combined inhibitory input (Erev = −75 mV). Conductance traces were pregenerated (see below) and applied to the point processes at the rate of the simulation (dt = 10 μs). Conductance templates were fitted to match the rise and decay time constants of recorded waveforms. Excitatory conductances corresponded to the EPSPs recorded in vivo (compare Kuenzel et al., 2011), inhibitory conductances were fitted to IPSC data from slice recordings presented in the study at hand. Temperature differences between the different sources of templates were not taken into account for the simulations.
Rate-dependent plasticity of the inhibitory synaptic mechanism was implemented as a dynamic state-variable that relaxed back to the resting state with double-exponential functions (compare Varela et al., 1997). Exponential parameters of the synaptic plasticity model were as follows: Afast 0.79, Aslow 0.99, τfast 31.1 ms, τslow 316 ms. Parameters of the decay τ change were as follows: Afast 8.7, Aslow 2.8, τfast 16.9 ms, τslow 151 ms, τmax 80 ms. Exponential parameters were extracted by fitting the model results to recorded data. The conductance templates were thereafter convolved with pregenerated spike times to create the conductance traces. Spike times were either single interactions or pulse trains at various frequencies, with the inhibitory input leading excitation by 1 ms. Such a latency was chosen to match the slice experiments investigating the inhibitory potency to suppress spikes evoked by current injection (compare Fig. 2 and 6). Given the nonprimary source of inhibition in the cochlear nucleus, it is conceivable that at some time point the excitatory input will precede the effect of inhibition. However, considering slow decay times of inhibition measured here, at input frequencies in physiological range it is likely that the inhibition overlaps the excitatory conductance. In all other models (see Figs. 8B and 9), inhibition was trailing excitation by 1 ms, to mimic an additional synapse in the inhibitory pathway reaching the SBC. This is based on studies showing that neurons in the deep dorsal cochlear nucleus, which are directly innervated through AN inputs, are probably the major source of inhibition to SBCs (Wickesberg and Oertel, 1990; Xie and Manis, 2013; Campagnola and Manis, 2014). To assess the model performance with in vivo like inputs, spontaneous (silence) and sound-driven (white-noise) spike trains were generated with the help of a simple AN model. This AN model was implemented as a gamma-tone filter bank driving leaky integrate-and-fire neurons with noise and refractoriness, using the tools provided by the spiking neural-network simulator “Brian” (Goodman and Brette, 2008). In the absence of sound stimuli, the AN model generates random spike trains with ISI that essentially follow a shifted-exponential distribution. During simulated sound stimulation, the model generates phase-locked, AN-like spike trains.
In the model, we also implemented a more realistic representation of the excitatory endbulb synapse. Numerous in vitro studies showed prominent short-term synaptic depression, occurring also in preparations from older animals (Wang and Manis, 2008; Yang and Xu-Friedman, 2009; Wang et al., 2010). Analogous to our approach for the inhibitory synaptic plasticity, a phenomenological model was fitted to pulse-train data presented in respective studies. Data for 100, 200, 300, and 333 Hz trains were used and the mean results of these fits were taken as model parameters. A depression state variable was reduced by every event and relaxed back to rest with a double exponential function (Afast 0.75, Aslow 0.94, τfast 17.2 ms, τslow 57.0 ms). The maximal endbulb conductance was multiplied with the state variable to derive the excitatory conductance active at a given time. Maximal endbulb conductance of 75 nS was used to counteract tonic depression by spontaneous activity in the model. As shown previously by Pliss et al. (2009), endbulb conductances measured in vitro can span a large range of values, with 75 nS being toward the lower end.
As the short-term depression at the endbulb synapse may not play such a prominent role in the adult auditory brainstem in vivo (Borst, 2010; Crins et al., 2011; Kuenzel et al., 2011), a nondepressing stochastic endbulb model synapse was also included. Here, the EPSG amplitude for every event was drawn from a normal distribution with a SD of ±0.2 around the mean, which led to an average EPSG of 48.9 ± 9.8 nS.
Interaction of the excitatory and inhibitory conductance in the SBC model was scrutinized by simulating the membrane potential response of the recorded SBC. All experiments were simulated at 35°C.
Results
Inhibition attunes the in vivo input–output rate at the endbulb of Held synapse
To evaluate the impact of acoustically evoked inhibition on AP generation in SBCs, we used pure tone stimulation within the unit's high-frequency inhibitory sideband. The presently recorded units were located in the rostral AVCN, had low CFs (median 2.0 [first quartile 1.39, third quartile 2.24] kHz, n = 16), primary-like response patterns, and three-component complex spike waveforms. Together, these characteristics identify them as large SBCs (Tolbert et al., 1982; Rouiller and Ryugo, 1984; Hackney et al., 1990; Ostapoff et al., 1994; Bazwinsky et al., 2008; Englitz et al., 2009). The average spontaneous firing rate was 53.5 ± 17.2 Hz and thresholds ranged from −8.1 to 63.3 dB SPL (median 12 [first quartile 5.9, third quartile 30.2] dB).
The juxtacellular recorded voltage signals from each SBC were divided into two categories, based on principal component analysis complemented by hierarchical clustering (see also Englitz et al., 2009, Typlt et al., 2010): (1) three components signals (P-A-B), i.e., signals containing a discernable postsynaptic AP (EPSPsucc); and (2) signals composed of the components ‘P’ and ‘A,’ but lacking ‘B,’ i.e., signals representing only the endbulb discharge and the postsynaptic EPSP (EPSPfail) (Fig. 1A).
Typically, the acoustically evoked response profiles of SBCs show inhibitory sidebands flanking the excitatory response areas (Fig. 1B). The respective excitatory and inhibitory domains were determined as frequency-intensity ranges in which the discharge rates were significantly increased/reduced compared with the spontaneous firing rate (p < 0.01, two-sample t test). The inhibitory sidebands spanned the frequency range 0.75–1.9 octaves above unit's CF (mean ± SD = 1.27 ± 0.4 O, n = 16). Acoustically evoked inhibition was elicited by a tone burst stimulation at 4.8 [first quartile 4, third quartile 4.9] kHz and 44 ± 15 dB SPL, n = 16 (Fig. 1B). Under these conditions, the incidence of failures of postsynaptic spikes was significantly increased (EPSPfail) with respect to spontaneous activity (Fig. 1C,D; EPSPfail/(EPSPfail + EPSPsucc) = 0.61 ± 0.19 vs 0.27 ± 0.13, p < 0.001, paired t test). Notably, such acoustic stimulation did not change the overall input to the SBC, as indicated by the constant overall EPSP rate. This excludes the possibility that a smaller incidence of postsynaptic APs is due to a reduced ANF input (Fig. 1D; EPSP = EPSPsucc + EPSPfail, 10–90 ms vs 150–230 ms, p = 0.34, paired t test) (further elaborated in Fig. 8A). These in vivo data demonstrate that acoustically evoked inhibition can control the fidelity of signal transmission at the endbulb of Held synapse by engaging postsynaptic mechanisms while leaving the presynaptic calyceal input unaffected.
Synaptic inhibition prevents the AP generation
Generally, GABAA-R and glycine-R mediated inhibition in mammalian neurons is achieved by both shunting (i.e., reduction of the amplitude of EPSPs due to a local reduction of the input resistance in the vicinity of excitatory synapses) and membrane hyperpolarization (i.e., moving the membrane potential away from the action potential threshold). Accordingly, increased occurrence of postsynaptic AP failures during acoustic stimulation might be attributable to the inhibitory transmission interfering with conductance generated by the calyceal excitatory inputs. To assess such potential interaction, the waveforms of SBC signals recorded in vivo during acoustic stimulation within the inhibitory sidebands were analyzed and compared with the respective signals recorded during spontaneous activity (Fig. 2A). Under both conditions, the EPSP-to-AP transition times were negatively correlated with the EPSP rising slopes, i.e., the more shallow the slopes the longer the EPSP-to-AP transition times (spontaneous rs = −0.56, p < 0.001; inh. rs = −0.55, p < 0.001) (Fig. 2B, left). Yet, the EPSP-to-AP transition time was significantly longer and rising slopes flattened during stimulation in the units' acoustic sidebands (Fig. 2A,B; delay: p < 0.01 in 16 of 16 neurons, slope: p < 0.05 in 10 of 16 neurons, Mann–Whitney U test). From these results, it can be concluded that the inhibition reduces the slope steepness of the membrane potential depolarization, which was shown to be crucial for successful AP generation in SBCs (McGinley and Oertel, 2006).
Postnatal development of the chloride homeostasis in SBCs renders inhibitory transmission hyperpolarizing before the onset of acoustic information processing (Milenkovic et al., 2007). To assess the potency of hyperpolarizing inhibition to impair the AP generation, gramicidin perforated-patch recordings were conducted in brainstem slices during the stimulation of synaptic inputs in combination with suprathreshold depolarizing current injections that, when presented alone, elicited a single action potential. Evoked IPSPs were hyperpolarizing from Vrest by −5.6 ± 1.0 mV (n = 6) (Fig. 2C, middle). Similar hyperpolarizations were observed in whole-cell recordings with Cl−pip = 5 mm, the estimated chloride concentration in SBCs around hearing onset (Milenkovic and Rübsamen, 2011) (Fig. 2D, middle; ΔVm = −6.4 ± 0.7 mV, n = 16; p = 0.54 for whole-cell vs gramicidin perforated patch, t test). Under both recording conditions, larger synaptically evoked hyperpolarizations correlated with AP failures (termed Depolfail-current injection that evokes no AP due to synaptic inhibition) (Fig. 2C,D, right). The amplitude of hyperpolarization (ΔVm) was linearly related to the depolarization rising slope (Fig. 2E; mean r = 0.94 ± 0.01, p < 0.001 for each cell, n = 8). Thus, the synaptically evoked inhibition flattens the depolarization slope rate and has, therefore, the comparable effect to acoustic stimulation within inhibitory sidebands in vivo. Together, our data obtained in vivo and in slice recordings are consistent with earlier studies by showing that the inhibitory inputs have the potency to suppress the AP generation in SBCs (Caspary et al., 1994; Gai and Carney, 2008; Dehmel et al., 2010; Xie and Manis, 2013).
Surprisingly slow kinetics of inhibitory transmission
To assess the properties of inhibitory synaptic transmission onto identified SBCs (Fig. 3A), eIPSCs were recorded at different input frequencies (Fig. 3B). The baseline-to-peak amplitude of single events ranged from −132 to −1164 pA, with a mean ± SD = −494.8 ± 270.2 pA, n = 43. Comparison across the cells yielded neither a correlation between the IPSC rise time and amplitude (rs = 0.04, p = 0.78, n = 43) nor between decay time and amplitude (rs = −0.13, p = 0.41, n = 43). Also, no systematic variation was seen across the ages (P22-P33) with respect to current decay times (rs = −0.07, p = 0.91), rise times (rs = −0.23, p = 0.14), and amplitudes (rs = 0.25, p = 0.11). The single-event eIPSCs exhibited remarkably slow synaptic decays (τw = 23.7 ± 5.3 ms; 10%–90% rise time = 0.46 ± 0.18 ms, n = 43), compared with decay rates measured from neurons in the superior olivary complex (Awatramani et al., 2004; Magnusson et al., 2005; Fischl et al., 2012; Stange et al., 2013; Kramer et al., 2014). Furthermore, the synaptically evoked hyperpolarization under gramicidin-perforated-patch condition had also a slow decay time constant (τw = 18.5 ± 4.0 ms, n = 6). Thus, repetitive stimulation caused a current summation in consecutive events. At 10 Hz stimulation, the first event had on average the highest conductance (calculated from Ibaseline), whereas for frequencies 50–333 Hz the maximal total conductance was reached between IPSC pulses 7–9 (Fig. 3C; maximum ginh 10 Hz = 22.1 ± 3.1 nS, n = 7; 50 Hz = 36.7 ± 3.3 nS, n = 7; 100 Hz = 51.3 ± 3.2 nS, n = 24; and 333 Hz = 68.0 ± 5.4 nS, n = 12). Because of frequency-dependent current summation, the average amplitude of the last three IPSCs in a 50-pulse train increased fourfold between 10 and 333 Hz (10 Hz = −0.26 ± 0.03 nA; 333 Hz = −1.01 ± 0.08 nA, p < 0.001).
Analysis of the phasic peak currents (Ipeak, Fig. 3B, inset) revealed a prominent frequency-dependent depression (Fig. 3D; normalized mean 40–50th Ipeak 10 Hz = 0.58 ± 0.05, n = 7; 50 Hz = 0.35 ± 0.03, n = 7; 100 Hz = 0.18 ± 0.02, n = 24; p < 0.001, ANOVA). The prolongation of consecutive synaptic currents is demonstrated by τw values of the last event in a 100 Hz train (Fig. 3E; τw 1st = 23.1 ± 0.7 ms, 5th = 44.3 ± 3.3 ms, 10th = 54.3 ± 4.6 ms, 20th = 62.9 ± 4.6 ms, 30th = 68.4 ± 6.0 ms, 40th = 73.2 ± 5.3 ms, 50th = 71.7 ± 3.7 ms; n = 6). Analogously, the current decay also showed a frequency-dependent prolongation of the last IPSC (Fig. 3F; mean τw 50th 0.1 Hz = 25.1 ± 2.2 ms; 10 Hz = 30.8 ± 2.6 ms; 20 Hz = 40.6 ± 1.5 ms; 50 Hz = 52.5 ± 3.6 ms; 100 Hz = 67.7 ± 2.7 ms; 200 Hz = 79.3 ± 2.7 ms; and 333 Hz = 77.6 ± 4.4 ms; n = 6). The maximum decay time was reached at 200 Hz with no further increase toward higher stimulus frequency. In summary, these data indicate activity-dependent buildup of inhibitory conductance due to slow current decay. Furthermore, a synaptic depression shapes the current profile during ongoing inhibitory activity.
GABA enhances and slows down the predominantly glycinergic inhibition
To determine the contribution of transmitters involved in synaptic inhibition, eIPSCs were pharmacologically characterized by superfusion of glycine-R and GABAA-R antagonists. First, the specificity of antagonists at the concentrations used was confirmed in control experiments. The currents evoked by focal puff application (5 ms) of either glycine or GABA were efficiently blocked by strychnine and SR95531, respectively (Fig. 4A; Iglycine 0.5 mm = −1.67 ± 0.12 nA, +strychnine 0.5 μm = −0.08 ± 0.02 nA, n = 6, p < 0.001, paired t test; IGABA 0.5 mm = −1.23 ± 0.36 nA, +SR95531 20 μm = −0.01 ± 0.007 nA, n = 5, p < 0.001, paired t test).
While SR95531 (20 μm) did not affect the Iglycine amplitude (Fig. 4C; p = 0.28, n = 6, ANOVA followed by pairwise comparisons), the current was progressively blocked by increasing strychnine concentrations, achieving maximal inhibition at 1 μm (98.8 ± 0.38%, n = 6) (Fig. 4B,C). The effect of strychnine was reversible, again reaching 93% of the initial Iglycine after a washout of 30 min (Fig. 4B). Similarly, strychnine concentrations of 0.5 and 1 μm had no effect on IGABA (Fig. 4D; p = 0.99 for 0.5 vs 1 μm, ANOVA). Hence, it was concluded that, at the concentrations used, strychnine and SR95531 specifically blocked the respective glycine and GABA evoked currents. Further corroboration was provided by showing equal inhibitory potency of increasing SR95531 concentrations on eIPSCs (Fig. 4E; eIPSC = −0.52 ± 0.05 nA, n = 17; +SR95531 10 μm = −0.40 ± 0.07 nA, 82 ± 2% of control, n = 6; +SR95531 20 μm = −0.43 ± 0.05 nA, 84 ± 3% of control, n = 11; p = 0.89 for 10 vs 20 μm SR95531, ANOVA). Similar effects of SR95531 were observed on spontaneous IPSCs, suggesting that the GABAA-R contribution is an intrinsic property of inhibitory synapses and not a result of synaptic stimulation (Fig. 4E; sIPSC = −40.6 ± 3.6 pA, +SR95531 20 μm = −33.8 ± 2.7 pA, 84 ± 4% of control, n = 5, p = 0.03, paired t test). Based on these data, it seems unlikely that, in our experiments, 20 μm SR95531 acted as a low-affinity antagonist of glycine receptors (Beato et al., 2007).
After validating the specificity of antagonists, the contribution of both glycine and GABAA receptors to IPSCs was quantified (Fig. 5). The currents were differentially affected by 0.5–1 μm strychnine and 20 μm SR95531 (Fig. 5A). The glycinergic component, measured under SR95531, dominated the peak amplitude (Fig. 5A,B; single or first eIPSC: control = −0.60 ± 0.09 nA, +SR95531 = −0.44 ± 0.05 nA, p < 0.01; average of 7–9th baseline IPSCs at 100 Hz: control = −1.06 ± 0.10 nA, +SR95531 = −0.89 ± 0.07 nA, p < 0.05; average of 48–50th baseline IPSCs at 100 Hz: control = −0.59 ± 0.08 nA, +SR95531 = −0.49 ± 0.06 nA, p = 0.12; n = 15; paired t test). Yet, the initially small GABAergic component conspicuously increased from 5 ± 1% to 12 ± 3% during ongoing 100 Hz activity (single: control = −0.65 ± 0.12 nA, +strychnine 0.5–1 μm = −0.04 ± 0.02 nA, p = 0.001; 7–9th baseline IPSCs: control = −1.11 ± 0.14 nA, +strychnine 0.5–1 μm = −0.09 ± 0.01 nA, p < 0.001; 48–50th baseline IPSCs: control = −0.55 ± 0.13 nA, +strychnine 0.5–1 μm = −0.06 ± 0.01 nA, p < 0.01; n = 7, paired t test).
Synaptically released GABA not only contributed to the IPSC amplitude, but it also prolonged the current kinetics in an activity-dependent manner (Fig. 5C). The decays of cotransmission, glycinergic, and GABAergic eIPSCs were similar for single events (τw: control = 25.7 ± 1.2 ms, n = 25; +SR95531 20 μm = 23.2 ± 1.2 ms, n = 13; +strychnine 1 μm = 26.2 ± 3.2 ms, n = 12, p > 0.5, two-way ANOVA followed by Holm-Sidak test). However, the differences became evident at higher input frequencies (interaction of frequency and drug superfusion p = 0.02). At 100 Hz, control IPSCs had time constants between the fast glycinergic and slow GABAergic currents (Fig. 5Ci; τw 50th: control = 64.2 ± 2.1 ms, n = 25 vs +SR95531 20 μm = 56.3 ± 2.8 ms, n = 14, p < 0.05; control vs +strychnine 1 μm = 81.1 ± 6.4 ms, n = 11, p < 0.001, two-way ANOVA followed by Holm–Sidak test). Notably, the kinetics of control events at 333 Hz was similar to isolated GABAergic IPSCs (τw 50th: control = 76.9 ± 2.7 ms, n = 21 vs +SR95531 20 μm = 65.6 ± 3.5 ms, n = 10, p < 0.05; control vs +strychnine 1 μm= 84.1 ± 5.3 ms, n = 11, p = 0.11, two-way ANOVA with post hoc Holm–Sidak test). Together, these results demonstrate that GABA progressively slows the inhibition with increasing rates of activity. GABA-mediated prolongation of inhibitory synaptic events was also shown for other central synapses (Jonas et al., 1998; Nabekura et al., 2004; Apostolides and Trussell, 2013). The eIPSC rise times were similar for dual glycine-GABA and for pharmacologically isolated events (10–90% rise time control = 0.42 ± 0.05 ms, +strychnine = 0.46 ± 0.05 ms, n = 7, p = 0.63; control = 0.49 ± 0.03 ms, +SR95531 = 0.48 ± 0.03 ms, n = 14, p = 0.60, paired t test for both). Thus, given the rapid and similar onset of all inhibitory currents, it is unlikely that the slow decay kinetics are due to dendritic filtering or space clamp issues.
GABA increases the potency of inhibition
The in vivo results demonstrate that acoustically evoked inhibition can constrain AP firing in SBCs. By means of slice recordings, we showed that the inhibitory potency is enhanced by the synergistic action of glycine and GABA. Next we investigated how each of the two inhibitory mechanisms interferes with AP generation. Under current-clamp, the APs were elicited by a 5 ms current injection set slightly over the threshold to reliably evoke APs during a 100 Hz train (Fig. 6A, AP probability of 1 in Fig. 6E). Concurrent stimulation of inhibitory synaptic input (10 pulses at 100 Hz, each pulse presented 5 ms before corresponding depolarizing current injection) reliably blocked the generation of APs. This block outlasted the duration of inhibitory stimulation by 45 ± 3 ms (n = 12), consistent with the prolonged efficacy of inhibition (Fig. 6B,F). The duration of this inhibitory after-effect is comparable with the τw of 54.3 ± 4.6 ms, measured for the 10th eIPSC at 100 Hz (compare Fig. 3E). This implies that the slow decay of inhibition can account for the failures in AP firing.
Application of 0.5 μm strychnine completely abolished the inhibition of APs (Fig. 6C). Blocking GABAA-R had a more subtle effect: The AP inhibition was delayed and lost its potency faster after the end of inhibitory stimulation (Fig. 6D). Figure 6E, F demonstrates the prolonged poststimulation inhibitory effect of the glycine-GABA coinhibition with respect to the sole glycinergic inhibition. Summary data in Figure 6E show that the time when the first AP occurs is significantly longer after full, compared with pharmacologically isolated glycinergic, inhibition. Although the GABAergic component had such a prominent effect on capacity of inhibition, it just moderately contributed to the membrane hyperpolarization evoked by stimulation of synaptic inputs. Figure 6G summarizes the relative contribution of each transmitter component to total hyperpolarization (ΔVm single: control=-6.6 ± 0.9 mV, +SR95531=-5.0 ± 0.6 mV, n = 9, p = 0.011; control=-6.1 ± 1.1 mV, +strychnine = −0.19 ± 0.1 mV, n = 7, p = 0.002; 10th at 100 Hz: control = −15.3 ± 1.3 mV, +SR95531 = −12.7 ± 1.1 mV, n = 13, p < 0.001; control = −9.9 ± 1.0 mV, +strychnine = 0.65 ± 0.3 mV, p < 0.001, n = 7; paired t test for all). These data are in agreement with the results from voltage-clamp experiments (Fig. 4E). Although synaptically released GABA seems not capable of preventing the AP generation per se, it appears necessary to adjust the duration of AP inhibition through a synergistic action with glycine.
Interaction of excitation and inhibition in silico
These results give rise to the question of how efficient the slow inhibition is in preventing calyceal inputs from triggering APs considering a wider range of physiologically relevant conductances. To get a closer insight into the inhibitory impact on the powerful and fast excitatory synaptic conductance (i.e., endbulb-like excitation), we performed in silico experiments. The model inhibitory synaptic conductance was matched to the data from slice recordings. Simulating synaptic currents without considering the plasticity and rate-dependent reduction of conductance (e.g., 23.8 nS initial conductance and fixed 24 ms decay time-constant) only poorly matched the data measured in vitro (Fig. 7A, compare Fig. 3C). A better match of the current profile was achieved after including the rate-dependent increase of τ decay (Fig. 3F), which counteracts depression during the stimulus train by facilitating summation. Conductance waveforms underlying the synaptic model mechanism are shown in Figure 7B. The respective conductance values for glycine and GABAergic components were implemented in the model, whereas other aspects, such as onset and decay time-constants, were simulated as identical for both parts.
Referring to the data from slice experiments, the interactions between inhibition and excitation were assessed for 100 Hz trains. Initially, a simplified excitatory input with fixed conductance was used for simulations (Figs. 7 and 8) despite the fact that short-term plasticity of the endbulb of Held synapse has been demonstrated in many studies (Bellingham and Walmsley, 1999; Wang and Manis, 2008; Yang and Xu-Friedman, 2008). As these particular experiments aimed to assess the interaction of the two inhibitory transmitters and their effect on AP firing, the model was simplified by investigating the dynamics of inhibition acting on a constant excitation. These data are directly comparable to the results in Figure 6, which also had a constant strength of excitation due to current stimulation (but see also the extended model in Fig. 9). Even with inhibition preceding the excitation, individual inhibitory events were too weak to overcome endbulb-like excitation of 50 nS. However, the strong summation of inhibitory conductance during pulse trains eventually interfered with the excitation (Fig. 7Ci). This was, however, only true for the synergistic glycine + GABA conductance (Fig. 7Cii, black trace). The GABAergic component, albeit small in amplitude, increases the total effectiveness of inhibition. This is illustrated in Figure 7D by plotting the conductance of excitatory events and the effective inhibitory conductance during each interaction in a 100 Hz train combined as coordinates (glycine + GABA: white circles; glycine only: black crosses). The trajectories of EI-coordinates are superimposed on a 2D plot illustrating the efficacy of inhibition over a wide range of inhibitory and excitatory conductances. Here, hotter colors indicate higher resulting membrane potentials (note that peak Vm of respective EPSP and AP events is plotted with corresponding color-coded waveforms on the right). These data demonstrate that, for the assumed endbulb-triggered conductance of 50 nS, the additional GABAergic component expands the effective inhibitory strength into the range of secure spike suppression.
We next assessed the impact of inhibitory conductance during simulated 100 Hz neuronal activity in analogy to the experiments shown in Figure 6. For this, a range of initial inhibitory conductances was regarded interacting with a fixed excitatory conductance. As measures of inhibitory efficacy, the fraction of spikes in the train (Fig. 7E) and the time of spike occurrence after the last inhibitory event were quantified (Fig. 7F). Figure 7E illustrates the inhibitory conductances necessary to suppress spikes in a 100 Hz train. The glycinergic component alone abolished spiking for excitation <42 nS but failed to suppress spiking when excitation was >57 nS. The glycine + GABA inhibition efficiently blocked APs when excitation was <45 nS and had no effect when gexc was >65 nS.
These results strongly corroborate the hypothesis that synergistic glycine-GABA action extends the inhibitory potency to a higher range of excitatory conductances. Figure 7F shows the duration of inhibitory effect in the model (compare Fig. 6E). In contrast to slice measurements, there was no difference between the sole glycinergic and the glycine-GABA inhibition. In both cases, the effective inhibitory conductance ended too fast to suppress the next event in the train. This indicates that our model is probably not taking into account all aspects of inhibitory synapse dynamics. The proportionally stronger prolongation of GABAergic decay rate, compared with glycinergic decay rate, was not included in the model. This most likely underlies the prolonged effect of inhibition on spike probability in slice recordings.
Together, the in silico data suggest that the slow inhibition found in the in vitro experiments has the potency to effectively suppress endbulb-triggered AP generation, whereby the inhibitory strength and the dynamics are finely tuned through interaction of glycinergic and GABAergic mechanisms.
Slow inhibitory dynamics in vivo
Excitatory acoustic stimulation at SBC's CF evokes a primary-like response, characterized by a rapid increase in AP spiking followed by a steady-state activity during ongoing stimulation (Blackburn and Sachs, 1989; Typlt et al., 2012). The response to excitatory auditory nerve (AN) fiber input, triggered by a tone burst at CF (20 dB above threshold), is characterized by the initial increase in firing rate (onset) and the spike rate reduction at the end of the stimulus (offset) (Fig. 8Aii, gray line; onset: 3.6 ms, offset: 1.8 ms). When the tone burst was presented in the unit's inhibitory sideband (Fig. 1B), the AP rate was reduced, although the excitatory input to the SBC remained constant, as seen from the steady EPSP rate shown as mean for 16 cells (Fig. 8Ai). The time course of AP inhibition in vivo, assessed as the dynamics of the increase in AP failure rate during inhibitory sideband stimulation, showed a slow onset (half-rise time from stimulus onset = 7.3 ms, n = 16) (Fig. 8Aii, red). Moreover, the failure fraction tapered down during ongoing inhibitory stimulation, being significantly lower at its end (15–30 ms vs 80–95 ms; 0.64 ± 0.19 vs 0.58 ± 0.18, n = 16, p < 0.01, paired t test). Also, the release from inhibition was remarkably slower than the offset of excitation (Fig. 8Aii; decay time constant computed from fitting the average trace for 16 units: τinhibition = 9.8 ms [95% confidence interval: 5.9, 14.1 ms], τexcitation = 1.4 [95% confidence interval: 0.7, 2.0 ms], black and gray, respectively). One conspicuous effect of inhibition, consistent across in vivo and slice experiments, was the change of the EPSP slope and prolongation of the EPSP-to-AP transition time (Fig. 2). Notably, the prolongation of the EPSP-to-AP delay was most prominent after the onset of inhibitory stimulation, reaching maximal delay values within initial 30 ms. Although the delay remained prolonged during ongoing inhibitory stimulation (Fig. 8Aiii; average delay for 150–230 ms spontaneous activity = 0.23 ± 0.02 ms; average delay during 10–90 ms of inhibitory stimulation = 0.25 ± 0.02 ms, n = 16, p < 0.001, paired t test), this effect became progressively weaker with time (EPSP-to-AP delay during 15–30 ms of inhibitory stimulation = 0.25 ± 0.02 ms; during 80–95 ms = 0.24 ± 0.02 ms; n = 16, p < 0.01, paired t test). After the release from inhibition, the delay decreased to initial values with a time constant of τ = 14.7 ms (95% confidence interval: 12.1, 17.4 ms), calculated from fitting the average trace from 16 units.
Finally, in silico experiments were used to explore whether the EPSP-to-AP delay prolongation could be reproduced with inhibitory parameters assessed from slice experiments. Indeed, during excitation–inhibition interaction, the AP peak delay strongly depended on the conductance strength revealing that AP initiation close to threshold (i.e., excitation at just marginally suprathreshold values) caused a later occurring peak (Fig. 8B). The dynamics of EPSP-to-AP prolongation showed a time course that was approximately comparable with in vivo recordings (Fig. 8C, compare Fig. 8Aiii). Increased delay values were reached within the first few events, but after reaching a maximum, the delay was progressively reduced depending on inhibitory conductance strength. However, small variations in the balance between excitatory and inhibitory conductance can have tremendous impact on the EPSP-to-AP delay in the model (Fig. 8B). Thus, in this context only a qualitative comparison can be made between the in vivo and in silico results. These data imply that the dynamic efficacy of inhibition, as observed in vivo, crucially depends on the inhibitory conductance profile characterized by the current depression and tau prolongation (Fig. 7A). Modeling the inhibitory conductance without τ plasticity and/or depression yields a time course of delay prolongation different from the in vivo data (data not shown).
To explore the properties of the inhibition during in vivo-like spontaneous and sound-evoked activities, the inhibitory conductance profiles measured in slice experiments were used. To this end, AN spike trains (2 s) were generated and implemented to drive excitatory as well as inhibitory synaptic mechanisms. The inhibition was trailing the excitation by 1 ms, assuming the engagement of an additional synapse activated by AN fibers (Wickesberg and Oertel, 1990; Xie and Manis, 2013; Campagnola and Manis, 2014). Spontaneous AN events (30 ± 5 spikes/s) reliably caused action potentials in the SBC although the inhibitory conductance showed activity-dependent buildup (Fig. 9Ai–Ci, Aii–Cii). Modeling of stochastic and dynamic endbulb conductances revealed failures occurring during spontaneous activity. Thus, in combination with variable excitatory conductances, tonic inhibition caused some failures during spontaneous activity, as shown by the analysis of 100 repetitions (Fig. 9Aiii–Ciii). However, failure incidence strongly increased during the response to simulated sound bursts (201 ± 37 spikes/s) (red, Fig. 9Aii–Cii). The quantification of spike probability during simulated acoustic responses for all three conditions yielded different results: in the fixed endbulb condition (black line), the probability dropped from p = 0.95 ± 0.11 to p = 0.23 ± 0.22, in the stochastic endbulb condition (orange line) from p = 0.88 ± 0.08 to p = 0.4 ± 0.17, and in the dynamic endbulb condition (blue line) from p = 0.74 ± 0.13 to p = 0.1 ± 0.21 (Fig. 9D). Notably, in the dynamic endbulb condition, the spike probability dropped to zero in the sustained part of the simulated acoustic responses due to strong depression at high input rates. In addition to increased refractoriness related to short interspike intervals (data not shown), the reduction in spike probability was mainly caused by a strong increase in tonically active inhibitory conductance from 21.5 ± 0.9 nS to 43.8 ± 7.8 nS (Fig. 9G). The spike probability measured in vivo changed from p = 0.73 to p = 0.36 during acoustic stimulation in the inhibitory sideband and the inhibition offset had the time constant of 9.8 ms (Fig. 8Aii). Compared with these data, the offset of inhibition computed by our models was in the same range being 5 ms for the fixed and 20 ms for the stochastic excitatory conductance (measured as the time constant of the monoexponential fit to the spike probability traces after the acoustic stimulation; Fig. 9D). Recovery from inhibition was much slower with dynamically depressing excitation.
The tonic inhibition also caused an increase in AP delay of successful events, quantified as the interval between the peaks of the EPSP and AP waveform components. This effect was most pronounced for the fixed endbulb condition (Fig. 9E, black line; EPSP-to-AP delay spontaneous = 0.20 ± 0.01 ms; stimulation = 0.24 ± 0.03 ms; maximum = 0.31 ms), and less strong for the stochastic endbulb condition (orange line; EPSP-to-AP delay spontaneous = 0.20 ± 0.01 ms; stimulation = 0.22 ± 0.02 ms; maximum = 0.25 ms). In the dynamic endbulb condition, EPSP-to-AP delay initially increased up to 0.26 ms but could not be determined throughout the ongoing simulated acoustic response, as only few successful events occurred. Nevertheless, the EPSP-to-AP delay was still increased to 0.26 ms at 34 ms after stimulus and returned to normal only after >50 ms. On the other hand, the effect on the delay prolongation decayed with a τ = 15 ms for fixed and τ = 48 ms for the stochastic synapse model, which is comparable to τ = 14.7 ms measured in vivo (Fig. 8Aiii).
The time course of the mean endbulb EPSG underlying the dynamic condition showed a tonic depression to 81.5% (61.1 ± 2.8 nS; Fig. 9F) during the spontaneous activity. Depression further increased during simulated sound evoked activity resulting in 43.5% EPSG of the initial value (32.6 ± 11.1 nS). At the end of the simulated sound-evoked response, EPSG was as low as 26.3% (19.7 nS). Consistent with the slice data, the activity-evoked inhibition showed a dynamic change of the phasic amplitude of individual events (Fig. 9H; spontaneous = 18.1 ± 1.3 nS; stimulation = 6.9 ± 4.3 nS) and in τ decay (Fig. 9I; spontaneous = 40.3 ± 6.5 ms; stimulation = 64.1 ± 13.4 ms). After acoustic-like stimulation at the average rate of 201 ± 37 spikes/s, the prolonged τ decay returned back to initial value with a time constant of 99 ms. Very similar decay rate was measured in slice experiments after the 200 Hz stimulation (weighted τ decay of ∼80 ms; Fig. 3F).
In summary, because of the rate-dependent dynamics of the inhibitory input, its efficacy against a fixed excitatory conductance gradually increased during the initial part of the stimulus response, reached a maximum, and thereafter declined despite the ongoing presentation of the inhibitory stimulus. Against the stochastically variant excitatory conductance, the overall potency of the inhibition was on average lower but still efficient in the ongoing part of the simulated stimulus response. Counteracting the dynamically depressing excitatory conductance, inhibition was overwhelmingly efficient, completely silencing all spiking activity in the model SBC for >50 ms. Together, the data obtained by using the fixed and stochastic endbulb model are compatible with the time course of inhibitory action measured in vivo (Fig. 8). Hence, the modeling results demonstrate that synergistic glycine/GABAergic transmission yields inhibitory input that effectively constrains AP generation at input rates characteristic for acoustically evoked activity.
Discussion
Synaptic inhibition in auditory brainstem circuits of mammals is mostly mediated by glycine, but few auditory synapses use both GABA and glycine beyond the early postnatal development (Balakrishnan and Trussell, 2008; Balakrishnan et al., 2009; Apostolides and Trussell, 2013). To date, it remained elusive why both glycine-R and GABAA-R are engaged with inhibition of SBCs and what is the inherent advantage of such a two-transmitter inhibitory mechanism. We addressed this question using the endbulb of Held-SBC synapse, which allows for an in vivo assessment of input–output function based on the analysis of systemic, acoustically evoked activity. This model synapse enables characterization of a combined glycine-GABA inhibition in slice experiments and further in vivo evaluation of its impact. We combined both approaches with in silico modeling to demonstrate that the acoustically evoked inhibition enables only strong or coincident excitatory inputs to generate APs. Hence, the slow glycine-GABA inhibition acts as a dynamic high-pass filter. The time course of acoustically evoked AP inhibition in vivo matches the inhibitory conductance profile characterized by a transient summation followed by depression. Comparable inhibitory action was observed in the models assuming fixed or stochastically changing excitatory conductances, whereas modeled inhibitory conductance was apparently too potent when interacting with dynamically depressing excitation. The dual glycine-GABA synaptic transmission endows each transmitter with a distinct function: whereas the glycinergic component dominates the inhibitory conductance of evoked and spontaneous IPSCs, the GABAergic fraction extends the offset of inhibition in an activity-dependent manner. By increasing the inhibitory strength and prolonging the decay time, signaling via GABAA-R effectively increases the capacity of inhibition to counteract the endbulb of Held excitation.
Sole glycinergic versus synergistic glycine-GABA inhibition
The inhibitory inputs of noncochlear origin (Smith and Rhode, 1989; Benson and Potashner, 1990; Saint Marie et al., 1991; Wenthold, 1991; Oertel and Wickesberg, 1993; Ostapoff et al., 1997; Arnott et al., 2004; Campagnola and Manis, 2014) terminate on SBCs in three classes of synapses: glycinergic, GABAergic, and synapses containing both (Kolston et al., 1992; Juiz et al., 1996; Mahendrasingam et al., 2004). The postnatal structural development of respective inputs spans the period of early auditory experience, reaching adult-like configuration of terminals on SBCs ∼P21 (Luján et al., 2008). These results are validated by the demonstration of glycine-R (Wenthold et al., 1988) and GABAA-R (Campos et al., 2001) on SBCs, as well as the respective uptake transporters in presynaptic terminals (Mahendrasingam et al., 2000). Our data demonstrate that spontaneously occurring (sIPSCs) and events induced by synaptic stimulation (eIPSCs) consisted of both glycinergic and GABAergic components. Although the used experimental protocols enabled us to quantify the physiological impact of the glycine-GABA transmission, we cannot draw a conclusion about the possible corelease of both transmitters from the same vesicle as shown for the spinal and brainstem motoneurons (Jonas et al., 1998; Russier et al., 2002). As indicated above, such an assumption seems plausible, but measurements of miniature IPSCs are necessary to verify this hypothesis.
Across inputs on SBCs, the GABA content in the presynaptic terminals is apparently rather variable, contrasting with the larger and more constant synaptic amount of glycine (Mahendrasingam et al., 2004). Accordingly, the inhibitory strength of a sole GABAA signaling was not sufficient to interfere with spike activity per se, as indicated by our slice and modeling data. Still, the GABA contribution to the inhibition cannot be considered a minor one because the impact of the GABAergic fraction can account for up to ∼12% of the total inhibitory conductance during ongoing activity. The activity-dependent depression of glycinergic currents is attributable to the relative increase in GABA component (Ipeak; compare Fig. 5), which progressively prolongs IPSC kinetics with increasing activity.
Our data from gerbils and results from mice (Xie and Manis, 2013) consistently show slow kinetics of inhibitory currents in SBCs. This seems, prima facie, surprising given the role of SBCs in the submillisecond processing of the temporal information of low-frequency sounds, which is then binaurally integrated by the coincidence detector neurons of the MSO (Smith et al., 1993; Grothe, 2000). The question arises: how can acoustically evoked inhibition, which is orders of magnitude slower than acoustically triggered excitation, contribute to the retention (or even improvement) of temporal resolution? Particularly at higher input frequencies, inhibition is remarkably prolonged by GABA. Yet, the glycinergic IPSCs are also summating due to decay time constants lasting tens of milliseconds. Possible mechanisms for inhibitory current summation are the change to asynchronous transmitter release at high frequencies and/or transmitter pooling due to poor clearance from the synaptic cleft (Lu and Trussell, 2000; Balakrishnan et al., 2009). Hence, depending on the activity, SBC can be controlled by a tonic inhibition arising through current summation (present study and Xie and Manis, 2013), as also shown for DCN granular cells (Balakrishnan and Trussell, 2008). As we presently found no correlation between the amplitudes and the rise times of IPSCs, the longer decay times and larger amplitudes of evoked versus spontaneous events might point to transmitter rebinding and multivesicular release. The fast glycinergic IPSCs in AVCN stellate cells, in MSO, and MNTB neurons, enable phasic inhibition (Awatramani et al., 2004; Magnusson et al., 2005; Chirila et al., 2007; Couchman et al., 2010; Xie and Manis, 2013). Contrary to this, the temporal fidelity in SBC spiking is controlled by a synergistic and rather slow glycine-GABA transmission through filtering out the weak, or poorly timed, excitatory inputs.
Contribution of inhibition to auditory processing
Biophysically, SBCs are characterized by a fast membrane time constant, low input resistance, and nonlinearity of voltage-gated conductances. These properties control the phasic AP responses to sharply rising, large EPSPs (Rothman and Manis, 2003; McGinley and Oertel, 2006; Cao et al., 2007). During spontaneous activity in vivo, postsynaptic AP failures occur with an incidence of up to ∼20% probably caused by stochastic fluctuations in excitatory synaptic strength (Englitz et al., 2009; Kuenzel et al., 2011). Our data support this assumption by demonstrating that AP failures are associated with smaller EPSPs, despite constant presynaptic spike amplitudes (Fig. 1). Leaving the excitatory input rate unchanged, the acoustic stimulation within the inhibitory sideband can increase the failure fraction up to ∼65%. In such cases, the steepness of the EPSP slope was reduced, supposedly below the integration time window required for AP generation (McGinley and Oertel, 2006). Hence, the dynamic adjustment of hyperpolarizing inhibitory strength mediated by glycine and GABA, as shown by our data, is likely to tune the fidelity of the endbulb of Held synapse to fast rising and large EPSPs.
The spiking activity of SBCs can be synchronized to a particular sine wave phase of a low-frequency tone burst (phase-coupling). Surprisingly, the precision of phase-coupling is improved compared with the AN input (Joris et al., 1994; Paolini et al., 2001; Joris and Smith, 2008). Two mechanisms possibly underlie this improvement: (1) jitter reduction through coincidence detection of convergent inputs (Kuhlmann et al., 2002; Xu-Friedman and Regehr, 2005); and (2) synaptic inhibition mediated by glycine and GABA. Iontophoretic administration of glycine/GABAA-R antagonists during in vivo recordings and acoustic stimulation at units CFs deteriorated phase-coupling and increased the firing rates at higher stimulus intensities (Dehmel et al., 2010). The latter is probably due to largely overlapping excitatory and inhibitory response areas (on-CF inhibition) (Winter and Palmer, 1990; Kopp-Scheinpflug et al., 2002; Kuenzel et al., 2011). According to our results, the activation of inhibition at higher sound pressure levels might induce larger inhibitory current summation (tonic inhibition due to higher firing rate or recruitment of inputs), leading to both nonmonotonic rate-level function and tightening of phase-coupling (Dehmel et al., 2010). As the sound-driven responses may cause variations in the EPSP slope steepness, we cannot estimate the potency of inhibition at CF, but the overall output would again depend on the actual excitatory conductance determined by the amount of depression, which correlates the input frequency and synchrony of excitatory inputs. In this respect, it is necessary to point out that the observed reduction in APs is not caused by cochlear suppression, as evidenced from the constant excitatory input activity during acoustic stimulation in the inhibitory sideband.
Most in vivo studies that reported a contribution of both glycine and GABA to acoustic signal processing in the CN could not specifically assign particular role to either transmitter (Caspary et al., 1994; Backoff et al., 1997; Dehmel et al., 2010). It remains an intriguing question whether the respective transmitters may be differentially associated with specific auditory mechanisms, as suggested previously (Ebert and Ostwald, 1995; Gai and Carney, 2008). Our results argue that the potent glycinergic inhibition on SBCs primarily determines the slow filtering properties. However, this study also reveals the particular role of GABAA receptors on SBCs. By mediating dynamic enhancement of inhibitory strength and shaping its duration, postsynaptic GABAergic transmission can contribute to an effective gain control.
Footnotes
This work was supported by the Deutsche Forschungsgemeinschaft Grant MI 954/2-1 to I.M., Grant MI 954/1-1 to I.M., Grant RU 390/19-1 to R.R. and J.N., Grant GRK 1097 to J.N. and C.K., and Grant KU 2529/2-1 to T.K., A.K. was supported by the DAAD project “Akademischer Neuaufbau Südosteuropa.” We thank Stefan Oline for providing the slope analysis routine for slice recordings.
The authors declare no competing financial interests.
- Correspondence should be addressed to Dr. Ivan Milenkovic, Carl Ludwig Institute for Physiology, Faculty of Medicine, University of Leipzig, Liebigstrasse 27, D-04103 Leipzig, Germany. Ivan.Milenkovic{at}medizin.uni-leipzig.de