|
|
||||||||
1Institute of Physiology, University of Bern, 3012 Bern, Switzerland; and 2Department of Biophysical and Electronic Engineering, University of Genoa, 16145 Genoa, Italy
Submitted 20 January 2004; accepted in final form 17 March 2004
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Even more important, the characterization of the single-neuron response properties, performed under appropriate and realistic conditions, is a key element to predict and understand how a population of neurons collectively interacts and processes information in the intact brain. This view is supported by several theoretical studies (Amit and Tsodyks 1992
; Brunel 2000
; Mattia and Del Giudice 2002
; Salinas 2003
), where the single-neuron response properties were used to make predictions about the collective phenomena, such as the global spontaneous irregular activity (Amit and Brunel, 1997b
), the emergence of network-driven oscillations (Brunel and Wang 2003
; Fuhrmann et al. 2002
), and of selective delay-activity states (Amit and Brunel 1997b
; Wang 2001
; Yakovlev et al. 1998
).
Under these perspectives, novel electrophysiological paradigms were recently proposed by several experimenters, focusing on the integrative (Destexhe and Paré 1999
, 2000
; Poliakov et al. 1997
; Rauch et al. 2003
), computational (Chance et al. 2002
; Mainen and Sejnowski 1995
; Protopapas and Bower 2001
), and adaptive properties of single neurons (Fuhrmann et al. 2002
; Paninski et al. 2003
) and synapses (Froemke and Dan 2002
; Sjöström et al. 2001
). In particular, Rauch et al. (2003)
proposed to mimic the conditions that neurons experience in the recurrent networks of the intact cortex. They reproduced a realistic background synaptic activity as a computer-synthesized noisy current that was injected into the soma, as shown in Fig. 1, A and B (see also Destexhe and Paré 2000
). Under such conditions, layer V pyramidal neurons in acute neocortical slices respond as integrate-and-fire (IF) point neurons (Rauch et al. 2003
). This suggests that an extremely simplified model neuron may be used to accurately describe the discharge properties of a single cell under simulated in vivo conditions and to predict and interpret collective emergent phenomena. Unfortunately, because it is not straightforward to observe, induce, or re-create coordinated network activity in neocortical slices (but see Giugliano and Lüscher 2003
; Giugliano et al. 2003
; and Shu et al. 2003
), there was no attempt at validating any network-level prediction in vitro, resulting from the proposed minimal characterization of the single cell. In addition, a direct confirmation of the same predictions in awake andbehaving animals is at present impossible or still requires extremely challenging technical issues to be overcome.
|
In this paper, we show that mature cultured neurons retain the features that characterize neurons in acute brain slices, investigated under the same experimental paradigm. Under these conditions, they respond as IF point neurons, not only with respect to the output mean firing rate but also to the subthreshold membrane voltage statistics, confirming and extending the work of Rauch et al. (2003)
. As a step toward the quantitative description of the mechanisms underlying population phenomena in cultured networks (Jimbo et al. 1999
; Kamioka et al. 1996
; Maeda et al. 1995
; Marom and Shahaf, 2002
), we computer-simulated the electrical activity of a network of model neurons, incorporating the single-neuron details identified in the experiments. We further interpret the results of our simulations by the extended mean-field theory (Renart et al. 2003
), and we compare them with in vitro experimental recordings, using the MEAs.
| METHODS |
|---|
|
|
|---|
Cultures of neocortical neurons were obtained from the somatosensory/motor cortex of newborn Wistar rats (P1-2), following standard procedures. Rats were anesthetized (0.4 ml Vetanarcol) and killed by decapitation. Transverse brain slices, 225 µm thick, were exposed to a 0.3% trypsin solution for 3 min at 37°C for enzymatic digestion. Cells were mechanically dissociated and plated on substrate arrays of planar microelectrodes (MEAs) or glass coverslips, at a density of 150,000 or 75,000/150 µl, respectively. MEAs were microfabricated as described previously (Tscherter et al. 2001
) and coated for 1 h with diluted (1:50) Matrigel (Falcon/Biocoat, Becton Dickinson AG, Basel, Switzerland). The glass coverslips were coated with polylysine (1 mg/ml overnight at 37°C or Matrigel (1:50). Cells were restricted to a small area (
50 mm2) using cloning glass cylinders attached to the MEAs or coverslips by silicone sealant. Cultures were incubated at 36.5°C in a 5% CO2 air atmosphere and maintained in 150 µl nutrient medium. The medium contained MEM Eagle (Sigma) supplemented with fetal bovine serum 10%, glucose 0.2%, B27, and Glutamax (Gibco BRL, Life Technologies AG, Basel, Switzerland) for the cultures on glass coverslips, and serum-free Neurobasal medium supplemented with B27 and Glutamax (Gibco BRL) for the cultures on the MEAs. Half of the medium was changed weekly.
Recordings were made in a chamber mounted on the stage of an upright microscope (Nikon, Tokyo, Japan). Patch-clamp experiments were carried out on cells in which a few processes could be visually identified (e.g., Fig. 1E), from cultures of 47 wk of in vitro age. Within this period, no systematic change of single-neuron and network properties was detected, suggesting a complete maturation of the neurons and indicating that the cultures' development had reached a steady state (Kamioka et al. 1996
; Marom and Shahaf 2002
).
Pharmacology
Several minutes before starting each recording session, the culture medium was replaced by an extracellular solution containing (in mM): NaCl 145, KCl 4, MgCl2 1, CaCl2 2, HEPES 5, Na-pyruvate 2, glucose 5, at pH 7.4. Recordings were performed either in the absence of a continuous flow of solution, with solution changes every 1015 min, or in the presence of continuous superfusion at 1 ml/min. No differences were detected between these protocols. All recordings were made at room temperature (2324°C).
As extensively described in the literature, mature cultures of dissociated neocortical neurons exhibit a spontaneous electrical activity, which results in simultaneous bursting network activity and mainly depends on the synaptic development and coupling between the neurons (Jimbo et al. 1999
; Kamioka et al. 1996
; Maeda et al. 1995
; Marom and Shahaf 2002
). In the patch-clamp experiments, given that our goal was to characterize the single-neuron response properties to a class of computer-synthesized stimulus waveforms, single-neuron recordings were performed while blocking glutamatergic synaptic transmission with D-APV (D-2-amino-5-phosphonovalerate 50 µM), a competitive antagonist of the NMDA (N-methyl-D-aspartate) receptors, and CNQX (6-cyano-7-nitroquinoxaline-2-3-dione, 10 µM) (both Sigma, Buchs, Switzerland), a competitive antagonist of non-NMDA receptors. These substances were bath applied and completely suppressed incoming synaptic activity and spontaneous spiking in the neurons. In the experiments involving multisite extracellular recordings of the neuronal electrical activity by the MEAs, an antagonist of GABAA receptors (bicuculline methochloride, 10 µM) (Tocris, Anawa Trading SA, Wangen, Switzerland) was bath applied to block synaptic inhibition, with the aim of focusing on the disinhibited pattern generation as well as establishing a more direct comparison to the computer-simulated networks, which included only excitatory model neurons.
Whole cell patch-clamp recording from single neurons
The patch-clamp technique was used, in the whole cell configuration (Hamill et al. 1981
), by using an Axoclamp-2B amplifier (Axon Instruments, Union City, CA). Signals were recorded in current clamp, filtered at 1 kHz, sampled at 5 kHz, and digitized by a 12-bit A/D converter (Digidata 1200) and pClamp 8 software (Axon Instruments). Electrodes were pulled from filamented borosilicate glass capillaries (GC150F, Harvard Apparatus GmbH, March-Hugstetten, Germany) on a horizontal puller (DMZ, Zeitz Instrumente GmbH, Munich, Germany) and their resistance was 5.5 ± 0.6 M
. Electrodes were filled with a solution containing (in mM: K-gluconate 100, KCl 20, HEPES 10, Mg-ATP 4, Na2-GTP 0.3, Na2 phosphocreatine 10; pH 7.3, 290 mOsm. Other (standard) pipette solutions were reported not to alter significantly the response properties of the cells, under the very same experimental protocol (Rauch et al. 2003
). High-resistance seals (24 G
) were formed and the whole cell configuration was achieved by the application of a negative pressure pulse. The bath application of D-APV and CNQX always followed the establishment of the whole cell patch configuration.
Conventional procedures, consisting of repetitive hyperpolarizing step currents, were used to obtain a direct estimate of the passive (cable) properties of cultured neurons (Iansek and Redman 1973
). Such estimates are related to the time constant
m and capacitance Cm of an equivalent lumped RC compartment (Abbott and Dayan 2001
).
Emulating a stationary realistic input from the network
Independent realizations of the stochastic OrnsteinUhlenbeck process were computer-synthesized and injected under current clamp (see Fig. 1C). (Interactive stimuli-synthesis and analysis software tools were developed and are available on request.) Such nondeterministic current stimuli mimic a realistic barrage of excitatory/inhibitory postsynaptic currents (EPSCs/IPSCs) for a cell embedded in a large in vivo network, spontaneously and randomly active at low rates (Destexhe et al. 2003
), as well as in a cultured network where spontaneous neurotransmitter release at individual synapses, as well as other sources of inhomogeneity and randomness determine an irregular background synaptic noise in vitro (Maeda et al. 1995
). Although a conductance injection would have appeared more appropriate (e.g., by means of the dynamic-clamp technique; see Destexhe and Pare 2000), it has been proven that noisy conductance-driven and noisy current-driven stimulations are equivalent with respect to the evoked neuronal steady-state mean firing rates, apart from 2 stimulus-independent scaling factors, for the input mean and variance (La Camera et al. 2003
; Rauch et al. 2003
). More important, the comparison of the experimental data to the predictions of mathematical model neurons has been carried out under the same current-driven conditions. Under these hypotheses and because the impact of a single EPSC/IPSC on the postsynaptic membrane voltage is weak in evoking suprathreshold responses in vitro (Nakanishi and Kukita 1998
; Nakanishi et al. 1999
), the overall current experienced by a generic postsynaptic neuron can be approximated as a diffusion stochastic process (Destexhe and Paré 2000
; Fourcaud and Brunel 2002
). Therefore, under extended mean-field hypotheses on the interactions between neurons in a population (Amit and Brunel 1997b
), such noisy stimuli may statistically account for a wide range of feedforward/recurrent network architectures and regimes. For instance, by indicating with Ne/i the number of excitatory/inhibitory afferents with stationary mean activation rates fe/i, under the hypothesis that synaptic inputs are approximately independent, the distribution of the resulting postsynaptic somatic current amplitude becomes Gaussian (i.e., by the central-limit theorem), with steady-state mean m and variance s2 given by the expressions reported below (Amit and Brunel 1997b
)
![]() |
e/i are the effective peak-amplitude and decay time constant at the soma, for individual excitatory and inhibitory postsynaptic currents, respectively. Throughout this work, we set
e =
i =
I
{1; 5} ms, thereby mimicking (AMPA- and GABAA-mediated) fast synaptic currents (Destexhe et al. 1994
We note that the features of the incoming synaptic current arising from the particular spike-timing precision and reliability of individual presynaptic neurons were not explicitly investigated (Mainen and Sejnowski 1995
; Jolivet et al. 2004). Consistently, the neuronal response to such noisy current input was routinely analyzed by characterizing the output mean spike rate at the steady state. However, provided that the hypothesis on the statistical independence of the presynaptic activity holds, the impact of neuronal precision and reliability is expected to weakly contribute to the network activity we discuss in the present work. Spontaneous synaptic release and other sources of network randomness will in fact introduce uncorrelated stochastic components to the overall synaptic inputs to any neuron. Therefore, under such an hypothesis, the overall resulting current incoming to a generic postsynaptic neuron of the network is still Gauss-distributed and completely characterized by mean m and variance s2.
Noisy current-clamp protocol
The stimulation protocol consisted of the repeated somatic injection of independent current realizations I(t), each lasting 20 s and interleaved by 3060 s of recovery time. For any pair (m, s2), the following iterative expression was used to synthesize a realization I(t) of the process (Cox and Miller 1965
)
![]() | (1) |
t is a unitary Gauss-distributed random variable (Press et al. 1992Analysis of the single-neuron response properties
Collected data consisted of the voltage responses to each noisy stimulus, lasting 20 s and associated with the pair (m, s2) (see Fig. 1D). To account for current-clamp offsets, the actual injected current was also monitored and the actual value of (m, s2) directly estimated from it. Raw voltage traces were processed in Matlab (The MathWorks, Natick, MA) by a peak-detection algorithm to extract individual spike times and shape. When no substantial change in the shape of the individual action potentials, or in the instantaneous firing rate at the beginning and at the end of the elicited spike train occurred (representative of possible nonstationarities), the trial was accepted. The steady-state mean spiking frequency f was estimated, discarding an initial transient (i.e., 25 s) and averaging across the 2 available repetitions. After a successful completion of the stimulation protocol, the experimental curve {fh = fh(mh, sh), h = 1, 2,... , M} (for at least 3 distinct values of s and typically M
20) was plotted and compared to the theoretical responses of 2 model neurons, driven by the same current statistics (Fig. 2). The recordings were further compared to the data available for pyramidal layer V neurons in acute slices, collected under the same experimental conditions (Rauch et al. 2003
).
|
The mathematical modeling and the procedures for the model fit to the collected response data are standard and closely followed Rauch et al. (2003)
. Briefly, 2 single-compartment standard mathematical descriptions of neuronal excitability were considered: the Lapicque's or leaky integrate-and-fire neuron LIF (Abbott and Dayan 2001
; Tuckwell 1988
) and the constant leakage integrate-and-fire neuron with a floor (CLIFF) (Fusi and Mattia 1999
; Mongillo and Amit 2001
; Rauch et al. 2003
) (see Fig. 2). These models have been studied in depth (Fourcaud and Brunel 2002
; Mattia and Del Giudice 2002
) and widely used in simulations of large-scale networks (Mattia and Del Giudice 2000
; Reutimann et al. 2003
) and hardware implementations (Chicca et al. 2003
). As opposed to the biophysically realistic conductance-based models (Abbott and Dayan 2001
), these descriptions are characterized by a single state variable V (i.e., the membrane potential) and by a reduced set of effective constant parameters: the absolute refractory period
arp and reset voltage H, the membrane capacitance C, the voltage threshold for spike emission
, and the subthreshold voltage decay rate C(E V)/
for the LIF (
being the membrane effective time constant and E the resting potential) or
for the CLIFF model. Finally, both models incorporated a simplified spike-frequencydependent adaptation, modeling the contribution of intracellular calcium- and/or sodium-activated outward currents to the net membrane current, modulated by a constant factor
and implemented as described in Liu and Wang (2001)
and van Vreeswijk and Hansel (2001)
. The stationary effect of such an adaptation is to reduce the gain of the frequencycurrent curve by a factor that does not depend on the adaptation time constant
a (see Eqs. 2 and 3). We refer as aLIF and aCLIFF, to the LIF and CLIFF models incorporating the spike-frequencydependent adaptation, respectively (see Fig. 2B).
Model parameters fit to the experimental data
The full analysis of the response of the spiking neuron models to the input current, specified by Eq. 1 (i.e., a colored noise), can be performed only under approximate treatments (Fourcaud and Brunel 2002
) but it was not considered in the present work. Instead, it is considerably easier to analyze the models' response to an idealized "equivalent" white current (i.e.,
-correlated), characterized by the infinitesimal mean µ and variance
2 determining an asymptotically equivalent effect on the subthreshold membrane voltage V (Rauch et al. 2003
). The last is a satisfactory approximation, when the evoked mean interspike interval is much larger than
I, as verified by computer simulations. Under these hypotheses, the steady-state current-to-rate response function
(m, s2) can be analytically determined (Fusi and Mattia 1999
; Tuckwell 1988
) (Fig. 2B) and the corresponding parameter space may be searched for the best fit to the stationary data points {fh(mh, sh), h = 1, 2,... , M}, collected for each cell. The implicit formulas corresponding to the mean firing rate
aLIF/aCLIFF under noisy current stimulation are
![]() | (2) |
![]() | (3) |
![]() |
![]() |
= RC the membrane time constant (Fig. 2A), whereas for the aCLIFF model
![]() |
was kept fixed at 20 mV, with respect to the resting membrane potential (i.e., E = 0 mV and
= E + 20 mV).
For each cell, simulated-annealing optimization techniques (Press et al. 1992
) were used to fit the above reported theoretical response functions to the data, by minimizing the following mean-square error with respect to the model parameters
![]() |
![]() |
![]() |
1, the resulting accuracy interval
h± for fh corresponds approximately to a confidence of at least 68%.
Finally,
2 and its minimum
min2 are random variables, known to be approximately distributed according to a
2 distribution (Press et al. 1992
). This makes it possible to refer to the probability Prob (
2 >
min2) as a standard indication of the goodness of the fit. Model fits were accepted when such a probability was >0.1. The minimal value of K of successful
2-test was chosen as a quantitative analogue measure of the quality of the fit, comparing the performances of the aLIF, LIF, aCLIFF, and CLIFF models over the entire data set. Finally, the parameters search was repeated twice for each experiment: under free-search conditions and introducing additive quadratic cost-penalty constraints to the simulated annealing energy (Press et al. 1992
). This aimed at discouraging the exploration of the parameter space in a limit-valued region (i.e., H
and
arp
+
). We refer to these conditions as free search, or no-penalty best fit, and as penalty best fit.
MEA recording and analysis of the network activity
In some experiments, MEAs were used as a substrate for the cultured networks, so that neuronal somata as well as axons were allowed to develop close to the individual MEA microelectrodes. Such a proximity makes possible the extracellular detection of the emission of action potentials by one or more neighboring cells (Streit et al. 2001
). MEAs contained 68 platinum planar electrodes, spaced at 200-µm intervals and laying out in the form of a rectangle. Recording channels showing activity were selected by eye and their recordings digitized at 6 kHz per channel and stored on a hard disk. The detection of extracellularly recorded action potentials (i.e., fast voltage transients) and further analysis were performed off-line in IGOR. (WaveMetrics, Lake Oswego, OR), as described previously (Streit et al. 2001
; Tscherter et al. 2001
). No attempt was made to sort spikes detected by the same MEA electrode. The electrical noise of individual channels was very stable. The selectivity of event detection was routinely checked by using recordings obtained in the presence of tetrodotoxin (TTX, 1.5 µM, Sigma), as a "zero" reference.
The spontaneous network activity, detected by the MEA, consisted of asynchronous activity and population bursts (PBs). The asynchronous activity was defined as one or a few spikes detected by a single or several electrodes, whereas the PBs, consisting of episodes of activity occurring simultaneously across several recording channels, spread over the network. The processed network activity was visualized in the form of event raster plots and of the instantaneous population mean firing rate (see Fig. 5A). The last was computed by counting the total number of detected events from all recording channels, within a sliding time window of 10 ms. With the aim of summarizing the features of the detected network activity, the mean and the coefficient of variation (CV) of the interburst intervals (IBI) distribution and of the burst durations (PBd) distribution were estimated over at least 10 min of continuous recording (15 PBs minimum).
|
We studied and computer-simulated the collective electrical activity of small networks of Ne = 1001,000 interacting excitatory aLIF identical neurons (van Vreeswijk and Hansel 2001
), using the single-neuron effective parameters identified in the experiments (Table 1). We chose to neglect network inhomogeneities to keep the interpretation of the results as simple as possible. Anyway, our interpretative framework can be extended to cover a similar situation, following the approach described in Amit and Brunel (1997a)
.
|
![]() |
![]() |
], the previous equations are complemented by the conditions
![]() |
![]() |
(t) the unitary step function (i.e.,
(t) = 0, t < 0 and
(t) = 1, t > 0), the total synaptic current Ii(t) into the neuron ith is given by
![]() |
Focusing on cultured neocortical networks, we chose an unstructured connection topology (Marom and Shahaf 2002
) with a probability Cee of synaptic connection between any 2 neurons of 0.30.4 (i.e. Prob {Cij = 1|i
j} = Cee) (Nakanishi and Kukita 1998
; Nakanishi et al. 2001
). In agreement with the hypotheses underlying the stimuli injected in the experiments, the synaptic interactions were described by currents rather than conductance changes. The adaptation currents were also described as current changes, as well. Describing these variables by conductance changes does not qualitatively affect the results reported. Synaptic interactions between 2 connected neurons were triggered by the presynaptic emission of action potentials, after an effective delay
of 1.5 ms, which included the axonal propagation delay and synaptic release latency (Nakanishi and Kukita 1998
). The resulting individual postsynaptic currents were characterized by an instantaneous rise to Je and by an exponential decay with a time constant
e = 5 ms, and no activity-dependent short- or long-term change. The effect of spontaneous synaptic release and of other sources of randomness (Maeda et al. 1995
) was incorporated into an activity-independent additional random synaptic drive. The effective value of
unavoidably included the superimposed stationary contribution of both fast and slower adaptation mechanisms, estimated in the patch-clamp experiments over 20 s of stimulation time. With the aim of focusing on the time scales characterizing the network bursting (i.e.,
0.11 s), sometimes we decreased the value of
that affected the PBd and not the IBIs statistics. For the sake of simplicity, the intrinsic cumulative inactivation experimentally measured was not included in the simulations (but see Giugliano et al. 2002
).
Under the same extended mean-field hypotheses that underlie the noisy currents that were injected into real neurons, the collective activity of the simulated network can be fully predicted in terms of its stationary mean firing rate f, by the knowledge of
aLIF(m, s2) and the details of the synaptic connectivity (Amit and Brunel 1997b
). In the present case, the statistics of the total recurrent synaptic current experienced by a generic neuron of the network is approximately Gauss-distributed and it can be fully described by its mean and variance
![]() | (4) |
Statistics
The Pearson's r linear correlation as well as the Kendall's Tau nonparametric (rank-order) test (Press et al. 1992
) were used to assess statistical correlations. The last provides a correlation measure together with an estimated significance level pK, which corresponds to the probability of obtaining the same correlation from statistically independent samples. Averages are expressed as means ± SE.
| RESULTS |
|---|
|
|
|---|
As opposed to a DC stimulation (see Fig. 1A), under noisy current injection, the neuronal membrane voltage evolves in time as in a random walk, leading to irregular spike emission (see Fig. 1, A and B). Its subthreshold amplitude distribution becomes bell-shaped, with mean and SD increasing with the steady-state mean m and variance s2 of the injected current, respectively (see Figs. 1D and 4).
|
) was indeed required by the models to fit the data points at the steady state (La Camera et al. 2002
No difference was found in the response properties of the neurons to noisy currents characterized by an autocorrelation time length
1 = 1 ms (n = 17) and 5 ms (n = 18). Thus, in both cases the responses could be captured by the white-noise input approximation, mentioned under METHODS.
However, collected data resembled the aLIF and not the aCLIFF responses (Fig. 3), as opposed to pyramidal cells in acute slices, where both models were reported to describe the data set equally well (Rauch et al. 2003
). The aLIF model was substantially better in fitting the responses of the cells over the entire set of experiments (n = 35) (see Fig. 3 and Table 1).
|
Best-fit effective parameters and passive membrane properties
As mentioned in METHODS, the passive membrane properties (
m and Cm) of the patched neurons were routinely measured. As opposed to the previous report by Rauch et al. (2003)
, the passive properties and the corresponding best-fit effective parameters of the IF models showed a significant cross-correlation, both for the best-fit free search and the penalty search. We report the results of the Pearson's r-test and, between parentheses, those of the Kendall's Tau test. In the best-fit free search, between Cm and C the correlation coefficient was 0.66 (0.54) for the aLIF and 0.73 for the aCLIFF model, whereas between
m and
it was 0.77 (0.57) for the aLIF. In the best-fit search with penalties, between Cm and C the coefficient was 0.56 (0.44) for the aLIF and 0.68 for the aCLIFF model, whereas between
m and
it was 0.81 (0.61) for the aLIF. Such results were validated by a very high significance of the nonparametric test (i.e., pK < 105; see METHODS).
Interspike interval variability and subthreshold voltage distribution
As stated in the previous sections, by using the set of best-fit effective parameters, the steady-state response of the aLIF model accurately matches the corresponding experimental mean firing rates, in the plane (m, s2). In layer V pyramidal neurons the best-fit parameters could sometimes account for the CV of the interspike intervals distribution as well (Rauch et al. 2003
), although the fit criterion involved the mean firing rates only. Carrying out a similar analysis on the aLIF model in cultured cells, it turned out that model responses poorly matched the CV, experimentally estimated at the steady state (not shown).
However, when the steady-state subthreshold distribution of the membrane voltage, recorded in the experiments, was estimated and compared to the aLIF model prediction, a good agreement was observed over the available set of pairs (m, s2) (see Fig. 4). To quantitatively compare the aLIF model behavior to the experimental voltage traces, individual action potentials were clipped and a unique offset and a scaling factor were required for the model internal state-variable V to optimally match the voltage distribution (i.e., V
V E'). For each cell, these additional 2 parameters do not affect the current-to-rate response function and they were the same over the pairs (m, s2) (Fig. 4).
Slow/cumulative inactivation and the stationary spike frequency
As indicated in METHODS, special care was taken in assessing the stationarity of the neuronal responses, with the aim of direct comparison with the model predictions, available at the steady state. Nevertheless the temporal dynamics of the output firing rate was characterized by fast (frequency-dependent) adaptation components, occurring over a time scale of several hundreds of milliseconds (see Fig. 1D), and by a slower component, occurring over a time scale of several seconds. Although the steady-state effect of both adaptation processes was captured by the model (i.e., by
; see METHODS), a cumulative (inactivation-related) component contributed to set an upper limit to the maximal stationary response frequency, sustained by the neuron over the entire duration of the stimulation (Fleidervish et al. 1996
; Powers et al. 1999
; Rauch et al. 2003
; Sanchez-Vives et al. 2000
; Sawczuk et al. 1997
). For instance, the injection of a noisy stimulus current for 20 s with m > 200300 pA or more, depending on the cell input resistance, resulted in a slowly decaying instantaneous firing rate, eventually turning into a cumulative inactivation of the action-potential generation (not shown). This was reminiscent of the slow cumulative sodium-current inactivation characterized by Fleidervish et al. (1996)
, and made it impossible to quantify the steady-state responses above about 30 Hz, with the current stimulation protocol.
The presence of such a cumulative inactivation was explicitly tested (n = 8, 2143 DIVs) by extending the protocol described in (Fleidervish et al. 1996
; Schwindt et al. 1989
), consisting in a repeated-pulse stimulation, lasting 1 s, with a very short recovery time. Under noisy current injection, by using the same current realization for each repetition, the same phenomenon occurred, although the voltage fluctuations, induced by the nondeterministic stimulus waveform, delayed the onset of the inactivation at parity of m, and sometimes transiently reversed the inactivation for a few tens of milliseconds, compared to DC stimuli. As expected, larger values of
1 induced slow modulations on the membrane voltage trajectory (see Svirskis and Rinzel 2000
), making the episodes of partial transient recovery more frequent, at parity of m and s (not shown).
The spontaneous emergence of patterned network activity in vitro
Networks of neocortical dissociated neurons show spontaneous collective patterned activity (Fig. 5A). Such an activity starts as asynchronous and spatially uncorrelated firing, toward the end of the first week in culture, and evolves into nonperiodic, synchronized, population bursting after 34 wk (Kamioka et al. 1996
). Such an activity is reported to stay unchanged for more than 8 wk and thus represents the mature state of the network (Marom and Shahaf 2002
). During such a period, single cells emit rare and irregular spikes or bursts of action potential, superimposed on spontaneous voltage fluctuations around a membrane potential of about 60 mV (Nakanishi and Kukita 1998
). A similar spontaneous network-driven activity constitutes a typical feature of dissociated neuronal networks (Marom and Shahaf 2002
) and was also observed in adult neocortex (Sanchez-Vives and McCormick 2000
).
The recording of the network activity by means of MEAs led to a quantitative characterization of the distributions of the interburst interval (IBI) and of the population burst durations (PBd), in 7 cortical cultures. Under control extracellular medium, the IBIs were characterized by a mean ranging from 4.6 ± 0.4 to 30.3 ± 4.8 s, and by a CV ranging from 44 to 70%. The PBs were characterized by a mean ranging from 54 ± 8 to 146 ± 8.5 ms, and by a CV ranging from 26 to 80%. Under pharmacological disinhibition (see METHODS), the PBs became longer with a mean ranging from 600 ± 70 ms to 1.7 ± 0.053 s, and a CV ranging from 7 to 30%, whereas IBIs were characterized by a mean ranging from 12.7 ± 1.3 to 51.8 ± 11.1 s, and by a CV ranging from 31 to 96%. Figure 5D summarizes graphically the results from 4 experiments, under bicuculline.
Simulations of networks of aLIF model neurons
By using the single-neuron effective parameters, we computer-simulated the collective electrical activity, emerging from a small homogeneous population of excitatory aLIF neurons (see METHODS). We considered estimates for the synaptic connectivity available from the literature (i.e., Cee; see METHODS), and we set (m0, s0) in Eq. 4 to match the spontaneous low-rate asynchronous background activity observed in the MEAs experiments (
1 Hz).
Similarly to previous theoretical reports (Segev et al. 2001
; Tateno 2002
; van Vreeswijk and Hansel 2001
; Wiedemann and Lüthi 2003
), but with a stronger original motivation for the use of the aLIF model, and a quantitative goal to compare computer simulations to the available MEAs recordings, we found that the simulated network activity consisted of asynchronous activity and of population bursts.
Although the individual neurons of the simulated network were not intrinsic burster cells and no pacemaker mechanisms had been introduced in the simulations, the emerging activity evolved into a collective network-driven bursting, depending on the strength of synaptic coupling Je (see Fig. 8). In agreement with the experiments of Maeda et al. (1995)
, the spatial origin of PBs varied randomly with each burst, consistent with the conclusions about the lack of a unique pacemaker mechanism, driving the network.
|
a, which was set in the simulations in the range of 0.72 s (see Fig. 5, B and C, gray shading). Transiently and following each PB, the network activity was almost completely suppressed by the adaptation currents in each neuron and it later recovered. On the other hand, the synaptic coupling as well as the strength of the spontaneous synaptic release or other sources of randomness (i.e., m0 and s0) were correlated to the bursting frequency.
Three different global regimes were observed. The first one corresponds to a situation in which the excitatory synaptic interactions between neurons are very weak or the network connectivity is very low. Under such conditions, the network of model neurons was exclusively dominated by low-rate asynchronous activity (
1 Hz). No PB occurred either spontaneously or evoked by any brief depolarizing stimulus (as opposed to Fig. 7A). Actually, in such a regime the global dynamics was dominated by a single low-rate stable state. These conditions approximate early developmental stages of cultured networks, described in the literature to display a very similar electrical activity pattern.
|
. Actually, in such a regime the global dynamics of the network is transiently bistable and once a PB is started, the adaptation slowly redefines the location and the existence of the stable states for the network dynamics, until there is suddenly only a single stable state at 0 Hz. This somehow corresponds to a reset for the collective network activity. In details, the adaptation hyperpolarizing contribution, which starts to build up in every neuron recruited by the PB, decreases the output firing rate until the recurrent synaptic inputs to any neuron stop. In such a regime, the statistics of the simulated IBIs and PBd matched the experiments, performed under pharmacological disinhibition (see Fig. 5, BD).
The third and last regime occurs for stronger excitatory synaptic interactions in the network. Under such conditions, the network is characterized by a high-rate asynchronous regime (
5580 Hz). In this regime, the firing of the individual neurons is more regular than in the previous ones. Moreover, any attempt at transiently silencing the activity of the network, by hyperpolarizing a large fraction of the neurons of the population, would not prevent the network to later recover its global state (as opposed to the simulations reported in Fig. 7). Although, under physiological conditions, it is probably not possible for a cortical network to sustain such a firing regime for a long time, an interesting phenomenon was observed for an intermediate synaptic interactions strength: occasional population breaks. This appears as rare and unpredictable simultaneous short interruptions of the global activity and is determined by the finite-size fluctuations of the activity.
The single-neuron response and the mean-field theory interpretation
In the present section, we show that the quantitative knowledge of the single-neuron response function, identified in the previous experiments, is very relevant to predict and interpret the emergence of the network activity described above. As widely discussed in the literature, a homogeneous network of synaptically interacting excitatory neurons may be regarded as a single dynamical system. Its stationary states can then be predicted and interpreted, in the limit of an infinite number of neurons Ne
, by using the mean-field Eq. 4 and studying
aLIF [m(f), s(f)] as a function of f (Amit and Brunel 1997b
).
We first consider the collective activity in the absence of spike-frequencydependent adaptation (i.e.,
= 0). Under such hypotheses, because of the simultaneous dependency of
aLIF on m and s2, the collective firing rate of the network may be characterized by 2 stable dynamic-equilibrium states, in a small range of average synaptic coupling Je (Figs. 6A and 8, left panel) (Amit and Brunel 1997b
; Fusi and Mattia 1999
). Such global activity configurations correspond to the solutions f* of the following self-consistent network equation, further satisfying a stability condition
![]() |
aLIF(f) with the unitary-slope line (see Figs. 6A and 7, insets).
|
m0 ±
m0) and resulting into an intrinsically bistable network activity (Fig. 7A). This has been already described as a possible neuronal correlate of the selective delay-activity states in vivo (Yakovlev et al. 1998The mechanisms of a PB
Considering the full network model, where individual neurons keep adapting their output rate as a function of the activity, it is possible to carry out a simple approximate analysis. Provided that the mechanisms responsible for the excitability reduction (e.g., the adaptation) act on a time scale (
a) that is longer compared to the single-neuron dynamics (
), an analysis of the quasi-stationary equilibria may tell us a lot about the collective activity of the network. In other words, by assuming that adaptation is delayed and transiently uncoupled from the neuronal dynamics, we m