The temporal structure of spike trains in the primate basal ganglia: afferent regulation of bursting demonstrated with precentral cerebral cortical ablation

We studied the temporal pattern of discharge of single units in the basal ganglia of awake primates sitting quietly. Bursting was studied with a procedure that identified individual bursts in a spike train, quantifying burst intensity (surprise), bursts per 1000 spikes, spikes per burst and burst length. Autocorrelation techniques were used to assess the dependencies of spike trains on the temporal order of intervals. Striatal units had a greater tendency to burst (79% of units) than pallidal units (50%). The caudate nucleus and putamen had nearly identical burst properties on all measures. In the pallidum, bursting was more prevalent in the external segment and bursts were more intense and more frequent than in the internal segment. The autocorrelation analysis revealed that the temporal structure of the spike train was more dependent on the order of intervals in the striatum than in the pallidum. Bursting units had an increased probability of discharge after each spike and the relative refractory period was shorter in bursting units than units without bursts. Very few units exhibited cyclic discharge properties. Ablations of areas 4 and 6 in the precentral cortex demonstrated that striatal bursting was under afferent control. The putamen, which receives more cortical afferents from areas 4 and 6 than the caudate nucleus, had fewer and less intense bursts after the afferents were lesioned. Burst intensity did not change in the pallidum after the lesion. The findings indicate that bursting properties contribute to discharge variability in the basal ganglia and suggest that information transfer in the striatum may utilize bursts. In contrast, rate coding may be a more important mechanism for units in the pallidum.     

Methods for characterizing interspike intervals and identifying bursts in neuronal activity

Neurons produce complex patterns of electrical spikes, which are often clustered in bursts. The patterns of spikes and bursts can change substantially when neurons are exposed to toxins and chemical agents. For that reason, characterization of these patterns is important for the development of neuron-based biosensors for environmental threat exposure. Here, we develop a quantitative approach to describe the distribution of interspike intervals, based on plotting histograms of the logarithm of the interspike interval. This approach provides a method for automatically classifying spikes into bursts, which does not depend on assumptions about the burst parameters. Furthermore, the approach provides a sensitive technique for detecting changes in spike and burst patterns induced by pharmacological exposure. Hence, it is suitable for use both as a research tool and for deployment in a neuron-based biosensor.     

Estimating the correlation between bursty spike trains and local field potentials

To further understand rhythmic neuronal synchronization, an increasingly useful method is to determine the relationship between the spiking activity of individual neurons and the local field potentials (LFPs) of neural ensembles. Spike field coherence (SFC) is a widely used method for measuring the synchronization between spike trains and LFPs. However, due to the strong dependency of SFC on the burst index, it is not suitable for analyzing the relationship between bursty spike trains and LFPs, particularly in high frequency bands. To address this issue, we developed a method called weighted spike field correlation (WSFC), which uses the first spike in each burst multiple times to estimate the relationship. In the calculation, the number of times that the first spike is used is equal to the spike count per burst. The performance of this method was demonstrated using simulated bursty spike trains and LFPs, which comprised sinusoids with different frequencies, amplitudes, and phases. This method was also used to estimate the correlation between pyramidal cells in the hippocampus and gamma oscillations in rats performing behaviors. Analyses using simulated and real data demonstrated that the WSFC method is a promising measure for estimating the correlation between bursty spike trains and high frequency LFPs.     

Identification of bursts in spike trains.

A computer algorithm to identify 'bursts' in trains of spikes is described. The algorithm works by constructing a histogram of interspike intervals, then analyzing the histogram to detect the critical interval value in the distribution that represents the break between short intervals within a burst and the longer intervals between bursts. When such a value is found, it is used as the 'threshold' to determine those intervals in the spike train that lie within a burst and those that lie between bursts and, thereby, to identify the beginning and end of each burst in the train. The validity of the bursts is evaluated with a chi-square test. The performance of the algorithm and how it can be assessed is discussed.     

Burst Discharge In Primary Sensory Neurons: Triggered By Subthreshold Oscillations, Maintained By Depolarizing Afterpotentials

Afferent discharge generated ectopically in the cell soma of dorsal root ganglion (DRG) neurons may play a role in normal sensation, and it contributes to paraesthesias and pain after nerve trauma. This activity is critically dependent on subthreshold membrane potential oscillations; oscillatory sinusoids that reach threshold trigger low-frequency trains of intermittent spikes. Ectopic firing may also enter a high-frequency bursting mode, however, particularly in the event of neuropathy. Bursting greatly amplifies the overall ectopic barrage. In the present report we show that subthreshold oscillations and burst discharge occur in vivo, as they do in vitro. We then show that although the first spike in each burst is triggered by an oscillatory sinusoid, firing within bursts is maintained by brief regenerative post-spike depolarizing afterpotentials (DAPs). Numerical simulations were used to identify the cellular process underlying rebound DAPs, and hence the mechanism of the spike bursts. Finally, we show that slow ramp and hold (tonic) depolarizations of the sort that occur in DRG neurons during physiologically relevant events are capable of triggering sustained ectopic bursting, but only in cells with subthreshold oscillatory behavior. Oscillations and DAPs are an essential substrate of ectopic burst discharge. Therefore, any consideration of the ways in which cellular regulation of ion channel synthesis and trafficking implement normal sensation and, when disrupted, bring about neuropathic pain must take into account the effects of this regulation on oscillations and bursting.     

Computer analysis of vestibular responses to putative neurotransmitters

An inexpensive microcomputer (Commodore 64K) based system was developed for the analysis of neural spike trains. The trains were recorded from single ampullary units in response to mechanical stimulation of the isolated semicircular canal of the bullfrog (Rana catesbeiana). A BASIC program provided a number of options while machine language subroutines generated interstimulus interval (ISI) and peristimulus time (PST) histograms. Up to thirty 5-s spike trains could be combined for analysis (0.1 ms resolution ISI, 100 ms bin width PST). Histograms and summary statistics were saved on floppy disks. The cost of adding this computer system to an existing neurophysiology laboratory is less than US $600 (printed, tape, and disk versions of these programs are available). The system was used to measure vestibular responses to putative vestibular neurotransmitters such as carbachol (an acetylcholine mimic) (Rossi et al., 1980) and gamma-aminobutyric acid (GABA) (Flock and Lam, 1974).     

Ontogenetic development of somatosensory thalum. II. Electrogenesis.

Electrodes were implanted into the brains of newborn kittens to study chronically the ontogenesis of unit spike and slow-wave phenomena in the somatosensory thalamus and adjacent structures. In general, the upper brain stem and midline (nonspecific) thalamus showed evidences of functional maturation appreciably sooner than the more laterally lying sensory nuclei. State-dependent changes in single-unit behavior began to develop relatively early (fourth postnatal day) and occasional, scattered spike bursts could be seen by the tenth to 14th day, especially during slow wave sleep. However, this scattered, low-density burst pattern, with total absence of synchrony across leads from the same electrode cluster, at a time when runs of thalamic slow waves or spindles are already very obvious, casts considerable doubt on any obligatory relationship between these two phenomena. Alternative functional models for the production of thalamic wave trains are considered.     

Dopamine modulation of spike dynamics in bursting neurons

The pyloric network of the lobster stomatogastric ganglion is a prime example of an oscillatory neural circuit. In our previous study on the firing patterns of pyloric neurons we observed characteristic temporal structures termed 'interspike interval (ISI) signatures' which were found to depend on the synaptic connectivity of the network. Dopamine, a well-known modulator of the pyloric network, is known to affect inhibitory synapses so it might also tune the fine temporal structure of intraburst spikes, a phenomenon not previously investigated. In the recent work we study the DA modulation of ISI patterns of two identified pyloric neurons in normal conditions and after blocking their glutamatergic synaptic connections. Dopamine (10-50 microM) strongly regularizes the firing of the lateral pyloric (LP) and pyloric dilator (PD) neurons by increasing the reliability of recurrent spike patterns. The most dramatic effect is observed in the LP, where precisely replicated spike multiplets appear in a normally 'noisy' neuron. The DA-induced regularization of intraburst spike patterns requires functional glutamatergic inputs to the LP neuron and this effect cannot be mimicked by simple intracellular depolarization. Inhibitory synaptic inputs arriving before the bursts are important factors in shaping the intraburst spike dynamics of both the PD and the LP neurons. Our data reveal a novel aspect of chemical neuromodulation in oscillatory neural networks. This effect sets in at concentrations lower than those affecting the overall burst pattern of the network. The sensitivity of intraburst spike dynamics to preceding synaptic inputs also suggests a novel method of temporal coding in neural bursters.     

Bursting properties of units in cat globus pallidus and entopeduncular nucleus: the effect of excitotoxic striatal lesions

The bursting properties of units recorded in globus pallidus and entopeduncular nucleus were studied in awake cats sitting quietly before and after ipsilateral excitotoxic striatal lesions. A computerized statistical procedure was used to identify and evaluate bursts in the recorded spike trains. Bursts were assigned a quantitative statistical measure of burst 'strength' (or improbability) - the surprise value. Before the lesion, 34% of units in the globus pallidus and 60% of units in the entopeduncular nucleus exhibited bursts. Burst units had a significantly slower discharge rate and a significantly greater variability of discharge than non-burst units. The mean length of the interspike intervals immediately preceding the bursts was significantly longer than the overall median intervals in burst units. After the lesion, 21% of units in the globus pallidus and 11% of the units in the entopeduncular nucleus exhibited bursts. Burst units had significantly higher discharge rates and lower discharge variability after the lesion. In contrast, the lesion had no significant effect on the rate or variability of non-burst units. The differences between bursting and non-bursting units in discharge rate and variability disappeared after the lesion. In globus pallidus, the lesion resulted in a significant reduction in the mean number of bursts per unit, surprise value per burst, mean length of bursts, and number of spikes per burst, and a significant increase in the mean discharge rate of burst units. In entopeduncular nucleus, the small number of bursts recorded after the lesion precluded a useful statistical comparison of the effect of striatal lesions on the properties of the bursts. This study demonstrates that removing striatal projections to globus pallidus and entopeduncular nucleus decreases bursting in these nuclei, indicating that intact striatal projections are necessary for the normal production of bursts in these regions.     

A simple indicator of nonstationarity of firing rate in spike trains

The assessment of stationarity of firing rate in neural spike trains is important but is often performed only visually. Facing the growing amount of neural data generated by multi-electrode recording, there is a need for an automatic method to identify and disqualify spike trains with highly nonstationary firing rates. In this report, we propose a simple test of nonstationarity, associated with an indicator quantifying the degree of nonstationary in a spike train. This method is compared to the Mann-Kendall test of trend detection and the Runs test on simulated and real spike trains.     

Burst firing of oxytocin neurons in male rat hypothalamic slices

Burst firing and single spike activity play different roles in the modulation of local neuronal circuit activity and neurosecretion. In hypothalamic oxytocin (OT) neurons in vivo, burst firing is associated with pulsatile secretion of OT in the milk ejection reflex, and can be observed in slices from both immature and lactating rats in vitro. Whether OT neurons from male rats also possess burst firing capability is still an open question. To examine this possibility, whole-cell patch clamp recordings were made in supraoptic nucleus OT neurons in brain slices from male rats. In low Ca(2+) medium, the alpha(1)-adrenoceptor agonist, phenylephrine evoked bursts that were highly similar to those from lactating rats in vivo and in vitro: explosive onset, short-duration, quickly reaching peak firing rate and displaying an exponential decay in returning to the pre-burst rate. During bursts, spike durations increased, and spike amplitudes decreased, while riding on an arc of depolarization around peak rate. In comparison to those from lactating rats in vitro, the rising phase of male bursts was more rapid, the decay phase was slower, and the rising phase of the spike after hyperpolarization was faster. No significant differences, however, were seen in burst characteristics that are most important in determining the amount of peptide release: burst amplitudes (the number of spikes in a burst), firing frequency within bursts or peak firing rate. Thus, we conclude that OT neurons in males are capable of burst firing highly similar to that seen in lactating rats.     

Spike detection in human muscle sympathetic nerve activity using the kurtosis of stationary wavelet transform coefficients

The accurate assessment of autonomic sympathetic function is important in the diagnosis and study of various autonomic and cardiovascular disorders. Sympathetic function in humans can be assessed by recording the muscle sympathetic nerve activity, which is characterized by synchronous neuronal discharges separated by periods of neural silence dominated by colored Gaussian noise. The raw nerve activity is generally rectified, integrated, and quantified using the integrated burst rate or area. We propose an alternative quantification involving spike detection using a two-stage stationary wavelet transform (SWT) de-noising method. The SWT coefficients are first separated into noise-related and burst-related coefficients on the basis of their local kurtosis. The noise-related coefficients are then used to establish a threshold to identify spikes within the bursts. This method demonstrated better detection performance than an unsupervised amplitude discriminator and similar wavelet-based methods when confronted with simulated data of varying burst rate and signal to noise ratio. Additional validation on data acquired during a graded head-up tilt protocol revealed a strong correlation between the mean spike rate and the mean integrate burst rate (r = 0.85) and burst area rate (r = 0.91). In conclusion, the kurtosis-based wavelet de-noising technique is a potentially useful method of studying sympathetic nerve activity in humans.     

Short-term plasticity optimizes synaptic information transmission

Short-term synaptic plasticity (STP) is widely thought to play an important role in information processing. This major function of STP has recently been challenged, however, by several computational studies indicating that transmission of information by dynamic synapses is broadband, i.e., frequency independent. Here we developed an analytical approach to quantify time- and rate-dependent synaptic information transfer during arbitrary spike trains using a realistic model of synaptic dynamics in excitatory hippocampal synapses. We found that STP indeed increases information transfer in a wide range of input rates, which corresponds well to the naturally occurring spike frequencies at these synapses. This increased information transfer is observed both during Poisson-distributed spike trains with a constant rate and during naturalistic spike trains recorded in hippocampal place cells in exploring rodents. Interestingly, we found that the presence of STP in low release probability excitatory synapses leads to optimization of information transfer specifically for short high-frequency bursts, which are indeed commonly observed in many excitatory hippocampal neurons. In contrast, more reliable high release probability synapses that express dominant short-term depression are predicted to have optimal information transmission for single spikes rather than bursts. This prediction is verified in analyses of experimental recordings from high release probability inhibitory synapses in mouse hippocampal slices and fits well with the observation that inhibitory hippocampal interneurons do not commonly fire spike bursts. We conclude that STP indeed contributes significantly to synaptic information transfer and may serve to maximize information transfer for specific firing patterns of the corresponding neurons.     

Phasic bursts in rat magnocellular neurosecretory cells are not intrinsically regenerative in vivo

Vasopressinergic hypothalamic magnocellular neurosecretory cells fire in phasic bursts. Burst initiation involves summation of postsynaptic potentials to generate action potentials. Action potentials are each followed by a nonsynaptic depolarizing after-potential that summates temporally to generate a plateau potential and so sustain activity throughout the burst. It is unknown whether this plateau potential exceeds spike threshold in vivo to cause intrinsic regenerative firing or simply approaches threshold to increase the probability that excitatory postsynaptic potentials will trigger further action potentials. Here we show that pharmacological blockade of ionotropic glutamatergic transmission by microdialysis application of kynurenic acid into the supraoptic nucleus of anaesthetized rats prevents spontaneous bursts and bursts (after-discharge) evoked by short trains of antidromically stimulated action potentials in magnocellular neurosecretory cells. Even during prolonged depolarization induced by 1 m NaCl infusion, kynurenic acid microdialysis application still blocked after-discharge. The ability of kynurenic acid to block after-discharge during osmotic stimulation was not caused by an unmasking of inhibitory postsynaptic potentials as kynurenic acid was equally effective in the presence of the ionotropic gamma-aminobutyric acid receptor antagonist, bicuculline, nor did it result from inhibition of plateau potential amplitude as this was unaffected by kynurenic acid and bicuculline in vitro, as was after-discharge evoked in vitro. We conclude that phasic bursts are nonregenerative in vivo but rather require continued excitatory synaptic input activity superimposed upon a subthreshold plateau potential to sustain burst activity.     

Learning-related patterns of CA1 spike trains parallel stimulation parameters optimal for inducing hippocampal long-term potentiation.

Recent studies have revealed 3 stimulation parameters that together comprise the temporal pattern of neuronal activation optimal for the induction of hippocampal LTP: high-frequency bursts, activity 100-200 ms prior to a burst, and burst delivery in phase with the ongoing hippocampal theta rhythm. The present paper reports that these 3 aspects of patterned neural activity, collectively referred to as "theta-bursting," are characteristic of the spike trains of CA1 pyramidal cells in rats during the sampling and analysis of learning cues in an odor discrimination task and during performances of a spatial memory task. In contrast, theta-bursting occurs relatively infrequently during behavioral events less directly related to task-relevant mnemonic processing. These findings suggest that the optimal conditions for the induction of LTP occur naturally in behaving animals, time-locked to behavioral events critical to learning.     

Classification of extracellularly recorded neurons by their discharge patterns and their correlates with intracellularly identified neuronal types in the frontal cortex of behaving monkeys

Neurons in the cerebral cortex are not homogeneous. However, neuronal types have been ignored in most previous work studying neuronal processes in behaving monkeys. We propose a new method to identify neuronal types in extracellular recording studies of behaving monkeys. We classified neurons as either bursting or non‐bursting, and then classified the bursting neurons into three types: (i) neurons displaying a burst of many spikes (maximum number of spikes within a burst; NSB max ≥ 8) at a high discharge rate (maximum interspike interval; ISI max < 5 ms); (ii) neurons displaying a burst of fewer spikes (NSB max ≤ 5) at a high discharge rate (ISI max < 5 ms); and (iii) neurons displaying a burst of a few spikes (NSB max ≤ 7) at relatively long ISIs (ISI max > 5 ms). We found that the discharge patterns of the four groups corresponded to those of regular spiking (RS), fast spiking (FS), fast rhythmic bursting (FRB) and intrinsic bursting (IB) neurons demonstrated in intracellular recording studies using in vitro slice preparations, respectively. In addition, we examined correlations with the task events for neurons recorded in the frontal eye field and neuronal interactions for pairs of neurons recorded simultaneously from a single electrode. We found that they were substantially different between RS and FS types. These results suggest that neurons in the frontal cortex of behaving monkeys can be classified into four types based on their discharge patterns, and that these four types contribute differentially to cortical operations.     

Physiological effects of sustained blockade of excitatory synaptic transmission on spontaneously active developing neuronal networks--an inquiry into the reciprocal linkage between intrinsic biorhythms and neuroplasticity in early ontogeny

Spontaneous bioelectric activity (SBA) taking the form of extracellularly recorded spike trains (SBA) has been quantitatively analyzed in organotypic neonatal rat visual cortex explants at different ages in vitro, and the effects investigated of both short- and long-term pharmacological suppression of glutamatergic synaptic transmission. In the presence of APV, a selective NMDA receptor blocker, 1-2- (but not 3-)week-old cultures recovered their previous SBA levels in a matter of hours, although in imitation of the acute effect of the GABAergic inhibitor picrotoxin (PTX), bursts of action potentials were abnormally short and intense. Cultures treated either overnight or chronically for 1-3 weeks with APV, the AMPA/kainate receptor blocker DNQX, or a combination of the two were found to display very different abnormalities in their firing patterns. NMDA receptor blockade for 3 weeks produced the most severe deviations from control SBA, consisting of greatly prolonged and intensified burst firing with a strong tendency to be broken up into trains of shorter spike clusters. This pattern was most closely approximated by acute GABAergic disinhibition in cultures of the same age, but this latter treatment also differed in several respects from the chronic-APV effect. In 2-week-old explants, in contrast, it was the APV+DNQX treated group which showed the most exaggerated spike bursts. Functional maturation of neocortical networks, therefore, may specifically require NMDA receptor activation (not merely a high level of neuronal firing) which initially is driven by endogenous rather than afferent evoked bioelectric activity. Putative cellular mechanisms are discussed in the context of a thorough review of the extensive but scattered literature relating activity-dependent brain development to spontaneous neuronal firing patterns.