|
|
||||||||
1Departamentos de Biofísica and 2Neurociencias, Instituto de Fisiología Celular, Universidad Nacional Autónoma de México, Mexico City, Mexico
Submitted 12 October 2007; accepted in final form 8 January 2008
|
|
ABSTRACT |
|---|
|
|
|
INTRODUCTION |
|---|
|
Basal ganglia (BG) contain CPGs that activate innate behavioral routines, procedural memories, and learned motor programs (Barnes et al. 2005
; Graybiel 1995
; Grillner et al. 2005a
,b
; Takakusaki et al. 2004a
). A major component of the BG is the striatum, which receives a widespread input from the cerebral cortex and thalamus. Striatal circuits process cortico-thalamic inputs to produce specific outputs consisting of bursts of action potentials during the execution of motor tasks, probably following voltage transitions to depolarized "up-states" (Hikosaka et al. 2000
; Kasanetz et al. 2006
; Mahon et al. 2006
; Romo et al. 1992
; Schultz et al. 1993
; Wilson 1993
). As in other isolated nervous tissue preparations known to contain CPGs (e.g., Guertin and Hounsgaard 1998
), addition of NMDA to neostriatal circuits in vitro (Vergara et al. 2003
) and in vivo (Herrling et al. 1983
) induces recurrent bursting and pattern generation in single neurons. Moreover, intrastriatal application of NMDA in vivo generates turning behavior when administered unilaterally in freely moving animals (Ossowska and Wolfarth 1995
); demonstrating that motor behaviors arise from the striatal processing of enhanced excitatory drives. This evidence suggests that the striatum posses the connectivity and intrinsic mechanisms to orchestrate pattern generation (Grillner 2006
).
To test this hypothesis, we used calcium imaging of neuronal populations in a corticostriatal slice preparation to monitor, with single-cell resolution, dozens of cells simultaneously and then discern pattern generation produced at the microcircuit level by the activity of cell assemblies. Simultaneous electrophysiological recordings from striatal neurons demonstrated that calcium transients result from neurons bursting on top of suprathreshold up-states (Kerr and Plenz 2002
).
|
|
METHODS |
|---|
|
Transverse corticostriatal slices (300 µm thickness) were obtained from PD14-29 Wistar rats as previously described (Kawaguchi et al. 1989
; Vergara et al. 2003
). All procedures conformed to the guidelines of the UNAM's Animals Scientific Procedures Committee. Slices were obtained with ice-cold saline (4°C) containing in mM: 123 NaCl, 3.5 KCl, 1 MgCl2, 1 CaCl2, 26 NaHCO3, and 11 glucose (25°C; saturated with 95% O2-5% CO2; pH = 7.4; 298 mosM/l). Slices were then transferred to saline at room temperature (21–25°C) where they remained for
1 h before recording. The cationic concentration of this saline favors the appearance of up states in vitro (Vergara et al. 2003
). Although the present results were performed mainly in young animals (PD14-21), basically the same network behavior was observed in older animals (PD25-29) (J. Bargas, unpublished data and see RESULTS).
Calcium imaging
Slices were incubated at room temperature in the dark for 20–30 min in the presence of 10–20 µM fluo 4-AM (Tef Labs, Austin, TX) in 0.1% dimethylsulphoxide (35°C), equilibrated with 95% O2-5% CO2. Slices were perfused with control saline (see preceding text) in a perfusion chamber on the stage of an upright microscope equipped with a x10 water-immersion objective (Eclipse E600FN; Nikon, Melville, NY). Excitation at 488 nm was performed with a Lambda LS illuminator (Sutter instruments, Novato CA). Experiments were performed at room temperature.
Images were acquired with a cooled digital camera (SenSys 1401E, Roper Scientific, Tucson, AZ) at 250–500 ms/frame. Imaging software used was RS Image (Photometrics; Roper Scientific). The imaged field was 800 x 600 µm in size. Short movies (100–250 s, 50- to 100-ms exposure, 2–4 image/s) were taken at time intervals of 5–20 min during 1 h.
The number of fluo 4-loaded neurons in the field was determined at the end of the experiment with a 5-s puff of 50 mM KCl. This maneuver disclosed all fluo-4-labeled neurons (either active or silent during the experiment). Cells active during the experiment were analyzed, and the ratio of active/silent cells was obtained. Spontaneous or evoked calcium transients together with voltage responses were recorded electrophysiologically in some cells, both in control saline and during the application of 5–12 µM NMDA (Sigma-Aldrich-RBI, St. Louis, MO). In some experiments, cortical sensory motor areas were stimulated with a concentric bipolar electrode (12 µm; FHC, Bowdoinham, ME). Stimuli consisted of different trains of 5–10 stimuli at 20 Hz. Each stimulus was 100–200 µs and 50–120 µA. In some experiments, we used the minimal stimulus intensity necessary to evoke peaks of synchrony with amplitudes above chance (P < 0.05). This allowed us to study changes of electrical evoked synchrony under different pharmacological conditions.
Immunohistochemistry
Sections were processed to fix fluo-4 fluorescence—or cells active during experiments—with N-(3-dimethylaminopropyl)-N'-ethylcarbodiimide hydrochloride (EDAC), and to perform conventional immunocytochemistry in fluorescent cells to demonstrate either substance P (SP) or enkephalin (ENK) on fluo-4-labeled cells using commercially available antisera (Peninsula Labs, San Carlos, CA) conjugated to CY3 or CY5. Slices were not processed for both antisera, but one was chosen in each case. Thus in each trial, either SP or ENK positive and negative neurons could be observed. Briefly, sections were rinsed in PBS and incubated for 18–24 h at 4°C with primary rabbit antibody against ENK or SP (diluted 1:200). Sections were mounted in an anti-quenching media (Vectashield, Vector Laboratories) and examined under a confocal microscope (MRC-1024; Bio-Rad, Natford, UK) equipped with a krypton–argon mixed-gas laser. Immunostained cells were studied in either single confocal images or reconstructed sections made by projecting z-series of 10–40 consecutive confocal images 1 µm apart collected throughout the thickness of the section. The background noise was reduced averaging three to six images. Digitized images were transferred to a personal computer (Confocal Assistant, T. C. Brelje). More than 80% of fluo-4-loaded cells were medium spiny neurons.
Drugs
Stock solutions were prepared before each experiment and added to the perfusion solution in the final concentration indicated. NMDA, APV, nicardipine, CNQX, biocytin, and bicuculline methiodide or hydrochloride were obtained from Sigma (St. Louis, MO).
Electrophysiology
Calcium imaging and simultaneous electrophysiological recordings were obtained from areas of the dorsal striatum previously shown as receiving numerous cortical fibers (Vergara et al. 2003
).
An Axoclamp 2B amplifier (Axon Instruments, Foster City, CA) was used to perform whole cell current- and voltage-clamp recordings. Signals were filtered at 1–3 kHz and digitized at 3–9 kHz with an AT-MIO-16E4 board (National Instruments, Austin, TX) in a PC computer. Data acquisition used a software designed in the LabView environment (Lemus-Aguilar et al. 2006
). Patch pipettes (3–6 M
) were filled with (in mM) 115 KH2PO4, 2 MgCl2, 10 HEPES, 0.5 EGTA, 0.2 Na2ATP, and 0.2 Na3GTP. In some experiments, biocytin 0.5%, and fluo-4 salt (20–30 µM) were added to the recording pipettes.
Image analysis
Image processing was carried out with Image J (v.1.36, National Institutes of Health), Multicell 2.0 (kindly supplied by Robert Froemke), and custom-made programs written in IDL (Cossart et al. 2003
; Mao et al. 2001
; Schwartz et al. 1998
) or MATLAB (The Math-Works, Natick, MA).
All active neurons in a field were semi-automatically identified, and their mean fluorescence was measured as a function of time. Single-pixel noise was discarded using a 5-pixel ratio mean filter. Calcium-dependent fluorescence signals were computed as (Fi –Fo)/Fo, where Fi: fluorescence intensity at any frame and Fo: resting fluorescence, i.e., average fluorescence of the first four frames of the movie.
Calcium signals elicited by action potentials were detected based on a threshold value given by their first time derivative (2.5 times the SD of the noise value). Thus we obtained a C x F binary matrix; were C represents the number of active cells and F the number of frames for each movie. Spike onsets were signaled by ones in the matrix representing transitions to the up states. Recordings were inspected manually to remove artifacts and slow calcium transients which are likely to correspond to glial cells (Ikegaya et al. 2005
; Sasaki et al. 2007
).
Statistical methods
To determine if calcium transients recorded from different cells were correlated, the numbers of simultaneous activations per trial (onset of signals occurring within 3 frame windows) were detected. To determine the P value of simultaneous transients occurring by chance, the distribution under the null hypothesis of independent transients using Monte Carlo simulations with 1,000 replications were computed (Mao et al. 2001
).
The degree of correlation between active cells was calculated with the Jaccard correlation coefficient. Nevertheless the magnitude of the correlations is difficult to discern when many lines are superimposed. Thus we constructed cross-correlation maps of Jaccard correlation coefficients to show the magnitude of the correlations between all cells pairs.
In addition, we identified the sets of cells activated simultaneously over time. To identify peaks of synchronous activity (i.e., that included more cells than those expected by chance), Monte Carlo simulations were also used to estimate the significance of their firing together. The threshold corresponded to a significance level of P < 0.05. We then examined peaks that were significant at P < 0.05. Peaks of synchronous and recurrent bursting activity that remained significant during the experiment were selected for further analysis.
To analyze the dynamics of cell assemblies over time (network dynamics), D x N matrices were constructed, where D represents the number of active neurons in a set of experiments, and N denotes the firing of cells during 250-ms to 1-s time bins (NMDA-induced up states last between 0.5 and 5 s) (Vergara et al. 2003
). Peaks of recurrent synchronous activity were vectorized so that bursting over time was associated with different neurons. Each vector element is formed by the sum of the number of calcium spikes displayed by a single neuron during the time bin, where peak derivatives denote the time onset of electrophysiologically recorded up states (see preceding text and RESULTS). Therefore the set of these vectors denote network activity as a function of time (Brown et al. 2005
; Sasaki et al. 2007
).
To measure the similarity index between the network vectors, we computed the normalized dot product of all possible vector pairs, which is equivalent to the cosine of the angle between the vectors (Sasaki et al. 2006
; Schreiber et al. 2003
). Then we plotted the similarity indexes as a pseudocolor matrix in which functional states sustained by significantly correlated or synchronized bursting cell assemblies appear as cluster-like structures (Sasaki et al. 2006
).
To reduce the dimensionality of the network vectors, locally linear embedding (LLE) was chosen. LLE is an unsupervised learning algorithm that discloses nonlinear structures from multi-dimensional data (Roweis and Saul 2000
). Multidimensional scaling (MDS) (Systat, Richmond, CA) gave results quantitatively similar from those obtained with LLE. Nevertheless LLE depicted more effectively the trajectories of the functional states in the network. The trajectories of our high-dimensional data were not well described by linear dimensionality reduction methods (see Brown et al. 2005
). To choose the optimal number of states depicted from the LLE reduction, we used hard and fuzzy clustering algorithms taking the Dunn's index as a validity function (Bezdek et al. 1997
; Sasaki et al. 2007
). To determine the number of neurons in each state, hierarchical cluster analysis was computed using Euclidean distances and the nearest neighbor single linkage method (Systat Software, San Jose, CA).
|
|
RESULTS |
|---|
|
To study striatal microcircuits in vitro, we used Ca2+ imaging (Fig. 1) to measure electrical activity in many cells simultaneously with single-cell resolution. Activity was measured indirectly, as changes in fluorescence. Fields from the dorsal striatum were imaged in slices loaded with the Ca2+ indicator fluo-4 AM in 174 experiments performed in 86 corticostriatal slices. Figure 1A shows all dye-loaded neurons (see METHODS). Contours of cells both active and inactive during the experiment are depicted in Fig. 1B (filled circles and empty contours, respectively). Only neurons active during an experiment were analyzed. Under control conditions (with no drugs added), only a few cells were active, and their firing was asynchronous (filled circles in Fig. 1B). Records of Ca2+ transients from three spontaneously active neurons are shown (Fig. 1C). These experiments confirm that the striatal circuitry is mostly quiescent (n = 40 slices) under control conditions.
|
F/F) corresponding to up states, had time derivatives [d(
F/F)/dt] that matched bursts duration (Cossart et al. 2003
|
In corticostriatal slices, electrical stimulation of the cortex evokes long-lasting depolarizations with overriding spikes in medium spiny neurons (Bargas et al. 1991
; Vergara et al. 2003
). Trains of cortical stimuli (Fig. 3A) also result in prolonged synaptic depolarizations with action potentials (Fig. 3B). Figure 3B illustrates an experiment where whole cell current-clamp recordings (Fig. 3B1) and calcium imaging (Fig. 3B2) were performed simultaneously from a medium spiny neuron (stimulus frequency: 0.1 Hz; stimuli are signaled with arrows at the bottom). A calcium transient accompanies each orthodromic response to the cortical stimulus (Fig. 3B2, arrows at the bottom). The time derivatives of the calcium responses are shown in the third row (Fig. 3B3); note similar duration of derivative positive peaks and voltage responses. Accordingly, each peak from the differentiated calcium signal generates a dot used to build raster plots as that illustrated in Fig. 3D. Each row in the raster plot represents an active neuron (filled circles in Fig. 3C). Filled circles in Fig. 3C indicate fluo-4-loaded neurons within the field of observation that responded to cortical stimulation in one representative experiment. Empty circles indicate loaded neurons that did not respond to cortical stimulation (see METHODS). Responsive cells are a minority of the total number of neurons which are scattered throughout the observational field among many unresponsive neurons (Fig. 3C).
|
10–20% of fluo-4-loaded cells. Only a few neurons were active before and after the stimulus train, and these cells present a lack of correlation (Fig. 3, D and E). At each trial, >90% of the responsive neurons followed the cortical stimulus (Fig. 3D, bottom). Some neurons spontaneously active before the stimulus did not follow the electrical stimulation (blue lines). Targeted electrophysiological recordings of cells responsive to cortical stimulation revealed that most active neurons were medium spiny neurons (Fig. 3B). Network dynamics set into action by NMDA
Ca2+ imaging allowed us to observe a population of striatal cells activated by cortical stimulation in control conditions (Fig. 3) and in the bath presence of NMDA (5–12 µM; Fig. 4). Simultaneous voltage recordings of imaged neurons (Fig. 4A) show that voltage transitions (Vergara et al. 2003
) were accompanied by Ca2+ transients during NMDA activation (Kerr and Plenz 2002
, 2004
). Cell activity in the circuit was followed by recording their Ca2+ transients and plotting them as raster plots (Fig. 4B). Similar NMDA-induced electrical activity has been observed in neurons that are part of a CPG (e.g., Grillner et al. 1981
; Guertin and Hounsgaard 1998
; Hsiao et al. 1998
; Takakusaki et al. 2004b
). Interestingly, the raster plot (Fig. 4B, top) shows that several neurons are involved in the activity. Moreover, the time histogram (Fig. 4B, bottom; asterisks, red lines) clearly shows periods of increased synchrony and correlated firing that occur spontaneously in different sets of neurons (Cossart et al. 2003
; Ikegaya et al. 2004
). Statistically significant (threshold P < 0.05) spontaneous peaks of synchrony occurred in sets of neurons from n = 38/45 slices in these conditions. On average, there were 2.5 ± 0.2 peaks of synchronous bursting per 201 ± 10-s time epoch. Peaks of recurrent synchronous activity showed no apparent periodicity, and occurred with a mean interval of 46 ± 10 s. Notably, neurons firing synchronously during NMDA treatment could be hundreds of microns apart intermingled with silent cells (Fig. 4C, filled circles). Spatial correlation maps (Fig. 4D) show the pairs of neurons that exhibit statistically significant (P < 0.05) correlated activity (line thickness is proportional to the degree of correlation). On average, 77 ± 4% of active cells had correlated activity (n = 38 slices) at any given moment. Cells active during synchrony peaks represent most of the correlations between cells. Correlation plots of this activity (Fig. 4E; P < 0.05) show that the degree of synchronization among active neurons is heterogeneous, thus the network activity is not due to chance correlation between any pair of neurons but to the spatiotemporal dynamics of several cells.
|
|
Visualizing functional states during network dynamics
We next investigated network dynamics in the striatum over long periods of time (n = 15 slices). Brief movies (100- to 250-s epochs) were taken at different intervals for up to an hour (e.g., Fig. 6C, 1–5, separated by lines). To analyze network dynamics, we binned the raster plots of all the neurons involved in the peaks of synchrony (Fig. 6). Time segments represent vectors coding coactive burst activity of individual neurons (see METHODS). Network dynamics consisted of a set of N vectors (bins) in D dimensions (active neurons). Vectorization allowed us to compare different states of network activity rigorously over time (Brown et al. 2005
; Sasaki et al. 2007
). The normalized inner product (see METHODS) of all possible vector pairs provides a measure of the similarity among states (Sasaki et al. 2006
; Schreiber et al. 2003
). To evaluate this similarity along network activity, we plotted all vector pairs as a matrix (Fig. 6A). Notably, abrupt transitions in the similarity index showed the presence of cluster-like structures (Sasaki et al. 2006
).
|
20% the size of those induced by cortical stimulation (cf. Fig. 6C, 1–5, bottom). Nonetheless the general behavior of the active network after a cortical stimulus (Fig. 6C, 4 and 5) remained essentially comparable to the activity before the stimulus (Fig. 6C2 and see following text). In fact, time histograms showed similar numbers of spontaneous peaks of synchronous activity (Fig. 6C, bottom; P < 0.05) before and after the cortical stimulus. Moreover, cortical stimulation (Fig. 6C3) elicited a trajectory looping back to the NMDA-induced dynamics (Fig. 6B), showing recurrent activation of the same assembly after a perturbation (Stopfer et al. 2003Pattern generation denoted by functional states sustained by cell assemblies
The preceding analyzed states included some with imposed synchronization, i.e.: with electrical cortical stimulation. If network dynamics induced by NMDA are mediated by cell assemblies, then one may expect the appearance of distinguishable functional states evoked by the excitatory tonic drive brought about by NMDA only (without electrical stimulation). To address this issue, we investigated the relationships between the peaks of synchronous activity induced only by NMDA (similar results were obtained in n = 14/15 slices that presented NMDA-induced peaks of synchrony).
To observe alternation between functional states, we investigated the NMDA-induced network dynamics in the striatum over long periods of time (Fig. 7, n = 15 slices). Brief image series (100–250 s in duration) were obtained at different intervals for
1 h (Figs. 7C, 1-,3 top, representative image series obtained at different times during the experiment are shown separated by vertical lines).
|
Figure 7D shows the spatial distribution of neurons generating the patterns and producing the different functional states. Network states share some elements but most of the neurons are dissimilar (Fig. 7E). Finally, hierarchical cluster analysis confirmed the recruitment of different cell assemblies during specific states (Fig. 7F). Nevertheless there is a small core of cells shared by all the network states (
8% of neurons involved in the synchrony peaks).
We concluded that functional states, denoted by spontaneous peaks of synchrony (Barnes et al. 2005
), are generated by the coordinated participation of cell assemblies, as it is the case of "unit CPGs" (Grillner 2006
). The transitions from one state to the other display long-term network dynamics. The return of the same cell assemblies after net traveling through different trajectories shows that the elements of these modules are robustly associated, allowing alternating participation. Also, a core of neurons is shared by all the states.
Mechanisms underlying network dynamics for pattern generation
For a preliminary investigation into the synaptic and intrinsic requirements for multistable network dynamics, we challenged network activity with both synaptic and intrinsic ion channel antagonists. Fast synaptic inhibition within the striatal circuit (Czubayko and Plenz 2002
; Koos et al. 2004
; Tepper et al. 2004
; Tunstall et al. 2002
) was tested with the GABAA receptor antagonist, bicuculline, which was applied once the network became active.
Figure 8 shows that on exposure to 10 µM bicuculline, the number of state transitions is drastically reduced in the NMDA-treated slices (n = 6 slices). The network becomes locked in a preferred state, which is recurrently revisited. Similarity indexes representing network dynamics suggest the absence of diverse cell assemblies (Fig. 8A). Only occasional transitions to another state occurred (Fig. 8B). Paradoxically, the raster plot of network activity showed an increased frequency of recurrent synchrony peaks: from 2.7 ± 0.2 peaks in the NMDA condition to 8.3 ± 2 peaks in the presence of bicuculline (P < 0.01). There was also an increase in the number of coactive cells supporting the preferred state (Fig. 8C). Most synchrony peaks, however, correspond to the same dominant state. Neurons involved in network sates are shown in Figs. 8, D and E. Hierarchical cluster analysis did not reveal cluster-like structures in spite of synchronization (Fig. 8F), meaning that all active neurons were embedded into the same state. Cortical stimulation in the presence of bicuculline recruited even more striatal neurons than those obtained in control conditions or in the presence of NMDA alone (data not shown). Thus fast GABAergic transmission is a necessary requirement for the network to exhibit multistable dynamics and different functional states.
|
Active networks were challenged with the AMPA receptor antagonist CNQX (10 µM) (Fig. 9A). As shown in the raster plot of Fig. 9A1, CNQX drastically reduced the number of active cells (cf., left and right) and disrupted the appearance of synchrony peaks in all NMDA-treated slices (Fig. 9A, bottom; n = 6 slices). Interestingly, a few neurons that were active during synchrony peaks maintained spontaneous activity in the absence of fast AMPA transmission (Fig. 9A, 1 and 2, gray circles) (Vergara et al. 2003
). Also in the presence of CNQX, cortical stimulation was unable to recruit striatal neurons (data not shown) and all correlated activity was suppressed (Fig. 9A).
|
Shared core of neurons
In most slices presenting multistable network dynamics (n = 20/25 slices), we found a core of neurons, which were shared by all states. So we aimed at recording electrophysiologically some of these neurons after their identification as a part of a core.
These neurons were 8 ± 2% of all active cells in the network states (n = 20 slices). Some of these core neurons exhibited periodic calcium transients during NMDA treatment or CNQX. Bicuculline disrupted the periodic firing of these cells. Cell-attached (Fig. 10A1) or whole cell current-clamp recordings (Fig. 10, C1 and D1) were performed on some of these cells, and we found that 40% of the core assembly (n = 4/10 cells) showed firing properties characteristic of GABAergic interneurons (Plenz and Aertsen 1996
; Tepper et al. 2004
), such as high-frequency bursts (Fig. 10A) and periodic pacemaking activity (Fig. 10D, 1–3). However, by fixing fluorescence of cells active during a given experiment (see METHODS), Fig. 11 shows that most active cells in the different assemblies were immunoreactive to either substance P or enkephalines.
|
|
|
|
DISCUSSION |
|---|
|
Ca2+ imaging of neuronal populations (Mao et al. 2001
; Sasaki et al. 2006
; Schwartz et al. 1998
; Stopfer et al. 2003
) and simultaneous whole cell recordings show that calcium transients correspond to burst firing in single cells, the synchronization and recurrence of which generates the activity patterns. Pattern generation during network states depend on synaptic and intrinsic properties as well as the activity of a core of active neurons (Berke et al. 2004
). In summary, this work shows how pattern generation and correlated activity is generated in the striatal microcircuitry in vitro identifying the neuronal elements involved. The connections and modulation of these elements deserve further study because they may reveal basic general properties of striatal microcircuitry and may further suggest in vivo experiments.
Voltage transitions in striatal neurons
Striatal neurons, both in vivo and in vitro, exhibit spontaneous voltage transitions that sustain burst firing. Transitions occur between a quiescent down state and a recurrent bursting up state (Herrling et al. 1983
; Kerr and Plenz 2002
, 2004
; Mahon et al. 2006
; Vergara et al. 2003
; Wilson 1993
). It is largely unknown how these transitions reflect network activity and how they propagate through the network. It is probable that different classes of bursts exist, for example, during different functional states (e.g., sleep vs. movement) (N. Vautrelle, personal communication). Here we took advantage of the tonic drive provided by NMDA to induce the transitions because in the striatum (Ossowska 1995
), as in other circuits (Gordon and Whelan 2006
; Grillner et al. 1981
; Guertin and Hounsgaard 1998
; Hsiao et al. 1998
; Kiehn 2006
; Takakusaki et al. 2004b
), NMDA sets into action patterned activity that generates movement. Accordingly, it was demonstrated that NMDA induces correlated activity that evolves in time. In contrast to control conditions in which the striatum has only a few active cells with no correlated activity, both cortical stimulation and the tonic excitatory drive provided by NMDA exposure (Vergara et al. 2003
) recruit neurons spread over a wide area; which, in the case of NMDA, spontaneously express particular spatiotemporal dynamics, indicating that the striatal microcircuit in vitro preserves a set of "unit CPGs" (Grillner 2006
), that is, even if most connections present in more intact preparations are severed, remaining sets of neurons, probably belonging to larger CPGs or modules in vivo, preserve some of the properties of these larger modules. In support of this inference, independent evidence in the cortex suggests that cell assemblies activity is preserved along different spatial scales (Plenz and Thiagarajan 2007
).
Mechanisms of network dynamics
The conditions that generate circuit dynamics reside in both the synaptic and intrinsic properties of striatal neurons. NMDA induces synchronous peaks of neuronal activity emerging from the co-activation of robust and recurrent cell assemblies that alternate in their patterns displaying a sequence of trajectories. These results suggest the presence of robust mechanisms that maintain and stabilize the cells participating in the assemblies. In particular, AMPA transmission appears to be necessary for turning on network dynamics. In fact, the number of active neurons drastically falls when AMPA transmission is blocked.
Because all glutamatergic synapses in the striatum originate from cortical or thalamic afferents, cortico-thalamic afferents are likely to be the source of tonic excitatory drive underlying synchronous activity. NMDA presumably generates an increase in the tonic excitatory drive conveyed by cortico-thalamic afferents (Grillner 2006
; Grillner et al. 1981
), and although some neurons maintain spontaneous activity in the presence of CNQX—or after dissecting away the cerebral cortex (Vergara et al. 2003
)—the analysis shows that this behavior is asynchronous so that the emergence of striatal cell assemblies appears to require fast AMPA glutamatergic transmission, presumably originating from cortical and thalamic neurons (Kasanetz et al. 2006
; Yuste et al. 2005
). The most parsimonious interpretation of these results is that cortico-thalamic afferents are essential for the appearance of multistate network dynamics. That is, network dynamics are commanded by the cortex and/or the thalamus (Kasanetz et al. 2006
; Magill et al. 2001
). However, the actual relationship between cortical or thalamic activation and cell assembly activity in the striatum needs further investigation. It is known that activation of NMDA receptors is enough to induce bistability in striatal neurons, amplifying the excitatory drive conveyed by cortical or thalamic inputs (Grillner 2006
; Grillner et al. 1981
) so that activity (e.g., from the cortex or thalamus) can be relied on. And although cell assembly activity has been recorded in the cortex using NMDA plus dopamine, it has not been recorded after NMDA only (Tseng and O'Donnell 2005
). Finally, striatal assemblies receive the convergence of widespread cortical areas, and slices cannot contain all these areas (Parthasarathy and Graybiel 1997
). However, cortical cell assemblies driving striatal cell assemblies are not excluded and may be possible in several conditions (e.g., with added dopamine receptor agonists). Therefore striatal assemblies may represent the coordinated activity between cortical areas, that once formed by development and learning, only need a cortical trigger or a subset of the original stimuli for activation.
In addition, our data show that the striatum does not simply follow cortical commands. The role of the striatal circuit proper in the management of cortical drives is revealed by the result of inhibiting GABAA receptors. Blockade of GABAA inhibition not only increases the peaks of synchronous activity and recruits more cells during synchrony peaks but also shifts the network toward a preferred state that recruits most active neurons. Without inhibition, the striatal circuit gets tied up in a preferred state, from which alternation and selection among different cell assemblies becomes unlikely. Therefore we conclude that inhibitory circuits in the striatum are responsible for transforming cortical commands into a sequential activity of striatal cell assemblies (Ossowska 1995
; Wickens and Oorschot 2000
). In addition, the blockade of L-type calcium channels disrupts the peaks of synchrony. This result confirms that intrinsic conductances, such as voltage-gated calcium channels, by generating plateau potentials may work as synchronization enablers of neuronal networks (Kiehn 2006
; Yuste et al. 2005
). Thus NMDA-induced network dynamics in the striatum requires the participation of both synaptic and intrinsic mechanisms.
Finally, electrophysiological recordings of neurons exhibiting Ca2+ transients showed that most of the active neurons in cell assemblies were medium spiny projection neurons. This was confirmed by immunocitochemistry (Fig. 11). Nevertheless, a subpopulation of neurons shared by all functional states (see Parthasarathy and Graybiel 1997
) exhibited periodic Ca2+ transients and the firing properties of pacemaking GABAergic interneurons (e.g., Berke et al. 2004
; Tepper et al. 2004
). The fact that striatal GABAergic interneurons comprise 2–5% of striatal cells and yet make up 40% of the neurons active in multiple assemblies suggests that striatal interneurons participate in the orchestration of multistate dynamics (Berke et al. 2004
). Activity of pacemaker neurons is a frequent finding in CPGs (Grillner 2006
; Yuste et al. 2005
).
Functional implications and perspectives
Recurrent and alternating bursting is characteristic of cell assemblies included in CPGs in vivo and in vitro (Barnes et al. 2005
; Grillner 2006
; Ikegaya et al. 2004
). Synchronized activity of these modules exhibit attractor dynamics (Cossart et al. 2003
), providing a unified description for circuits encoding the storage and retrieval of long-term and working memory. Attractors are seen as memory traces retrieved through excitatory tonic driving or partial cues useful for executing motor programs. The persistence of cell assemblies along time is due to recurrent connectivity (Barnes et al. 2005
; Tsodyks 2005
). The alternation of their activity is under neuromodulatory control (Yuste et al. 2005
). Apparently, the properties of these microcircuits are at the interface between small networks and global brain functions (Plenz and Thiagarajan 2007
). Their disturbance provokes abnormal processes of synchrony, associated with different disorders such as schizophrenia and Parkinson Disease (Schnitzler and Gross 2005
; Uhlhaas and Singer 2006
). We show that an isolated circuit set into action, in vitro, can exhibit synchronized states in specific cell assemblies emerging and returning during relatively long periods of time.
|
|
GRANTS |
|---|
|
|
|
ACKNOWLEDGMENTS |
|---|
|
|
|
FOOTNOTES |
|---|
Address for reprint requests and other correspondence: J. Bargas, Instituto de Fisiología Celular UNAM, PO Box 70-253, Mexico City, DF 04510 Mexico (E-mail: jbargas{at}ifc.unam.mx)
|
|
REFERENCES |
|---|
|
Barnes TD, Kubota Y, Hu D, Jin DZ, Graybiel AM. Activity of striatal neurons reflects dynamic encoding and recoding of procedural memories. Nature 437: 1158–1161, 2005.[CrossRef][Medline]
Berke JD, Okatan M, Skurski J, Eichenbaum HB. Oscillatory entrainment of striatal neurons in freely moving rats. Neuron 43: 883–896, 2004.[CrossRef][Web of Science][Medline]
Bezdek JC, Li WQ, Attikiouzel Y, Windham M. A geometric approach to cluster validity for normal mixtures. Soft Comput 1: 166–179, 1997.
Brown SL, Joseph J, Stopfer M. Encoding a temporally structured stimulus with a temporally structured neural representation. Nat Neurosci 8: 1568–1576, 2005.[CrossRef][Web of Science][Medline]
Cossart R, Aronov D, Yuste R. Attractor dynamics of network UP states in the neocortex. Nature 423: 283–288, 2003.[CrossRef][Medline]
Czubayko U, Plenz D. Fast synaptic transmission between striatal spiny projection neurons. Proc Natl Acad Sci USA 99: 15764–15769, 2002.
Froemke RC, Kumar VS, Czkwianianc P, Yuste R. Analysis of multineuronal activation patterns from calcium-imaging experiments in brain slices. Trends Cardiovasc Med 12: 247–52, 2002.[CrossRef][Web of Science][Medline]
Gordon IT, Whelan PJ. Deciphering the organization and modulation of spinal locomotor central pattern generators. J Exp Biol 209: 2007–2014, 2006.
Graybiel AM. Building action repertoires: memory and learning functions of the basal ganglia. Curr Opin Neurobiol 5: 733–741, 1995.[CrossRef][Web of Science][Medline]
Grillner S. Biological pattern generation: the cellular and computational logic of networks in motion. Neuron 52: 751–766, 2006.[CrossRef][Web of Science][Medline]
Grillner S, Hellgren J, Menard A, Saitoh K, Wikstrom MA. Mechanisms for selection of basic motor programs—roles for the striatum and pallidum. Trends Neurosci 28: 364–370, 2005a.[CrossRef][Web of Science][Medline]
Grillner S, Markram H, De Schutter E, Silberberg G, LeBeau FE. Microcircuits in action-from CPGs to neocortex. Trends Neurosci 28: 525–533, 2005b.[CrossRef][Web of Science][Medline]
Grillner S, McClellan A, Sigvardt K, Wallen P, Wilen M. Activation of NMDA-receptors elicits "fictive locomotion" in lamprey spinal cord in vitro. Acta Physiol Scand 113: 549–551, 1981.[Web of Science][Medline]
Guertin PA, Hounsgaard J. NMDA-induced intrinsic voltage oscillations depend on L-type calcium channels in spinal motoneurons of adult turtles. J Neurophysiol 80: 3380–3382, 1998.
Herrling PL, Morris R, Salt TE. Effects of excitatory amino acids and their antagonists on membrane and action potentials of cat caudate neurons. J Physiol 339: 207–222, 1983.
Hikosaka O, Takikawa Y, Kawagoe R. Role of the basal ganglia in the control of purposive saccadic eye movements. Physiol Revs 80: 953–978, 2000.
Hsiao C, del Negro CA, Trueblood PR, Chandler SH. Ionic basis for serotonin-induced bistable membrane properties in guinea pig trigeminal motoneurons. J Neurophysiol 79: 2847–2856, 1998.
Hutchinson WD, Dostrovsky JO, Walters JR, Courtemanche R, Boraud T, Goldberg J, Brown P. Neuronal oscillations in the basal ganglia and movement disorders: evidence from whole animal and human recordings. J Neurosci 24: 9240–9243, 2004.
Ikegaya Y, Le Bon-Jego M, Yuste R. Large-scale imaging of cortical network activity with calcium indicators. Neurosci Res 52: 132–138, 2005.[CrossRef][Web of Science][Medline]
Ikegaya Y, Aaron G, Cossart R, Aronov D, Lampl I, Ferster D, Yuste R. Synfire chains and cortical songs: temporal modules of cortical activity. Science 304: 559–564, 2004.
Izhikevich EM. Dynamical Systems In Neuroscience. Cambridge, MA: MIT Press, 2007.
Kasanetz F, Riquelme LA, O'Donnell P, Murer MG. Turning off cortical ensembles stops striatal up states and elicits phase perturbations in cortical and striatal slow oscillations in rat in vivo. J Physiol 577: 97–113, 2006.
Kawaguchi Y, Wilson CJ, Emson PC. Intracellular recording of identified neostriatal patch and matrix spiny cells in a slice preparation preserving cortical inputs. J Neurophysiol 62: 1052–1068, 1989.
Kerr JN, Plenz D. Dendritic calcium encodes striatal neuron output during up-states. J Neurosci 22: 1499–1512, 2002.
Kerr JN, Plenz D. Action potential timing determines dendritic calcium during striatal up-states. J Neurosci 24: 877–885, 2004.
Kiehn O. Locomotor circuits in the mammalian spinal cord. Annu Rev Neurosci 29: 279–306, 2006.[CrossRef][Web of Science][Medline]
Koos T, Tepper JM, Wilson CJ. Comparison of IPSCs evoked by spiny and fast-spiking neurons in the neostriatum. J Neurosci 24: 7916–7922, 2004.
Lemus-Aguilar I, Bargas J, Tecuapetla F, Galarraga E, Carrillo-Reid L. Diseño modular de instrumentación virtual para la manipulación y el análisis de señales eletrotisiológicas. Rev Mex Ing Biomed 27: 82–92, 2006.
Magill PJ, Bolam JP, Bevan MD. Dopamine regulates the impact of the cerebral cortex on the subthalamic nucleus-globus pallidus network. Neuroscience 106: 313–330, 2001.[CrossRef][Web of Science][Medline]
Mahon S, Vautrelle N, Pezard L, Slaght SJ, Deniau JM, Chouvet G, Charpier S. Distinct patterns of striatal medium spiny neuron activity during the natural sleep-wake cycle. J Neurosci 26: 12587–12595, 2006.
Mao BQ, Hamzei-Sichani F, Aronov D, Froemke RC, Yuste R. Dynamics of spontaneous activity in neocortical slices. Neuron 32: 883–898, 2001.[CrossRef][Web of Science][Medline]
Ossowska K. Interaction between striatal excitatory amino acid and gamma-aminobutyric acid (GABA) receptors in the turning behaviour of rats. Neurosci Lett 202: 57–60, 1995.[CrossRef][Web of Science][Medline]
Ossowska K, Wolfarth S. Stimulation of glutamate receptors in the intermediate/caudal striatum induces contralateral turning. Eur J Pharmacol 273: 89–97, 1995.[CrossRef][Web of Science][Medline]
Parthasarathy HB, Graybiel AM. Cortically driven immediate-early gene expression reflects modular influence of sensorimotor cortex on identified striatal neurons in the squirrel monkey. J Neurosci 17: 2477–2491, 1997.
Peterlin ZA, Kozloski J, Mao BQ, Tsiola A, Yuste R. Optical probing of neuronal circuits with calcium indicators. Proc Natl Acad Sci USA 97: 3619–3624, 2000.
Plenz D, Aertsen A. Neural dynamics in cortex-striatum co-cultures. II. Spatiotemporal characteristics of neuronal activity. Neuroscience 70: 893–924, 1996.[CrossRef][Web of Science][Medline]
Plenz D, Thiagarajan TC. The organizing principles of neuronal avalanches: cell assemblies in the cortex? Trends Neurosci 30: 101–110, 2007.[CrossRef][Web of Science][Medline]
Romo R, Scarnati E, Schultz W. Role of primate basal ganglia and frontal cortex in the internal generation of movements. II. Movement-related activity in the anterior striatum. Exp Brain Res 91: 385–395, 1992.[Web of Science][Medline]
Roweis S, Saul LK. Nonlinear dimensionality reduction by locally linear embedding. Science 290: 2323–2326, 2000.
Sasaki T, Kimura R, Tsukamoto M, Matsuki N, Ikegaya Y. Integrative spike dynamics of rat CA1 neurons: a multineuronal imaging study. J Physiol 574: 195–208, 2006.
Sasaki T, Matsuki N, Ikegaya Y. Metastability of active CA3 networks. J Neurosci 27: 517–528, 2007.
Schnitzler A, Gross J. Normal and pathological oscillatory communication in the brain. Nat Rev Neurosci 6: 285–296, 2005.[CrossRef][Web of Science][Medline]
Schreiber S, Fellous JM, Whitmer D, Tiesinga P, Sejnowski TJ. A new correlation-based measure of spike timing reliability. Neurocomputing 52–54: 925–931, 2003.
Schultz W, Apicella P, Ljungberg T, Romo R, Scarnati E. Reward-related activity in the monkey striatum and substantia nigra. Prog Brain Res 99: 227–235, 1993.[Web of Science][Medline]
Schwartz TH, Rabinowitz D, Unni V, Kumar VS, Smetters DK, Tsiola A, Yuste R. Networks of coactive neurons in developing layer 1. Neuron 20: 541–552, 1998.[CrossRef][Web of Science][Medline]
Stopfer M, Jayaraman V, Laurent G. Intensity versus identity coding in an olfactory system. Neuron 39: 991–1004, 2003.[CrossRef][Web of Science][Medline]
Takakusaki K, Oohinata-Sugimoto J, Saitoh K, Habaguchi T. Role of basal ganglia-brainstem systems in the control of postural muscle tone and locomotion. Prog Brain Res 143: 231–237, 2004a.[Web of Science][Medline]
Takakusaki K, Saitoh K, Harada H, Kashiwayanagi M. Role of basal ganglia-brain stem pathways in the control of motor behaviors. Neurosci Res 50: 137–151, 2004b.[CrossRef][Web of Science][Medline]
Tepper JM, Koos T, Wilson CJ. GABAergic microcircuits in the neostriatum. Trends Neurosci 27: 662–669, 2004.[CrossRef][Web of Science][Medline]
Tsodyks M. Attractor neural networks and spatial maps in hippocampus. Neuron 48: 168–169, 2005.[CrossRef][Web of Science][Medline]
Tunstall MJ, Oorschot DE, Kean A, Wickens JR. Inhibitory interactions between spiny projection neurons in the rat striatum. J Neurophysiol 88: 1263–1269, 2002.
Uhlhaas PJ, Singer W. Neural synchrony in brain disorders: relevance for cognitive dysfunctions and pathophysiology. Neuron 52: 155–168, 2006.[CrossRef][Web of Science][Medline]
Vergara R, Rick C, Hernandez-Lopez S, Laville JA, Guzman JN, Galarraga E, Surmeier DJ, Bargas J. Spontaneous voltage oscillations in striatal projection neurons in a rat corticostriatal slice. J Physiol 553: 169–182, 2003.
Tseng KY, O'Donnell P. Post-pubertal emergence of prefrontal cortical up states induced by D1-NMDA co-activation. Cereb Cortex 15: 49–57, 2005.
Wickens JR, Oorschot DE. Neural dynamics and surround inhibition in the neostriatum: a possible connection. In: Brain Dynamics and the Striatal Complex, edited by Miller R, Wickens JR. Australia: Harwood Acad, 2000, p. 141–149.
Wilson CJ. The generation of natural firing patterns in neostriatal neurons. Prog Brain Res 99: 277–297, 1993.[Web of Science][Medline]
Yuste R, MacLean JN, Smith J, Lansner A. The cortex as a central pattern generator. Nature Reviews 6: 477–483, 2005.
This article has been cited by other articles:
![]() |
L. Carrillo-Reid, F. Tecuapetla, N. Vautrelle, A. Hernandez, R. Vergara, E. Galarraga, and J. Bargas Muscarinic Enhancement of Persistent Sodium Current Synchronizes Striatal Medium Spiny Neurons J Neurophysiol, August 1, 2009; 102(2): 682 - 690. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. E. Rueda-Orozco, E. Mendoza, R. Hernandez, J. J. Aceves, O. Ibanez-Sandoval, E. Galarraga, and J. Bargas Diversity in long-term synaptic plasticity at inhibitory synapses of striatal spiny neurons Learn. Mem., July 24, 2009; 16(8): 474 - 478. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. M. Yu, J. P. Cunningham, G. Santhanam, S. I. Ryu, K. V. Shenoy, and M. Sahani Gaussian-Process Factor Analysis for Low-Dimensional Single-Trial Analysis of Neural Population Activity J Neurophysiol, July 1, 2009; 102(1): 614 - 635. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. J. Gruber, E. M. Powell, and P. O'Donnell Cortically Activated Interneurons Shape Spatial Aspects of Cortico-Accumbens Processing J Neurophysiol, April 1, 2009; 101(4): 1876 - 1882. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Carrillo-Reid, F. Tecuapetla, O. Ibanez-Sandoval, A. Hernandez-Cruz, E. Galarraga, and J. Bargas Activation of the Cholinergic System Endows Compositional Properties to Striatal Cell Assemblies J Neurophysiol, February 1, 2009; 101(2): 737 - 749. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. E. Pomata, M. A. Belluscio, L. A. Riquelme, and M. G. Murer NMDA Receptor Gating of Information Flow through the Striatum In Vivo J. Neurosci., December 10, 2008; 28(50): 13384 - 13389. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. R. Miller, A. G. Walker, A. S. Shah, S. J. Barton, and G. V. Rebec Dysregulated Information Processing by Medium Spiny Neurons in Striatum of Freely Behaving Mouse Models of Huntington's Disease J Neurophysiol, October 1, 2008; 100(4): 2205 - 2216. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Rayport Visualizing Striatal Networks. Focus on: "Encoding Network States by Striatal Cell Assemblies" J Neurophysiol, March 1, 2008; 99(3): 1053 - 1054. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |