Tetrodes markedly improve the reliability and yield of multiple single-unit isolation from multi-unit recordings in cat striate cortex

Spike sorting based on discrete wavelet transform coefficients

Using the novel mathematical technique known as wavelet analysis, a new method (WSC) is presented to sort spikes according to a decomposition of neural signals in the time-frequency space. The WSC method is implemented by a pyramidal algorithm that acts upon neural signals as a bank of quadrature mirror filters. This algorithm is clearly explained and an overview of the mathematical background of wavelet analysis is given. An artificial spike train, especially designed to test the specificity and sensibility of sorting procedures, was used to assess the performance of the WSC method as well as of methods based on principal component analysis (PCA) and reduced feature set (RFS). The WSC method outperformed the other two methods. Its superior performance was largely due to the fact that spike profiles that could not be separated by previous methods (because of the similarity of their temporal profile and the masking action of noise) were separable by the WSC method. The WSC method is particularly noise resistant, as it implicitly eliminates the irrelevant information contained in the noise frequency range. But the main advantage of the WSC method is its use of parameters that describe the joint time-frequency localization of spike features to build a fast and unspecialized pattern recognition procedure.

Relationship among discharges of neighboring neurons in the rat prefrontal cortex during spatial working memory tasks

The relationship among discharges of neurons that were recorded simultaneously with tetrodes in the rat medial prefrontal cortex was analyzed. Spatial working memory tasks were divided into several distinct stages based on the behavioral correlates of individual neurons, and interneuronal correlation of signal (mean discharge rate at each stage) and noise (trial-to-trial deviation from the signal) was calculated. Behavioral correlates of neighboring neurons were quite heterogeneous and, accordingly, average signal correlation was relatively low ( approximately 0.16). Noise correlation was even lower ( approximately 0.06), but neuronal noise was more correlated among the neurons with similar signals. Spikes underlying the signal and noise correlation among the prefrontal cortical neurons were loosely synchronized over a few hundred milliseconds. These results suggest that neighboring prefrontal cortical neurons process largely independent information and have weakly correlated noise and that precisely synchronized spikes play a relatively minor role in producing the correlated signal and noise among these neurons.

Cross-correlation measures of unresolved multi-neuron recordings

An increasing number of laboratories are studying population properties of the nervous system using data where the spike activity of more than one neuron is recorded on each electrode and where, accidentally or deliberately, these activities are not resolved into single unit spike trains. We have previously examined the consequences for measurement of cross-correlation between two such electrodes in the limited case where all individual distant (between electrode) correlations are the same and all individual close (on a single electrode) correlations are the same [Bedenbaugh, P.H., and Gerstein, G.L. (1997). Multiunit normalized cross correlation differs from the average single-unit normalized correlation. Neural Computation 9, 1265-1275]. Here, we lift these unrealistic restrictions to allow all values of individual correlation, and examine explicitly the cases of two or three unresolved neurons on each electrode. In these situations, the cross-correlation coefficient measured between the electrodes is a linear sum of the distant correlations, divided by a non-linear function of the close correlations. We then examine in detail the case of a single direct distant correlation and take account of all relevant indirect correlations. The measured interelectrode correlation shows a reduction of this actual distant correlation by a non-linear function of the close correlations on each electrode over most of their possible values. Finally, we examine the consequences of poor waveform sorting for correlation measures; here a supposedly isolated spike train is contaminated by some fraction of the activity of another train, a situation that unfortunately is all too common in experiments. All these distortions become far more serious in the more realistic situation of dynamic firing rates and correlations. This paper is intended as a cautionary note for those who want to draw inferences about neuronal organization and/or coding or representation by using cross-correlation analysis of unresolved recordings.

Intracellular features predicted by extracellular recordings in the hippocampus in vivo

Multichannel tetrode array recording in awake behaving animals provides a powerful method to record the activity of large numbers of neurons. The power of this method could be extended if further information concerning the intracellular state of the neurons could be extracted from the extracellularly recorded signals. Toward this end, we have simultaneously recorded intracellular and extracellular signals from hippocampal CA1 pyramidal cells and interneurons in the anesthetized rat. We found that several intracellular parameters can be deduced from extracellular spike waveforms. The width of the intracellular action potential is defined precisely by distinct points on the extracellular spike. Amplitude changes of the intracellular action potential are reflected by changes in the amplitude of the initial negative phase of the extracellular spike, and these amplitude changes are dependent on the state of the network. In addition, intracellular recordings from dendrites with simultaneous extracellular recordings from the soma indicate that, on average, action potentials are initiated in the perisomatic region and propagate to the dendrites at 1.68 m/s. Finally we determined that a tetrode in hippocampal area CA1 theoretically should be able to record electrical signals from approximately 1, 000 neurons. Of these, 60-100 neurons should generate spikes of sufficient amplitude to be detectable from the noise and to allow for their separation using current spatial clustering methods. This theoretical maximum is in contrast to the approximately six units that are usually detected per tetrode. From this, we conclude that a large percentage of hippocampal CA1 pyramidal cells are silent in any given behavioral condition.

Accuracy of tetrode spike separation as determined by simultaneous intracellular and extracellular measurements

Simultaneous recording from large numbers of neurons is a prerequisite for understanding their cooperative behavior. Various recording techniques and spike separation methods are being used toward this goal. However, the error rates involved in spike separation have not yet been quantified. We studied the separation reliability of "tetrode" (4-wire electrode)-recorded spikes by monitoring simultaneously from the same cell intracellularly with a glass pipette and extracellularly with a tetrode. With manual spike sorting, we found a trade-off between Type I and Type II errors, with errors typically ranging from 0 to 30% depending on the amplitude and firing pattern of the cell, the similarity of the waveshapes of neighboring neurons, and the experience of the operator. Performance using only a single wire was markedly lower, indicating the advantages of multiple-site monitoring techniques over single-wire recordings. For tetrode recordings, error rates were increased by burst activity and during periods of cellular synchrony. The lowest possible separation error rates were estimated by a search for the best ellipsoidal cluster shape. Human operator performance was significantly below the estimated optimum. Investigation of error distributions indicated that suboptimal performance was caused by inability of the operators to mark cluster boundaries accurately in a high-dimensional feature space. We therefore hypothesized that automatic spike-sorting algorithms have the potential to significantly lower error rates. Implementation of a semi-automatic classification system confirms this suggestion, reducing errors close to the estimated optimum, in the range 0-8%.

Multiple site silicon-based probes for chronic recordings in freely moving rats: implantation, recording and histological verification

This paper describes the procedure of assembling a miniature microdrive and silicon probe system for surgical implantation into the adult rat brain. Successful recordings of single and multiunit activity with parallel depth profiles of spontaneous and evoked field potentials are shown. The procedure for histological verification of the position of the silicon probe is described.

Fast oscillations display sharper orientation tuning than slower components of the same recordings in striate cortex of the awake monkey

We wanted to know whether fast oscillations ( approximately 30-80 Hz) in striate cortex of awake monkeys show sharper orientation selectivity than (i) slower components, including spike rate modulations, and (ii) broad-band signals of the same recordings. As fast oscillations are probably of cortical origin this may further clarify whether cortical network mechanisms are substantially involved in generating orientation selectivity. We recorded multi unit activity (MUA) and local field potentials (LFP, 1-140 Hz) by the same microelectrodes from upper layers of macaque striate cortex during visual stimulation with grating textures of different orientations. An orientation index (OI) was derived from the cortical responses in three frequency ranges (low, 0-11.7 Hz; medium, 11.7-31.3 Hz; and fast oscillations, 31.3-62.5 Hz) and for the broad-band LFP and MUA power. (i) Both LFP and MUA fast oscillations reveal a higher orientation index than signal components in the low and medium frequency ranges. (ii) For MUA the orientation index was significantly higher with fast oscillations than for the lower frequency ranges and the initial broad-band transient responses. (iii) LFPs show a significantly higher orientation index only for the fast oscillations during sustained activation compared with their broad-band power during the transient responses. Thus, our main result is the sharper orientation tuning of fast oscillations in spike activities of local populations compared with slower components of the same broad-band recordings. As fast oscillations occur synchronized in the awake monkey's striate cortex we assume that they have enhanced probability of activating successive stages of visual processing and hence contribute to the perception of orientation.

Functional coupling shows stronger stimulus dependency for fast oscillations than for low-frequency components in striate cortex of awake monkey

It has been argued that coupling among the neural signals activated by a visual object supports binding of local features into a coherent object perception. During visual stimulation by a grating texture we studied functional coupling by calculating spectral coherence among pairs of signals recorded in the striate cortex of awake monkeys. Multiple unit activity (MUA) and local field potentials (LFP, 1-140 Hz) were extracted from seven parallel broad band recordings. Spectral coherence was dominated by high-frequency oscillations in the range 35-50 Hz and often by additional low-frequency components (0-12 Hz). Functional coupling among separate cortical sites was more stimulus specific for MUA than for LFP: MUA coherence at high and low frequencies depended highly significantly on: (i) the similarity of the preferred orientations at the two sites - the more similar the higher the coherence; (ii) the orientation of the stimulus grating - with highest coherence at half angle between the preferred orientations at the two sites; (iii) cortical distance - coherence decreases to noise levels at approximately 3 mm (MUA) and 6 mm (LFP). Coherence of fast oscillations did not depend on the degree of coaxiality of the orientation-sensitive receptive fields, whereas low frequencies showed significant dependency. This indicates that different frequency components can engage different coupling networks in the striate cortex which probably support different coding tasks. Changes in average oscillation frequency with stimulus orientation were highly significant for fast oscillations while there was no dependency for low frequencies. Finally, stimulus-related spectral power and coherence of fast oscillations were considerably higher than of low frequency components. Fast oscillations may therefore contribute more to feature binding and coding of object continuity than low-frequency components, at least for texture surfaces as analysed here.

Trial-to-trial variability and state-dependent modulation of auditory-evoked responses in cortex

Recent experimental work has provided evidence that trial-to-trial variability of sensory-evoked responses in cortex can be explained as a linear superposition of random ongoing background activity and a stationary response. While studying single trial variability and state-dependent modulation of evoked responses in auditory cortex of ketamine/xylazine-anesthetized rats, we have observed an apparent violation of this model. Local field potential and unit spike trains were recorded and analyzed during different anesthesia depths-deep, medium, and light-which were defined by the pattern of ongoing cortical activity. Estimation of single trial evoked response was achieved by considering whole waveforms, rather than just one or two peak values from each wave. Principal components analysis was used to quantitatively classify waveforms on the basis of their time courses (i.e., shapes). We found that not only average response but also response variability is modulated by depth of anesthesia. Trial-to-trial variability is highest under medium levels of anesthesia, during which ongoing cortical activity exhibits rhythmic population bursting activity. By triggering the occurrence of stimuli from the spontaneously occurring burst events, we show that the observed variability can be accounted for by the background activity. In particular, the ongoing activity was found to modulate both amplitude and shape (including latency) of evoked local field potentials and evoked unit activity in a manner not predicted by linear superposition of background activity and a stereotyped evoked response. This breakdown of the linear model is likely attributable to rapid transitions between different levels of thalamocortical excitability (e.g., spike-wave discharges), although brain "state" is relatively fixed.

Responses of macaque perirhinal neurons during and after visual stimulus association learning

Recent lesion studies have implicated the perirhinal cortex in learning that two objects are associated, i.e., visual association learning. In this experiment we tested whether neuronal responses to associated stimuli in perirhinal cortex are altered over the course of learning. Neurons were recorded from monkeys during performance of a visual discrimination task in which a predictor stimulus was followed, after a delay, by a GO or NO-GO choice stimulus. Association learning had two major influences on neuronal responses. First, responses to frequently paired predictor-choice stimuli were more similar to one another than was the case with infrequently paired stimuli. Second, the magnitude of activity during the delay was correlated with the magnitude of responses to both the predictor and choice stimuli. Both of these learning effects were found only for stimulus pairs that had been associated on at least 2 d of training. Early in training, the delay activity was correlated only with the response to the predictor stimuli. Thus, with long-term training, perirhinal neurons tend to link the representations of temporally associated stimuli.

Interaction between spike waveform classification and temporal sequence detection

In vivo extracellular recordings have allowed researchers to study the response properties of neurons to behaviorally relevant stimuli. In this paper we use multiple tetrode recordings from the hippocampus of the freely behaving rat to show that the action potential amplitude of a given cell can vary in a systematic and activity dependent manner over behaviorally relevant time scales. Since the discrimination algorithms used by experimenters to isolate cells from extracellular recordings are based on differences in waveforms, we show how these systematic changes in waveform shape can lead to non-random errors in single cell isolation. We further demonstrate that these non-random errors can lead to apparent temporal ordering effects between neurons in the absence of any specific temporal relationship. A firm understanding of these naturally occurring physiological changes is therefore critical for the evaluation of higher order phenomena such as the temporally correlated firing of ensembles of neurons.

Coding of locomotor phase in populations of neurons in rostral and caudal segments of the neonatal rat lumbar spinal cord

Several experiments have demonstrated that rostral segments of the vertebrate lumbar spinal cord produce a rhythmic motor output more readily and of better quality than caudal segments. Here we examine how this rostrocaudal gradient of rhythmogenic capability is reflected in the spike activity of neurons in the rostral (L(2)) and caudal (L(5)) lumbar spinal cord of the neonatal rat. The spike activity of interneurons in the ventromedial cord, a region necessary for the production of locomotion, was recorded intracellularly with patch electrodes and extracellularly with tetrodes during pharmacologically induced locomotion. Both L(2) and L(5) neurons tended to be active in phase with their homologous ventral root. L(5) neurons, however, had a wider distribution of their preferred phases of activity throughout the locomotor cycle than L(2) neurons. The strength of modulation of the activity of individual L(2) neurons was also larger than that of L(5) neurons. These differences resulted in a stronger rhythmic signal from the L(2) neuronal population than from the L(5) population. These results demonstrate that the rhythmogenic capability of each spinal segment was reflected in the activity of interneurons located in the same segment. In addition to paralleling the rostrocaudal gradient of rhythmogenic capability, these results further suggest a colocalization of motoneurons and their associated interneurons involved in the production of locomotion.

Long-term neural recording characteristics of wire microelectrode arrays implanted in cerebral cortex

This paper describes a detailed protocol for obtaining chronic, multi-site unit recordings in cerebral cortex of awake animals for periods of three months or more. The protocol includes details for making relatively simple and inexpensive implantable multichannel electrodes that consist of arrays of separate microwires. The results reported in this paper suggest that a viable implant will have discriminable unit activity on about 80% of the electrodes, resulting in, on average, the simultaneous unit recording of upwards of 60 units during a daily recording session. The active electrodes during one recording session tend to remain active in subsequent recording sessions for several weeks. Using the methods described here, implants have been constructed which incorporate several different electrode materials, coatings, sizes, and electrode separation within a single array. These microwire electrode arrays provide the basic technology for obtaining unit recordings for several months. This provides a model system for studying biocompatibility of neural implants, which is a critical component for the development of neural implants that have an indefinite working span.

Building neural representations of habits

Memories for habits and skills ("implicit or procedural memory") and memories for facts ("explicit or episodic memory") are built up in different brain systems and are vulnerable to different neurodegenerative disorders in humans. So that the striatum-based mechanisms underlying habit formation could be studied, chronic recordings from ensembles of striatal neurons were made with multiple tetrodes as rats learned a T-maze procedural task. Large and widely distributed changes in the neuronal activity patterns occurred in the sensorimotor striatum during behavioral acquisition, culminating in task-related activity emphasizing the beginning and end of the automatized procedure. The new ensemble patterns remained stable during weeks of subsequent performance of the same task. These results suggest that the encoding of action in the sensorimotor striatum undergoes dynamic reorganization as habit learning proceeds.

Hippocampal encoding of non-spatial trace conditioning

Trace eyeblink classical conditioning is a non-spatial learning paradigm that requires an intact hippocampus. This task is hippocampus-dependent because the auditory tone conditioned stimulus (CS) is temporally separated from the corneal airpuff unconditioned stimulus (US) by a 500-ms trace interval. Our laboratory has performed a series of neurophysiological experiments that have examined the activity of pyramidal cells in the CA1 area of the hippocampus during trace eyeblink conditioning. We have found that the non-spatial stimuli involved in this paradigm are encoded in the hippocampus in a logical order that is necessary for their association and the subsequent expression of behavioral learning. Although there were many profiles of single neurons responding to the CS-US trial during training, the majority of the neurons showed an increase in activity to the airpuff-US. Prior to learning, it appears that hippocampal cells and ensembles of cells were preferentially attending to the stimulus with immediate behavioral importance, the US. Hippocampal cells then began to respond to the associated neutral stimulus, the CS. Shortly thereafter, animals began to show increases in the behavioral expression of CRs. In some experiments, hippocampal neurons from aged animals exhibited impairments in the encoding of CS and US information. These aged animals were not able to associate these stimuli and acquire trace eyeblink CRs. Our findings along with the findings of other spatial learning studies, suggest that the hippocampus is involved in encoding information about discontiguous sets of stimuli, either spatial or nonspatial, especially early in the learning process.

Examination of the spatial and temporal distribution of sensory cortical activity using a 100-electrode array

This paper introduces improved techniques for multichannel extracellular electrophysiological recordings of neurons distributed across a single layer of topographically mapped cortex. We describe the electrode array, the surgical implant techniques, and the procedures for data collection and analysis. Neural events are acquired through an array of 25 or 100 microelectrodes with a 400-microm inter-electrode spacing. One advantage of the new methodology is that implantation is achieved through transdural penetration, thereby reducing the disruption of the cortical tissue. The overall cortical territory sampled by the 25-electrode array is 1.6 x 1.6 mm (2.56 mm2) and by the 100-electrode array 3.6 x 3.6 mm (12.96 mm2). Using a recording system with 100 channels available, neural activity is simultaneously acquired on all electrodes, amplified, digitized, and stored on computer. In our data, average peak-to-peak signal/noise ratio was 11.5 and off-line waveform analysis typically allowed the separation of at least one well-discriminated single-unit per channel. The reported technique permits analysis of cortical function with high temporal and spatial resolution. We use the technique to create an 'image' of neural activity distributed across the whisker representation of rat somatosensory (barrel) cortex.

Stimulation of the paraventricular nucleus modulates firing of neurons in the nucleus of the solitary tract

The present study assessed whether the baroreflex inhibition elicited by electrical stimulation of the hypothalamic paraventricular nucleus (PVN) involves altered activity in the nucleus of the solitary tract (NTS). Unit recordings were made from 107 neurons in the NTS in anesthetized rabbits. Intravenous phenylephrine was used to induce a pressor response and to activate baroreflexes. Of the neurons that responded to pressor responses, two-thirds were excited and one-third was inhibited. Stimulation of the PVN inhibited 70% of the phenylephrine-responsive NTS neurons, with or without concurrent baroreceptor stimulation. When PVN stimulation was delivered concurrently with phenylephrine injection, more NTS neuronal inhibition and less excitation occurred than with phenylephrine alone. Usually PVN stimulation inhibited NTS neurons that were excited by pressor responses; less commonly, PVN stimulation excited NTS neurons that were inhibited by pressor responses. The findings are consistent with the view that PVN activation during the defense reaction inhibits baroreflexes by altering firing of NTS neurons.

Neural mechanisms for encoding binocular disparity: receptive field position versus phase

The visual system uses binocular disparity to discriminate the relative depth of objects in space. Because the striate cortex is the first site along the central visual pathways at which signals from the left and right eyes converge onto a single neuron, encoding of binocular disparity is thought to begin in this region. There are two possible mechanisms for encoding binocular disparity through simple cells in the striate cortex: a difference in receptive field (RF) position between the two eyes (RF position disparity) and a difference in RF profiles between the two eyes (RF phase disparity). Although there is evidence that supports each of these schemes, both mechanisms have not been examined in a single study to determine their relative roles. In this study, we have measured RF position and phase disparities of individual simple cells in the cat's striate cortex to address this issue. Using a sophisticated RF mapping technique that employs binary m-sequences, we have obtained left and right eye RF profiles of two or more cells recorded simultaneously. A version of the reference-cell method was used to estimate RF position disparity. We find that RF position disparities generally are limited to values that are not sufficient to encode large binocular disparities. In contrast, RF phase disparities cover a wide range of binocular disparities and exhibit dependencies on RF orientation and spatial frequency in a manner expected for a mechanism that encodes binocular disparity. These results suggest that binocular disparity is encoded mainly through RF phase disparity. However, RF position disparity may play a significant role for cells with high spatial frequency selectivity that are constrained to have only small RF phase disparities.

Ultra-miniature headstage with 6-channel drive and vacuum-assisted micro-wire implantation for chronic recording from the neocortex

We describe a head-stage, with precision microtranslators for the chronic placement of micro-wire electrodes in the neocortex, that minimizes compressive damage to the brain. The head-stage has a diameter of 5.8 mm and allows six electrodes, separated by 450 microm on a hexagonal grid, to be individually and continuously positioned throughout a depth of approximately 3 mm. Suction is used to transiently support the dura against a curved array of tubes that guide and stabilize the electrodes as a means to prevent compression of the neocortex as the electrodes breach the dura. With this headstage we recorded extracellular signals in a rat immediately after surgery. Single-unit waveforms at a given electrode position were stable for at least several hours in the freely behaving animal and were obtained throughout the depth of the neocortex for at least 2 months. Electrophysiological records and histological examination showed that the upper layers of the neocortex were intact and minimally damaged after the implantation.

Cellular mechanisms contributing to response variability of cortical neurons in vivo

Cortical neurons recorded in vivo exhibit highly variable responses to the repeated presentation of the same stimulus. To further understand the cellular mechanisms underlying this phenomenon, we performed intracellular recordings from neurons in cat striate cortex in vivo and examined the relationships between spontaneous activity and visually evoked responses. Activity was assessed on a trial-by-trial basis by measuring the membrane potential (Vm) fluctuations and spike activity during brief epochs immediately before and after the onset of an evoked response. We found that the response magnitude, expressed as a change in Vm relative to baseline, was linearly correlated with the preceding spontaneous Vm. This correlation was enhanced when the cells were hyperpolarized to reduce the activation of voltage-gated conductances. The output of the cells, expressed as spike counts and latencies, was only moderately correlated with fluctuations in the preceding spontaneous Vm. Spike-triggered averaging of Vm revealed that visually evoked action potentials arise from transient depolarizations having a rise time of approximately 10 msec. Consistent with this, evoked spike count was found to be linearly correlated with the magnitude of Vm fluctuations in the gamma (20-70 Hz) frequency band. We also found that the threshold of visually evoked action potentials varied over a range of approximately 10 mV. Examination of simultaneously recorded intracellular and extracellular activity revealed a correlation between Vm depolarization and spike discharges in adjacent cells. Together these results demonstrate that response variability is attributable largely to coherent fluctuations in cortical activity preceding the onset of a stimulus, but also to variations in action potential threshold and the magnitude of high-frequency fluctuations evoked by the stimulus.

Dynamics of tremor-related oscillations in the human globus pallidus: a single case study

Physiological evidence indicates that the resting tremor of Parkinson's disease originates in oscillatory neural activity in the forebrain, but it is unknown whether that activity is globally synchronized or consists of parallel, independently oscillating circuits. In the present study, we used dual microelectrodes to record tremor-related neuronal activity from eight sites in the internal segment of the globus pallidus (GPi) from an awake Parkinson's disease patient undergoing stereotaxic pallidotomy. We utilized spectral analysis to evaluate the temporal correlations between multiunit activity at spatially separated sites and between neural and limb electromyographic activity. We observed that some GPi neural pairs oscillated synchronously at the tremor frequency, whereas other neural pairs oscillated independently. Additionally, we found that GPi tremor-related activity at a given site could fluctuate between states of synchronization and independence with respect to upper limb tremor. Consistent with this finding, some paired recording sites within GPi showed periods of transient synchronization. These observations support the hypothesis of independent tremor-generating circuits whose coupling can fluctuate over time.

Oscillatory coupling of hippocampal pyramidal cells and interneurons in the behaving Rat

We examined whether excitation and inhibition are balanced in hippocampal cortical networks. Extracellular field and single-unit activity were recorded by multiple tetrodes and multisite silicon probes to reveal the timing of the activity of hippocampal CA1 pyramidal cells and classes of interneurons during theta waves and sharp wave burst (SPW)-associated field ripples. The somatic and dendritic inhibition of pyramidal cells was deduced from the activity of interneurons in the pyramidal layer [int(p)] and in the alveus and st. oriens [int(a/o)], respectively. Int(p) and int(a/o) discharged an average of 60 and 20 degrees before the population discharge of pyramidal cells during the theta cycle, respectively. SPW ripples were associated with a 2.5-fold net increase of excitation. The discharge frequency of int(a/o) increased, decreased ("anti-SPW" cells), or did not change ("SPW-independent" cells) during SPW, suggesting that not all interneurons are innervated by pyramidal cells. Int(p) either fired together with (unimodal cells) or both before and after (bimodal cells) the pyramidal cell burst. During fast-ripple oscillation, the activity of interneurons in both the int(p) and int(a/o) groups lagged the maximum discharge probability of pyramidal neurons by 1-2 msec. Network state changes, as reflected by field activity, covaried with changes in the spike train dynamics of single cells and their interactions. Summed activity of parallel-recorded interneurons, but not of pyramidal cells, reliably predicted theta cycles, whereas the reverse was true for the ripple cycles of SPWs. We suggest that network-driven excitability changes provide temporal windows of opportunity for single pyramidal cells to suppress, enable, or facilitate selective synaptic inputs.

Sustained activation of hippocampal pyramidal cells by 'space clamping' in a running wheel

In contrast to sensory cortical areas of the brain, the relevant physiological inputs to the hippocampus, leading to selective activation of pyramidal cells, are largely unknown. Pyramidal cells are thought to be phasically activated by spatial cues and a variety of sensory and motor stimuli. Here, we used a behavioural 'space clamp' method, which involved the confinement of the actively running animal in a defined position in space (running wheel) and kept sensory inputs constant. Twelve percent of the recorded CA1 pyramidal cells were selectively active while the rat was running in the wheel. Cell firing was specific to the direction of running and disappeared after rotating the recording apparatus. The discharge frequency of pyramidal cells and interneurons was sustained as long as the rat ran continuously in the wheel. Furthermore, the discharge frequency of pyramidal cells and interneurons increased with increasing running velocity, even though the frequency of hippocampal theta waves remained constant. The discharge frequency of some 'wheel-related' pyramidal cells could increase more than 10-fold between 10 and 100 cm/s, whereas the firing rate of 'non-wheel' cells remained constantly low. We hypothesize that: (i) a necessary condition for place-specific discharge of hippocampal pyramidal cells is the presence of theta oscillation; and (ii) relevant stimuli can tonically and selectively activate hippocampal pyramidal cells as long as theta activity is present.

Azimuth coding in primary auditory cortex of the cat. I. Spike synchrony versus spike count representations

The neural representation of sound azimuth in auditory cortex most often is considered to be average firing rate, and azimuth tuning curves based thereupon appear to be rather broad. Coincident firings of simultaneously recorded neurons could provide an improved representation of sound azimuth compared with that contained in the firing rate in either of the units. In the present study, a comparison was made between local field potentials and several measures based on unit firing rate and coincident firing with respect to their azimuth-tuning curve bandwidth. Noise bursts, covering a 60-dB intensity range, were presented from nine speakers arranged in a semicircular array with a radius of 55 cm in the animal's frontal half field. At threshold intensities, all local field potential (LFP) recordings showed preferences for contralateral azimuths. Multiunit recordings showed in 74% a threshold for contralateral azimuths, in 16% for frontal azimuths, and in only 5% showed an ipsilateral threshold. The remaining 5% were not spatially tuned. Representations for directionally sensitive units based on coincident firings provided significantly sharper tuning (50-60 degrees bandwidth at 25 dB above the lowest threshold) than those based on firing rate (bandwidths of 80-90 degrees). The ability to predict sound azimuth from the directional information contained in the neural population activity was simulated by combining the responses of the 102 single units. Peak firing rates and coincident firings with LFPs at the preferred azimuth for each unit were used to construct a population vector. At stimulus levels of >/=40 dB SPL, the prediction function was sigmoidal with the predicted frontal azimuth coinciding with the frontal speaker position. Sound azimuths >45 degrees from the midline all resulted in predicted values of -90 or 90 degrees, respectively. No differences were observed in the performance of the prediction based on firing rate or coincident firings for these intensities. This suggests that although coincident firings produce narrower azimuth tuning curves, the information contained in the overall neural population does not increase compared with that contained in a firing rate representation. The relatively poor performance of the population vector further suggests that primary auditory cortex does not code sound azimuth by a globally distributed measure of peak firing rate or coincident firing.

Analysing functional connectivity in brain slices by a combination of infrared video microscopy, flash photolysis of caged compounds and scanning methods

We evaluate a novel set-up for scanning functional connectivity in brain slices from the somatosensory cortex of the rat. Upright infrared video microscopy for targeted placement of electrodes is combined with rapid photolysis of bath-applied caged neurotransmitter induced by a xenon flash lamp. Flash photolysis of caged glutamate and electrical stimulation produce comparable field potential responses and demonstrate that the viability of the submerged slices exceeds several hours. Glutamate release leads to field potential responses whose two phases are differentially affected by selective blockade of N-methyl-D-aspartate- and alpha-amino-3-hydroxy-5-methyl-4-isoxazolepropionate-type glutamate receptors with DL-2-amino-5-phosphonovaleric acid and 1,2,3,4-tetrahydro-6-nitro-2,3-dioxo-benzo[f]quinoxaline-7-sulphonamide, respectively. Rapid computer-controlled scanning of hundreds of distinct stimulation sites with simultaneous recordings at a fixed reference site allows construction of functional input maps from peak amplitudes and delays to peak of field potential responses. Selective laminar expansion of the functional input maps after bicuculline application demonstrates that the combination of this conveniently assembled set-up with pharmacological and physical manipulations can provide insights into the determinants of functional connectivity in brain slices.

Reliability and state dependence of pyramidal cell-interneuron synapses in the hippocampus: an ensemble approach in the behaving rat

Spike transmission probability between pyramidal cells and interneurons in the CA1 pyramidal layer was investigated in the behaving rat by the simultaneous recording of neuronal ensembles. Population synchrony was strongest during sharp wave (SPW) bursts. However, the increase was three times larger for pyramidal cells than for interneurons. The contribution of single pyramidal cells to the discharge of interneurons was often large (up to 0.6 probability), as assessed by the presence of significant (<3 ms) peaks in the cross-correlogram. Complex-spike bursts were more effective than single spikes. Single cell contribution was higher between SPW bursts than during SPWs or theta activity. Hence, single pyramidal cells can reliably discharge interneurons, and the probability of spike transmission is behavior dependent.

The 'Ideal Homunculus': decoding neural population signals

Information processing in the nervous system involves the activity of large populations of neurons. It is possible, however, to interpret the activity of relatively small numbers of cells in terms of meaningful aspects of the environment. 'Bayesian inference' provides a systematic and effective method of combining information from multiple cells to accomplish this. It is not a model of a neural mechanism (neither are alternative methods, such as the population vector approach) but a tool for analysing neural signals. It does not require difficult assumptions about the nature of the dimensions underlying cell selectivity, about the distribution and tuning of cell responses or about the way in which information is transmitted and processed. It can be applied to any parameter of neural activity (for example, firing rate or temporal pattern). In this review, we demonstrate the power of Bayesian analysis using examples of visual responses of neurons in primary visual and temporal cortices. We show that interaction between correlation in mean responses to different stimuli (signal) and correlation in response variability within stimuli (noise) can lead to marked improvement of stimulus discrimination using population responses.

Identification of synaptic connections in neural ensembles by graphical models

A method for the identification of direct synaptic connections in a larger neural net is presented. It is based on a conditional correlation graph for multivariate point processes. The connections are identified via the partial spectral coherence of two neurons, given all others. It is shown how these coherences can be calculated by inversion of the spectral density matrix. In simulations with GENESIS, we discuss the relevance of the method for identifying different neural ensembles including an excitatory feedback loop and networks with lateral inhibitions.

Coding for auditory space in the nucleus of the brachium of the inferior colliculus in the ferret

Coding for auditory space in the nucleus of the brachium of the inferior colliculus in the ferret. J. Neurophysiol. 78: 2717-2731, 1997. The nucleus of the brachium of the inferior colliculus (BIN) projects topographically to the deeper layers of the superior colliculus (SC), which contain a two-dimensional map of auditory space. In this study, we have used broadband stimuli presented in the free field to investigate how auditory space is represented in the BIN of the ferret. Response latencies and temporal firing patterns were comparable with those in the SC, and both properties showed some variation with stimulus location. We obtained spatial response profiles at two sound levels (5-15 and 25-35 dB above unit threshold). A large proportion of azimuth profiles (41% in the suprathreshold condition, 80% in the near-threshold condition) presented a single peak, indicating that they were tuned to single regions in space. For some of these units, the preferred speaker position varied considerably with sound level. The remaining units showed predominantly either broad "hemifield" or spatially ambiguous "bilobed" response profiles. At suprathreshold sound levels, the preferred azimuths of the tuned cells were ordered topographically along the rostrocaudal axis of the BIN, although this representation is considerably more scattered than that in the SC. In contrast to the SC, we observed no systematic variation in the distribution of near-threshold best azimuths, which were instead concentrated around the interaural axis in the contralateral hemifield. The azimuth tuning of individual units in the BIN was generally broader at both sound levels than that in the SC. Many units also were tuned for the elevation of the sound source (48% for supra-, 77% for near-threshold stimulation), but there was no evidence for topographic order in the distribution of preferred elevations within the BIN. These results suggest that the BIN sends inputs to the SC that are already selective for sound azimuth and elevation and that show some degree of topographic order for sound azimuth. These inputs then presumably are sharpened and their topography refined by a mechanism that is likely to involve convergence of other inputs and activity-dependent fine tuning of terminal connections, to result in a precise two-dimensional map of auditory space in the SC.

Memory reprocessing in corticocortical and hippocampocortical neuronal ensembles

Hippocampal cells that fire together during behaviour exhibit enhanced activity correlations during subsequent sleep, with some preservation of temporal order information. Thus, information reflecting experiences during behaviour is re-expressed in hippocampal circuits during subsequent 'offline' periods, as postulated by some theories of memory consolidation. If the hippocampus orchestrates the reinstatement of experience-specific activity patterns in the neocortex, as also postulated by such theories, then correlation patterns both within the neocortex and between hippocampus and neocortex should also re-emerge during sleep. Ensemble recordings were made in the posterior parietal neocortex, in CA1, and simultaneously in both areas, in seven rats. Each session involved an initial sleep episode (S1), behaviour on a simple maze (M), and subsequent sleep (S2). The ensemble activity-correlation structure within and between areas during S2 resembled that of M more closely than did the correlation pattern of S1. Temporal order (i.e. the asymmetry of the cross-correlogram) was also preserved within, but not between, structures. Thus, traces of recent experience are re-expressed in both hippocampal and neocortical circuits during sleep, and the representations in the two areas tend to correspond to the same experience. The poorer preservation of temporal firing biases between neurons in the different regions may reflect the less direct synaptic coupling between regions than within them. Alternatively, it could result from a shift, between behavioural states, in the relative dominance relations in the corticohippocampal dialogue. Between-structure order will be disrupted, for example, if, during behaviour, neocortical patterns tend to drive corresponding hippocampal patterns, whereas during sleep the reverse occurs. This possibility remains to be investigated.

The effect of aging on experience-dependent plasticity of hippocampal place cells

The firing characteristics of 1437 CA1 pyramidal neurons were studied in relation to both spatial location and the phase of the theta rhythm in healthy young and old rats performing a simple spatial task on a rectangular track. The old rats had previously been found to be deficient on the Morris spatial learning task. Age effects on the theta rhythm per se were minimal. Theta amplitude and frequency during rapid eye movement sleep were virtually identical. During behavior, theta frequency was slightly reduced with age. In both groups, cell firing occurred at progressively earlier phases of the theta rhythm as the rat traversed the place field of the cell (i. e., there was "phase precession," as reported by others). The net phase shift did not differ between age groups. The main finding of the study was a loss of experience-dependent plasticity in the place fields of old rats. During the first lap around the track on each day, the initial sizes of the place fields were the same between ages; however, place fields of young rats, but not old, expanded significantly during the first few laps around the track in a given recording session. As the place fields expanded, the rate of change of firing with phase slowed accordingly, so that the net phase change remained constant. Thus changes in field size and phase precession are coupled. A deficit in plasticity of place fields in old rats may lead to a less accurate population code for spatial location.

Sequence of single neuron changes in CA1 hippocampus of rabbits during acquisition of trace eyeblink conditioned responses

The sequence of changes in single neuron activity in the CA1 area of the rabbit hippocampus was examined during daily sessions (80 trials/session) of hippocampally dependent nonspatial trace eyeblink (i.e., nictitating membrane response) conditioning. Each trial for trace conditioned animals (n = 7) consisted of a tone conditioned stimulus (CS; 6 kHz; 90 dB, 100 ms) followed by a 500-ms silent trace period, then a corneal airpuff unconditioned stimulus (US; 3.0 psi; 150 ms). Control animals (n = 5) received unpaired CSs and USs. Most pyramidal (n = 309) and theta (n = 21) cells were recorded for a single day of training. The activity of cells for each animal were grouped according to: the day of training that CRs began to increase and the day of training that CR performance became asymptotic. Pyramidal cells from trace conditioned animals demonstrated several stages of learning-related activity: large increases in activity after both the CS and US early in conditioning on the day of training when CRs began to increase, smaller moderate increases in activity on the following days of training, and decreases in activity after the US during asymptotic CRs. Pyramidal cell-increases declined significantly across the trials of each daily session. Theta cells showed an activity pattern opposite to the pyramidal cells, consistent with the notion that theta cells have an inhibitory influence on pyramidal cells. Single pyramidal cells also were categorized into response profiles. Most pyramidal response profiles showed increases in activity specific to the day of initial CRs. Two of the pyramidal response profiles may be involved in assessing the temporal properties of the CS-US trace conditioning trial.

Orientation selectivity in pinwheel centers in cat striate cortex

In primary visual cortex of higher mammals neurons are grouped according to their orientation preference, forming "pinwheels" around "orientation centers." Although the general structure of orientation maps is largely resolved, the microscopic arrangement of neuronal response properties in the orientation centers has remained elusive. The tetrode technique, enabling multiple single-unit recordings, in combination with intrinsic signal imaging was used to reveal the fine-grain structure of orientation maps in these locations. The results show that orientation centers represent locations where orientation columns converge containing normal, sharply tuned neurons of different orientation preference lying in close proximity.

Stimulus-dependent neuronal oscillations and local synchronization in striate cortex of the alert cat

Neuronal responses to visual stimuli that are correlated on a millisecond time scale are well documented in several areas of the mammalian visual cortex. This coherent activity often takes the form of synchronous rhythmic discharges ranging in frequency from 20 to 70 Hz. We performed experiments to determine the incidence and properties of this rhythmic activity in the striate cortex of alert cats and to compare this activity to similar data collected in the striate cortex of anesthetized cats. The results demonstrate that optimal visual stimuli evoke robust, locally synchronous, 20-70 Hz oscillatory responses in the striate cortex of cats that are fully alert and performing a visual fixation task. The oscillatory activity is stimulus dependent, largely absent during periods of spontaneous activity, and shows a systematic increase in frequency with increasing stimulus velocity. Thus, the synchronous oscillatory activity observed in this and earlier studies cannot be explained as an artifact of anesthesia nor as a phenomenon that occurs independent of visual stimulation. Rather, it is a robust process that is present in the alert state and is dependent on the presence and specific properties of visual stimuli.

Reconstructing the engram: simultaneous, multisite, many single neuron recordings

Little is known about the physiological principles that govern large-scale neuronal interactions in the mammalian brain. Here, we describe an electrophysiological paradigm capable of simultaneously recording the extracellular activity of large populations of single neurons, distributed across multiple cortical and subcortical structures in behaving and anesthetized animals. Up to 100 neurons were simultaneously recorded after 48 microwires were implanted in the brain stem, thalamus, and somatosensory cortex of rats. Overall, 86% of the implanted microwires yielded single neurons, and an average of 2.3 neurons were discriminated per microwire. Our population recordings remained stable for weeks, demonstrating that this method can be employed to investigate the dynamic and distributed neuronal ensemble interactions that underlie processes such as sensory perception, motor control, and sensorimotor learning in freely behaving animals.

Analysis of neuro-behavioral Discovery data on the Macintosh computer

This paper describes a suite of routines using IgorPro, a powerful analysis and graphing software package for the Macintosh computer, to enhance the ability to analyze, manipulate, and display data recorded with the Discovery acquisition software marketed by DataWave Technologies. The routines are able to time-align fast and slow data channels, and are especially useful for analyses that involve both neural and behavioral data. The software was designed for eyeblink conditioning and vocalization experiments, but it can easily be used for analyzing other types of neurobehavioral data. The data are first prepared on the PC with routines that inspect the header of the data file and translate the data file into a compact binary format that can be read by IgorPro. An option is also available to splice out data from unnecessary portions of an intertrial interval. The new file is then put on the Macintosh computer for display and analysis by IgorPro. These routines enable both neural and behavioral data to be quickly and easily reduced, manipulated, and statistically and graphically summarized.

Heterogeneity in local distributions of orientation-selective neurons in the cat primary visual cortex

We have employed the tetrode technique, which allows accurate discrimination of individual neuronal spike trains from multiunit recordings, in order to examine the variation of orientation selectivity among local groups of neurons. We recorded a total of 321 cells from 62 sites in area 17 of halothane-anesthetized cats; each site contained between three to ten neurons that were estimated to be less than 65 microns away from the tetrode tip. For each cell, we determined the orientation tuning in response to moving bars. Of the cells tested, 8.4% were unresponsive, 22.7% had no preferential response to any particular orientation, while 68.8% were tuned. The average difference in preferred orientation between cell pairs recorded at the same site was 10.7 deg, but the variance in preferred orientation differences differed significantly among sites. Some clusters of cells exhibited the same or nearly the same orientation preference, while others had orientation preferences that differed by as much as 90 deg. Our data demonstrate that the tuning for orientation is more heterogeneously distributed at a local level than previous studies have suggested.