# Models of Metaplasticity: A Review of Concepts

Pierre Yger, Matthieu Gilson

► To cite this version:

Pierre Yger, Matthieu Gilson. Models of Metaplasticity: A Review of Concepts. Frontiers in Computational Neuroscience, 2015, 9, pp.138. 10.3389/fncom.2015.00138. hal-01263771

HAL Id: hal-01263771

https://hal.sorbonne-universite.fr/hal-01263771

Submitted on 28 Jan 2016

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

CrossMark

# Models of Metaplasticity: A Review of Concepts

Pierre Yger¹,²* and Matthieu Gilson³

¹ Sorbonne Université, UPMC Univ Paris06 UMRS968, Paris, France, ² Institut de la Vision, INSERM, U968, Centre National de la Recherche Scientifique, UMR7210, Paris, France, ³ Computational Neurosciences Group, Departament de Tecnologies de la Informació i les Comunicacions, Universitat Pompeu Fabra, Barcelona, Spain

Part of hippocampal and cortical plasticity is characterized by synaptic modifications that depend on the joint activity of the pre- and post-synaptic neurons. To which extent those changes are determined by the exact timing and the average firing rates is still a matter of debate; this may vary from brain area to brain area, as well as across neuron types. However, it has been robustly observed both *in vitro* and *in vivo* that plasticity itself slowly adapts as a function of the dynamical context, a phenomena commonly referred to as metaplasticity. An alternative concept considers the regulation of groups of synapses with an objective at the neuronal level, for example, maintaining a given average firing rate. In that case, the change in the strength of a particular synapse of the group (e.g., due to Hebbian learning) affects others' strengths, which has been coined as heterosynaptic plasticity. Classically, Hebbian synaptic plasticity is paired in neuron network models with such mechanisms in order to stabilize the activity and/or the weight structure. Here, we present an oriented review that brings together various concepts from heterosynaptic plasticity to metaplasticity, and show how they interact with Hebbian-type learning. We focus on approaches that are nowadays used to incorporate those mechanisms to state-of-the-art models of spiking plasticity inspired by experimental observations in the hippocampus and cortex. Making the point that metaplasticity is an ubiquitous mechanism acting on top of classical Hebbian learning and promoting the stability of neural function over multiple timescales, we stress the need for incorporating it as a key element in the framework of plasticity models. Bridging theoretical and experimental results suggests a more functional role for metaplasticity mechanisms than simply stabilizing neural activity.

**Keywords:** synaptic plasticity, metaplasticity, Hebbian learning, homeostasis, STDP

## 1. INTRODUCTION

The brain is made of billions of neurons able to efficiently process the huge flow of information impinging continuously on sensory modalities, extracting relevant data, and producing appropriately timed responses. Even during development (Corlew et al., 2007; Wang et al., 2012) or when lesioned (Young et al., 2007; Beck and Yaari, 2008), the brain has the striking capability to adapt in order to maintain the stability of neural functions. Importantly, this slow adaptation, acting at a timescale of hours or days (Turrigiano and Nelson, 2000; Davis, 2006) is performed in conjunction with fast changes often observed in the so called Hebbian learning (Hebb, 1949). Understanding the mechanisms leading to the dynamical organization of neuronal network via the fine interactions of those two competing processes is therefore a crucial step toward analyzing the stability of the computations performed by cerebral activity.

**OPEN ACCESS**

**Edited by:**

Friedemann Zenke,
Stanford University, USA

**Reviewed by:**

Harel Z. Shouval,
University of Texas Medical School at
Houston, USA
Christian Tetzlaff,
Max Planck Institute for Dynamics and
Self-Organization, Germany

*Correspondence:**

Pierre Yger
pierre.yger@inserm.fr

**Received:** 02 July 2015

**Accepted:** 27 October 2015

**Published:** 10 November 2015

**Citation:**

Yger P and Gilson M (2015) Models of Metaplasticity: A Review of Concepts. *Front. Comput. Neurosci.* 9:138. doi: 10.3389/fncom.2015.00138
Following the seminal idea that neurons firing together should wire together (Hebb, 1949), numerous experimental studies have been conducted to unravel part of the links between plasticity and neuronal activity. Nowadays, this so-called Hebbian form of plasticity in the brain has been characterized experimentally in many areas, involving multiple but still misunderstood molecular pathways (see Abbott and Nelson, 2000; Caporale and Dan, 2008, for reviews). While it is commonly assumed that NMDA receptors are the primary actors in long-term potentiation, or LTP (Feldman, 2012), the biochemical pathways for long-term depression (LTD) seem to differ in cortex and in hippocampus (Wang et al., 2005; Bender et al., 2006; Nevian and Sakmann, 2006). In controlled *in vitro* experiments, it has also been shown that LTP and LTD depend on the precise timing of pre- and post-synaptic spikes (Markram et al., 1997; Bi and Poo, 1998), leading to the concept of timing-LTP/LTD or spike-timing-dependent plasticity (STDP).

By acting independently at each synapse without spatial or temporal crosstalk among synapses, Hebbian learning is a form of homosynaptic plasticity that is intrinsically unstable. In point of fact, provided synapses are reinforced when both the pre- and post-synaptic neurons are active, nothing prevents the synapses from strengthening themselves boundlessly, which causes the post-synaptic activity to explode (Rochester et al., 1956; von der Malsburg, 1973; Miller, 1996). While this instability can be avoided by artificially imposing hard boundaries onto the synaptic weights, several learning models came with intrinsic mechanisms regulating the synaptic efficacies (Bienenstock et al., 1982; Oja, 1982) in order to solve this issue in a less fine-tuned manner.

The present paper reviews such mechanisms that aim to tame the positive feedback provided by Hebbian plasticity. In the biology, some homeostatic mechanisms can be viewed as independent from the Hebbian learning that they counterbalance. For example, the sum of synaptic strengths may be up or down regulated to maintain the average post-synaptic firing rate; see Vitureira and Goda (2013) for a review of the biophysics of such mechanisms. In contrast, other processes directly modulate the learning rule itself as a function of the dynamical context, which is referred to as metaplasticity. This concept is the plasticity of the synaptic plasticity itself (Abraham and Bear, 1996; Abraham, 2008), and it is tightly related to the notion of homeostasis (O'Leary and Wyllie, 2011). To ensure the overall stability of the neuronal system, a key role for metaplasticity is to regulate the synaptic update rules in terms of the past history of the activity at the whole neuronal level. Many experiments have demonstrated metaplasticity using distinct protocols (Abraham, 2008). Quite often, it also involves some form of heterosynaptic plasticity, in the sense that the local changes affecting a particular synapse onto a post-synaptic neuron influence the plasticity for neighboring synapses.

The study of the dynamical implications of the interaction between homeostatic mechanisms and Hebbian plasticity requires the integration of experimental data in model studies (Marder and Goaillard, 2006). From a modeler's point of view, interactions between Hebbian learning and its regulating counterpart, either by homeostatic mechanisms or

by metaplasticity, is problematic. The principal reason being that those two distinct forms of plasticity do not act on similar timescales. Following experimental results, it is commonly assumed that synaptic changes triggered by Hebbian plasticity protocols are rather fast (Bliss and Lomo, 1973; Sjöström et al. 2001, 2003; Wang et al., 2005), occurring in the timescale of minutes or faster, while metaplasticity or homeostatic changes are much slower (Abraham and Bear, 1996), in the order of days. The present paper provides a theoretical framework to analyze the interaction between Hebbian and homeostatic plasticities at different timescales. In this way it gives an overarching view of different methods used in the literature to solve the above-mentioned instability issue of Hebbian plasticity. Maintaining the stability only being one of the requirements for proper behavior, we will discuss how homeostatic constraints can also be used to adjust the function implemented by the neural circuits.

## 2. THE APPARENT ANTAGONISM BETWEEN HEBBIAN AND HOMEOSTATIC PLASTICITY

### 2.1. Two Divergent Goals

As it has already been observed (Turrigiano and Nelson, 2000; Watt and Desai, 2010; Vitureira and Goda, 2013), Hebbian and homeostatic plasticities are two apparently opposing processes, which compete at the synaptic level to fulfill different goals. Hebbian learning promotes strong or synchronous firing among neurons, which is hypothesized to be a building block for memory storage (Nabavi et al., 2014). In contrast, homeostatic processes counterbalance such intense spiking activity to maintain the global stability in neuronal networks (Turrigiano and Nelson, 2000; Turrigiano, 2008; Pozo and Goda, 2010). Several types of homeostatic processes have been observed at the neuronal level in many brain areas, such as synaptic scaling (Turrigiano et al., 1998) and intrinsic plasticity (Zhang and Linden, 2003).

It has been long known that Hebbian plasticity *alone* is intrinsically unstable (Rochester et al., 1956; von der Malsburg, 1973; Miller, 1996). The entrainment between synapses often force all to grow boundlessly or to a maximal set value; in other cases, they may all become silent. To circumvent these issues of traditional rate-based Hebbian learning, weight normalization can be introduced to prevent the runaway of synapses (Oja, 1982; Miller, 1996). In the context of spiking activity, STDP has been termed “temporally Hebbian” when it promotes synchronous firing. Weight-dependent STDP update rules, which induces more LTD than LTP for strong synapses, provide a fixed point in the learning dynamics (van Rossum et al., 2000; Gütig et al., 2003). Although this ensures some stability, it may dramatically change the weight distribution from being bimodal to being unimodal. In the case of a narrow unimodal weight distribution, competition induced by STDP among synapses is weakened between pathways with distinct characteristics (e.g., rate, correlation), which is not functionally interesting. For weight-dependent STDP, this trade-off compromise is only fulfilled in a given parameter range. In recurrent networks
especially, the synaptic specialization by competition may be severely impaired without fine tuning (Morrison et al., 2007; Gilson and Fukai, 2011).

## 2.2. Two Different Timescales

Most of the plasticity protocols performed *in vitro* are based on either input stimulation at a high/low frequency leading to LTP/LTD (Bliss and Lomo, 1973) or STDP-type pairings of pre-post spikes (Markram et al., 1997; Bi and Poo, 1998; Sjöström et al., 2001; Froemke and Dan, 2002; Wang et al., 2005). The typical protocol used in cortical or hippocampal slices to elicit STDP *in vitro* using spike pairs is represented in **Figure 1A**: a spike is triggered at the pre-synaptic neuron and another at the post-synaptic neuron with time difference $\delta t = t_{\text{pre}} - t_{\text{post}}$. This pairing is repeated approximately 60 times with frequency $f_{\text{pairing}} = 1$ Hz in order to see a robust change in the weight: it has been shown that after an induction phase, the total weight change evolves non-linearly up to a saturation plateau, at around 60–100 pairings (Froemke et al., 2006), which corresponds to the number of protocol repetitions in most studies.

For modelers, this STDP protocol leads to the simplified view of the time-difference window in **Figure 1B**, where a single pre spike followed by a post spike will trigger LTP, whereas post followed by pre causes LTD. This is clearly an over-simplification of a much more complex phenomenon. Just to mention some limitations of this simplified view, it has been shown that if the frequency $f_{\text{pairing}}$ of the pairing is changed, the typical STDP curve with LTP for $\delta t < 0$ and LTD for $\delta t > 0$ is dramatically modified (Sjöström et al., 2001). Depression is only visible for low

frequency pairings, when pairings are performed with $\delta t < 0$ and $f_{\text{pairing}} < 20$ Hz. For $f_{\text{pairing}} > 20$ Hz, however, synapses undergo LTP irrespective of the sign for $\delta t$. Moreover, several *in vitro* studies on cortical pyramidal neurons showed that the canonical shape of the STDP curve for such pre-post pairings strongly depends on the position of the synapse along the dendritic tree (Froemke et al., 2005; Letzkus et al., 2006; Kampa et al., 2007), as well as the post-synaptic voltage (Artola et al., 1990). Those experimental findings led to the refinements of initial STDP models based on the curve, in order to incorporate the observed effects for triplets of spikes, spike bursts, clamping the post-synaptic membrane potential and so on (Pfister and Gerstner, 2006; Clopath et al., 2010; El Boustani et al., 2012; Graupner and Brunel, 2012; Yger and Harris, 2013).

Despite those efforts, there is a point that is almost never considered: STDP changes are not instantaneous. In most experiments, when plasticity protocols are performed, the resulting weight is recorded up to 30 min later. The curve in **Figure 1B** corresponds to the corresponding weight change divided by the number of pairings. In models of classical (van Rossum et al., 2000; Song and Abbott, 2001) and weight-dependent (van Rossum et al., 2000; Gütig et al., 2003; Morrison et al., 2007; Gilson and Fukai, 2011) STDP, its final value is the results of additive instantaneous and independent weight updates following each pairing. In fact, even elaborate models consider the linear summation of weight updates, even when contributions are restricted to neighboring spikes (Burkitt et al., 2004). Only a few attempts have been done to change this property that is convenient for theory, such as probabilistic

<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | Intrinsic timescale of Hebbian learning. (A) The classical STDP pairing protocols widely used in the literature. (B) Synaptic modification for one pair of pre- and post-synaptic spikes, as a function of their relative timing. (C) Evolution as a function of time of a single synaptic weight, after an STDP protocol, for various papers taken from the literature, both for LTP of LTD protocols [dash-dotted thin black line is the null-line for Sjöström et al. (2001)]. (D) Adapted from Keck et al. (2013), Normalized mEPSC amplitude in a layer 5 cell in the mice visual cortex following a lesion in the retina. (E) Adapted from Huang et al. (1992), Prior synaptic activity triggered during the red shaded area (LTP priming, red curve) reduces LTP in CA1 hippocampus compared to control without pre-activation (black curve). (F) Adapted from Mockett et al. (2002), Low frequency stimulation (LFS, red shaded areas) influences non-linearly the amount of LTD in CA1 hippocampus: black curve, (control with only one LFS), red curve (two consecutive LFS).</figcaption></figure>
models of STDP (Appleby and Elliott, 2005). By re-examining the weight traces found in the STDP literature (Bi and Poo, 1998; Sjöström et al., 2001; Froemke and Dan, 2002; Froemke et al., 2006) and reproduced in Figure 1C, it can be seen that the weights actually evolve continuously in vitro. Therefore, plasticity should better be seen as a phenomenon that is triggered by a stimulation event and evolves toward a new equilibrium with a time constant $τ_{Hebb} \simeq 10$ min.

Now considering that Hebbian plasticity induces such a transient synaptic change, the question arises about its interaction with homeostatic plasticity. Those processes, either intrinsic or synaptic, are assumed to be much slower. For example synaptic scaling, one of the numerous mechanisms of homeostasis, takes place in vitro with a time constant $τ_{homeo}$ of the order of a day (Turrigiano and Nelson, 2000), and in vivo during the 2-3 days after an abrupt change, as observed for neurons in the visual cortex following visual deprivation (Hengen et al., 2013; Keck et al., 2013). Figure 1D, adapted from Keck et al. (2013), shows the amplitude of miniature EPSC in V1 neurons after a bilateral lesion in the adult retina: after an initial period of about a day, amplitudes are scaled up to compensate for the reduced inputs. Together, these results stress the fact that Hebbian and homeostatic processes have distinct timescales. Understanding the biological mechanisms responsible for those changes at the molecular level is necessary to gain a better insight on the interaction between them, especially in vivo where synapses are constantly bombarded by spikes.

## 2.3. Primings as an Evidence for Metaplasticity

Although on a first approximation it may appear that $τ_{homeo} \gg τ_{Hebb}$, several experiments show that those two timescales may be more interleaved. In hippocampal slices, it has been shown in so-called priming experiments that the activation of a synapse before its reactivation modulates the plasticity triggered later at that particular synapse (see Figures 1E,F). In Figure 1E that is adapted from Huang et al. (1992), weak tetanic priming stimulations can reduce the amount of LTP obtained during a strong subsequent tetanic stimulation; note that the effects last more than 1 h. On the contrary, the LTD pathway seems to be facilitated when the synapse is preactivated a few hours before the plasticity protocol (Christie and Abraham, 1992; Wang, 1998; Mockett et al., 2002). This is illustrated in Figure 1F, adapted from Mockett et al. (2002), where the effect lasts at least 2 h. Those primings experiments suggest the existence of long-lasting regulation mechanisms, acting over large time constants, which counteracts the effect of Hebbian learning. This modulation of the Hebbian plasticity by preactivation of the synaptic pathway is a direct application of the so-called metaplasticity (Abraham and Bear, 1996), i.e., the plasticity of the learning rules themselves.

## 3. MATHEMATICAL FORMALISM

To formally study the interactions between Hebbian and homeostatic plasticity, we use the following mathematical formalism. We consider a Poisson neuron (Kempter et al., 1999) with N synapses indexed by i, corresponding to the input firing

rates $r_i$; for STDP examples, we also define the input cross-covariances $c_{ij}$ between neurons i and j. The equations for the output firing rate $r_{post}$ and pre-post covariances $c_{i-post}$ between synapse i and the post-synaptic spike train in a feedforward scenario are given by

$$ 
\begin{aligned}
r_{\text{post}} &= \sum_{1 \le j \le N} w_j r_j \\
c_{i-\text{post}} &= \sum_{1 \le j \le N} w_j c_{ij}
\end{aligned} 
\quad (1) $$

In order to compare several learning rules in the context of metaplasticity, we consider the following general equations for the evolution of a given weight $w_i$ and a modulation parameter $θ$:

$$ 
\begin{aligned}
\dot{w}_i &= \frac{1}{\tau_{\text{Hebb}}} \Phi(w_i, r_i, r_{\text{post}}, c_{i-\text{post}}, \theta) \\
\dot{\theta} &= \frac{1}{\tau_{\text{homeo}}} [\Psi(r_{\text{post}}) - \theta]
\end{aligned} 
\quad (2) $$

The motivation for these expressions is to model the two timescales explicitly, as previously done for the BCM rule (Bienenstock et al., 1982) and for an extension of the triplet STDP rule (Zenke et al., 2013): $τ_{Hebb}$ and $τ_{homeo}$ are the two time constants at which both Hebbian and homeostatic changes are propagated onto the synapses. The Hebbian plasticity update is embodied in $Φ$, which also depends on $r_{pre}$, $c_{pre-post}$, etc. The parameter $θ$ is global for all synapses of a neuron and interacts or modulates the corresponding weight updates. Typically, it is used to implement a homeostatic mechanism, as we will see for several models of synaptic plasticity that are commonly used in the literature. The present framework could be extended to incorporate other non-linearities in the firing mechanism (e.g., LIF neuron), adaptation or intrinsic plasticity.

### 3.1. Stability Analysis for the Mean-field Dynamical System

Ignoring correlations and inhomogeneities across synapses, we focus on the analysis of the mean weight $\bar{w} = \sum_j w_j/N > 0$ with mean input rate $r_{pre}$. The rate Equation (1) simply becomes

$$ r_{\text{post}} = \bar{w} r_{\text{pre}} \quad (3) $$

This allows for an easy comparison of the weight dynamics based on polynomial expressions in $\bar{w}$. Other neuron models usually give more complex mapping between input and output rate/correlations, but the common trend is that they are monotonically increasing function of the weight $\bar{w}$. This property is the cause for the instability of Hebbian learning, as it increases $\bar{w}$ all the more as $r_{post}$ is large. Therefore, we will review through the example of the Poisson neuron how stabilization mechanisms interact with the Hebbian component.

In order to examine the stability of the mean-field dynamical system (Equation 2) where $w_i$ is replaced by $\bar{w}$, we consider its Jacobian matrix.

$$ \left( \begin{array}{c} \frac{1}{\tau_{\text{Hebb}}} \left[ \frac{\partial \Phi}{\partial w} + \frac{\partial \Phi}{\partial r_{\text{post}}} r_{\text{pre}} \right] \frac{1}{\tau_{\text{Hebb}}} \frac{\partial \Phi}{\partial \theta} \\ \frac{1}{\tau_{\text{homeo}}} \frac{\partial \Psi}{\partial r_{\text{post}}} r_{\text{pre}} \end{array} \right) = \left( \begin{array}{c} a \\ c \end{array} \right) \quad (4) $$
For the top-left term in the Jacobian, we have used the following equality for the feedforward architecture corresponding to Equation (3): $\frac{\partial r_{post}}{\partial w} = r_{pre}$. The eigenvalues of the Jacobian matrix are given by

$$x_{\pm} = \frac{1}{2}(T \pm \sqrt{T^2 - 4D}) \quad (5)$$

where $T = a + d$ is the trace and $D = ad - bc$ is the determinant. To ensure stability for this 2-dimensional dynamical system, these eigenvalues must be real negative. This requires that the following relationships are satisfied.

$$\begin{array}{l} T < 0 \\ 0 < 4D < T^2 \end{array} \quad (6)$$

If, however, the discriminant is positive with the trace still negative ($T < 0$ and $T^2 < 4D$) the system exhibits damped oscillations related to the imaginary eigenvalues. With purely imaginary eigenvalues, we may obtain a limit cycle. Finally, when $D < 0$ or $T > 0$, at least one eigenvalue is positive and can lead to an explosion of the mean weight.

## 3.2. Competition between Input Pathways

Following Kempter et al. (1999); Gütig et al. (2003); Gilson et al. (2009), we can use and rewrite Equation (1) to study the competition between learning weights.

$$r_{post} = w_1 r_1 + w_2 r_2 \quad (7)$$

Again ignoring correlations, we obtain the following 3-dimensional learning system.

$$\begin{align} \dot{w}_1 &= \frac{1}{\tau_{\text{Hebb}}} \Phi(w_1, r_1, r_{\text{post}}, \theta) \notag \\ \dot{w}_2 &= \frac{1}{\tau_{\text{Hebb}}} \Phi(w_2, r_2, r_{\text{post}}, \theta) \tag{8} \\ \dot{\theta} &= \frac{1}{\tau_{\text{homeo}}} [\Psi(r_{\text{post}}) - \theta] \notag \end{align}$$

Considering the equilibrium for the mean weight $\bar{w} = (w_1 + w_2)/2$ to be satisfied, the competition between the two input pathways can be studied for what is called “symmetry breaking,” namely the divergence of $w_1$ and $w_2$. This relates to the following differential equation for the weight difference $\Delta w = w_1 - w_2$, which quantifies the tendency for splitting

$$\begin{align} \dot{\Delta w} &= \frac{1}{\tau_{\text{Hebb}}} \left[ \Phi(w_1, r_1, r_{\text{post}}, \theta) - \Phi(w_2, r_2, r_{\text{post}}, \theta) \right] \quad (9) \\ &\approx \frac{1}{\tau_{\text{Hebb}}} \frac{\partial \Phi}{\partial w} (\bar{w}, \bar{r}_{\text{pre}}, r_{\text{post}}, \theta) \Delta w \notag \\ &\quad + \frac{1}{\tau_{\text{Hebb}}} \frac{\partial \Phi}{\partial r_{\text{pre}}} (\bar{w}, \bar{r}_{\text{pre}}, r_{\text{post}}, \theta) \Delta r \notag \end{align}$$

where $\bar{r}_{pre} = (r_1 + r_2)/2$ and $\Delta r = r_1 - r_2$ is assumed to be small here. The larger positive $\frac{\partial \Phi}{\partial w}(\bar{w}, \bar{r}_{pre}, r_{post}, \theta)$ is, the more strongly the weights $w_1$ and $w_2$ will move apart from each other.

## 3.3. Conditions for Joint Stability and Competition for Hebbian Learning with Synaptic Scaling

In general, the equations for stability and competition may turn out to be quite complex, even for the mean-field dynamical system. The ambition here is to describe the general trends for the influence of $\tau_{Hebb}$ and $\tau_{homeo}$ on the behavior of the dynamical learning system. To illustrate this, we examine the “simple” case of an arbitrary Hebbian-type learning rule with additional synaptic scaling. Inspired by experimental results (Turrigiano and Nelson, 2000) and used in previous studies (van Rossum et al., 2000; Yger and Harris, 2013; Zenke et al., 2013), synaptic scaling is used as a homeostatic mechanism that increases or decreases homogeneously the synaptic weights in order to reach a given firing rate $r_{target}$. In our generic formulation in Equation (2), this is equivalent to including an additive scaling term $\Gamma$ in the expression of $\Phi$ in addition to the Hebbian contribution $H$, while $\theta$ tracks the post-synaptic firing rate with a timescale $\tau_{homeo}$.

$$\begin{align} \Phi(\bar{w}, r_{\text{pre}}, r_{\text{post}}, \theta) &= H(\bar{w}, r_{\text{pre}}, r_{\text{post}}) + \Gamma(\bar{w}, \theta) \quad (10) \\ \Gamma(\bar{w}, \theta) &= \alpha \bar{w} (r_{\text{target}} - \theta) \notag \\ \Psi(r_{\text{post}}) &= r_{\text{post}} \notag \end{align}$$

For simplicity, we rewrite the Hebbian contribution using Equation (3) in terms of $\bar{w}$ only: $\tilde{H}(\bar{w}) = H(\bar{w}, r_{pre}, r_{post})$. This yields the following expression for the Jacobian in Equation (4):

$$\left( \frac{\frac{\tilde{H}'(\bar{w}) + \alpha(r_{\text{target}} - \theta)}{\tau_{\text{Hebb}}} - \frac{\alpha \bar{w}}{\tau_{\text{Hebb}}}}{\frac{\tau_{\text{Hebb}}}{\tau_{\text{homeo}}} - \frac{1}{\tau_{\text{homeo}}}} \right) \quad (11)$$

The equilibrium corresponds to the fixed point(s) where $\dot{\bar{w}} = 0$ and $\dot{\theta} = 0$, which implies that $r_{post} = \theta$ and $\tilde{H}(\bar{w}) = -\alpha \bar{w}(r_{target} - \theta)$. The trace and determinant of the Jacobian matrix are given by

$$\begin{align} T &= \frac{\frac{\tilde{H}'(\bar{w}) - \tilde{H}(\bar{w})}{\bar{w}}}{\tau_{\text{Hebb}}} - \frac{1}{\tau_{\text{homeo}}} \quad (12) \\ D &= \frac{-[\tilde{H}'(\bar{w}) - \tilde{H}(\bar{w})/\bar{w}] + \alpha r_{\text{post}}}{\tau_{\text{Hebb}} \tau_{\text{homeo}}} \notag \end{align}$$

As explained above, stability is ensured when the necessary conditions $T < 0$ and $0 < 4D < T^2$ in Equation (6) are met. These three conditions read.

$$\tilde{H}'(\bar{w}) - \frac{\tilde{H}(\bar{w})}{\bar{w}} < \frac{\tau_{\text{Hebb}}}{\tau_{\text{homeo}}} \quad (13)$$

$$\tilde{H}'(\bar{w}) - \frac{\tilde{H}(\bar{w})}{\bar{w}} < \alpha r_{\text{post}} \quad (14)$$

$$\alpha r_{\text{post}} < \frac{1}{4} \left\{ \left[ \tilde{H}'(\bar{w}) - \frac{\tilde{H}(\bar{w})}{\bar{w}} \right] \sqrt{\frac{\tau_{\text{homeo}}}{\tau_{\text{Hebb}}}} + \sqrt{\frac{\tau_{\text{Hebb}}}{\tau_{\text{homeo}}}} \right\}^2 \quad (15)$$

The term $\tilde{H}'(\bar{w}) - \tilde{H}(\bar{w})/\bar{w}$ corresponds to the sub/super-linearity of the effective weight update $\tilde{H}$ at the equilibrium $\bar{w}$, including
the effects of the neuron model. For the simplest Hebbian rule $H(\bar{w}, r_{\text{pre}}, r_{\text{post}}) = r_{\text{pre}} r_{\text{post}}, \dot{H}(\bar{w}) = r_{\text{pre}}^2 \bar{w}$ is linear and we always have $\dot{H}(\bar{w}) - \dot{H}(\bar{w})/\bar{w} = 0$. This implies that the first two conditions Equations (13) and (14) are always true, while the third condition Equation (15) reduces to $\alpha r_{\text{post}} < \tau_{\text{Hebb}}/4\tau_{\text{homeo}}$. For the synaptic scaling mechanism, $\alpha$ should be chosen sufficiently large in order to keep the output rate $r_{\text{post}}$ close to its target $r_{\text{target}}$. It follows that the third condition may be violated depending on the details of the parameters, in particular when $\tau_{\text{Hebb}} \ll \tau_{\text{homeo}}$. This corresponds to non-real eigenvalues, synonymous with oscillatory dynamics in the weights.

As a second example related to the BCM rule and triplet- STDP as will be detailed later, when $\dot{H}$ is a quadratic polynomial in $\bar{w}$ with positive second-order coefficient, we have $\dot{H}'(\bar{w}) - \dot{H}(\bar{w})/\bar{w} > 0$ for large weights. According to Equation (9), a large positive value for $\dot{H}'(\bar{w})$ implies strong competition as desired. However, the condition for the negativity of the trace in Equation (13) implies that $\tau_{\text{homeo}}$ should not be much larger than $\tau_{\text{Hebb}}$, as shown previously (Zenke et al., 2013). Then, assuming Hebbian learning to be relatively fast, Equations (14) and (15) define a limited range for the choice of $\alpha$, out of which divergence or oscillations may occur. As a conclusion, those stability and competition conditions oppose each other and make the fine tuning of the parameters necessary.

# 4. THE FAMILY OF STDP LEARNING RULES

## 4.1. Need for Regulation with Classical STDP

As a first example of learning rules, we consider the family of STDP rules to illustrate the interplay between Hebbian learning and synaptic scaling. We show that they fall into the mathematical framework developed in Section 3. To start, without any additional homeostatic regulation based on $\theta$, we recall that the convergence of the weight depends on the fixed points of $\Phi$ only. The original version of STDP simply describes the effect for pairs of input-output spikes using the well-known temporal window in Figure 1B, which determines the weight update as a function of spike-time difference. All contributions are then summed over time to obtain the total weight update. The net effect denoted by $H$ here can be decomposed into two terms, for the neuronal firing rates and covariances, respectively (Kempter et al., 1999; Gilson et al., 2009). In our framework based on the Poisson neuron (see Section 3), this gives the following differential equation for the mean weight $\bar{w}$

$$ 
\begin{aligned}
\dot{\bar{w}} &= \frac{1}{\tau_{\text{Hebb}}} H(\bar{w}, r_{\text{pre}}, r_{\text{post}}, c_{\text{pre-post}}) \\
&= \frac{1}{\tau_{\text{Hebb}}} (Ar_{\text{pre}}r_{\text{post}} + Bc_{\text{pre-post}}) \\
&= \frac{1}{\tau_{\text{Hebb}}} (Ar_{\text{pre}}^2 + Bc_{\text{pre}})\bar{w}
\end{aligned} 
\quad (16) $$

where the typical area under the curve $A < 0$ corresponds to more LTD than LTP for the rate contribution, while $B > 0$ describes LTP due to the temporal interaction for correlated

inputs. The last line is obtained using Equation (3), where the mean weight update can be expressed as a linear function of the weight from a macroscopic point of view. We obtain a first-order polynomial similar to that for classical Hebbian learning, where the coefficient depends on the input correlation. Two behaviors can occur for this system: for sufficiently strong input correlations $c_{\text{pre}}$, the factor for $\bar{w}$ becomes positive and the fixed point unstable, so positive weights are potentiated in a Hebbian fashion and diverge; otherwise weights are depressed and converge to the fixed point $\bar{w} = 0$. For a pool of synapses, competition is ensured provided $\frac{\partial H}{\partial \bar{w}} = (Ar_{\text{pre}}^2 + Bc_{\text{pre}})/\tau_{\text{Hebb}} > 0$, which occurs for sufficiently strong input correlation here. In that case, the diverging learning dynamics can result in a bimodal distribution when a positive upper bound is set (Kempter et al., 1999; Song and Abbott, 2001).

To change the fixed-point structure and enforce stability, one can add a penalty term on the weight update based on the current value of the weight (Oja, 1982). A usual example found in the literature uses a polynomial in $\bar{w}$, which leads to the following expression for $\Phi$.

$$ \dot{\bar{w}} = \frac{1}{\tau_{\text{Hebb}}} (Ar_{\text{pre}}^2 + Bc_{\text{pre}})\bar{w} - \alpha \bar{w}^n \quad (17) $$

The key point here is that $\dot{w}$ is a first-order polynomial in $w$ for classical STDP, so $n \ge 2$ stabilizes the system (Tetzlaff et al., 2011). Synaptic scaling maintains the synaptic competition while preventing weights from taking too high values, at the cost of not being able to control the post-synaptic firing rate, and without having any relationship to the real homeostatic timescale. Although that previous work studied in depth the interaction of synaptic scaling with more complex Hebbian learning rule, the temporal dynamics when the two processes are not acting on the same timescale is still poorly understood.

## 4.2. Synaptic Scaling Mechanism Targeting a Fixed Firing Rate Requires Fine Tuning

In order to target a fixed firing rate, weight normalization as previously defined is not sufficient. One must add a constraint enforcing the post-synaptic neuron to scale all its input weights such that, on average, a desired firing rate is maintained. Following previous studies (van Rossum et al., 2000; Yger and Harris, 2013; Zenke et al., 2013), it can be implemented by the term $\Gamma$ as in Equation (10), which depends on the difference between a running estimate of the post-synaptic firing rate and a desired firing rate $r_{\text{target}}$. The expression for $\Phi$ with the STDP contribution $H$ and $\Psi$ then read

$$ 
\begin{aligned}
\Phi(\bar{w}, r_{\text{post}}, c_{\text{pre-post}}, \theta) &= H(\bar{w}, r_{\text{post}}, c_{\text{pre-post}}) \\
&\quad + \alpha \bar{w}(r_{\text{target}} - \theta)
\end{aligned} 
\quad (18) $$

The constant $\alpha$ defines the strength of the homeostasis on the mean weight $\bar{w}$, while $\tau_{\text{homeo}}$ determines the timescale of the smoothing of the $r_{\text{post}}$ estimate tracked by $\theta$.

The analysis in Section 3 states that $\tau_{\text{Hebb}}$, $\tau_{\text{homeo}}$ and $\alpha$ must be chosen so as to avoid instability and trivial solutions where all
weights become silent. As it has been shown for other learning rules (Cooper et al., 2004; Zenke et al., 2013), the running estimate $\theta$ of the post-synaptic firing rate have to be rather fast, otherwise the system is subject to strong oscillations. To illustrate the problem, suppose we have a neurons targeting $r_{\text{target}} = 1$ Hz, with $r_{\text{pre}} = 0.9$ Hz, $c_{\text{pre}} = 0.1$, $A = -0.1$ and $B = 1$ (see Section 4.1). The value of $\alpha$ is varied between 0.01, 0.1, and 1. As we can see on **Figure 2A**, the convergence to the fixed point can be pretty fast if $\tau_{\text{Hebb}} = \tau_{\text{homeo}}$, and if $\alpha$ is strong enough to counterbalanced the Hebbian force that depresses synapses here; see panels with $\alpha \ge 0.1$, insets show the trajectory in the phase space $(\bar{w}, \theta)$ as function of time. However, we can see on **Figure 2B** that when $\tau_{\text{homeo}} \gg \tau_{\text{Hebb}}$, as it is found in vivo (Keck et al., 2013), strong oscillations emerge for strong value $\alpha = 1$. There is a fine tuning required between those two competing forces. To circumvent the problem, the use of a Proportional-Integral (PI) controller was incorporated in some study (van Rossum et al., 2000; Yger and Harris, 2013), but even when it prevents some oscillations from occuring, it does not abolish the requirement that $\tau_{\text{Hebb}}$ and $\tau_{\text{homeo}}$ should not be order of magnitudes apart.

## 4.3. Similar Stability Issues Occur for Weight-dependent and Triplet STDP

The analysis and the observations performed previously can be extended to several STDP-like learning rules. For example, a simple version of the weight-dependent STPD learning rule (van Rossum et al., 2000; Morrison et al., 2007; Gilson and Fukai, 2011) with linearly increasing LTD as a function of the weight and constant LTP gives

$$ 
\begin{aligned}
H(\bar{w}, r_{\text{post}}, c_{\text{pre-post}}) &= (A_+ + A_- \bar{w}) r_{\text{pre}} r_{\text{post}} \\
&\quad + B c_{\text{pre-post}} \\
&= (A_+ r_{\text{pre}}^2 + B c_{\text{pre}}) \bar{w} + A_- r_{\text{pre}}^2 \bar{w}^2
\end{aligned} 
\quad (19) $$

Again Equation (3) was used to obtain the second-order polynomial in $\bar{w}$. In **Figure 3B** that depicts the convergence of the system in a similar fashion to **Figure 2** with typical values for the parameters ($A_+ = 0.1$, $A_- = -0.3$, $B = 1$, $c_{\text{pre}} = 0.1$), the convergence is achieved if homeostatic coupling is weak ($\alpha = 0.1$). However, large oscillations arise for strong coupling ($\alpha = 1$) and when the ratio between the homeostatic and Hebbian timescales is large.

Likewise, the triplet STDP model (Pfister and Gerstner, 2006) corresponds to

$$ 
\begin{aligned}
H(\bar{w}, r_{\text{post}}, c_{\text{pre-post}}) &= (A_+ r_{\text{post}} + A_-) r_{\text{pre}} r_{\text{post}} \\
&\quad + B c_{\text{pre-post}} \\
&= (A_- r_{\text{pre}}^2 + B c_{\text{pre}}) \bar{w} + A_+ r_{\text{pre}}^3 \bar{w}^2
\end{aligned} 
\quad (20) $$

where $A_+ > 0$, $A_- < 0$ for the LTP and LTD rate contributions, respectively, as well as $B > 0$ for the correlation contribution. Again, for standard values of the parameters $A_+ = 0.05$, $A_- = -0.2$, $B = 1$, $c_{\text{pre}} = 0.1$ (Pfister and Gerstner, 2006), we see in **Figure 3B** the same qualitative behavior as with weight-dependent STDP.

The similarity can be explained by the fact that both Equations (19) and (20) are quadratic polynomials in $\bar{w}$. The difference between the two rules lies in the signs of the coefficients. Nevertheless, we have for the scaling term $\Gamma(\bar{w}, \theta) = \alpha \bar{w}(r_{\text{target}} - \theta) \simeq \alpha r_{\text{target}} \bar{w} - \alpha r_{\text{pre}} \bar{w}^2$, where we have used $\theta \simeq r_{\text{post}} = r_{\text{pre}} \bar{w}$. This means that, when $\Gamma$ overpowers the STDP contribution to enforce stability with a large $\alpha$, the coefficient for $\bar{w}$ in $\Phi$ is negative in both cases. It ensures stability, but generates similar oscillations for large $\tau_{\text{homeo}}$. The intuitive explanation is that large values for $\tau_{\text{homeo}}$ cause the gradient to have a strong horizontal component in the phase space $(\bar{w}, \theta)$ of **Figure 3**, which often implies oscillations around the fixed point.

## 4.4. Trade-off between Stability and Competition

While we analyzed the dynamical behavior of the learning rules for the mean weight to assess their implications for stability, we now examine the situation for two inputs in order to study how competition can be affected by this interaction between homeostatic and Hebbian learning. This yields an extra differential equation as explained in 3.2. **Figure 4** illustrates the evolution of the weights $w_1$ and $w_2$, as well as $\theta$, for the three learning rules previously mentioned combined with synaptic scaling, and show how competition can take place. We consider two input pathways with the same input rates $r_{1/2}$, but different levels of correlation: $c_{\text{pre}}^1 = 0.1$ and $c_{\text{pre}}^2 = 0.05$. The homeostatic mechanism targets the fixed firing rate $r_{\text{target}} = 2$ Hz. As shown in **Figures 4A,C**, strong competition is observed for both classical STDP and triplet STDP, leading to $w_2 = 0$ for the pathway with weaker correlation $c_2 < c_1$ (Kempter et al., 1999; van Rossum et al., 2000; Song and Abbott, 2001). For weight-dependent STDP, the competition is much weaker in **Figure 4B**. Nevertheless, in all cases, increasing the ratio $\tau_{\text{homeo}}/\tau_{\text{Hebb}}$ introduces oscillations of the weights during competition, exactly as previously observed for the mean weight $\bar{w}$. We also see that an increased strength for the homeostatic force ($\alpha = 0.5$ in the bottom row of **Figure 4**) does not solve the stability issue when $\tau_{\text{homeo}} \gg \tau_{\text{Hebb}}$, but causes larger fluctuations.

## 5. METAPLASTIC LEARNING RULES

The previous section showed the common trend for STDP learning rules paired with synaptic scaling targeting a desired firing rate: a large time constant to estimate the post-synaptic firing rate gives rise to instability or potentially large oscillations in the weights. Now we examine a second category of stabilizing mechanisms, where the homeostatic mechanism is implemented directly in the metaplastic learning rule; see Yeung et al. (2004) for an example for calcium-based regulation. Metaplasticity is often used to enforce a homeostatic behavior on the neural system and we will stick to this function here. Without loss of generality, we ignore correlations in the learning rules and focus on rate-based rules.
## 5.1. The Bienenstock-Cooper-Munro (BCM) Learning Rule

In order to extend rules based on correlations of rates (Oja, 1982) and approach the problem of synaptic competition via weight normalization, Bienenstock et al. (1982) designed a model of synaptic plasticity that was able to reproduce phenomenologically several observations made in vivo. Their so-called BCM rule is a physical theory of learning in the visual cortex; see Cooper and Bear (2012) for a review. The mechanism

consists in an efficient way to balance and regulate the amount of plasticity according to past activity by means of a heterosynaptic process.

Practically, a sliding threshold determines the boundary between LTP above and LTD below, and evolves according to the square of the postsynaptic firing rate (Bienenstock et al., 1982). In our formalism, this can be taken care of by a temporal tracking of $r_{post}^2$ using $\theta$ as the threshold variable with $\tau_{homeo} \gg \tau_{Hebb}$, such that $\theta \approx \langle r_{post}^2 \rangle$ with the angular brackets indicating the

<figure><img src="image_4.png" /><figcaption>**FIGURE 2** | Interplay between Hebbian and homeostatic timescales for pairwise STDP. (A) Evolution of the weight $w$ and the running estimate $\theta$ of the post-synaptic firing rate as function of time, for $\tau_{Hebb} = \tau_{homeo} = 10$ min and various gain $\alpha$ for the heterosynaptic scaling. Insets shows the trajectory in the phase space $(w, \theta)$. (B) Same as (A) with a slower homeostatic scaling: $\tau_{homeo} = 10 \tau_{Hebb} = 100$ min.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**FIGURE 3** | Interplay between Hebbian and homeostatic timescales for different learning rules and homeostatic forces. (A) Pairwise STDP with weight dependent modification in Equation (19). Left column: convergence in the phase space $(w, \theta)$ for a fast homeostatic force ($\tau_{Hebb} = \tau_{homeo} = 10$ min, upper row), or for a slow homeostatic force ($\tau_{homeo} = 10 \tau_{Hebb}$, lower row). Right column is the same, but with a stronger drive $\alpha = 1$ for the homeostatic force. (B) Same as (A) for the triplet learning rule (Pfister and Gerstner, 2006), see Equation (20).</figcaption></figure>
average over the randomness. The expression for $\Phi$ is a second-order polynomial in $r_{post}$ (Bienenstock et al., 1982), which finally gives

$$ \Phi(\bar{w}, r_{post}, \theta) = r_{pre}r_{post}(r_{post} - \theta) \\ \Psi(r_{post}) = r_{post}^2 $$

Here $\Phi$ has a similar form to that for the triplet rule in Equation (20), but the boundary between potentiation and depression is now given by $\theta$.

It is known that the BCM formalism can be subject to strong oscillations, when the timescales for the two differential equations are too far apart (Cooper et al., 2004; Toyoizumi et al., 2014). In Figure 5A, even when $\tau_{Hebb} = \tau_{homeo}$, weight oscillations are present. Moreover, for a slightly larger ratio $\tau_{Hebb}/\tau_{homeo}$, the oscillations can destroy the convergence of the system when the weights hit the lower bound 0, as illustrated in Figure 5B.

## 5.2. Modulation of STDP Depending on the Post-synaptic Firing Rate

In order to stabilize the triplet STDP rule (Pfister and Gerstner, 2006) in recurrent networks, further studies (Clopath et al., 2010;

Zenke et al., 2013) scaled the amount of LTD in terms of a smoothed average of the firing of the post-synaptic neuron. This modulation of LTD actually brings the triplet STDP rule closer to the BCM rule, by implementing a regulation of the threshold between effective LTP and LTD. In our formalism, the rule used by Zenke et al. (2013) can be implemented for rates as

$$ \Phi(\bar{w}, r_{post}, \theta) = r_{pre}r_{post}(A_+r_{post} + A_-\theta^2/r_{target}) \\ \Psi(r_{post}) = r_{post} $$

Note that the difference here compared to BCM is that $\theta$ tracks $r_{post}$ and not $r_{post}^2$, and the limit between depression and potentiation is related to $\theta^2$. As in Equation (20), we have $A_+ > 0$ and $A_- < 0$.

Figure 6 compares the evolution for this metaplastic triplet STDP rule with classical STDP combined with synaptic scaling: we can clearly see that the resulting dynamics is strongly affected by the ratio between Hebbian and homeostatic time constants in both cases. The trajectories of $w$ and $\theta$ in the same phase space as before show several types of instability, from weight (and rate) explosion for slow tracking with large $\tau_{homeo}$ (Zenke et al., 2013) to oscillations when $\tau_{homeo} = \tau_{Hebb}$. As before, slow tracking yields a gradient with a strong horizontal component,

<figure><img src="image_6.png" /><figcaption>FIGURE 4 | Competition for several plasticity rules with different timescales for Hebbian and homeostatic forces. (A) Pairwise STDP with weight-independent update. Convergence of two synaptic weights $w_{1/2}$ with different correlation inputs and the estimate of the post-synaptic firing rate, $\theta$ as function of time, for a fast homeostatic force ($\tau_{Hebb} = \tau_{homeo} = 10$ min, top row), for a slow homeostatic force ($\tau_{homeo} = 10\tau_{Hebb}$, middle row), or for a slow and stronger homeostatic force ($\alpha = 0.5$). (B) Same as A for the weight-dependent STDP learning rule (van Rossum et al., 2000). (C) Same as (A) for triplet learning rule (Pfister and Gerstner, 2006).</figcaption></figure>
hence oscillations. The limit cycle in the top panel of Figure 6B only happens for some limited range of the parameters, but this illustrates the severe instability issues even for the simple dynamical system considered here.

### 5.3. Non-linearly Gated STDP Rules

Another direction of research (Senn et al., 2001; El Boustani et al., 2012) introduced non-linearity in the effect of the Hebbian term, by turning it on and off depending on the past pre- and post-synaptic activity of the neuron. Taking a simplified version with a similar mechanism for both LTP and LTD, we consider

$$ 
\begin{aligned}
\Phi(\bar{w}, r_{\text{post}}, \theta) &= \|r_{\text{pre}}r_{\text{post}} - f_{+}(\theta)\|_{+} - \|r_{\text{pre}}r_{\text{post}} - f_{-}(\theta)\|_{+} \\
\Psi(r_{\text{post}}) &= r_{\text{pre}}r_{\text{post}}
\end{aligned}
\quad (21) $$

where $\|x\|_{+}$ is a non linear function equal to $x$ if $x > 0$ and $0$ otherwise. Now $\Psi$ is such that $\theta$ embodies a smoothed average of the pre-post correlations with the time constant $\tau_{homeo}$. When instantaneous correlations are higher than thresholds $f_{\pm}(\theta)$, for potentiation or depression respectively, plasticity effectively occurs. In the general case, $f_{\pm}$ could be any non-linear functions, and do not even need to rely on the same timescales (El Boustani et al., 2012). The simulations in Figure 6C correspond the simple case where $f_{\pm}(\theta) = a_{\pm} = \pm 0.4$ are constant. The problem with those non-linearities is that it becomes hard to perform an mathematical analysis of the equilibrium. As with other rules, we observe the same effect

of a large $\tau_{homeo}$ on the gradient and the same qualitative conclusion that slow tracking implies the slow convergence of the system.

### 5.4. Toward More Complex Models

The stability problem arises because, at the equilibrium state, the Hebbian and homeostatic mechanisms compete to balance each other, but they do not act on the same timescale. As pointed out recently (Toyoizumi et al., 2014), a solution can be found when considering that both do not interact linearly, i.e., summing their effects at the synapses, but rather work in a multiplicative manner to determine the synaptic weight. To be more precise, the model developed by Toyoizumi et al. (2014) can be integrated within our framework modulo a slightly more generic formulation for the equation in $\theta$. The model states that $w = \rho H$, where those two quantities are governed by the following system of differential equations.

$$ 
\begin{aligned}
\dot{\rho} &= \frac{1}{\tau_{\text{hebb}}} [(\rho_{\text{max}} - \rho)\|r_{\text{pre}}r_{\text{post}} - A_{+}\|_{+} \\
&\quad - (\rho - \rho_{\text{min}})\|A_{-} - r_{\text{pre}}r_{\text{post}}\|_{+}] \\
\dot{H} &= \frac{1}{\tau_{\text{homeo}}} H(1 - r_{\text{post}})
\end{aligned}
$$

Even if the lower and upper weight bounds $\rho_{min}$ and $\rho_{max}$ depend on $H$, the model can be written in a generalized version of Equation (2), using $\dot{w} = H\dot{\rho} + \dot{H}\rho$ with simply $\theta = H$. The final expression resembles non-linearly gated plasticity with an additional synaptic scaling, but involves further refinements compared to Equation (21).

<figure><img src="image_7.png" /><figcaption>FIGURE 5 | Interplay between Hebbian and homeostatic timescales for the BCM learning rule. (A) Upper row: evolution of the weight and the running estimate $\Psi$ of the post-synaptic firing rate as function of time, for $\tau_{\text{Hebb}} = \tau_{\text{homeo}} = 10$ min. Lower row: same but in the phase space $(w, \theta)$. (B) Same as (A) with $\tau_{\text{homeo}} = 2\tau_{\text{Hebb}}$.</figcaption></figure>
# 6. DISCUSSION AND PERSPECTIVES

In this paper, we have reviewed various homeostatic mechanisms that are used in recent state-of-the-art plasticity models to regulate Hebbian-type learning. We have focused on two main categories of models: (1) homeostatic synaptic scaling as an independent process that competes with the Hebbian force via an additive term, and (2) metaplastic rules, for which the Hebbian contribution is modulated in an homeostatic fashion. In both cases, the regulation is performed via an estimate of the neural activity (often the post-synaptic firing rate $r_{\text{post}}$) smoothed with a timescale $\tau_{\text{homeo}}$, whereas the Hebbian update corresponds to another timescale $\tau_{\text{Hebb}}$. We have shown for most models that, when $\tau_{\text{homeo}} \gg \tau_{\text{Hebb}}$, undesired behaviors such as oscillations in the synaptic weights occur, in particular in the case where the homeostatic force is strong. Moreover, competition and stability correspond to conflicting constraints on the parameters, which requires fine-tuning. There is thus a trade-off between the strength of the homeostatic regulation that must compete with the Hebbian drive without perturbing the convergence to a fixed point for the weights. Stability in the weights at a macroscopic level is necessary to ensure stability of the neural functions; note that we have not considered noise in the dynamics of individual weights here, but rather their mean for given pathways.

This constraint on the timescales $\tau_{\text{Hebb}}$ and $\tau_{\text{homeo}}$ is problematic in regards of available experimental data, as many of them point to slow homeostatic processes (Turrigiano et al., 1998) in comparison with Hebbian processes for which typically $\tau_{\text{Hebb}} \simeq 10$ min. Other models not considered here exhibit similar behavior, for example a homeostatic regulation obtained via intrinsic plasticity (see Zheng et al., 2013) for an example based on spike-threshold adaptation. As a conclusion, the control of the firing rate of the post-synaptic neuron should be taken care of by a mechanisms at a fast timescale, say few minutes at the maximum. Conversely, we point out that homeostatic mechanisms operating on a much slower timescale should be

related to other functions than maintaining the neural activity in a given range.

This claim is supported by several experimental and theoretical findings. Spiking activity of neurons *in vivo* is known to be sparse and highly irregular. Most V1 neurons display Poissonnian or supra Poisson spike-count variability in response to low dimensional stimuli such as bars and gratings (Dean, 1981). Even *in vitro*, they fire as Poisson sources, irregularly, with a coefficient of variation for their inter-spike intervals close to 1 (Nawrot et al., 2008). The origin of this irregular activity observed in the sub-threshold voltage and/or in spiking activity is linked to synaptic activity (Paré et al., 1998; Destexhe and Paré, 1999), and because it has been observed experimentally that excitatory and inhibitory conductances are closely balanced (Froemke et al., 2007; Okun and Lampl, 2008), such a fine balance has to be maintained by the system (Renart et al., 2010). Therefore, there is a crucial need for compensatory mechanisms that may interfere or act in concert with Hebbian learning to not only keep the neuron's firing rate within a certain range, but also guarantee this balance (Vogels et al., 2011), or the irregularity of the spiking discharge (Pozzorini et al., 2013). Weight normalization has also been studied in depth in the context of emergence of ocular dominance in order to adjust the competition between synaptic pathways, switching from winner-take-all to winner-share-all behaviors for example (Miller, 1996).

We should discuss several limitations of our study related to the proposed mathematical framework. We have focused on very simple and canonical models of synaptic plasticity, ignoring the fine morphological structure of the neurons. It was shown that the shape of the temporal learning window represented in **Figure 1B** depends on the synaptic position on the dendritic tree (Letzkus et al., 2006; Kampa et al., 2007). More importantly, homeostatic regulation or plasticity thresholds exhibit variability and affect predominantly neighboring synapses *in vivo* (Harvey and Svoboda, 2007). Therefore, we only address the temporal crosstalk between Hebbian and homeostatic plasticity at the

<figure><img src="image_8.png" /><figcaption>**FIGURE 6** | Convergence of metaplastic learning rules. (A) Convergence in the phase space $(w, \theta)$ of a classical STDP plasticity, either for a fast homeostatic time constant ($\tau_{\text{homeo}} = \tau_{\text{Hebb}}$, upper row) or for a slow one ($\tau_{\text{homeo}} = 10 \tau_{\text{Hebb}}$, lower row). (B) Same as (A) for a metaplastic learning rule combining the triplet learning rule and scaling of the LTD term (Zenke et al., 2013). (C) Same as (A) for a non-linear metaplastic learning rule including thresholds (El Boustani et al., 2012).</figcaption></figure>
<figure><img src="image_9.png" /><figcaption>**FIGURE 7** | Memory retention problem and timescales. (A) Illustration of metaplastic thresholds stabilizing the learning. Synapses are stable at in the ongoing regime, then a "plasticity trace" builds up during presentation of new sensory inputs, but this will eventually be stopped by a sliding activation threshold, allowing the synapse to adapt to those novel stimuli. (B) Illustration of the multiple timescales involved in plasticity, from the membrane time constant τm to the homeostatic one τhomeo, ranging from ms to days.</figcaption></figure>

largest scale and the question of defining the spatial extent for heterosynaptic mechanisms remains open. Nevertheless, we expect our conclusions to hold locally for groups of synapses that can be isolated and experience homogeneous processes.

In the general dynamical system in Equations (2) considered here, the timescales are explicitly defined via τHebb and τhomeo. In more complex dynamical systems involving noise and attractors, implicit time constants can emerge in a population of synapses (Tetzlaff et al., 2013). Usually, they are slow time constants though, and cannot be used for fast control of the rate, but rather to implement long-lasting memory patterns in the synaptic weights. Another limitation of our conclusions is that we only consider a feedforward model. To extend those to networks with plastic recurrent connections, the mathematical formalism should be modified to account for the case of synapses with the same pre- and post-synaptic firing rates $r_{post} = r_{pre} = r_{rec}$, and likewise the correlations $c_{pre-post} = c_{rec}$. Those quantities follow the consistency equations.

$$r_{rec} = \frac{r_{in}}{1 - w} \quad \text{and} \quad c_{rec} = \frac{c_{in}}{(1 - w)^2}$$

A similar analysis can be done to predict the behavior of learning rules and compare them. The difference compared to the feedforward case is that rates and correlations contributions to the weight updates are not of the same order. This implies that oscillations or other instability effects induced by spike synchrony are more likely to be amplified in recurrent networks than those due to firing rates. It remains that stability can be studied similarly via the Jacobian matrix. Note also that noise in firing and learning dynamics, as well as heterogeneity in neuron and network parameters, may help to prevent “pathological” weight trajectories such as limit cycles, as they smooth the dynamical landscape and degenerate too stereotypical situations.

Beyond those technical details, the puzzling question with plasticity is how synapses can store relevant information while neurons are constantly bombarded by spiking activity in vivo. This ongoing input stimulation is quite often considered to be noise in models, which impairs stability of dynamical systems over long time-scales. Although this issue has been addressed theoretically for various models (Clopath et al., 2008; Billings and van Rossum, 2009; Gilson and Fukai, 2011; Tetzlaff et al., 2013; Zenke et al., 2013), it suggests that additional timescales are necessary to properly combine short-term and long-term properties such that the system learns fast and slowly forgets. **Figure 7A** recapitulates several timescales involved in learning and memory. In essence, for the neural system to retain memories, synaptic plasticity should only be turned on by metaplasticity when “new” incoming stimuli impinge neurons. Once this novelty has been learnt, metaplasticity should stop synaptic changes. Then a selection process should trim all newly formed memories to keep only appropriate ones (Frey and Morris, 1997). This is illustrated in **Figure 7B**, where several interleaved timescales interact to bridge all mechanisms, from the effective membrane time constant τm (order of ms) that interacts with STDP to the homeostatic time constants τhomeo, which can range from hours to days (Turrigiano et al., 1998; Turrigiano and Nelson, 2004). Calcium signals can act as activity buffers at a timescale τCa²⁺ (Artola et al., 1990; Shouval et al., 2002; Yeung et al., 2004; Graupner and Brunel, 2012), whereas reward signals or neuromodulation would affect plasticity at a larger timescale τreward (Izhikevich, 2007), comparable to the one observed for Hebbian changes (τHebbian). There is also evidence for a control of intrinsic excitability, synaptic scaling at the post-synaptic density, adaptation of the pre-synaptic neurotransmitter release (Davis, 2006). As most models incorporate only a few of those at a time, we stress the need for a better understanding of the complex interactions that may arise when bringing together those experimentally observed mechanisms.
REFERENCES

Abbott, L. F., and Nelson, S. B. (2000). Synaptic plasticity: taming the beast. *Nat. Neurosci.* **3**, 1178–1183. doi: 10.1038/81453

Abraham, W. C. (2008). Metaplasticity: tuning synapses and networks for plasticity. *Nat. Rev. Neurosci.* **9**, 387. doi: 10.1038/nrn2356

Abraham, W. C., and Bear, M. F. (1996). Metaplasticity: the plasticity of synaptic plasticity. *Trends Neurosci.* **19**, 126–130.

Appleby, P. A., and Elliott, T. (2005). Synaptic and temporal ensemble interpretation of spike-timing-dependent plasticity. *Neural Comput.* **17**, 2316–2336. doi: 10.1162/0899766054796879

Artola, A., Bröcher, S., and Singer, W. (1990). Different voltage-dependent thresholds for inducing long-term depression and long-term potentiation in slices of rat visual cortex. *Nature* **347**, 69–72.

Beck, H., and Yaari, Y. (2008). Plasticity of intrinsic neuronal properties in CNS disorders. *Nat. Rev. Neurosci.* **9**, 357–369. doi: 10.1038/nrn2371

Bender, V. A., Bender, K. J., Brasier, D. J., and Feldman, D. E. (2006). Two coincidence detectors for spike timing-dependent plasticity in somatosensory cortex. *J. Neurosci.* **26**, 4166–4177. doi: 10.1523/JNEUROSCI.0176-06.2006

Bi, G. Q., and Poo, M. M. (1998). Synaptic modifications in cultured hippocampal neurons: dependence on spike timing, synaptic strength, and postsynaptic cell type. *J. Neurosci.* **18**, 10464–10472.

Bienenstock, E. L., Cooper, L. N., and Munro, P. W. (1982). Theory for the development of neuron selectivity: orientation specificity and binocular interaction in visual cortex. *J. Neurosci.* **2**, 32–48.

Billings, G., and van Rossum, M. C. W. (2009). Memory retention and spike-timing-dependent plasticity. *J. Neurophysiol.* **101**, 2775–2788. doi: 10.1152/jn.91007.2008

Bliss, T. V. P., and Lomo, T. (1973). Long-lasting potentiation of synaptic transmission in the dentate area of the anaesthetized rabbit following stimulation of the preforamt path. *J. Physiol.* **232**, 331–356.

Burkitt, A. N., Meffin, H., and Grayden, D. B. (2004). Spike-timing-dependent plasticity: the relationship to rate-based learning for models with weight dynamics determined by a stable fixed point. *Neural Comput.* **16**, 885–940. doi: 10.1162/089976604773135041

Caporale, N., and Dan, Y. (2008). Spike timing-dependent plasticity: a Hebbian learning rule. *Annu. Rev. Neurosci.* **31**, 25–46. doi: 10.1146/annurev.neuro.31.060407.125639

Christie, B. R., and Abraham, W. C. (1992). Priming of associative long-term depression in the dentate gyrus by theta frequency synaptic activity. *Neuron* **9**, 79–84.

Clopath, C., Büsing, L., Vasilaki, E., and Gerstner, W. (2010). Connectivity reflects coding: a model of voltage-based STDP with homeostasis. *Nat. Neurosci.* **13**, 344–352. doi: 10.1038/nn.2479

Clopath, C., Ziegler, L., Vasilaki, E., Büsing, L., and Gerstner, W. (2008). Tag-trigger-consolidation: a model of early and late long-term-potentiation and depression. *PLoS Comput. Biol.* **4**:e1000248. doi: 10.1371/journal.pcbi.1000248

Cooper, L. N., and Bear, M. F. (2012). The BCM theory of synapse modification at 30: interaction of theory with experiment. *Nat. Rev. Neurosci.* **13**, 798–810. doi: 10.1038/nrn3353

Cooper, L. N., Intrator, N., Blais, B. S., and Shouval, H. Z. (2004). *Theory of Cortical Plasticity*. New Jersey, NJ: World Scientific.

Corlew, R., Wang, Y., Ghermazien, H., Erisir, A., and Philpot, B. D. (2007). Developmental switch in the contribution of presynaptic and postsynaptic NMDA receptors to long-term depression. *J. Neurosci.* **27**, 9835–9845. doi: 10.1523/jneurosci.5494-06.2007

Davis, G. W. (2006). Homeostatic control of neural activity: from phenomenology to molecular design. *Annu. Rev. Neurosci.* **29**, 307–323. doi: 10.1146/annurev.neuro.28.061604.135751

Dean, A. F. (1981). The relationship between response amplitude and contrast for cat striate cortical neurones. *J. Physiol.* **318**, 413–427.

Destexhe, A., and Paré, D. (1999). Impact of network activity on the integrative properties of neocortical pyramidal neurons *in vivo*. *J. Neurophysiol.* **81**, 1531–1547.

El Boustani, S., Yger, P., Frégnac, Y., and Destexhe, A. (2012). Stable learning in stochastic network states. *J. Neurosci.* **32**, 194–214. doi: 10.1523/JNEUROSCI.2496-11.2012

Feldman, D. E. (2012). The spike-timing dependence of plasticity. *Neuron* **75**, 556–571. doi: 10.1016/j.neuron.2012.08.001

Frey, U., and Morris, R. G. (1997). Synaptic tagging and long-term potentiation. *Nature* **385**, 533–536.

Froemke, R. C., Merzenich, M. M., and Schreiner, C. E. (2007). A synaptic memory trace for cortical receptive field plasticity. *Nature* **450**, 425–429. doi: 10.1038/nature06289

Froemke, R. C., Poo, M.-M., and Dan, Y. (2005). Spike-timing-dependent synaptic plasticity depends on dendritic location. *Nature* **434**, 221–225. doi: 10.1038/nature03366

Froemke, R. C., Tsay, I. A., Raad, M., Long, J. D., and Dan, Y. (2006). Contribution of individual spikes in burst-induced long-term synaptic modification. *J. Neurophysiol.* **95**, 1620–1629. doi: 10.1152/jn.00910.2005

Froemke, R. C. R., and Dan, Y. (2002). Spike-timing-dependent synaptic modification induced by natural spike trains. *Nature* **416**, 433–438. doi: 10.1038/416433a

Gilson, M., Burkitt, A. N., Grayden, D. B., Thomas, D. A., and van Hemmen, J. L. (2009). Emergence of network structure due to spike-timing-dependent plasticity in recurrent neuronal networks. I. Input selectivity-strengthening correlated input pathways. *Biol. Cybern.* **101**, 81–102. doi: 10.1007/s00422-009-0319-4

Gilson, M., and Fukai, T. (2011). Stability versus neuronal specialization for STDP: long-tail weight distributions solve the dilemma. *PLoS ONE* 6:e25339. doi: 10.1371/journal.pone.0025339

Graupner, M., and Brunel, N. (2012). Calcium-based plasticity model explains sensitivity of synaptic changes to spike pattern, rate, and dendritic location. *Proc. Natl. Acad. Sci. U.S.A.* **109**, 21551–21551. doi: 10.1073/pnas.1109359109

Gütig, R., Aharonov, R., Rotter, S., and Sompolinsky, H. (2003). Learning input correlations through nonlinear temporally asymmetric Hebbian plasticity. *J. Neurosci.* **23**, 3697–3714.

Harvey, C. D., and Svoboda, K. (2007). Locally dynamic synaptic learning rules in pyramidal neuron dendrites. *Nature* **450**, 1195–1200. doi: 10.1038/nature06416

Hebb, D. O. (1949). *The Organization of Behavior: A Neuropsychological Theory*, Vol. 44. New York, NY: Science Education.

Hengen, K. B., Lambo, M. E., Van Hoooser, S. D., Katz, D. B., and Turrigiano, G. G. (2013). Firing rate homeostasis in visual cortex of freely behaving rodents. *Neuron* **80**, 335–342. doi: 10.1016/j.neuron.2013.08.038

Huang, Y.-Y., Colino, A., Selig, D. K. D., and Malenka, R. C. R. (1992). The influence of prior synaptic activity on the induction of long-term potentiation. *Science* **255**, 730–733.

Izhikevich, E. M. (2007). Solving the distal reward problem through linkage of STDP and dopamine signaling. *Cereb. Cortex* **17**, 2443–2452. doi: 10.1126/science.1346729

Kampa, B. M., Letzkus, J. J., and Stuart, G. J. (2007). Dendritic mechanisms controlling spike-timing-dependent synaptic plasticity. *Trends Neurosci.* **30**, 456–463. doi: 10.1016/j.tins.2007.06.010

Keck, T., Keller, G. B., Jacobsen, R. I., Eysel, U. T., Bonhoeffer, T., and Hübener, M. (2013). Synaptic scaling and homeostatic plasticity in the mouse visual cortex *in vivo*. *Neuron* **80**, 327–334. doi: 10.1016/j.neuron.2013.08.018

Kempter, R., Gerstner, W., Van Hemmen, J. L., and Hemmen, J. V. (1999). Hebbian learning and spiking neurons. *Phys. Rev. E* **59**, 4498–4514.

Letzkus, J. J., Kampa, B. M., and Stuart, G. J. (2006). Learning rules for spike timing-dependent plasticity depend on dendritic synapse location. *J. Neurosci.* **26**, 10420–10429. doi: 10.1523/JNEUROSCI.2650-06.2006

Marder, E., and Goaillard, J.-M. (2006). Variability, compensation and homeostasis in neuron and network function. *Nat. Rev. Neurosci.* **7**, 563–574. doi: 10.1038/nrn1949

Markram, H., Lübke, J., Frotscher, M., and Sakmann, B. (1997). Regulation of synaptic efficacy by coincidence of postsynaptic APs and EPSPs. *Science* **275**, 213–215.

Miller, K. (1996). Synaptic economics: competition and cooperation in synaptic plasticity. *Neuron* **17**, 371–374.

Mockett, B., Coussens, C., and Abraham, W. C. (2002). NMDA receptor-mediated metaplasticity during the induction of long-term depression by low-frequency stimulation. *Eur. J. Neurosci.* **15**, 1819–1826. doi: 10.1046/j.1460-9568.2002.02008.x
Morrison, A., Aertsen, A., and Diesmann, M. (2007). Spike-timing-dependent plasticity in balanced random networks. *Neural Comput.* 19, 1437–1467. doi: 10.1162/neco.2007.19.6.1437

Nabavi, S., Fox, R., Proulx, C. D., Lin, J. Y., Tsien, R. Y., and Malinow, R. (2014). Engineering a memory with LTD and LTP. *Nature* 511, 348–352. doi: 10.1038/nature13294

Nawrot, M. P., Boucsein, C., Rodriguez Molina, V., Riehle, A., Aertsen, A., and Rotter, S. (2008). Measurement of variability dynamics in cortical spike trains. *J. Neurosci.* Methods 169, 374–390. doi: 10.1016/j.jneumeth.2007.10.013

Nevian, T., and Sakmann, B. (2006). Spine Ca2+ signaling in spike-timing-dependent plasticity. *J. Neurosci.* 26, 11001–11013. doi: 10.1523/JNEUROSCI.1749-06.2006

Oja, E. (1982). Simplified neuron model as a principal component analyzer. *J. Math. Biol.* 15, 267–273.

Okun, M., and Lampl, I. (2008). Instantaneous correlation of excitation and inhibition during ongoing and sensory-evoked activities. *Nat. Neurosci.* 11, 535–537. doi: 10.1038/nn.2105

O'Leary, T., and Wyllie, D. J. (2011). Neuronal homeostasis: time for a change? *J. Physiol.* 589, 4811–4826. doi: 10.1113/jphysiol.2011.210179

Paré, D., Shink, E., Gaudreau, H., Destexhe, A., and Lang, E. J. (1998). Impact of spontaneous synaptic activity on the resting properties of cat neocortical pyramidal neurons in vivo. *J. Neurophysiol.* 79, 1450–60.

Pfister, J.-P., and Gerstner, W. (2006). Triplets of spikes in a model of spike timing-dependent plasticity. *J. Neurosci.* 26, 9673–9682. doi: 10.1523/JNEUROSCI.1425-06.2006

Pozo, K., and Goda, Y. (2010). Unraveling mechanisms of homeostatic synaptic plasticity. *Neuron*, 66, 337–351. doi: 10.1016/j.neuron.2010.04.028

Pozzorini, C., Naud, R., Mensi, S., and Gerstner, W. (2013). Temporal whitening by power-law adaptation in neocortical neurons. *Nat. Neurosci.* 16, 942–948. doi: 10.1038/nn.3431

Renart, A., de la Rocha, J., Bartho, P., Hollender, L., Parga, N., Reyes, A., and Harris, K. D. (2010). The asynchronous state in cortical circuits. *Science* 327, 587–590. doi: 10.1126/science.1179850

Rochester, N., Holland, J., Haibt, L., and Duda, W. (1956). Tests on a cell assembly theory of the action of the brain, using a large digital computer. *IRE Trans. Inf. Theory* 2, 80–93. doi: 10.1109/TIT.1956.1056810

Senn, W., Markram, H., and Tsodyks, M. (2001). An algorithm for modifying neurotransmitter release probability based on pre- and postsynaptic spike timing. *Neural Comput.* 13, 35–67. doi: 10.1162/089976601300014628

Shouval, H. Z., Bear, M. F., and Cooper, L. N. (2002). A unified model of NMDA receptor-dependent bidirectional synaptic plasticity. *Proc. Natl. Acad. Sci. U.S.A.* 99, 10831–10836. doi: 10.1073/pnas.152343099

Sjöström, P. J., Turrigiano, G. G., and Nelson, S. B. (2001). Rate, timing, and cooperativity jointly determine cortical synaptic plasticity. *Neuron* 32, 1149–1164. doi: 10.1016/S0896-6273(01)00542-6

Sjöström, P. J., Turrigiano, G. G., and Nelson, S. B. (2003). Neocortical LTD via coincident activation of presynaptic NMDA and cannabinoid receptors. *Neuron* 39, 641–654. doi: 10.1016/S0896-6273(03)00476-8

Song, S., and Abbott, L. F. (2001). Cortical development and remapping through spike timing-dependent plasticity. *Neuron* 32, 339–350. doi: 10.1016/S0896-6273(01)00451-2

Tetzlaff, C., Kolodziejski, C., Timme, M., Tsodyks, M., and Wörgötter, F. (2013). Synaptic scaling enables dynamically distinct short- and long-term memory formation. *PLoS Comput. Biol.* 9:e1003307. doi: 10.1371/journal.pcbi.10.03307

Tetzlaff, C., Kolodziejski, C., Timme, M., and Wörgötter, F. (2011). Synaptic scaling in combination with many generic plasticity mechanisms stabilizes circuit connectivity. *Front. Comput. Neurosci.* 5:47. doi: 10.3389/fncom.2011.00047

Toyoizumi, T., Kaneko, M., Stryker, M. P., and Miller, K. D. (2014). Modeling the dynamic interaction of Hebbian and homeostatic plasticity. *Neuron* 84, 497–510. doi: 10.1016/j.neuron.2014.09.036

Turrigiano, G. G. (2008). The self-tuning neuron: synaptic scaling of excitatory synapses. *Cell* 135, 422–435. doi: 10.1016/j.cell.2008.10.008

Turrigiano, G. G., Leslie, K. R., Desai, N. S., Rutherford, L. C., and Nelson, S. B. (1998). Activity-dependent scaling of quantal amplitude in neocortical neurons. *Nature* 391, 892–896.

Turrigiano, G. G., and Nelson, S. B. (2000). Hebb and homeostasis in neuronal plasticity. *Curr. Opin. Neurobiol.* 10, 358–364. doi: 10.1016/S0959-4388(00)00091-X

Turrigiano, G. G., and Nelson, S. B. (2004). Homeostatic plasticity in the developing nervous system. *Nat. Rev. Neurosci.* 5, 97–107. doi: 10.1038/nrn1327 van Rossum, M. C., Bi, G. Q., and Turrigiano, G. G. (2000). Stable Hebbian learning from spike timing-dependent plasticity. *J. Neurosci.* 20, 8812–8821.

Vitureira, N., and Goda, Y. (2013). The interplay between Hebbian and homeostatic synaptic plasticity. *J. Cell Biol.* 203, 175–186. doi: 10.1083/jcb.201306030

Vogels, T. P., Sprekeler, H., Zenke, F., Clopath, C., and Gerstner, W. (2011). Inhibitory plasticity balances excitation and inhibition in sensory pathways and memory networks. *Science* 334, 1569–1573. doi: 10.1126/science.1211095

von der Malsburg, C. (1973). Self-organization of orientation sensitive cells in the striate cortex. *Kybernetik* 14, 85–100.

Wang, H.-X., Gerkin, R. C., Nauen, D. W., and Bi, G.-Q. (2005). Coactivation and timing-dependent integration of synaptic potentiation and depression. *Nat. Neurosci.* 8, 187–193. doi: 10.1038/nn1387

Wang, L., Fontanini, A., and Maffei, A. (2012). Experience-dependent switch in sign and mechanisms for plasticity in layer 4 of primary visual cortex. *J. Neurosci.* 32, 10562–10573. doi: 10.1523/JNEUROSCI.0622-12.2012

Wang, X. J. (1998). Calcium coding and adaptive temporal computation in cortical pyramidal neurons. *J. Neurophysiol.* 79, 1549–1566.

Watt, A. J., and Desai, N. S. (2010). Homeostatic plasticity and STDP: keeping a neuron's cool in a fluctuating world. *Front. Synaptic Neurosci.* 2:5. doi: 10.3389/fnsyn.2010.00005

Yeung, L. C., Shouval, H. Z., Blais, B. S., and Cooper, L. N. (2004). Synaptic homeostasis and input selectivity follow from a calcium-dependent plasticity model. *Proc. Natl. Acad. Sci. U.S.A.* 101, 14943–14948. doi: 10.1073/pnas.0405555101

Yger, P., and Harris, K. D. (2013). The Convallis rule for unsupervised learning in cortical networks. *PLoS Comput. Biol.* 9:e1003272. doi: 10.1371/journal.pcbi.1003272

Young, J. M., Waleszczyk, W. J., Wang, C., Calford, M. B., Dreher, B., and Obermayer, K. (2007). Cortical reorganization consistent with spike timing-but not correlation-dependent plasticity. *Nat. Neurosci.* 10, 887–895. doi: 10.1038/nn1913

Zenke, F., Hennequin, G., and Gerstner, W. (2013). Synaptic plasticity in neural networks needs homeostasis with a fast rate detector. *PLoS Comput. Biol.* 9:e1003330. doi: 10.1371/journal.pcbi.1003330

Zhang, W., and Linden, D. J. (2003). The other side of the engram: experience-driven changes in neuronal intrinsic excitability. *Nat. Rev. Neurosci.* 4, 885–900. doi: 10.1038/nrn1248

Zheng, P., Dimitrakakis, C., and Triesch, J. (2013). Network self-organization explains the statistics and dynamics of synaptic connection strengths in cortex. *PLoS Comput. Biol.* 9:e1002848. doi: 10.1371/journal.pcbi.1002848

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2015 Yger and Gilson. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.