Gangs of Boomtown provides a platform for certain missions. These missions are based on two interface lines. These two inerface lines are:-
  • Reward obtaining interface line
  • Storyline Proceedings

Rewards obtained during missions have often proved beneficial to a player. Therefore, it becomes mandatory to first sort out these missions and then proceed to personal play (free play). On the other hand, Certain missions have also rewards but are attained after reaching an appropriate level. They always appear in a substabntial form. These missions proceeds the storyline of the game, that is, killing all the members of the Assassins Gang in maintained intervals. This interface of the missions in Gangs of Boomtown often co-relates with certain requirements of multiple targets. These requirements are greatly associated with the main target, the Assassin, also called the PRIMARY TARGET.

Besides considering the missions on the criteria of thier interface, they are also categorised on the criteria of different stages. There are two categories of missions when considering this criteria. These are:-
  • Tutorial Missions (missions available at the very initial stages of the game), also referred to as beginner's missions. There are 6 tutorial missions.
  • Moderate Missions (missions that are followed after the tutorial missions are completed. These missions are also called normal game missions). There are 105 moderate missions.

Collectively, Taking all the missions as a whole, there are 111 missions.

Moreover, There are certain missions that are subjected to the buildings of the game and certain missions that are with regard to the duels. When these missions are met, rewards can be collected. Another aspect, when missions that are directly related to the upgrade of Guns of Boomtown Building are also known as exclusive missions as they are always followed after the missions subjected to the storyline (i.e., killing the assassin) are completed. These exclusive missions provide a heavy reward.

Missions are always stacked along the left corner of the game display called mission stack. Here missions that can be viewed normally are limited to four. When missions are more than this normal limit, that is four missions for display, there will be scroll buttons available. Missions that generally show "NEW MISSION" icon near them indicates a player that the arrival of the mission is new and thus not accepted. When the missions are accepted, they present the requirements for completion. However, the programmers and the developers of the game have updated "information upon scroll". This works as follows:-
When you scroll your mouse (computers) or mouse pads (laptops commonly), the information, actually a hint, will appear. This does not only happens with missions but also other requirements of the game that are not related to the missions such as the requirements of a particular buildings (Guns of Boomtown Building, Population subjected building and Housing subjected builings).
At times, generally observed, missions come new but they aren't successful enough to catch the attention of the player. This is because, when there are already four missions in the mission stack and the new mission joins the queue, it takes place as the fifth mission. There will appear a scroll button. Upon scrolling itself the fifth mission catches the attention of the player. If missions that are stacked in the mission stack of the player are only 3, then the new missions that arrive catches the attention of the player immediately as it arrives in the mission stack taking the place as the fourth mission.

There is another fact that is normally considered among players. This fact is: WHEN THE ENERGY IS FULL (25 complete energy unit), AND THE MISSION REWARDS YOU WITH 7 EXTRA ENERGY UNITS, THEN THE ENERGY UNITS IN THE ENERGY BAR REMAINS THE SAME AND THE 7 EXTRA ENERGY UNITS GO WASTE?

Answer to this is quite simple. Infact, the game is designed in such a way that it corresponds to the commands where the missions confer it. When you have 25 energy units present in your energy bar and you then receive seven addtional energy units, the energy bar wont remain the same. The missions confer the system commands that the energy wont proceed with respond to regeneration time when the energy is to the limit (25 energy units), by adding extra 7 energy units to the already existing ones. In this case, you will have a total of 32 energy units filled in your energy bar. If you loose even one of these units, the regeneration time will not occur because the regeneration time only occurs when the energy unit has reached 24.

Missions are also based on social aspects. Missions that relate to the requirements regarding neighbour help are in this category.


Science fiction convention used to circumvent the fact that, according to conventional science, FTL space travel is probably impossible. A generation ship is a vast space ship that travels at sub-light speeds and is capable of sustaining human life for centuries or more. The idea is that if it takes generations to reach another star system then you do just that: the crew live their lives aboard the ship, breed, and then their children take control of the ship and then their children and so on until the generation ship actually arrives at its destination.


Separation Agreement

Is it legal to draw up our own separation agreement and file it with the courts with some sort of clause that states we will divide assets at a later time?

Separation Agreements are not filed with the court, and you must divide the assets prior to divorce via separation agreement, or file a claim for ED at that time.


Journal of Computational Neuroscience

, Volume 39, Issue 3, pp 255–270 | Cite as

Impact of sub and supra-threshold adaptation currents in networks of spiking neurons

  • David Colliaux
  • Pierre Yger
  • Kunihiko Kaneko
Open Access

Neuronal adaptation is the intrinsic capacity of the brain to change, by various mechanisms, its dynamical responses as a function of the context. Such a phenomena, widely observed in vivo and in vitro, is known to be crucial in homeostatic regulation of the activity and gain control. The effects of adaptation have already been studied at the single-cell level, resulting from either voltage or calcium gated channels both activated by the spiking activity and modulating the dynamical responses of the neurons. In this study, by disentangling those effects into a linear (sub-threshold) and a non-linear (supra-threshold) part, we focus on the the functional role of those two distinct components of adaptation onto the neuronal activity at various scales, starting from single-cell responses up to recurrent networks dynamics, and under stationary or non-stationary stimulations. The effects of slow currents on collective dynamics, like modulation of population oscillation and reliability of spike patterns, is quantified for various types of adaptation in sparse recurrent networks.

Adaptation Neuronal dynamics

1 Introduction

Most neurons in primary sensory areas tend to change the strength of their dynamical responses over time for sustained and constant inputs, in a so-called adaptation process. The detailed mechanisms of this adaptation are still not clearly understood, and can result from various phenomena that might be combined: homeostasis or intracellular mechanisms (Turrigiano and Nelson 2004; Benda and Herz 2003), short term plasticity (Tsodyks et al. 2000), or even network-wide effect originating from lateral connections (Haider et al. 2010). In this paper, we propose to investigate computationally some effects of neuronal adaptation, from both a single-cell and a network point of view. By using a phenomenological model for neurons based on an integrate-and-fire model with intrinsic adaptation and its macroscopic counterparts, we studied the effects of slow adaptation currents on neuronal dynamics at the network level.

Detailed models of intrinsic plasticity as a source of homeostasis and neuronal adaptation have already been investigated. Such biological models are used to explore the biomechanistic effects of slow ionic channels on the microscopic cortical activity of single cells (Benda and Herz 2003; Gigante et al. 2007). On a more phenomenological level, they can be used to study the firing rates dynamics (Treves 1993) or constrain their distribution (Benucci et al. 2013) in macroscopic models of cortical networks. All those models, however, are mostly based on two types of current involved in slow adaptation and related to the flow of potassium and calcium ions through the membrane. The interaction of the adaptation currents with the neuronal dynamics is complex, impacting both the firing threshold and the behaviour of the cell (Benda et al. 2010). On one hand, voltage-gated potassium currents, like Kv1, are activated by a mild depolarization of the membrane potential and control the propagation of spikes by modulating the spike threshold (Higgs and Spain 2011), a phenomena also referred to as accommodation. On the other hand, calcium-gated potassium currents are activated only at high levels of depolarization and result in after-hyper-polarization (AHP), a drop of the membrane potential and long lasting decrease of excitability after a spike is emitted (Andrade 2011). The time scale of these adaptation mechanisms can expand over a wide range up to minutes (Pozzorini et al. 2013) and are very heterogeneously distributed among cells but especially pronounced in pyramidal cells (Nowak et al. 2003).

Combined all together, those two distinct adaptation mechanisms can prevent the saturation in the spike generation process of a neuron or modulate its synchronization properties, both acting in a different manner as shown in previous experimental (Deemyad et al. 2012) or theoretical works (Ermentrout 1998; Ermentrout et al. 2001; Prescott and Sejnowski 2008; La Camera et al. 2004). The voltage-gated sub threshold current shifts the input threshold for triggering spike depending on the basal activity: it can thus be seen as a good mechanism to explain that some response properties of a neuron are independent of the background inputs. It has also been shown to play a role in the homeostasis after sensory deprivation (Nataraj et al. 2010). Conversely, the calcium-gated current triggered by a spike does not only affect the threshold but also increase the dynamic range of the neuronal response and thus avoids saturation. Based on those observations, we will refer to those two different mechanisms of adaptation as a linear or sub-threshold one (voltage-gated channels) versus a non-linear or supra-threshold one (calcium-gated channels). The role of adaptation in enhancing reliability of spike-timing of neurons stimulated with periodic inputs was also studied in Schreiber et al. (2004). Adaptation could also act functionally as a decorrelation machine (Wang 1998).

The dynamics of cortical networks are considerably enriched with adaptation currents, inducing bursting activity, slow oscillations, and homeostasis (Tabak et al. 2000; Giugliano et al. 2004; Gigante et al. 2007). Rhythmic activity of central pattern generators involved in locomotor behaviour (Grillner 1997), or multi-stable dynamics of cortical networks when ambiguous stimulus is presented (Wilson 2003) were related to slow potassium currents. More recently it has been shown that such currents should be included to explain the rich repertoire of ongoing activity observed in vivo (Mattia and Sanchez-Vives 2012) and which could be the microscopic substrate of the resting state activity recorded in fMRI (Gusnard and Raichle 2001). It is therefore crucial to have a better understanding of the effect of this adaptation at a single-cell or at a population level.

In this study, we investigated those physiological observations on adaptation in a model of cortical dynamics simple enough, so that large scale simulations can be performed. To do so, we used an adaptive exponential integrate-and-fire neuron model suited for large-scale simulations of cortical networks (Brette and Gerstner 2005) and known to be complex enough to reproduce all the dynamical repertoire recorded in vitro in various cell types (Izhikevich 2001; Brette and Gerstner 2005). Indeed, the dynamics of networks of such units have been recently investigated in Ladenbauer et al. (2014) and Farkhooi et al. (2011) and this model have been successful in capturing more diverse dynamics by generating a slow inhibitory feedback, reflecting the fact that neurons tends to adapt when stimulated with a constant inputs. While a classical model would provide a sustained response, models with adaptation will have response closer to what is observed in biological recordings. By studying the sub-threshold (linear) and the supra-threshold (non-linear) effects of the adaptation on single-cell response or in a neuronal network, we were able to disentangle the functional role of those two components on aspects of the neuronal dynamics, like oscillations or the reliability of spike patterns.

2 Materials and methods

Neuron model

Simulations of the spiking neurons were performed using a custom version of the NEST simulator (Gewaltig and Diesmann 2007) and the PyNN interface (Davison et al. 2008), with a fixed time step of 0.1ms. In all simulations, we use a planar integrate and fire neuron with exponential non-linearity as introduced in Brette and Gerstner (2005). The dynamics of the membrane potential is composed of a capacitive current \(C_{\mathrm {m}} \frac {dV_{\mathrm {m}}(t)}{dt}\) and a leak current −g L(V m(t)−E L ), with leak conductance g L and leak reversal potential E L . The ratio τ m = C m/g L gives the relaxation time constant for the membrane voltage. Spikes are generated quasi-instantaneously by active conductances rendered by an exponential non-linearity for the current \(\psi (V_{\mathrm {m}}(t))={\Delta }_{T} e^{\frac {V_{\mathrm {m}}(t)-V_{T}}{{\Delta }_{T}}}\). Thus, a spike is initiated when the membrane potential V m goes over V T (with V T > E L ). The spike is cut when the voltage reaches V spike, the membrane potential is then reset to V reset.

An additional slow variable u, with timescale τ u, accounts for the effects of adaptation currents resulting from channels with slow dynamics. The coupling parameter a between V m and u is a linear approximation of hyper-polarizing (a>0) ionic conductances such as those associated with voltage gated potassium channels. Finally, u is increased by an amount b after each spike, which models the effects of highly non-linear conductances such as those associated with calcium gated potassium channels. This results in the following system for (V m, u):
$$\begin{array}{@{}rcl@{}} C_{\mathrm{m}} \frac{dV_{\mathrm{m}}(t)}{dt}&=&-g_{\mathrm{L}} (V_{\mathrm{m}}(t)-E_{\mathrm{L}}) \end{array} $$
$$\begin{array}{@{}rcl@{}} &&+g_{\mathrm{L}}\psi(V_{\mathrm{m}}(t)) +I_{\text{syn}}-u \\ \tau_{\mathrm{u}} \frac{du}{dt}&=&a(V_{\mathrm{m}}(t)-E_{\mathrm{L}})-u \end{array} $$
with the spike condition:
$$\begin{array}{@{}rcl@{}} V_{\mathrm{m}}(t)>V_{\text{spike}}&\rightarrow& V_{\mathrm{m}}(t^{+})=V_{\text{reset}}\\ && u(t^{+})=u(t)+b \end{array} $$
The details of all numerical values for cell properties can be found in Table 1. These values were chosen according to those found in the literature for cortical neurons (Pospischil et al. 2011; Rossant et al. 2011; Hertäg et al. 2014). Parameters are identical for excitatory and inhibitory neurons except when specified, for example adaptation parameters are set to 0 for inhibitory cells. For every simulation, initial values of V m are drawn from a distribution uniform in [V rest, V T].
Table 1

Parameters of the adaptive exponential and fire neuron used in all the simulations

C m

Membrane capacitance

g L

Total leak conductance

τ refrac

Duration of refractory period

Spike initiation threshold

Resting membrane potential

Slope factor

τ m

Membrane time constant

V reset

Reset potential

100 (E) 0 (I)

adaptation conductance

1 (E) 0 (I)

Spike-triggered adaptation current

τ u

Adaptation time constant

E rev

(E) 0 (I) −80

Synaptic reversal potential

V spike

Spike threshold

For parameters specific to a neuron type, E denotes excitatory cells and I denotes inhibitory cells. Note that adaptation parameters a and b are varied widely across simulations

Changes in synaptic conductances triggered by incoming spikes from excitatory and inhibitory neurons are modeled such that the total synaptic current to a neuron can be written as
$$ I_{\text{syn}}(t)=\sum\limits_{s \in \{exc, inh\}} (V(t)-E_{\text{rev}}^{\mathrm{s}})\sum\limits_{k} g_{s}(t-{t_{k}^{s}}) $$
The times \({t_{k}^{s}}\) (s∈{exc, inh}) are the times of the incoming spikes, respectively at excitatory and inhibitory synapses. The dynamics g s∈{exc, inh}(t) after a spike is described by an alpha-function, from Rall synapse model (Bard Ermentrout and Terman 2010), so that the synaptic inputs may be rewritten as the convolution of the spike trains with kernels K s with \(K^{s}(t)=[t]_{+} e^{-\frac {t}{\tau _{s}}}/\tau _{s}\), [.]+ representing the Heaviside function. If the maximal conductance for a synapse type s is written \({g_{s}^{m}}\), we have \(g_{s}(t)={\sum }_{k}{g_{s}^{m}} K^{s}(t-{t_{k}^{s}})\), with k running over incoming spikes and s∈{exc, inh}. Figure 1 is illustrating all those concepts, in a condensed form. In all the paper, we took τ exc = 2ms and τ inh = 3ms for the synaptic time constants.

Diffusion approximation

The spiking activity for a neuron receiving only excitatory inputs at a rate ν exc is controlled by the total effective conductance input \(\mu =g^{m}_{exc}\nu _{\text {exc}}\tau _{\text {exc}}\). The minimum effective conductance input for spike to be triggered, related to the rheobase current, is the one that brings the membrane potential just above its threshold value and it is defined in the following by \(g^{\text {rheo}}=\frac {g_{L}(V_{T}-E_{L})}{V_{T}}\). The deviation of the membrane potential from the diffusion approximation is measured by the Kullback-Leibler divergence between the simulated membrane potential P mc and a theoretical Gaussian distribution P th having mean and variance as predicted from the diffusion approximation:
$$\mathcal{D}(P_{mc}||P_{th})=\int\limits_{-\infty}^{\infty} P_{mc}(v) ln \frac{P_{mc}(v)}{P_{th}(v)} dv. $$

In all the manuscript, we consider the sub-threshold or linear part of the adaptation as the one controlled by the a parameter in the equation of the adaptive exponential neuron, and the supra-threshold or non-linear part of the adaptation as the one controlled by the b parameter. Therefore, a neuron with only linear adaptation is one with b = 0, and one with only non-linear adaptation has a = 0. Physiological interpretation of these parameters is discussed in the neuron model description.

Cortical column

A column is composed of two populations of excitatory and inhibitory neurons connected in a random manner (Erdös-Renyi wiring) with excitatory and inhibitory weights g exc and g inh, and receiving external input ν ext from spike source generated through Poisson processes with weights g ext. This is equivalent to a so-called sparse balanced network (Brunel 2000). Neuron parameters are the same as those described in Table 1. A schematic drawing of the column can be found on Fig. 4a. More precisely, the sizes of the populations are N exc = 800, N inh = 200 and N ext = 200 (a classical 4:1 ratio between excitatory and inhibitory cells). Connections among neurons are drawn randomly with probabilities p AB , A (and B) being populations from which input (and output) neurons are selected (E for excitatory, I for inhibitory, ext for external input). In and out degrees of neurons are thus distributed according to a Poisson law with parameter p AB N A N B . In all the simulations, we use a value of p EE = p EI = 0.05, p II = p IE = 0.05 and p ext = 0.01. Excitatory weights are fixed to g exc = 6nS and inhibitory weights are modified through the control parameter g defined by g inh = g.g exc. For external connections, we have g ext = g exc/6. Delays of the connection are all equal to the simulation time step, i.e. 0.1ms.

Classification of dynamical regimes

The column is considered to be in a Synchronous regime if the pairwise spike correlations 〈CC(0)〉 are over a 0.026 threshold value. Pairwise spike correlations are computed as the mean Pearson correlation coefficient averaged over N = 10000 pairs of randomly chosen cells. The area for saturated regime with Synchronous regular dynamics correspond to an average firing rate over 75 Hz and a mean coefficient of variation for the interspike intervals (CV ISI) less than 0.2. Silent regimes correspond to firing rates lower than 0.2 Hz. Regions of the diagram not detected by these criteria are denoted as the Asynchronous Irregular regime.

Reliability of responses

The reliability of the response is assessed by considering repeated input spike trains from a population of 2000 neurons connected with probability p ext = 0.01 to the excitatory cells of the column. For inputs as homogeneous Poisson process, we consider spike trains with firing rate, r cst = 85Hz and for inputs as inhomogeneous Poisson processes, the firing rate is modulated by a sine function, r var = f 0+f 1 cos(2πωt) with f 0 = 100Hz, f 1 = 40Hz and ω = 5Hz. The spiking response and the average membrane voltage is recorded for 40 times and measures are averaged on these 40 repetitions of the same input spike trains.

The reliability of the spiking response is monitored by the pairwise spike correlation averaged over pairs of repetitions. To assess the reliability of membrane potentials, we monitor the Coefficient of Variation (CV) of the mean voltage V c of excitatory neurons, that is the standard deviation divided by the mean over the 40 repeated trials, and is unitless. The spike and voltage reliability measures are averaged over 10 different inputs and reported with standard deviations for various values of a and b in Fig. 9.

Mean field

The coarse-grained dynamics of a column can be captured by the firing rates (r E , r I ) averaged over excitatory and inhibitory populations. The evolution equation for the firing rate dynamics is then reduced to the Wilson-Cowan system:
$$ \left\{\begin{array}{l} \tau_{E} \frac{dr_{E}}{dt} = -r_{E}+f_{E}(w_{EE} r_{E} \tau_{E} -w_{IE} r_{I} \tau_{I})\\ \tau_{I} \frac{dr_{I}}{dt} = -r_{I}+f_{I}(w_{EI} r_{E} \tau_{E} -w_{II} r_{I} \tau_{I}) \end{array} \right. $$
with w AB standing for the effective coupling from the sub-population A to the sub-population B (\(w_{AB}=p_{AB}N_{A} N_{B} {g_{A}^{m}}\), p AB probability of connection from A to B, N A number of neurons in population A, \({g_{A}^{m}}\) conductance level at a synaptic contact), τ A the synaptic time constant of the sub-population A mentionned in the neuron model and f A the firing rate response of sub-population A. The function f A (I) is the response curve, that is the mean firing rate of neurons in sub-population A when stimulated by incoming inputs I. These response curves and their approximations are discussed in the results Section 3.1. We are using standard numerical methods of bifurcation theory as implemented in the AUTO library with the python interface PyDSTools (Clewley 2012). The local stability analysis for the fixed points of the dynamical system consists in the study of the set of parameters where some eigenvalues of the Jacobian are zero. Further study of the Jacobian are indicative of the type of instability at these points (Kuznetsov 1998). A fold (or saddle-node) bifurcation curve separates regions with one stable fixed point from bistable regions with two stable fixed points separated by an unstable fixed point. Along those curves, the Jacobian has a zero eigenvalue and a stable fixed point and an unstable fixed point annihilate. When an eigenvalue is non-zero but has a zero real part, there is a Hopf bifurcation curve separating a region with one stable fixed point from a region with a limit cycle. On those curves, eigenvalues of the Jacobian are a pair of pure conjugated numbers. If 2 saddle-node branches collide in the parameter plane, the resulting bifurcation point is a so-called cusp. Other bifurcations points are related to Hopf bifurcations when the linear field cancels (Bogdanov-Takens point) or when the quadratic contribution to the vector field cancels (Bautin point). Note that Bautin points separate the part of the Hopf curve where the limit cycle arising has infinitesimal amplitude from the part where it has finite amplitude. Simulations of the stochastic dynamics where implemented in C++ using the Milstein scheme.

3 Results

3.1 Dynamics of a single neuron without adaptation

First, we consider the case of a single adaptive exponential integrate-and-fire neuron neuron simply bombarded with excitatory and inhibitory homogeneous Poisson inputs. Both types of inputs are triggering conductances changes at the soma level as illustrated on Fig. 1a, and the spiking activity of the neuron is modified by the adaptation mechanisms regulating the activity. In all the following, we assume that “linear” (or sub-threshold) adaptation is mediated by the parameter a of the model (see Section 2), while “non-linear” (or supra-threshold) adaptation is related to the b parameter. To gain a better understanding of the responses, we consider the spiking activity of the neuron without any adaptation (a = 0, b = 0, see Section 2). Without a loss of generality, we can restrict the analysis to the case where only excitatory inputs are impinging the cell (ν inh = 0). In this case, the average synaptic input onto the cell for filtered Poisson spike train of rate ν exc with ν exc τ exc>>1 can simply be approximated by a Gaussian process of distribution G exc = N(μ, σ 2) with a mean μ and a variance σ 2 (Papoulis 1965; Ricciardi and Sacerdote 1979; Lánský and Lánská 1987; Richardson and Gerstner 2005; Hertäg et al. 2014) (see Section 2 for the synapse model):
$$ \left\{\begin{array}{rcr} \mu &=& \nu_{\text{exc}} {\int}_{0}^{\infty} K^{\text{exc}}(\tau)d\tau = g^{m}_{\text{exc}} \tau_{\text{exc}} \nu_{\text{exc}}\\ \sigma^{2} &=& \nu_{\text{exc}} {\int}_{0}^{\infty} (K^{\text{exc}}(\tau))^{2}d\tau = (g^{m}_{\text{exc}})^{2} \tau_{\text{exc}} \nu_{\text{exc}} \end{array}\right. $$
Fig. 1

Illustrationof synaptic current and diffusion approximation. a. Schematic view of a neuron with the various components included in the model. Excitatory and Inhibitory inputs are modelled as α-shaped conductances (inset), and examples of the V m dynamics can be seen in different cases: without adaptation or with non-linear adaptation. b.Plain line: KL-divergence between input simulated conductance distribution and theoretical estimate \(\mathcal {N}(\mu , \sigma ^{2})\) in the diffusion limit (ν exc> > τ exc) for various synaptic time constants, assuming independent Poisson inputs at excitatory synapses. Dashed line KL-divergence between the recorded membrane potential distribution and its closest Gaussian distribution with same mean and variance. Insets are example of simulated conductance traces and recorded V m distribution (in mV), compared to theoretical predictions

The validity of such a diffusion approximation can be checked by computing the Kullback-Leibler divergence (KL, see Section 2) between the theoretical distribution of the synaptic input conductances (G exc) and numerical simulations for various values of the maximal conductance \(g^{m}_{\text {exc}}\) and input firing rate ν exc, keeping the total mean effective input \(\mu =g^{m}_{\text {exc}} \tau _{\text {exc}} \nu _{\text {exc}}\) constant. In all the numerical simulations, we use the planar adaptive integrate and fire neuron described in Brette and Gerstner (2005) and further studied in Touboul and Brette (2008).

As we can see from Fig. 1b, in absence of adaptation and for our particular model, this diffusion approximation on G exc becomes valid for a synaptic time constant of τ exc = 2 ms., and for input rates of approximately 1000Hz (KL ≃ 0). This can be seen in the insets by the convergence of the conductance distribution toward a Gaussian distribution, while incoming synaptic conductances are shown as black traces. Similarly, the KL-divergence between the empirically measured distribution of the V m and a Gaussian distribution having mean μ V and variance \({\sigma _{V}^{2}}\) of the real V m distribution (dashed line in Fig. 1b) is converging toward zero. The voltage being filtered with a slower timescale (in the order of τ m≈10 ms, see Section 2), converges faster than for the conductance distribution. Therefore, this diffusion approximation is valid for a neuron receiving input spikes at realistic rate (1–10Hz) from a population as small a 1000 neurons. A reasonable size of a population from which a mean field model can be used is needed to describe its dynamics and to justify the size of the cortical columns we will use in the following.

Then, it is known that for a stationary Poisson input at a firing rate ν exc in the diffusion limit, i.e for low values of \(g^{m}_{\text {exc}}\) and approximated by a Gaussian process \(\mathcal {N}(\mu , \sigma ^{2})\), individual fluctuations in V m are small compared to V T . Therefore, the synaptic current \(I_{\text {exc}}(t)=G_{\text {exc}}(t)(V_{\mathrm {m}}(t)-E^{exc}_{\text {rev}})\) can also be approximated by a Gaussian process with mean \(\hat {\mu } = \mu (V_{T}-E^{exc}_{\text {rev}})\) and variance \(\hat {\sigma }^{2}=\sigma ^{2} (V_{T}-E^{exc}_{\text {rev}})^{2}\) (Destexhe et al. 2001; Richardson and Gerstner 2005). The dynamics for the exponential integrate-and-fire neuron can be summarized by the adaptive integrate-and-fire Ornstein-Uhlenbeck process defined as follows
$$\begin{array}{@{}rcl@{}} \tau_{\mathrm{m}} \frac{dV_{\mathrm{m}}(t)}{dt} &=& -(V_{\mathrm{m}}(t)-E_{\mathrm{L}})+ \psi(V_{\mathrm{m}}(t)) \end{array} $$
$$\begin{array}{@{}rcl@{}} && +\frac{\hat{\mu}+\hat{\sigma}\sqrt{\tau_{m}}\eta(t)}{g_{\mathrm{L}}} \end{array} $$
where η(t) is drawn from a Gaussian distribution \(\mathcal {N}(0, 1)\).
For this model, as already studied in La Camera et al. (2004) the rheobase effective input g rheo, i.e. the minimum conductance input triggering a spike, is a fixed value that can be compared with estimations obtained from numerical simulations. As we can see on Fig. 2a, for low input firing rate ν exc, the neuron tends to fire spikes at lower values h rheo of input threshold than the theoretical one but as the input firing rate is increased, the fluctuations in synaptic inputs vanish and the neuron tends to behave as its deterministic limit g rheo. From Fig. 2a, we conclude that this deterministic limit, reached at around 5000Hz, is more restrictive than the one given by the diffusion approximation alone, but remains valid for the networks considered in the following.
Fig. 2

a. Difference between the real rheobase conductance g rheo and its estimated value h rheo from diffusion approximation as a function of the rate of excitatory Poissonian input. b. Firing rate of the neuron without adaptation subject to excitatory inputs from a Poisson source as a function of the input μ, for several maximal conductances \(g^{max}_{exc}\). Red curve is for low membrane conductances values and high firing rates, while green curve is obtained for high conductances and low firing rates. Plain lines are theoretical approximations obtained for the situations (low rate, high conductance) (green) and (high rate, low conductance) (red), dots are obtained from numerical simulations. Inset shows simulations of couples (rate, conductance) for several values of g m ., color-coded

The two limiting cases of large and small conductance g exc can be studied while maintaining the mean effective input μ = g exc ν exc τ exc constant. When the maximal conductance is close to \(g^{max}_{\text {exc}}\), response is linear at low inputs rates and saturates to the maximal firing rate imposed by the refractory period, with \(\nu _{\text {out}}=\frac {\nu _{\text {exc}}}{1 + \tau _{\text {ref}} \nu _{\text {exc}}}\) (see green curve on Fig. 2b). On the opposite, in the limit of high firing rates, ν exc≫1/τ exc, the sub-threshold dynamics is well approximated by the Ornstein-Uhlenbeck process shown above, and it is known Ricciardi and Sacerdote (1979) and Sacerdote and Giraudo (2013) that the mean first passage time is given by
$$T_{OU} = \sqrt{\frac{\pi \tau_{\mathrm{m}}}{\hat{\sigma}^{2}}} {\int}_{V_{\text{reset}}-\hat{\mu} \tau_{\mathrm{m}}}^{V_{\text{spike}}-\hat{\mu} \tau_{\mathrm{m}}} \{1+Erf(z/\hat{\sigma} \sqrt{\tau_{\mathrm{m}}})\} e^{\frac{z^{2}}{\hat{\sigma}^{2}}}dz $$
where Erf is the error function and other parameters are described in the Methods (see Richardson and Gerstner 2005, Ostojic and Brunel 2011 for a fast method to compute this quantity). The firing rate for the adaptive exponential neuron without adaptation is then simply ν out = 1/(τ ref +T OU ) (see red curve on Fig. 2b). Note that in the inset of Fig. 2b firing rate response are plotted as a function of g exc ranging from 0.05 to 50 nS, color-coded. Response behaves linearly for \(g_{exc} \simeq g^{max}_{\text {exc}}\). Interestingly, for values of g exc in between, response curves cross almost all in the same region around which the response firing rate is not dependent on the size of post-synaptic potential for a given effective input. This could be interesting for networks including heterogeneities of synaptic weights as their output rate would be similar for identical effective input when input rates are scaled appropriately.

3.2 Dynamics of a single neuron with adaptation

Now if we consider adaptation, its dynamics (τ u = 100ms) is slower than timescales involved in spike generation or membrane relaxation (τ m = 10ms) so for the analysis of response properties it is possible to assume timescale separation (van Vreeswijk and Hansel 2001; Benda and Herz 2003; Ermentrout 1998). Let us first study the effect of linear (sub-threshold) adaptation current on the response to Poisson input spike train. At slow time scale, the fast membrane dynamics can be averaged so that \(\langle (V_{\mathrm {m}}(t)-E_{L})\rangle =\frac {\hat {\mu }}{g_{L}}\) and injecting this value in the dynamics for the adaptation current gives the stationary current \(u=a \frac {\hat {\mu }}{g_{L}}\) resulting in the effective dynamics for V m(t)
$$\begin{array}{@{}rcl@{}} \tau_{m} \frac{dV_{\mathrm{m}}(t)}{dt}&=&-g_{\mathrm{L}} (V_{\mathrm{m}}(t)-E_{\mathrm{L}})+g_{L} \psi(V_{m}) \\ &&+(1-a)\frac{\hat{\mu}}{g_{L}} + \frac{\hat{\sigma}}{g_{L}} \eta(t) \end{array} $$

The effect of linear adaptation is thus to reduce the effective input and the firing rate response for all inputs, shifting the response curve. This negative feedback could be useful to tune the threshold depending on the basal level of computation, and a similar mechanism has been shown to enable contrast invariant computations in the visual system (Carandini and Ferster 1997).

If we now consider the non-linear effect of adaptation, we can notice that between two consecutive spikes t k+1t k = T, a discrete map describes the dynamics of the adaptation current. At each spike time, we have \(u(t) \rightarrow u(t)e^{-T/\tau _{u}} + b\), so that between 0 and T we have \(u(t)=b(1-e^{-t/\tau _{u}})/(1-e^{-T/\tau _{u}})\). For a high firing rate ν exc compared to the adaptation time scale τ u , the average adaptation current can therefore be estimated as 〈u〉=τ u exc . Figure 2b (red curve) showed that for ν exc far below saturation, the response curves of the neuron without adaptation behave linearly. In this linear regime of the response, the firing rate can therefore be described with a linear relationship \(\nu _{\text {out}}=k(\hat {\mu }-\hat {\mu }_{c})\), with \(k, \hat {\mu }_{c}\) being constants and \(\hat {\mu }\) the total current to the membrane \(\hat {\mu }=I_{syn}-u\). By replacing u(t) by its mean value estimated above, we have \(\nu _{\text {out}}= k (I_{syn}-\hat {\mu }_{c})/(1+kb\tau _{u})\), such that the net effect of non-linear adaptation is to change the slope of the response to slow input variations, as also found in Ermentrout (1998) and Benda and Herz (2003).

Figure 3a shows examples of the dynamics of the exact same neuron, in different conditions of adaptation. We compare single-cell responses to a particular input without any adaptation (a = 0, b = 0) to a regime with full adaptation (a, b with standard values as defined in Section 2), one with only linear adaptation (b = 0), and one with only non-linear spike frequency adaptation (a = 0). In Fig. 3b, the firing rate responses as functions of the input μ are plotted for those different adaptation conditions. The effect of linear adaptation is a shift in the spike threshold (see Fig. 3b, the fact that when b = 0 (green curve), the curve is rising for higher values of μ compared to the response curve with no adaptation (in red)). Similarly, the effect of non-linear adaptation (a = 0) is equivalent to a change in the slope of the response curve to slow input variations (see blue curve). Schematically, we can fit with a least square fitting procedure those response curves to sigmoidal functions \(f(\nu _{\text {ext}})=1/(1+e^{-\alpha (\nu _{\text {ext}} -\beta )})\) (see Fig. 3c, top) and quantify this behavior.
Fig. 3

Influence of linear and non linear adaptation a. Comparison of the response to the same excitatory input from a Poisson spiking process for a neuron without adaptation and neuron with (up) linear adaptation (middle) non linear adaptation (bottom) combined linear and non linear adaptation conditions (see Section 2). b. Firing rate response as a function of the input rate in the 4 conditions previously listed. c.Top: Sigmoidal rate function used to estimate the effect of adaptation. Middle: Effect of linear adaptation on threshold (β, unfilled square markers) and gain (α, filled square markers) of the sigmoid. Bottom: same as middle, but when the non-linear adaptation is varied

Using the sigmoidal fits, we study the variation of the thresholds β and the gains α, as functions of the adaptation parameters of the neuron, a and b. We can see on Fig. 3c, middle, that the linear adaptation (controlled by a) affects mostly the threshold β shifting the response curve to the right with higher values of a resulting in higher thresholds β. This dependence is linear as expected from the analysis via time scales separation. On the other hand, the non-linear adaptation (controlled by the parameter b) affects mostly the gain α of the sigmoidal response curve by decreasing the slope of the response function, with higher values of b resulting in lower α (see Fig. 3c, bottom). Again, as expected from former considerations, this dependence is non-linear. The decrease in gain also results in a wider dynamic range. Note that the change in α is necessary to keep the rheobase constant so that non-linear adaptation changes both α and β. We note that having lower gain, the input range for which linear approximation of the response is valid is increased and this can be summarized by stating that the net effect of non-linear adaptation is to linearize the response curve.

3.3 Dynamics of a cortical column with adaptation

In the neocortex, neurons are arranged in complex microcircuits affecting their response properties and giving rise to internal dynamics. The cortical column (Horton and Adams 2005) is a good example of such canonical circuits encountered in the brain and has been shown to give rise to oscillatory rhythms and self-sustained irregular activity (Brunel and Hakim 1999; Vogels et al. 2005). If the fine connectivity details of such a column are still poorly understood (Binzegger et al. 2007), we can still understand the dynamical properties of a homogeneous balanced network of adapting neurons, arranged in a columnar fashion. To this aim, we simulated what we called a cortical column composed of two populations of 800 excitatory and 200 inhibitory neurons connected in a random manner with an averaged probability of connection of 5 % (Renart et al. 2010; Brunel 2000), and receiving external input from spike source generated through Poisson processes. Such a generic and classical network is often referred to in the literature as a random balanced network (Brunel 2000; Vreeswijk et al. 1996), and has been well used as a good model of in vivo activity in sensory cortices. More details could be found in the corresponding section of the Section 2. Schematically, the structure of the network is represented in Fig. 4a. Excitatory weights are fixed to g exc = 6nS and inhibitory weights are modified through the control parameter g such that g inh = g.g exc. Note that on Fig. 4a, insets shows the distribution of the indegree for all the different connections.
Fig. 4

Structure and dynamics of the cortical column model a. Schematics of the wiring connections in the cortical column. Two populations of 800 excitatory (red) and 200 inhibitory (blue) neurons are reciprocally connected with a probability of 5 %. Excitatory cells are also stimulated by external inputs at a rate ν ext. Insets are showing the distribution of the indegree for all the different connections between populations involved in the column. b. Phase diagram as a function of external firing rate ν ext and inhibitory maximal conductance (when g is varied), without any adaptation. We can see four distinct regimes of activity. c. Typical spike rasters for the three non-silent regions of the phase diagram

In Fig. 4b, we show the phase diagram of such a cortical column without any adaptation when the external rate ν ext and the balance g between excitation and inhibition are varied. We can see that the dynamics of the cortical column can be maintained in an Asynchronous Irregular regime (AI, Section 2 for details on the classification) with low firing rate when external input is small enough, as already found (Brunel 2000). This large region is termed Asynchronous Irregular because neurons are firing irregularly (high coefficient of variation of their inter spikes intervals), and the network’s firing rate displays no clear oscillations. Similarly, we can see in the diagram a Synchronous Regular (SR) or a Synchronous Irregular (SI) regime. It is commonly assumed that the AI regime is a good candidate for describing cortical dynamics observed in vivo and it has been identified as the operating point of the brain with transient perturbation leading to amplified response with fast recovery to the resting operating point (Renart et al. 2010; Stimberg et al. 2009). Canonical raster plots of the three main dynamical regimes identified in the phase diagram are shown in Fig. 4c. In such AI regime, when input is increased, the dynamics of the network depends on the dominant polarity of the synaptic current. For an excitation dominated column, increasing input destabilizes the asynchronous irregular state so activity is amplified to a persistent saturated state with all neurons firing at their maximal frequency. For an inhibition dominated network, the destabilization of the irregular state through increased stimulation leads to oscillatory dynamics with neurons firing synchronously.

We can see what is the direct effect of the two previously discussed components of adaptation (linear v.s non-linear) in Fig. 5. As we can notice in Fig. 5a, when linear adaptation is turned on, the Synchronous Irregular region of the phase diagram tends to disappear, while the overall structure of the diagram is preserved in case of non-linear adaptation. In both cases, because the average firing rates in those regions is similar (data not shown), it is really the influence of single-cell response curves that is affecting the phase diagram. Linear adaptation, responsible for a shift in the response threshold is decreasing the synchronous activity, thus reducing the amount of correlations in the network.
Fig. 5

Dynamics of the cortical column model with adaptation. a. Phase diagram as a function ν ext and g, and raster plots of the three non-silent regions with linear adaptation only b. Same as in a., but with only non-linear adaptation

3.4 Mean field description of the cortical column

To get a better understanding on how adaptation is impacting the phase diagram, we used the fact that the coarse-grained dynamics of the cortical column can be captured, at a macroscopic level, by describing the firing rates (r E , r I ) averaged over homogeneous excitatory and inhibitory neurons. Using a mean field approach, as in Tabak et al. (2000) and Giugliano et al. (2004), the evolution equation for the firing rate dynamics can be reduced to the Wilson-Cowan system:
$$\begin{array}{@{}rcl@{}} \tau_{E} \frac{dr_{E}}{dt}&=&-r_{E}+f_{E}(w_{EE} r_{E} \tau_{E} -w_{IE} r_{I} \tau_{I}) \\ \tau_{I} \frac{dr_{I}}{dt}&=&-r_{I}+f_{I}(w_{EI} r_{E} \tau_{E} -w_{II} r_{I} \tau_{I}) \end{array} $$
For more details, see Section 2. To describe the dynamics, we keep the sigmoidal approximation discussed in the previous section for the response curve \(f_{A}(x)=(1+e^{-\alpha (\nu _{\mathrm {A}} -\beta )})^{-1}\). Beyond a fixed point attractor, this system is known to exhibit multistability and limit cycle. The transition among these regimes can be studied through a bifurcation analysis (see Section 2 for details about the related numerical tools). For this matter, as linear and non-linear adaptation were shown to affect the threshold β and the gain α of the response curves, these parameters are taken as bifurcation parameters in Fig. 6 to study the various possible instabilities.
Fig. 6

Mean Field dynamic of the cortical column with adaptation a. Diagram showing bifurcations of codimension 1 and 2 as a function of threshold β and gain α of the excitatory response function. The stationary distribution for the activity of the excitatory population under weak stochastic perturbation is also depicted at various points close to the instability curve. b. Dynamics, as a function of time, of the excitatory population at various points in the phase space

Figure 6 depicts the distributions and the temporal traces of the activity for the excitatory sub population in the system under stochastic perturbation close to various bifurcation points (see Methods for a description of the numerical tools used for bifurcation analysis). Close to saddle-node points, the distribution is bimodal with a bias toward up or down state depending on which side of the bistable region is considered (low or high α). Similar bistability is observed between a limit cycle and a fixed point of high, low or moderate level of activity close to the Hopf curve. A prediction of the model is that if sufficient linear adaptation, affecting the threshold parameter β in the coarse-grained model, is at work in the excitatory cell dynamics, a column which would be initially oscillatory or bistable would be silenced to a low activity level in an abrupt fashion. The effect of non-linear adaptation, affecting the gain parameter α, is more subtle with possible transitions along the vertical axis from bistable regime in the upper cusp to various oscillatory regimes, changes in the oscillatory pattern or transition from oscillation to fixed point.

From Fig. 6 we see that the crossing of the Hopf curve is reached with a reasonable change of the threshold value β (related to linear adaptation) whereas a large change in gain α (related to nonlinear adaptation) is necessary to escape from the oscillatory region. This is consistent with previous simulations of spiking neurons where we observed that linear adaptation reduces greatly the synchronous irregular region of the diagram (see Fig. 5).

3.5 Dynamical response to external stimuli with or without adaptation

3.5.1 Oscillatory dynamics

To study the effect of adaptation on columnar dynamics modeled with spiking neurons (see Section 2), the network conductances were set so that the mean firing rate of the network is oscillatory when no adaptation is considered (a = 0, b = 0) and adaptation parameters were then varied from this setting. The measures of the activity of the excitatory sub population (firing rate, coefficient of variation of interspike) depending on those parameters are shown on the phase diagrams of Fig. 7a, b, c when the two components a and b of the adaptation are varied. As expected, when linear or non-linear adaptation is increased the mean firing rate activity is decreased but the effects of adaptation is different on the second order statistics for the dynamics at network level. The pairwise spike correlations in the network drop abruptly when linear adaptation is varied (see panel c) whereas much smoother transition is observed when non-linear adaptation is varied. Furthermore, when b is increased, the average firing rate in the column is oscillating at half the period of the original network. We can see on Fig. 7d, e, f examples of the activity in the three canonical regimes: without adaptation, with only linear adaptation (b = 0), or with only non-linear adaptation (a = 0).
Fig. 7

Influence of adaptation on the columnar dynamics. a. Measures of the network activity depending on linear and non linear adaptation parameters (a, b). Lower point in the network is set to the parameters of purple sample from Fig. 4 d, e, f. Samples of network activity (spike raster and mean rate) in no adaptation, linear adaptation and non linear adaptation conditions

For a better grasp at the effects of adaptation on the dynamics of a column, we consider the time course of the activity when a stimulus is presented starting from quiescent state. The resulting responses are reported on Fig. 8. As expected from time averaged measures, after the stimulus was presented, the mean firing rate of the column oscillates when no adaptation is considered. Strong adaptation currents result in relaxation to a stationary firing rate. This relaxation occurs at the slow time scale imposed by the adaptation current for linear adaptation whereas it occurs much faster, within an oscillation cycle, for non-linear adaptation currents. To summarize, with subthreshold adaptation a, oscillatory fluctuations first decrease in amplitudes and then for higher values of a the column firing rate also decreases. In contrary, with suprathreshold adaptation b, both oscillatory fluctuations, oscillation frequency and mean firing rate of the columns are decreasing simultaneously.
Fig. 8

Columnar dynamics under sustained Poissonian input at constant rate. a. Time course of the mean firing rate of of cortical column for various levels of adaptation currents coded in colors when linear adaptation only is varied (b = 0). The shaded areas represent two standard deviation around the mean. b. Firing rate for various levels of adaptation currents: solid line is the stationary mean and shaded areas represent two standard deviations. The dotted line is the peak response right after the stimulus is presented. c., d. Same as a., b. but when only non-linear adaptation is varied (a = 0)

We showed through theoretical arguments that effects of adaptation on the neuronal dynamics may be related to α and β parameters in the macroscopic model under the diffusion approximation and slow adaptation. The study of the 2 variables macroscopic model is then useful to interpret some aspects of the dynamics like oscillations. There are aspects of changes induced by adaptation which are not discussed in this way because they would go beyond the 2 variables mean field model we consider to discuss qualitative dynamics. For instance chaotic dynamics where found in a macroscopic model including short term plasticity in Cortes et al. (2013).

3.5.2 Reliability of spike patterns

In previous sections, we analyzed the stationary responses of a neuronal network with adaptation, and we distinguished the linear v.s the non-linear part of that mechanism, either with simulations of spiking neurons or with a macroscopic approximation. Now we consider the dynamical responses of the same cortical column, but when stimulated with a Poisson process at constant rate, or with one with rate modulated in a sinusoidal manner (see Methods). For both scenarios, we observed the effect of adaptation on the reproducibility of the responses, when the exact same realization of the input spike patterns were replayed to the exact same network. At each trial, initial conditions of the network were different (see Methods), and to visualize the mean response of the column, we computed and plotted the mean voltage over all N neurons within the column, l(t)=<V m (t)> N .

On Fig. 9a, b, we can see the temporal dynamics of the mean membrane potential, averaged over 40 repetitions, and when stimulated with constant Poisson input (Fig. 9a) or time-varying inhomogeneous Poisson input (Fig. 9b) (see Methods for a description of the input stimulus). In both constant and varying rate cases, the variance over trials of the responses are reduced by the adaptation currents. This can be viewed more clearly in the raster plots below Fig. 9a and b, showing the trial-to-trial spiking responses of 3 sampled neurons taken in the excitatory population of the column for the 40 repetitions of the same input.
Fig. 9

Reliability of columnar dynamics under repeated stimulus. We consider neurons with no adaptation current in black, with only linear (a = 500nS) in red and with only nonlinear adaptation current (b = 1500nA). a Mean membrane potential ± std averaged over repetitions to a constant Poisson input at 850Hz, without any adaptation (black), with linear only (red), or non-linear only (green). Below are the trial-to-trial raster plots of three representative cells in the column each spiking over 40 repetitions, in those three conditions. b. Same as in a. but in response to an Poisson process with rate varying as a sinus during a period of 200 ms (see Methods). c. Measures of the Coefficient of Variation (CV) of the mean membrane potential in the period for the 40 repetitions of the stimulus averaged across 10 realizations of the network. The measures was computed for constant (solid) or varying (dashed) input rate as a (red) or b (green) were varied. d. Normalized pairwise cross-correlations between repeated response of a neuron, averaged across 100 non-silent cells and over 10 realizations. Line and color codes are same as in c

To quantify this increase in reliability between trials, we used two measures. The first one is the CV of the mean voltage averaged over all N e neurons within the column. We can see on Fig. 9c that this measure is affected similarly by both adaptation currents in the case of the stimulation at constant rate (dashed lines), while it has a minimum at a finite value of nonlinear adaptation current b in the case of a temporally varying stimulation (solid green line). For high values of a, the CV is almost 0, meaning that the reliability is very high. To rule out the fact that this will depend on the average activity, we also computed the normalized averaged cross-correlation coefficient within trials and among neurons spike trains as described in Methods (see Fig. 9d). On Fig. 9d, we can see similar trends as for the average membrane potential: the normalized correlation coefficient is close to 1 for high values of linear adaptation (red curves) and there is an optimal value of b at which maximal reliability is achieved with time-varying input rate.

4 Discussion

In this paper, we analyzed the effects of adaptation in the adaptive exponential integrate-and-fire neuron model, distinguishing linear and non-linear mechanisms, related biologically to voltage-gated and calcium-gated channels. Using mathematical observations and comparisons with simulations, we showed how those two distinct components could lead to different changes of the firing rate response curve for this particular neuron model. In single-cell simulations, we found that linear adaptation affects mostly the threshold at which a neuron starts to fire, while non-linear adaptation tends to lower the slope of the response curve. Extending those observations to cortical networks, we studied the role of adaptation onto the dynamics of a cortical column (a so-called random balanced network) with an activity similar to what can be observed in vivo. We found that linear adaptation introduces a switch from cortical oscillations to a fixed point of stationary low firing rate while non-linear adaptation preserves cortical oscillations but shifts their frequency to lower values. These are thus two possible ways to modulate synchrony in a neuronal network.

From a functional point of view, we showed that the dynamic range of the response of a cortical column is increased when non linear adaptation currents are included and that there is a shift of the threshold and response to lower values when linear adaptation current is included. These aspects may be important to understand the contrast invariant computations at the neuronal level and predictive aspects of the neuronal response (Deneve 2008; Boerlin et al. 2013). We also demonstrate that when oscillatory dynamics arise in a cortical column, both adaptation currents attenuate these oscillations but in different manners. Our study thus sheds lights on the possible ways to modulate oscillatory dynamics with slow currents which is of crucial interest for cognitive neuroscience (Buzsáki 2010). Moreover, we found that adaptation was also able to greatly increase the reliability of the neuronal responses and as observed in vitro (Mainen and Sejnowski 1995) and in vivo (Haider et al. 2010), responses are reliable on a trial-to-trial basis. The linear part of the adaptation mechanism used in the integrate-and-fire model considered here was the one mostly involved in this process. Because this linear adaptation is related to voltage-gated channels, we can establish a link between the temporal precision of the responses and the distribution of those channels in neurons. Interestingly, it is known that the repartition of such channels can also be activity-dependent (Lu et al. 2004), so neurons may have mechanisms to adapt their precision as a function of their inputs.

Adaptation is an ubiquitous phenomena in the brain that can spawn multiple time scales: from time constants of several minutes or even hours, it has been shown to be crucial for homeostasis and stability in neuronal network. In this work, we investigated only adaptation mechanisms in the order of hundreds of milliseconds, so relatively close to the membrane time constant. This form of adaptation is known to rely on the kinetics of the voltage or calcium gated channels, and for a more in depth knowledge, see Marder (2012). Therefore, our study focused only on dynamical responses and equilibrium that could be reached within those time scales. Results are observed for random balanced networks with sparse connectivity (Vogels et al. 2005), and it is likely that differences may be found in denser networks, or when the balance between excitatory and inhibitory conductances is reduced (Vreeswijk et al. 1996; Ostojic 2014). All together, unravelling the effect of adaptation in spiking networks is crucial to understand the computation that can be performed by such dynamical systems. Functionally, one could hypothesize that linear (sub-threshold) adaptation, responsible for a shift in the threshold of the response curve is useful to implement contrast invariant responses, while the change in the slope of the response curve induced by non-linear (supra-threshold) adaptation increases the dynamic range of the neuron and softly modulates the oscillatory dynamics. This findings remains to be tested in vivo.

This work is supported by the Japan Society for the Promotion of Science. This work is also supported in part by the Platform for Dynamic Approaches to Living System from the Ministry of Education, Culture, Sports, Science, and Technology of Japan. P.Y was supported by the Fondation pour la Recherche Médicale (FRM).

Conflict of interests

The authors declare that they have no conflict of interest.

  1. Andrade, R. (2011). Serotonergic regulation of neuronal excitability in the prefrontal cortex. Neuropharmacology, 61(3), 382–6.PubMedCentralCrossRefPubMedGoogle Scholar
  2. Benda, J., & Herz, A.V.M. (2003). A universal model for spike-frequency adaptation. Neural Computation, 15(11), 2523–2564.CrossRefPubMedGoogle Scholar
  3. Benda, J., Maler, L., & Longtin, A. (2010). Linear versus nonlinear signal transmission in neuron models with adaptation currents or dynamic thresholds. Journal of Neurophysiology, 104, 2806–2820.CrossRefPubMedGoogle Scholar
  4. Benucci, A., Saleem, A.B., & Carandini, M. (2013). Adaptation maintains population homeostasis in primary visual cortex. Nature Neuroscience, 16(6), 724–9.PubMedCentralCrossRefPubMedGoogle Scholar
  5. Binzegger, T., Douglas, R.J., & Martin, K.A.C. (2007). Stereotypical bouton clustering of individual neurons in cat primary visual cortex. Journal of Neuroscience, 27(45), 12242–54.CrossRefPubMedGoogle Scholar
  6. Boerlin, M., Machens, C.K., & Denève, S. (2013). Predictive coding of dynamical variables in balanced spiking networks. PLoS Computational Biology, 9.Google Scholar
  7. Brette, R., & Gerstner, W. (2005). Adaptive exponential integrate-and-fire model as an effective description of neuronal activity. Journal of Neurophysiology, 94(5), 3637–42.CrossRefPubMedGoogle Scholar
  8. Brunel, N, & Hakim, V (1999). Fast global oscillations in networks of integrate-and-fire neurons with low firing rates. Neural Computation, 11(7), 1621–1671.CrossRefPubMedGoogle Scholar
  9. Brunel, N. (2000). Dynamics of sparsely connected networks of excitatory and inhibitory spiking neurons. Journal of Computational Neuroscience, 8(8), 183–208.CrossRefPubMedGoogle Scholar
  10. Buzsáki, G. (2010). Neural syntax: cell assemblies, synapsembles, and readers. Neuron, 4;68(3), 362–85.CrossRefGoogle Scholar
  11. Carandini, M., & Ferster, D. (1997). A tonic hyperpolarization underlying contrast adaptation in cat visual cortex. Science, 276(5314), 949–52.Google Scholar
  12. Clewley, R. (2012). Hybrid models and biological model reduction with PyDSTool. PLoS Computational Biology, 8.Google Scholar
  13. Cortes, J.M., Desroches, M., Rodrigues, S., Veltz, R., a Muñoz, M., & Sejnowski, T.J. (2013). Short-term synaptic plasticity in the deterministic Tsodyks-Markram model leads to unpredictable network dynamics. Proceedings of the National Academy of Sciences of the United States of America, 110, 16610–5.PubMedCentralCrossRefPubMedGoogle Scholar
  14. Davison, A., Brüderle, D., Eppler, J., Kremkow, J., Muller, E., Pecevski, D., Perrinet, L., & Yger, P. (2008). PyNN: a common interface for neuronal network simulators. Frontiers in Neuroinformatics, 2, 11.Google Scholar
  15. Deemyad, T., Kroeger, J., & Chacron, M.J. (2012). Sub- and suprathreshold adaptation currents have opposite effects on frequency tuning. The Journal of Physiology, 590, 4839–4858.PubMedCentralCrossRefPubMedGoogle Scholar
  16. Deneve, S. (2008). Bayesian spiking neurons I: inference. Neural Computation, 20, 91–117.CrossRefPubMedGoogle Scholar
  17. Destexhe, A., Rudolph, M., Fellous, J.M., & Sejnowski, T.J. (2001). Fluctuating synaptic conductances recreate in vivo-like activity in neocortical neurons. Neuroscience, 107(1), 13–24.PubMedCentralCrossRefPubMedGoogle Scholar
  18. Ermentrout, B. (1998). Linearization of F-I curves by adaptation. Neural Computation, 10, 1721–1729.CrossRefPubMedGoogle Scholar
  19. Ermentrout, B., Pascal, M., & Gutkin, B. (2001). The effects of spike frequency adaptation and negative feedback on the synchronization of neural oscillators. Neural Computation, 13, 1285–1310.CrossRefPubMedGoogle Scholar
  20. Bard Ermentrout, G., & Terman, D.H. (2010). Mathematical foundations of neuroscience. Springer.Google Scholar
  21. Farkhooi, F., Muller, E., & Nawrot, M.P. (2011). Adaptation reduces variability of the neuronal population code. Physical Review E, 83.Google Scholar
  22. Gewaltig, M.-O., & Diesmann, M. (2007). NEST (NEural Simulation Tool).Google Scholar
  23. Gigante, G., Del Giudice, P., & Mattia, M. (2007). Frequency-dependent response properties of adapting spiking neurons. Mathematical Biosciences, 207, 336–351.CrossRefPubMedGoogle Scholar
  24. Gigante, G., Mattia, M., & Del Giudice, P. (2007). Diverse population-bursting modes of adapting spiking neurons. Physical Review Letters, 98.Google Scholar
  25. Giugliano, M., Darbon, P., Arsiero, M., Lüscher, H.-R., & Streit, J. (2004). Single-neuron discharge properties and network activity in dissociated cultures of neocortex. Journal of Neurophysiology, 92, 977–996.CrossRefPubMedGoogle Scholar
  26. Grillner, S. (1997). Ion channels and locomotion. Science, 278(5340), 1087–1088.CrossRefPubMedGoogle Scholar
  27. Gusnard, D.A., & Raichle, M.E. (2001). Searching for a baseline: functional imaging and the resting human brain. Nature Reviews Neuroscience, 2(10), 685–94.CrossRefPubMedGoogle Scholar
  28. Haider, B., Krause, M.R., Duque, A., Yu, Y., Touryan, J., Mazer, J.A., & McCormick, D.A. (2010). Synaptic and network mechanisms of sparse and reliable visual cortical activity during nonclassical receptive field stimulation. Neuron, 65(1), 107–21.PubMedCentralCrossRefPubMedGoogle Scholar
  29. Hertäg, L., Durstewitz, D., & Brunel, N. (2014). Analytical approximations of the firing rate of an adaptive exponential integrate-and-fire neuron in the presence of synaptic noise. Frontiers in Computational Neuroscience, 8, 116.PubMedCentralCrossRefPubMedGoogle Scholar
  30. Higgs, M.H., & Spain, W.J. (2011). Kv1 channels control spike threshold dynamics and spike timing in cortical pyramidal neurones. Journal of Physiology, 589(21), 5125–5142.PubMedCentralCrossRefPubMedGoogle Scholar
  31. Horton, J.C., & Adams, D.L. (2005). The cortical column: a structure without a function. Philosophical Transactions Of The Royal Society Of London. Series B, Biological Sciences, 360(1456), 837–62.PubMedCentralCrossRefPubMedGoogle Scholar
  32. Izhikevich, E.M. (2001). Resonate-and-fire neurons. Neural networks: The Official Journal of the International Neural Network Society, 14(6–7), 883–94.CrossRefGoogle Scholar
  33. Kuznetsov, Y.A. (1998). Elements of applied bifurcation theory(p. 591). Springer.Google Scholar
  34. La Camera, G., Rauch, A., Lüscher, H.-R., Senn, W., & Fusi, S. (2004). Minimal models of adapted neuronal response to in vivo-like input currents. Neural Computation, 16, 2101–2124.CrossRefPubMedGoogle Scholar
  35. Ladenbauer, J., Augustin, M., & Obermayer, K. (2014). How adaptation currents change threshold, gain and variability of neuronal spiking. Journal of Neurophysiology, 111, 939–953.CrossRefPubMedGoogle Scholar
  36. Lánský, P., & Lánská, V. (1987). Diffusion approximation of the neuronal model with synaptic reversal potentials. Biological Cybernetics, 56, 19–26.CrossRefPubMedGoogle Scholar
  37. Lu, Y., Monsivais, P., Tempel, B.L., & Rubel, E.W. (2004). Activity-Dependent regulation of the potassium channel subunits Kv1.1 and Kv3.1. Journal of Comparative Neurology, 470, 93–106.CrossRefPubMedGoogle Scholar
  38. Mainen, Z.F., & Sejnowski, T.J. (1995). Reliability of spike timing in neocortical neurons. Science, 268, 1503–1506.CrossRefPubMedGoogle Scholar
  39. Marder, E. (2012). Neuromodulation of neuronal circuits: back to the future. Neuron, 76(1), 1–11.PubMedCentralCrossRefPubMedGoogle Scholar
  40. Mattia, M., & Sanchez-Vives, M.V. (2012). Exploring the spectrum of dynamical regimes and timescales in spontaneous cortical activity. Cognitive Neurodynamics, 6, 239–250.PubMedCentralCrossRefPubMedGoogle Scholar
  41. Nataraj, K., Le Roux, N., Nahmani, M., Lefort, S., & Turrigiano, G. (2010). Visual deprivation suppresses L5 pyramidal neuron excitability by preventing the induction of intrinsic plasticity. Neuron, 68(4), 750–62.PubMedCentralCrossRefPubMedGoogle Scholar
  42. Nowak, L.G., Azouz, R., Sanchez-Vives, M.V., Gray, C.M., & McCormick, D.A. (2003). Electrophysiological classes of cat primary visual cortical neurons in vivo as revealed by quantitative analyses. Journal of Neurophysiology, 89(3), 1541–66.CrossRefPubMedGoogle Scholar
  43. Ostojic, S. (2014). Two types of asynchronous activity in networks of excitatory and inhibitory spiking neurons. Nature Neuroscience, (October 2013).Google Scholar
  44. Ostojic, S., & Brunel, N. (2011). From spiking neuron models to linear-nonlinear models. PLoS Computational Biology, 7(1), e1001056.PubMedCentralCrossRefPubMedGoogle Scholar
  45. Papoulis, A. (1965). Probability, random variables, and stochastic processes. McGraw-Hill Kogakusha.Google Scholar
  46. Pospischil, M., Piwkowska, Z., Bal, T., & Alain, D. (2011). Comparison of different neuron models to conductance-based post-stimulus time histograms obtained in cortical pyramidal cells using dynamic-clamp in vitro. Biological Cybernetics, 105(2), 167–80.CrossRefPubMedGoogle Scholar
  47. Pozzorini, C., Naud, R., Mensi, S., & Gerstner, W. (2013). Temporal whitening by power-law adaptation in neocortical neurons. Nature Neuroscience, 16(7), 942–8.CrossRefPubMedGoogle Scholar
  48. Prescott, S.A., & Sejnowski, T.J. (2008). Spike-rate coding and spike-time coding are affected oppositely by different adaptation mechanisms. Journal of Neuroscience, 28, 13649–13661.PubMedCentralCrossRefPubMedGoogle Scholar
  49. Renart, A., de la Rocha, J., Bartho, P., Hollender, L., Parga, N., Reyes, A., & Harris, K.D. (2010). The asynchronous state in cortical circuits. Science, 327(5965), 587–90.PubMedCentralCrossRefPubMedGoogle Scholar
  50. Ricciardi, L.M., & Sacerdote, L. (1979). The Ornstein-Uhlenbeck process as a model for neuronal activity. I. Mean and variance of the firing time. Biological Cybernetics, 35(1), 1–9.CrossRefPubMedGoogle Scholar
  51. Richardson, M.J.E., & Gerstner, W. (2005). Synaptic shot noise and conductance fluctuations affect the membrane voltage with equal significance. Neural Computation, 17(4), 923–947.CrossRefPubMedGoogle Scholar
  52. Rossant, C., Goodman, D.F.M., Fontaine, B., Platkiewicz, J., Magnusson, A.K., & Brette, R. (2011). Fitting neuron models to spike trains. Frontiers in Neuroscience, 5, 9.PubMedCentralCrossRefPubMedGoogle Scholar
  53. Sacerdote, L., & Giraudo, M.T. (2013). Stochastic integrate and fire models: a review on mathematical methods and their applications. Stochastic Biomathematical Models,Lecture Notes in Mathematics (Vol. 2058).Google Scholar
  54. Schreiber, S., Fellous, J.-M., Tiesinga, P., & Sejnowski, T.J. (2004). Influence of ionic conductances on spike timing reliability of cortical neurons for suprathreshold rhythmic inputs. Journal of Neurophysiology, 91, 194–205.PubMedCentralCrossRefPubMedGoogle Scholar
  55. Stimberg, M., Wimmer, K., Martin, R., Schwabe, L., Mariño, J., Schummers, J., Lyon, D.C., Sur, M., & Obermayer, K. (2009). The operating regime of local computations in primary visual cortex. Cerebral Cortex, 19(9), 2166–2180.PubMedCentralCrossRefPubMedGoogle Scholar
  56. Tabak, J., Senn, W., O’Donovan, M.J., & Rinzel, J. (2000). Modeling of spontaneous activity in developing spinal cord using activity-dependent depression in an excitatory network. Journal of Neuroscience, 20, 3041–3056.PubMedGoogle Scholar
  57. Touboul, J., & Brette, R. (2008). Dynamics and bifurcations of the adaptive exponential integrate-and-fire model. Biological Cybernetics, 99(4–5), 319–334.CrossRefPubMedGoogle Scholar
  58. Treves, A. (1993). Mean-field analysis of neuronal spike dynamics. Network, 4, 259–284.CrossRefGoogle Scholar
  59. Tsodyks, M., Uziel, A., & Markram, H. (2000). Synchrony generation in recurrent networks with frequency-dependent synapses. Journal of Neuroscience, 20, RC50.PubMedGoogle Scholar
  60. Turrigiano, G.G., & Nelson, S.B. (2004). Homeostatic plasticity in the developing nervous system. Nature Reviews. Neuroscience, 5(2), 97–107.CrossRefPubMedGoogle Scholar
  61. van Vreeswijk, C., & Hansel, D. (2001). Patterns of synchrony in neural networks with spike adaptation. Neural Computation, 13(5), 959–992.CrossRefPubMedGoogle Scholar
  62. Vogels, T.P., Rajan, K., & Abbott, L.F. (2005). Neural network dynamics. Annual Review of Neuroscience, 28, 357–376.CrossRefPubMedGoogle Scholar
  63. Vreeswijk, C.v., Sompolinsky, H., & van Vreeswijk, C. (1996). Chaos in neuronal networks with balanced excitatory and inhibitory activity. Science, 274(5293), 1724–1726.CrossRefPubMedGoogle Scholar
  64. Wang, X.J. (1998). Calcium coding and adaptive temporal computation in cortical pyramidal neurons. Journal of Neurophysiology, 79, 1549–1566.PubMedGoogle Scholar
  65. Wilson, H.R. (2003). Computational evidence for a rivalry hierarchy in vision. Proceedings of the National Academy of Sciences of the United States of America, 100, 14499–14503.PubMedCentralCrossRefPubMedGoogle Scholar

Copyright information

© The Author(s) 2015

Authors and Affiliations

  1. 1.Institut des Systèmes Intelligents et de Robotique (ISIR)CNRS UMR 7222, UPMC University ParisParisFrance
  2. 2.Institut d’Etudes de la Cognition, ENSParisFrance
  3. 3.Sorbonne Université, UPMC University Paris06 UMRS968, Insititut de la VisionParisFrance
  4. 4.INSERM, U968ParisFrance
  5. 5.CNRS, UMR7210ParisFrance
  6. 6.Department of Basic ScienceThe University of TokyoMeguro-kuJapan

Personalised recommendations


Tests have been carried out by Highways England to raise the upper limits from 50mph to 55mph or 60mph to help improve traffic flow and congestion.

Now Highways England has announced that drivers could be allowed to travel at increased speed through motorway roadworks depending on when they are travelling while further tests are planned to see if varying speed limits could safely be operated within a set of roadworks.

Motorway average speed cameras have also been tweaked where trials have taken place to allow for the new limits.

But Joshua Harris, director of campaigns for Brake, said: “Speed is the critical factor in the severity of collisions and so any move to increase limits is clearly a step in the wrong direction.

"Strict enforcement of safe speed limits is needed to assure the safety of workers and drivers, not an increase in speed. Average speed cameras are a proven and effective deterrent to speeding and we want to see their further rollout across the motorway network to limit dangerous driving behaviour.

“A recent Brake survey has found that a quarter of drivers think it likely that they will be involved in a fatal or serious crash on a motorway or dual carriageway at some point in the future.

"This shocking statistic makes clear that safety should be the overriding priority for investment on our motorway network, and we hope Highways England listens to the people and reconsider their plans."

There are also average speed cameras on the busy motorway for those travelling southbound from Leicester.

They start from around junction 19 and last for several miles.

Queues on the M1
Queues on the M1

It's part of a £15 billion road programme with around 150 smaller projects, such as work to install new slip roads, under way.

Traffic at Junction 21 on the M1
File picture of traffic queuing on the M1


Are you a member? Register / Log in
Does baby brain exist?

If you’ve got a little one on the way, you may find yourself blaming baby brain every time you lose your keys. But why do some women feel more forgetful during pregnancy?

Along with a whole list of things that happen to our bodies when we’re pregnant, so called ‘baby brain’ is one of the more bizarre. Your once perfectly functioning brain seems to have turned to mush overnight and now you can’t remember your pin code, phone number or even your partner’s name. But is this ‘mumnesia’ due to actual changes in your brain or simply a side effect of being more tired, emotional and preoccupied during pregnancy? Here are some theories…

Theory 1: Changes in the brain

One study done at Royal Holloway, University of London, found that pregnant women show more activity in the emotional, right side of the brain as they prepare to bond with their babies. This could explain why some women become more forgetful as they are, in turn, underusing the part of the brain that deals with memory.

Dr Victoria Bourne, from the Department of Psychology at Royal Holloway said, “The results suggest that during pregnancy, there are changes in how the brain processes facial emotions so that mothers are neurologically prepared to bond with their babies at birth.” So, it’s great news that your brain is giving you a helping hand to bond with your baby – not so great that you’ve lost your phone. Again.

Theory 2: Tiredness

During pregnancy a lot of your energy is being used to nurture your unborn baby. Couple this with raging hormones and lack of sleep and it’s hardly surprising you’re feeling exhausted. So, one explanation for baby brain could be that we’re just too wiped out to be on top of things. With this in mind, make sure you take every opportunity to rest and give your mind and body time to cope with all the changes. Don’t put yourself under pressure to do too much and listen to your body when it’s telling you to slow down. The HSE has more advice about about staying healthy and well during pregnancy.

Theory 3: Concentrating on other things

Pregnancy is a big deal, especially if it’s your first child. Not only do you have to cope with enormous changes happening to your body, you also have to prepare emotionally for how your life is about to change. So it’s little wonder that you might forget to put the washing on, or that you need to buy a pint of milk on your way home from work. Brenda McLackland, a consultant clinical psychologist, explains, “During most pregnancies women prioritise information about both their pregnancy and baby. Less relevant information may be missed and therefore will not be remembered. Tiredness is likely to exacerbate this further.” So don’t beat yourself up if you suffer from the occasional memory lapse. Maybe it’s just your brain doing what it’s designed to do – concentrating on your baby.

Theory 4: Expecting it to happen

Another theory is that women are so aware of the baby brain phenomenon that they blame pregnancy for normal memory lapses because it’s what’s expected of them. A study into cognitive changes in pregnancy had some interesting results. The study asked women about changes they’d noticed during pregnancy, and only three out of 198 women mentioned cognitive changes like forgetfulness. But when the same women were given a checklist of possible changes including options like memory, more women said they had experienced negative changes – for example, 32% reported that their memory was worse.

The study concluded that women would only report negative changes to their memory if they were asked specifically about it, not if they were thinking about how pregnancy affected them in general. So it could be that pregnant women are just attributing their usual forgetfulness to their pregnancy because they are expecting their memory to get worse.

It would seem that no one is quite sure what causes baby brain, or if it actually exists in the first place – there are plenty of theories but nothing conclusive. But that doesn’t make it any the less real for those who are experiencing it. So, whichever theory you want to believe, next time you find yourself at the shops still wearing your slippers or closing your front door with your keys still inside, try not to beat yourself up about it – you’ve got a much more important task ahead of you!

Like this? Read the eight embarrassing pregnancy side effects your pharmacist won’t blush at.


Compiled by Dr. Beth Gott of the School of Biological Sciences, Monash University.

All over Australia, Aborigines had their own local yearly calendars. Just as the climate in Kakadu is very different from that of Melbourne, so the Wurundjeri had their own way of marking the changing seasons. The division of the year into four seasons comes from Northern Europe, and does not fit Melbourne. We still think of winter as an unfavourable season for plants, when northern European trees drop their leaves and become dormant, but for our native plants, especially the small tuberous herbs, winter is a season of growth. At this time the bush is green, and the temperatures are rarely low enough to stop growth. The unfavourable season is high summer, when water is scarce, and much of the ground flora becomes brown and dies off. Water-plants such as Cumbungi are usually green during the summer, they die off during the winter.

Alan Reid originally suggested a calendar for the Middle Yarra region which has six seasons. Autumn, Winter, Pre-spring, True Spring, Early Summer and Late Summer. Glen Jameson (Victorian Naturalist 1996, Vol. 1 13 pp.26,67,123,269,313. 1997, Vol. 1 14, p.4S.) has modified this to: Late Summer, Early Winter, Deep Winter, Early Spring, True Spring, and High Summer, and has given lists of many plants and animals to be found in each of these seasons. Alan Reid agrees with Jameson, but emphasises that the seasons vary from place to place - for example the coastal seasons of the Bunwurrung would have been different.

Jones, Mackay and Pisani, from the University of Adelaide (Jones, D., Mackay, S. & Pisani, A. 1997 Patterns in the Valley of the Christmas Bush: a seasonal calendar for the upper Yarra Valley. Victorian Naturalist 114(5):246-249.) have proposed a seven season calendar.

This has Kangaroo Apple Season (December), Dry Season (Jan-Feb), Eel Season (March), Wombat Season (April-August), Orchid Season (September), Tadpole Season (October), Grass-flowering Season (November). The Museum of Victoria has chosen to follow this pattern in its displays in Melbourne and at Coranderrk. The Museum have named the seasons after the Woiwurrung language names for eels, wombats, and orchid, tadpoles and grass, but it is not known if these names were used by the Woiwurrung.

All the above are attempts to find the natural seasonal divisions of the area based on modern observations. Although William Thomas gave some seasonal happenings, that also did not represent the Aboriginal view. In the seasonal descriptions which follow we have combined elements of Jameson and Thomas, and have tried to include the activities of the people.

Lifestyle depended very much on the rhythm of the seasons. People moved over their extensive territories in regular cycles, maximising the availability of food and shelter. When food was plentiful there were opportunities to carry out ritual responsibilities and opportunities for ceremonial gatherings.

They were always mindful of the seasons in selecting the localities in which to spend their time, taking into account not only the natural features of the ground, but the facilities for obtaining food.

Much was known about the stars and the seasons. People could read the sky. For example, they would know that when a particular constellation appeared it was time to hunt a particular animal, or dig up a particular plant. See Thomas v.21, p98 "they are great observers of the starry heavens" for details. The flowering of particular plants was often used to mark seasonal events. Even today, fisherman use the flowering of the Coast Tea-tree in early November to mark the entry of Snapper into Port Phillip Bay.

We will start our survey in March, because this marked the end of the dry summer.

Late Summer, February - Mid March.
After the dry hot summer, the Autumn rains arrived and the days became cooler. People started burning those parts of the land where the scrub or tussock grass had become too dense during the summer, but they were advised by their elders where and when to burn, taking the weather into account so that the fires did not spread too far. It was important to clear the undergrowth and provide fertilising ash so that the small tuberous food plants could grow well after the rains came in March. Burning also made it easier to catch animals. Plants which had suffered from lack of water during the summer were now able to renew growth. DJAAK, Wattle gum wattle gum, was plentiful, and in the middle of this month the WARRAK Banksia or Honeysuckle, Long-leaf Box and Silver-leaf Stringybark came into blossom, providing sweet nectar, and attracting birds. March has been called the Eel Season, Eel because the female Short-finned eels were moving down the streams to the sea; the male eels had been leaving in smaller numbers during the spring and summer. These were an important food, and among the vegetables there were the starchy roots of the water plants, which began to die down after their summer growth. Some late summer fruits such as mistletoe berries were also available. Birds started to flock before heading north for the winter, to be replaced by other birds which will soon start to arrive from Tasmania.

Early Winter, April & May.
All sorts of fungi appeared with the rains, while the ground was still warm. BUNJIL, the Eagle, was building his nest, and the Brush-tail and Ringtail Possums were mating. Bolin Bolin billabong started to fill. Many different moths emerged, and were food for birds during the day and for Sugar and Feathertail gliders at night. Eastern Grey Kangaroos and Wallabies fed on the new growth.

Deep Winter- June, mid July
This cold time of the year slowed down but did not stop plant growth. Animals such as Echidnas were breeding, birds nesting. The flats near the rivers and creeks were often flooded; and the low lands generally were wet and cold, and unsuitable for camping, so people moved to the best sheltered spots on the uplands, where they were able to catch koalas, possums, and wombats, and to find grubs in the trees The leaves of the water plants had become dry and brown, but the small tuberous herbs were green and growing; the roots of both were good food. Fragrant nectar came from BURGILBURGIL, Honey-pots, Acrotriche serrulata, a small shrub which hid its flowers close to the ground. BULAIT- Cherry Ballart formed fruit. People constructed good bark WILLAMS (shelters) and kept fires burning for warmth. They wrapped themselves in rugs made from possum skins.

Early Spring - Mid July, August
As early as mid-July, MUYAN, Silver Wattle, started to flower, the first of the wattles to do so. It earned the name of 'Barak's Wattle' because when he died at Coranderrk on August 15th 1903, MUYAN was in full golden bloom. Yellow Box also flowered, providing much nectar. Early Nancy was the first of the small food plants to flower, and by late August the MURNONG was budding.Pacific Black Duck People moved slowly towards the lower lands as the spring temperatures rose, there they were able to snare ducks, to catch other kinds of wild-fowl, and, as the season advanced, to get eggs from the nests of all kinds of birds.

True Spring- September, October.
This was a time of plenty. Lilies, Orchids and MURNONG floweredmurnong, and still provided root vegetables. Greens were consumed in large quantities. Flowers were everywhere -Wattles, Hop Goodenia, Burgan, Kangaroo Apple, as well as orchids and small lilies which had been building their tubers over the winter. Snakes and Lizards became active, young Kangaroos came out of the pouch. Migrant birds - the Sacred Kingfisher for example, returned from the north. Tadpoles appeared in ponds, and the river, fed by melting snows from the mountains, flowed into the flood-plains and replenished the billabongs. Water-plants put on green leaves. Nowadays this flooding is prevented.

High Summer November, December, January
WARRA WARRAP/GARRONG, late Black Wattle, with pale yellow blossoms, flowered in November. As the summer advanced, the land began to dry, and people congregated around the reliable water-sources, the creeks, rivers and billabongs. Fish was an important food - Galaxias moved up the river from the sea. Where rocky falls blocked the river, as in the Prince's Bridge area and at Dight's Fails, fish would accumulate in large numbers, and could be easily taken. Eels started to come downriver. Fish traps were set. Water sources were important for the wildlife, so large animals such as Kangaroos and Emus would come to drink and could be caught. Lizards and snakes were active. Grasses flowered - Kangaroo Grass, Wallaby Grass, Spear Grass, Tussock Grass and the Common Reed. Fruits ripened -, MORR - Currant-bush, GARRAWANG - Apple-berry, White Elderberry, Kangaroo Apples and sweet LAAP - Manna, could be collected beneath the WURUN - Manna Gums The small tuberous plants died back, but the women still knew where they could dig for their roots, which at this time were at their best. . When people went up into the mountain gullies to get firedrills, they ate the pith from the centre of the treeferns. In the warm weather, big shelters were not needed unless it rained.

As food was plentiful, large gatherings of the tribes and clans took place. With permission from Bunurong clans, people went to the sea-coast to swim and gather shellfish and the fruits of Pigface and Coast Beard-heath. Flounder and Flathead could be speared or netted in the shallows, and shellfish were gathered.

The Dandenong Ranges were the hunting grounds for both the Bunurong (Western Port tribe) whose land lay to the south, and the Wurundjeri (Yarra Yarra) tribe whose land lay to the north and west. When the first pastoralists came: Blacks from the Western Port and Yarra Yarra tribes were frequently seen during the summer months, hunting in the forest for wallaby, possum and koalas.

Our knowledge of pre-European life in Victoria is fragmentary and sometimes confusing. Below is listed some of the names of Aboriginal seasons as recorded by three writers working separately during the nineteenth century.
WinterPer-ring-nger-wein, Perrin
(perein = no more sun)
wudawiin = cold weather)
Summer 1Nerrem-nre-wein
Nger-wein = sun)
Bullarto n'yoweenth
(plenty sun)
Summer 2 -Wygabil-ny-ewin
(Old Man sun)
AutumnBerripManemit (?good) birrin

Graphic by Damian Curtain & Stanley Barker.