(C) PLOS One [1]. This unaltered content originally appeared in journals.plosone.org.
Licensed under Creative Commons Attribution (CC BY) license.
url:https://journals.plos.org/plosone/s/licenses-and-copyright

------------



STDP and the distribution of preferred phases in the whisker system

['Nimrod Sherf', 'Physics Department', 'Ben-Gurion University Of The Negev', 'Beer-Sheva', 'Zlotowski Center For Neuroscience', 'Maoz Shamir', 'Department Of Physiology', 'Cell Biology Faculty Of Health Sciences']

Date: 2021-10

Rats and mice use their whiskers to probe the environment. By rhythmically swiping their whiskers back and forth they can detect the existence of an object, locate it, and identify its texture. Localization can be accomplished by inferring the whisker’s position. Rhythmic neurons that track the phase of the whisking cycle encode information about the azimuthal location of the whisker. These neurons are characterized by preferred phases of firing that are narrowly distributed. Consequently, pooling the rhythmic signal from several upstream neurons is expected to result in a much narrower distribution of preferred phases in the downstream population, which however has not been observed empirically. Here, we show how spike timing dependent plasticity (STDP) can provide a solution to this conundrum. We investigated the effect of STDP on the utility of a neural population to transmit rhythmic information downstream using the framework of a modeling study. We found that under a wide range of parameters, STDP facilitated the transfer of rhythmic information despite the fact that all the synaptic weights remained dynamic. As a result, the preferred phase of the downstream neuron was not fixed, but rather drifted in time at a drift velocity that depended on the preferred phase, thus inducing a distribution of preferred phases. We further analyzed how the STDP rule governs the distribution of preferred phases in the downstream population. This link between the STDP rule and the distribution of preferred phases constitutes a natural test for our theory.

The distribution of preferred phases of whisking neurons in the somatosensory system of rats and mice presents a conundrum: a simple pooling model predicts a distribution that is an order of magnitude narrower than what is observed empirically. Here, we suggest that this non-trivial distribution may result from activity-dependent plasticity in the form of spike timing dependent plasticity (STDP). We show that under STDP, the synaptic weights do not converge to a fixed value, but rather remain dynamic. As a result, the preferred phases of the whisking neurons vary in time, hence inducing a non-trivial distribution of preferred phases, which is governed by the STDP rule. Our results imply that the considerable synaptic volatility which has long been viewed as a difficulty that needs to be overcome, may actually be an underlying principle of the organization of the central nervous system.

Copyright: © 2021 Sherf, Shamir. This is an open access article distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Here we investigated the hypothesis that the distribution of phases in the downstream layer is governed by the interplay of the distribution in the upstream layer and the STDP rule. The remainder of this article is organized as follows. First, we define the network architecture and the STDP learning rule. We then derive a mean-field approximation for the STDP dynamics in the limit of a slow learning rate for a threshold-linear Poisson downstream neuron model. Next, we show that STDP dynamics can generate non-trivial distributions of preferred phases in the downstream population and analyze how the parameters characterizing the STDP govern this distribution. Finally, we summarize the results, discuss the limitations of this study and suggest essential empirical predictions deriving from our theory.

Recently, the effects of rhythmic activity on the spike timing dependent plasticity (STDP) dynamics of feed-forward synaptic connections have been examined [ 22 , 23 ]. It was shown that in this case the synaptic weights remain dynamic. As a result, the phase of the downstream neuron is arbitrary and drifts in time; thus, effectively, inducing a distribution of preferred phases in the downstream population. However, if the phases of the downstream population are arbitrary and drift in time, how can information about the whisking phase be transmitted?

Fig 2b shows the distribution of preferred phases in a naive model that pools the responses of N = 100 VPM neurons. The source of variability in the preferred phases of the downstream neuron is the random choice of N VPM neurons, each characterized by a preferred phase that is distributed in an i.i.d. manner following Eq (1) . As a result, one would expect that the number, N, of pooled thalamic neurons would govern the width of the distribution. Fig 2c shows the expected distribution width, κ L4I , as a function of N for uniform (squares) and random (circles) pooling. As can be seen from the figure, even a random pooling of N = 10 results in a distribution of preferred phases that is considerably narrower than empirically observed. This result is particularly surprising since the number of thalamic neurons synapsing onto a single L4 neuron was estimated to be on the order of 100, see e.g. [ 16 , 21 ]. Moreover, in the naive pooling model, the mean preferred phase of the downstream population is given by the mean phase of the upstream population. The delay in the response of downstream neurons to their input should also be added. It is estimated to be around 1-10ms; hence, for whisking at 10Hz translates into ψ VPM − ψ L4I = 0.06-0.6 [rad]. This result also runs counter empirical findings reporting ψ VPM − ψ L4I ≈ 2.2 [rad]. Thus, the synaptic weight profile in a model that lacks activity-dependent plasticity is expected to converge to a random pooling (due to constant synaptic remodelling), which is inconsistent with the observed distribution of preferred phases.

Throughout this work we ignore the possible contribution of recurrent connections to the shaping of phase selectivity in layer 4 (but see Discussion ). Assuming the rhythmic input to L4I neurons originates solely from the VPM, the distribution of preferred phases of L4I neurons is determined by the distribution in the thalamus and the profile of thalamo-cortical synaptic weights. However, synaptic weights have been reported to be highly volatile [ 17 – 20 ]. If synaptic plasticity is activity independent, then one would expect the synaptic weight profile to become unstructured. In the naive pooling model, which lacks activity-dependent plasticity, we consider two types of unstructured connectivity; namely, uniform pooling and random pooling.

(a) The von-Mises distribution, Eq (1) , with ψ = 0 is shown for different values of κ as depicted by color. (b) The distribution of preferred phases in the uniform pooling model. The distribution was estimated from 10000 repetitions of drawing N = 100 neurons from the upstream population. The preferred phases in the upstream population were drawn in an i.i.d. manner following Eq (1) , with κ UpStream = 1 and ψ UpStream = 0. The dashed red line depicts the mean phase. (c) Distribution width in the naive uniform/random pooling model. The width of the distribution of preferred phases, κ L4I , in the downstream layer (L4I) is shown as a function of the number of pooled VPM neurons, N, for the uniform and random pooling in squares and circles, respectively. The width of the distribution, κ L4I , was estimated from 10000 repetitions of drawing N preferred phases of the upstream population with κ = 1 (blue) and κ = 2 (red).

(a) Mice and rats can infer the azimuthal location of an object by touch. (b) The angular position of a whisker, β, (during whisking) is shown as a function of time. The angle is often modeled as β(t) = β midpoint (t) + Δβ(t) cos(ϕ(t)), where β midpoint (t) and Δβ(t) are the midpoint and the whisking amplitude, respectively. (c) The whisking phase ϕ as a function of time is ϕ(t) = (νt) mod2π , where ν is the angular frequency of the whisking. (d) & (e) Raster plot and normalized tuning curve of a neuron with a preferred phase near maximal retraction.

During whisking, the animal moves its vibrissae back and forth in a rhythmic manner, Fig 1a–1c . Neurons that track the azimuthal position of the whisker by firing in a preferential manner to the phase of the whisking cycle are termed whisking neurons. Whisking neurons in the ventral posteromedial nucleus (VPM) of the thalamus as well as inhibitory whisking neurons in layer 4 of the barrel cortex are characterized by a preferred phase at which they fire with the highest rate during the whisking cycle [ 1 , 6 , 12 – 16 ], Fig 1d and 1e . The distribution of preferred phases is non-uniform and can be approximated by the circular normal (Von-Mises) distribution (1) where ψ is the mean phase and I 0 (κ) is the modified Bessel function of order 0, Fig 2a . The parameter κ quantifies the width of the distribution; κ = 0 yields a uniform (flat) distribution, whereas in the limit of κ → ∞ the distribution converges to a delta function. Typical values for κ in the thalamus and for layer 4 inhibitory (L4I) whisking neurons are κ VPM ≈ κ L4I ≈ 1 where ψ VPM ≈ 5π/6 [rad] and κ L4I ≈ 0.5 [rad] [ 2 , 13 ].

The whisker system is used by rats and mice to actively gather information about their proximal environment [ 1 – 4 ]. Information about whisker position, touch events, and texture is relayed downstream the somatosensory system via several tracks; in particular, the lemniscal pathway that relays information about both whisking and touch [ 5 – 11 ].

The distribution of preferred phases of L4I neurons, κ L4I , is shown as a function of the characteristic timescales of the STDP: τ − in (a) & (c) and τ − + in (b) & (d), for the temporally asymmetric rule in (a) & (b) and the symmetric rule in (c) & (d). The dashed black lines depict the mean phase. Unless stated otherwise, the parameters used here were: κ VPM = 1, γ = 1, D = 10hz, τ − = 50ms, τ + = 22ms, d = 10ms, ψ VPM = 5π/6, μ = 0.01, α = 1.1 and λ = 0.01.

Fig 11 shows the effect of the temporal structure of the STDP on the distribution of preferred phases in the downstream population. As can be seen from the figure, varying the characteristic timescales of potentiation and depression, τ + and τ − , induces quantitative effects on the distribution of the preferred phases. However, qualitative effects result from changing the nature of the STDP rule. Comparing the temporally asymmetric rule, Fig 11a and 11b , to the temporally symmetric rule, Fig 11c and 11d , reveals a dramatic difference in the mean preferred phase of L4I neurons, dashed black lines.

(a) The mean input phase κ L4I as a function of α. Probability of occurrences is depicted by color as shown in the right color bar. (b) The order parameters and are shown as a function of α. Here, we used: κ VPM = 1, ψ 0 = 5π/6, γ = 1, D = 10hz, and d = 3ms. The temporally asymmetric STDP rule, Eq (9) , was used with: τ − = 50ms, τ + = 22ms, μ = 0.01 and λ = 0.01.

As one might expect, decreasing α beyond a certain value will result in a potentiation dominated STDP dynamics, saturating all synapses close to their maximal value, thus, approaching the uniform solution, which is characterized by a narrow preferred phase distribution centered around the mean preferred phase of the VPM neurons shifted by the delay, and low values of γ L4I , Fig 10 . Increasing α strengthens the competitive winner-take-all like nature of the STDP dynamics. Initially, this competition will generate a fixed point with non-uniform synaptic weights, thus increasing γ L4I . A further increase of α results in a limit cycle solution to the STDP dynamics that widens the distribution of the preferred phases. Increasing α beyond a certain critical value will result in depression dominated STDP dynamics, thus driving the synaptic weights to zero, Fig 10 .

(a) The distribution of preferred phases of L4I neurons, κ L4I , is shown by color as a function of μ. (b) The modulation depth, γ L4I , is depicted as a function of μ. The parameters used here were: κ VPM = 1, ψ VPM = 5π/6, γ = 0.9, D = 10hz, and d = 3ms. The temporally asymmetric STDP rule, Eq (9) , was used with: τ − = 50ms, τ + = 22ms, α = 1.1 and λ = 0.01.

(a), (c), (e) & (f) show the synaptic weight dynamics for different values of μ. The synaptic weights are depicted by different traces, colored according to the preferred phase of the pre-synaptic neuron, see legend. (b), (d), (f) & (h) show the preferred phase of the downstream L4I neuron, κ L4I , as a function of time. The parameters used in these simulation were: N = 150, κ = 1, , γ = 0.9, D = 10Hz, and d = 3ms. We used here the temporally asymmetric learning rule, Eq (9) , with: τ − = 50ms, τ + = 22ms, and α = 1.1. For the non-linearity parameter we used: μ = 0.01 in (a) & (b), μ = 0.06. in (c) & (d), μ = 0.07 in (e) & (f). and μ = 0.1 in (g) & (h).

(a) The (non-normalized) distribution of L4I neuron phases, κ L4I , is depicted by color as a function, . (b) The width of the distribution, κ L4I , of L4I neurons, is shown as a function of . Unless stated otherwise, the parameters used here were: κ VPM = 1, ψ VPM = 5π/6, γ = 0.9, D = 10hz, and d = 10ms. We used the temporally asymmetric exponential learning rule, Eq (10) , with: τ − = 50ms, τ + = 22ms, μ = 0.01, α = 1.1, and λ = 0.01.

(a) The (non-normalized) distribution of L4I neuron phases, κ L4I , is depicted by color as a function, . (b) The width of the distribution, κ L4I , of L4I neurons is shown as a function of . The parameters used in these graphs were: κ VPM = 1, ψ VPM = 5π/6 γ = 0.9, D = 10hz, and d = 3ms. We used the temporally asymmetric exponential learning rule, Eq (9) , with τ − = 50ms, τ + = 22ms, μ = 0.01, α = 1.1, and λ = 0.01.

The effect of the whisking frequency is shown in Fig 6a . For moderate rhythms that are on a similar timescale to that of the STDP rule, STDP dynamics can generate a wide distribution of preferred phases. However, in the high frequency limit, ν → ∞, the synaptic weights converge to a uniform solution with w(ϕ) = (1 + α 1/μ ) −1 , ∀ϕ (see [ 22 ]). In this limit, due to the uniform pooling, there is no selective amplification of phases and the whisking signal is transmitted downstream due to the selectivity of the thalamic population, κ VPM > 0. Consequently, at high frequencies, the distribution of preferred phases will be extremely narrow, because of the quenched disorder, and the rhythmic signal will be attenuated, γ L4I = γ VPM × I 1 (κ)/I 0 (κ) (where I j (κ) is the modified Bessel function of order j), Fig 6b . The rate of convergence to the high frequency limit is governed by the smoothness of the STDP rule. Discontinuity in the STPD rule, such as in our choice of a temporally asymmetric rule, will induce algebraic convergence in ν to the high frequency limit, whereas a smooth rule, such as our choice of a temporally symmetric rule, will manifest in exponential convergence; compare Fig 6 with Fig 7 , see also [ 22 ]. In our choice of parameters, τ + ≈ 20ms and τ − ≈ 50ms, STDP dynamics induced a wide distribution for frequencies in the α, β and low γ bands in the case of the asymmetric learning rule Fig 6 , and around 5–15hz in the case of the symmetric learning rule Fig 7 .

(a) The (non-normalized) distribution of the preferred phases of L4I neurons, κ L4I , is shown by color as a function of the width of the distribution of preferred phases in VPM, κ VPM . (b) The distribution width in layer 4, κ L4I , is shown as a function of κ VPM (blue). The identity line is shown (dashed black) for comparison. (c) The modulation depth in the downstream population, γ L4I , is shown as a function of κ VPM (blue). For comparison, the modulation depth of the uniform pooling model is also presented (dashed red). The parameters used here are: , κ VPM = 1, γ = 1, D = 10hz, d = 3ms, ϕ 0 = 5π/6. For the STDP we used the temporally exponential asymmetric kernel, Eq (9) , with τ − = 50ms, τ + = 22ms, μ = 0.01, α = 1.1, and λ = 0.01.

Fig 5a depicts the distribution of preferred phases as a function of the distribution width of their thalamic input, κ VPM . For a uniform input distribution, κ VPM = 0, the downstream distribution is also uniform, κ L4I = 0, see [ 22 ]. As the distribution in the VPM becomes narrower, so does the distribution in the L4I population. If the distribution of the thalamic population is narrower than a certain critical value, STDP will converge to a fixed point and κ L4I will diverge. Typically, we find that the width of the cortical distribution, κ L4I , is similar to or larger than that of the upstream distribution, κ VPM , Fig 5b . This sharpening is obtained via STDP by selectively amplifying certain phases while attenuating others. Consequently the rhythmic component, in terms of the modulation depth, γ L4I , is also typically amplified relative to the uniform pooling model, Fig 5c .

For small N, N ⪅ 70, STDP induces a point measure distribution over time, Fig 4a . This is due to pinning in the noiseless STDP dynamics in the mean field limit, λ → 0. Stochastic dynamics, due to noisy neuronal responses, can overcome pinning. Fig 4c depicts the distribution of the preferred phases for different values of N, which results from both STDP dynamics and quenched averaging over different realizations of the preferred phases of the thalamic neurons. For small values of N, the distribution is dominated by the quenched statistics, in terms of its narrow width and mean that is dominated by the mean phase of the thalamic population. As N increases, the distribution widens and is centered around a preferred phase that is determined by the STDP. Thus, activity dependent plasticity helps to shape the distribution of preferred phases in the downstream population. Below, we examine how these different parameters affect the ways in which STDP shapes the distribution; hence, we will not average over the quenched disorder.

The distribution of the preferred phases of L4I neurons, κ L4I , is shown as a function of N for different sources of variability. In each column (value of N), the non-normalized distribution: Pr(ψ)/max{Pr(ψ)}, is presented by color. (a) The temporal distribution of the preferred phase of a single L4I neuron as a result of STDP dynamics without quenched disorder or averaging (see Modeling pre-synaptic phase distributions) is presented. (b) Distribution due to quenched disorder without STDP. The distribution of L4I neurons in the randomly pooling model was estimated from 1000 trials of drawing the preferred phases of N VPM neurons in an i.i.d. manner from Eq (1) . (c) The distribution due to both quenched disorder and STDP dynamics is shown. Unless stated otherwise, the parameters used in these graphs are as follows: , κ VPM = 1, γ = 0.9, D = 10hz, d = 3ms, ϕ 0 = 5π/6. For the STDP we used the temporally exponential asymmetric kernel, Eq (9) , with τ − = 50ms, τ + = 22ms, μ = 0.01, α = 1.1, and λ = 0.01.

Fig 4a shows the (non-normalized) distribution of preferred phases induced by STDP as a function of the number of pooled thalamic neurons, N. For large N, the STDP dynamics converge to the continuum limit of Eq (11) , and the distribution converges to a limit that is independent of N. This contrasts with the naive pooling model lacking plasticity (cf. Fig 4b ). In that model, the distribution of preferred phases in the cortical population results from a random process of pooling phases from the thalamic population. We shall refer to this type of variability as quenched disorder, since this randomness is frozen and does not fluctuate over time. This process is characterized by a distribution of preferred phases with vanishing width, κ L4I → ∞, in the large N limit, Fig 4b . In addition to the different widths of the distribution, the mean of the preferred phases also differs greatly. Whereas in Fig 4a the mean preferred phase is determined by the STDP rule, without STDP, Fig 4b , the mean phases of L4I neurons is given by the mean preferred phase of the VPM population shifted by the delay, κ L4I = ψ VPM + dν.

The continuous remodelling of the synaptic weights causes the preferred phase of the downstream neuron, κ L4I (see Eq (4) ) to drift in time, Fig 3b . As can be seen from the figure, the drift velocity is not constant. Consequently, the downstream (L4I) neuron ‘spends’ more time in certain phases than others. Thus, the STDP dynamics induce a distribution over time for the preferred phases of the downstream neuron. One can estimate the distribution of the preferred phases of L4I neurons by tracking the phase of a single neuron over time. Alternatively, since our model is purely feed-forward, the preferred phases of different L4I neurons are independent; hence, this distribution can also be estimated by sampling the preferred phases of different L4I neurons at the same time. Fig 3d and 3e show the distribution of preferred phases of thalamic and L4I whisking neurons, respectively. Thus, STDP induces a distribution of preferred phases in the L4I population, which is linked to the temporal distribution of single L4I neurons.

(a) Synaptic weight dynamics. Each trace depicts the time evolution of a single synaptic weight, differentiated by color according to its preferred phase, see legend. (b) & (c) Dynamics of the order parameters. The preferred phase of the downstream neuron, κ L4I , (in (b)), and the mean, , and the magnitude of the first Fourier component, , (in red and black, respectively, in (c)) are shown as a function of time. (d) The distribution of preferred phases in the thalamic population that served as input to the downstream L4I neuron is presented as a polar histogram. The distribution followed Eq (1) with κ VPM = 1 and ψ VPM = 5π/6 rad. (e) The temporal distribution of preferred phases of the downstream L4I neuron is shown in a polar plot. Fitting the von-Mises distribution yielded κ L4I = 1.1 and κ L4I = 0.8 rad. In this simulation we used the following parameters: N = 150, , γ VPM = 1, D VPM = 10hz, and d = 3ms. The temporally asymmetric STDP rule, Eq (9) , was used with τ − = 50ms, τ + = 22ms, μ = 0.01, α = 1.1, and λ = 0.01.

Fig 3a shows the temporal evolution of the synaptic weights (color-coded by their preferred phases). As can be seen from the figure, the synaptic weights do not relax to a fixed point; instead there is a continuous remodelling of the entire synaptic population. Examining the order parameters, Fig 3b and 3c , reveals that the STDP dynamics converge to a limit cycle.

In the limit of a slow learning rate, λ → 0, we obtain deterministic dynamics for the mean synaptic weights (see [ 32 ] for a detailed derivation) (11) with (12) where Γ j, L4I (Δ) is the cross-correlation function between the jth thalamic pre-synaptic neuron and the L4I post-synaptic neuron; see detailed derivation in Temporal correlations.

For the temporally symmetric learning rule we used a difference of Gaussians model, (10) where τ ± are the temporal widths. In this case, the order of firing is not important, only the absolute time difference.

Empirical studies portray a wide variety of temporal kernels [ 34 – 42 ]. Specifically, in our work, we used two families of STDP rules: 1. A temporally asymmetric kernel [ 35 , 37 – 39 ]; 2. A temporally symmetric kernel [ 38 , 40 – 42 ]. Both of these rules have been observed in the the barrel system of mice, at least for some developmental period [ 43 – 45 ]. For the temporally asymmetric kernel we used the exponential model, (9) where Θ(x) is the Heaviside function, and τ ± are the characteristic timescales of the potentiation (+) or depression (−). We take τ − > τ + as typically reported.

We model the modification of the synaptic weight, Δw, following either a pre- or post-synaptic spike as the sum of two processes: potentiation (+) and depression (-) [ 31 – 33 ], as (7) The parameter λ denotes the learning rate. We further assume a separation of variables and write each term (potentiation and depression) as the product of the function of the synaptic weight, f ± (w), and the temporal kernel of the STDP rule, K ± (Δt), where Δt = t post − t pre is the time difference between pre- and post-synaptic spike times. Following Gütig et al. [ 33 ] the weight dependence functions, f ± (w), were chosen to be: (8a) (8b) where α > 0 is the relative strength of depression and μ ∈ [0, 1] controls the non-linearity of the learning rule.

Due to the rhythmic nature of the thalamic inputs, Eq (2) , and the linearity of the L4I neuron, Eq (3) , the L4I neuron will exhibit rhythmic activity: (4) with a mean, D L4I , a modulation depth, γ L4I , and a preferred phase ψ L4I , that depend on global order parameters that characterize the thalamocortical synaptic weights population. For large N these order parameters are given by: (5) and (6) where is the mean synaptic weight and is its first Fourier component. The phase ψ is determined by the condition that is real and non-negative. Consequently, the L4I neurons in our model respond to whisking with a mean , a modulation depth of , and a preferred phase ψ L4I = ψ + νd.

To facilitate the analysis, we model the response of the downstream layer 4 inhibitory (L4I) neuron to its thalamic inputs by the linear Poisson model, which has been used frequently in the past [ 22 , 24 – 31 ]. Given the thalamic responses, the firing of the L4I neuron follows inhomogeneous Poisson process statistics with an instantaneous firing rate (3) where d > 0 represents the characteristic delay, and w k is the synaptic weight of the kth thalamic neuron.

We model a system of N thalamic excitatory whisking neurons, synapsing in a feed-forward manner onto a single inhibitory L4 barrel cortical neuron. Unless stated otherwise, following [ 16 ], in our numerical simulations we used N = 150. The spiking activity of the thalamic neurons is modelled by independent inhomogeneous Poisson processes with an instantaneous firing rate that follows the whisking cycle: (2) where ρ k (t) = ∑ i δ(t − t k,i ), k ∈ {1, …N}, is the spike train of the kth thalamic neuron, with denoting its spike times. The parameter D is the mean firing rate during whisking (averaged over one cycle), γ is the modulation depth, ν is the angular frequency of the whisking, and ϕ k is the preferred phase of firing of the kth thalamic neuron. We further assume that the preferred phases in the thalamic population, , are distributed i.i.d. according to the von-Mises distribution, Eq (1) .

Discussion

We examined the possible contribution of STDP to the diversity of preferred phases of whisking neurons. Whisking neurons can be found along different stations throughout the somatosensory information processing pathway [1, 4–6, 9, 10, 13, 16, 47]. Here we focused on L4I whisking neurons that receive their whisking input mainly via excitatory inputs from the VPM, and suggested that the non-trivial distribution of preferred phases of L4I neurons results from a continuous process of STDP.

STDP has been reported in thalamocortical connections in the barrel system. However, STDP has only been observed during the early stages of development [43–45]. Is it possible that the contribution of STDP to shaping thalamocortical connectivity is restricted to the early developmental stages? This is an empirical question that can only be answered experimentally. Nevertheless, several comments can be made.

First, Inglebert and colleagues recently showed that pairing single pre- and post-synaptic spikes, under normal physiological conditions, does not induce plastic changes [48]. On the other hand, activity dependent plasticity was observed when stronger activity was induced. Thus, whisking activity, which is a strong collective and synchronized activity, may induce STDP of thalamocortical connectivity.

Second, in light of the considerable volatility of synaptic connections observed in the central nervous system [17–20], it is hard to imagine that thalamocortical connectivity remains unchanged throughout the lifetime of an animal. Third, if activity independent plasticity alone underlies synaptic volatility, thalamocortical synaptic weights would be expected to be random. In this case, thalamic whisking input to layer 4 should be characterized by an extremely narrow distribution centered around the delayed mean thalamic preferred phase, cf. Fig 4b. Consequently, neglecting the possible contributions due to recurrent interactions with layer 4, a non-trivial distribution of L4I phases, with κ L4I ∼ 1, can be obtained either via STDP or by pooling the whisking signal from an extremely small VPM population of N < 10 neurons. However, in the latter scenario (without STDP) the mean preferred phase of L4I neurons is expected to be determined by the (delayed) mean phase of VPM neurons (see e.g. Fig 4b), thus raising serious doubts as to the viability of the latter solution.

An alternative hypothesis to our proposed mechanism is that the recurrent dynamics within layer 4 generate the empirically observed wide distribution of preferred phases centered on κ L4I ≈ ψ VPM − 2.2 [rad]. As layer 4 excitatory neurons have been reported to exhibit very weak rhythmic activity [13, 16], such a mechanism would mainly be based on recurrent interaction within the L4I population. The putative ways in which inhibitory interactions widen the distribution remain unclear. Nevertheless, this hypothesis should be explored further.

Our hypothesis views STDP as a continuous process. As a result of STDP, the preferred phase of a L4I neuron, in our model, is a stochastic process with a non-uniform drift, Fig 3b, which, in turn, induces a distribution over time for the phase, Fig 3e. As our model lacks recurrent interactions, phases of different L4I neurons will be i.i.d. random processes, albeit with different random initial conditions. Consequently, simultaneous sampling of the distribution of preferred phases of a population of L4I neurons will yield the same distribution as obtained by sampling the preferred phase of a single neuron over time. Thus, functionality, in terms of transmission of the whisking signal and retaining the stationary distribution of the preferred phases in the downstream population is obtained as a result of continuous remodelling of the entire population of synaptic weights.

Several key features of our hypothesis provide clear empirical predictions. Our theory links two empirical observations. One is the distribution of preferred phases in up- and downstream populations. The other is the STDP rule. The literature can provide a good assessment of the distribution of preferred phases of thalamic and of cortical whisking neurons. Further effort is required to estimate the STDP rule of the thalamocortical connections to L4I neurons. Thus, we can reject our hypothesis if the predicted distribution, based on our theory and the estimated learning rule, is not consistent with the empirically observed distribution.

Our theory also posits that the distribution of preferred phases results from their stochastic dynamics. Our theory also incorporates a direct relationship between the drift velocity and the distribution of preferred phases. Essentially, the system will spend more time around phases characterized by slower velocity; hence, a higher probability of preferred phases should correlate with slower drift velocity. To the best of our knowledge, the dynamics of the preferred phases of whisking neurons has not been studied. Monitoring the preferred phases of single L4I neurons over long periods of time (days and weeks) could serve as additional tests of our theory. Specifically, it would address the question of whether the preferred phases drift in time, and if so, whether the drift velocity and the distribution of preferred phases are consistent with our theory.

Further empirical exploration can test the degree in which our approximation of a purely feed-forward architecture (i.e., neglecting recurrent interactions) holds. Simultaneous recordings of a large population of fast spiking layer 4 whisking neurons can be harnessed to compare the distribution of preferred phases over space and time, which are expected to be identical in the absence of recurrent interactions. In addition, in the visual system, Lien and Scanziani combined intracellular recordings with optogenetic methods to assess the contribution of thalamic inputs to the selectivity of cortical neurons [49]. In the context of the current work, we explored the thalamic inputs to whisking neurons in the barrel cortex. Thus, applying ideas and approaches from the work of Lien and Scanziani could (i) provide a measurement of the thalamic input, and (ii) determine the relative contribution of feed-forward and recurrent connections to the selectivity to the phase of whisking in L4I neurons.

More broadly, our findings raise the question of whether the constant remodelling of synaptic efficacies is an artefact of the ‘biological hardware’ the brain must overcome, or whether it reflects an underlying principle of the central nervous system. Though we cannot provide an answer to this question at this point in time, we may speculate on the possible advantages of a dynamic solution. One interesting hypothesis to consider is that the constant remodelling makes the neuronal network more flexible and consequently enables a faster response to a changing environment.

In the current study, we made several simplifying assumptions. The spiking activity of the thalamic population was modeled as a rhythmic activity with a well-defined frequency. However, whisking activity spans a frequency band of several Hertz. Moreover, the thalamic input relays additional signals, such as touch and texture. These signals modify the cross-correlation structure and add ‘noise’ to the dynamics of the preferred phase of the downstream neuron. As a result, the distribution of preferred phases in the downstream population is likely to be wider. In addition, our analysis used a purely feed-forward architecture ignoring recurrent connections in layer 4, which may also affect the preferred phases in layer 4. A quantitative application of our theory to the whisker system should consider all these effects. Nevertheless, the theory presented provides the basic foundation to address these effects.

[END]

[1] Url: https://journals.plos.org/ploscompbiol/article?id=10.1371/journal.pcbi.1009353

(C) Plos One. "Accelerating the publication of peer-reviewed science."
Licensed under Creative Commons Attribution (CC BY 4.0)
URL: https://creativecommons.org/licenses/by/4.0/


via Magical.Fish Gopher News Feeds:
gopher://magical.fish/1/feeds/news/plosone/

You are viewing proxied material from magical.fish. The copyright of proxied material belongs to its original authors. Any comments or complaints in relation to proxied material should be directed to the original authors of the content concerned. Please see the disclaimer for more details.