Hippocampal "time cells" bridge the gap in memory for discontiguous events

The hippocampus is critical to remembering the flow of events in distinct experiences and, in doing so, bridges temporal gaps between discontiguous events. Here, we report a robust hippocampal representation of sequence memories, highlighted by "time cells" that encode successive moments during an empty temporal gap between the key events, while also encoding location and ongoing behavior. Furthermore, just as most place cells "remap" when a salient spatial cue is altered, most time cells form qualitatively different representations ("retime") when the main temporal parameter is altered. Hippocampal neurons also differentially encode the key events and disambiguate different event sequences to compose unique, temporally organized representations of specific experiences. These findings suggest that hippocampal neural ensembles segment temporally organized memories much the same as they represent locations of important events in spatially defined environments.

A point process framework for relating neural spiking activity to spiking history, neural ensemble, and extrinsic covariate effects

Multiple factors simultaneously affect the spiking activity of individual neurons. Determining the effects and relative importance of these factors is a challenging problem in neurophysiology. We propose a statistical framework based on the point process likelihood function to relate a neuron's spiking probability to three typical covariates: the neuron's own spiking history, concurrent ensemble activity, and extrinsic covariates such as stimuli or behavior. The framework uses parametric models of the conditional intensity function to define a neuron's spiking probability in terms of the covariates. The discrete time likelihood function for point processes is used to carry out model fitting and model analysis. We show that, by modeling the logarithm of the conditional intensity function as a linear combination of functions of the covariates, the discrete time point process likelihood function is readily analyzed in the generalized linear model (GLM) framework. We illustrate our approach for both GLM and non-GLM likelihood functions using simulated data and multivariate single-unit activity data simultaneously recorded from the motor cortex of a monkey performing a visuomotor pursuit-tracking task. The point process framework provides a flexible, computationally efficient approach for maximum likelihood estimation, goodness-of-fit assessment, residual analysis, model selection, and neural decoding. The framework thus allows for the formulation and analysis of point process models of neural spiking activity that readily capture the simultaneous effects of multiple covariates and enables the assessment of their relative importance.

Dynamic analysis of neural encoding by point process adaptive filtering

Neural receptive fields are dynamic in that with experience, neurons change their spiking responses to relevant stimuli. To understand how neural systems adapt their representations of biological information, analyses of receptive field plasticity from experimental measurements are crucial. Adaptive signal processing, the well-established engineering discipline for characterizing the temporal evolution of system parameters, suggests a framework for studying the plasticity of receptive fields. We use the Bayes' rule Chapman-Kolmogorov paradigm with a linear state equation and point process observation models to derive adaptive filters appropriate for estimation from neural spike trains. We derive point process filter analogues of the Kalman filter, recursive least squares, and steepest-descent algorithms and describe the properties of these new filters. We illustrate our algorithms in two simulated data examples. The first is a study of slow and rapid evolution of spatial receptive fields in hippocampal neurons. The second is an adaptive decoding study in which a signal is decoded from ensemble neural spiking activity as the receptive fields of the neurons in the ensemble evolve. Our results provide a paradigm for adaptive estimation for point process observations and suggest a practical approach for constructing filtering algorithms to track neural receptive field dynamics on a millisecond timescale.

Biophysical foundations underlying TMS: setting the stage for an effective use of neurostimulation in the cognitive neurosciences

Transcranial Magnetic Stimulation (TMS) induces electrical currents in the brain to stimulate neural tissue. This article reviews our present understanding of TMS methodology, focusing on its biophysical foundations. We concentrate on how the laws of electromagnetic induction apply to TMS; addressing issues such as the location, area (i.e., focality), depth, and mechanism of TMS. We also present a review of the present limitations and future potential of the technique.

Hippocampal replay reflects specific past experiences rather than a plan for subsequent choice

Executing memory-guided behavior requires storage of information about experience and later recall of that information to inform choices. Awake hippocampal replay, when hippocampal neural ensembles briefly reactivate a representation related to prior experience, has been proposed to critically contribute to these memory-related processes. However, it remains unclear whether awake replay contributes to memory function by promoting the storage of past experiences, facilitating planning based on evaluation of those experiences, or both. We designed a dynamic spatial task that promotes replay before a memory-based choice and assessed how the content of replay related to past and future behavior. We found that replay content was decoupled from subsequent choice and instead was enriched for representations of previously rewarded locations and places that had not been visited recently, indicating a role in memory storage rather than in directly guiding subsequent behavior.

Transcranial magnetic stimulation and stroke: A computer-based human model study

This paper explores how transcranial magnetic stimulation (TMS) induced currents in the brain are perturbed by electrical and anatomical changes following a stroke in its chronic stage. Multiple MRI derived finite element head models were constructed and evaluated to address the effects that strokes can have on the induced stimulating TMS currents by comparing stroke models of various sizes and geometries to a healthy head model under a number of stimulation conditions. The TMS induced currents were significantly altered for stimulation proximal to the lesion site in all of the models analyzed. The current density distributions were modified in magnitude, location, and orientation such that the population of neural elements that are stimulated will be correspondingly altered. The current perturbations were minimized for conditions tested where the coil was far removed from the lesion site, including models of stimulation contralateral to the lesioned hemisphere. The present limitations of TMS to the peri-lesional cortex are explored, ultimately concluding that conventional clinical standards for stimulation are unreliable and potentially dangerous predictors of the site and degree of stimulation when TMS is applied proximal to infarction site.

Network inference with confidence from multivariate time series

Networks--collections of interacting elements or nodes--abound in the natural and manmade worlds. For many networks, complex spatiotemporal dynamics stem from patterns of physical interactions unknown to us. To infer these interactions, it is common to include edges between those nodes whose time series exhibit sufficient functional connectivity, typically defined as a measure of coupling exceeding a predetermined threshold. However, when uncertainty exists in the original network measurements, uncertainty in the inferred network is likely, and hence a statistical propagation of error is needed. In this manuscript, we describe a principled and systematic procedure for the inference of functional connectivity networks from multivariate time series data. Our procedure yields as output both the inferred network and a quantification of uncertainty of the most fundamental interest: uncertainty in the number of edges. To illustrate this approach, we apply a measure of linear coupling to simulated data and electrocorticogram data recorded from a human subject during an epileptic seizure. We demonstrate that the procedure is accurate and robust in both the determination of edges and the reporting of uncertainty associated with that determination.

Analysis of between-trial and within-trial neural spiking dynamics

Recording single-neuron activity from a specific brain region across multiple trials in response to the same stimulus or execution of the same behavioral task is a common neurophysiology protocol. The raster plots of the spike trains often show strong between-trial and within-trial dynamics, yet the standard analysis of these data with the peristimulus time histogram (PSTH) and ANOVA do not consider between-trial dynamics. By itself, the PSTH does not provide a framework for statistical inference. We present a state-space generalized linear model (SS-GLM) to formulate a point process representation of between-trial and within-trial neural spiking dynamics. Our model has the PSTH as a special case. We provide a framework for model estimation, model selection, goodness-of-fit analysis, and inference. In an analysis of hippocampal neural activity recorded from a monkey performing a location-scene association task, we demonstrate how the SS-GLM may be used to answer frequently posed neurophysiological questions including, What is the nature of the between-trial and within-trial task-specific modulation of the neural spiking activity? How can we characterize learning-related neural dynamics? What are the timescales and characteristics of the neuron's biophysical properties? Our results demonstrate that the SS-GLM is a more informative tool than the PSTH and ANOVA for analysis of multiple trial neural responses and that it provides a quantitative characterization of the between-trial and within-trial neural dynamics readily visible in raster plots, as well as the less apparent fast (1-10 ms), intermediate (11-20 ms), and longer (>20 ms) timescale features of the neuron's biophysical properties.

A state-space analysis for reconstruction of goal-directed movements using neural signals

The execution of reaching movements involves the coordinated activity of multiple brain regions that relate variously to the desired target and a path of arm states to achieve that target. These arm states may represent positions, velocities, torques, or other quantities. Estimation has been previously applied to neural activity in reconstructing the target separately from the path. However, the target and path are not independent. Because arm movements are limited by finite muscle contractility, knowledge of the target constrains the path of states that leads to the target. In this letter, we derive and illustrate a state equation to capture this basic dependency between target and path. The solution is described for discrete-time linear systems and gaussian increments with known target arrival time. The resulting analysis enables the use of estimation to study how brain regions that relate variously to target and path together specify a trajectory. The corresponding reconstruction procedure may also be useful in brain-driven prosthetic devices to generate control signals for goal-directed movements.

General-Purpose Filter Design for Neural Prosthetic Devices

Brain-driven interfaces depend on estimation procedures to convert neural signals to inputs for prosthetic devices that can assist individuals with severe motor deficits. Previous estimation procedures were developed on an application-specific basis. Here we report a coherent estimation framework that unifies these procedures and motivates new applications of prosthetic devices driven by action potentials, local field potentials (LFPs), electrocorticography (ECoG), electroencephalography (EEG), electromyography (EMG), or optical methods. The brain-driven interface is described as a probabilistic relationship between neural activity and components of a prosthetic device that may take on discrete or continuous values. A new estimation procedure is developed for action potentials, and a corresponding procedure is described for field potentials and optical measurements. We test our framework against dominant approaches in an arm reaching task using simulated traces of ensemble spiking activity from primary motor cortex (MI) and a wheelchair navigation task using simulated traces of EEG-band power. Adaptive filtering is incorporated to demonstrate performance under neuron death and discovery. Finally, we characterize performance under model misspecification using physiologically realistic history dependence in MI spiking. These simulated results predict that the unified framework outperforms previous approaches under various conditions, in the control of position and velocity, based on trajectory and endpoint mean squared errors.

Using point process models to compare neural spiking activity in the subthalamic nucleus of Parkinson's patients and a healthy primate

Placement of deep brain stimulating electrodes in the subthalamic nucleus (STN) to treat Parkinson's disease (PD) also allows the recording of single neuron spiking activity. Analyses of these unique data offer an important opportunity to better understand the pathophysiology of PD. Despite the point process nature of PD neural spiking activity, point process methods are rarely used to analyze these recordings. We develop a point process representation of PD neural spiking activity using a generalized linear model to describe long- and short-term temporal dependencies in the spiking activity of 28 STN neurons from seven PD patients and 35 neurons from one healthy primate (surrogate control) recorded, while the subjects executed a directed-hand movement task. We used the point process model to characterize each neuron's bursting, oscillatory, and directional tuning properties during key periods in the task trial. Relative to the control neurons, the PD neurons showed increased bursting, increased 10-30 Hz oscillations, and increased fluctuations in directional tuning. These features, which traditional methods failed to capture accurately, were efficiently summarized in a single model in the point process analysis of each neuron. The point process framework suggests a useful approach for developing quantitative neural correlates that may be related directly to the movement and behavioral disorders characteristic of PD.

Transcranial magnetic stimulation and brain atrophy: a computer-based human brain model study

This paper is aimed at exploring the effect of cortical brain atrophy on the currents induced by transcranial magnetic stimulation (TMS). We compared the currents induced by various TMS conditions on several different MRI derived finite element head models of brain atrophy, incorporating both decreasing cortical volume and widened sulci. The current densities induced in the cortex were dependent upon the degree and type of cortical atrophy and were altered in magnitude, location, and orientation when compared to healthy head models. Predictive models of the degree of current density attenuation as a function of the scalp-to-cortex distance were analyzed, concluding that those which ignore the electromagnetic field-tissue interactions lead to inaccurate conclusions. Ultimately, the precise site and population of neural elements stimulated by TMS in an atrophic brain cannot be predicted based on healthy head models which ignore the effects of the altered cortex on the stimulating currents. Clinical applications of TMS should be carefully considered in light of these findings.

Scalp recorded spike ripples predict seizure risk in childhood epilepsy better than spikes

In the past decade, brief bursts of fast oscillations in the ripple range have been identified in the scalp EEG as a promising non-invasive biomarker for epilepsy. However, investigation and clinical application of this biomarker have been limited because standard approaches to identify these brief, low amplitude events are difficult, time consuming, and subjective. Recent studies have demonstrated that ripples co-occurring with epileptiform discharges ('spike ripple events') are easier to detect than ripples alone and have greater pathological significance. Here, we used objective techniques to quantify spike ripples and test whether this biomarker predicts seizure risk in childhood epilepsy. We evaluated spike ripples in scalp EEG recordings from a prospective cohort of children with a self-limited epilepsy syndrome, benign epilepsy with centrotemporal spikes, and healthy control children. We compared the rate of spike ripples between children with epilepsy and healthy controls, and between children with epilepsy during periods of active disease (active, within 1 year of seizure) and after a period of sustained seizure-freedom (seizure-free, >1 year without seizure), using semi-automated and automated detection techniques. Spike ripple rate was higher in subjects with active epilepsy compared to healthy controls (P = 0.0018) or subjects with epilepsy who were seizure-free ON or OFF medication (P = 0.0018). Among epilepsy subjects with spike ripples, each month seizure-free decreased the odds of a spike ripple by a factor of 0.66 [95% confidence interval (0.47, 0.91), P = 0.021]. Comparing the diagnostic accuracy of the presence of at least one spike ripple versus a classic spike event to identify group, we found comparable sensitivity and negative predictive value, but greater specificity and positive predictive value of spike ripples compared to spikes (P = 0.016 and P = 0.006, respectively). We found qualitatively consistent results using a fully automated spike ripple detector, including comparison with an automated spike detector. We conclude that scalp spike ripple events identify disease and track with seizure risk in this epilepsy population, using both semi-automated and fully automated detection methods, and that this biomarker outperforms analysis of spikes alone in categorizing seizure risk. These data provide evidence that spike ripples are a specific non-invasive biomarker for seizure risk in benign epilepsy with centrotemporal spikes and support future work to evaluate the utility of this biomarker to guide medication trials and tapers in these children and predict seizure risk in other at-risk populations.

Successful reconstruction of a physiological circuit with known connectivity from spiking activity alone

Identifying the structure and dynamics of synaptic interactions between neurons is the first step to understanding neural network dynamics. The presence of synaptic connections is traditionally inferred through the use of targeted stimulation and paired recordings or by post-hoc histology. More recently, causal network inference algorithms have been proposed to deduce connectivity directly from electrophysiological signals, such as extracellularly recorded spiking activity. Usually, these algorithms have not been validated on a neurophysiological data set for which the actual circuitry is known. Recent work has shown that traditional network inference algorithms based on linear models typically fail to identify the correct coupling of a small central pattern generating circuit in the stomatogastric ganglion of the crab Cancer borealis. In this work, we show that point process models of observed spike trains can guide inference of relative connectivity estimates that match the known physiological connectivity of the central pattern generator up to a choice of threshold. We elucidate the necessary steps to derive faithful connectivity estimates from a model that incorporates the spike train nature of the data. We then apply the model to measure changes in the effective connectivity pattern in response to two pharmacological interventions, which affect both intrinsic neural dynamics and synaptic transmission. Our results provide the first successful application of a network inference algorithm to a circuit for which the actual physiological synapses between neurons are known. The point process methodology presented here generalizes well to larger networks and can describe the statistics of neural populations. In general we show that advanced statistical models allow for the characterization of effective network structure, deciphering underlying network dynamics and estimating information-processing capabilities.

To appreciate how neural circuits control behaviors, we must understand two things. First, how the neurons comprising the circuit are connected, and second, how neurons and their connections change after learning or in response to neuromodulators. Neuronal connectivity is difficult to determine experimentally, whereas neuronal activity can often be readily measured. We describe a statistical model to estimate circuit connectivity directly from measured activity patterns. We use the timing relationships between observed spikes to predict synaptic interactions between simultaneously observed neurons. The model estimate provides each predicted connection with a curve that represents how strongly, and at which temporal delays, one circuit element effectively influences another. These curves are analogous to synaptic interactions of the level of the membrane potential of biological neurons and share some of their features such as being inhibitory or excitatory. We test our method on recordings from the pyloric circuit in the crab stomatogastric ganglion, a small circuit whose connectivity is completely known beforehand, and find that the predicted circuit matches the biological one — a result other techniques failed to achieve. In addition, we show that drug manipulations impacting the circuit are revealed by this technique. These results illustrate the utility of our analysis approach for inferring connections from neural spiking activity.

Stable encoding of task structure coexists with flexible coding of task events in sensorimotor striatum

The sensorimotor striatum, as part of the brain's habit circuitry, has been suggested to store fixed action values as a result of stimulus-response learning and has been contrasted with a more flexible system that conditionally assigns values to behaviors. The stability of neural activity in the sensorimotor striatum is thought to underlie not only normal habits but also addiction and clinical syndromes characterized by behavioral fixity. By recording in the sensorimotor striatum of mice, we asked whether neuronal activity acquired during procedural learning would be stable even if the sensory stimuli triggering the habitual behavior were altered. Contrary to expectation, both fixed and flexible activity patterns appeared. One, representing the global structure of the acquired behavior, was stable across changes in task cuing. The second, a fine-grain representation of task events, adjusted rapidly. Such dual forms of representation may be critical to allow motor and cognitive flexibility despite habitual performance.

Drawing inferences from Fano factor calculations

An important characterization of neural spiking is the ratio of the variance to the mean of the spike counts in a set of intervals--the Fano factor. For a Poisson process, the theoretical Fano factor is exactly one. For simulated or experimental neural data, the sample Fano factor is never exactly one, but often appears close to one. In this short communication, we characterize the distribution of the Fano factor for a Poisson process, allowing us to compute probability bounds and perform hypothesis tests on the distribution of recorded neural spike counts. We show that for a Poisson process the Fano factor asymptotically follows a gamma distribution with dependence on the number of observations of spike counts, and that convergence to this asymptotic distribution is fast. The analysis provides a simple method to determine how close to 1 the computed Fano factor should be and to formally test whether the observed variability in the spiking is likely to arise in data generated by a Poisson process.

Clusterless Decoding of Position from Multiunit Activity Using a Marked Point Process Filter

Point process filters have been applied successfully to decode neural signals and track neural dynamics. Traditionally these methods assume that multiunit spiking activity has already been correctly spike-sorted. As a result, these methods are not appropriate for situations where sorting cannot be performed with high precision, such as real-time decoding for brain-computer interfaces. Because the unsupervised spike-sorting problem remains unsolved, we took an alternative approach that takes advantage of recent insights into clusterless decoding. Here we present a new point process decoding algorithm that does not require multiunit signals to be sorted into individual units. We use the theory of marked point processes to construct a function that characterizes the relationship between a covariate of interest (in this case, the location of a rat on a track) and features of the spike waveforms. In our example, we use tetrode recordings, and the marks represent a four-dimensional vector of the maximum amplitudes of the spike waveform on each of the four electrodes. In general, the marks may represent any features of the spike waveform. We then use Bayes's rule to estimate spatial location from hippocampal neural activity. We validate our approach with a simulation study and experimental data recorded in the hippocampus of a rat moving through a linear environment. Our decoding algorithm accurately reconstructs the rat's position from unsorted multiunit spiking activity. We then compare the quality of our decoding algorithm to that of a traditional spike-sorting and decoding algorithm. Our analyses show that the proposed decoding algorithm performs equivalent to or better than algorithms based on sorted single-unit activity. These results provide a path toward accurate real-time decoding of spiking patterns that could be used to carry out content-specific manipulations of population activity in hippocampus or elsewhere in the brain.

Construction of point process adaptive filter algorithms for neural systems using sequential Monte Carlo methods

The stochastic state point process filter (SSPPF) and steepest descent point process filter (SDPPF) are adaptive filter algorithms for state estimation from point process observations that have been used to track neural receptive field plasticity and to decode the representations of biological signals in ensemble neural spiking activity. The SSPPF and SDPPF are constructed using, respectively, Gaussian and steepest descent approximations to the standard Bayes and Chapman-Kolmogorov (BCK) system of filter equations. To extend these approaches for constructing point process adaptive filters, we develop sequential Monte Carlo (SMC) approximations to the BCK equations in which the SSPPF and SDPPF serve as the proposal densities. We term the two new SMC point process filters SMC-PPFs and SMC-PPFD, respectively. We illustrate the new filter algorithms by decoding the wind stimulus magnitude from simulated neural spiking activity in the cricket cercal system. The SMC-PPFs and SMC-PPFD provide more accurate state estimates at low number of particles than a conventional bootstrap SMC filter algorithm in which the state transition probability density is the proposal density. We also use the SMC-PPFs algorithm to track the temporal evolution of a spatial receptive field of a rat hippocampal neuron recorded while the animal foraged in an open environment. Our results suggest an approach for constructing point process adaptive filters using SMC methods.

Impact of brain tissue filtering on neurostimulation fields: a modeling study

Electrical neurostimulation techniques, such as deep brain stimulation (DBS) and transcranial magnetic stimulation (TMS), are increasingly used in the neurosciences, e.g., for studying brain function, and for neurotherapeutics, e.g., for treating depression, epilepsy, and Parkinson's disease. The characterization of electrical properties of brain tissue has guided our fundamental understanding and application of these methods, from electrophysiologic theory to clinical dosing-metrics. Nonetheless, prior computational models have primarily relied on ex-vivo impedance measurements. We recorded the in-vivo impedances of brain tissues during neurosurgical procedures and used these results to construct MRI guided computational models of TMS and DBS neurostimulatory fields and conductance-based models of neurons exposed to stimulation. We demonstrated that tissues carry neurostimulation currents through frequency dependent resistive and capacitive properties not typically accounted for by past neurostimulation modeling work. We show that these fundamental brain tissue properties can have significant effects on the neurostimulatory-fields (capacitive and resistive current composition and spatial/temporal dynamics) and neural responses (stimulation threshold, ionic currents, and membrane dynamics). These findings highlight the importance of tissue impedance properties on neurostimulation and impact our understanding of the biological mechanisms and technological potential of neurostimulatory methods.

Characterizing learning by simultaneous analysis of continuous and binary measures of performance

Continuous observations, such as reaction and run times, and binary observations, such as correct/incorrect responses, are recorded routinely in behavioral learning experiments. Although both types of performance measures are often recorded simultaneously, the two have not been used in combination to evaluate learning. We present a state-space model of learning in which the observation process has simultaneously recorded continuous and binary measures of performance. We use these performance measures simultaneously to estimate the model parameters and the unobserved cognitive state process by maximum likelihood using an approximate expectation maximization (EM) algorithm. We introduce the concept of a reaction-time curve and reformulate our previous definitions of the learning curve, the ideal observer curve, the learning trial and between-trial comparisons of performance in terms of the new model. We illustrate the properties of the new model in an analysis of a simulated learning experiment. In the simulated data analysis, simultaneous use of the two measures of performance provided more credible and accurate estimates of the learning than either measure analyzed separately. We also analyze two actual learning experiments in which the performance of rats and of monkeys was tracked across trials by simultaneously recorded reaction and run times and the correct and incorrect responses. In the analysis of the actual experiments, our algorithm gave a straightforward, efficient way to characterize learning by combining continuous and binary measures of performance. This analysis paradigm has implications for characterizing learning and for the more general problem of combining different data types to characterize the properties of a neural system.

The dependence of spike field coherence on expected intensity

The coherence between neural spike trains and local-field potential recordings, called spike-field coherence, is of key importance in many neuroscience studies. In this work, aside from questions of estimator performance, we demonstrate that theoretical spike-field coherence for a broad class of spiking models depends on the expected rate of spiking. This rate dependence confounds the phase locking of spike events to field-potential oscillations with overall neuron activity and is demonstrated analytically, for a large class of stochastic models, and in simulation. Finally, the relationship between the spike-field coherence and the intensity field coherence is detailed analytically. This latter quantity is independent of neuron firing rate and, under commonly found conditions, is proportional to the probability that a neuron spikes at a specific phase of field oscillation. Hence, intensity field coherence is a rate-independent measure and a candidate on which to base the appropriate statistical inference of spike field synchrony.

Hippocampal replay of experience at real-world speeds

Representations related to past experiences play a critical role in memory and decision-making processes. The rat hippocampus expresses these types of representations during sharp-wave ripple (SWR) events, and previous work identified a minority of SWRs that contain ‘replay’ of spatial trajectories at ∼20x the movement speed of the animal. Efforts to understand replay typically make multiple assumptions about which events to examine and what sorts of representations constitute replay. We therefore lack a clear understanding of both the prevalence and the range of representational dynamics associated with replay. Here, we develop a state space model that uses a combination of movement dynamics of different speeds to capture the spatial content and time evolution of replay during SWRs. Using this model, we find that the large majority of replay events contain spatially coherent, interpretable content. Furthermore, many events progress at real-world, rather than accelerated, movement speeds, consistent with actual experiences.

Decoding Movement Trajectories Through a T-Maze Using Point Process Filters Applied to Place Field Data from Rat Hippocampal Region CA1

Firing activity from neural ensembles in rat hippocampus has been previously used to determine an animal's position in an open environment and separately to predict future behavioral decisions. However, a unified statistical procedure to combine information about position and behavior in environments with complex topological features from ensemble hippocampal activity has yet to be described. Here we present a two-stage computational framework that uses point process filters to simultaneously estimate the animal's location and predict future behavior from ensemble neural spiking activity. First, in the encoding stage, we linearized a two-dimensional T-maze, and used spline-based generalized linear models to characterize the place-field structure of different neurons. All of these neurons displayed highly specific position-dependent firing, which frequently had several peaks at multiple locations along the maze. When the rat was at the stem of the T-maze, the firing activity of several of these neurons also varied significantly as a function of the direction it would turn at the decision point, as detected by ANOVA. Second, in the decoding stage, we developed a state-space model for the animal's movement along a T-maze and used point process filters to accurately reconstruct both the location of the animal and the probability of the next decision. The filter yielded exact full posterior densities that were highly nongaussian and often multimodal. Our computational framework provides a reliable approach for characterizing and extracting information from ensembles of neurons with spatially specific context or task-dependent firing activity.