Research Articles, Development/Plasticity/Repair

Predictable Fluctuations in Excitatory Synaptic Strength Due to Natural Variation in Presynaptic Firing Rate

Naixin Ren, Ganchao Wei, Abed Ghanbari and Ian H. Stevenson
Journal of Neuroscience 16 November 2022, 42 (46) 8608-8620; https://doi.org/10.1523/JNEUROSCI.0808-22.2022

Abstract

Many controlled in vitro studies have demonstrated how postsynaptic responses to presynaptic spikes are not constant but depend on short-term synaptic plasticity (STP) and the detailed timing of presynaptic spikes. However, the effects of short-term plasticity (depression and facilitation) are not limited to short, subsecond timescales. The effects of STP appear on long timescales as changes in presynaptic firing rates lead to changes in steady-state synaptic transmission. Here, we examine the relationship between natural variations in the presynaptic firing rates and spike transmission in vivo. Using large-scale spike recordings in awake male and female mice from the Allen Institute Neuropixels dataset, we first detect putative excitatory synaptic connections based on cross-correlations between the spike trains of millions of pairs of neurons. For the subset of pairs where a transient, excitatory effect was detected, we use a model-based approach to track fluctuations in synaptic efficacy and find that efficacy varies substantially on slow (∼1 min) timescales over the course of these recordings. For many connections, the efficacy fluctuations are correlated with fluctuations in the presynaptic firing rate. To understand the potential mechanisms underlying this relationship, we then model the detailed probability of postsynaptic spiking on a millisecond timescale, including both slow changes in postsynaptic excitability and monosynaptic inputs with short-term plasticity. The detailed model reproduces the slow efficacy fluctuations observed with many putative excitatory connections, suggesting that these fluctuations can be both directly predicted based on the time-varying presynaptic firing rate and, at least partly, explained by the cumulative effects of STP.

SIGNIFICANCE STATEMENT The firing rates of individual neurons naturally vary because of stimuli, movement, and brain state. Models of synaptic transmission predict that these variations in firing rates should be accompanied by slow fluctuations in synaptic strength because of short-term depression and facilitation. Here, we characterize the magnitude and predictability of fluctuations in synaptic strength in vivo using large-scale spike recordings. For putative excitatory connections from a wide range of brain areas, we find that typical synaptic efficacy varies as much as ∼70%, and in many cases the fluctuations are well described by models of short-term synaptic plasticity. These results highlight the dynamic nature of in vivo synaptic transmission and the interplay between synaptic strength and firing rates in awake animals.

Introduction

The firing activity of individual neurons can vary over time because of external stimuli (Hubel and Wiesel, 1959), the behavior of the animal (Churchland et al., 2012; Shenoy et al., 2013), and the internal state of the network (Poulet and Petersen, 2008; McGinley et al., 2015). Although many statistical models assume that coupling between neurons has a fixed strength (Harris et al., 2003; Truccolo et al., 2005; Pillow et al., 2008), firing rate variation can directly affect synaptic strength because of short-term synaptic plasticity (STP). For instance, synapses with short-term depression should be weaker when presynaptic firing rates are higher and stronger when the rates are lower. On the other hand, synapses with short-term facilitation should be stronger when presynaptic firing rates are higher and weaker when the rates are lower. In many studies, STP is induced in vitro by modifying the presynaptic firing rate using external electrical stimulation at different frequencies (Abbott et al., 1997; Tsodyks and Markram, 1997; Fortune and Rose, 2001; Chung et al., 2002; Zucker and Regehr, 2002). These studies illustrate how STP can act to high pass or low pass filter the inputs of neurons and could potentially underlie working memory and decision-making processes (Deng and Klyachko, 2011).

One prediction from these in vitro studies of STP using simple, controlled stimuli is that synaptic strengths should also vary in response to the natural changes in firing rate that occur in vivo. Indeed, slice studies have shown that for neurons receiving naturalistic spike trains as input, synaptic strength varies substantially over time, with patterns matching those predicted by STP (Klyachko and Stevens, 2006; Kandaswamy et al., 2010; Costa et al., 2013). Studies in awake animals have also shown that postsynaptic responses to presynaptic spikes vary depending on the preceding presynaptic interspike intervals (ISIs; Swadlow and Gusev, 2001; Stoelzel et al., 2008, 2009; English et al., 2017) and presynaptic firing rates (Fujisawa et al., 2008; Stoelzel et al., 2015; McKenzie et al., 2021). Together, these findings suggest a general conclusion that the natural variation in presynaptic firing rates in vivo may be accompanied by large fluctuations in synaptic strength and that these changes may be because of the effects of STP on longer timescales.

Quantifying synaptic strength in vivo is a difficult statistical problem. Many studies of synaptic transmission have used intracellular recordings to study the dynamics in individual postsynaptic potentials or currents (Pala and Petersen, 2015, 2018; Sedigh-Sarvestani et al., 2017). Long-term extracellular spike recordings may allow us to track slow fluctuations more easily in vivo. Unlike intracellular signals that contain postsynaptic responses to single presynaptic spikes, extracellular spikes are sparse binary events, and studies of synapses often rely on the cross-correlograms (CCGs) between the presynaptic and postsynaptic spike trains. If two neurons are connected with an excitatory synapse, there is often a fast-onset, short-latency peak in the CCG (Perkel et al., 1967; Fetz et al., 1991), and several techniques have been developed to automatically detect these connected pairs (Barthó et al., 2004; Amarasingham et al., 2012; Kobayashi et al., 2019; Ren et al., 2020).

Synaptic transmission with spikes is often measured by synaptic efficacy using the CCG, which is defined as the excess probability of postsynaptic spiking in the short transmission window following presynaptic spikes (Levick et al., 1972). By focusing on specific presynaptic spike patterns, such as pairs with specific ISIs, previous studies have shown how synaptic efficacy can vary depending on features of presynaptic spike timing (Usrey et al., 2000; Swadlow and Gusev, 2001; English et al., 2017). Modeling individual spike transmission probabilities for the entire spike train can also reveal how presynaptic spike features, as well as different types of plasticity, interact (Ghanbari et al., 2017, 2020; Song et al., 2018; Wei and Stevenson, 2021). Together, descriptive and model-based approaches can track synaptic dynamics in vivo (Swadlow and Gusev, 2001; English et al., 2017; Ghanbari et al., 2017, 2020) and reveal changes associated with learning, behavior, or stimuli (Fujisawa et al., 2008; Ghanbari et al., 2020; McKenzie et al., 2021). Here, we focus on a general description of the ongoing changes in synaptic strength that may be merely because of changes in presynaptic rate. We divide long-term spike recordings into windows to measure the fluctuations in synaptic efficacy on a timescale of minutes and model the changes as a function of presynaptic ISI to build the potential link between these slow changes and STP.

Using in vivo large-scale spike recordings from the Allen Institute Neuropixels dataset (Siegle et al., 2021), we examine natural fluctuations in synaptic efficacy and firing rates and quantify the relationships between them. These recordings from awake, behaving mice contain the spiking of hundreds of neurons from many brain regions while the mice view natural and artificial visual stimuli (Siegle et al., 2021). Because neurons in different brain regions have different firing rates and firing regularity (Mochizuki et al., 2016), we also examine whether putative synapses from different regions show different patterns of efficacy fluctuations. We find that for individual connections the synaptic efficacy often varies substantially over the course of the recording. Interestingly, the efficacy is often correlated with the presynaptic firing rate, with the postsynaptic firing rate, or with both. We then test whether the time-varying efficacy is predictable based on the presynaptic firing. We model the time-varying synaptic effect as a function of presynaptic ISI by fitting a generalized bilinear model (GBLM) to the postsynaptic spike trains. For many connections, the GBLM reproduces the observed efficacy fluctuations, suggesting that these fluctuations can be, at least partly, explained by a cumulative effect of STP. Altogether, these results illustrate how natural fluctuations in firing rate could lead to large, but predictable, changes in synaptic strength in vivo and in many different brain areas.

Materials and Methods

We use the Visual Coding–Neuropixels dataset from the Allen Institute (https://portal.brain-map.org/explore/circuits). During the recordings, head-fixed mice viewed a standardized set of visual stimuli (including Gabor patches, full-field drifting gratings, moving dots, and natural images and movies) while they were free to run on a wheel. We analyze 20 recordings here (electrophysiology sessions 715093703–746083955 and 819186360–847657808), 6 from female mice (125 ± 11.6 d), 14 from male mice (117.8 ± 8.3 d). Nine recordings were from wild-type mice, and the remainder were transgenic Cre lines expressing channelrhodopsin in inhibitory interneuron subtypes (Sst-IRES-Cre; Ai32, Pvalb-IRES-Cre; Ai32 and Vip-IRES-Cre; Ai32). Each recording lasts ∼2.5 h and contains 365–893 well-isolated single units from five or six Neuropixels probes from mouse visual cortex, thalamus, and hippocampus. In 10 recordings, the mice had a resting period that lasted for ∼30 min. This period is referred to as spontaneous activity period for further analysis. Additional information about this dataset can be found in Siegle et al. (2021).

Detection of putative excitatory synaptic connections

We detect putative synaptic connections using the CCGs between pairs of neurons—a source neuron (presynaptic) i and target neuron (postsynaptic) j. For the binned spike trains of the two neurons, ni and nj (1 when there is a spike and 0 otherwise), the CCG is given by yij(m)=∑tni(t)nj(t−m), where m denotes the interval between presynaptic and postsynaptic spikes, and yij(m) is the number of the times spikes in ni and nj that are separated by an interval [m−Δtbin2,m+Δtbin2] for bin size Δtbin (1 ms used here).

Given the CCG, we detect putative excitatory synapses in two stages. First, we use a hypothesis test based on the jittered spikes to detect fast changes in the CCG, and then we use a model-based approach to estimate and screen for putative synaptic effects at specific latencies and timescales. Our initial screening of neuron pairs is based on comparing the observed CCG to a null model of the CCG generated using jittered spike times. An exact null distribution can be approximated by directly sampling jittered spike trains. However, here, for computational convenience, we approximate the null distribution using the simplifying assumption that the counts in each bin of the CCG yij(m) are independent of each other (Amarasingham et al., 2012). We assume that yij(m)∼Binomial(numpre,μ(m)/numpre), where numpre is the total number of presynaptic spikes and μ(m)=(yij * f)(m), the convolution of the CCG with a jitter function. We use a 10 ms uniform jitter f(t)=1 for |t|<5ms, 0 otherwise. We then compute p values for each bin of the CCG using the cumulative distribution function of the binomial distribution and evaluate statistical significance correcting for multiple comparisons by controlling the false discovery rate using a Benjamini–Hochberg procedure (α=0.00001). We then identify pairs of neurons where two criteria are satisfied, (1) there are two neighboring statistically significant bins in the CCG (where m > 0) and (2) the center bin (m=0) is not statistically significant. This test yields 0.5% of pairs of neurons that have fast, transient changes in their CCGs.

For each pair of neurons that satisfies the testing criteria we then fit a parametric model that aims to estimate the latency, timescale, and strength of a putative synaptic effect. We fit an extended GLM based on the model used in Ren et al. (2020). Briefly, we describe the CCG using two components, (1) a slow fluctuation because of background fluctuations and (2) a fast, transient effect because of the synaptic effect. The rate of spike counts λ at bin m in a CCG is given by the following: λ(m)=exp(sXs(m) + wα(m,τ,Δt)).

Here, sXs describes the slow fluctuation, where Xs represents a set of smooth basis functions learned by applying a low-rank, nonlinear matrix factorization to all the cross-correlograms in the dataset (Ren et al., 2020); and wα(m,τ,Δt) represents the fast synaptic effect, where w is the synaptic strength and α(m,τ,Δt) is an alpha function with a time constant τ and a latency Δt that is then convolved with the autocorrelogram of the presynaptic neuron to account for possible bursting in the neuron. After fitting the parameters {s,w,τ,Δt} using a penalized maximum likelihood method, we select the neuron pairs where the CCG has a strong sharp peak (w>0.3,τ<0.8ms) with a short latency (Δt<10 ms) and a relatively flat slow fluctuation (the coefficient of variation of the slow fluctuation < 0.15). In the end, we detect 1383 putative excitatory connections (∼0.17% of all pairs of neurons). This screening procedure returns a set of neuron pairs whose CCGs are consistent with the biophysics of excitatory synaptic transmission, which we then characterize in more detail.

Estimating time-varying synaptic efficacy

Synaptic efficacy reflects the (excess) probability of observing a postsynaptic spike following a presynaptic spike. Here, we estimate synaptic efficacy using the model-based approach described in the previous section. Namely, we calculate the efficacy by comparing the model with and without the synaptic effect as follows: eff=∑m(λ(m)−λslow(m))/numpre, where λslow is λ evaluated with θ={s,w=0} and numpre is the number of presynaptic spikes. To estimate changes in observed efficacy over time, we divide the presynaptic and postsynaptic spike trains into 5 min windows with 80% overlap and estimate efficacy for each window. We fit the cross-correlogram for each window using the extended GLM with a fixed alpha function α(m,τ̂,Δt̂), where τ̂ and Δt̂ are the estimated parameters from the overall cross-correlogram of the connection. We calculate the coefficient of variation (CV; the ratio of the SD to the mean) for time-varying efficacy to measure the fluctuation.

For comparison, we also shuffle the postsynaptic spike timings to create surrogate data and then calculate the fluctuations in efficacy that would be expected by chance. With this shuffling, the overall average efficacy is preserved, but any fluctuations in the efficacy over long time periods should be because of chance. To do this shuffling, we find all pairs of presynaptic and postsynaptic spikes with an interval of (0, 25] ms. Then for each presynaptic spike in this subset of spike pairs, we resample the subsequent postsynaptic spikes. Any postsynaptic spikes that do not occur following a presynaptic spike within the interval are not changed. This way the postsynaptic spikes following each presynaptic spike in the CCG are shuffled across time, and any systematic, slow changes in the CCG are removed. For each putative synapse, we do the shuffling 100 times to generate 100 surrogate time-varying efficacies. We can then create a surrogate distribution of the CV for efficacy and calculate a z score for the observed CV to measure how far away it is from the surrogate distribution. In addition, to examine the relationship between the time-varying efficacy and presynaptic firing rate, we calculate Spearman's correlation coefficient, and then similarly calculate a z score for the observed correlation using the surrogate data.

Estimating synaptic efficacy with different presynaptic ISI ranges

To further examine the relationship between efficacy and the presynaptic firing, we also estimate the efficacy for presynaptic spikes with different ISIs. Here, we group the presynaptic spikes by the ISI range their preceding spike fell within and calculate a separate CCG for each subset of presynaptic spikes. We then estimate the efficacy using the extended GLM for each separate ISI range.

Modeling short-term synaptic plasticity with spike train observations

In addition to quantifying the observed efficacy from noisy correlograms, we also aim to assess to what extent changes in efficacy can be explained by STP. Here, we model the time-varying synaptic effect as a function of presynaptic ISI by fitting a GBLM to the postsynaptic spike trains. This model has been previously described (Ghanbari et al., 2017; Wei and Stevenson, 2021) and extends previous models of coupled GLMs (Harris et al., 2003; Truccolo et al., 2005; Pillow et al., 2008; Rebesco et al., 2010) to account for fluctuating excitability and plasticity. Briefly, for each connection, we model the postsynaptic firing rate γ at time t as follows: γ(t|β0(t),βs,βc,w(t))=exp(β0(t) + βsyslow(t) + βcycoup(t)w(t)) npost(t)∼Poisson(γ(t|β0(t),βs,βc,w(t))) yslow=npre * xslow ycoup=npre * xcoup, where npre and npost are the presynaptic and postsynaptic spike trains, respectively.

Here, β0(t) represents a time-varying baseline firing rate, which is estimated using an adaptive filtering algorithm (Wei and Stevenson, 2021). βsyslow(t) accounts for slow fluctuations that are shared by the presynaptic neuron and postsynaptic neuron, where yslow denotes a set of covariates generated by filtering the presynaptic spike train (cubic B spline functions with four equally spaced knots over 150 ms xslow). βcycoup(t)w(t) represents the contribution from the synaptic effect. Here, βcycoup(t) describes the stationary coupling effect, where ycoup is the presynaptic spike train convolved with the alpha function we learned from the extended GLM α(m,τ,Δt), and w(t) describes the effects of STP. Namely, we model w(t) as a function of presynaptic ISIs Δtk as follows: w(t)=1+∑kqk(Δtk)∑l=0Tδ(tk−l)exp(−lτstp), where qk(Δtk) is a modification function that determines how the synaptic strength increases/decreases following each presynaptic spike at time tk. We model the modifications using a group of raised cosine basis functions B(Δtk) and qk(Δtk)=θB(Δtk) and assume that the effects of STP decay with τstp=200 ms. We fit the free parameters {β0(t),βs,βc,θ} by maximizing the Poisson log-likelihood as follows: LLHGBLM=∑t(npost(t)logγ(t)−γ(t)).

To estimate the parameters, we use an alternating optimization (coordinate ascent); we fit β0(t) using fixed {βs,βc,θ} and then fit {βs,βc,θ} using fixed β0(t). We repeat this alternating pattern until convergence.

After optimization we then have a model-based estimate of the impact of each individual presynaptic spike and can predict how the efficacy would change over time under a model of short-term plasticity. We calculate GBLM estimated synaptic efficacy in a similar way as we estimate efficacy using CCGs as follows: effGBLM=∑t(γ(t)−γ0(t))/numpre, where γ0 is the postsynaptic rate with the synaptic effect removed βc=0; effGBLM then reflects the excess probability of observing a postsynaptic spike because of the occurrence of a presynaptic spike within the model.

For comparison, we also fit a static-synapse model where the synaptic coupling effect is βcycoup(t) instead of βcycoup(t)w(t). Under the static-synapse model, the coupling effect stays the same over time and will allow us to assess whether STP, specifically, meaningfully improves the description of spike transmission. To quantify how well a model predicts whether there is a postsynaptic spike after each presynaptic spike, we use receiver operating characteristic (ROC) curves, and the area under the curve (AUC) reflects the performance. For each presynaptic spike, we identify a transmission interval (ta,tb) where the coupling filter (alpha function in GBLM) is >0.5, and then we use the model fit within the transmission interval 1−∏t∈(ta,tb)pp(npost=0|γ(t)) as the score for predicting each individual postsynaptic spike.

Data availability

The MATLAB code for data analysis is available from GitHub at https://github.com/NaixinRen/predictable-fluctuations-in-synaptic-efficacy.

Results

Here, we characterize fluctuations in putative synapses from 20 recordings from the Allen Institute Neuropixels dataset (Siegle et al., 2021). We detect putative excitatory synaptic connections based on the CCGs of pairs of neurons using a two-stage synapse detection method. In stage 1, we use a hypothesis test based on the jittered spikes to detect fast changes in the CCG, and in stage 2, we use a model-based approach to estimate and screen for putative synaptic effects at specific latencies and timescales (see above, Materials and Methods). Figure 1A shows the CCGs of example connections. By using this two-stage screening procedure, our goal is not to detect as many putative synapses as possible but rather to generate a subset of strong putative synapses with CCGs that are consistent with the biophysics of excitatory synaptic transmission, including synaptic latencies (1.29 ± 1.23 SD ms; Fig. 1B) and time constants (0.55 ± 0.13 SD ms; Fig. 1C).

Figure 1.
Figure 1.

Detection of putative excitatory synaptic connections. A, Examples of CCGs of putative excitatory synaptic connections. The bin size is 0.5 ms here for better visualization. In the analysis, we use 1 ms. B, The histogram of synaptic latency of all the detected synapses. The distribution tends to peak at 1, 2, 3 … [ms] because the bin size of the CCG is 1 ms. C, The histogram of the synaptic time constant of all the detected synapses. D, Two measurements of the duration of waveforms. Each dot represents one neuron. Most putative presynaptic neurons (yellow dots) show broader spike waveforms compared with all the neurons in the dataset (gray dots). E, The histogram of the synaptic efficacy of all the detected synapses. Left, Example CCGs with low or high efficacy and the model fits of the extended GLM (red line, model fits with synaptic effect; yellow line, model fits without synaptic effect). F, Connectivity matrix across different areas (y-axis, presynaptic area; x-axis, postsynaptic area). Color represents the number of connections detected (log transformed). The gridlines separate the brain regions of visual cortex, thalamus, hippocampus, midbrain, and unclassified. Visual cortex (VIS), Primary Visual Area (VISp), Lateral visual area (VISl), Rostrolateral visual area (VISrl), Anteromedial visual area (VISal), Posteromedial visual area (VISpm), Anteromedial visual area (VISam), Mediomedial posterior visual area (VISmmp), Mediomedial anterior visual area (VISmma), thalamus (TH), dorsal/ventral lateral/medial geniculate (LGd/LGv/MGd/MGv), Lateral posterior nucleus (LP), Intergeniculate leaflet of the lateral geniculate complex (IGL), Posterior nucleus (PO), Posterior limiting nucleus (POL), suprageniculate nucleus (SGN), Ventral posteromedial nucleus (VPM), Dentate gyrus (DG), Subiculum (SUB), Prosubiculum (ProS), Anterior pretectal nucleus (APN), Nucleus of the optic tract (NOT), Posterior pretectla nucleus (PPT), Superior colliculus, intermediate grey matter (SCig), lateral ventricle (VL). G, Box plots for the mean firing rate (top) and efficacy (bottom) for brain areas. The areas are grouped by brain regions mentioned above. Only areas with ≥10 efferent synapses are shown here. The line inside each box is the median. The top and bottom edges of the box are the upper and lower quartiles, respectively. The whiskers show the nonoutlier maximum and minimum. The outliers are not shown here.

The presynaptic neurons for the majority of putative synapses that we detect have broad waveforms (Fig. 1D), which previous work has found tend to correspond to excitatory neurons (Connors and Gutnick, 1990; Barthó et al., 2004; Trainito et al., 2019). For each putative synapse, we estimate the overall synaptic efficacy using an extended GLM (see above, Materials and Methods; Fig. 1E). The distribution of synaptic efficacy is approximately log normal (0.037 ± 0.035, mean ± SD; Fig. 1E), consistent in shape with the previous findings on postsynaptic potentials (Buzsáki and Mizuseki, 2014) and efficacy (English et al., 2017). Here, we detect 1383 putative excitatory synapses in total, most of which were between neurons recorded on the same probe (1174 of 1383) and within the same brain area (1090 of 1383; Fig. 1F).

Across the 20 recordings, we identify putative synapses from many labeled brain regions. The regions with >10 efferent synapses included six visual cortical areas (VISp, VISl, VISrl, VISal, VISpm, VISam), three thalamic areas (LP, MGd, MGv), three hippocampal areas (CA1, CA3, DG), and in locations that were not precisely identified (grey or VIS). Synapses from different brain regions have distinct presynaptic firing rates (5.75 ± 9.93 SD spikes/s overall average) and synaptic efficacies (Fig. 1G). Notably, putative synapses from thalamic areas appear to have higher presynaptic rates and lower efficacy on average. Firing rates are expected to fluctuate because of varying stimuli, behavior, and internal state variables. Here, we aim to measure the extent to which the efficacy of putative synapses also fluctuates and to determine whether these fluctuations are predictable.

Presynaptic firing rate and efficacy both fluctuate substantially over time

To examine the natural fluctuations in presynaptic firing rate and synaptic efficacy, we divide the presynaptic and postsynaptic spike trains into 5 min windows with 80% overlap and estimate the presynaptic firing rate and efficacy for each window separately (Fig. 2A,B, example of fluctuations). Here, we calculate the CV (the ratio of the SD to the mean) for both variables for each connection to measure the fluctuation. We find that both presynaptic firing rate and synaptic efficacy fluctuate substantially over time in all areas (Fig. 2C,D; CV for firing rate, 0.55 ± 0.29 SD; CV for efficacy, 0.69 ± 0.46 SD). Although large fluctuations in firing rate are expected, large fluctuations in synaptic efficacy of this scale are somewhat unexpected during ongoing behavior. When firing rates are low, efficacy estimates are noisy, and some variation in efficacy is expected purely by chance. However, the observed fluctuations are not simply explainable by measurement error. Here, we compare the observed fluctuations in efficacy to surrogate data where postsynaptic spike trains are shuffled to remove any systematic fluctuations over long timescales but preserve the overall efficacy (see above, Materials and Methods, Estimating time-varying synaptic efficacy). For each putative synapse, we shuffle 100 times to create a null distribution for CV of efficacy and calculate a z score for the observed CV using the mean and SD of the null distribution. The average z score across all putative synapses is 3.26, suggesting that the observed efficacy fluctuations are substantially higher in most cases than would be expected by chance. In addition, we find that the CV for firing rate is correlated with CV for efficacy for individual synapses (r=0.24,p<0.001). Averaging within brain regions, we also find that the CV of firing rate and the CV of efficacy are positively correlated across regions (r=0.77,p<0.01; n = 14 regions with ≥10 putative synapses). Regions with greater fluctuations in the presynaptic firing rates tend to have higher fluctuations in the efficacies of their putative synapses.

Figure 2.
Figure 2.

Fluctuations in presynaptic firing rate and synaptic efficacy. A, B, Three examples for fluctuations in presynaptic firing rate (A) and synaptic efficacy (B), respectively. Top, Low CV. Middle, Moderate CV. Bottom, High CV. The firing rates in A are estimated using optimal time bins that are generated based on Shimazaki and Shinomoto (2010) for better visualization on finer timescales. The efficacy in B is estimated based on the 5 min window described in Materials and Methods. C, D, The box plots for CVs of presynaptic firing rate (C) and synaptic efficacy (D) for each area. Only areas with ≥10 efferent synapses are shown here. The notations are the same as in Figure 1.

Efficacy fluctuations mirror fluctuations in the presynaptic firing rate for individual neuron pairs

The efficacies of putative synapses appear to fluctuate substantially in vivo across multiple brain regions. To examine these fluctuations in more detail, we calculate correlations between the time-varying presynaptic firing rate and time-varying synaptic efficacy for individual connections. The average absolute Spearman's correlation is 0.32 ± 0.22 SD (983 of 1383 connections are statistically significant; p < 0.05; Fig. 3A,B). As with the CV of efficacy fluctuations, these correlations do not appear to be simply because of chance. By shuffling the postsynaptic spike times as described above we can also measure to what extent the observed Spearman's correlation for each synapse differs from that expected by chance (see above, Materials and Methods, Estimating time-varying synaptic efficacy). Here, the average z score of the synapses with positive correlation coefficient is 2.34, and the average z score of the synapses with negative correlation coefficient is −1.83 (Fig. 4A–C, third column). This suggests that efficacy fluctuations may be predicted based on the presynaptic firing rate, with positive or negative correlations larger than expected by chance, depending on the neuron pair. For comparison, we also examine the correlation between the fluctuations in the postsynaptic firing rate and synaptic efficacy; the average absolute Spearman's correlation in this case is 0.38 ± 0.23 SD (1095 of 1383 connections are significant; p < 0.05; Fig. 3C,D). Positive correlations between the fluctuations in postsynaptic firing rate and synaptic efficacy are somewhat expected as the excitability of the postsynaptic neuron may affect both and increased synaptic efficacy might causally affect the firing rate of the postsynaptic neuron (>75% of the synapses show a positive correlation in Fig. 3C). However, the correlation between the fluctuations in the presynaptic firing rate and efficacy is less predictable and worth exploration. To further determine the influence of the presynaptic firing rate on synaptic efficacy, for each putative synapse, we run a linear regression predicting the time-varying efficacy from the log-transformed time-varying postsynaptic firing rate and the log-transformed time-varying presynaptic firing rate. The Wald test shows that for 67.1% of the putative synapses, the presynaptic firing rate makes a statistically significant contribution to the prediction (928 of 1383, p < 0.05). For many synapses, the fluctuations in the synaptic efficacy mirror the natural fluctuations in presynaptic firing rate.

Figure 3.
Figure 3.

Correlation between the fluctuation in synaptic efficacy and that in the presynaptic and postsynaptic firing rate. A, C, Histogram of the Spearman's correlation between the fluctuations in synaptic efficacy and the natural fluctuations in the presynaptic firing rate (A, blue) or the postsynaptic firing rate (C, green) and the corresponding distributions generated using postsynaptic spikes jittered on a uniform interval of [−1,1] s (gray). B, D, Violin plot of the Spearman's correlation for presynaptic firing rate (B) or postsynaptic firing rate (D) for each area. Only areas with ≥10 efferent synapses are shown.

Figure 4.
Figure 4.

Three kinds of relationships between the fluctuations in synaptic efficacy and presynaptic firing rate. A, An example putative synapse (DG–CA3), where there is a positive correlation between the synaptic efficacy and the presynaptic firing rate. B, An example putative synapse (APN–APN), where there is a negative correlation between the synaptic efficacy and the presynaptic firing rate. C, An example putative synapse (CA1–VIAam), where the relationship between the synaptic efficacy and the presynaptic firing rate is not clear. First column (A–C), CCG (top) and autocorrelogram (ACG, bottom). Second column (A–C), The fluctuations in the presynaptic firing rate (top) and the synaptic efficacy (bottom). The shaded area represents 1 SD. Third column (A–C), The distribution of Spearman's correlation from surrogate, shuffled data (gray, 1000 samples shown for visualization). The red bar shows the observed correlation from the original data. Note that in C the observed correlation is close to zero. This does not mean the fluctuation in synaptic efficacy is random, it just indicates that the fluctuation in synaptic efficacy is not correlated with presynaptic firing rate. Fourth column, The distribution of the presynaptic ISI (top) and the efficacy estimated for presynaptic spikes with different ISIs using the extended GLM.

One concern with this analysis may be that some of the putative synapses are false positives. To address this concern, we consider jittered data where the fast changes in the CCGs are removed. We calculate efficacy using the same neuron pairs where the postsynaptic spikes are jittered on a uniform interval of [−1, 1] s. We find that the Spearman's correlation between firing rate and efficacy in this case is substantially smaller compared with the original spikes (Fig. 3A,C, gray).

One potential explanation for the correlation between presynaptic rate and efficacy is short-term STP. If a synapse exhibits short-term synaptic depression, higher presynaptic firing rates lead to a depletion of resources and lower efficacy, leading to a negative correlation. On the other hand, if a synapse exhibits short-term synaptic facilitation, higher presynaptic rates lead to an increase in release probability and higher efficacy, leading to a positive correlation. In other cases, facilitation and depression may coexist, leading to nonmonotonic relationships between presynaptic rate and efficacy. Here, we find a range of relationships (Fig. 4A,B) as well as cases where there are clear fluctuations but no clear relationship between presynaptic rate and efficacy (Fig. 4C). Previous studies have found short-term changes in synaptic efficacy in vivo by showing how synaptic efficacy varies depending on the preceding presynaptic ISIs (Swadlow and Gusev, 2001; Stoelzel et al., 2008, 2009; English et al., 2017; Ghanbari et al., 2020). As previous studies have done, here, we estimate efficacies at different presynaptic ISI ranges (see above, Materials and Methods, Estimating synaptic efficacy with different presynaptic ISI ranges). We find that the relationships between efficacy and presynaptic ISI (Fig. 4A–C, fourth column) are often consistent with the correlation between efficacy and presynaptic firing rates we find on a minute timescale. For example, the putative synapse in Figure 4A exhibits short-term facilitation, and there is a consistent efficacy–ISI relationship where short presynaptic ISI leads to high efficacy. Overall, this suggests that fluctuations in synaptic efficacy can potentially be predicted in detail by modeling STP as a function of presynaptic ISI.

A model of short-term plasticity reproduces efficacy fluctuations

Because many synapses show a facilitating or depressing relationship between the presynaptic firing rate and synaptic efficacy, we examine explanatory models based on STP. Here, we build a GBLM that aims to describe the detailed (1 ms bins) postsynaptic spiking by modeling STP, slow fluctuations in postsynaptic excitability, and slow common input to both the presynaptic and postsynaptic neurons (see above, Materials and Methods). Here, we estimate a modification function that describes how the synaptic strength changes following the occurrence of presynaptic spikes at a given ISI, and we use adaptive smoothing to track slower fluctuations in the postsynaptic rate. Altogether, this model captures the probability of postsynaptic spiking following every single presynaptic spike and allows us to track how the time-varying synaptic strength increases (facilitation) or decreases (depression) following presynaptic spikes with specific intervals.

To quantify how well the models predict the postsynaptic spiking after a presynaptic spike, we use ROC curves, where the AUC reflects performance (see above, Materials and Methods). We find that the GBLM predicts postsynaptic spikes better than the chance level (AUC, 0.68 ± 0.05 SD). In the 10 recordings where there is a 30 min period of spontaneous activity, we find that the GBLM predicts postsynaptic spikes better than the chance level even when only examining this spontaneous period where we calculate the AUC using only the spikes during the spontaneous activity (AUC, 0.67 ± 0.06 SD). Although many putative synapses involve neurons in nonvisual regions, this comparison illustrates how fluctuations in efficacy are predictable even in the absence of specific visual stimuli and are not merely an artifact of fast stimulus correlations.

The GBLM provides a model of the time-varying impact of each individual presynaptic spike. However, to illustrate how slow fluctuations in efficacy might be explainable by STP more directly, we evaluate the predictions of the GBLM across 5 min windows as before. We find that the median coefficient of determination (R2) between the GBLM estimated efficacy and the observed efficacy is 0.22 [0.07–0.43 interquartile range (IQR); Fig. 5A; 86.6% statistically significant, p<0.05], and the results do not vary substantially across different brain regions (Fig. 5B). Moreover, the modification function of presynaptic ISI from the GBLM tends to be consistent with the relationship between presynaptic firing rate and observed synaptic efficacy (Fig. 5C–E). As a comparison, we also build a static-synapse model without STP, where the synaptic coupling stays the same over time (see above, Materials and Methods). In this case, the median coefficient of determination between the static model estimated efficacy and the observed efficacy 0.17 (0.05–0.37 IQR; Fig. 5A; 83.3% of them show p<0.05). The static model can reproduce some fluctuations in the synaptic efficacy, indicating that some of the fluctuations can be attributed to non-STP components, such as the excitability of the postsynaptic neuron (Gal et al., 2010) or slow common inputs shared by the presynaptic and postsynaptic neurons. However, there is 30% improvement in the median R2 when comparing the GBLM to the static model, which suggests that STP plays a role in the efficacy fluctuations as well. For many connections, the GBLM can reproduce fluctuations in the synaptic efficacy with high fidelity, suggesting that STP can, at least partially, explain the relationship between fluctuations in the presynaptic firing rate and synaptic efficacy.

Figure 5.
Figure 5.

GBLM of STP reproduces the fluctuation in synaptic efficacy. A, The cumulative probability of R2 between the model estimated efficacy and the observed efficacy. The dotted line shows the median (50% cumulative probability). B, Box plot of R2 for different areas. Only areas with ≥10 efferent synapses are shown. The notations are the same as in Figures 1–4. C–E, The synapses are the same as those shown in Figure 4. Left, The fluctuations in the presynaptic firing rate (top) and the synaptic efficacy (bottom). The blue line is the observed efficacy. The orange line is the GBLM estimated efficacy. The yellow line is the static model estimated efficacy. Right, The distribution of the presynaptic ISI (top) and the modification function estimated from the GBLM. The shaded orange area represents 1 SE.

To characterize the diversity of STP patterns in the putative excitatory connections examined here, we further analyze the modification functions estimated from the GBLM. The modification function gives a simple description of how the synaptic effect varies as a function of the presynaptic ISI. When the modification is positive, synaptic strength is facilitated, and when it is negative it is depressed. However, the influence of the modification function, overall, depends on the ISI distribution. Although the modification function may be negative (for ISI < 10 ms in Fig. 5C, e.g.), this will only influence the synaptic strength when spikes occur with those particular ISIs. Here, we align the modification function by ISI distributions. We sample the modification function at the 5th–95th percentiles of the ISI distribution (with 5% increments) and normalize by dividing the ISI-scaled function by the absolute value of its average. To see if there is any structure in the modification functions for all the putative synapses, we use hierarchical clustering on the ISI-scaled modification function (Euclidean distance with unweighted average linkage). We find the functions are broadly split into two clusters with some putative synapses mostly facilitating and others mostly depressing. However, the specific ISI range where facilitation or depression occurs varies widely across putative synapses (Fig. 6A). Summarizing the ISI-scaled modification functions with their average, we find that the overall level of facilitation or depression appears to vary across brain regions (Fig. 6B) in a similar trend as the relationship between the minute-level fluctuations in presynaptic firing rate and efficacy (Fig. 3B). These results suggest that differences in the association between efficacy and presynaptic rate that we observe on ∼1 min timescales may be related to differences in STP on subsecond timescales.

Figure 6.
Figure 6.

Hierarchical clustering for modification functions. A, Right, The modification functions of all putative synapses sorted based on the clustering results; each row represents one modification function. Left, The dendrogram of the hierarchical clustering. Colors represent the top two clusters. B, The violin plot of the mean of the modification function for each brain area. Only areas with ≥10 efferent synapses are shown. C, Schematic illustrating how the same modification rule and presynaptic ISI distribution can result in different efficacy fluctuations. Top, Two example presynaptic spike trains with the same number of spikes and the same ISI distribution but different ISI sequences. Bottom, The resulting coupling effects for both spike trains simulated using one facilitation modification function. Note that the presynaptic spike train on the left results in a synaptic strength that is stronger on average than the train on the right. D, The GBLM estimated efficacy from the original data (the orange lines in Fig. 6D are the same as the orange lines in Fig. 5C,D (bottom), and from the postsynaptic activity reconstructed using shuffled presynaptic spikes (gray).

Although the sign of the modification functions estimated with the GBLM typically reflects the positive or negative correlation between the minute-level fluctuations in the presynaptic firing rate and synaptic efficacy, it is also important to note that this relationship is not one to one. The modification function summarizes how the putative synaptic weight changes based on the single preceding ISI, whereas the overall efficacy in a time window is a result of not only the preceding ISI but also the sequence of presynaptic ISIs within the window. The same ISI distribution can give rise to different efficacy fluctuations depending on the specific order of the presynaptic spikes (Fig. 6C). To illustrate the influence of specific presynaptic spike patterns on the observed efficacy fluctuations, we compared postsynaptic responses to the original spikes with responses to shuffled presynaptic ISIs with the same modification function. For each example synapse (same as those in Fig. 5C,D), we shuffle the presynaptic ISI within a 1 min window, so the ISI distribution and minute-level fluctuation in the presynaptic firing rate are preserved but the ISI sequence is randomized. We then reconstruct the postsynaptic activities using the GBLM parameters estimated from the original data where the coupling effect changes as a function of the shuffled presynaptic spike train. For both synapses, the predicted efficacy fluctuations with shuffled presynaptic ISIs (Fig. 6D, gray lines) are weaker compared with the predicted fluctuations from the original presynaptic spikes (Fig. 6D, orange lines) even under the same modification functions. These differences demonstrate how, in addition to the modification function and the overall presynaptic rate, the detailed sequence of presynaptic spikes plays a role in determining the correlation between presynaptic firing rate and synaptic efficacy.

Discussion

In this study we examine associations between fluctuations in firing rates and efficacies of putative excitatory synapses using large-scale extracellular spike recordings from awake, behaving mice. Using both descriptive and model-based analysis we identify three main findings. First, we found that for many putative synapses there are substantial fluctuations in the synaptic efficacy that occur on a slow timescale of minutes alongside the expected variations in presynaptic and postsynaptic firing rates. Second, we found that these slow fluctuations in synaptic efficacy are often correlated (either positively or negatively) with fluctuations in presynaptic rate. And third, by fitting a detailed 1 ms resolution model of postsynaptic spiking probability where synaptic effects have short-term plasticity, we show that these slow fluctuations in synaptic efficacy are consistent with the cumulative effects of STP in many cases.

Although previous studies have shown that synaptic efficacy varies in vivo with learning, behavior, and stimuli (Fedulov et al., 2007; Fujisawa et al., 2008; Ghanbari et al., 2020; McKenzie et al., 2021), differences in efficacy across conditions may partially reflect differences in presynaptic firing rates and spike patterns across conditions. Firing rates depend both on external variables, such as stimuli and movement (Ben-Yishai et al., 1995; Kang et al., 2004; Vogels et al., 2005), as well as internal variables, such as arousal, attention, adaption, and excitability (Gilbert and Sigman, 2007; Gal et al., 2010; Goris et al., 2014). Firing rate fluctuations can occur spontaneously (Teich et al., 1990; Yaksi and Friedrich, 2006; Faisal et al., 2008) and are influenced by brain state and the population activity of the surrounding neurons (Tsodyks et al., 1999; Poulet and Petersen, 2008; McGinley et al., 2015). These many sources of firing rate variability occur simultaneously in vivo and may have complex effects on synaptic efficacy because of diversity in both cell types and synaptic dynamics. Understanding how synaptic efficacy is affected by presynaptic firing rates first may allow us to better understand the influence of specific external or internal variables.

Here, we find that the brain regions with larger fluctuations in the presynaptic firing rate (measured by CV) also show larger fluctuations in efficacy of putative excitatory synapses. It is important to note that this observation is just based on the putative excitatory synapses we detected. Synapse detection from spikes is biased toward strong synapses between neurons with high firing rates. Moreover, most synapses analyzed here are between neurons within the same brain area and detected on the same probe. These findings thus represent only a partial picture of the overall synaptic dynamics in these regions. However, the correlations we observe are consistent with our core hypothesis that firing rate variation leads directly to fluctuations in synaptic efficacy because of STP. Here, we find such patterns not just across brain areas, but for many individual synapses, fluctuations in synaptic efficacy mirror fluctuations in the presynaptic firing rate. Synaptic efficacy can be affected by the postsynaptic activity as well (Abbott and Nelson, 2000; Ibata et al., 2008; Sjöström et al., 2008), but here we find that the presynaptic firing rate is still predictive, even after accounting for the postsynaptic firing rate. Previous in vitro experiments have shown that natural spike patterns lead to changes in synaptic strength consistent with STP on short timescales (Klyachko and Stevens, 2006; Kandaswamy et al., 2010). Here, we find evidence that synaptic efficacy fluctuates with presynaptic firing rate in vivo, during behavior, with large fluctuations on timescales of minutes.

Using a model-based approach to predict postsynaptic spiking on 1 ms timescales (GBLM), we also examined to what extent the slow fluctuations in synaptic efficacy can be explained by STP. Here, we modeled the time-varying synaptic effect using a modification function where the synaptic effect changes as a function of presynaptic ISI. We find the GBLM can predict when individual postsynaptic spikes will occur after a presynaptic spike better than chance level, as well as accurately describe the slow fluctuations in efficacy on timescales of minutes. The modification functions estimated by the model show a range of ISI–efficacy relationships, both depressing and facilitating, and the distribution of this relationship varies in different brain areas. Previous intracellular studies have shown that short-term synaptic dynamics often depend on the presynaptic and postsynaptic neuron type and brain area (Thomson and Lamy, 2007; Blackman et al., 2013; Lee et al., 2019; Campagnola et al., 2022), and the presynaptic spike statistics that drive STP also vary with cell type and brain area (Mochizuki et al., 2016). A more detailed comparison of in vivo results with previous in vitro results (Seeman et al., 2018) may help determine to what extent in vitro results generalize to naturalistic settings. Similarly, comparisons to previous in vivo studies, such as those finding depression in thalamocortical synapses, (Chung et al., 2002; Swadlow and Gusev, 2002; Bruno and Sakmann, 2006), can help determine to what extent these results generalize across species and experimental settings.

It should be noted that by using the GBLM we are only able to partially reproduce the fluctuations in the synaptic efficacy using the presynaptic firing activity, and many other factors could contribute to variation in synaptic efficacy. First, our model does not account for the long-term changes in the synaptic efficacy. Previous work has found that long-term plasticity can be induced in the visual cortex after repeated visual stimulation (Heynen and Bear, 2001; Frenkel et al., 2006; Cooke and Bear, 2010; Sale et al., 2011), as well as in hippocampus after behavioral training (Gruart et al., 2006; Whitlock et al., 2006; Fedulov et al., 2007). A model that accounts for the long-term plasticity may give a better prediction for the fluctuations in synaptic efficacy (Stevenson and Kording, 2011; Linderman et al., 2014; Song et al., 2018; Wei and Stevenson, 2021). Second, in our model, the STP rules themselves are assumed to be fixed for the whole recording. However, experiments have shown that the induction of long-term potentiation and long-term depression tend to shift STP itself to be more depressing or more facilitating, respectively (Markram and Tsodyks, 1996; Hardingham et al., 2007; Costa et al., 2015, 2017), which may allow for more flexibility in signal transmission (Carvalho and Buonomano, 2011). Third, although we have a term in our model to describe the slow fluctuation in network shared by the presynaptic and postsynaptic neurons, it does not account for the fast synaptic effect of other potential presynaptic neurons (Harris et al., 2003). A postsynaptic neuron can receive thousands of synaptic inputs, and correlated presynaptic inputs can affect the synaptic responses (Salinas and Sejnowski, 2000; Reinartz et al., 2014; Lee et al., 2016). Finally, we model STP as a function of presynaptic ISI to directly examine the influence of presynaptic firing, but there are also many alternative models for short-term synaptic dynamics (Hennig, 2013). Our results find that ∼20% of the variance in efficacy explained by the GBLM, but better models may be able to improve on this result and more accurately predict fluctuations in synaptic efficacy in vivo.

Using extracellular spike recordings from behaving mice, we find that the substantial fluctuations in the efficacy of putative excitatory synapses can be partially predicted by modeling short-term synaptic plasticity. Large-scale spike recordings have the advantage of allowing us to study putative connections between many brain areas for long periods of time (>2 h). Data from a wide range of synapses across the entire brain, may drive richer theoretical and normative models of synaptic transmission (Pfister et al., 2010; Hennig, 2013). However, studying synapses from spikes has the disadvantage that it requires sufficient numbers of presynaptic spikes to accurately estimate efficacy. Advances in stimulation and recording techniques may make it possible to study fluctuations on faster timescales (Kwan and Dan, 2012; English et al., 2017; Walker et al., 2021; Hage et al., 2022), and systematic fluctuations in efficacy for inhibitory synapses may also occur (Reyes et al., 1998; Beierlein et al., 2003; Ma et al., 2012). Although disentangling the complex relationships between synaptic efficacy and external or internal variables in vivo is an ongoing challenge in systems neuroscience, our results here highlight the potentially important role of STP in driving slow fluctuations in synaptic efficacy. In putative excitatory synapses from a wide range of brain areas, we find that short-term synaptic dynamics may lead to large, predictable fluctuations on longer timescales.

Footnotes

  • This work was supported by National Science Foundation Grant 1651396. We thank Harvey Swadlow and Monty Escabi for discussions, and the staff at Allen Institute for Brain Science for sharing the Visual Coding–Neuropixels dataset and for supporting open science.

  • The authors declare no competing financial interests.

  • Correspondence should be addressed to Naixin Ren at naixin.ren{at}uconn.edu

References

  1. ↵
    2. Abbott LF,
    3. Varela JA,
    4. Sen K,
    5. Nelson SB
    (1997) Synaptic depression and cortical gain control. Science 275:220–224. doi:10.1126/science.275.5297.221 pmid:8985017
    OpenUrlCrossRefPubMed
  2. ↵
    2. Abbott LF,
    3. Nelson SB
    (2000) Synaptic plasticity: taming the beast. Nat Neurosci 3:1178–1183. doi:10.1038/81453
    OpenUrlCrossRefPubMed
  3. ↵
    2. Amarasingham A,
    3. Harrison MT,
    4. Hatsopoulos NG,
    5. Geman S
    (2012) Conditional modeling and the jitter method of spike resampling. J Neurophysiol 107:517–531. doi:10.1152/jn.00633.2011 pmid:22031767
    OpenUrlCrossRefPubMed
  4. ↵
    2. Barthó P,
    3. Hirase H,
    4. Monconduit L,
    5. Zugaro M,
    6. Harris KD,
    7. Buzsáki G
    (2004) Characterization of neocortical principal cells and interneurons by network interactions and extracellular features. J Neurophysiol 92:600–608. doi:10.1152/jn.01170.2003
    OpenUrlCrossRefPubMed
  5. ↵
    2. Beierlein M,
    3. Gibson JR,
    4. Connors BW
    (2003) Two dynamically distinct inhibitory networks in layer 4 of the neocortex. J Neurophysiol 90:2987–3000. doi:10.1152/jn.00283.2003 pmid:12815025
    OpenUrlCrossRefPubMed
  6. ↵
    2. Ben-Yishai R,
    3. Bar-Or RL,
    4. Sompolinsky H
    (1995) Theory of orientation tuning in visual cortex. Proc Natl Acad Sci U S A 92:3844–3848. doi:10.1073/pnas.92.9.3844 pmid:7731993
    OpenUrlAbstract/FREE Full Text
  7. ↵
    2. Blackman AV,
    3. Abrahamsson T,
    4. Costa RP,
    5. Lalanne T,
    6. Sjöström PJ
    (2013) Target-cell-specific short-term plasticity in local circuits. Front Synaptic Neurosci 5:11.
    OpenUrlCrossRefPubMed
  8. ↵
    2. Bruno RM,
    3. Sakmann B
    (2006) Cortex is driven by weak but synchronously active thalamocortical synapses. Science 312:1622–1627. doi:10.1126/science.1124593 pmid:16778049
    OpenUrlAbstract/FREE Full Text
  9. ↵
    2. Buzsáki G,
    3. Mizuseki K
    (2014) The log-dynamic brain: how skewed distributions affect network operations. Nat Rev Neurosci 15:264–278. doi:10.1038/nrn3687 pmid:24569488
    OpenUrlCrossRefPubMed
  10. ↵
    2. Campagnola L, et al
    . (2022) Local connectivity and synaptic dynamics in mouse and human neocortex. Science 375:eabj5861.doi:10.1126/science.abj5861
    OpenUrlCrossRef
  11. ↵
    2. Carvalho TP,
    3. Buonomano DV
    (2011) A novel learning rule for long-term plasticity of short-term synaptic plasticity enhances temporal processing. Front Integr Neurosci 5:20. doi:10.3389/fnint.2011.00020 pmid:21660100
    OpenUrlCrossRefPubMed
  12. ↵
    2. Chung S,
    3. Li X,
    4. Nelson SB
    (2002) Short-term depression at thalamocortical synapses contributes to rapid adaptation of cortical sensory responses in vivo. Neuron 34:437–446. doi:10.1016/s0896-6273(02)00659-1 pmid:11988174
    OpenUrlCrossRefPubMed
  13. ↵
    2. Churchland MM,
    3. Cunningham JP,
    4. Kaufman MT,
    5. Foster JD,
    6. Nuyujukian P,
    7. Ryu SI,
    8. Shenoy KV
    (2012) Neural population dynamics during reaching. Nature 487:51–56. doi:10.1038/nature11129 pmid:22722855
    OpenUrlCrossRefPubMed
  14. ↵
    2. Connors BW,
    3. Gutnick MJ
    (1990) Intrinsic firing patterns of diverse neocortical neurons. Trends Neurosci 13:99–104. doi:10.1016/0166-2236(90)90185-d pmid:1691879
    OpenUrlCrossRefPubMed
  15. ↵
    2. Cooke SF,
    3. Bear MF
    (2010) Visual experience induces long-term potentiation in the primary visual cortex. J Neurosci 30:16304–16313. doi:10.1523/JNEUROSCI.4333-10.2010 pmid:21123576
    OpenUrlAbstract/FREE Full Text
  16. ↵
    2. Costa RP,
    3. Jesper Sjöström P,
    4. van Rossum MCW
    (2013) Probabilistic inference of short-term synaptic plasticity in neocortical microcircuits. Front Comput Neurosci 7:75. pmid:23761760
    OpenUrlCrossRefPubMed
  17. ↵
    2. Costa RP,
    3. Froemke RC,
    4. Sjöström PJ,
    5. van Rossum MCW
    (2015) Unified pre- and postsynaptic long-term plasticity enables reliable and flexible learning. Elife 4:e09457. doi:10.7554/eLife.09457
    OpenUrlCrossRefPubMed
  18. ↵
    2. Costa RP,
    3. Mizusaki BEP,
    4. Sjöström PJ,
    5. van Rossum MCW
    (2017) Functional consequences of pre- and postsynaptic expression of synaptic plasticity. Phil Trans R Soc B 372:20160153. doi:10.1098/rstb.2016.0153
    OpenUrlCrossRefPubMed
  19. ↵
    2. Deng PY,
    3. Klyachko VA
    (2011) The diverse functions of short-term plasticity components in synaptic computations. Commun Integr Biol 4:543–548. doi:10.4161/cib.4.5.15870 pmid:22046457
    OpenUrlCrossRefPubMed
  20. ↵
    2. English DF,
    3. Mckenzie S,
    4. Evans T,
    5. Kim K,
    6. Yoon E,
    7. Buzsáki G
    (2017) Pyramidal cell-interneuron circuit architecture and dynamics in hippocampal networks. Neuron 96:505–520.e7. doi:10.1016/j.neuron.2017.09.033
    OpenUrlCrossRefPubMed
  21. ↵
    2. Faisal AA,
    3. Selen LPJ,
    4. Wolpert DM
    (2008) Noise in the nervous system. Nat Rev Neurosci 9:292–303. doi:10.1038/nrn2258 pmid:18319728
    OpenUrlCrossRefPubMed
  22. ↵
    2. Fedulov V,
    3. Rex CS,
    4. Simmons DA,
    5. Palmer L,
    6. Gall CM,
    7. Lynch G
    (2007) Evidence that long-term potentiation occurs within individual hippocampal synapses during learning. J Neurosci 27:8031–8039. doi:10.1523/JNEUROSCI.2003-07.2007 pmid:17652593
    OpenUrlAbstract/FREE Full Text
  23. ↵
    2. Fetz E,
    3. Toyama K,
    4. Smith W
    (1991) Synaptic interactions between cortical neurons. In: Normal and altered states of function. Cerebral cortex (Peters A, Jones EG, eds), pp 1–47. Boston: Springer.
  24. ↵
    2. Fortune ES,
    3. Rose GJ
    (2001) Short-term synaptic plasticity as a temporal filter. Trends Neurosci 24:381–385. doi:10.1016/s0166-2236(00)01835-x pmid:11410267
    OpenUrlCrossRefPubMed
  25. ↵
    2. Frenkel MY,
    3. Sawtell NB,
    4. Diogo ACM,
    5. Yoon B,
    6. Neve RL,
    7. Bear MF
    (2006) Instructive effect of visual experience in mouse visual cortex. Neuron 51:339–349. doi:10.1016/j.neuron.2006.06.026 pmid:16880128
    OpenUrlCrossRefPubMed
  26. ↵
    2. Fujisawa S,
    3. Amarasingham A,
    4. Harrison MT,
    5. Buzsáki G
    (2008) Behavior-dependent short-term assembly dynamics in the medial prefrontal cortex. Nat Neurosci 11:823–833. doi:10.1038/nn.2134
    OpenUrlCrossRefPubMed
  27. ↵
    2. Gal A,
    3. Eytan D,
    4. Wallach A,
    5. Sandler M,
    6. Schiller J,
    7. Marom S
    (2010) Dynamics of excitability over extended timescales in cultured cortical neurons. J Neurosci 30:16332–16342. doi:10.1523/JNEUROSCI.4859-10.2010 pmid:21123579
    OpenUrlAbstract/FREE Full Text
  28. ↵
    2. Ghanbari A,
    3. Malyshev A,
    4. Volgushev M,
    5. Stevenson IH
    (2017) Estimating short-term synaptic plasticity from pre- and postsynaptic spiking. PLOS Comput Biol 13:e1005738. doi:10.1371/journal.pcbi.1005738 pmid:28873406
    OpenUrlCrossRefPubMed
  29. ↵
    2. Ghanbari A,
    3. Ren N,
    4. Keine C,
    5. Stoelzel C,
    6. Englitz B,
    7. Swadlow HA,
    8. Stevenson IH
    (2020) Modeling the short-term dynamics of in vivo excitatory spike transmission. J Neurosci 40:4185–4202. doi:10.1523/JNEUROSCI.1482-19.2020 pmid:32303648
    OpenUrlAbstract/FREE Full Text
  30. ↵
    2. Gilbert CD,
    3. Sigman M
    (2007) Brain states: top-down influences in sensory processing. Neuron 54:677–696. doi:10.1016/j.neuron.2007.05.019 pmid:17553419
    OpenUrlCrossRefPubMed
  31. ↵
    2. Goris RLT,
    3. Movshon JA,
    4. Simoncelli EP
    (2014) Partitioning neuronal variability. Nat Neurosci 17:858–865. doi:10.1038/nn.3711 pmid:24777419
    OpenUrlCrossRefPubMed
  32. ↵
    2. Gruart A,
    3. Muñoz MD,
    4. Delgado-García JM
    (2006) Involvement of the CA3–CA1 synapse in the acquisition of associative learning in behaving mice. J Neurosci 26:1077–1087. doi:10.1523/JNEUROSCI.2834-05.2006 pmid:16436593
    OpenUrlAbstract/FREE Full Text
  33. ↵
    2. Hage TA,
    3. Bosma-Moody A,
    4. Baker CA,
    5. Kratz M,
    6. Campagnola L,
    7. Jarsky T,
    8. Zeng H,
    9. Murphy GJ
    (2022) Synaptic connectivity to L2/3 of primary visual cortex measured by two-photon optogenetic stimulation. Elife 11:e71103.doi:10.7554/eLife.71103
    OpenUrlCrossRef
  34. ↵
    2. Hardingham NR,
    3. Hardingham GE,
    4. Fox KD,
    5. Jack JJB
    (2007) Presynaptic efficacy directs normalization of synaptic strength in layer 2/3 rat neocortex after paired activity. J Neurophysiol 97:2965–2975. doi:10.1152/jn.01352.2006 pmid:17267749
    OpenUrlCrossRefPubMed
  35. ↵
    2. Harris KD,
    3. Csicsvari J,
    4. Hirase H,
    5. Dragoi G,
    6. Buzsáki G
    (2003) Organization of cell assemblies in the hippocampus. Nature 424:552–556. doi:10.1038/nature01834 pmid:12891358
    OpenUrlCrossRefPubMed
  36. ↵
    2. Hennig MH
    (2013) Theoretical models of synaptic short term plasticity. Front Comput Neurosci 7:45. doi:10.3389/fncom.2013.00045 pmid:23626536
    OpenUrlCrossRefPubMed
  37. ↵
    2. Heynen AJ,
    3. Bear MF
    (2001) Long-term potentiation of thalamocortical transmission in the adult visual cortex in vivo. J Neurosci 21:9801–9813. pmid:11739588
    OpenUrlAbstract/FREE Full Text
  38. ↵
    2. Hubel DH,
    3. Wiesel TN
    (1959) Receptive fields of single neurones in the cat's striate cortex. J Physiol 148:574–591. doi:10.1113/jphysiol.1959.sp006308 pmid:14403679
    OpenUrlCrossRefPubMed
  39. ↵
    2. Ibata K,
    3. Sun Q,
    4. Turrigiano GG
    (2008) Rapid synaptic scaling induced by changes in postsynaptic firing. Neuron 57:819–826. doi:10.1016/j.neuron.2008.02.031 pmid:18367083
    OpenUrlCrossRefPubMed
  40. ↵
    2. Kandaswamy U,
    3. Deng PY,
    4. Stevens CF,
    5. Klyachko VA
    (2010) The role of presynaptic dynamics in processing of natural spike trains in hippocampal synapses. J Neurosci 30:15904–15914. doi:10.1523/JNEUROSCI.4050-10.2010 pmid:21106829
    OpenUrlAbstract/FREE Full Text
  41. ↵
    2. Kang K,
    3. Shapley RM,
    4. Sompolinsky H
    (2004) Information tuning of populations of neurons in primary visual cortex. J Neurosci 24:3726–3735. doi:10.1523/JNEUROSCI.4272-03.2004 pmid:15084652
    OpenUrlAbstract/FREE Full Text
  42. ↵
    2. Klyachko VA,
    3. Stevens CF
    (2006) Excitatory and feed-forward inhibitory hippocampal synapses work synergistically as an adaptive filter of natural spike trains. PLoS Biol 4:e207. doi:10.1371/journal.pbio.0040207
    OpenUrlCrossRefPubMed
  43. ↵
    2. Kobayashi R,
    3. Kurita S,
    4. Kurth A,
    5. Kitano K,
    6. Mizuseki K,
    7. Diesmann M,
    8. Richmond BJ,
    9. Shinomoto S
    (2019) Reconstructing neuronal circuitry from parallel spike trains. Nat Commun 10:4468. doi:10.1038/s41467-019-12225-2
    OpenUrlCrossRef
  44. ↵
    2. Kwan AC,
    3. Dan Y
    (2012) Dissection of cortical microcircuits by single-neuron stimulation in vivo. Curr Biol 22:1459–1467. doi:10.1016/j.cub.2012.06.007 pmid:22748320
    OpenUrlCrossRefPubMed
  45. ↵
    2. Lee JH,
    3. Campagnola L,
    4. Seeman SC,
    5. Jarsky T,
    6. Mihalas S
    (2019) Functional synapse types via characterization of short-term synaptic plasticity. bioRxiv 648725. doi:10.1101/648725
    OpenUrlCrossRef
  46. ↵
    2. Lee KFH,
    3. Soares C,
    4. Thivierge JP,
    5. Béïque JC
    (2016) Correlated synaptic inputs drive dendritic calcium amplification and cooperative plasticity during clustered synapse development. Neuron 89:784–799. doi:10.1016/j.neuron.2016.01.012 pmid:26853305
    OpenUrlCrossRefPubMed
  47. ↵
    2. Levick WR,
    3. Cleland BG,
    4. Dubin MW
    (1972) Lateral geniculate neurons of cat: retinal inputs and physiology. Invest Ophthalmol 11:302–311. pmid:5028229
    OpenUrlAbstract/FREE Full Text
  48. ↵
    2. Linderman SW,
    3. Stock CH,
    4. Adams RP
    (2014) A framework for studying synaptic plasticity with neural spike train data. arXiv:1411.4077. doi:10.48550/arXiv.1411.4077
  49. ↵
    2. Ma Y,
    3. Hu H,
    4. Agmon A
    (2012) Short-term plasticity of unitary inhibitory-to-inhibitory synapses depends on the presynaptic interneuron subtype. J Neurosci 32:983–988. doi:10.1523/JNEUROSCI.5007-11.2012 pmid:22262896
    OpenUrlAbstract/FREE Full Text
  50. ↵
    2. Markram H,
    3. Tsodyks M
    (1996) Redistribution of synaptic efficacy between neocortical pyramidal neurons. Nature 382:807–810. doi:10.1038/382807a0 pmid:8752273
    OpenUrlCrossRefPubMed
  51. ↵
    2. McGinley MJ,
    3. Vinck M,
    4. Reimer J,
    5. Batista-Brito R,
    6. Zagha E,
    7. Cadwell CR,
    8. Tolias AS,
    9. Cardin JA,
    10. McCormick DA
    (2015) Waking state: rapid variations modulate neural and behavioral responses. Neuron 87:1143–1161. doi:10.1016/j.neuron.2015.09.012 pmid:26402600
    OpenUrlCrossRefPubMed
  52. ↵
    2. McKenzie S,
    3. Huszár R,
    4. English DF,
    5. Kim K,
    6. Christensen F,
    7. Yoon E,
    8. Buzsáki G
    (2021) Preexisting hippocampal network dynamics constrain optogenetically induced place fields. Neuron 109:1040–1054.e7. doi:10.1016/j.neuron.2021.01.011 pmid:33539763
    OpenUrlCrossRefPubMed
  53. ↵
    2. Mochizuki Y, et al
    . (2016) Similarity in neuronal firing regimes across mammalian species. J Neurosci 36:5736–5747. doi:10.1523/JNEUROSCI.0230-16.2016 pmid:27225764
    OpenUrlAbstract/FREE Full Text
  54. ↵
    2. Pala A,
    3. Petersen CCH
    (2015) In vivo measurement of cell-type-specific synaptic connectivity and synaptic transmission in layer 2/3 mouse barrel cortex. Neuron 85:68–75. doi:10.1016/j.neuron.2014.11.025 pmid:25543458
    OpenUrlCrossRefPubMed
  55. ↵
    2. Pala A,
    3. Petersen CCH
    (2018) State-dependent cell-type-specific membrane potential dynamics and unitary synaptic inputs in awake mice. Elife 7:e35869. doi:10.7554/eLife.35869
    OpenUrlCrossRef
  56. ↵
    2. Perkel DH,
    3. Gerstein GL,
    4. Moore GP
    (1967) Neuronal spike trains and stochastic point processes: II. Simultaneous spike trains. Biophys J 7:419–440. doi:10.1016/S0006-3495(67)86597-4 pmid:4292792
    OpenUrlCrossRefPubMed
  57. ↵
    2. Pfister JP,
    3. Dayan P,
    4. Lengyel M
    (2010) Synapses with short-term plasticity are optimal estimators of presynaptic membrane potentials. Nat Neurosci 13:1271–1275. doi:10.1038/nn.2640 pmid:20852625
    OpenUrlCrossRefPubMed
  58. ↵
    2. Pillow JW,
    3. Shlens J,
    4. Paninski L,
    5. Sher A,
    6. Litke AM,
    7. Chichilnisky EJ,
    8. Simoncelli EP
    (2008) Spatio-temporal correlations and visual signalling in a complete neuronal population. Nature 454:995–999. doi:10.1038/nature07140 pmid:18650810
    OpenUrlCrossRefPubMed
  59. ↵
    2. Poulet JFA,
    3. Petersen CCH
    (2008) Internal brain state regulates membrane potential synchrony in barrel cortex of behaving mice. Nature 454:881–885. doi:10.1038/nature07150 pmid:18633351
    OpenUrlCrossRefPubMed
  60. ↵
    2. Rebesco JM,
    3. Stevenson IH,
    4. Körding KP,
    5. Solla SA,
    6. Miller LE
    (2010) Rewiring neural interactions by micro-stimulation. Front Syst Neurosci 4:39.
    OpenUrlCrossRefPubMed
  61. ↵
    2. Reinartz S,
    3. Biro I,
    4. Gal A,
    5. Giugliano M,
    6. Marom S
    (2014) Synaptic dynamics contribute to long-term single neuron response fluctuations. Front Neural Circuits 8:71. doi:10.3389/fncir.2014.00071 pmid:25071452
    OpenUrlCrossRefPubMed
  62. ↵
    2. Ren N,
    3. Ito S,
    4. Hafizi H,
    5. Beggs JM,
    6. Stevenson IH
    (2020) Model-based detection of putative synaptic connections from spike recordings with latency and type constraints. J Neurophysiol 124:1588–1604.
    OpenUrl
  63. ↵
    2. Reyes A,
    3. Lujan R,
    4. Rozov A,
    5. Burnashev N,
    6. Somogyi P,
    7. Sakmann B
    (1998) Target-cell-specific facilitation and depression in neocortical circuits. Nat Neurosci 1:279–85285. doi:10.1038/1092 pmid:10195160
    OpenUrlCrossRefPubMed
  64. ↵
    2. Sale A,
    3. De Pasquale R,
    4. Bonaccorsi J,
    5. Pietra G,
    6. Olivieri D,
    7. Berardi N,
    8. Maffei L
    (2011) Visual perceptual learning induces long-term potentiation in the visual cortex. Neuroscience 172:219–225. doi:10.1016/j.neuroscience.2010.10.078 pmid:21056088
    OpenUrlCrossRefPubMed
  65. ↵
    2. Salinas E,
    3. Sejnowski TJ
    (2000) Impact of correlated synaptic input on output firing rate and variability in simple neuronal models. J Neurosci 20:6193–6209. pmid:10934269
    OpenUrlAbstract/FREE Full Text
  66. ↵
    2. Sedigh-Sarvestani M,
    3. Vigeland L,
    4. Fernandez-Lamo I,
    5. Taylor MM,
    6. Palmer LA,
    7. Contreras D
    (2017) Intracellular, in vivo, dynamics of thalamocortical synapses in visual cortex. J Neurosci 37:5250–5262. doi:10.1523/JNEUROSCI.3370-16.2017 pmid:28438969
    OpenUrlAbstract/FREE Full Text
  67. ↵
    2. Seeman SC, et al
    . (2018) Sparse recurrent excitatory connectivity in the microcircuit of the adult mouse and human cortex. Elife 7:doi:10.7554/eLife.37349
    OpenUrlCrossRefPubMed
  68. ↵
    2. Shenoy KV,
    3. Sahani M,
    4. Churchland MM
    (2013) Cortical control of arm movements: a dynamical systems perspective. Annu Rev Neurosci 36:337–359. doi:10.1146/annurev-neuro-062111-150509 pmid:23725001
    OpenUrlCrossRefPubMed
  69. ↵
    2. Shimazaki H,
    3. Shinomoto S
    (2010) Kernel bandwidth optimization in spike rate estimation. J Comput Neurosci 29:171–182. doi:10.1007/s10827-009-0180-4 pmid:19655238
    OpenUrlCrossRefPubMed
  70. ↵
    2. Siegle JH, et al
    . (2021) Survey of spiking in the mouse visual system reveals functional hierarchy. Nature 592:86–92. doi:10.1038/s41586-020-03171-x pmid:33473216
    OpenUrlCrossRefPubMed
  71. ↵
    2. Sjöström PJ,
    3. Rancz EA,
    4. Roth A,
    5. Häusser M
    (2008) Dendritic excitability and synaptic plasticity. Physiol Rev 88:769–840. doi:10.1152/physrev.00016.2007 pmid:18391179
    OpenUrlCrossRefPubMed
  72. ↵
    2. Song D,
    3. Robinson BS,
    4. Berger TW
    (2018) Identification of short-term and long-term functional synaptic plasticity from spiking activities. In: Adaptive learning methods for nonlinear system modeling (Comminiello D, Principe JC, eds), pp 289–312. Oxford: Elsevier.
  73. ↵
    2. Stevenson IH,
    3. Kording KP
    (2011) Inferring spike-timing-dependent plasticity from spike train data. Advances in neural information processing systems, 24.
  74. ↵
    2. Stoelzel CR,
    3. Bereshpolova Y,
    4. Gusev AG,
    5. Swadlow HA
    (2008) The impact of an LGNd impulse on the awake visual cortex: synaptic dynamics and the sustained/transient distinction. J Neurosci 28:5018–5028. doi:10.1523/JNEUROSCI.4726-07.2008 pmid:18463255
    OpenUrlAbstract/FREE Full Text
  75. ↵
    2. Stoelzel CR,
    3. Bereshpolova Y,
    4. Swadlow HA
    (2009) Stability of thalamocortical synaptic transmission across awake brain states. J Neurosci 29:6851–6859. doi:10.1523/JNEUROSCI.5983-08.2009 pmid:19474312
    OpenUrlAbstract/FREE Full Text
  76. ↵
    2. Stoelzel CR,
    3. Huff JM,
    4. Bereshpolova Y,
    5. Zhuang J,
    6. Hei X,
    7. Alonso JM,
    8. Swadlow HA
    (2015) Hour-long adaptation in the awake early visual system. J Neurophysiol 114:1172–1182. doi:10.1152/jn.00116.2015 pmid:26108950
    OpenUrlCrossRefPubMed
  77. ↵
    2. Swadlow HA,
    3. Gusev AG
    (2001) The impact of “bursting” thalamic impulses at a neocortical synapse. Nat Neurosci 4:402–408. doi:10.1038/86054 pmid:11276231
    OpenUrlCrossRefPubMed
  78. ↵
    2. Swadlow HA,
    3. Gusev AG
    (2002) Receptive-field construction in cortical inhibitory interneurons. Nat Neurosci 5:403–404. doi:10.1038/nn847 pmid:11967546
    OpenUrlCrossRefPubMed
  79. ↵
    2. Teich MC,
    3. Johnson DH,
    4. Kumar AR,
    5. Turcott RG
    (1990) Rate fluctuations and fractional power-law noise recorded from cells in the lower auditory pathway of the cat. Hear Res 46:41–52. doi:10.1016/0378-5955(90)90138-f pmid:2380126
    OpenUrlCrossRefPubMed
  80. ↵
    2. Thomson AM,
    3. Lamy C
    (2007) Functional maps of neocortical local circuitry. Front Neurosci 1:19–42. doi:10.3389/neuro.01.1.1.002.2007 pmid:18982117
    OpenUrlCrossRefPubMed
  81. ↵
    2. Trainito C,
    3. von Nicolai C,
    4. Miller EK,
    5. Siegel M
    (2019) Extracellular spike waveform dissociates four functionally distinct cell classes in primate cortex. Curr Biol 29:2973–2982.e5. doi:10.1016/j.cub.2019.07.051 pmid:31447374
    OpenUrlCrossRefPubMed
  82. ↵
    2. Truccolo W,
    3. Eden UT,
    4. Fellows MR,
    5. Donoghue JP,
    6. Brown EN
    (2005) A point process framework for relating neural spiking activity to spiking history, neural ensemble, and extrinsic covariate effects. J Neurophysiol 93:1074–1089. doi:10.1152/jn.00697.2004 pmid:15356183
    OpenUrlCrossRefPubMed
  83. ↵
    2. Tsodyks M,
    3. Kenet T,
    4. Grinvald A,
    5. Arieli A
    (1999) Linking spontaneous activity of single cortical neurons and the underlying functional architecture. Science 286:1943–1946. doi:10.1126/science.286.5446.1943 pmid:10583955
    OpenUrlAbstract/FREE Full Text
  84. ↵
    2. Tsodyks MV,
    3. Markram H
    (1997) The neural code between neocortical pyramidal neurons depends on neurotransmitter release probability. Proc Natl Acad Sci U S A 94:719–723. doi:10.1073/pnas.94.2.719 pmid:9012851
    OpenUrlAbstract/FREE Full Text
  85. ↵
    2. Usrey WM,
    3. Alonso JM,
    4. Reid RC
    (2000) Synaptic interactions between thalamic inputs to simple cells in cat visual cortex. J Neurosci 20:5461–5467. pmid:10884329
    OpenUrlAbstract/FREE Full Text
  86. ↵
    2. Vogels TP,
    3. Rajan K,
    4. Abbott LF
    (2005) Neural network dynamics. Annu Rev Neurosci 28:357–376. doi:10.1146/annurev.neuro.28.061604.135637 pmid:16022600
    OpenUrlCrossRefPubMed
  87. ↵
    2. Walker AS,
    3. Raliski BK,
    4. Karbasi K,
    5. Zhang P,
    6. Sanders K,
    7. Miller EW
    (2021) Optical spike detection and connectivity analysis with a far-red voltage-sensitive fluorophore reveals changes to network connectivity in development and disease. Front Neurosci 15:643859.
    OpenUrl
  88. ↵
    2. Wei G,
    3. Stevenson IH
    (2021) Tracking fast and slow changes in synaptic weights from simultaneously observed pre- and postsynaptic spiking. Neural Comput 33:2682–2709. doi:10.1162/neco_a_01426 pmid:34530452
    OpenUrlCrossRefPubMed
  89. ↵
    2. Whitlock JR,
    3. Heynen AJ,
    4. Shuler MG,
    5. Bear MF
    (2006) Learning induces long-term potentiation in the hippocampus. Science 313:1093–1097. doi:10.1126/science.1128134 pmid:16931756
    OpenUrlAbstract/FREE Full Text
  90. ↵
    2. Yaksi E,
    3. Friedrich RW
    (2006) Reconstruction of firing rate changes across neuronal populations by temporally deconvolved Ca2+ imaging. Nat Methods 3:377–383. doi:10.1038/nmeth874 pmid:16628208
    OpenUrlCrossRefPubMed
  91. ↵
    2. Zucker RSRS,
    3. Regehr WGWG
    (2002) Short-term synaptic plasticity. Ann Rev Physiol 64:355–405.
    OpenUrlCrossRefPubMed
Back to top

In this issue

Email
Print
View Full Page PDF
Citation Tools
Respond to this article
Predictable Fluctuations in Excitatory Synaptic Strength Due to Natural Variation in Presynaptic Firing Rate
Naixin Ren, Ganchao Wei, Abed Ghanbari, Ian H. Stevenson
Journal of Neuroscience 16 November 2022, 42 (46) 8608-8620; DOI: 10.1523/JNEUROSCI.0808-22.2022
Twitter logo Facebook logo Mendeley logo

Jump to section

Keywords

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.