Oscillations and synchrony in basal ganglia circuits may play a key role in the organization of voluntary actions and habits. We recorded single units and local field potentials from multiple striatal and cortical locations simultaneously, over a range of behavioral states. We observed opposite gradients of oscillatory entrainment, with dorsal/lateral striatal neurons entrained to high-voltage spindle oscillations ("spike wave discharges") and ventral/medial striatal neurons entrained to the hippocampal theta rhythm. While the majority of units were likely medium-spiny projection neurons, a second neuronal population showed characteristic features of fast-spiking GABAergic interneurons, including tonic activity, brief waveforms, and high-frequency bursts. These fired at an earlier spindle phase than the main neuronal population, and their density within striatum corresponded closely to the intensity of spindle oscillations. The orchestration of oscillatory activity by networks of striatal interneurons may be an important mechanism in the pathophysiology of neurological disorders such as Parkinson's disease.

Analyzing the dependencies between spike trains is an important step in understanding how neurons work in concert to represent biological signals. Usually this is done for pairs of neurons at a time using correlation-based techniques. Chornoboy, Schramm, and Karr (1988) proposed maximum likelihood methods for the simultaneous analysis of multiple pair-wise interactions among an ensemble of neurons. One of these methods is an iterative, continuous-time estimation algorithm for a network likelihood model formulated in terms of multiplicative conditional intensity functions. We devised a discrete-time version of this algorithm that includes a new, efficient computational strategy, a principled method to compute starting values, and a principled stopping criterion. In an analysis of simulated neural spike trains from ensembles of interacting neurons, the algorithm recovered the correct connectivity matrices and interaction parameters. In the analysis of spike trains from an ensemble of rat hippocampal place cells, the algorithm identified a connectivity matrix and interaction parameters consistent with the pattern of conjoined firing predicted by the overlap of the neurons' spatial receptive fields. These results suggest that the network likelihood model can be an efficient tool for the analysis of ensemble spiking activity.

No abstract

Temporal difference learning (TD) is a popular algorithm in machine learning. Two learning signals that are derived from this algorithm, the predictive value and the prediction error, have been shown to explain changes in neural activity and behavior during learning across species. Here, the predictive value signal is used to explain the time course of learning-related changes in the activity of hippocampal neurons in monkeys performing an associative learning task. The TD algorithm serves as the centerpiece of a joint probability model for the learning-related neural activity and the behavioral responses recorded during the task. The neural component of the model consists of spiking neurons that compete and learn the reward-predictive value of task-relevant input signals. The predictive-value signaled by these neurons influences the behavioral response generated by a stochastic decision stage, which constitutes the behavioral component of the model. It is shown that the time course of the changes in neural activity and behavioral performance generated by the model exhibits key features of the experimental data. The results suggest that information about correct associations may be expressed in the hippocampus before it is detected in the behavior of a subject. In this way, the hippocampus may be among the earliest brain areas to express learning and drive the behavioral changes associated with learning. Correlates of reward-predictive value may be expressed in the hippocampus through rate remapping within spatial memory representations, they may represent reward-related aspects of a declarative or explicit relational memory representation of task contingencies, or they may correspond to reward-related components of episodic memory representations. These potential functions are discussed in connection with hippocampal cell assembly sequences and their reverse reactivation during the awake state. The results provide further support for the proposal that neural processes underlying learning may be implementing a temporal difference-like algorithm.

Abstract:We describe a method for computing a pair of spike detection thresholds, called 'truncation thresholds', using truncated probability distributions, for extracellular recordings. In existing methods the threshold is usually set to a multiple of an estimate of the standard deviation of the noise in the recording, with the multiplication factor being chosen between 3 and 5 according to the researcher's preferences. Our method has the following advantages over these methods. First, because the standard deviation is usually estimated from the entire recording, which includes the spikes, it increases with firing rate. By contrast, truncation thresholds decrease in absolute value with increasing firing rate, thereby capturing more of the signal. Second, the parameters of the selected noise distribution are estimated more accurately by maximum likelihood fitting of the truncated distribution to the data delimited by the truncation thresholds. Third, the computation of the truncation thresholds is completely data-driven. It does not involve a userdefined multiplication factor. Fourth, methods that use a threshold that is proportional to the estimated standard deviation of the noise assume that the noise distribution is symmetrical around the mean. By contrast, truncation thresholds are not linked to each other by an assumption of symmetry about some axis. Fifth, existing methods do not verify that subthreshold data obey a noise distribution. Truncation thresholds, however, are defined by the fact that the distribution of the data they delimit is statistically indistinguishable, according to the Kolmogorov-Smirnov test, from a selected distribution, truncated at those thresholds. Application of the method is illustrated using recordings from cortical area M1 in awake behaving rats, as well as in simulated recordings. Source code and executables of a software suite that computes the truncation thresholds are provided for the case when the noise distribution is modeled as truncated normal.