Axonal Na+ channels detect and transmit levels of input synchrony in local brain circuits

See allHide authors and affiliations

Science Advances  06 May 2020:
Vol. 6, no. 19, eaay4313
DOI: 10.1126/sciadv.aay4313


Sensory processing requires mechanisms of fast coincidence detection to discriminate synchronous from asynchronous inputs. Spike threshold adaptation enables such a discrimination but is ineffective in transmitting this information to the network. We show here that presynaptic axonal sodium channels read and transmit precise levels of input synchrony to the postsynaptic cell by modulating the presynaptic action potential (AP) amplitude. As a consequence, synaptic transmission is facilitated at cortical synapses when the presynaptic spike is produced by synchronous inputs. Using dual soma-axon recordings, imaging, and modeling, we show that this facilitation results from enhanced AP amplitude in the axon due to minimized inactivation of axonal sodium channels. Quantifying local circuit activity and using network modeling, we found that spikes induced by synchronous inputs produced a larger effect on network activity than spikes induced by asynchronous inputs. Therefore, this input synchrony–dependent facilitation may constitute a powerful mechanism, regulating synaptic transmission at proximal synapses.


Higher brain function such as sensory processing or encoding of memory requires biological mechanisms to be able to discriminate highly synchronized inputs from weakly synchronized or asynchronous inputs with a temporal accuracy in the millisecond range (14). This property is classically thought to be achieved by mechanisms of coincidence detection occurring at the postsynaptic side such as activation of NMDA receptors at glutamatergic synapses (5), dendritic amplification of synaptic responses (6), and spike threshold dynamics (7). Spike threshold adaptation represents the most accurate mechanism known for detection of fast input synchrony in central neurons (710). The spike threshold recorded in cortical neurons in vivo depends on the rate of depolarization and the nature of the stimulus (7, 11). For example, in the barrel cortex, the spike threshold is hyperpolarized for the preferred stimulus and depolarized for the nonpreferred stimulus (11). Theoretical work has underlined the importance of Na+ channel inactivation in the adaptation of spike threshold (9, 12). However, it is not known whether the readout of input synchrony by adaptation of spike threshold can be transmitted within the network. This amounts to asking the question: Does an action potential (AP) produced by synchronous inputs convey more information than a spike produced by asynchronous inputs? Axonal ion channels have been shown to determine presynaptic spike waveform and hence synaptic transmission in a context-dependent manner that corresponds to analog-digital modulation of synaptic strength (1318). We therefore examined whether levels of input synchrony could be read and transmitted through the network by axonal channels.

We show here that axonal sodium channels detect levels of input synchrony with millisecond precision and transmit this information to the network through modulation of AP amplitude and presynaptic release. Synaptic transmission is facilitated at cortical connections when the presynaptic AP is produced by synchronous inputs. Using dual soma-axon recordings, calcium imaging, and computer modeling, we show that this facilitation results from the enhancement of axonal spike amplitude caused by minimized inactivation of axonal sodium channels. In addition, we reveal that synchronous input-mediated spikes provoke a larger increase in local circuit activity than asynchronous input-mediated spikes. This suggests that this form of context-dependent enhancement of synaptic transmission may represent a robust way for information coding.


Input synchrony–dependent facilitation

Pairs of monosynaptically connected L5 pyramidal neurons were recorded in acute slices of P13-49 rats. The effect of synchronous activation on synaptic transmission was tested by injecting AMPA-like excitatory presynaptic potentials (EPreSPs) into the presynaptic neuron using dynamic clamp. Postsynaptic responses were found to be larger when the presynaptic AP was produced by highly synchronized inputs (0-ms delay between inputs) than when elicited by intermediate (5-ms delay) or weakly synchronous (10-ms delay) inputs (140 ± 14%, n = 8; Wilcoxon test, P < 0.01; Friedman-test, P < 0.01; Fig. 1A). Therefore, the synchrony level of the inputs determines output synaptic strength. We named this phenomenon input synchrony–dependent facilitation (ISF). As previously reported, the AP threshold of the presynaptic neuron was found to be lower as EPreSP synchrony increased (fig. S1).

Fig. 1 ISF at L5-L5 synapses.

(A) Induction of ISF with physiological stimuli. Left: Schematic representation of the system used to inject a dynamic current mimicking five glutamatergic inputs into the presynaptic neuron. Middle: Examples of electrophysiological recordings from a connected pair of L5 neurons. The five EPreSGs were injected with a delay of 10 ms (black traces), 5 ms (orange traces), or synchronously (red traces). Note the increase in the postsynaptic response amplitude with the increase in EPreSGs synchrony. Resting membrane potential of the presynaptic neuron: −78 mV. Right: Statistics (open circles, individual pairs; closed circles, averages). Willcoxon test, **P < 0.01. (B) Induction of ISF with current pulses. The presynaptic AP was induced by a brief step of current with or without preceding 50-ms prepulse to mimic weakly synchronized or highly synchronized inputs, respectively. Note the increased excitatory postsynaptic current (EPSC) amplitude induced by a presynaptic AP induced by a synchronous-like input (si-AP) compared to that produced by an asynchronous-like input (ai-AP). Resting membrane potential of the presynaptic neuron: −77 mV. Right: Plot of ISF observed in 33 cell pairs. Wilcoxon test, ***P < 0.001.

ISF could also be induced by a simple protocol based on two ways of inducing the presynaptic AP (i.e., directly from resting potential to mimic synchronous inputs or after a 50-ms depolarizing prepulse to mimic asynchronous inputs). Each stimulation type was evoked in alternation. Consistent with our preceding results (Fig. 1A), APs emitted by a synchronous-like input (si-APs) produced larger postsynaptic responses than APs emitted by asynchronous-like input (ai-APs) (117 ± 3%, n = 33; Wilcoxon test, P < 0.001; Fig. 1B). During ISF, the EPSC2/EPSC1 ratio diminished (from 64 ± 4% to 47 ± 3%, n = 24; Wilcoxon test, P < 0.001; Fig. 2A), and the CV−2 of EPSC amplitude increased (181 ± 22% of the control, n = 26; Wilcoxon test, P < 0.001; Fig. 2B), indicating the presynaptic origin of ISF. Magnitude of ISF was not related to age (fig. S2A) but followed an inverse function of synaptic strength (fig. S2B). ISF was also observed at local excitatory synapses in area CA3 of the hippocampus (117 ± 3% of the control EPSC, n = 20; fig. S2C), indicating that ISF may be a general feature in excitatory brain circuits.

Fig. 2 Properties of ISF.

(A) Reduced EPSC2/EPSC1 ratio after ISF. Left traces, paired-pulse ratio (PPR) with an asynchronous-like weak (black) and synchronous-like (red) input. Note the decrease in EPSC2/EPSC1 ratio (here from 76 to 59%). Right: Group data (Wilcoxon test, ***P < 0.001). (B) Reduced EPSC variability during ISF. Left: Examples of 30 EPSCs evoked by an asynchronous-like (back) or synchronous-like (red) input. Right: Plot of normalized CV−2 as a function of ISF (Wilcoxon test, ***P < 0.001). (C) Time course of ISF. Presynaptic spikes were evoked with prepulse duration (PD) ranging from 50 to 2 ms. Note that evoked postsynaptic responses increased when the PD decreased. Right: Plot of EPSC amplitude normalized to that produced by a PD of 50 ms (black line, exponential regression y = 39.6 * ex + 99.98 with τ = 5.1 ms; R2 = 0.3) Wilcoxon test, *P < 0.05, **P < 0.01.

We next determined the time course of ISF by measuring ISF magnitude as a function of the pre-pulse duration (PD). Presynaptic spikes were induced following depolarizations ranging between 50 and 2 ms (corresponding to input synchrony increase), and the resulting postsynaptic responses were recorded. ISF magnitude followed an exponential decay with increasing PD (time constant of 5.1 ms; Fig. 2C). Synaptic facilitation observed with a PD of 2 ms was found to be substantially larger than that observed for a PD of 5 ms (Wilcoxon test, P < 0.01). We conclude that ISF is highly time discriminative.

Mechanisms of ISF: Modulation of axonal spike amplitude

Since ISF has a presynaptic origin, we examined whether ISF was associated with an elevation in presynaptic Ca2+ influx. L5 neurons were recorded with a pipette filled with 50 μM Alexa Fluor 594 and 250 μM Fluo-4, and spike-evoked Ca2+ signals were measured in putative en passant boutons (Fig. 3A). The amplitude of the spike-evoked Ca2+ transient was increased when the spike was evoked by a synchronous-like input (144 ± 8%, n = 6, Wilcoxon test, P < 0.05; Fig. 3A). The ratio between si-AP–evoked Ca2+ entry and ai-AP–evoked Ca2+ entry declined with the distance from the soma with a space constant of 232 μm (Fig. 3B), a value compatible with previous estimation of axonal space constant (14, 15).

Fig. 3 Ca2+ imaging.

(A) Left: Ensemble view and detail of axonal arborization of an L5 pyramidal neuron filled with Alexa 594 and Fluo-4. The white arrow indicates the line scan made at ~200 μm from the cell body. Scale bars, 50 and 5 μm. Middle traces: Simultaneous electrophysiological (top) and Ca2+ imaging (bottom) recordings during ISF. Note the increase in Ca2+ fluorescence when the spike is evoked by a synchronous-like input compared with a spike evoked by an asynchronous-like input. Resting membrane potential: −76 mV. Right: Normalized amplitude of the spike-evoked Ca2+ signal (Wilcoxon test, *P < 0.05). (B) Attenuation of the enhanced spike-evoked Ca2+ influx produced by si-APs along the axon. Plot of the normalized amplitude of Ca2+ signal as a function of the distance from the soma. The data have been fitted with an exponential regression (y = 89.6 * ex + 100, with λ = 232 μm, R2 = 0.65).

To examine the presynaptic spike shape, L5 neurons were recorded simultaneously from the soma with a pipette filled with Alexa 488 (50 μM) to visualize the axon and from the axon cut end (i.e., bleb) at distances ranging from 50 to 191 μm (mean, 92 ± 19 μm; n = 9). The subthreshold depolarization attenuated along the axon according to an exponential decay with a space constant of 360 μm (fig. S3A). In accordance with the increased Ca2+ influx during ISF, amplitude of si-AP was found to be increased in the axon (by 12.4 ± 0.8 mV, n = 9; Wilcoxon test, P < 0.05; Fig. 4A) but not in the soma (difference, 1.0 ± 1.2 mV; n = 9; Wilcoxon test, P > 0.1). Moreover, the axonal spike amplitude followed an exponential decay as a function of PD that is similar to the decay of synaptic strength found in Fig. 1C (time constant of 4.0 ms; fig. S3B). However, the half-width of the axonal si-AP was not notably broader than the ai-AP (105 ± 3%, n = 8; Wilcoxon test, P > 0.1; fig. S3C). The spike threshold measured in the axon was hyperpolarized by ~4 mV during ISF (from −43.5 ± 1.8 to −48.9 ± 2.2 mV, n = 9; Wilcoxon test, P < 0.01; Fig. 4B), and the maximal slope of the rising phase of the spike was increased (from 210 ± 12 to 259 ± 15 mV/ms; n = 9; Wilcoxon test, P < 0.01; Fig. 4B), indicating that voltage-dependent sodium (Nav) channel inactivation is minimized when the spike is induced by synchronous inputs. Together, these data indicate that modulation of spike amplitude is at the origin of ISF through a Ca2+-dependent mechanism (19).

Fig. 4 Modulation of axonal spike amplitude.

(A) ISF is associated with a modulation of axonal spike amplitude. Left: Recording configuration. Simultaneous whole-cell recordings from the soma and cut end axon in L5 pyramidal cells were obtained. Middle: Comparison of somatic and axonal AP amplitudes evoked by asynchronous-like (black traces) and synchronous-like (red traces) inputs. Note the increased amplitude in the axon but not in the soma. Resting membrane potential: −78 mV. Right: Quantitative data (Wilcoxon test, *P < 0.05). (B) Hyperpolarization of spike threshold and increase in the slope of the spike rising-phase during ISF (Wilcoxon test, **P < 0.01). Left: Comparison of axonal spikes induced by asynchronous and synchronous inputs. Resting membrane potential: −80 mV. Middle: Phase-plot of the APs. Note the hyperpolarization of the threshold during and the increase of maximal slope during ISF. Right: Quantitative data.

Role of sodium channel inactivation

The differential sensitivity of the spike amplitude in the axonal and the somatic compartment to depolarization may result from the difference in Nav channel inactivation in these two compartments. Axonal Nav channels display a hyperpolarized inactivation curve (half inactivation of −80 mV) and are highly sensitive to inactivation by subthreshold depolarization (20). We used a NEURON model of a fully reconstructed L5 neuron (20) to determine the role of axonal Nav channel inactivation in ISF. Synaptic conductance was injected in the soma in the form of five inputs whose delay varied between 0 and 10 ms to simulate different input synchrony. Presynaptic Nav channel inactivation, spike-evoked Na+ current, and AP overshoot were measured under each condition in the presynaptic terminal located at 150 μm from the soma. Increased input synchrony leads to a decrease in Nav channel inactivation and therefore to larger Na+ current and spike overshoot (Fig. 5A). Next, we used the model of an L5 neuron connected to a postsynaptic cell by a synapse to determine the incidence of Nav channel availability on ISF. In agreement with the experiments, a si-AP had a larger amplitude and evoked a larger postsynaptic response than an ai-AP in control conditions [i.e., when half inactivation of axonal Nav channels was set to −80 mV (20); Fig. 5B]. The difference between si-AP and ai-AP amplitude was enhanced if inactivation of axonal Nav channels was increased by hyperpolarizing their basal inactivation curve (half inactivation at −85 mV; Fig. 5B), leading to enhanced ISF. We verified experimentally the model prediction by using carbamazepine (CBZ; 100 μM), an anti-epileptic drug that hyperpolarizes Nav half-inactivation by ~6 mV (21, 22). In L5 connections, addition of CBZ reduced synaptic transmission to 56 ± 9% of the control amplitude (n = 7; Wilcoxon test, P < 0.05; Fig. 5C) and increased the magnitude of ISF by 23 ± 8% (n = 7; Wilcoxon test, P < 0.05; Fig. 5C). Thus, hyperpolarizing Nav half-inactivation by CBZ enhances ISF possibly through the increase of the dynamic range of the axonal AP amplitude.

Fig. 5 Minimization of Nav channel inactivation.

(A) Reduced Nav channel inactivation during ISF in a model L5 pyramidal cell. Left: Voltage (top line), Nav channel inactivation (middle line), and presynaptic Na+ current (bottom line) traces for asynchronous (black), intermediate (orange), or synchronous (red) inputs. Note the increase in spike and Nav current amplitude and the decreased Nav channel inactivation measured before the spike. Resting membrane potential of the presynaptic bouton: −78 mV. Right: Plot of presynaptic Nav channel inactivation and presynaptic spike overshoot as a function of synaptic input synchrony. (B) Increasing presynaptic Nav channel inactivation enhances ISF in a Hodgkin-Huxley model of an L5-L5 connection. Top left: Control conditions (V1/2 inactiv = −80 mV). A presynaptic AP evoked by a synchronous input is 13.8 mV larger than that produced by an asynchronous input. This larger presynaptic spike produces a larger EPSP (ISF = 115%). Bottom left: Test conditions (V1/2 inactiv = −85 mV). Because Nav channel inactivation is larger, the modulation of presynaptic spike amplitude (20 mV) and ISF (217%) is enhanced. Resting membrane potential of the presynaptic bouton: −78 mV. Right: Plot of the amplitude of presynaptic spike overshoot as a function of Nav channel availability. (C) Experimental verification of the role of Nav channel inactivation during ISF. Top left: Example of L5 connection showing no ISF in control conditions. Bottom left: Same connection in the presence of carbamazepine (CBZ). The evoked EPSC is ~40% smaller and ISF is now clearly visible. Resting membrane potential of the presynaptic neuron: −83 mV. Middle: Plot of EPSC reduction caused by CBZ. Right: Plot of ISF in control (red) and in the presence of CBZ (light blue). Wilcoxon test, *P < 0.05.

To confirm the effect of changing Nav channel availability on ISF in the model, we next reduced the density of Nav channels in the axon from 300 to 150 pS/μm2 while keeping all other parameters constant. In this condition, ISF was enhanced from 116 to 165% (fig. S4A). We tested experimentally the second prediction of the model by partially blocking Nav channels with a low concentration of tetrodotoxin (TTX; 20 to 40 nM). TTX reduced synaptic transmission to 75 ± 9% of the control value (n = 7; fig. S4B) and increased ISF by 10 ± 2% (n = 7; Wilcoxon test, P < 0.05). We conclude that ISF results from the minimization of Nav channel inactivation by synchronous inputs that, in turn, leads to the increase in presynaptic AP amplitude, calcium influx, and glutamate release.

ISF induced in a single neuron affects network activity

Next, we tested whether ISF triggered in a single neuron had a substantial impact on spontaneous activity in cortical circuits. For this, we approached the physiological conditions of extracellular calcium [i.e., 1.3 mM external calcium instead of 3 mM (23, 24)]. First, we measured the effect of reducing external calcium on ISF. We found that decreasing extracellular calcium concentration from 3 to 1.3 mM decreased synaptic strength (to 40 ± 8% of the control EPSC, n = 8; Wilcoxon test, P < 0.01) but increased ISF amplitude (from 106 ± 4% to 118 ± 6%, n = 8; Wilcoxon test, P < 0.01; fig. S5, A and B). In a second step, paired recordings from adjacent L5 pyramidal neurons were obtained, and a spike was elicited in one of the two cells. The spontaneous activity was measured in the nonspiking cell by counting the number of postsynaptic events occurring 200 ms before the spike and 200 ms after the spike (Fig. 6A). To verify the integrity of the axon, neurons were filled with Alexa 488, and recordings from cells with damaged axons were rejected. We found that spontaneous activity in the nonspiking cell increased during the first 200 ms following a si-AP but not after an ai-AP (127 ± 6% versus 96 ± 5%, n = 6; Wilcoxon test, P < 0.03; Fig. 6B). We conclude that ISF induced in a single neuron substantially affects spontaneous activity in local cortical circuits, suggesting that transmission of input synchrony through ISF would represent a robust way for information coding.

Fig. 6 Effects of ISF on network activity.

(A) Incidence of ISF on cortical network activity. Spontaneous activity was measured in cell 2 before and after the spike triggered in cell 1 and the event frequency after spike was normalized to that measured 200 ms before the spike. Black traces, no change in network activity is induced by ai-APs. Red traces, transient increase in network activity induced by si-APs (127%). (B) Plots of normalized EPSC frequency triggered by ai-APs (black circles) and by si-APs (red circles) in different time windows (0 to 200, 200 to 400, and 400 to 600 ms). Note that the increase in network activity induced by si-APs is limited to the first 200-ms postspike. Wilcoxon test, *P < 0.05, n.s., not significant. (C) Effect of ISF in a balanced network model. (i) Averaged firing rate (over 100 repeats on 100 neurons and over five networks) just before and after the injection of an ai-AP spike (black curve) and a si-AP spike (red curve). (ii) Individual normalized changes in spiking activity, for the 100 neurons in a given network, when injecting an ai-AP (black) and a si-AP (red).

To confirm this result, we used a computational approach. We built a network model made of 4000 excitatory cells and 1000 inhibitory cells randomly connected with a probability of 10%. Neurons present a low spiking activity around 0.2 Hz (fig. S5C). The network was balanced for excitatory and inhibitory synaptic weights (see Materials and Methods for details; fig. S5, D to H). We reproduced the protocol used in the experiments to check the effect of ISF on the network. One excitatory neuron was randomly chosen and was forced to fire one extra spike. To model a si-AP, EPSPs induced by this extra spike were increased by 30% in target neurons. For an ai-AP, EPSPs were left unchanged. The spiking rate of the network was measured 100 ms before and 100 ms after the extra spike. This protocol was repeated for 100 ai-APs and 100 si-APs for a given neuron, and results over five random networks were averaged. Overall, we found that an ai-AP produced, on average, an increase of 8.7% of the network spiking activity, while a si-AP induced, on average, an increase of 25.2% (t test, P < 0.00001; Fig. 6C). This result confirms our experimental observation that ISF induced in a single neuron can affect local circuit activity.


We show here that sodium channels read and transmit levels of input synchrony with a high temporal precision to local cortical circuits through the modification of presynaptic AP amplitude. Synaptic transmission was found to be facilitated at L5-L5 or at CA3-CA3 connections when the presynaptic AP was triggered by highly synchronous inputs or by a short pulse of current that produced in all cases a steep depolarization from resting membrane potential. Using dual soma-axon recordings, calcium imaging, pharmacological tools, and computer modeling, we demonstrate that this facilitation results from an elevation in presynaptic calcium due to the enhancement of axonal spike amplitude caused by minimized inactivation of axonal Nav channels. This ISF was found to enhance local circuit activity, suggesting that it may represent a robust way of transmitting synchrony in the brain.

So far, spike threshold adaptation reported in cortical neurons was thought to represent the most accurate mechanism for detection of fast input synchrony in central neurons (710). We suggest here that it may represent only the first visible step of a longer process because spike threshold dynamics appears to be only an indicator of Nav channel inactivation in the presynaptic cell, and the consequences on the network can only be evaluated by paired intracellular recording from connected neurons, a particular challenging technique in vivo. The most important consequence of input synchrony–dependent adaptation in spike threshold is an enhancement of presynaptic release caused by modification of presynaptic spike amplitude. As a consequence, input synchrony increases both cell excitability by lowering the spike threshold and synaptic transmission via the enhancement of the spike amplitude in the axon.

An AP evoked by synchronous AMPA-like synaptic inputs (0-ms delay) produced a larger postsynaptic response than an AP evoked by asynchronous AMPA-like synaptic inputs (5- or 10-ms delay between inputs). We have not tested in this study the impact of dendritic input location on ISF. However, one can expect that asynchronous distal inputs, which present slower rise time and decay time, may maximize both the duration and magnitude of depolarization preceding spike emission. This should lead to a greater amount of Nav inactivation and a larger decrease of spike overshoot and synaptic release. Therefore, distal inputs may provoke larger ISF than proximal ones.

Synaptic change during ISF was found to vary from no change to +70% at L5-L5 connections. This variation can be attributed to the size of unitary connection as ISF follows an inverse function of EPSC size (fig. S2B). However, one cannot exclude the contribution of the axonal length that can also vary from connection to connection and that may attenuate propagation of voltage along the axon.

The contribution of axonal Nav channels in input synchrony–dependent change in presynaptic release is shown by four lines of evidence. First, whole-cell recording from the axon showed that the spike amplitude in the axonal compartment was highly modulated. In particular, the maximal slope of the rising phase of the spike was increased. Second, computer simulation showed that Nav channel inactivation caused by subthreshold depolarization preceding the spike decreased substantially when input synchrony increased. Third, enhancement of Nav channel inactivation by CBZ was found to enhance the magnitude of ISF. Last, reduction of Nav channel availability by TTX increased ISF. The implication of Nav channels had been shown in another context-dependent enhancement of synaptic transmission that depends on presynaptic hyperpolarization (17). This facilitation was due to recovery of axonal sodium channels from inactivation. Therefore, availability of presynaptic Nav channels represents a major determinant of synaptic transmission at local excitatory connections. ISF relies on inactivation of axonal Nav channels. Several studies report that axonal Nav channels display a more hyperpolarized inactivation curve than somatic Nav channels. Nav1.6 channels are present in the axon from the AIS to synaptic terminals, while Nav1.2 channels are present in the somato-dendritic region of the neuron including the proximal region of the AIS (20). Inactivation of sodium channels found in the soma (i.e., Nav1.2) is ~15 mV more depolarized than sodium channels found in the axon (Nav1.6). An accessory channel subunit of Nav channels, A-type fibroblast growth factor homologous factor (FHF), is specifically localized in the proximal part of the axon and could be responsible for this difference. FHF shifts Nav channel inactivation by ~15 mV in the depolarized direction (25, 26). The hyperpolarized inactivation curve of axonal Nav channels makes them more sensitive to subthreshold variations of membrane potential and therefore more suited to underlie ISF.

The modification of release we report here is due to a Nav channel–dependent enhancement of presynaptic AP amplitude and calcium influx. The fact that increasing presynaptic spike amplitude enhances calcium entry and neurotransmitter release has been shown in many axon types. For example, in cultured hippocampal neurons, voltage imaging of presynaptic boutons indicates that APs display a very small overshoot (27), a condition favorable for the modulation of Ca2+ entry by changes in AP amplitude. Moreover, voltage-clamp recordings from presynaptic boutons in cerebellar cultured neurons show that enhanced spike amplitude results in larger presynaptic Ca2+ currents and in enhanced release (19). In addition, partial blockade of Nav channels with low concentration of TTX reduces both AP amplitude and synaptic release in connected pairs of L5 pyramidal neurons (28). To optimize detection of small unitary synaptic currents, external calcium concentration was elevated to 3 mM under our control conditions. Under these conditions, release probability is already high, and upward modulations of axonal Nav spike and calcium entry may probably not produce a large increase in glutamate release. Supporting this view, when release probability was reduced by recording L5-L5 connection in 1.3 mM external calcium, ISF was found to be enhanced.

ISF enters into the group of context-dependent modifications of synaptic release that depends on axonal ion channels. To date, three following types have been identified: (i) a Kv1 channel–dependent form of synaptic facilitation induced by long presynaptic depolarizations preceding the presynaptic AP (14, 15, 29, 30), (ii) a Cav2.1 channel–dependent accumulation of presynaptic Ca2+ induced by slow presynaptic depolarizations that facilitates presynaptic release in short cerebellar axons (31, 32), and (iii) a Nav-dependent facilitation induced by a short presynaptic hyperpolarization (17). Among these context-dependent modifications of synaptic transmission, only the last one may have a substantial impact on the activity of fast cortical circuits, the two other processes being too slow and probably involved in slower network activity (14, 15, 29, 30). While both ISF and Kv1-dependent modulation of synaptic release are observed on the same synapses, their respective time windows are distinct (in the millisecond range for ISF and in the second range for Kv1-dependent modulation). In addition, Kv1-dependent facilitation induced by a long depolarization and Nav-dependent facilitation induced by transient hyperpolarization have been shown to be additive (17).

From a functional point of view, ISF can be seen as a time coding process in which APs elicited with a short latency would produce stronger responses than those induced with a long latency. Theoretical work suggests that coding based on absolute spike time is by far the most powerful because of its high capacity for information transfer compared to rate coding or rank coding (33). In addition, several experimental studies indicate that the precise timing of the first spike carries critical information about the stimulus in sensory responses (3437). Our data indicate that ISF conveys presynaptic spike timing information to the postsynaptic cell with high temporal precision and thus provides strong biological support for time coding based on spike latency at cortical synapses (33, 38). ISF defines a rule in which highly relevant neuronal information comes early and is transmitted in the circuit, while poorly relevant information comes later and is filtered. This rule is reminiscent of the principle governing the release-dependent variation in synaptic delay at unitary L5-L5 synapses where large synaptic responses display the shortest synaptic delays and small responses, the longest ones (39).

Does ISF also occur in vivo? Although a direct answer to this question is technically extremely challenging, we assume that this might be the case. Visually evoked inputs differentially inactivate Nav channels and therefore highly modulate spike threshold in primary visual cortical neurons (7). In the barrel cortex, preferred direction elicits steeper EPSPs that, in turn, evoke spikes with a lower threshold (11). Thus, spikes evoked by highly synchronous inputs (or fast rates of depolarization) are likely to better transmit information in the network than spikes evoked by weakly synchronous inputs (or slow rates of depolarization). Our experimental results show that this might be the case. Spikes elicited by synchronous-like inputs have a stronger impact on network activity compared to spikes elicited by asynchronous-like inputs.

The magnitude of ISF reported in the present study is relatively modest (20 to 40%). There are, however, good reasons to believe that ISF could be substantially larger in vivo. First, we show that ISF is notably enhanced in physiological extracellular Ca2+. Second, the presence of neuromodulators such as acetylcholine may promote ISF. Acetylcholine specifically depolarizes axon membrane potential by blocking Kv7 channels (40) and reduces sodium current (41), which should increase ISF (see TTX experiment, fig. S6). In addition, spike threshold adaptation (the prerequisite for ISF) has been reported in at least three independent in vivo studies (7, 10, 11). In addition, ISF might be larger at 37°C than at 30°C (experimental temperature) mainly because as kinetics is faster, sodium channels would be fully inactivated when an ai-AP is generated. Although ISF is present in mature cortical neurons (P42-P49), the contribution of ISF to network dynamics might be even larger during early development. The low density of Nav channels, together with the slow kinetics of EPSPs observed during development, represents two factors that may optimize ISF. Further experimental investigations will be needed to determine the importance of ISF in the cortex in vivo.


Acute slices of rat neocortex

Our experiments were conducted according to the European and Institutional guidelines [Council Directive 86/609/EEC and French National Research Council and approved by the local health authority (Préfecture des Bouches-du-Rhône, Marseille)]. Neocortical slices (350 to 400 μm) were obtained from 14- to 49-day-old Wistar rats of both sexes. Rats were deeply anesthetized with chloral hydrate (intraperitoneal, 200 mg/kg) and killed by decapitation. For animals older than P40, rats were deeply anesthetized with isoflurane and ice-cold perfused with the slicing solution (see composition below) and killed by decapitation. Slices were cut in an ice-cold solution containing 92 mM n-methyl-d-glutamine, 30 mM NaHCO3, 25 mM d-glucose, 10 mM MgCl2, 2.5 mM KCl, 0.5 mM CaCl2, 1.2 mM NaH2PO4, 20 mM Hepes, 5 mM sodium ascorbate, 2 mM thiourea, and 3 mM sodium pyruvate and were bubbled with 95% O2–5% CO2 (pH 7.4). Slices recovered (1 hour) in a solution containing 125 mM NaCl, 26 mM NaHCO3, 2 mM CaCl2, 2.5 mM KCl, 2 mM MgCl2, 0.8 mM NaH2PO4, and 10 mM d-glucose and were equilibrated with 95% O2–5% CO2. Each slice was transferred to a submerged chamber mounted on an upright microscope (Olympus BX51WI or Zeiss Axio-Examiner Z1), and neurons were visualized using differential interference contrast infrared videomicroscopy.

Organotypic slices of rat hippocampus

Hippocampal slice cultures were prepared using an interface technique (42). Briefly, postnatal day 5 to 7 Wistar rats were deeply anesthetized by intraperitoneal injection of chloral hydrate, the brain was quickly removed, and each hippocampus was individually dissected. Hippocampal slices (350 μm) were placed on 20-mm latex membranes (Millicell) inserted into 35-mm petri dishes containing 1 ml of culture medium and maintained for up to 21 days in an incubator at 34°C, 95% O2–5% CO2. The culture medium contained 25 ml of Eagle’s Minimum Essential Medium (MEM), 12.5 ml of Hank’s balanced salt solution, 12.5 ml of horse serum, 0.5 ml of penicillin/streptomycin, 0.8 ml of glucose (1 M), 0.1 ml of ascorbic acid (1 mg/ml), 0.4 ml of Hepes (1 M), 0.5 ml of B27, and 8.95 ml of sterile H2O. To avoid glial proliferation, 5 μM Ara-C was added to the medium at 3 days in vitro (DIV) for one night. Pyramidal cells from the CA3 area were recorded at the age of 8 to 12 DIV.

Paired recordings and data analysis

Dual recording from pairs of L5 neurons or CA3 neurons were obtained as previously described (42). The external saline contained 125 mM NaCl, 26 mM NaHCO3, 3 mM CaCl2 (1.3 mM in specific cases), 2.5 mM KCl, 2 MgCl2 (1 mM in specific cases), 0.8 mM NaH2PO4, and 10 mM d-glucose and was equilibrated with 95% O2–5% CO2. Synaptic inhibition was kept intact in both cortical and hippocampal recordings. Patch pipettes (5 to 10 megohm) were pulled from borosilicate glass and filled with an intracellular solution containing 120 mM K gluconate, 20 mM KCl, 10 mM Hepes, 0.5 mM EGTA, 2 mM MgCl2, 2 mM Na2ATP, and 0.3 mM NaGTP (pH 7.4). Recordings were performed with 2 Axopatch-200B (Axon Instruments, Molecular Devices) or a MultiClamp-700B (Molecular Devices) at 30°C in a temperature-controlled recording chamber (Luigs & Neumann, Ratingen, Germany). Usually, the presynaptic neuron was recorded in current clamp and the postsynaptic cell in voltage clamp. The membrane potential was corrected for the liquid junction potential (−13 mV). Both pre- and postsynaptic cells were held at their resting membrane potential (approximately −77 mV). Presynaptic APs were generated by injecting brief (3 nA for 2 ms) depolarizing pulses of current at a frequency of 0.1 Hz. For si-APs, the depolarizing pulse was evoked directly from the resting membrane potential. For ai-APs, the depolarizing pulse was evoked after a 50-ms depolarization of the presynaptic cell to ~−60/−55 mV. Synaptic transmission rundown can occur during the recordings. To prevent the potential rundown that could mask or artificially enhance ISF, si-APs and ai-APs were collected in alternation all along the recordings (si-AP, ai-AP, si-AP, ai-AP, …). EPSC2/EPSC1 ratio was assessed with two presynaptic stimulations (50-ms interval). Voltage and current signals were low pass–filtered (3 kHz), and sequences (200 to 500 ms) were acquired at 10 to 20 kHz with pClamp 10 (Axon Instruments, Molecular Devices). In all recordings, access resistance was fully compensated using bridge balance and capacitance neutralization (>70%). Electrophysiological signals were analyzed with ClampFit (Axon Instruments). Postsynaptic responses were averaged over 30 to 50 trials following alignment of the presynaptic APs using automatic peak detection. The presence or absence of a synaptic connection was determined on the basis of averages of 30 individual traces. Spike threshold was determined on AP phase plots as the voltage value exceeding 10 mV/ms.

Calcium imaging

Thin axon collaterals from L5 pyramidal neurons were imaged with a LSM-710 Zeiss confocal microscope. Putative boutons were identified as swellings on thin axon branches. For imaging calcium in the axon of L5 pyramidal neurons, 50 μM Alexa 594 and 250 μM Fluo-4 (Invitrogen) were added to the pipette solution, as previously described (17). Alexa 594 fluorescence was used to reveal neuronal morphology, whereas fluorescence signals emitted by Fluo-4 were used for calcium imaging. L5 neurons with a long axon in the plane of the slice were selected. Laser sources for fluorescence excitation were set at 488 nm for Fluo-4 and 543 nm for Alexa 594. Calcium signals were measured in the linear range of Fluo-4 detection. Emitted fluorescence was collected between 500 and 580 nm for Fluo-4 and between 620 and 750 nm for Alexa 594. After whole-cell access, the dyes were allowed to diffuse for at least 15 to 20 min before acquisition of fluorescence signals. APs were evoked by a brief (2 ms) and strong (2 to 3 nA) pulse of depolarizing current with or without a preceding subthreshold depolarization (50 ms, 50 to 100 pA). Electrophysiological signals were synchronized with calcium imaging in line scan mode (speed, 0.5 to 2 kHz). Acquired Fluo-4 signals were converted to ΔF/F values, and peak amplitudes were measured with a custom-made analysis program (LabView, National Instruments). Rundown of calcium signals may occur during imaging. To prevent the potential rundown that could mask or artificially enhance ISF, calcium signals resulting from si- and ai-APs were collected in alternation all along the recordings (si-AP, ai-AP, si-AP, ai-AP, etc.).

Axonal recordings

Simultaneous recordings from the soma and axonal cut ends (blebs, at 50 to 152 μm from the cell body) were obtained in L5 pyramidal neurons, as previously described (17). The axon was recorded 10 to 15 min after whole-cell access of the somatic compartment and the full labeling of the axon by Alexa 488 and visualized with the LSM-710 confocal microscope.

Network activity

To test the incidence of spiking activity on cortical network activity, paired recordings from adjacent L5 pyramidal neurons were obtained in physiological extracellular calcium conditions (1.3 mM). A spike was elicited in one of the two cells after having verified the integrity of the axon in the spiking neuron. For this, neurons were recorded with 50 μM Alexa 488, and only recordings in which the spiking cell displayed an axon longer than 100 μm were kept for the final analysis. The spontaneous excitatory postsynaptic activity was measured in the nonspiking cell by counting the number of postsynaptic events occurring 200 ms before and 200 ms after the spike with the use of a custom-made software programmed in LabView.


TTX and CBZ were provided by Tocris. All drugs were bath applied.

Dynamic clamp

Excitatory postsynaptic conductance (Fig. 1A) was injected into the presynaptic neuron using a dynamic clamp system (SM-1; Cambridge Conductance, Cambridge, UK), and a computer conductance profile was delivered via an analog-digital interface (Digidata 1440A; Molecular Devices). The conductance profile was derived from the Exp2syn Point Process from the NEURON 7.4 software (rise time, 1.5 ms; decay time, 10 ms). The reversal potential for AMPA receptor–like component was set at 0 mV. Amplitude of AMPA-like conductance was adjusted to produce in each case only one AP in the presynaptic neuron per episode of stimulation. In the synchronous configuration, the latency of the spike from the EPSP onset was 4.7 ± 0.7 ms (n = 8). In the intermediate configuration, the latency of the spike from the EPSP onset was 24.9 ± 1.3 ms (n = 8). In the asynchronous configuration, the latency of the spike from the EPSP onset was 43.7 ± 0.7 ms (n = 8).

Hodgkin-Huxley modeling

A multicompartment model of L5 pyramidal neuron was simulated with NEURON 7.4. The neuronal morphology was taken from Hu and coworkers (20). All simulations were run with 10-μs time steps, and the nominal temperature of simulation was 30°C. The voltage dependence of activation and inactivation of Hodgkin-Huxley–based conductance models were taken from Hu et al. (20) for gNav1.2 and gNav1.6, from Mainen et al. (43) for gKDR, from Shu et al. (44) for gKv1, and from Bischofberger et al. (45) for gCa. The equilibrium potentials for Na+, K+, Ca2+, and passive channels were set to +60, −90, +140, and −70 mV, respectively. The conductance density is provided in table S1, and the biophysics of Nav channels is provided in fig. S6. The resting membrane potential was set to −78 mV. The current was injected in the soma to produce APs. EPSPs were simulated in Fig. 5A with synaptic conductance whose waveform was determined by two exponential functions (EPSC rise time, 1.5 ms; and EPSC decay, 10 ms) and matched experimentally determined mEPSCs (46). In other modeling figures (Fig. 5B and fig. S4A), APs evoked by synchronous-like inputs were produced by injection of 1.6 nA during 2 ms directly from the resting membrane potential. APs evoked by asynchronous-like inputs were produced by injection of 1.6 nA during 2 ms after a 50-ms subthreshold depolarization via injection of 0.3 nA. To observe the postsynaptic response, a presynaptic terminal of the L5 neuron located at 150 μm from the soma and located on an axonal collateral was synaptically connected to a postsynaptic neuron. The postsynaptic neuron was a single-compartment model, containing only passive conductance. The level of release was dependent on the Ca2+ concentration in the presynaptic terminal. The synapse model was taken from Destexhe et al. (47).

Network modeling

Neuron model. Simulations were performed using current-based leaky integrate-and-fire neurons with membrane time constant τm = 20 ms and resting membrane potential Vrest = − 70 mV. When the membrane potential Vm reaches the spiking Vthresh = − 55 mV, a spike is generated and the membrane potential is clamped to the reset potential Vreset = − 70 mV during a refractory period of duration τref = 5 ms. Each presynaptic spike produces an exponentially decaying current with time constant τexc = 5 ms for excitation and τinh = 10 ms for inhibition. All neurons are stimulated with a zero-mean white noise current ξ(t), resulting in voltage fluctuation of standard deviation σ = 4 mV. The model equations are thusτmdVm(t)dt=(VrestVm(t))+gexc(t)+ginh(t)+σ(2τm)ξ(t)τsyndgsyn(t)dt=gsyn(t)where syn ∈ {exc, inh} and each presynaptic spike trigger an instantaneous increase: gsyngsyn + λsynwi, where wi is the synaptic weight and λsyn is a scaling factor calculated so that a synaptic weight of 1 mV produces PSPs of peak size 1 mV. Note that gsyn is in units of volt, i.e., the membrane resistance is implicitly included in the variable. For exponential synapses, we can compute analyticallyλsyn=(τmτsyn)τmτmτsyn

Network model. We used a random balanced network (48) composed of 4000 excitatory neurons and 1000 inhibitory neurons. Every neuron in the network is connected to 10% of the others, with delays drawn from a uniform distribution between 0 and 5 ms. To set the network in a balanced regime, we neglected the contribution of spikes and used Campbell’s theorems, which give the mean and variance of a shot noise (corresponding to sums of PSPs with Poisson statistics)<Vm>=FNexcwexc<EPSP>+FNinhwinh<IPSP>+Vrestwhere F is in the firing rate of the network, and Nsyn is the average number of incoming synapses per neuron. In the balanced regime, we want to have <Vm > = Vrest, so we should havewexcNexc<EPSP>=winhNinh<IPSP>

Since the connection probability is uniform, we have Nexc = 4Ninh. Using the fact that <∫EPSP> = λexcτexc and <∫IPSP> = λinhτinh, we have the following equation for the balance g = winh/wexcg=4τexcλexcτinhλinh=2λexcλinh

Initial synaptic weights for recurrent excitatory connections are drawn from a log-normal distribution such that <wexc > = 1 mV, to match experimental findings (see fig. S5D) (49). To be in a strong balanced regime, recurrent inhibitory connections are drawn from a Gaussian distribution such that <winh > = − 1.5g < wexc>. SD of the Gaussian distributions is chosen to be a fourth of the means.

In the spontaneous regime, the network displays an asynchronous irregular regime characterized by a low firing rate (0.2 Hz) (see fig. S5C) and an average coefficient of variation for the interspike intervals close to 1 (see fig. S5G).

Protocol. To mimic the experimental protocol, we randomly selected 100 excitatory neurons within the network and forced them to spikes at controlled times. In half of the trials (100 repeats), the spike had a “normal” effect (ai-AP), while in the other half of the trials (100 repeat), to mimic the fact that this particular spike was considered to be a si-AP spike, all the effective weights of the spiking neuron were artificially increased by 30%. Injections of consecutive spikes were spaced by a short duration of 200 ms, and the si-AP spike effect was limited to the effective weights onto excitatory neurons. Simulations where all effective weights were increased (both to excitatory and inhibitory neurons) showed that the effect was strongly reduced.

Simulator. All simulations were performed using the Brian simulator version 2 (50) with a fixed time step of 0.1 ms.

Quantification and statistical analysis

Membrane potentials are specified with correction for liquid junction potentials. Values indicate mean value ± SEM. Significance of difference between two groups of paired data was assessed by nonparametric Wilcoxon signed-rank test. For comparison between unpaired groups of data, Mann-Whitney U test was used. Differences were considered significant if P < 0.05.


Supplementary material for this article is available at

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.


Acknowledgments: We thank A. Bialowas and V. Marra for sharing preliminary data; U. Gerber, G. Laurent, E. Neher, and W. Singer for critically reading a preliminary version of the manuscript; and R. Cossart and M. Scanziani for helpful discussion. Funding: This work was supported by INSERM, CNRS (to D.D.), Ecole Normale Supérieure (doctoral grant to M.Z.), Agence Nationale de la Recherche (grant ANR-14-CE13-0003 to D.D., R.B., and P.Y.), and Fondation pour la Recherche Médicale (DVS-2013-1228768 to D.D. and FDT-2015-0532147 to M.Z.). Author contributions: M.Z., S.R., and D.D. conceived the project. M.Z., S.R., Y.I., N.B.-G., L.F.-M., and M.R. performed the experiments. M.Z., S.R., Y.I., and D.D. analyzed the data. M.Z., P.Y., and R.B. performed modeling and M.Z., S.R., P.Y., and D.D. wrote the manuscript. D.D. supervised the project. Competing interests: The authors declare that they have no competing interests. Data and material availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.
View Abstract

Stay Connected to Science Advances

Navigate This Article