Sequences Of Action Potentials Research Articles

Steady-state approximations for Hodgkin-Huxley cell models: Reduction of order for uterine smooth muscle cell model.

Multi-scale mathematical bioelectrical models of organs such as the uterus, stomach or heart present challenges both for accuracy and computational tractability. These multi-scale models are typically founded on models of biological cells derived from the classic Hodkgin-Huxley (HH) formalism. Ion channel behaviour is tracked with dynamical variables representing activation or inactivation of currents that relax to steady-state dependencies on cellular membrane voltage. Timescales for relaxation may be orders of magnitude faster than companion ion channel variables or phenomena of physiological interest for the entire cell (such as bursting sequences of action potentials) or the entire organ (such as electromechanical coordination). Exploiting these time scales with steady-state approximations for relatively fast-acting systems is a well-known but often overlooked approach as evidenced by recent published models. We thus investigate feasibility of an extensive reduction of order for an HH-type cell model with steady-state approximations to the full dynamical activation and inactivation ion channel variables. Our effort utilises a published comprehensive uterine smooth muscle cell model that encompasses 19 ordinary differential equations and 105 formulations overall. The numerous ion channel submodels in the published model exhibit relaxation times ranging from order 10-1 to 105 milliseconds. Substitution of the faster dynamic variables with steady-state formulations demonstrates both an accurate reproduction of the full model and substantial improvements in time-to-solve, for test cases performed. Our demonstration here of an effective and relatively straightforward reduction method underlines the particular importance of considering time scales for model simplification before embarking on large-scale computations or parameter sweeps. As a preliminary complement to more intensive reduction of order methods such as parameter sensitivity and bifurcation analysis, this approach can rapidly and accurately improve computational tractability for challenging multi-scale organ modelling efforts.

A biophysical and statistical modeling paradigm for connecting neural physiology and function

To understand single neuron computation, it is necessary to know how specific physiological parameters affect neural spiking patterns that emerge in response to specific stimuli. Here we present a computational pipeline combining biophysical and statistical models that provides a link between variation in functional ion channel expression and changes in single neuron stimulus encoding. More specifically, we create a mapping from biophysical model parameters to stimulus encoding statistical model parameters. Biophysical models provide mechanistic insight, whereas statistical models can identify associations between spiking patterns and the stimuli they encode. We used public biophysical models of two morphologically and functionally distinct projection neuron cell types: mitral cells (MCs) of the main olfactory bulb, and layer V cortical pyramidal cells (PCs). We first simulated sequences of action potentials according to certain stimuli while scaling individual ion channel conductances. We then fitted point process generalized linear models (PP-GLMs), and we constructed a mapping between the parameters in the two types of models. This framework lets us detect effects on stimulus encoding of changing an ion channel conductance. The computational pipeline combines models across scales and can be applied as a screen of channels, in any cell type of interest, to identify ways that channel properties influence single neuron computation.

Signal Fluctuations and the Information Transmission Rates in Binary Communication Channels.

In the nervous system, information is conveyed by sequence of action potentials, called spikes-trains. As MacKay and McCulloch suggested, spike-trains can be represented as bits sequences coming from Information Sources (). Previously, we studied relations between spikes’ Information Transmission Rates and their correlations, and frequencies. Now, I concentrate on the problem of how spikes fluctuations affect . The are typically modeled as stationary stochastic processes, which I consider here as two-state Markov processes. As a spike-trains’ fluctuation measure, I assume the standard deviation , which measures the average fluctuation of spikes around the average spike frequency. I found that the character of and signal fluctuations relation strongly depends on the parameter s being a sum of transitions probabilities from a no spike state to spike state. The estimate of the Information Transmission Rate was found by expressions depending on the values of signal fluctuations and parameter s. It turned out that for smaller , the quotient has a maximum and can tend to zero depending on transition probabilities, while for , the is separated from 0. Additionally, it was also shown that quotient by variance behaves in a completely different way. Similar behavior was observed when classical Shannon entropy terms in the Markov entropy formula are replaced by their approximation with polynomials. My results suggest that in a noisier environment , to get appropriate reliability and efficiency of transmission, with higher tendency of transition from the no spike to spike state should be applied. Such selection of appropriate parameters plays an important role in designing learning mechanisms to obtain networks with higher performance.

Broad-Spectrum Antiviral Entry Inhibition by Interfacially Active Peptides.

Numerous peptides inhibit the entry of enveloped viruses into cells. Some of these peptides have been shown to inhibit multiple unrelated viruses. We have suggested that such broad-spectrum antiviral peptides share a property called interfacial activity; they are somewhat hydrophobic and amphipathic, with a propensity to interact with the interfacial zones of lipid bilayer membranes. In this study, we further tested the hypothesis that such interfacial activity is a correlate of broad-spectrum antiviral activity. In this study, several families of peptides, selected for the ability to partition into and disrupt membrane integrity but with no known antiviral activity, were tested for the ability to inhibit multiple diverse enveloped viruses. These include Lassa pseudovirus, influenza virus, dengue virus type 2, herpes simplex virus 1, and nonenveloped human adenovirus 5. Various families of interfacially active peptides caused potent inhibition of all enveloped viruses tested at low and submicromolar concentrations, well below the range in which they are toxic to mammalian cells. These membrane-active peptides block uptake and fusion with the host cell by rapidly and directly interacting with virions, destabilizing the viral envelope, and driving virus aggregation and/or intervirion envelope fusion. We speculate that the molecular characteristics shared by these peptides can be exploited to enable the design, optimization, or molecular evolution of novel broad-spectrum antiviral therapeutics.IMPORTANCE New classes of antiviral drugs are needed to treat the ever-changing viral disease landscape. Current antiviral drugs treat only a small number of viral diseases, leaving many patients with established or emerging infections to be treated solely with supportive care. Recent antiviral peptide research has produced numerous membrane-interacting peptides that inhibit diverse enveloped viruses in vitro and in vivo Peptide therapeutics are becoming more common, with over 60 FDA-approved peptides for clinical use. Included in this class of therapeutics is enfuvirtide, a 36-residue peptide drug that inhibits HIV entry/fusion. Due to their broad-spectrum mechanism of action and enormous potential sequence diversity, peptides that inhibit virus entry could potentially fulfill the need for new antiviral therapeutics; however, a better understanding of their mechanism is needed for the optimization or evolution of sequence design to combat the wide landscape of viral disease.

Blind deconvolution for spike inference from fluorescence recordings

The parallel developments of genetically-encoded calcium indicators and fast fluorescence imaging techniques allows one to simultaneously record neural activity of extended neuronal populations in vivo. To fully harness the potential of functional imaging, one needs to infer the sequence of action potentials from fluorescence traces. Here we build on recently proposed computational approaches to develop a blind sparse deconvolution (BSD) algorithm based on a generative model for inferring spike trains from fluorescence traces. BSD features, (1) automatic (fully unsupervised) estimation of the hyperparameters, such as spike amplitude, noise level and rise and decay time constants, (2) a novel analytical estimate of the sparsity prior, which yields enhanced robustness and computational speed with respect to existing methods, (3) automatic thresholding for binarizing spikes that maximizes the precision-recall performance, (4) super-resolution capabilities increasing the temporal resolution beyond the fluorescence signal acquisition rate. BSD also uniquely provides theoretically-grounded estimates of the expected performance of the spike reconstruction in terms of precision-recall and temporal accuracy for each recording. The performance of the algorithm is established using synthetic data and through the SpikeFinder challenge, a community-based initiative for spike-rate inference benchmarking based on a collection of joint electrophysiological and fluorescence recordings. Our method outperforms classical sparse deconvolution algorithms in terms of robustness, speed and/or accuracy and performs competitively in the SpikeFinder challenge. This algorithm is modular, easy-to-use and made freely available. Its novel features can thus be incorporated in a straightforward way into existing calcium imaging packages.

Neuronal Transmission of Subthreshold Periodic Stimuli Via Symbolic Spike Patterns.

We study how sensory neurons detect and transmit a weak external stimulus. We use the FitzHugh–Nagumo model to simulate the neuronal activity. We consider a sub-threshold stimulus, i.e., the stimulus is below the threshold needed for triggering action potentials (spikes). However, in the presence of noise the neuron that perceives the stimulus fires a sequence of action potentials (a spike train) that carries the stimulus’ information. To yield light on how the stimulus’ information can be encoded and transmitted, we consider the simplest case of two coupled neurons, such that one neuron (referred to as neuron 1) perceives a subthreshold periodic signal but the second neuron (neuron 2) does not perceive the signal. We show that, for appropriate coupling and noise strengths, both neurons fire spike trains that have symbolic patterns (defined by the temporal structure of the inter-spike intervals), whose frequencies of occurrence depend on the signal’s amplitude and period, and are similar for both neurons. In this way, the signal information encoded in the spike train of neuron 1 propagates to the spike train of neuron 2. Our results suggest that sensory neurons can exploit the presence of neural noise to fire spike trains where the information of a subthreshold stimulus is encoded in over expressed and/or in less expressed symbolic patterns.

Precise timing is ubiquitous, consistent, and coordinated across a comprehensive, spike-resolved flight motor program

Sequences of action potentials, or spikes, carry information in the number of spikes and their timing. Spike timing codes are critical in many sensory systems, but there is now growing evidence that millisecond-scale changes in timing also carry information in motor brain regions, descending decision-making circuits, and individual motor units. Across all of the many signals that control a behavior, how ubiquitous, consistent, and coordinated are spike timing codes? Assessing these open questions ideally involves recording across the whole motor program with spike-level resolution. To do this, we took advantage of the relatively few motor units controlling the wings of a hawk moth, Manduca sexta. We simultaneously recorded nearly every action potential from all major wing muscles and the resulting forces in tethered flight. We found that timing encodes more information about turning behavior than spike count in every motor unit, even though there is sufficient variation in count alone. Flight muscles vary broadly in function as well as in the number and timing of spikes. Nonetheless, each muscle with multiple spikes consistently blends spike timing and count information in a 3:1 ratio. Coding strategies are consistent. Finally, we assess the coordination of muscles using pairwise redundancy measured through interaction information. Surprisingly, not only are all muscle pairs coordinated, but all coordination is accomplished almost exclusively through spike timing, not spike count. Spike timing codes are ubiquitous, consistent, and essential for coordination.

Restitution and Stability of Human Ventricular Action Potential at High and Variable Pacing Rate

Despite the key role of beat-to-beat action potential (AP) variability in the onset of ventricular arrhythmias at high pacing rate, the knowledge of the involved dynamics and of effective prognostic parameters is largely incomplete. Electrical restitution (ER), the way AP duration (APD) senses changes in preceding cycle length (CL), has been used to monitor transition to arrhythmias. The use of standard ER (sER), though, is controversial, not always suitable for in vivo and only rarely for clinical applications. By means of simulations on a human ventricular AP model, I investigate the dynamics of APD at high pacing rate under sinusoidally, saw-tooth, and randomly variable pacing CLs. AP sequences were compared in terms of beat-to-beat restitution (btb-ER) and of the collections of sER curves generated from each beat. A definition of APD stability is also proposed, based on successive APD changes introduced in an AP sequence by a premature beat. The explored CL range includes values leading to APD alternans under constant pacing. Three different types of response to CL variability were found, corresponding to progressively higher rate of beat-to-beat CL changes. Low rates (∼1 ms/beat) generate a btb-ER dominated by steady-state rate dependence of APD (type 1), intermediate rates (∼5 ms/beat) lead to a btb-ER similar to a single sER (type 2), and high rates (∼20 ms/beat) to hysteretic btb-ER under periodic pacing and to a vertically spread btb-ER in the case of random pacing (type 3). Stability of AP repolarization always increases with the rate of CL changes. Thus, rather than looking at sER slope, which requires additional interventions during the recording of cardiac electrical activity, this study provides rationale for the use of btb-ER representations as predictors of repolarization stability under extreme pacing conditions, known to be critical for the arrhythmia development.

Synaptic Information Transmission in a Two-State Model of Short-Term Facilitation

Action potentials (spikes) can trigger the release of a neurotransmitter at chemical synapses between neurons. Such release is uncertain, as it occurs only with a certain probability. Moreover, synaptic release can occur independently of an action potential (asynchronous release) and depends on the history of synaptic activity. We focus here on short-term synaptic facilitation, in which a sequence of action potentials can temporarily increase the release probability of the synapse. In contrast to the phenomenon of short-term depression, quantifying the information transmission in facilitating synapses remains to be done. We find rigorous lower and upper bounds for the rate of information transmission in a model of synaptic facilitation. We treat the synapse as a two-state binary asymmetric channel, in which the arrival of an action potential shifts the synapse to a facilitated state, while in the absence of a spike, the synapse returns to its baseline state. The information bounds are functions of both the asynchronous and synchronous release parameters. If synchronous release facilitates more than asynchronous release, the mutual information rate increases. In contrast, short-term facilitation degrades information transmission when the synchronous release probability is intrinsically high. As synaptic release is energetically expensive, we exploit the information bounds to determine the energy–information trade-off in facilitating synapses. We show that unlike information rate, the energy-normalized information rate is robust with respect to variations in the strength of facilitation.

Simulation Study of Intermittent Axonal Block and Desynchronization Effect Induced by High-Frequency Stimulation of Electrical Pulses.

Deep brain stimulation (DBS) has been successfully used in treating neural disorders in brain, such as Parkinson’s disease and epilepsy. However, the precise mechanisms of DBS remain unclear. Regular DBS therapy utilizes high-frequency stimulation (HFS) of electrical pulses. Among all of neuronal elements, axons are mostly inclined to be activated by electrical pulses. Therefore, the response of axons may play an important role in DBS treatment. To study the axonal responses during HFS, we developed a computational model of myelinated axon to simulate sequences of action potentials generated in single and multiple axons (an axon bundle) by stimulations. The stimulations are applied extracellularly by a point source of current pulses with a frequency of 50–200 Hz. Additionally, our model takes into account the accumulation of potassium ions in the peri-axonal spaces. Results show that the increase of extracellular potassium generates intermittent depolarization block in the axons during HFS. Under the state of alternate block and recovery, axons fire action potentials at a rate far lower than the frequency of stimulation pulses. In addition, the degree of axonal block is highly related to the distance between the axons and the stimulation point. The differences in the degree of block for individual axons in a bundle result in desynchronized firing among the axons. Stimulations with higher frequency and/or greater intensity can induce axonal block faster and increase the desynchronization effect on axonal firing. Presumably, the desynchronized axonal activity induced by HFS could generate asynchronous activity in the population of target neurons downstream thereby suppressing over-synchronized firing of neurons in pathological conditions. The desynchronization effect generated by intermittent activation of axons may be crucial for DBS therapy. The present study provides new insights into the mechanisms of DBS, which is significant for advancing the application of DBS.

Estimation of neuronal firing rate using Bayesian Adaptive Kernel Smoother (BAKS).

Neurons use sequences of action potentials (spikes) to convey information across neuronal networks. In neurophysiology experiments, information about external stimuli or behavioral tasks has been frequently characterized in term of neuronal firing rate. The firing rate is conventionally estimated by averaging spiking responses across multiple similar experiments (or trials). However, there exist a number of applications in neuroscience research that require firing rate to be estimated on a single trial basis. Estimating firing rate from a single trial is a challenging problem and current state-of-the-art methods do not perform well. To address this issue, we develop a new method for estimating firing rate based on a kernel smoothing technique that considers the bandwidth as a random variable with prior distribution that is adaptively updated under an empirical Bayesian framework. By carefully selecting the prior distribution together with Gaussian kernel function, an analytical expression can be achieved for the kernel bandwidth. We refer to the proposed method as Bayesian Adaptive Kernel Smoother (BAKS). We evaluate the performance of BAKS using synthetic spike train data generated by biologically plausible models: inhomogeneous Gamma (IG) and inhomogeneous inverse Gaussian (IIG). We also apply BAKS to real spike train data from non-human primate (NHP) motor and visual cortex. We benchmark the proposed method against established and previously reported methods. These include: optimized kernel smoother (OKS), variable kernel smoother (VKS), local polynomial fit (Locfit), and Bayesian adaptive regression splines (BARS). Results using both synthetic and real data demonstrate that the proposed method achieves better performance compared to competing methods. This suggests that the proposed method could be useful for understanding the encoding mechanism of neurons in cognitive-related tasks. The proposed method could also potentially improve the performance of brain-machine interface (BMI) decoder that relies on estimated firing rate as the input.

Slow gain fluctuations limit benefits of temporal integration in visual cortex.

Sensory neurons represent stimulus information with sequences of action potentials that differ across repeated measurements. This variability limits the information that can be extracted from momentary observations of a neuron's response. It is often assumed that integrating responses over time mitigates this limitation. However, temporal response correlations can reduce the benefits of temporal integration. We examined responses of individual orientation-selective neurons in the primary visual cortex of two macaque monkeys performing an orientation-discrimination task. The signal-to-noise ratio of temporally integrated responses increased for durations up to a few hundred milliseconds but saturated for longer durations. This was true even when cells exhibited little or no adaptation in their response levels. These observations are well explained by a statistical response model in which spikes arise from a Poisson process whose stimulus-dependent rate is modulated by slow, stimulus-independent fluctuations in gain. The response variability arising from the Poisson process is reduced by temporal integration, but the slow modulatory nature of variability due to gain fluctuations is not. Slow gain fluctuations therefore impose a fundamental limit on the benefits of temporal integration.

Dendritic Properties Control Energy Efficiency of Action Potentials in Cortical Pyramidal Cells.

Neural computation is performed by transforming input signals into sequences of action potentials (APs), which is metabolically expensive and limited by the energy available to the brain. The metabolic efficiency of single AP has important consequences for the computational power of the cell, which is determined by its biophysical properties and morphologies. Here we adopt biophysically-based two-compartment models to investigate how dendrites affect energy efficiency of APs in cortical pyramidal neurons. We measure the Na+ entry during the spike and examine how it is efficiently used for generating AP depolarization. We show that increasing the proportion of dendritic area or coupling conductance between two chambers decreases Na+ entry efficiency of somatic AP. Activating inward Ca2+ current in dendrites results in dendritic spike, which increases AP efficiency. Activating Ca2+-activated outward K+ current in dendrites, however, decreases Na+ entry efficiency. We demonstrate that the active and passive dendrites take effects by altering the overlap between Na+ influx and internal current flowing from soma to dendrite. We explain a fundamental link between dendritic properties and AP efficiency, which is essential to interpret how neural computation consumes metabolic energy and how biophysics and morphologies contribute to such consumption.

Emergent stochastic oscillations and signal detection in tree networks of excitable elements

We study the stochastic dynamics of strongly-coupled excitable elements on a tree network. The peripheral nodes receive independent random inputs which may induce large spiking events propagating through the branches of the tree and leading to global coherent oscillations in the network. This scenario may be relevant to action potential generation in certain sensory neurons, which possess myelinated distal dendritic tree-like arbors with excitable nodes of Ranvier at peripheral and branching nodes and exhibit noisy periodic sequences of action potentials. We focus on the spiking statistics of the central node, which fires in response to a noisy input at peripheral nodes. We show that, in the strong coupling regime, relevant to myelinated dendritic trees, the spike train statistics can be predicted from an isolated excitable element with rescaled parameters according to the network topology. Furthermore, we show that by varying the network topology the spike train statistics of the central node can be tuned to have a certain firing rate and variability, or to allow for an optimal discrimination of inputs applied at the peripheral nodes.

Motor control by precisely timed spike patterns

A fundamental problem in neuroscience is understanding how sequences of action potentials ("spikes") encode information about sensory signals and motor outputs. Although traditional theories assume that this information is conveyed by the total number of spikes fired within a specified time interval (spike rate), recent studies have shown that additional information is carried by the millisecond-scale timing patterns of action potentials (spike timing). However, it is unknown whether or how subtle differences in spike timing drive differences in perception or behavior, leaving it unclear whether the information in spike timing actually plays a role in brain function. By examining the activity of individual motor units (the muscle fibers innervated by a single motor neuron) and manipulating patterns of activation of these neurons, we provide both correlative and causal evidence that the nervous system uses millisecond-scale variations in the timing of spikes within multispike patterns to control a vertebrate behavior-namely, respiration in the Bengalese finch, a songbird. These findings suggest that a fundamental assumption of current theories of motor coding requires revision.

Kernel Analysis for Estimating the Connectivity of a Network with Event Sequences

AbstractEstimating the connectivity of a network from events observed at each node has many applications. One prominent example is found in neuroscience, where spike trains (sequences of action potentials) are observed at each neuron, but the way in which these neurons are connected is unknown. This paper introduces a novel method for estimating connections between nodes using a similarity measure between sequences of event times. Specifically, a normalized positive definite kernel defined on spike trains was used. The proposed method was evaluated using synthetic and real data, by comparing with methods using transfer entropy and the Victor-Purpura distance. Synthetic data was generated using CERM (Coupled Escape-Rate Model), a model that generates various spike trains. Real data recorded from the visual cortex of an anaesthetized cat was analyzed as well. The results showed that the proposed method provides an effective way of estimating the connectivity of a network when the time sequences of events are the only available information.

Metabolic Energy of Action Potentials Modulated by Spike Frequency Adaptation.

Spike frequency adaptation (SFA) exists in many types of neurons, which has been demonstrated to improve their abilities to process incoming information by synapses. The major carrier used by a neuron to convey synaptic signals is the sequences of action potentials (APs), which have to consume substantial metabolic energies to initiate and propagate. Here we use conductance-based models to investigate how SFA modulates the AP-related energy of neurons. The SFA is attributed to either calcium-activated K+ (IAHP) or voltage-activated K+ (IM) current. We observe that the activation of IAHP or IM increases the Na+ load used for depolarizing membrane, while produces few effects on the falling phase of AP. Then, the metabolic energy involved in Na+ current significantly increases from one AP to the next, while for K+ current it is less affected. As a consequence, the total energy cost by each AP gets larger as firing rate decays down. It is also shown that the minimum Na+ charge needed for the depolarization of each AP is unaffected during the course of SFA. This indicates that the activation of either adaptation current makes APs become less efficient to use Na+ influx for their depolarization. Further, our simulations demonstrate that the different biophysical properties of IM and IAHP result in distinct modulations of metabolic energy usage for APs. These investigations provide a fundamental link between adaptation currents and neuronal energetics, which could facilitate to interpret how SFA participates in neuronal information processing.

A minimum-error, energy-constrained neural code is an instantaneous-rate code.

Sensory neurons code information about stimuli in their sequence of action potentials (spikes). Intuitively, the spikes should represent stimuli with high fidelity. However, generating and propagating spikes is a metabolically expensive process. It is therefore likely that neural codes have been selected to balance energy expenditure against encoding error. Our recently proposed optimal, energy-constrained neural coder (Jones et al. Frontiers in Computational Neuroscience, 9, 61 2015) postulates that neurons time spikes to minimize the trade-off between stimulus reconstruction error and expended energy by adjusting the spike threshold using a simple dynamic threshold. Here, we show that this proposed coding scheme is related to existing coding schemes, such as rate and temporal codes. We derive an instantaneous rate coder and show that the spike-rate depends on the signal and its derivative. In the limit of high spike rates the spike train maximizes fidelity given an energy constraint (average spike-rate), and the predicted interspike intervals are identical to those generated by our existing optimal coding neuron. The instantaneous rate coder is shown to closely match the spike-rates recorded from P-type primary afferents in weakly electric fish. In particular, the coder is a predictor of the peristimulus time histogram (PSTH). When tested against in vitro cortical pyramidal neuron recordings, the instantaneous spike-rate approximates DC step inputs, matching both the average spike-rate and the time-to-first-spike (a simple temporal code). Overall, the instantaneous rate coder relates optimal, energy-constrained encoding to the concepts of rate-coding and temporal-coding, suggesting a possible unifying principle of neural encoding of sensory signals.

Statistical structure of neural spiking under non-Poissonian or other non-white stimulation.

Nerve cells in the brain generate sequences of action potentials with a complex statistics. Theoretical attempts to understand this statistics were largely limited to the case of a temporally uncorrelated input (Poissonian shot noise) from the neurons in the surrounding network. However, the stimulation from thousands of other neurons has various sorts of temporal structure. Firstly, input spike trains are temporally correlated because their firing rates can carry complex signals and because of cell-intrinsic properties like neural refractoriness, bursting, or adaptation. Secondly, at the connections between neurons, the synapses, usage-dependent changes in the synaptic weight (short-term plasticity) further shape the correlation structure of the effective input to the cell. From the theoretical side, it is poorly understood how these correlated stimuli, so-called colored noise, affect the spike train statistics. In particular, no standard method exists to solve the associated first-passage-time problem for the interspike-interval statistics with an arbitrarily colored noise. Assuming that input fluctuations are weaker than the mean neuronal drive, we derive simple formulas for the essential interspike-interval statistics for a canonical model of a tonically firing neuron subjected to arbitrarily correlated input from the network. We verify our theory by numerical simulations for three paradigmatic situations that lead to input correlations: (i) rate-coded naturalistic stimuli in presynaptic spike trains; (ii) presynaptic refractoriness or bursting; (iii) synaptic short-term plasticity. In all cases, we find severe effects on interval statistics. Our results provide a framework for the interpretation of firing statistics measured in vivo in the brain.

Independent theta phase coding accounts for CA1 population sequences and enables flexible remapping.

Hippocampal place cells encode an animal's past, current, and future location through sequences of action potentials generated within each cycle of the network theta rhythm. These sequential representations have been suggested to result from temporally coordinated synaptic interactions within and between cell assemblies. Instead, we find through simulations and analysis of experimental data that rate and phase coding in independent neurons is sufficient to explain the organization of CA1 population activity during theta states. We show that CA1 population activity can be described as an evolving traveling wave that exhibits phase coding, rate coding, spike sequences and that generates an emergent population theta rhythm. We identify measures of global remapping and intracellular theta dynamics as critical for distinguishing mechanisms for pacemaking and coordination of sequential population activity. Our analysis suggests that, unlike synaptically coupled assemblies, independent neurons flexibly generate sequential population activity within the duration of a single theta cycle.