Elsevier

Biosystems

Volume 79, Issues 1–3, January–March 2005, Pages 11-20
Biosystems

Dynamics of pruning in simulated large-scale spiking neural networks

https://doi.org/10.1016/j.biosystems.2004.09.016Get rights and content

Abstract

Massive synaptic pruning following over-growth is a general feature of mammalian brain maturation. This article studies the synaptic pruning that occurs in large networks of simulated spiking neurons in the absence of specific input patterns of activity. The evolution of connections between neurons were governed by an original bioinspired spike-timing-dependent synaptic plasticity (STDP) modification rule which included a slow decay term. The network reached a steady state with a bimodal distribution of the synaptic weights that were either incremented to the maximum value or decremented to the lowest value. After 1×106 time steps the final number of synapses that remained active was below 10% of the number of initially active synapses independently of network size. The synaptic modification rule did not introduce spurious biases in the geometrical distribution of the remaining active projections. The results show that, under certain conditions, the model is capable of generating spontaneously emergent cell assemblies.

Introduction

Massive synaptic pruning following over-growth is a general feature of mammalian brain maturation Rakic et al., 1986, Zecevic and Rakic, 1991. Pruning starts near time of birth and is completed by time of sexual maturation. Biological mechanisms that regulate pruning mechanisms involve complex neurochemical pathways of cell signaling and are not intended to be reviewed here. Trigger signals able to induce synaptic pruning could be related to dynamic functions that depend on the timing of action potentials. Spike-timing-dependent synaptic plasticity (STDP) is a change in the synaptic strength based on the ordering of pre- and post-synaptic spikes. This mechanism has been proposed to explain the origin of long-term potentiation (LTP), i.e. a mechanism for reinforcement of synapses repeatedly activated shortly before the occurrence of a post-synaptic spike Kelso et al., 1986, Bi and Poo, 1998, Froemke and Dan, 2002, Kepecs et al., 2002, Markram et al., 1997. STDP has also been proposed to explain long-term depression (LTD), which corresponds to the weakening of synapses strength whenever the pre-synaptic cell is repeatedly activated shortly after the occurrence of a post-synaptic spike (Karmarkar and Buonomano, 2002).

The glutamatergic NMDA receptors were initially identified as the receptor site with all biological features compatible with LTP induced by coincident pre- and post-synaptic cell discharges (Wigstrom and Gustafsson, 1986). The involvement of NMDA receptors in timing-dependent long-term depression (tLTD) has been recently described (Sjstrm et al., 2003). Recent investigations suggest that glutamatergic receptors with AMPA channels and GABAergic receptors may also undergo modifications of the corresponding post-synaptic potentials as a function of the timing of pre- and post-synaptic activities Engel et al., 2001, Woodin et al., 2003. These studies suggest that several mechanisms mediated by several neurotransmitters may exist at the synaptic level for changing the post-synaptic potential, either excitatory or inhibitory, as a function of the relative timing of pre- and post-synaptic spikes.

The important consequences that changes in synaptic strength may produce for information transmission, and subsequently for synaptic pruning, have raised an interest to simulate the activity of neural networks with embedded synapses characterized by STDP Lumer et al., 1997, Fusi et al., 2000, Hopfield and Brody, 2004 The relation between synaptic efficacy and synaptic pruning Chechik et al., 1999, Mimura et al., 2003, suggest that the weak synapses may be modified and removed through competitive “learning” rules. Competitive synaptic modification rules maintain the average neuronal input to a post-synaptic neuron, but provoke selective synaptic pruning in the sense that converging synapses are competing for control of the timing of post-synaptic action potentials Song et al., 2000, Song and Abbott, 2001.

This article studies the synaptic pruning that occurs in a large network of simulated spiking neurons in the absence of specific input patterns. The originality of our study stands on the size of the network, up to 10,000 units, the duration of the experiment, 1,000,000 time units (one time unit corresponding to the duration of a spike), and the application of an original bioinspired STDP modification rule compatible with hardware implementation Eriksson et al., 2003, Tyrrell et al., 2003. The network is composed of a mixture of excitatory and inhibitory connections that maintain a balanced input locally connected in a random way.

STDP is considered an important mechanism that modifies the gain of several types of synapses in the brain. In this study the synaptic modification rule was applied only to the excitatory–excitatory connections. This plasticity rule might produce the strengthening of the connections among neurons that belong to cell assemblies characterized by recurrent patterns of firing. Conversely, those connections that are not recurrently activated might decrease in efficiency and eventually be eliminated. The main goal of our study is to determine whether or not, and under which conditions, such cell assemblies may emerge from a large neural network receiving background noise and content-related input organized in both temporal and spatial dimension. In order to reach this goal the first step consisted in characterizing the dynamics of synaptic pruning in absence of content-related input. This first step is described here.

Section snippets

Network connectivity

The network is a 2D lattice folded as a torus to limit the edge effect where the units near the boundary received less input. The size of the network varies between 10×10 to 100×100 units. Several types of units may be defined. In this study we define two types, q{1,2} , 80% of Type I (q=1) units and 20% of Type II (q=2) units are uniformly distributed over the network according to a space-filling quasi-random Sobol distribution (Press et al., 1992, Fig. 7.7.1). A unit of either type may

Results

All synapses of type [1,1] , i.e. (excexc) were initialized with Aij(t=0)=[A2] . In presence of background activity only, most synapses were characterized by a decrement of the activation level. After a long time, t=tsteady , the network activity is stabilized and STDP does not modify any more the activation level of the synapses. At time t=tsteady most modifiable synapses were eliminated and almost all remaining active synapses were characterized by the highest possible activation level, i.e.

Discussion

We assumed a number of simplified hypotheses about the presence of only two types of units, their leaky integrate-and-fire dynamics, their distribution and the dynamics of the transfer functions of the synapses that connect these units. With all these assumptions we observed that the network reached a steady state when the synaptic weights were either incremented to the maximum value or decremented to the lowest value. Our result is in agreement with the bimodal distribution of synaptic

Acknowledgements

The authors thank Dr. Yoshiyuki Asai for discussions and comments on the manuscript. This work is partially funded by the Future and Emerging Technologies program (ISTFET) of the European Community, under grant IST-2000-28027 (POETIC), and under grant OFES 00.0529-2 by the Swiss Government.

References (32)

  • R.C. Froemke et al.

    Spike-timing-dependent synaptic modification induced by natural spike trains

    Nature

    (2002)
  • S. Fusi et al.

    Spike-driven synaptic plasticity: theory, simulation, VLSI implementation

    Neural Comput.

    (2000)
  • S.L. Hill et al.

    Dynamic transitions in global network activity influenced by the balance of excitation and inhibition

    Network: Computat. Neural Syst.

    (1997)
  • J.J. Hopfield et al.

    What is a moment? “Cortical” sensory integration over a brief interval

    Proc. Natl. Acad. Sci. USA

    (2000)
  • J.J. Hopfield et al.

    Learning rules and network repair in spike-timing-based computation networks

    Proc. Natl. Acad. Sci. USA

    (2004)
  • U.R. Karmarkar et al.

    A model of spike-timing dependent plasticity: one or two coincidence detectors?

    J. Neurophysiol.

    (2002)
  • Cited by (91)

    • Growth strategy determines the memory and structural properties of brain networks

      2021, Neural Networks
      Citation Excerpt :

      The mammalian brain is initially formed through an initial rapid proliferation of synapses. Synaptic density thus reaches a peak during early infancy and from then on it begins a steady decline down to about half this value later in life, in a process known as synaptic pruning (Chechik, Meilijson, & Ruppin, 1998; Iglesias, Eriksson, Grize, Tomassini, & Villa, 2005). It is believed that the reason for reducing synaptic density is becoming more energetically efficient (Chechik, Meilijson, & Ruppin, 1999; Stepanyants, Hof, & Chklovskii, 2002).

    • SNAVA—A real-time multi-FPGA multi-model spiking neural network simulation architecture

      2018, Neural Networks
      Citation Excerpt :

      The majority of instructions are carried out in a single clock cycle in all Processing Elements. The following evaluations were carried out to obtain the SNAVA performance figures in terms of processing speed and spike distribution time, by considering the simulation of Iglesias and Villa model (Iglesias et al., 2005), Izhikevich model (Izhikevich, 2003) and Leaky integrate-and-fire model in 16-bit fixed point arithmetic operations. The algorithms for the simulation of the mentioned models have been programmed in assembler code in order to achieve the maximum efficiency in terms of the execution time.

    View all citing articles on Scopus
    View full text