Abstract

The degree to which individual neuronal input-output functions are modulated by voltage fluctuations, however, is not well established across different cortical areas. Additionally, the extent and mechanisms of input-output modulation through fluctuations have been explored largely in simplified models of spike generation, and with limited consideration for the role of nonlinear and voltage-dependent membrane properties. To address these issues, we studied fluctuation-based modulation of input-output responses in medial entorhinal cortical (MEC) stellate cells of rats, which express strong subthreshold nonlinear membrane properties. Using in vitro recordings, dynamic clamp and modeling, we show that the modulation of input-output responses by random voltage fluctuations in stellate cells is significantly limited. In stellate cells, a voltage-dependent increase in membrane resistance at subthreshold voltages mediated by Na+ conductance activation limits the ability of fluctuations to elicit spikes. Similarly, in exponential leaky integrate-and-fire models using a shallow voltage-dependence for the exponential term that matches stellate cell membrane properties, a low degree of fluctu-ation-based modulation of input-output responses can be attained. These results demonstrate that fluctuation-based modulation of input-output responses is not a universal feature of neurons and can be significantly limited by subthreshold voltage-gated conductances.

Author Summary

It has been speculated that membrane voltage fluctuations in neurons play an important role in scaling the relationship between input amplitude and spike rate response. For this to be true, neuronal spike input-output behavior must be sensitive to physiological membrane voltage fluctuations. Using a combination of single cell recordings and modeling, we investigated the mechanisms through Which voltage fluctuations modulate neuronal input-output responses. We find that neurons that express an increase in membrane input resistance With depolarization show low levels of noise-mediated modulation of input-output responses due, in part, to voltage trajectories that suppress the likelihood of generating a spike in response to random current input fluctuations. Hence, nonlinear membrane properties arising from certain types of voltage-gated conductances limit noise-based modulation of neuronal input-output responses.

Introduction

Because probabilistic threshold crossings associated with fluctuations lower spike threshold, enabling spike response to otherwise subthreshold inputs [5,6] , it has been hypothesized that background activity amplifies neuronal sensitivity, and in doing so permits fluctuations to modify the input-output functions of neurons [7—12]. Consistent with this hypothesis, recordings in vivo often show a large variance in interspike intervals [13,14]. Spectral properties of voltage fluctuations are also correlated with different cognitive states, lending support to the idea that fluctuations play an important role in modulating spike output [4,15—17]. Finally, computational models suggest that neurons are sensitive to transient inputs and modulate their input-output function in response to changes in the size of membrane voltage fluctuations [10,18—20].

First, data supporting a strong relationship come from only a few types of neurons [8,11,21—23]. Second, even these restricted studies have shown significant variability in the magnitude of the effect [21,23—25]. These observations indicate a possible complex relationship between membrane voltage fluctuations and neuronal input-output modulation. Modulation of input-output responses is likely influenced by numerous factors, including subthreshold voltage-dependent properties present in neurons. For example, the negative slope conductance associated with Na+ current, which increases membrane resistance in close proximity to spike threshold [26] , has been shown to reduce neuronal responsiveness to high frequency voltage fluctuations in model neurons [27].

These neurons express strong nonlinear membrane properties at subthreshold voltages and are characterized by a voltage-dependent change in membrane resistance [28—30]. Like other cortical neurons, in vivo recordings of stellate cells have established the presence of large membrane voltage fluctuations that have the potential to influence input-output responses [31,32].

We find that nonlinear membrane properties associated with increased membrane resistivity at sub and peri-threshold voltages reduce fluctuation-based modulation of input-output responses. Overall, our results indicate that fluc-tuation-based modulation of neuronal input-output responses can be very low, with limited scaling of spike output via changes in noise and conductance levels.

Results

Stellate cell input-output functions are modulated weakly by membrane voltage fluctuations

These include the slope (gain) and rheobase of frequency-current (f-I) and spike-probability curves. Voltage fluctuations were generated using current-based fluctuations that were constructed using low-pass filtered white noise (see Methods). For each cell, we recorded a short trial period in which the current input fluctuation amplitude was adjusted to maintain a standard deviation (SD) in output voltage at rest (-75 mV, corrected for the electrodes junction potential) of approximately 2.5 mV (2.41 i 0.1 mV), a value commonly observed in vivo [33]. Controlling for the SD of voltage fluctuations was essential since the intrinsic properties of neurons are voltage-dependent and a fair comparison across different cells, conditions and models require that the SD of membrane voltage be constant. Furthermore, previous work addressing similar issues has controlled fluctuation sizes in terms of the SD of membrane voltage and used similar values [8,Q,Q,§,fi,fi,3j] , For f-I curves, spike frequency was determined using only the first three inter-spike intervals in order to avoid complications arising from the interaction between the time scale of voltage fluctuations and spike frequency adaptation [21,23]. Gain was calculated individually for each cell using the slope of a linear fit (r2: 0.64 to 0.96, mean: 0.89 i 0.02) to the f-I relationship, while rheobase was measured as the current required to elicit a minimum of 3 inter-spike intervals from a holding voltage of -75 mV.

The addition of fluctuations generated a small, but non-sig-nificant, leftward shift in rheobase (Fig 1Ci; 199 i 20 pA vs. 158 i 20 pA, p = 0.11, n = 20, 18). As with rheobase, the addition of voltage fluctuations did not generate a significant reduction in the f-I curve gain (Fig 1Cii; 0.161 1 0.10 spikes/ pA s. vs. 0.148 i 0.08 spikes/pA s, P = 0.31, n = 20, 18). To quantify potential changes in the f-I curve more carefully, we also measured the effects of voltage fluctuations on firing rate within discrete regions of the f-I curve (Fig 1D; low, mid and high). Previous modeling and experimental work has shown that random voltage fluctuations induce the largest increase in firing rate in the low spike rate region of the f-I curve, near the transition between rest and firing [8—10,19,20,24]. For our data, the low region was defined individually for each cell as the frequency of spike discharge at rheobase, while the mid and high regions corresponded to current values eliciting 15 spikes/s and 25 spikes/s more than the initial frequency, respectively. For each cell, we measured the change in firing rate brought about by voltage fluctuations for low, mid and high current input regions relative to the same cell’sf-I curve acquired without fluctuations. Differences in spike rate were small but changed significantly in the low and mid regions of the f-I curve. For the low region of the f-I curve, fluctuations induced an increase of 3.3 i 0.34 spikes/s (Fig 1D; P <0.001, n = 18), while in the mid region, these values were 3.2 i 1.4 spikes/s (P = 0.03, n = 18).

Unlike the f-I curve, which requires repetitive spike generation at each current step size, a spike-probabili-ty curve quantifies the probability of generating just a single spike within a given time window (100 ms in our case, Fig 1E and 1F). In the absence of membrane voltage fluctuations, the transition from a probability of zero to one in spike discharge occurred almost always within a single current step (Fig 1F ; vertical line). For experiments without artificial voltage fluctuations, the current amplitude associated with this transition point was defined as rheobase. The addition of voltage fluctuations smoothed the relationship between spike probabilities and current steps such that data points could be fit with a sigmoid function (Fig 1F; Boltzmann fit, r2 >096). In the presence of voltage fluctuations, rheobase was defined as the current needed to establish a 0.2 probability (P02) in generating a single spike, while the slope was quantified using the “k” term in the Boltzmann function (see Methods). As shown, voltage fluctuations induced a leftward shift in rheobase relative to conditions without fluctuations (Fig 1F; 56 i5.6 pA, P <0.001, Student t-test, n = 20) and resulted in spike-probability curves with an average slope of 26 i 8.4 pA (mean i s.e.m).

In comparison, past work using similarly sized voltage fluctuations has reported reductions in gain of up to 50% [8,9]. Nevertheless, stellate cells do show some degree of fluctuation-mediated modulation of input-output responses as indicated by a shiftand smoothing of the spike-probability curve.

Stellate cells express significant nonlinear membrane properties leading up to spike threshold

As a potential cause for the limited fluctuation-based modulation of input-output responses, we considered the role of nonlinear membrane properties leading up to spike threshold. In simple models, realistic spike generation dynamics have been shown to reduce the likelihood of spike response to rapid voltage fluctuations [27,36]. We hypothesized that an extension of this effect over a much larger subthreshold voltage region than has been previously considered could significantly reduce fluctu-ation-based modulation of input-output responses.

At each holding voltage, membrane resistance was measured in voltage-clamp using a 5 mV voltage step of 100 ms duration. Depolarizing stellate cells led to a progressive increase in steady-state membrane input resistance (Fig 2A; one-way ANOVA, P <0.001, n = 12). Membrane resistance nearly tripled over a 20 mV range, increasing from 51.8 i 3.9 M!) at -85 mV to 151.4 i 15.5 M!) at -65 mV (Fig 2A). We should note that resistance also kept increasing with additional depolarization to the extent that very small voltage steps (~ 1 mV) often elicited spikes at levels more depolarized than -65 mV and prevented an accurate measure of membrane resistance at these voltage values.

For current inputs eliciting small changes in voltage (5 mV), the resulting voltage trajectory was fit accurately with an exponential function and used to extract the membrane time constant near -75 mV (12.0 i 0.9 ms, n = 19). In contrast, the membrane voltage trajectory to the first spike (for a 50 ms first spike latency) from a holding voltage of -75 mV (~ 15 mV range) could not be fit with an exponential function due to the more linear profile of the trajectory (Fig 2B). Similarly, the average voltage trajectory between the trough of the afterhyperpolarization (AHP) and spike threshold during continuous f1ring (~ 4 Hz) was not exponential (Fig 2C). Thus, the voltage trajectories to spike threshold starting either from resting voltages or the trough of the AHP were relatively linear compared to that expected from our measures of the membrane time constant at—75 mV. To assess how changes in spike rate relate to changes in the average voltage associated with spike trajectories at different spike frequencies, we quantified the relationship between spike Am 9‘“

a co g-n —stellate v,, trajectory

Stellate cells express nonlinear membrane properties that shape the voltage trajectory leading up to spike threshold. (A) Average steady-state input resistance as a function of membrane voltage. Measures were taken in voltage-clamp using a 100 ms duration 5 mV step. (B) Average membrane voltage trajectory (grey lines indicate SEM) associated with approach from rest to first spike. For comparison, the exponential approximation using the membrane time constant measures taken at -75 mV is also shown. (C) Average membrane voltage trajectory (grey lines indicate SEM) associated with the interspike interval. (D) Representative example of an f-V curve from a stellate cell (solid line). Line indicates fit to a power-law function of the form shown in panel inset. Inset also shows box-plot of mean exponent (p) value forthe power-Iaw function fit.

As With the f-I curves, spike rate for the f-V curve was taken from the first three inter-spike intervals, With membrane voltage values calculated as the mean during the same period of time. Surprisingly, we found that f-V curves were nonlinear and could be fit With a power-law function (Fig 2D; mean r2 = 0.95 i 0.02, range: 07—099) using an exponent (p) of 1.45 i 0.08 (Fig 2D; n = 19, range: 0.45—2.0, 17/19 had p values >1). Cell f-V curves were also shallow, With an average slope across the firing range of 4.5 i 0.4 spikes/mV s (n = 19). Contrary to previous assumptions [12,37], therefore, neuronal f-V curves can eXpress significant power-law scaling in the absence of any fluctuation-based smoothing. In summary, our measures of membrane resistance, voltage trajectories, and f-V curves indicate that stellate cells eXpress significant sub and peri-threshold nonlinearities.

A gradual increase in membrane resistance is critical to reduced fluctuation-based modulation of input-output responses in an eLlF model

To understand the biophysical mechanisms and consequences of a voltage-dependent membrane resistance over a large region of subthreshold voltage, we started by studying the effects of fluctuation-based modulation of input-output behavior in a simplified model of spike generation in the form of an exponential leaky integrate-and-fire (eLIF) model [27]. This model has the advantage of incorporating important nonlinear membrane properties associated with spike threshold using a small set of parameters that are related to physiological measures (e.g. voltage-dependent subthreshold membrane resistance).

The exponential term models the increase in membrane resistance associated with an increase in Na+ conductance activation and spike generation. With small AT values, the I—V curve slope (membrane resistance) changes abruptly as the system approaches VT (Fig 3A). Conversely, with large AT values, the change in slope is more gradual. Thus, for a small AT value (e.g. 2 mV), the I -V curve is largely linear with the exception of a small voltage range in the immediate vicinity of VT. As AT becomes larger, however, membrane resistance increases gradually over a relatively large span of subthreshold membrane voltage values. We found that a AT value of 15 mV best matched experimental values of steady-state membrane input resistance observed in stellate cells (Fig 3B). We should note that a value 15 mV for AT is large compared to that implemented in previous work (0.8 m V to 6 mV) using an eLIF model to study cortical neurons [27,38,39] .

For comparison, we also included a completely passive model (i.e. standard leaky integrate-and fire model-LIF) consisting of a linear conductance term (15 n8, -75 mV reversal) and an artificial threshold (-55 mV). For all models, membrane voltage was reset (V,) to -65 mV after crossing threshold.

Likewise, for the interspike interval voltage trajectory between the AHP trough and spike threshold during repetitive spike discharge (~4 HZ), a AT value of 15 mV best captures the linear trajectory leading up to spike threshold observed in stellate cells (Fig 3Cii). As AT is reduced, the voltage trajectories approaching spike threshold become more exponential and qualitatively more like the passive model (Fig 3C).

For the passive model, the f-V curve is steep and has a negative slope. In the active models, both the slope and shape changes considerably with different AT values. As the value of AT is increased from 2 mV to 15 mV, the f-chrve slope goes from negative and steep to shallow and positive (Fig 3D). Furthermore, over the range of 1—60 spikes/ s, the f-chrve of the eLIF using a value of AT of 15 mV can be accurately fit with a power-law function with an exponent of 1.78, which is within the range of values observed experimentally for stellate cells (Fig 3F).

We delivered the same current input fluctuations as in stellate cells, and maintained voltage fluctuations with a SD of 2.5 mV at -75 mV. Previous studies have also established that increasing membrane conductance via a variety of mechanisms, which include shunting inhibition, balanced synaptic conductances or simply increasing membrane leak, facilitates modulation of the f-I curve using similarly sized voltage fluctuations [8—10,22]. In addition, as a result of the different slopes of the exponential terms, the eLIF models using AT values of 2 mV and 15 mV have different input resistance values at -75 mV (Fig 3B), a property that may account for potential differences in fluctuation-based modulation of input-output responses. For these reasons, a separate shunt or leak conductance (g) of 15 n8 was introduced in each of the models to test the effect of increasing membrane conductance on input-output modulation. Fluctuation-based modulation of both the f-I and spike-probability curves are substantially larger using a value for AT of 2 mV than With 15 mV (Fig 4). Consistent With previous computational and experimental results [7—10], increasing membrane conductance using shunting inhibition increases fluctuation-induced modulation of input-output responses in both models, 0100200300400500000 100200300400500000 current amplitude (pA) current amplitude (pA)

(D) Changes in firing rate induced by membrane voltage fluctuations for low, mid and high regions of the H curves for each model. For comparison, average stellate cell values are also shown. (E) Plots of changes in rheobase (i) and sigmoid slope factor (ii) in spike-probability measures associated with the introduction of membrane voltage fluctuations for each model.

With AT 2 2 mV, fluctuation-induced increases in initial spike firing rates are 23 spikes/s and 30 spikes/s under baseline and with increased membrane conductance, respectively (Fig 4Bii—4D). In comparison, with AT 2 15 mV, these value are only 4.2 spikes/s and 5.8 spikes/s (Fig 4Bi—4D). Similarly, for spike-probability curves, voltage fluctuations result in larger leftward shifts and smoothing under both conductance conditions with AT 2 2 mV (Fig 4C—4E). As a result, the changes in firing rate, rheobase and gain induced by membrane voltage fluctuations correspond more closely to those observed in stellate cells when AT is set to 15 mV (Fig 4D and 4E).

Consequently, the increase in fluctuation-induced modulation of the f-I curves with lower AT values could result from an intrinsic higher gain value that provides greater sensitivity to changes in current input fluctuations. Prior work has established that decreasing the Vr value has a divisive effect on gain [40]. To eliminate the influence of higher gain values, therefore, we incrementally decreased the Vr value to decrease gain when AT was small (Fig 5A). Using this approach, we could maintain gain approximately equal across different AT values and test if differences in fluctuation-induced smoothing of the f-I curve are entirely due to changes in intrinsic gain (Fig 5A and 5B). Although compensating for gain through changes in V, decreases fluctuation-based modulation of the f-I curve, lower AT values still result in greater smoothing and increases in the initial firing rate of the f-I curve (Fig 5C).

For these measures the choice of V, has no impact. As shown, a gradual decrease in AT from 15 mV to 1 mV results in a gradual increase in both the ability for fluctuations to shift rheobase and smooth the spike-probability curves (Fig 5E).

Limited fluctuation-based modulation of input-output responses can be reproduced in a conductance-based model using Hodgkin and Huxley formulism

For our H-H-based model, we started with a non-inactivating Na+ conductance (INap) that generates the gradual increase in membrane input resistance at subthreshold membrane voltages observed in stellate cells (Fig 6A). For spiking currents, we used a standard transient Na+ current coupled with a slower K+ current (see Methods). Note that the rheobase values in the conductance-based model differed from the eLIF model because the voltage threshold for spiking was more depolarized than for the eLIF.

Further, the reduction of INap increases the slope of the f-chrve (Fig 6B). Consequently, the presence of IN,p drives the model neuron to generate an f-Vmore similar to the eLIF model using a AT 2 15 mV, while without INap, thef-V curve is more similar to when AT 2 2 mV.

As before, we added current input fluctuations so as to generate voltage fluctuations With a 2.5 mV STD at -75 mV. For both f-I (Fig 6C and 6D) and spike probability (Fig 6E and 6F) measures, the reduction of INap leads to a progressive increase in the ability for membrane voltage fluctuations to modulate the input-output responses. Thus, a standard H-H conductance-based model can reproduce experimental results and observations from the eLIF model.

Large AT values reduce modulation of input-output responses through voltage fluctuations by slowing membrane voltage

In the phase-plane plot representation, the dashed lines indicate the membrane derivative function (dV/dt). The farther the value of the dashed line is from the X-aXis at zero, the faster the membrane voltage changes. Through its effect on the shape of the dV/dt line, the AT parameter determines the rate at Which membrane voltage reaches spike threshold (Fig 7A and 7B). With AT 2 15 mV, the dV/dt function is shallow and generates small values, as indicated by the close proximity to the zero X-aXis, for a large portion of the trajectory leading up to spike threshold; this results in a trajectory that changes more slowly and spends a small fraction of time in close proximity to spike threshold (Fig 7A). Reducing AT increases dV/dt values for all voltages leading up to threshold. Under these conditions, the voltage trajectory only begins to slow in the immediate vicinity of spike threshold and, consequently, spends a large fraction of the interspike interval in close proximity to threshold (Fig 7A and 7B; insets).

As AT becomes smaller, voltage trajectories spend an increasing fraction of time near the voltage for spike threshold, whereby small fluctuations can lead to spike events very early in the evolution of the trajectory. Conversely, by linearizing the trajectory, large AT values limit fluctuation-induced spikes to time points later in the evolution of the voltage trajectory. To illustrate this, we quantified the probability that voltage fluctuations cause a spike at different time and voltage points during a single inter-spike interval trajectory (Fig 7C). Thus, each point along the trajectory provided initial conditions from which to calculate the likelihood of generating a spike at that given voltage and time point. From each of these points, we ran the models for a 50 ms period of time and calculated the probability of spike discharge in response to voltage fluctuations (SD = 2.5 mV) within this time period using 1000 trials (see Methods). As indicated in Fig 7C, for the majority of the interval, the likelihood of generating a spike in response to voltage fluctuations is substantially higher for AT 2 2 mV than that for AT 2 15 mV. Although a value of AT 2 15 mV leads to higher spike discharge probabilities towards the end of the interspike interval trajectory, the effect is limited in time and occurs late in the cycle, whereby the impact on spike rate relative to the deterministic case is small.

As indicated above, changes in AT alter the dV/dt line, with a larger value of AT generating a shallower dV/dt line. With the addition of positive (depolarizing) current, the dV/dt line moves upwards and away from X-aXis, which results in a faster voltage trajectory (i.e. shorter interspike intervals).

With AT 2 2 mV, the membrane derivative only slows down near the inflection point of the dV/dt line. Consequently, as the dV/dt line shifts upwards, the mean of the voltage trajectory remains largely unchanged because the majority of the trajectory is represented by the value near the inflection point of the dV/dt line (Fig 7E and 7F). In essence, a AT 2 2 mV compresses the voltage trajectory to a value approximately equal to VT_ With AT 2 15 mV, however, a shift upwards in the dV/dt line significantly accelerates the voltage variable for a much larger fraction of the trajectory profile, especially for time points early in the trajectory that are far from the inflection point of the dV/dt line. As a result, the mean of value of the voltage variable depolarizes as the applied current magnitude increases and the dV/dt line shifts upwards (Fig 7E and 7F ). By generating a shallow f- V, a large AT limits the ability for a change in voltage brought about through random fluctuations to increase spike firing rate. Conversely, when AT is small and the f-Vrelationship is steep, voltage fluctuations can give rise to a large change in spike firing rate (Fig 7D). With large AT values, the increase in membrane resistance, and the resulting increase in voltage fluctuations, could potentially overcome the reduction in fluctuation-based modulation established by a shallow f-V relationship. Although the SD of membrane voltage fluctuations is larger in the peri-threshold region with AT 2 15 mV compared to AT 2 2 mV, the difference is small compared to the changes in the f-V relationship established by increasing the value of AT

Thus, at the mean voltage Where fluctuations first induce spiking in both models (~-68 mV), the SD of voltage fluctuations only increases from 2.64 mV to 3.06 mV When AT is changed from 2 mV to 15 mV. In comparison, changes in the f-V relationship are much greater. With AT 2 2 mV, a 60 spikes/s range occurs entirely Within a 0.5 mV range, but increases to a 10 mV range When AT 2 15 mV (Fig 7D). As a result, increasing AT from 2 mV to 15 mV has a much larger effect on the scaling of the f-V relationship than the size of membrane voltage fluctuations. Even under conditions Where the SD of membrane voltage fluctuations is increased to 3.6 mV or 7.2 mV (by increasing the current-input fluctuations), values much larger than those used With AT 2 2 mV, the modulation of the f-V and f-I relationships are still significantly less With AT 2 15 mV than With AT 2 2 mV (Fig 8B and 8C). In summary, a large AT generates voltage trajectories that spend a smaller fraction of time in close proximity to spike threshold and Whose mean changes significantly With spike frequency. Both characteristics are the result of a gradual increase in membrane input resistance and help reduce modulation of input-output responses by voltage fluctuations.

Reducing voltage-dependence of membrane resistance reduces fluctuation-based modulation of input-output curves in stellate cells

Reducing the increase in input resistance associated with depolarization over the subthreshold region should lead to an increase in fluctuation-based modulation of input-output responses. This manipulation is akin to reducing the value of AT. Second, manipulations of membrane resistance using negative and positive sloped conductances should result in a reduction and increase, respectively, in fluctuation-based modulation of input-output responses, by manipulating the influence of the endogenous negative slope conductance and altering the voltage trajectories associated with the approach to spike threshold.

To establish the role of Na+ conductance in determining subthreshold input resistance in stellate cells, we used a small concentration of TTX (10 nM) that was able to significantly alter subthreshold membrane resistance but also maintain the ability to generate at least one spike. Application of TTX significantly reduced the gradual increase in input resistance across different voltages (2-way ANOVA, P <0.001,

Because the membrane voltage trajectory to the first spike from a holding voltage of -75 mV could not be fit with an exponential function (Fig 2B and Fig 9B), we quantif1ed changes in the first spike voltage trajectory by measuring the fraction of time spent above the midpoint of the trajectory. A trajectory with a large fraction above the midpoint is more similar to an exponential approach. The midpoint to threshold was calculated for each individual cell’s response to a square current step eliciting a ~50 ms latency to first spike. In the presence of TTX, the voltage trajectory spent a significantly larger fraction of time above the midpoint compared to control (Fig 9B; TTX: 0.84 i 0.03 vs. control: 0.68 i 0.01, P <0.001, n = 6—18). Unfortunately, 10 nM TTX eliminated the ability to generate continuous spike discharge that is required for f-I curve measures. As a consequence, we limited our analysis of cell output in the presence of TTX to spike-proba-bility curves. As indicated, membrane voltage fluctuations were more effective at shifting the rheobase and smoothing the spike-probability curve in the presence of TTX (Fig 9D and 9E). The leftward shift in rheobase increased from 64.4 i 7.7 pA to 134 i 25 pA (Fig 9Ei; paired Student t-test, P = 0.001, n = 10), while the slope factor increased from 18.0 i 1.3 pA to 29.6 i 4.7 pA (Fig 9Eii; paired Student t-test, P = 0.02, n = 10). These results indicate that reducing the amount of subthreshold Na+ conductance generates more exponential-like trajectories that result in increased fluctuation-based modulation of the spike-probability curves.

Manipulation of membrane conductance using dynamic clamp alters voltage trajectories and modulation of input-output responses by voltage fluctuations

To further investigate the role of membrane resistance and negative slope conductance in shaping modulation of input-output response by voltage fluctuations, we manipulated membrane resistance in stellate cells using dynamic clamp by introducing artificial negative or positive slope conductances. In particular, we were interested in the effects of negative slope conductance since the ability for Na+ conductance to increase input resistance and slow the rate of change of membrane voltage is related to the negative slope associated with the Na+ I -V relationship [26]. Hence, the introduction of an artificial negative slope conductance should further decrease fluctuation-based modulation of input-output responses. Conversely, the addition of a positive slope conductance should increase modulation via voltage fluctuations by decreasing the influence of the endogenously eXpressed negative slope conductance.

The positive conductance was set to 15 n8 and decreased membrane resistance to 34.1 i 2.7 M52. For the positive conductance, values greater than 15 n8 led to a loss of continuous spiking in fashion similar to what has been reported in CA1 pyramidal cells using this manipulation [42,43]. Both the negative and positive conductances were linear with a reversal potential of -75 mV. As with TTX experiments, we quantified changes in the first spike voltage trajectories using the fraction of time above the midpoint. Changes in membrane conductance had a significant impact on the trajectory leading up to spike threshold (Fig 10A; one-way ANOVA, P <0.001,

Changes in initial spike and interspike interval voltage trajectory using dynamic clamp. (A) Average membrane voltage trajectories for initial spike approach in stellate cells under -5 n8 (red), control (black) and 15 n8 of artificial membrane conductance added with dynamic clamp. Finer lines indicate sem. (B) Plot of average fraction above midpoint for initial spike voltage trajectories using -5 n8, control and 15 n8 levels of artificial membrane conductance. (C) Average membrane voltage trajectories in stellate cells using dynamic clamp with -5 n8 (red), control (black) and 15 n8 (blue) of artificial membrane conductance. (D) Plot of average AHP half duration forthe interspike interval voltage trajectories under -5 n8, control and 15 n8 of artificial membrane conductance. n = 8—18). Negative conductance decreased the fraction of time above the midpoint to 0.59 i 0.02, While adding positive conductance increased this this value to 0.82 i 0.02 (Fig 10B; P <0.001, Tukey’s test).

Negative conductance led to a significant increase in the AHP half-duration (time from trough to midpoint voltage between trough and spike threshold, Fig 10D; 98.1 i 5.3 ms, P <0.001, Tukey’s test). Although positive conductance did not significantly alter the AHP half duration relative to control When taking repeated measures into account, there was a decrease in mean values (43.4 i 4.5 ms, P = 0.07, Tukey’s test). Hence, overall, changes in membrane conductance altered the duration of membrane voltage trajectories leading to spike threshold in a form consistent With our analysis of the eLIF model.

For each conductance level, we compared changes in the slope and rheobase of f-I and spike probabilities curves induced by the introduction of membrane voltage fluctuations. As before, membrane fluctuations were kept at a SD of~ 2.5 mV (-5 n8, 2.3 i 0.1 mV, control: 2.41 i 0.1 mV; 15 n8: 2.35 i 0.01 mV, n = 5, 19, 18). Analysis of f-I curves indicated that gain was significantly modulated by changes in membrane conductance, but not by the introduction of membrane voltage fluctuations (Fig 11A and 11B; 2-Way ANOVA, P <0.001 for conductance, P = 0.36 for voltage fluctuations). Contrary to expectations [8,9,44,45] ,f-I curve gain can be modulated, albeit modestly, by changes in membrane conductance only. More importantly, stellate cells f-I curves maintained a low degree of fluctuation-based modulation at all three conductance levels. We should note that eLIF model using a AT 2 15 mV, but not AT 2 2 mV, also generates a small decrease in gain When membrane conductance is increased. This result is related to changes in mean interspike interval membrane voltage induced by changes in conductance When AT 2 15 mV, but not AT 2 2 mV. We proceeded to measure changes in rheobase and initial firing rate associated Withf-I curves resulting from voltage fluctuations under each of the conductance conditions. Changes in membrane conductance significantly impacted the ability of artificial fluctuations to change rheobase (Fig 11Ci; P <0.02, one-way ANOVA, P = 0.01, Tukey’s test). Further, voltage fluctuations were not able to significantly reduce rheobase values in the presence of -5 nS conductance (Fig 11Di; one sample Student t-test, P = 0.58), while with 15 nS rheobase significantly increased (Fig 11Di; one sample Student t-test, P <0.001). As with rheobase, changes in the initial spike firing rate were also significantly changed by membrane conductance (Fig 11Dii; P = 0.001, one-way ANOVA), with each conductance level generating significant differences in initial spike firing rates (Fig 11Cii; P <0.01, Tukey’s test). With -5 nS conductance voltage fluctuations were not able to significantly increase firing rate from zero (Fig 11Cii; one sample Student t-test, P = 0.79), while these changes were significant under control and with 15 nS (Fig 11Cii; one sample Student t-test, P <0.001).

For both rheobase and slope factors, conductance had a significant impact (Fig 11Ei and 11Eii; P <0.001, one-way ANOVA). Increasing membrane conductance by 15 nS led to both an increase in the ability for voltage fluctuations to shift the rheobase and smooth the spike-probability curves (Fig 11E). With a decrease in membrane conductance (-5 nS), rheobase changes and smoothing were less pronounced (Fig 11D and 11E).

Additionally, the linearization of voltage trajectories is fundamentally related to the negative slope conductance associated with Na+ current activation. Our data show that nonlinear membrane properties shape the potential for voltage fluctuations to modulate the input-output responses of neurons.

Discussion

First, stellate cell input-output responses are modulated to a small degree by membrane voltage fluctuations under control conditions, as well as with increased membrane conductance. Second, in both models and stellate cells, a voltage-dependent inward current resulting from a negative slope conductance that activates over a significant sub and peri-threshold region is responsible for low levels of fluctuation-mediated modulation of input-output responses. Finally, negative slope conductance reduced fluctuation-mediated changes in input-output response via a line-arization and slowing-down of membrane voltage trajectories that also established a shallow and non-linearf-V curve.

Controlling neuronal spike output with membrane voltage fluctuations and changes in membrane conductance

A generally accepted role of voltage fluctuations is to lower spike threshold and amplify weak inputs [5,6]. Thus, the size and spectral content of voltage fluctuations can be used to potentially gate inputs, as suggested by some network models [46]. In addition, fluctuating current input can divisively scale the f-I curve by disproportionately increasing spike firing rate near threshold [5,7—10]. Moreover, under fluctuation-driven spiking, it has been shown that spike firing rate can scale with subthreshold membrane voltage for a meaningful portion of a cell’s dynamic range [9,10]. For this reason, in the presence of significant voltage fluctuations, changes in the subthreshold I—V relationship brought about through shunting inhibition, balancing excitatory and inhibitory conductances, or simply introducing a leak conductance are expected to translate into a further reduction in the slope of the f-I curve. Changes in the slope of the f-I curve are believed to be critical for setting the tuning curve of individual cells and have been proposed to play a critical role in sensory processing, particularly in the visual system With regards to setting the neuronal spike response to contrast [3,12,37,45]. Synaptic-mediated voltage fluctuations have therefore been implicated in setting both neuronal spike threshold and the overall scaling of the input-output relationship.

In granule cells, for example, fluctuations permit spike generation over a range of more than 100 spikes/s across What would otherWise be sub-thresh-old voltage values in the absence of fluctuations. Our data, however, suggest that this scenario is not generalizable, at least Within the limits of physiological voltage fluctuations. In the case of stellate cells and our eLIF model (AT 2 15 mV), spike frequency scaling With subthreshold membrane voltage is severely constrained due to a very shallow f-V relationship, Which limits the spike frequency range in Which voltage fluctuations can generate spikes. Consistent With our interpretation, granule cells have very steep f-V curves in the absence of significant membrane voltage fluctuations, With a 250 spikes/s range occurring over less than 2.5 mV [9]. Similarly, simulations supporting this mechanism [10] have been carried out in a compartmental model also eXpressing a steep f-Vrelationship (80 spikes/s over ~ 2.5 mV). In contrast, stellate cells and the eLIF model (AT 2 15 mV) generate f-V curves With slopes in the range of 4—5 spikes/mVs. Results presented here, therefore, indicate that the characteristics of the f-V curve and its relation to voltage trajectories leading up to spikes are crucial to understanding the degree that neuronal input-output functions are modulated by voltage fluctuations.

Power-law scaling of stellate cell f-V curve without voltage fluctuations

Modeling has shown that a power-law scaling with an exponent near 2 between spike firing rate and voltage can arise from the combination of an intrinsic, steep and linear f-Vrelationship, and smoothing through Gaussian-distributed voltage fluctuations [3,19,34,47]. In contrast with past assumptions, our data indicate that the f-V curve is not well approximated by a steep linear function, and that part of the nonlinear scaling between spike rate and voltage can result from intrinsic voltage-dependent membrane properties.

By activating gradually, depolarization results in a change of mean voltage at different spike discharge rates. This is because the rate of change in voltage is heavily influenced by the activation of the negative slope conductance. As greater amounts of the negative slope conductance are activated incrementally, the shape and mean of the interspike-interval voltage trajectory change considerably. This is not the case with AT 2 2 mV because the voltage trajectory at different frequencies is largely set by the passive properties of the membrane, which result in an exponential approach that does not experience a change in mean with increasing levels of depolarization and firing rate.

Implications for MEC network function

Stellate cells have been implicated in spatial navigation Via their grid-like spatial firing fields [48]. It is possible that the low degree of modulation of input-output responses by membrane voltage fluctuations in stellate cells helps maintain stable firing patterns With respect to spatial position. By reducing the influence of rapid voltage fluctuations on firing rate and input-output responses, stellate cell behavior likely promotes increased reliability to inputs associated directly With relevant network activity [31,32]. Consistent With this interpretation, our past work on spike-phase locking in stellate cells demonstrated a very high degree of spike-phase locking to slow (1—10 HZ) oscillatory inputs in the presence of random voltage fluctuations and increased membrane conductance [49].

Role of Na+ conductance and threshold dynamics in determining input-output modulation by voltage fluctuations

By slowing the approach to spike threshold, the gradual activation of Na+ conductance results in membrane voltage being farther from threshold for a large fraction of the trajectory. Current injections associated with f-I curve measures lead to a graded change in mean voltage and small, incremental increases in firing rate, resulting in a shallow f-Vrelationship. Both an LIF model and an eLIF model implemented with small AT values eXpress linear subthreshold membrane properties. This leads to an eXponential approach to threshold in which membrane voltage plateaus early and a large fraction of time is spent near spike threshold. As a result, small changes in current input and membrane voltage result in large changes in firing rate. In addition to stellate cells, cortical neurons in the visual corteX [39], cerebellar Purkinje cells [50] and striatal interneurons [51] show a gradual increase in membrane resistance with depolarization that is mediated by Na+ current. Hence, the behavior observed in stellate cells is likely applicable to a wide range of different neurons.

This bifurcation is often associated with type I characteristics present in cortical pyramidal cells [52]. On the other hand, fast-firing interneurons in cortex have been classified as type II, with threshold behavior often modeled using a Hopf bifurcation and hence not requiring an increase in membrane resistance [52,53]. Previous modeling and experimental studies have suggested that type I behavior promotes a high sensitivity to membrane voltage fluctuations [18,52,54]. Unfortunately, drawing a clear relationship between the degree of modulation of input-output responses by voltage fluctuations and the type of threshold bifurcation is difficult. Our mechanism requires a graded increase in membrane resistivity over a 20 mV range, while a determination between a saddle-node and a Hopf bifurcation is established in the immediate vicinity of spike threshold, which is often less than 1 mV. In vivo voltage fluctuations, however, can span more than a 10 mV range such that the integration behavior of a cell over large regions of subthreshold voltage, which are well outside the immediate vicinity of spike threshold, become crucial to understanding how a cell reacts to voltage fluctuations. For these reasons, we believe that either form of bifurcation and type can give rise to low or high sensitivity in input-output responses to random voltage fluctuations.

Materials and Methods

Ethics statement All experimental protocols were approved by the University of Utah Institutional Animal Care and Use Committee.

Tissue preparation

All chemicals were obtained from Sigma-Aldrich (St. Louis, MO) unless otherwise noted. After anesthetization with isoflurane and decapitation, brains were removed and immersed in 0°C artificial cerebrospinal fluid (ACSF) solution consisting of the following: (in mM): NaCl (125), NaHCO3 (25), D-glucose (25), KCl (2), CaClz (2), NaHzPO4 (1.25), MgC12 (1), and buffered to pH 7.4 with 95/ 5% 02/C02 gas. Horizontal slices were cut to a thickness of 400 um (Leica VT 1200, Leica Microsystems; Wetzlar, Germany). After the cutting procedure, slices were incubated in ACSF at 30°C for 20 minutes before being cooled to room temperature (20°C). After the incubation period, slices were moved to the stage of an infrared differential interference contrast-equipped microscope (AXioscope 2+; Zeiss, Oberkochen, Germany). All recordings were conducted between 32 and 34°C.

Electrophysiology

Stellate cell identity was established using the following criteria: 1) presence of a hyperpolar-ization-mediated membrane voltage sag, 2) impedance and resonance measures indicating a steady-state input resistance at rest between 35 MS) and 80 MS) and the presence of a ~5 HZ resonance peak determined using methods described previously [49] and 3) the location and cell morphology under DIC-IR optics (i.e. in layer II of MEC and with a non-pyramidal cell body shape). For voltage clamp experiments, we held cells at each corresponding voltage (-85 to -65 mV) and used a small step (5 mV) and measured the change in current. The ratio of the change in voltage and current was used to measure input resistance at each corresponding holding voltage. For dynamic clamp experiments, the current-clamp amplifier was driven by an analog signal from an x86 personal computer running Real-Time Application Interface Linux and Real-Time eXperimental Interface (RTXI) [57,58]. Shunting inhibition (link) was implemented using RTXI with the following equation:

were set to -75 mV and 15 n8, respectively. In the case of negative conductance (Fig 10 and 11), the ginh term was set to -5 n8. For all experiments, the sample rate of the dynamic clamp was set to 10 kHz. A measured junction potential of approximately 10 mV was subtracted from all recordings. Data were collected at 10 kHz and filtered at 3 kHz. Current input fluctuations were implemented with filtered white noise using a low pass implement an inverse Fourier transform and generate the time series from signals constructed in the frequency domain. For each cell, we recorded a short trial period in which the current input fluctuations were adjusted to maintain a standard deviation (SD) in membrane voltage fluctuations at rest (-75 mV) of ~25 mV. For Na+ channel block, tetrodotoxin (TTX, Tocris, Bristol, UK) was bath applied at a concentration of 10 nM. Recordings with TTX were carried out approximately 15 minutes after bath application of the drug.

Simulations and models

For the exponential leaky integrate-and-f1re (eLIF) model [27] , membrane voltage dynamics were governed by the following differential equation: where C = 170 pF, VT 2 -60 mV, gL = 25 n8, EL 2 -75 mV and AT 2 15 mV (default). Note, for the passive version of the model (i.e. standard leaky integrate-and-fire), the first gL term was set to zero and a separate leak term was used with the same reversal potential and using a conductance value of 15 n8. This value of conductance generates a passive model with the same input resistance as the eLIF model at -75 mV. Because of the eXponential term in the eLIF, membrane voltage diverges to infinity upon crossing VT. For the eLIF simulations, membrane voltage was reset to VR (-65 mV) upon reaching a value of 0 mV. The passive model lacks a true threshold phenomena, therefore an artificial threshold was set at -55 mV and membrane voltage was reset to -65 mV after crossing this threshold value. All models were simulated in Matlab and solved with a forward Euler method using a time step of 0.01 ms. We also tested model solutions with a time step of 0.001 ms and found the same results. For the H-H formulism-based model (Fig 6), membrane voltage was governed by the following equations: As With experiments, current fluctuations were generated using filtered White noise generated using the same equation and cutoff frequency (100 HZ) as in experiments. To ensure a SD

Current fluctuations were added to the DC current term (16) in the above equation (i.e. additive noise). To increase model membrane conductance, a separate therm was added using a reversal potential of -75 mV.

Analysis and statistics

For power-law and Boltzmann fits, we also confirmed fits in Origin 8.5 (OriginLab, Northampton, MA). Spike times were determined using a threshold crossing for membrane voltage. Spike frequency was determined using the mean inverse of the first three inter-spike intervals elicited from one or two second current steps. Average gain values were determined by averaging the individual slope values attained using a linear regression analysis of the f-I relationship. Rheo-base values were calculated as the minimal current required to elicit 4 spikes from a holding voltage of -75 mV. For illustrative purposes concerning the average f-I curves shown in Fig 1 and 8, we calibrated the starting point of the f-I curve such that each cell’s initial value was near the average rheobase calculated for a given condition. For spike-probability curves, the current step duration was 100 ms and each current step size was repeated 15 to 25 times, with spike-probability defined as the number of steps that evoked spikes divided by the total number of steps. Individual spike-probability (P(I)) curves were fit with a Boltzmann function as follows: where P is the probability of spike discharge,Ihalfis the current step size value required to elicit 0.5 probability (P05) in spike discharge and k is the slope (larger values denote a shallower slope) of the curve. Fits were used to calculate the slope factor (k), while the rheobase for probability curves was defined as the current step size required to elicit a P02 in spike discharge. For the f-V curve, experimental and modeling results were fit using a power-law function: where f is the firing rate, p is the exponent of the fit reported in the results section, a and b are positive constants and VC is the minimal voltage required to elicit spike generation. Note that the term a was bounded such that only values of 1 or greater were possible. All fits used a least-squares method.

At every given fraction of the interval, we took the deterministic state of the model (value for voltage and all other variables) and used it as the initial condition for a set of 1000 test simulations, each 50 ms long in which fluctuations were then added. We defined probability as the number of simulations generating a spike divided by the total number of simulations at each given interval fraction. Thus, our probability measures how likely it is that voltage fluctuations will introduce one spike at a given phase and provides a measure for the ability of voltage fluctuations to change the spike rate associated with the low current region of the f—I curves. For multiple comparisons, statistical significance was determined using either a one-way or two-way ANOVA. For repeated measures of means, statistical difference was determined using Tukey’s honestly significant criteria, while a Student t-test (one or two-sample) was used for comparison of one or two values. Means are presented along With the standard error of the mean (s.e.m. ).

Acknowledgments

The authors would like to thank Michael N. Economo, W. Hamish Mehaffey and Richard D. Rabbitt for extensive discussions and comments on the manuscript.

Author Contributions

Performed the experiments: FRF. Analyzed the data: FRF PM. Contributed reagents/materials/ analysis tools: FRF PM. Wrote the paper: FRF.