1. Neuroscience
Download icon

Gradients in the biophysical properties of neonatal auditory neurons align with synaptic contact position and the intensity coding map of inner hair cells

  1. Alexander L Markowitz
  2. Radha Kalluri  Is a corresponding author
  1. Neuroscience Graduate Program, University of Southern California, United States
  2. Department of Otolaryngology, Keck School of Medicine, University of Southern California, United States
  3. Zilkha Neurogenetic Institute, Keck School of Medicine, University of Southern California, United States
Research Article
Cite this article as: eLife 2020;9:e55378 doi: 10.7554/eLife.55378
7 figures, 5 tables and 1 additional file


Figure 1 with 1 supplement
Semi-intact preparation combines SGN biophysics with identification of pattern of connectivity with hair cells.

(A) Schematic of the innervation at auditory hair cell along with key anatomical structures provided with the semi-intact preparation presented in the XZ plane. (B) Confocal image of two type I spiral ganglion neurons (SGN) injected with biocytin in the acute semi-intact preparation of rat middle-turn cochlea. Hair-cells and biocytin-labeled fibers are visualized by immunolabeling for Myosin VI (red; hair cells) and streptavidin conjugated to Alexa Flour 488 (green, fiber). (C) Myelinating Schwann cells are mechanically removed by suction pipettes to make SGN somata accessible to recording electrodes. D.1-F.1. Confocal images in the XY plane show type I SGN afferent fibers making exclusive connections with inner hair cells. D.2-F.2 The synaptic position of the afferent fiber onto the inner hair cell is assessed by examining the cross-section of the hair cell in XZ plane. (G) Type II SGN projecting radially and turn within the outer hair cell region. (H) A schematic of an individual inner hair cell with guides showing how the normalized basal position (NBP) is measured for each type I SGN. NBP is measured by the synaptic position along the radial axis of the hair cell,c, divided by the maximum length of the inner hair cell, L. Positive NBP values indicate that the synaptic position is on the modiolar-side of the inner hair cell, while negative NBP values indicate that the synaptic position is on the pillar-side of the inner hair cell. (I) The distribution of type I SGN contacts as a function of NBP. NBP values averaged at 0.11 +/- 0.23 with a range from −0.27 to +0.52 (n = 38).

Figure 1—figure supplement 1
Polarization vector shows that branching fibers have a preferred polarization to either modiolar or pillar face.

The net polarization of the fiber was computed by taking the sum of the individual polarizations divided by the total number of branches (Figure 1).

SGN can be qualitatively grouped by firing patterns in response to step current injections.

(A) Rapidly-accommodating neurons fire a single action potential. (B) Intermediately-accommodating neurons fire more than one action potential and eventually reach a steady-state potential. (C) Slowly-accommodating neurons fire multiple actions and do not reach a steady-state potential. (D) Graded-response neurons do not produce an action potential, but produce graded-depolarization in response to incrementing current steps. These neurons are capable of firing action potentials after their membrane potential are held at −80 mV (inset). (E) Action potential features including current threshold (Ithreshold), voltage threshold (Vthreshold), and after hyperpolarization potential time constant (AHP tau) vary among spiking neurons. The p-values indicate when two means were assessed to be significantly different (see Materials and methods); asterisks indicate p<0.001. (F) Action potential features were significantly correlated with normalized basal position (NBP). As NBP values increased, current thresholds increased (r(20) = 0.68, p=0.0007), voltage thresholds hyperpolarized (r(20) = −0.45, p=0.045), and AHP time constants became faster (r(20) = 0.62, p<0.0001). Best fit lines with 95% confidence intervals are plotted to show the gradients along the normalized basal position axis. (G) The distribution of firing pattern subgroups as a function of NBP shows a significant relationship between the firing-pattern subgroup and NBP (p=0.0097, t-test), with rapidly-accommodating neurons found primarily on the modiolar-face and intermediately-accommodating neurons found on the pillar-face of the inner hair cell. Graded-firing neurons are found throughout the NBP scale on both the modiolar and pillar faces of the hair cell.

Figure 2—source data 1

Current-clamp features versus normalized basal position data set for Figure 2F.

Response latencies are alternate predictors for normalized basal position.

There is a significant inverse relationship between first-spike latency and current thresholds. Black line shows a fit Ithresh = constant/latency. Color transition from red to blue indicates position on NBP scale. Type II SGN are colored black. Unlabeled and therefore unclassified SGN are in gray. (A) First spike latencies compared between three spiking SGN with different NBP values. Latencies are faster for more positive NBP (B) In non-spiking neurons, response latency is computed by averaging the time from the onset of the stimulus until peak depolarization potential over a series of current steps. Similar to the dependence of first-spike latency, average response latency is faster for modiolar-contacting SGN (NBP >0, red) than for pillar-contacting SGN (NBP <0, blue). (C) Response latencies for all spiking (black dot) and non-spiking (triangle) neurons plotted against NBP. In both spiking- and non-spiking SGN, response latencies become faster as NBP values increase (R2 = 0.55, p<0.001). Latencies are highly variable across type II fibers (gray symbols).

Figure 4 with 1 supplement
Net whole-cell conductances increase as SGN synaptic position moves from pillar to modiolar-face of the inner hair cell.

(A–C) Voltage-clamp responses of a modiolar-contacting (A), pillar-contacting type I (B), and a type II SGN (C). For each cell, we held the potential at −60 mV, followed by a series of voltage 400 ms steps from −120 mV to 70 mV, and then finally back to the holding potential. For each cell, we measured current-voltage relationships at approximate steady-state (400 ms indicated by the arrow). We measured the time course for outward current inactivation by fitting a single exponential curve to the current evoked by a +70 mV voltage step (inset). (D) Steady-state conductance in microSiemens as a function of command voltage of all SGN recorded with three measurements: the magnitude of conductance at −30 mV (g-30), the maximum conductance (gmax), and the half-activation potential (V1/2) via Boltzmann function. (E) Steady-state conductances of identified type I SGN as a function of normalized basal position. (F–I) Voltage clamp properties as a function of normalized basal position (NBP) (Spiking, black dot; non-spiking, triangle). Linear regression models were fitted to each parameter with statistical significance test of p<0.05. All parameters were significant except for V1/2.

Figure 4—source data 1

Voltage-clamp features versus normalized basal position for Figure 4F through 4I.

Conductance values are shown in units of micro Siemens.

Figure 4—figure supplement 1
Considering the influence of space-clamp and series resistance errors.

(A) Circuit diagram of an acute SGN in the semi-intact preparation clamped in whole-cell configuration. This diagram shows the series resistance (Rpip), membrane resistance (Rm), axial resistance (Ra), and a secondary compartment resistance from the neurite (Rcomp) with corresponding capacitors (C). (B) Voltage-clamp response of the transient capacitance current of an SGN in the semi-intact preparation fitted with a bi-exponential function (green trace). The fast (pink trace) and slow (blue trace) components of the bi-exponential function are plotted to show the contribution of the current coming for the somata (fast/pink) and neurite (slow/blue). (C) The relative size of the magnitude of the fast (A2) and slow (A1) components of the trace, indicating that the fast component dominates the fit and currents are measured primarily from the soma. (D) The fast component of the bi-exponential (A2) relative proportion to the sum of the bi-exponential function (A1+A2). A2 relative proportion dominates in both modiolar(red) and pillar (blue) contacting SGN. (E) Cell capacitance, measured as the integral of the capacitance current, is not correlated with basal position (r = 0.05, p=0.85). (F) Series resistance as a function of normalized basal with the mean trend line is shown in red. (G) After performing series resistance error corrections per cell, we see that the position correlation of maximum outward current (Imax) and NBP does not collapse, and remains correlated (R2 = 0.23, p=0.0003). If we simulate increasing levels of systematic series resistance errors by three (blue trace) and six (black trace) times the amount we measured, trends along NBP do not collapse.

Combining voltage-clamp and current-clamp features in multiple-variable regression models to predict normalized basal position (NBP) based on biophysics.

(A) Prediction model for spiking neurons (P3–P10) relies on Ithresh, gmax and tinact. (B) Prediction model that pools spiking and non-spiking neurons (P3–P10) is successful by using latency in place of current threshold (C,D), respectively. Biophysical gradients are still successful at predicting spatial position even after filtering data into narrow age bins. P3-P5 in C and P6-P8 in D. Predictor variables used in each model are presented at the bottom right of the predicted vs. actual plot. (E) Current threshold increase with increases in g-30. (F) Response latency as a function of g-30. Fit is a function of the form latency = constant/g-30.

Biophysical properties of SGN change with age but remain diverse throughout 3 weeks of post-natal development.

(A) 1-A.3 Steady-state current/voltage curves show that net outward currents grow in three age bins (P1-P5; P6-9; and P10-P16). Small and large currents are evident in all age ranges. (B,D and F) Log-normal plots showing the distribution of maximum conductance (gmax), current threshold (Ithreshold) and latency in three age bins, respectively. Points are color coded to indicate normalized basal position (red to blue gradient), unlabeled cells (gray) and type II neurons (black). Within each age bin, biophysical properties of type I SGN change as a gradient along the NBP scale. The relative dispersion in each distribution is quantified by the CV (coefficient of variation indicated above each distribution) and is not systematically dependent on age. C,E and F Plots of same three features, gmax, current threshold, and latency plotted with age as a continuous variable but with spatial position binned as ‘modiolar-contacting’ (NBP >0; red symbols) and ‘pillar-contacting’ (NBP <= 0, blue symbols). Solid lines in C and E show the best regression fits through modiolar-contacting and pillar-contacting groups (gmaxmod=4.93Age+10;gmaxpillar=2.29Age+10). Solid lines in F plot curves defined by the inverting the regression plots in C; L(mod,pil)=420/gmax(mod,pil).

Morphological development of spiral ganglion terminal branches.

(A) Three-dimensional reconstructions of SGN terminals are made to quantify the number of branches at each terminal and the approximate diameter of the fibers. (B) The number of branches exponentially decreases over the course of development. One-to-one connections with the inner hair cell can be seen at P5. In some cases, terminals have multiple branches at P10. (C) Pillar-contacting SGN are not observed to have the mature one-connection morphology. (D) Fiber diameters do not differ between type I SGN along the NBP scale. Type II SGN have significantly smaller fiber diameters than type I SGN.


Table 1
Analysis of covariance.
FigureSourced.f.Sum of squaresF-scorep-value
3D LatencySpike
Spike * NBP14.140.130.71
4F gmaxSpike1150.991.370.25
Spike * NBP10.0380.00040.98
4 G g-30Spike141.271.970.17
Spike * NBP112.140.580.45
4 HV1/2Spike1632.643.530.068
Spike * NBP1107.860.600.44
4I τinactSpike118715.33.640.068
Spike * NBP19.850.00190.96
6C. gmaxClass (Modiolar/Pillar)11508.387.320.0107
Class* Age11232.785.980.0199
6E IthresholdClass (M/P)13827.398.870.0087
Class * Age1280.750.640.4334
6G LatencyClass (M/P)1577.73236.97<0.0001
Class * Age1131.2478.40.0066
Table 2
Correlations with normalized basal position.
Spiking neurons
Spiking plus non-spiking
Spiking plus non-spiking
Current-clamp features*IThreshold0.680.00070.810.00800.440.0772
*First-spike Latency−0.650.0013−0.760.0168−0.450.22
*Voltage Threshold−0.390.0169−0.610.08370.150.7390
*AHP Tau−0.610.0031−0.650.0217−0.210.5616
*AP Height0.310.170.030.56340.630.0069
Average Response Latency−0.89<0.0001−0.760.0010−0.770.0003
Resting Potential−0.540.0112−0.310.229−0.070.8048−0.480.0487
Voltage-clamp featuresgmax0.520.01660.580.01770.780.87990.620.0083
Table 3
Multiple variable regression models.
Table 4
Table of regression slopes.
FigureTermEstimateStandard errorTp-value
3D LatencySlope−22.84.02−5.66<0.0001
4F gmaxSlope27.872.5510.89<0.0001
4 Gg-30Slope17.460.8010.6<0.0001
4 HV1/2Slope−19.989.8−2.040.0495
4I τinactSlope178.962.32.870.0084
6C gmaxSlope2.801.062.630.0127
6E IthresholdSlope4.682.262.070.054
6G LatencySlope−1.7520.29−5.97<0.0001
Key resources table
Reagent type
(species) or resource
DesignationSource or referenceIdentifiersAdditional information
Strain, strain background (species)Long-Evans (rat)USC Vivarium and Charles RiverRRID:RGD_2308852Freshly isolated cochlear middle turns
Antibodyanti-myosin-VI (rabbit polyclonal)Proteus BiosciencesCat# 25–6791, RRID:AB_10013626IF(1:1000)
Antibodyanti-peripherin (rabbit polyclonal)MilliporeCat# 183 AB1530, RRID:AB_90725IF(1:500)
Antibodystreptavidin Alexa Fluor 488 conjugateMolecular ProbesCat# S32354, 182 RRID:AB_2315383IF(1:200)
AntibodyAlexa Fluor 594 (goat anti-rabbit, polyclonal)Thermo Fisher ScientificCat# A-11080, 185 RRID:AB_2534124IF(1:200)
OtherVecta-shieldVector LaboratoriesCat# H-191 1000, RRID:AB_2336789
Software, algorithmJMPSAS InstituteRRID:SCR_002865
Software, algorithmMatlabMathworksRRID:SCR_001622
Software, algorithmpClampMolecular DevicesRRID:SCR_011323
Software, algorithmOriginProOrigin LabsRRID:SCR_015636
Software, algorithmImarisOxford InstrumentsRRID:SCR_007370
Software, algorithmPrism 8Graph PadRRID:SCR_002798
Software, algorithmFilamentTracer in ImarisOxford InstrumentsRRID:SCR_007366

Data availability

Source data files are provided as an excel spreadsheet.

Additional files

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Download citations (links to download the citations from this article in formats compatible with various reference manager tools)

Open citations (links to open the citations from this article in various online reference manager services)