Abstract
The subiculum, which has a strategic position in controlling hippocampal activity, is receiving significant attention in epilepsy research. However, the functional organization of subicular circuits remains unknown. Here, we combined different recording and analytical methods to study focal and widespread population activity in the isolated subiculum in zero Mg2+ media. Patch and field recordings were combined to examine the contribution of different cell types to population activity. The properties of cells leading field activity were examined. Predictive factors for a cell to behave as leader included exhibiting the bursting phenotype, displaying a low firing threshold, and having more distal apical dendrites. A subset of bursting cells constituted the first glutamatergic type that led a recruitment process that subsequently activated additional excitatory as well as inhibitory cells. This defined a sequence of synaptic excitation and inhibition that was studied by measuring the associated conductance changes and the evolution of the composite reversal potential. It is shown that inhibition was time-locked to excitation, which shunted excitatory inputs and suppressed firing during focal activity. This was recorded extracellularly as a multi-unit ensemble of active cells, the spatial boundaries of which were controlled by inhibition in contrast to widespread epileptiform activity. Focal activity was not dependent on the preparation or the developmental state because it was also recorded under 5 mm [K+]o and in adult tissue. Our data indicate that the subicular networks can be spontaneously organized as leader-follower local circuits in which excitation is mainly driven by a subset of bursting cells and inhibition controls spatiotemporal firing.
Introduction
Neuronal circuits produce different patterns of activity depending on input configuration and physiological conditions. Among these patterns, the one that results from the spontaneous correlated activity in the absence of external inputs is likely to reflect the intrinsic architecture of the underlying network. Examples include ensemble activity of cortical circuits (Mao et al., 2001; Cossart et al., 2003), low-frequency cortical oscillations (Steriade et al., 1993; Sanchez-Vives and McCormick, 2000), and persistent activity (Goldman-Rakic, 1995). In all these cases, there is an intimate relationship between the electrophysiological cellular properties, the connectivity patterns, and the correlated spatiotemporal activity they produce.
The subiculum, which has a strategic position in controlling hippocampal activity, is currently receiving much experimental attention. This interest derives from an emerging awareness of its role in spatial encoding (Sharp and Green, 1994; Taube, 1995; Ekstrom et al., 2003) and human epilepsy (Cohen et al., 2002). The subiculum is a single-layered network exhibiting a large electrophysiological diversity of neuronal types, including bursting and non-bursting glutamatergic cells and fast-spiking (FS) inter-neurons (Greene and Totterdell, 1997; Menendez de la Prida et al., 2003). The differential role of these cell types in processing is poorly understood, although it has been shown that bursting cells participate in sustaining hippocampal-generated gamma rhythms (Stanford et al., 1998) and epileptiform activity (Harris and Stewart, 2001; Cohen et al., 2002). Similarly, little is known about the organization of glutamatergic subicular networks, although a laminar and columnar architecture has been suggested (Harris et al., 2001). Inhibitory circuits, which are functionally recruited in a feedforward manner (Colino and De Molina, 1986; Finch et al., 1988), exert a control of bursting cell firing (Menendez de la Prida, 2003). However, despite this knowledge, a large experimental gap exists between what we know about subicular architecture and its function. The importance of this problem is twofold, first because the internal operation of the subiculum would shape its response properties during spatial orientation tasks (Sharp, 1999; Taube and Basset, 2003), and second because by understanding how the subiculum generates activity, the role of each cell type will be uncovered (Harris and Stewart, 2001; Kawaguchi, 2001).
To better understand the functional organization of subicular circuits, we examined the patterning of activity induced by bathing subicular minislices in zero Mg2+ media. Using this approach (Walther et al., 1986; Traub et al., 1994; Behr and Heinemann, 1996; Harris and Stewart, 2001), we could identify the cellular and synaptic initiators of population activity. We combined electrophysiological recordings at several levels [field potential, multi-unit activity (MUA), and patch recordings] with computational analysis of field-associated synaptic conductances, to build a comprehensive view of the operation of subicular circuits. Our results explored several aspects of neuronal processing such as: (1) the leading role of a subset of cells; (2) the role of shunting inhibition; and (3) the spatiotemporal patterning of population activity and its excitatory-inhibitory control. Our data indicate that the subicular networks can be spontaneously organized as leader-follower circuits in which excitation is mainly driven by a subset of bursting cells and inhibition controls spatiotemporal firing.
Materials and Methods
Minislice preparation, stimulation, and recordings. Horizontal slices (hippocampus/entorhinal cortex, 350 μm) were prepared from juvenile (17-20 d old) and adult (2 months old) Wistar and Sprague Dawley rats, as described previously (Menendez de la Prida et al., 2003). No differences were found in data from both strains. Extracellular solution [artificial CSF (ACSF)] contained (in mm): 125 NaCl, 3 KCl, 1 MgCl2, 1.2 NaH2PO4, 2 CaCl2, 22 NaHCO3, and 10 glucose, pH 7.3 when bubbled with 95% O2-5% CO2. Population activity was enhanced either by omitting MgCl2 from the ACSF or by increasing [K+]o to5mm. The subiculum was isolated from the hippocampus, and the entorhinal cortex using a dissection scissor. Subicular minislices of different sizes were obtained, and their width was measured using a calibrated eyepiece (Fig. 1 A). Slices were maintained in ACSF at room temperature for 1 hr before being transferred to a submerged-type chamber (flow rate, ∼1.5 ml/min; 32-34°C) attached to an upright microscope (Olympus, Tokyo, Japan).
Experimental preparation. A, Isolated subicular minislices of different sizes (widths) were prepared from juvenile and adult rats. The proximal (near CA1) and distal (near pre-subiculum) borders are shown. Simultaneous field (fp1 and fp2) and whole-cell (wc) recordings were performed. B, Epileptiform field events were simultaneously recorded at the proximal (fp1) and distal (fp2) borders. Epileptiform events are depicted at high (left) and low (right) temporal resolution. C, Pharmacological dependency of epileptiform activity. The number of slices is indicated. D, Consecutive epileptiform events exhibited a large variability, with no clear leading region. Two consecutive events from a representative minislice are shown. E, Cross-correlation analysis of dual field recordings from a minislice ∼1000 μm wide (data from 4 consecutive events). Interelectrode distance, ∼300 μm. F, Frequency of epileptiform events versus minislice size (width). Data are from n = 29 minislices. G, Cross-correlation analysis of dual recordings (4 consecutive events) from the minislice shown in E, after re-sectioning (∼650 μm wide). Same interelectrode distance before and after cutting, ∼300 μm.
Extracellular field potentials were recorded using patch-clamp pipettes (borosilicate glass capillaries; outer diameter, 1.2 mm; inner diameter, 0.69 mm; Harvard Apparatus) filled with ACSF (tip resistance, ∼7 MΩ). Signals were amplified and bandpass filtered at 0.1 Hz to 5 kHz or at 0.1-300 Hz. MUA was recorded using a coated tungsten wire of 100 μm diameter etched to a point of 10-15 μm diameter. The tip was then re-coated using an insulating resin to obtain a final resistance of 0.7-0.8 MΩ (uninsulated tip length, ∼100-150 μm). MUA signals were amplified using two independent low-noise amplifiers (P15; Grass Instruments, Quincy, MA) and bandpass filtered at 100 Hz to 5 kHz. For paired MUA recordings, the electrodes were separated >200 μm to avoid over-sampling the same cellular population (Cohen and Miles, 2000). At the end of the recording session, a mark was made in the tissue and the interelectrode distance was measured.
Somatic patch-clamp recordings were made under visual control (60× immersion lens) using an Axoclamp 2B amplifier (Axon Instruments, Foster City, CA). Patch recordings were performed using both the cell-attached configuration in the voltage-clamp mode and the whole-cell configuration in the current- and voltage-clamp modes. Patch pipettes were filled with intracellular solution containing (in mm) 113 K-gluconate, 6 KCl, 1 MgCl2, 1 NaCl, 1 EGTA, 5 HEPES, 2 K2ATP, 0.3 NaGTP, and Neurobiotin (0.1-0.5%; Vector Laboratories, Burlingame, CA), pH 7.3 adjusted with KOH and osmolarity between 290 and 300 mOsm. For some voltage-clamp experiments, the following solution was used (in mm): 115 Cs-gluconate, 4 MgCl2, 1 EGTA, 5 HEPES, 2 Na2ATP, and 0.3 NaGTP, pH 7.3 adjusted with CsOH and osmolarity between 290 and 300 mOsm. Patch pipettes had a resistance of 4-6 MΩ. Capacitance compensation and bridge balance were performed before patching. The access resistance was calculated off-line from the cell response to subthreshold current pulses (double exponential fitting). The junction potential (∼-7 mV for K-based solution and ∼-10 mV for Cs-based solution) was corrected (Neher, 1994). All recordings were simultaneously digitized (Digidata 1200B; Axon Instruments) and stored on disk at a sampling frequency >10 kHz. Monopolar stimulation (square-wave; 0.1 msec duration; 0.5-1 mA amplitude) was delivered via tungsten electrodes (tip diameter, 10-20 μm).
Morphological reconstruction and analysis. After each experiment, the minislice was fixed overnight in 4% paraformaldehyde PBS (0.1 m; pH 7.4). Because of their small size, minislices were not dehydrated or resliced. After H2O2 (0.3%) and Triton X-100 (0.6%) pretreatment, minislices were processed by incubation in a 1:100 dilution of an ABC kit (Vector Laboratories) for 3 hr. Visualization was performed using a solution of 0.03% of 3,3-diaminobenzidine and 0.005% H2O2. The reaction was followed under visual control and stopped after 5-10 min depending on labeling intensity. Minislices were washed in TBS (0.5 m) and mounted on gelatin-coated glass slides using glycerol (50% in TBS) and Eukitt mounting medium (Fluka, Neu-Ulm, Germany). More than 60% of Neurobiotin-filled cells were successfully recovered. Cells were drawn using 40× and 60× (dendritic tree) and 100× (axon collaterals) objectives and a camera lucida attachment. The dendritic arborization was quantified using the Scholl concentric ring analysis (Scholl, 1953). For this purpose, soma-centered concentric rings of increasing diameters of 25 μm were drawn on the Camera Lucida reconstruction, and the number of dendritic crossings was counted (basal, negative numbers; apical, positive numbers). Axonal boutons were identified using a 100× oil immersion objective. No difference was found in the density of axonal boutons among different glutamatergic cell types. However, we cannot totally exclude this possibility because of a possible underestimation attributable to minislice thickness.
Electrophysiological classification of cell types. Subicular cells were electrophysiologically classified as weak (IB-) and strong (IB+) bursting, regular-spiking (RS) and FS, from their responses to depolarizing current pulses of 500 msec duration (Stewart and Wong, 1993; Staff et al., 2000; Menendez de la Prida et al., 2002, 2003). Intrinsic properties were determined in Mg2+-free solution. Bursting cells firing only one burst in response to several suprathreshold current pulses were classified as IB-cells, whereas those firing more than one burst were classified as IB+. Previous studies have shown that bursting and RS neurons are projection glutamatergic cells, whereas the FS type corresponds to inhibitory inter-neurons (Harris et al., 2001; Menendez de la Prida, 2003; Menendez de la Prida et al., 2003).
MUA activity: spike detection and firing rate data. Multiple spikes recorded using tungsten electrodes were detected off-line (Cohen and Miles, 2000). A Gaussian convolution was used to derive firing rate data (SD for convolution, 20 msec). In MUA recordings, collision of spikes may occur. This situation is more frequent while recording population activity that results from coordinated firing of several cells. In preliminary simulations, we estimated that pairs of spikes separated by intervals of 0.6-0.8 msec (1250-1666 Hz) could be resolved successfully.
Estimation of input conductance and composite reversal potential. Input conductance increases (ΔG) during focal field activity were measured using routines written in Matlab (MathWorks, Natick, MA) following the method of Borg-Graham et al. (1998). Conductance changes during the buildup period before epileptiform events could not be efficiently measured and are not discussed. Spontaneous synaptic currents associated with focal field activity were recorded in voltage clamp at three different holding potentials between -55 and -75 mV (transmembrane potential changes during field activity, <1 mV). The holding potential was corrected for the drop across access resistance, which was lower than 30 MΩ [18 ± 7MΩ (n = 19) and 19 ± 9MΩ (n = 6) for K-gluconate and Cs-gluconate electrodes, respectively]. For each holding potential, four to six field-triggered spontaneous increases of synaptic activity were recorded and averaged. G(t) was then derived from the slope of the linear regression at the I-V plot. The regression coefficient (r) ranged from 0.98 to 0.99 (see examples in supplementary Fig. 1; available at www.jneurosci.org), which suggest good voltage-clamp conditions and that voltage-gated channels were not typically activated during voltage-clamp experiments. The relative conductance increase ΔG(t) was given as [G(t)-Grest]/Grest in percentage, where Grest is the resting conductance. Similarly, the composite reversal potential of the synaptic input Vrev(t) was calculated as the voltage of the intersection between the I-V curve at time t and at rest. Vrev(t) is also known as the null-current potential. Because the intersection method is noisy for similar slopes, we only considered those cases in which synaptic conductance was >5% than resting conductance.
To dissect the contribution of glutamatergic and GABAergic synaptic currents to the composite Vrev(t), we estimated the reversal potential of evoked EPSCs and IPSCs. CNQX- and AP-5-sensitive glutamatergic currents reversed at -6 ± 4 mV (n = 4; see supplementary Fig. 1). The reversal potential of GABAA-mediated currents was -72 ± 2 mV (n = 4), which coincides with the theoretical value predicted for our pipette solution. Because a use-dependent shift of Vrev has been described for these currents (Thompson and Gähwiler, 1989; Staley and Proctor, 1999), we also estimated changes induced by repetitive activation of GABAA receptors. Repetitive stimulation of 25-100 Hz at 250-400 msec shifted the GABAA-mediated reversal potential to -64 ± 4 mV (n = 4; p < 0.05), as measured 40 msec after the train outset (see supplementary Fig. 1). As described previously (Thompson and Gähwiler, 1989; Staley and Proctor, 1999), this shift was associated with an increase of intracellular Cl- concentration, because it was induced at holding potentials more positive than the Cl- equilibrium potential (i.e., holding potential, -50 mV) but not at negative holding potentials (-80 mV). The reversal potential of GABAB-mediated currents was -88 ± 2mV(n = 5; data not shown).
We failed to obtain confident estimates of Gexc(t) and Ginh(t) using methods described previously (Borg-Graham, 2001). The reason for that derived from the need to place constraints to solve algebraic equations, which substantially distorted results. This is especially critical when more than two components (NMDA, non-NMDA, GABAA, and GABAB) are contributing together and in the case when GABAA-mediated currents exhibit a use-dependent shift of Vrev. Another potential limitation arises from possible differences in the spatial distribution of excitatory and inhibitory inputs. Therefore, we used the temporal evolution of the Vrev(t) to characterize the differential contribution of excitation and inhibition.
Data analysis. Field potential and MUA firing rate data were analyzed using routines written in Matlab. Normalized cross-correlation functions were computed using fast Fourier transform algorithms. Latencies between dual field recordings were calculated from the peak offset of cross-correlation. The irregular nature of population activity was quantified by calculating the coefficient of variation (CV = SD/mean) of a given measure.
The intrinsic firing pattern and membrane properties of cells were tested in the whole-cell current-clamp mode, as described previously (Menendez de la Prida, 2003). Input resistance was determined from cell responses to hyperpolarizing current pulses of 500 msec duration. To determine the relative [relative to resting membrane potential (RMP)] and absolute firing threshold, depolarizing pulses of increasing amplitude (50 msec duration) were injected at RMP. Firing threshold was defined at the inflection point of the voltage trace.
To examine the temporal correlation between single-cell firing and population activity, field-triggered data segments were used. For this purpose, 8-13 data segments were triggered by the first negative field deflection detected over noise level (∼20 μV). The mean amplitude of triggered data segments from field recordings had a CV that lies within 5-15%. The latency to the first spike recorded in the cell-attached mode was calculated, and histograms for each cell were obtained. To evaluate the role of distinct cell types, histograms from similar electrophysiological categories (IB+, IB-, RS, and FS) were pooled. Histograms from individual cells showed variability because of their irregular firing pattern. It was therefore important to statistically average the contribution of each cell over several data segments. For this purpose, we introduced the driving index (σ) as a measure of the probability of a given cell to start firing before the field:
A driving index σ = 0 identified those cells that fire after the field in all data segments (followers), whereas a driving index σ = 1 identified those cells that only fire before the field (leaders). A driving index σ = 0.5 means that it is not possible to statistically define a leader/follower firing pattern for that cell.
Results are given as means ± SD, with the number of minislices (or cells) indicated in every case. Results were compared using Student's t test or the ANOVA test. Simple linear regression was performed with a confidence interval of 95%.
Drugs. Picrotoxin (PTX), 2,3-dioxo-6-nitro-1,2,3,4-tetrahydro-benzo[f]quinoxaline-7-sulfonamide (NBQX), CNQX, AP-5, atropine, mecamylamine, and salts were purchased from Sigma (St. Louis, MO). Hydroxy-saclofen, α-methyl-4-carboxyphenylglycine (MCPG), and α-methyl-4-phosphonophenylglycine (MPPG) were obtained from Tocris Cookson.
Results
Temporal dynamics of population activity in the isolated subiculum in vitro
Lowering extracellular Mg2+ induced spontaneous epileptiform activity in the isolated subiculum after 10-20 min (∼70% minislices) (Behr and Heinemann, 1996; Harris and Stewart, 2001). Recordings were performed once activity stabilized, between 30 min and 2 hr after washout of Mg2+. Epileptiform activity was recorded at both the proximal (near CA1) and distal (near pre-subiculum) borders (Fig. 1A,B), indicating the widespread activation of all subicular cells. Local field activity was also detected (see below). The intervals between epileptiform events were widely irregular (mean frequency, 0.26 ± 0.19 Hz; CV, 48 ± 15%; n = 13) (Fig. 1B, right traces). The GABAA receptor antagonist PTX (100 μm) increased 1.3 ± 0.4 times the frequency of epileptiform activity (0.21 ± 0.07 Hz; CV, 46 ± 23%; n = 9), which was blocked by AP-5 (100 μm; n = 8) and strongly reduced by NBQX (n = 2) or CNQX (n = 6) at 10 μm (Fig. 1C).
To investigate whether there is a clear leading region of epileptiform activity, we performed cross-correlation analysis of dual field recordings. We looked for the appropriate slice width (Fig. 1A) for which epileptiform activity was spontaneously generated with minimal temporal variations. Minislices 900-1300 μm wide showed large variability of latencies between dual field recordings separated ∼300 μm (CV, 76 ± 33%) (Fig. 1E). Epileptiform activity was not recorded in minislices <500 μm wide (n = 9) (Fig. 1F). Minislices 500-750 μm wide (n = 15) were less excitable than minislices >750 μm wide (n = 15) (Fig. 1F). Interestingly, latency variations were significantly reduced in smaller minislices, and mean latency significantly decreased to -1.4 ± 1.9 msec (versus -14.6 ± 4.4 before cutting; same interelectrode distance before and after cutting; n = 4; p < 0.005) (Fig. 1G).
Therefore, similar to the hippocampus (Miles et al., 1984), the site of origin of population activity in the subiculum is spatially constrained in minislices 500-750 μm wide (latency, ∼0 msec). This also guaranteed that latencies are conserved among different episodes of field activity. We used this approach to better understand how different cell types interact to generate field activity in the isolated subiculum.
Focal and widespread field activity in the isolated subiculum
To examine the relationship between single-cell firing and population activity, cell-attached recordings were initially obtained. Whole-cell configuration in current and voltage clamp was then established to characterize membrane properties and synaptic activity. We took special care to avoid artificial delays between patch and field recordings by minimizing their relative distance (∼20-50 μm). The maximal error introduced was estimated from the delay between the cell and field spike evoked by local stimulation (0.1 ± 0.8 msec; n = 4 cells). A total of 54 cells were recorded: 17 were classified as IB+, 17 were IB-, 12 were RS, and 8 were of the FS type (Table 1; see Materials and Methods).
Passive and active properties of the different cell types
We classified field activity as desynchronized, widespread, and focal (Fig. 2). During periods of desynchronized activity, irregular cell firing was present mostly in FS and bursting cells (Fig. 2A; Table 1). The nature of this firing could be intrinsic or synaptic, because spontaneous EPSCs were evident in voltage-clamp configuration (Fig. 2A, vc trace) and EPSP-driven firing was present in current-clamp configuration (Fig. 2A, cc trace). Periods of desynchronized field activity were characterized by a flat cross-correlation of dual field recordings separated 450 ± 69 μm (n = 7) (Fig. 2D).
Simultaneous field (fp) and patch recordings. Cell responses were examined using cell-attached (ca) and whole-cell (wc) configurations in the current-clamp (cc) and voltage-clamp (vc) modes. Field recordings corresponding to ca, cc, and vctraces are depicted. A, Simultaneous field and patch recordings during desynchronized field activity. This cell was classified as weak bursting (IB-). B, Simultaneous field and patch recordings during widespread field activity. C, Simultaneous field and patch recordings during focal field activity. D, Cross-correlation analysis of simultaneous field recordings. Results are shown from three consecutive episodes of desynchronized field activity. E, Widespread activity. F, Focal activity. Insets show examples of dual field recordings (fp1 and fp2).
During widespread population activity (Fig. 2E), cell firing was recorded simultaneously with epileptiform field events (Fig. 2B). Interestingly, some cells (mostly bursting cells) increased their firing 50-100 msec before the onset of field events. Recordings under current and voltage clamp showed that this was associated with an increase of synaptic activity (Fig. 2B, cc and vc traces). This buildup period has been described in other in vitro models of epilepsy (Chamberlin et al., 1990; Traub and Dingledine, 1990) and in immature population activity (Menendez de la Prida and Sanchez-Andres, 1999; Wenner and O'Donovan, 2001).
The combination of cell-attached and field recordings facilitated the identification of focal activity (Fig. 2C), which was characterized by low-amplitude changes in the field that could be associated with cell firing. The mean, but not the variance, of the distribution of field event amplitude was different in focal (74 ± 28 μV) versus widespread episodes (131 ± 54 μV; data from 21 events from n = 7 minislices; p = 0.05). Once in whole-cell configuration, barrages of synaptic activity were evident in current- and voltage-clamp simultaneous to focal field activity (Fig. 2C, cc and vc traces). Not all the cells recorded exhibited firing in correlation with focal activity. Those that did included 5 of 17 IB+ cells, 5 of 17 IB-cells, 5 of 12 RS cells, and 6 of 8 FS interneurons. This suggests that only a fraction of cells contribute to the averaged field changes (see below). Focal activity occurs erratically, with no clear periodicity (mean frequency, 0.36 ± 0.16 Hz; CV, 58 ± 18%; n = 22 minislices). The absence of peaks in the auto-correlation (data not shown) and cross-correlation functions (Fig. 2F) from dual field recordings confirmed its local non-oscillatory nature and suggested that focal activity was not related with waves of activity across the tissue.
Role of bursting and FS cells in the generation of focal activity
We first examined the contribution of each cell type to focal activity using field-triggered data segments from cell-attached recordings (see Materials and Methods). Both IB+ (n = 5) and IB-(n = 5) cells that participated in focal activity typically began to fire before the field activity [-13 ± 35 msec (n = 41 episodes from five IB+ cells) and -12 ± 41 msec (n = 40 episodes from five IB+ cells)] (Fig. 3A,B). In contrast, RS cells were activated after the field activity (12 ± 25 msec; n = 40 episodes from five cells) (Fig. 3A,B), and a same cell could remain silent during some episodes (Fig. 3A, RS cell). In addition, bursting cells typically fired bursts, whereas RS cells fired single spikes (Fig. 3A). We also examined the behavior of FS interneurons. Similar to RS cells, FS interneurons began to fire 10-12 msec after focal field activity (11 ± 15 msec; n = 64 episodes from six cells) (Fig. 3B), but in contrast to the glutamatergic cell type, this was followed by high-frequency firing (Fig. 3A,C).
Cellular and synaptic behavior underlying focal activity. A, Cell responses in cell-attached (ca) configuration during focal activity were examined using field-triggered data segments (fp trace). Two consecutive episodes are shown for each cell (RS, IB+, IB-, and FS cells). B, Histograms of cell firing from each cell type. C, Representative response from another FS interneuron (ca trace) triggered by the onset of focal field activity (fp trace). Note the different time scale compared with A. D, Field-triggered synaptic currents from an IB-cell. E, Field-triggered synaptic currents from the cell shown in D at three different holding potentials (-58, -80, and -95 mV). F, Field-triggered synaptic currents from a different IB-cell at two holding potentials. G, Field-triggered composite reversal potential from cells shown in E (black trace) and F (gray trace). E-G share the same time scale (shown in G).
These data suggest that both IB+ and IB-cells constituted the first glutamatergic type that initiates a focal recruitment process that subsequently activates RS and FS cells. FS interneurons would in turn suppress activity. To test this idea, we looked at the field-triggered synaptic changes in 20 cells (5 IB+, 5 IB-, 5 RS, and 5 FS) and examined the temporal correlation between excitation and inhibition using the composite reversal potential (see Materials and Methods).
Glutamatergic- and GABAergic-mediated synaptic activity during focal activity
In most cells (12 of 19: 4 IB+, 4 IB-, and 4 RS), excitatory synaptic currents were recorded at the onset of focal field activity (Fig. 3D,E). The composite Vrev(t) typically showed a depolarized shift during the first 50-70 msec (peak value, -31 ± 11 mV; n = 12) (Fig. 3G, black trace), which suggests a dominance of glutamatergic excitatory currents. No differences were found between IB+ (-28 ± 17 mV; n = 4), IB-(-34 ± 17 mV; n = 4), and RS cells (-27 ± 5 mV; n = 4). All FS cells examined also exhibited a peak value of Vrev(t) that was even more depolarized than in glutamatergic cells (-7 ± 5 mV; n = 5; p < 0.005). Because both IB+ and IB-cells were the major glutamatergic type activated before the onset of focal activity, this analysis supports their role in the recruitment process. In addition, it also supports that FS cells experience strong glutamatergic activation at the onset of field activity.
Some bursting (IB+, one of five; IB-, one of five) and RS cells (one of five) exhibited inhibitory synaptic currents at the onset of focal activity (Fig. 3F, current trace at -58 mV holding potential). In these cases, the composite Vrev(t) showed a drop toward negative potentials ∼-75 mV, near the reversal potential of control GABAA-mediated chloride currents (Fig. 3F,G, gray trace). Because most FS inter-neurons began to fire 12 msec after the field (Fig. 3B), this probably reflects the effect of a fast inhibitory surround. A second barrage of excitatory synaptic activity was evident in these cells (Fig. 3F), as showed by the shift of Vrev(t) toward positive values (Fig. 3G, gray trace). This glutamatergic component could be related to the secondary recruitment of RS cells or with sustained activity of bursting cells (Fig. 3B).
In some neurons (5 of 20), long-lasting synaptic currents were recorded (Fig. 3D, gray shading). Their GABAergic nature was suggested by comparison with the temporal course of interneuron firing during field-triggered data (Fig. 3A,C,D) and by its reversal potential (Fig. 3E). The composite Vrev(t) showed a negative shift 150-200 msec after the onset of focal activity (-88 ± 6 mV; n = 5) (Fig. 3G, black trace). This is consistent with a late prevalence of GABAB-mediated currents.
Role of different cell types in the initiation of widespread epileptiform activity
We also examined the contribution of each cell type to the onset of widespread epileptiform activity using field-triggered averages. We found that both IB+ and IB-cells began to fire 50-100 msec before the onset of widespread epileptiform activity (Fig. 4A,B). In contrast, RS cells tended to fire after the onset of epileptiform activity (Fig. 4C). Most FS interneurons fired at the onset, although some firing increases were detected before (Fig. 4D).
Cell behavior underlying widespread activity. Field-triggered (fp) cell responses (5 consecutive events are depicted) in a cell-attached configuration (ca) from representative IB+ (A), IB-(B), RS (C), and FS (D) cells. Firing histograms from all cells in each group are shown at the bottom.
A second observation was evident when comparing histograms of single-cell firing during widespread and focal activity. During focal activity, bursting cells had a peak of firing before the field activity and subsequently decreased firing as the RS and FS firing increased (Fig. 3B). In contrast, generation of widespread activity was associated with a preceding monotonic increase of bursting cell firing that reached a peak at the onset of the widespread activity (Fig. 4A,B). Both RS and FS cells had similar spike time histograms in focal versus widespread field activity. This suggests that bursting cells are critical in initiating widespread forms of population activity and that FS cell firing, which is time-locked to leading cell firing, plays a role in controlling runaway excitation.
However, the temporal correlation between cell firing and the field activity varied among different episodes of focal and widespread activity. For a given cell, firing could occur before or after the onset of field activity (Fig. 5A; see also Fig. 3A, IB+ cell), suggesting that they variably participated in correlated firing. We defined the driving index (σ) as a statistical measure of the probability of a given cell to fire before the field activity by examining its behavior in 10-13 consecutive episodes of population activity (see Materials and Methods).
Electrophysiological properties and the driving index. A, Example of firing histograms from one cell (IB+ cell) in cell-attached recordings during field-triggered activity (data from 13 consecutive episodes). The driving index measures the probability of a given cell to fire before the field (see Materials and Methods). B, Statistically significant correlation (p < 0.05; r =-0.58) between the firing threshold and the driving index (n = 46 glutamatergic cells). C, Field-triggered averages (black traces) of synaptic currents during focal activity at three holding potentials. SDs are shown in gray. D, Conductance changes (up) and composite reversal potential (bottom) computed from data in C. E, Vrev-ΔG phase plots (see Materials and Methods and Borg-Graham et al. 1998). ΔG(t) is plotted as a function of the Vrev(t) during the first milliseconds before reaching the conductance peak. Arrows mark this temporal evolution. Gray circles denote the Vrev,ΔG values at conductance peak. Open circles denote the Vrev,ΔG values at the Vrev peak. Note the shunting effect of peak conductance increases (gray dots).
Both IB+ and IB-bursting cells exhibited a higher driving index than RS cells (Table 1) (p < 0.05). Nevertheless, the mean σ from bursting cells was <0.5, a statistical limit separating the leader and follower firing pattern (see Materials and Methods). This suggests that in statistical terms, bursting cells did not define a glutamatergic leader population (although 15 of 34 cells had σ > 0.5). The RS group, in contrast, behaved as a follower (mean σ = 0.08 ± 0.13; n = 12). All but one FS interneuron (seven of eight) had a driving index < 0.5, indicating that they were typically activated afterward.
The driving capacity of a cell is mainly determined by electrophysiological and morphological factors
Intrinsic electrophysiological properties such as resting membrane potential, input resistance, membrane time constant, and firing threshold could affect the firing capacity of a cell on synaptic activation. We examined these electrophysiological parameters among different cell types and found the firing threshold to be the most critical factor (Fig. 5B) (correlation coefficient r = -0.58; p < 0.05). Interestingly, bursting cells had a lower firing threshold (relative) than RS cells (Table 1). A higher input resistance in RS cells would favor their synaptic recruitment, but these cells exhibited the highest firing threshold. FS cells had high input resistance and a lower firing threshold than RS cells.
To address the possibility that those cells displaying σ > 0.5 received stronger synaptic activation, we examined the temporal evolution of synaptic conductance changes during field-triggered focal activity in a subset of cells (n = 17; see Materials and Methods) (Fig. 5C,D). The composite Vrev(t) initially evolved toward depolarized values and then hyperpolarized as the conductance reached a peak (Fig. 5D). This evolution is better depicted in Fig. 5E, in which ΔG(t) is plotted as a function of the Vrev(t) during the first milliseconds before reaching the conductance peak (arrows). From these plots, clearly the conductance increases associated with the Vrev peak (open circles; -27 ± 15 mV) were lower than the conductance peaks (gray circles; p < 0.055), which are associated with a Vrev of -56 ± 8 mV. This suggests that during conductance peaks concurrent contribution of inhibition and excitation existed. Values of conductance peaks and conductance increases at Vrev peak were similar in the three glutamatergic cell types (Table 2). FS interneurons were subject to larger increases than RS and bursting cells (Table 2) (p < 0.005). We did not find a good correlation between conductance increases in glutamatergic cells and their driving index (r = 0.13; p > 0.05). Additional experiments were performed using a cesium-based patch solution (n = 6). The driving index of these cells (ranging from 0 to 1) was established in the cell-attached configuration. No correlation between conductance changes and the driving index was documented (r = 0.12; p > 0.05).
Conductance changes associated with focal field activity
The subtle depolarized Vrev at peak conductance suggests a shunting GABAergic effect over excitatory inputs (Fig. 5D). A use-dependent shift of the GABAA reversal potential could also contribute (see Materials and Methods). Therefore, it could be possible that proximal GABAergic inputs mask distal excitatory ones (Koch et al., 1990). Space-clamp errors of distal inputs could also contaminate these results (Spruston et al., 1993). We therefore explored whether differences in apical/basal dendritic arborization were correlated with the driving index in a group of morphologically reconstructed cells (n = 15).
As reported previously, no difference in dendritic morphology was found between bursting and RS subicular cells (Staff et al., 2000; Harris et al., 2001). However, when cells were grouped by the driving index (σ = 0, n = 5 cells; σ > 0.5, n = 4 cells), the difference in the number of distal dendritic crossings was evident at 250-450 μm from the soma (Fig. 6A, arrow). The mean number of distal crossings from all cells showed a significant correlation with the driving index (Fig. 6C) (r = 0.76; p < 0.05). No difference in the spatial distribution of cells was present (Fig. 6D). Therefore, the difference in distal crossings suggests that driving cells may receive more inputs at secondary distal dendrites than follower cells. The good linear regression of I-V data (regression coefficient, 0.98-0.99) suggests good voltage-clamp condition and that voltage-gated channels were not activated. However, although reliable estimates of the Vrev of synaptic currents located 250-400 μm from the soma can be obtained (2-3 mV error), conductance values are significantly reduced by >20% because of space-clamp errors (Spruston et al., 1993). Therefore, the lack of correlation with somatic conductance changes could be related to significant electrotonic filtering.
Morphological properties and the driving index. A, Scholl diagram of apical and basal dendritic arbors of cells exhibiting low (σ = 0.0 ± 0.0; n = 5) and high (σ = 0.7 ± 0.2; n = 4) driving index. Statistical difference was evident at distal apical dendritic crossings (arrow). B, Representative examples of cells in each category. C, Correlation between the mean apical dendritic crossings at 250-450 μm from the soma and the driving index in all cells examined (n = 15). D, Spatial distribution of cells in each category. Positions are normalized.
MUA recordings during desynchronized, focal, and widespread population activity
Together, these results suggest that although some bursting cells play a leader role, network activation is a probabilistic process that needs the local recruitment of several cell types. Periods of focal activity seem to reflect the functional architecture of the underlying network with bursting, RS and FS cells differentially sculpting the spatiotemporal patterning. To gain in spatiotemporal resolution of the network activity, we combined field and multi-unit recordings. We used pharmacological tools to dissect the underlying circuit mechanism.
Data collected using MUA recordings were consistent with the picture that emerged from cell-attached recordings. Desynchronized, widespread, and focal field activities were correlated with three different MUA patterns (Fig. 7). During desynchronized field states, multi-unit firing was recorded (Fig. 7A). Both single spikes and bursts (Fig. 7A, inset a) of different amplitudes were erratically detected, suggesting the spontaneous activation of different cells. During widespread epileptiform activity, high-frequency firing was recorded in correlation with large field epileptiform events (Fig. 7B). Interestingly, firing increases were detected 50-100 msec before the onset of field events, which matched the buildup period. Cell ensemble activity was evident in MUA recordings during focal field activity (Fig. 7C). The amplitude, number, and sequence of spikes recorded within these ensembles varied substantially, confirming that many cells were variable recruited during consecutive episodes. This was also manifested in the variability of ensemble duration (124 ± 34 msec; CV, 57 ± 16%; n = 11 minislices). Ensemble activity during focal activity was similar to MUA activity detected before the onset of widespread field activity (Fig. 7, compare B, C).
Simultaneous field and MUA recordings during desynchronized (desynch; A), widespread (B), and focal (C) field activity. Insets a-c show additional details of firing activity.
Pharmacological dissection of synaptic factors contributing to focal field activity
We took advantage of the higher spatiotemporal resolution of MUA recordings to pharmacologically dissect the synaptic contributions to focal field activity. For this purpose, we initially addressed the question of what are the sources of cell firing during desynchronized activity by examining changes of firing rate under different receptor antagonists (Fig. 8A). We first derived confidence intervals of a firing rate from MUA recordings of 1 hr duration. Changes of 14% were not statistically significant (n = 5 minislices).
Pharmacological control of MUA firing during desynchronized activity. A, Representative MUA traces under control, zero Mg2+, and zero Mg2+ plus CNQX/AP5/PTX. B, Summary of results from pharmacological control of MUA firing in juvenile (left) and adult (right) rats. The number of slices is indicated. Statistical differences are shown: **p < 0.001; +p < 0.05.
Lowering extracellular Mg2+ induced a subtle increase of firing during desynchronized field activity (n = 7) (Fig. 8B). Most of this firing had a glutamatergic origin because it was reduced to 29 ± 18% of control by suppressing glutamatergic transmission with CNQX (or NBQX) and AP-5 (n = 6) (Fig. 8B). This was confirmed by adding a mixture of antagonists containing MCPG and MPPG (1 mm), atropine (1 μm), mecamylamine (20 μm), and saclofen (200 μm; n = 4) (Fig. 8B).
We next wondered about the synaptic components of focal field activity. Data from patch recordings suggested an initial glutamatergic activation, followed by a composite GABAergic component. Using pharmacological tools, we confirmed that NMDA, but not non-NMDA, transmission was critical in the generation of focal activity that was blocked by AP-5 (Fig. 9A,C).
Pharmacological dissection of synaptic factors mediating focal activity. A, Field-triggered MUA ensembles during focal activity from juvenile tissue in zero Mg2+ media and after bath application of PTX, CNQX, and AP-5. B, Field-triggered MUA ensembles during focal activity from adult tissue in zero Mg2+ media. The inset shows additional details of firing activity. C, Summary of pharmacological results. The number of slices is indicated. D, Nissl-stained subicular minislice. A mark denoted the position of tungsten electrodes. E, Simultaneous MUA recordings during focal activity. F, Cross-correlation analysis of simultaneously recorded MUA firing rate. Results from five consecutive episodes are shown in gray; the mean is shown in black. G, Same as in F after 4-7 min of PTX application. Note the significant increase of cross-correlation. H, Summary of the cross-correlation analysis under different conditions: zero Mg2+, PTX (4-7 min), saclofen, and CNQX. The number of slices is indicated. I, Pharmacological dependence of the frequency of focal and widespread field events.
Consistent with a role of inhibition in firing cessation, both saclofen and PTX induced an increase of ensemble duration (n = 5; p < 0.05) (Fig. 9C). Interestingly, this effect was also present under CNQX and NBQX (Fig. 9C)(p < 0.05). The prolongation of ensemble activity under PTX was less pronounced than under saclofen and CNQX (Fig. 9C). To better understand the network mechanisms, we analyzed data from dual MUA recordings (Fig. 9D).
Correlated ensembles were detected in <35% cases from paired MUA recordings separated 285 ± 30 μm (n = 9) (Fig. 9E). PTX significantly favored the spread of ensemble activity for 4-7 min after application (Fig. 9F-H). After ∼7 min, local activity became unstable, and widespread epileptiform events were more frequent (Fig. 9I). The spread of local activity was also favored by saclofen (Fig. 9H) (p < 0.05), but in contrast to PTX, this did not change neither focal nor widespread event frequency (Fig. 9I). CNQX (or NBQX) did not affect the local spread of activity, although it substantially increased ensemble duration (Fig. 9C,H). Because non-NMDA receptors mediate activation of FS interneurons (Menendez de la Prida, 2003), this probably reflected a disinhibitory effect as opposed to firing cessation. CNQX did reduce the frequency of both local ensembles and epileptiform events (Fig. 9I). Together, these results suggest that the spatial specificity of focal activity is strongly determined by local inhibition because impairment of GABAA,B synaptic transmission favored merging of formerly non-overlapping ensembles.
Model and developmental dependence of focal activity
Finally, we tested whether focal activity was specific of the zero Mg2+ model or whether it could be induced by other manipulations such as increasing cell excitability using high K+ media. Increasing [K+]o to5mm induced focal ensemble activity in four of five minislices examined (data not shown). Focal events in 5 mm K+ were less frequent (0.08 ± 0.06 Hz; n = 4) than in the zero Mg2+ model (p < 0.05).
We also examined tissue from adult rats to discriminate possible developmental changes. Focal (Fig. 9B) and widespread (data not shown) population activity were also recorded in the adult isolated subiculum in zero Mg2+ (n = 11). Although the frequency of these events, and the duration of focal activity, was lower than in juvenile tissue, their pharmacology was similar (Fig. 9C). This probably reflects the lower susceptibility of adult tissue to epilepsy. Interestingly, the synaptic component of single-cell firing during desynchronized states was lower in adult rats than in juvenile rats (Fig. 8B). This further suggests developmental differences in local connectivity between juvenile and adult tissue.
Discussion
The present study combines different recording and analytical methods to dissect the functional organization of rat subicular circuits using the zero Mg2+ model. The main results showed that focal and widespread field activity in the isolated subiculum is regulated by the firing correlation between leader and follower cells. Predictive factors for a cell to lead field activity include exhibiting the bursting phenotype, displaying a low firing threshold, and having more distal apical dendrites. A group of bursting cells constitutes the first glutamatergic type that initiates a local recruitment process that subsequently activates follower excitatory and inhibitory cells. This is recorded extracellularly as a multi-unit local ensemble. The underlying synaptic events define a sequence of excitation and inhibition that controls the spatiotemporal firing. We found that with this functional design, inhibition mediates firing cessation by a feedback mechanism that is time-locked to the firing of leader cells. There is also a fast inhibitory control of firing propagation that regulates transition from focal to widespread activity.
Functional organization of local subicular circuits
To better understand the local operation of neuronal circuits, it is essential to dissect their cellular and synaptic components and to reveal their patterning of activity. One approach is to artificially enhance glutamatergic currents by using low Mg2+ media. There are several computational and technical reasons for this. First, by increasing time constant of excitation through NMDA receptors, the window for summation is enlarged, which facilitates the recruitment process (Herron et al., 1985; Jones and Heinemann, 1988). Second, it has been theoretically demonstrated that slow excitatory synapses favor sustained network activity (Wang, 1999). This would allow detection of formerly hidden functional circuits. Additionally, in the zero Mg2+ model (1) recurrent inhibition is not directly impaired (Ling and Benardo, 1995), and (2) NMDA receptors contribute a small tonic depolarizing current that inject firing into the network (LoTurco et al., 1990; Myme et al., 2003). We also tested the effect of increasing cell excitability under 5 mm K+ and the possible developmental changes.
We found that focal activity is sculpted by the dynamic interactions between some bursting cells leading other bursting cells, RS cells, and FS interneurons. Because the activation of FS inter-neurons is time-locked to the firing of leader cells, inhibition is always adjusted by excitation. This is consistent with previous studies in the subiculum using extracellular stimulation (Finch et al., 1988; Taube, 1993; Behr et al., 1998; Menendez de la Prida, 2003). Follower bursting and RS cells generate a secondary front of excitation that could further increase cell firing and compensate with additional inhibition. As a result, local ensembles of active cells develop. Ensemble activity is very variable, suggesting that it does not always involve the same cellular groups. Using cell-attached recordings, we showed that a given cell variably participates in different focal events and that the definition of leader/follower has to be based on statistical average.
A classical theoretical constraint for the generation of focal activity is short-range excitation and long-range inhibition (Hartline et al., 1956; Gutkin et al., 2001). No morphological data are yet available to draw the picture in the subiculum, although broad excitatory arbors have been described (Berger et al., 1980; Finch et al., 1983; Harris et al., 2001). Computer simulations have shown that, even with excitation wider than inhibition, focal activity could still develop (Xing and Gerstein, 1996; Kang et al., 2003). Under this configuration, however, this activity could easily become unstable and spread throughout (Golomb and Amitai, 1997).
This picture is in agreement with our results: focal activity could generate widespread population activity. This transition is mainly controlled by fast GABAergic transmission. During the first minutes under PTX, cross-correlation significantly increased and widespread events emerged more frequently. Transition from focal to widespread activity is also regulated by glutamatergic components because CNQX decreased the frequency of focal and widespread activity in parallel. Possible factors underlying this transition include runaway excitation driven by follower glutamatergic cells. Because the circuit will compensate for that increase with more inhibition, additional mechanisms are needed to explain this transition. For example, it is possible that plasticity of synaptic activity at moderate and high-frequency firing could differentially affect glutamatergic and GABAergic synapses (Thomson and Deuchars, 1994; Galarreta and Hestrin, 1998), or significantly shift reversal potential for GABAA-mediated chloride currents (Thompson and Gähwiler, 1989; Staley and Proctor, 1999).
Comparison with previous studies
Focal field activity in the subiculum differs from other types of correlated activity described in the neocortex and the hippocampus. First, subicular focal activity is local and non-oscillatory as revealed by cross-correlation analysis. This is in contrast to many types of correlated activity that typically show oscillatory and propagating patterns (Steriade et al., 1993; Sanchez-Vives and McCormick, 2000; Wu et al., 2001) or to sharp-wave activity (Chrobak and Buzsaki, 1996). Instead, subicular focal activity is similar to spatially organized cortical ensembles (Mao et al., 2001; Cossart et al., 2003), although on a shorter time scale. Second, during focal activity, bursting cells lead population activity, whereas FS interneurons constitute a follower group. Such a temporal correlation is preserved during the buildup of widespread epileptiform activity. This is an important computational difference in comparison with models of slow propagating activity in which inhibitory firing precedes excitatory firing (Golomb and Ermentrout, 2002; Compte et al., 2003). Finally, repetitions of a precise firing pattern are not a feature of subicular focal activity (Cossart et al., 2003). Instead, the configuration of focal ensembles greatly varies with cells being recruited dynamically.
The firing correlation between driver and follower cells defines a succession of excitation and inhibition. This sequence is common in many preparations and represents the fact that local inhibition is secondarily recruited by excitatory circuits (Buzsaki, 1984). In contrast to other situations, in which excitation and inhibition were balanced continuously (Gutkin et al., 2001; Shu et al., 2003), subicular focal activity exhibited an initial dominance of excitation and a late prevalence of inhibition. Our analysis of synaptic conductances associated with focal activity revealed that conductance peaks were related with a slightly depolarized potential that subsequently hyperpolarized. This strongly suggests that the largest conductance increases (100-400%) are dominated by GABA-mediated currents, which in concomitance with a net hyperpolarizing driving force would shunt excitatory inputs and suppress cell firing. This mechanism is similar to that reported in the visual (Borg-Graham et al., 1998; Hirsch et al., 1998; Anderson et al., 2000) and auditory cortices (Zhang et al., 2003; Wehr and Zador, 2003).
The role of bursting and RS cells
The leading role of bursting cells has been reported in different acute preparations (Miles and Wong, 1983; Changac-Amitai and Connors, 1989; Jensen and Yaari, 1997) and in chronic epileptic tissue (Sanabria et al., 2001; Cohen et al., 2002). Indeed, it has been shown that a permanent upregulation of intrinsic bursting contributes to hyperexcitability in chronic epilepsy (Wellmer et al., 2002). In contrast to single spikes, bursts may enhance the security of spike transmission, depending on the release probability and the facilitating nature of synapses (Miles and Wong, 1986; Thomson and Deuchars, 1994). Our results support a role of some bursting cells leading field activity, in agreement with previous reports in the subiculum (Harris and Stewart, 2001). However, we approached this problem using cell-attached recordings that do not alter the firing capability of cells. In addition, we examined the firing patterns in statistical terms and evaluated their impact in local firing control.
We found that bursting cells are not a leader population per se, although many bursting cells behaved like leaders (σ > 0.5). In contrast, RS cells mostly behaved as followers. However, focal and widespread activities constitute a probabilistic process, and a given cell could lead or follow the field during different episodes. Therefore, it is the cooperative interactions between different cell types that generate a population event. The positive correlation between the driving index and the firing threshold and the number of distal apical dendrites suggests that leader cells would be those that are more easily recruited. Also, cells having a statistical leader role were of the IB+ and IB-bursting type, which confirms the critical contribution of these firing patterns in recurrent activity (Miles and Wong, 1983). Therefore, bursting subicular cells having more apical dendrites better contribute to collect and distribute activity through the network. The distal dendrites of these neurons can be also the site of the initiation of regenerative calcium spikes, which favor cell recruitment (Yuste et al., 1994; Kamondi et al., 1998; Larkum et al., 1999).
The role of followers in this scheme should not be considered secondary. This is particularly the case of RS cells, which displayed a low firing rate and high firing threshold. This determines that RS cells were sharply recruited and, therefore, the contrast between response and no-response is enhanced at the onset of field activity. However, RS cells remained silent during some episodes of focal activity, which possibly leaves them available for being recruited by other ensembles. In contrast, because follower cells inject firing into the circuit at the onset of field activity, they can participate in runaway excitation.
Footnotes
This work was supported by grants from the Spanish Ministry of Science and Technology (Programa Ramón y Cajal and BFI2003-04305). B.G. is supported by Convenio UEM OTRI 200301CII from the Universidad Europea de Madrid. We thank Drs. L. Graham, A. Araque, S. Canals, J. Tabak, A. Compte, and O. Herreras for helpful discussions and critical reading of this manuscript. We also thank Dr. R. Miles for discussion, Dr. I. Cohen for Labview routines for spike detection, and J. M. Ibarz for technical help.
Correspondence should be addressed to Dr. L. Menendez de la Prida, Departamento de Neurobiología-Investigación (planta-1D), Hospital Ramón y Cajal, Ctra Colmenar Km 9, Madrid 28034, Spain. E-mail: liset.m.prida{at}hrc.es.
Copyright © 2004 Society for Neuroscience 0270-6474/04/245525-12$15.00/0