Spiking Activity Of Individual Neurons Research Articles

Neuronal Spiking Responses to Direct Electrical Microstimulation in the Human Cortex.

Microstimulation can modulate the activity of individual neurons to affect behavior, but the effects of stimulation on neuronal spiking are complex and remain poorly understood. This is especially challenging in the human brain where the response properties of individual neurons are sparse and heterogeneous. Here we use microelectrode arrays in the human anterior temporal lobe in 6 participants (3 female) to examine the spiking responses of individual neurons to microstimulation delivered through multiple distinct stimulation sites. We demonstrate that individual neurons can be driven with excitation or inhibition using different stimulation sites, which suggests an approach for providing direct control of spiking activity at the single-neuron level. Spiking responses are inhibitory in neurons that are close to the site of stimulation, while excitatory responses are more spatially distributed. Together, our data demonstrate that spiking responses of individual neurons can be reliably identified and manipulated in the human cortex.SIGNIFICANCE STATEMENT One of the major limitations in our ability to interface directly with the human brain is that the effects of stimulation on the activity of individual neurons remain poorly understood. This study examines the spiking responses of neurons in the human temporal cortex in response to pulses of microstimulation. This study finds that individual neurons can either be excited or inhibited depending on the site of stimulation. These data suggest an approach for modulating the spiking activity of individual neurons in the human brain.

Predicting in vitro single-neuron firing rates upon pharmacological perturbation using Graph Neural Networks.

Modern Graph Neural Networks (GNNs) provide opportunities to study the determinants underlying the complex activity patterns of biological neuronal networks. In this study, we applied GNNs to a large-scale electrophysiological dataset of rodent primary neuronal networks obtained by means of high-density microelectrode arrays (HD-MEAs). HD-MEAs allow for long-term recording of extracellular spiking activity of individual neurons and networks and enable the extraction of physiologically relevant features at the single-neuron and population level. We employed established GNNs to generate a combined representation of single-neuron and connectivity features obtained from HD-MEA data, with the ultimate goal of predicting changes in single-neuron firing rate induced by a pharmacological perturbation. The aim of the main prediction task was to assess whether single-neuron and functional connectivity features, inferred under baseline conditions, were informative for predicting changes in neuronal activity in response to a perturbation with Bicuculline, a GABA A receptor antagonist. Our results suggest that the joint representation of node features and functional connectivity, extracted from a baseline recording, was informative for predicting firing rate changes of individual neurons after the perturbation. Specifically, our implementation of a GNN model with inductive learning capability (GraphSAGE) outperformed other prediction models that relied only on single-neuron features. We tested the generalizability of the results on two additional datasets of HD-MEA recordings-a second dataset with cultures perturbed with Bicuculline and a dataset perturbed with the GABA A receptor antagonist Gabazine. GraphSAGE models showed improved prediction accuracy over other prediction models. Our results demonstrate the added value of taking into account the functional connectivity between neurons and the potential of GNNs to study complex interactions between neurons.

Parameter Estimation in Multiple Dynamic Synaptic Coupling Model Using Bayesian Point Process State-Space Modeling Framework.

It is of great interest to characterize the spiking activity of individual neurons in a cell ensemble. Many different mechanisms, such as synaptic coupling and the spiking activity of itself and its neighbors, drive a cell's firing properties. Though this is a widely studied modeling problem, there is still room to develop modeling solutions by simplifications embedded in previous models. The first shortcut is that synaptic coupling mechanisms in previous models do not replicate the complex dynamics of the synaptic response. The second is that the number of synaptic connections in these models is an order of magnitude smaller than in an actual neuron. In this research, we push this barrier by incorporating a more accurate model of the synapse and propose a system identification solution that can scale to a network incorporating hundreds of synaptic connections. Although a neuron has hundreds of synaptic connections, only a subset of these connections significantly contributes to its spiking activity. As a result, we assume the synaptic connections are sparse, and to characterize these dynamics, we propose a Bayesian point-process state-space model that lets us incorporate the sparsity of synaptic connections within the regularization technique into our framework. We develop an extended expectation-maximization. algorithm to estimate the free parameters of the proposed model and demonstrate the application of this methodology to the problem of estimating the parameters of many dynamic synaptic connections. We then go through a simulation example consisting of the dynamic synapses across a range of parameter values and show that the model parameters can be estimated using our method. We also show the application of the proposed algorithm in the intracellular data that contains 96 presynaptic connections and assess the estimation accuracy of our method using a combination of goodness-of-fit measures.

Spatial Information Based OSort for Real-Time Spike Sorting Using FPGA

Spiking activity of individual neurons can be separated from the acquired multi-unit activity with spike sorting methods. Processing the recorded high-dimensional neural data can take a large amount of time when performed on general-purpose computers. In this paper, an FPGA-based real-time spike sorting system is presented which takes into account the spatial correlation between the electrical signals recorded with closely-packed recording sites to cluster multi-channel neural data. The system uses a spatial window-based version of the Online Sorting algorithm, which uses unsupervised template-matching for clustering. The test results show that the proposed system can reach an average accuracy of 86% using simulated data (16-32 neurons, 4-10dB Signal-to-Noise Ratio), while the single-channel clustering version achieves only 74% average accuracy in the same cases on a 128-channel electrode array. The developed system was also tested on in vivo cortical recordings obtained from an anesthetized rat. The proposed FPGA-based spike sorting system can process more than 11000 spikes/second, so it can be used during in vivo experiments providing real-time feedback on the location and electrophysiological properties of well-separable single units. The proposed spike sorting system could be used to reduce the positioning error of the closely-packed recording site during a neural measurement.

Microelectrophysiological Studies of the Ratio of Excitatory to Inhibitory Synaptic Processes in the Corticonigral Projection in a Model of Parkinson’s

Experiments on 23 white male rats (250 g) analyzed the spike activity of individual neurons in the substantia nigra pars compacta (SNc, 242 neurons, n = 11) and substantia nigra pars reticulata (SNr, 289 neurons, n = 12) during high-frequency stimulation of the primary motor cortex (M1)in normal animals and in animals with a rotenone model of Parkinson’s disease (BP). SNc neurons in the model of PD showed a complete absence of depressor effects induced by stimulation, though tetanic potentiation was accompanied by posttetanic potentiation and depression at levels 1.65 and 2.02 times greater than in normal animals. In SNr neurons in normal animals, tetanic potentiation, accompanied by post-tetanic potentiation and depression, was 2.37 times greater than tetanic depression, while in the model of PD the levels of both depressor and excitatory activity induced by stimulation were below normal. Spike activity frequency in SNc and SNr neurons preceding and accompanying stimulation was significantly greater than normal in the model of PD. This is evidence for excitotoxicity accompanying neurodegenerative damage, which is completed by neuron apoptosis and death. In SNr neurons, both depressor and excitatory reactions accompanying stimulation were markedly dominant over those in SNc neurons, which is evidence for more extensive cortical projections to the SNr. Furthermore, SNc neurons demonstrated greater susceptibility to pathological changes due to poststimulus depressor effects than SNr neurons, with formation of more marked excitatory effects, which is evidence of a greater involvement of the SNc in PD. In the model of PD, lacking stimulation-induced depressor effects and more marked excitatory effects in SNc neurons, SNr neurons retained their depressor reactions and relatively decreased excitatory reactions, which is evidence of a lower level of susceptibility of SNr neurons to excitotoxicity, extreme increases in the excitability of surviving neurons compensating for the lack of excitation of dead cells.

Spike-field Granger causality for hybrid neural data analysis.

Neurotechnological innovations allow for simultaneous recording at various scales, ranging from spiking activity of individual neurons to large neural populations' local field potentials (LFPs). This capability necessitates developing multiscale analysis of spike-field activity. A joint analysis of the hybrid neural data is crucial for bridging the scales between single neurons and local networks. Granger causality is a fundamental measure to evaluate directional influences among neural signals. However, it is mainly limited to inferring causal influence between the same type of signals-either LFPs or spike trains-and not well developed between two different signal types. Here we propose a model-free, nonparametric spike-field Granger causality measure for hybrid data analysis. Our measure is distinct from existing methods in that we use "binless" spikes (precise spike timing) rather than "binned" spikes (spike counts within small consecutive time windows). The latter clearly distort the information in the mixed analysis of spikes and LFP. Therefore, our spectral estimate of spike trains is directly applied to the neural point process itself, i.e., sequences of spike times rather than spike counts. Our measure is validated by an extensive set of simulated data. When the measure is applied to LFPs and spiking activity simultaneously recorded from visual areas V1 and V4 of monkeys performing a contour detection task, we are able to confirm computationally the long-standing experimental finding of the input-output relationship between LFPs and spikes. Importantly, we demonstrate that spike-field Granger causality can be used to reveal the modulatory effects that are inaccessible by traditional methods, such that spike→LFP Granger causality is modulated by the behavioral task, whereas LFP→spike Granger causality is mainly related to the average synaptic input.NEW & NOTEWORTHY It is a pressing question to study the directional interactions between local field potential (LFP) and spiking activity. In this report, we propose a model-free, nonparametric spike-field Granger causality measure that can be used to reveal directional influences between spikes and LFPs. This new measure is crucial for bridging the scales between single neurons and neural networks; hence it represents an important step to explicate how the brain orchestrates information processing.

Rich-club connectivity, diverse population coupling, and dynamical activity patterns emerging from local cortical circuits.

Experimental studies have begun revealing essential properties of the structural connectivity and the spatiotemporal activity dynamics of cortical circuits. To integrate these properties from anatomy and physiology, and to elucidate the links between them, we develop a novel cortical circuit model that captures a range of realistic features of synaptic connectivity. We show that the model accounts for the emergence of higher-order connectivity structures, including highly connected hub neurons that form an interconnected rich-club. The circuit model exhibits a rich repertoire of dynamical activity states, ranging from asynchronous to localized and global propagating wave states. We find that around the transition between asynchronous and localized propagating wave states, our model quantitatively reproduces a variety of major empirical findings regarding neural spatiotemporal dynamics, which otherwise remain disjointed in existing studies. These dynamics include diverse coupling (correlation) between spiking activity of individual neurons and the population, dynamical wave patterns with variable speeds and precise temporal structures of neural spikes. We further illustrate how these neural dynamics are related to the connectivity properties by analysing structural contributions to variable spiking dynamics and by showing that the rich-club structure is related to the diverse population coupling. These findings establish an integrated account of structural connectivity and activity dynamics of local cortical circuits, and provide new insights into understanding their working mechanisms.

From dynamics to links: a sparse reconstruction of the topology of a neural network

Abstract One major challenge in neuroscience is the identification of interrelations between signals reflecting neural activity and how information processing occurs in the neural circuits. At the cellular and molecular level, mechanisms of signal transduction have been studied intensively and a better knowledge and understanding of some basic processes of information handling by neurons has been achieved. In contrast, little is known about the organization and function of complex neuronal networks. Experimental methods are now available to simultaneously monitor electrical activity of a large number of neurons in real time. Then, the qualitative and quantitative analysis of the spiking activity of individual neurons is a very valuable tool for the study of the dynamics and architecture of the neural networks. Such activity is not due to the sole intrinsic properties of the individual neural cells but it is mostly the consequence of the direct influence of other neurons. The deduction of the effective connectivity between neurons, whose experimental spike trains are observed, is of crucial importance in neuroscience: first for the correct interpretation of the electro-physiological activity of the involved neurons and neural networks, and, for correctly relating the electrophysiological activity to the functional tasks accomplished by the network. In this work, we propose a novel method for the identification of connectivity of neural networks using recorded voltages. Our approach is based on the assumption that the network has a topology with sparse connections. After a brief description of our method, we will report the performances and compare it to the cross-correlation computed on the spike trains, which represents a gold standard method in the field.

Persistent neural activity in auditory cortex is related to auditory working memory in humans and nonhuman primates.

Working memory is the cognitive capacity of short-term storage of information for goal-directed behaviors. Where and how this capacity is implemented in the brain are unresolved questions. We show that auditory cortex stores information by persistent changes of neural activity. We separated activity related to working memory from activity related to other mental processes by having humans and monkeys perform different tasks with varying working memory demands on the same sound sequences. Working memory was reflected in the spiking activity of individual neurons in auditory cortex and in the activity of neuronal populations, that is, in local field potentials and magnetic fields. Our results provide direct support for the idea that temporary storage of information recruits the same brain areas that also process the information. Because similar activity was observed in the two species, the cellular bases of some auditory working memory processes in humans can be studied in monkeys.

Unbiased and robust quantification of synchronization between spikes and local field potential

BackgroundIn neuroscience, relating the spiking activity of individual neurons to the local field potential (LFP) of neural ensembles is an increasingly useful approach for studying rhythmic neuronal synchronization. Many methods have been proposed to measure the strength of the association between spikes and rhythms in the LFP recordings, and most existing measures are dependent upon the total number of spikes. New methodIn the present work, we introduce a robust approach for quantifying spike–LFP synchronization which performs reliably for limited samples of data. The measure is termed as spike-triggered correlation matrix synchronization (SCMS), which takes LFP segments centered on each spike as multi-channel signals and calculates the index of spike–LFP synchronization by constructing a correlation matrix. ResultsThe simulation based on artificial data shows that the SCMS output almost does not change with the sample size. This property is of crucial importance when making comparisons between different experimental conditions. When applied to actual neuronal data recorded from the monkey primary visual cortex, it is found that the spike–LFP synchronization strength shows orientation selectivity to drifting gratings. Comparison with existing methodsIn comparison to another unbiased method, pairwise phase consistency (PPC), the proposed SCMS behaves better for noisy spike trains by means of numerical simulations. ConclusionsThis study demonstrates the basic idea and calculating process of the SCMS method. Considering its unbiasedness and robustness, the measure is of great advantage to characterize the synchronization between spike trains and rhythms present in LFP.

Copula regression analysis of simultaneously recorded frontal eye field and inferotemporal spiking activity during object-based working memory.

Inferotemporal (IT) neurons are known to exhibit persistent, stimulus-selective activity during the delay period of object-based working memory tasks. Frontal eye field (FEF) neurons show robust, spatially selective delay period activity during memory-guided saccade tasks. We present a copula regression paradigm to examine neural interaction of these two types of signals between areas IT and FEF of the monkey during a working memory task. This paradigm is based on copula models that can account for both marginal distribution over spiking activity of individual neurons within each area and joint distribution over ensemble activity of neurons between areas. Considering the popular GLMs as marginal models, we developed a general and flexible likelihood framework that uses the copula to integrate separate GLMs into a joint regression analysis. Such joint analysis essentially leads to a multivariate analog of the marginal GLM theory and hence efficient model estimation. In addition, we show that Granger causality between spike trains can be readily assessed via the likelihood ratio statistic. The performance of this method is validated by extensive simulations, and compared favorably to the widely used GLMs. When applied to spiking activity of simultaneously recorded FEF and IT neurons during working memory task, we observed significant Granger causality influence from FEF to IT, but not in the opposite direction, suggesting the role of the FEF in the selection and retention of visual information during working memory. The copula model has the potential to provide unique neurophysiological insights about network properties of the brain.

Sparse sampling: theory, methods and an application in neuroscience.

The current methods used to convert analogue signals into discrete-time sequences have been deeply influenced by the classical Shannon–Whittaker–Kotelnikov sampling theorem. This approach restricts the class of signals that can be sampled and perfectly reconstructed to bandlimited signals. During the last few years, a new framework has emerged that overcomes these limitations and extends sampling theory to a broader class of signals named signals with finite rate of innovation (FRI). Instead of characterising a signal by its frequency content, FRI theory describes it in terms of the innovation parameters per unit of time. Bandlimited signals are thus a subset of this more general definition. In this paper, we provide an overview of this new framework and present the tools required to apply this theory in neuroscience. Specifically, we show how to monitor and infer the spiking activity of individual neurons from two-photon imaging of calcium signals. In this scenario, the problem is reduced to reconstructing a stream of decaying exponentials.

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.

Task-Dependent Changes in Cross-Level Coupling between Single Neurons and Oscillatory Activity in Multiscale Networks

Understanding the principles governing the dynamic coordination of functional brain networks remains an important unmet goal within neuroscience. How do distributed ensembles of neurons transiently coordinate their activity across a variety of spatial and temporal scales? While a complete mechanistic account of this process remains elusive, evidence suggests that neuronal oscillations may play a key role in this process, with different rhythms influencing both local computation and long-range communication. To investigate this question, we recorded multiple single unit and local field potential (LFP) activity from microelectrode arrays implanted bilaterally in macaque motor areas. Monkeys performed a delayed center-out reach task either manually using their natural arm (Manual Control, MC) or under direct neural control through a brain-machine interface (Brain Control, BC). In accord with prior work, we found that the spiking activity of individual neurons is coupled to multiple aspects of the ongoing motor beta rhythm (10–45 Hz) during both MC and BC, with neurons exhibiting a diversity of coupling preferences. However, here we show that for identified single neurons, this beta-to-rate mapping can change in a reversible and task-dependent way. For example, as beta power increases, a given neuron may increase spiking during MC but decrease spiking during BC, or exhibit a reversible shift in the preferred phase of firing. The within-task stability of coupling, combined with the reversible cross-task changes in coupling, suggest that task-dependent changes in the beta-to-rate mapping play a role in the transient functional reorganization of neural ensembles. We characterize the range of task-dependent changes in the mapping from beta amplitude, phase, and inter-hemispheric phase differences to the spike rates of an ensemble of simultaneously-recorded neurons, and discuss the potential implications that dynamic remapping from oscillatory activity to spike rate and timing may hold for models of computation and communication in distributed functional brain networks.

Automatic onlinespike sorting with singular value decomposition and fuzzy C-mean clustering

BackgroundUnderstanding how neurons contribute to perception, motor functions and cognition requires the reliable detection of spiking activity of individual neurons during a number of different experimental conditions. An important problem in computational neuroscience is thus to develop algorithms to automatically detect and sort the spiking activity of individual neurons from extracellular recordings. While many algorithms for spike sorting exist, the problem of accurate and fast online sorting still remains a challenging issue.ResultsHere we present a novel software tool, called FSPS (Fuzzy SPike Sorting), which is designed to optimize: (i) fast and accurate detection, (ii) offline sorting and (iii) online classification of neuronal spikes with very limited or null human intervention. The method is based on a combination of Singular Value Decomposition for fast and highly accurate pre-processing of spike shapes, unsupervised Fuzzy C-mean, high-resolution alignment of extracted spike waveforms, optimal selection of the number of features to retain, automatic identification the number of clusters, and quantitative quality assessment of resulting clusters independent on their size. After being trained on a short testing data stream, the method can reliably perform supervised online classification and monitoring of single neuron activity. The generalized procedure has been implemented in our FSPS spike sorting software (available free for non-commercial academic applications at the address: http://www.spikesorting.com) using LabVIEW (National Instruments, USA). We evaluated the performance of our algorithm both on benchmark simulated datasets with different levels of background noise and on real extracellular recordings from premotor cortex of Macaque monkeys. The results of these tests showed an excellent accuracy in discriminating low-amplitude and overlapping spikes under strong background noise. The performance of our method is competitive with respect to other robust spike sorting algorithms.ConclusionsThis new software provides neuroscience laboratories with a new tool for fast and robust online classification of single neuron activity. This feature could become crucial in situations when online spike detection from multiple electrodes is paramount, such as in human clinical recordings or in brain-computer interfaces.

Adaptive Changes in Neuronal Synchronization in Macaque V4

A fundamental property of cortical neurons is the capacity to exhibit adaptive changes or plasticity. Whether adaptive changes in cortical responses are accompanied by changes in synchrony between individual neurons and local population activity in sensory cortex is unclear. This issue is important as synchronized neural activity is hypothesized to play an important role in propagating information in neuronal circuits. Here, we show that rapid adaptation (300 ms) to a stimulus of fixed orientation modulates the strength of oscillatory neuronal synchronization in macaque visual cortex (area V4) and influences the ability of neurons to distinguish small changes in stimulus orientation. Specifically, rapid adaptation increases the synchronization of individual neuronal responses with local population activity in the gamma frequency band (30-80 Hz). In contrast to previous reports that gamma synchronization is associated with an increase in firing rates in V4, we found that the postadaptation increase in gamma synchronization is associated with a decrease in neuronal responses. The increase in gamma-band synchronization after adaptation is functionally significant as it is correlated with an improvement in neuronal orientation discrimination performance. Thus, adaptive synchronization between the spiking activity of individual neurons and their local population can enhance temporally insensitive, rate-based-coding schemes for sensory discrimination.

Adaptation-induced synchronization in laminar cortical circuits

A fundamental feature of information processing in neocortex is the ability of individual neurons to adapt to changes in incoming stimuli. It is increasingly being understood that cortical adaptation is a phenomenon that requires network interactions. The fact that the structure of local networks depends critically on cortical layer raises the possibility that adaptation could induce specific effects in different layers. Here we show that brief exposure (300 ms) to a stimulus of fixed orientation modulates the strength of synchronization between individual neurons and local population activity in the gamma-band frequency (30-80 Hz) in macaque primary visual cortex (V1) and influences the ability of individual neurons to encode stimulus orientation. Using laminar probes, we found that although stimulus presentation elicits a large increase in the gamma synchronization of rhythmic neuronal activity in the input (granular) layers of V1, adaptation caused a pronounced increase in synchronization in the cortical output (supragranular) layers. The increase in gamma synchronization after adaptation was significantly correlated with an improvement in neuronal orientation discrimination performance only in the supragranular layers. Thus, synchronization between the spiking activity of individual neurons and their local population may enhance sensory coding to optimize network processing across laminar circuits.

Interactions and Firing Patterns of Rabbit Amygdalar Neurons in Unconditioned Fear Manifested as Freezing

The characteristics of the operation of the amygdalar neural network in unconditioned fear were identified by studying the interactions and nature of the spike activity of individual neurons in the basal and central nuclei of the amygdala in rabbits during freezing (fear), active unconditioned motor reactions (absence of fear), and in response to emotionally significant stimuli, as well as in calm waking. When rabbits froze, there were specific changes in the interactions of close-lying amygdalar neurons as compared with other states; these changes were not seen in the spike activity of individual neurons. Freezing increased the numbers of short-latency (up to 100 msec) excitatory connections and decreased the numbers of long-latency (250-450 msec) inhibitory connections. Neuron interactions were seen at frequencies in the delta2 range (2-4 Hz) more often in this state than in others. When the animals made active motor responses to the stimulus, there were decreases in the numbers of interacting neurons and increases in the numbers of longlatency (200-250 msec) excitatory and short-latency (50-200 msec) inhibitory connections, and a greater proportion of interactions occurred at frequencies of the theta1 range more frequently than in other states. Thus, the balance between the excitatory and inhibitory components of the amygdalar neural network is important for the occurrence of fear. Language: en

Neural basis of visual selective attention.

Attentional modulation along the object-recognition pathway of the cortical visual system of primates has been shown to consist of enhanced representation of the retinal input at a specific location in space, or of objects located anywhere in the visual field which possess a critical object feature. Moreover, selective attention mechanisms allow the visual system to resolve competition among multiple objects in a crowded scene in favor of the object that is relevant for the current behavior. Finally, selective attention affects the spontaneous activity of neurons as well as their visually driven responses, and it does so not only by modulating the spiking activity of individual neurons, but also by modulating the degree of coherent firing within the critical neuronal populations. WIREs Cogni Sci 2011 2 392-407 DOI: 10.1002/wcs.117 For further resources related to this article, please visit the WIREs website.

A comparison of spike time prediction and receptive field mapping with point process generalized linear models, Wiener-Voltera kernels, and spike-triggered averaging methods

Characterizing neuronal encoding is essential for understanding information processing in the brain. Three methods are commonly used to characterize the relationship between neural spiking activity and the features of putative stimuli. These methods include: Wiener-Volterra kernel methods (WVK), the spike-triggered average (STA), and more recently, the point process generalized linear model (GLM). We compared the performance of these three approaches in estimating receptive field properties and orientation tuning of 251 V1 neurons recorded from 2 monkeys during a fixation period in response to a moving bar. The GLM consisted of two formulations of the conditional intensity function for a point process characterization of the spiking activity: one with a stimulus only component and one with the stimulus and spike history. We fit the GLMs by maximum likelihood using GLMfit in Matlab. Goodness-of-fit was assessed using cross-validation with Kolmogorov-Smirnov (KS) tests based on the time-rescaling theorem to evaluate the accuracy with which each model predicts the spiking activity of individual neurons and for each movement direction (4016 models in total, for 251 neurons and 16 different directions).