Oscillatory Waveform Shape and Temporal Spike Correlations Differ across Bat Frontal and Auditory Cortex.

Neural oscillations are associated with diverse computations in the mammalian brain. The waveform shape of oscillatory activity measured in the cortex relates to local physiology and can be informative about aberrant or dynamically changing states. However, how waveform shape differs across distant yet functionally and anatomically related cortical regions is largely unknown. In this study, we capitalize on simultaneous recordings of local field potentials (LFPs) in the auditory and frontal cortices of awake, male Carollia perspicillata bats to examine, on a cycle-by-cycle basis, waveform shape differences across cortical regions. We find that waveform shape differs markedly in the fronto-auditory circuit even for temporally correlated rhythmic activity in comparable frequency ranges (i.e., in the delta and gamma bands) during spontaneous activity. In addition, we report consistent differences between areas in the variability of waveform shape across individual cycles. A conceptual model predicts higher spike-spike and spike-LFP correlations in regions with more asymmetric shapes, a phenomenon that was observed in the data: spike-spike and spike-LFP correlations were higher in the frontal cortex. The model suggests a relationship between waveform shape differences and differences in spike correlations across cortical areas. Altogether, these results indicate that oscillatory activity in the frontal and auditory cortex possesses distinct dynamics related to the anatomical and functional diversity of the fronto-auditory circuit.

Noise or signal? Spontaneous activity of dorsal horn neurons: patterns and function in health and disease.

Spontaneous activity refers to the firing of action potentials by neurons in the absence of external stimulation. Initially considered an artifact or "noise" in the nervous system, it is now recognized as a potential feature of neural function. Spontaneous activity has been observed in various brain areas, in experimental preparations from different animal species, and in live animals and humans using non-invasive imaging techniques. In this review, we specifically focus on the spontaneous activity of dorsal horn neurons of the spinal cord. We use a historical perspective to set the basis for a novel classification of the different patterns of spontaneous activity exhibited by dorsal horn neurons. Then we examine the origins of this activity and propose a model circuit to explain how the activity is generated and transmitted to the dorsal horn. Finally, we discuss possible roles of this activity during development and during signal processing under physiological conditions and pain states. By analyzing recent studies on the spontaneous activity of dorsal horn neurons, we aim to shed light on its significance in sensory processing. Understanding the different patterns of activity, the origins of this activity, and the potential roles it may play, will contribute to our knowledge of sensory mechanisms, including pain, to facilitate the modeling of spinal circuits and hopefully to explore novel strategies for pain treatment.

Estimating Mutual Information for Spike Trains: A Bird Song Example.

Zebra finches are a model animal used in the study of audition. They are adept at recognizing zebra finch songs, and the neural pathway involved in song recognition is well studied. Here, this example is used to illustrate the estimation of mutual information between stimuli and responses using a Kozachenko-Leonenko estimator. The challenge in calculating mutual information for spike trains is that there are no obvious coordinates for the data. The Kozachenko-Leonenko estimator does not require coordinates; it relies only on the distance between data points. In the case of bird songs, estimating the mutual information demonstrates that the information content of spiking does not diminish as the song progresses.

Deficits in Behavioral and Neuronal Pattern Separation in Temporal Lobe Epilepsy.

In temporal lobe epilepsy, the ability of the dentate gyrus to limit excitatory cortical input to the hippocampus breaks down, leading to seizures. The dentate gyrus is also thought to help discriminate between similar memories by performing pattern separation, but whether epilepsy leads to a breakdown in this neural computation, and thus to mnemonic discrimination impairments, remains unknown. Here we show that temporal lobe epilepsy is characterized by behavioral deficits in mnemonic discrimination tasks, in both humans (females and males) and mice (C57Bl6 males, systemic low-dose kainate model). Using a recently developed assay in brain slices of the same epileptic mice, we reveal a decreased ability of the dentate gyrus to perform certain forms of pattern separation. This is because of a subset of granule cells with abnormal bursting that can develop independently of early EEG abnormalities. Overall, our results linking physiology, computation, and cognition in the same mice advance our understanding of episodic memory mechanisms and their dysfunction in epilepsy.SIGNIFICANCE STATEMENT People with temporal lobe epilepsy (TLE) often have learning and memory impairments, sometimes occurring earlier than the first seizure, but those symptoms and their biological underpinnings are poorly understood. We focused on the dentate gyrus, a brain region that is critical to avoid confusion between similar memories and is anatomically disorganized in TLE. We show that both humans and mice with TLE experience confusion between similar situations. This impairment coincides with a failure of the dentate gyrus to disambiguate similar input signals because of pathologic bursting in a subset of neurons. Our work bridges seizure-oriented and memory-oriented views of the dentate gyrus function, suggests a mechanism for cognitive symptoms in TLE, and supports a long-standing hypothesis of episodic memory theories.

Population burst propagation across interacting areas of the brain.

For many perceptual and behavioral tasks, a prominent feature of neural spike trains involves high firing rates across relatively short intervals of time. We call these events "population bursts." Because during a population burst information is, presumably, transmitted from one part of the brain to another, burst timing should reveal activity related to the flow of information across neural circuits. We developed a statistical method (based on a point process model) of determining, accurately, the time of the maximum (peak) population firing rate on a trial-by-trial basis and used it to characterize burst propagation across areas. We then examined the tendency of peak firing rates in distinct brain areas to shift earlier or later in time, together, across repeated trials, and found this trial-to-trial coupling of peak times to be a sensitive indicator of interaction across populations. In the data we examined, from the Allen Brain Observatory, we found many very strong correlations (95% confidence intervals above 0.75) in cases where standard methods were unable to demonstrate cross-area correlation. The statistical model introduced cross-area covariation only through population-level trial-dependent time shifts and gain constants (values of which were learned from the data), yet it provided very good fits to data histograms, including histograms of spike count correlations within and across visual areas. Our results demonstrate the utility of carefully assessing timing and propagation, across brain regions, of transient bursts in neural population activity, based on multiple spike train recordings.NEW & NOTEWORTHY We developed a novel statistical method for identifying coordinated propagation of activity across populations of spiking neurons, with high temporal accuracy. Using simultaneous recordings from three visual areas we document precise timing relationships on a trial-by-trial basis, and we show how previously existing techniques can fail to discover coordinated activity in cases where the new approach finds very strong cross-area correlation.

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 (IS). Previously, we studied relations between spikes' Information Transmission Rates (ITR) and their correlations, and frequencies. Now, I concentrate on the problem of how spikes fluctuations affect ITR. The IS 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 ITR 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 s<1, the quotient ITRσ has a maximum and can tend to zero depending on transition probabilities, while for 1

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.

Dynamical response properties of neocortical neurons to conductance-driven time-varying inputs.

Ensembles of cortical neurons can track fast-varying inputs and relay them in their spike trains, far beyond the cut-off imposed by membrane passive electrical properties and mean firing rates. Initially explored in silico and later demonstrated experimentally, investigating how neurons respond to sinusoidally modulated stimuli provides a deeper insight into spike initiation mechanisms and information processing than conventional F-I curve methodologies. Besides net membrane currents, physiological synaptic inputs can also induce a stimulus-dependent modulation of the total membrane conductance, which is not reproduced by standard current-clamp protocols. Here, we investigated whether rat cortical neurons can track fast temporal modulations over a noisy conductance background. We also determined input-output transfer properties over a range of conditions, including: distinct presynaptic activation rates, postsynaptic firing rates and variability and type of temporal modulations. We found a very broad signal transfer bandwidth across all conditions, similar large cut-off frequencies and power-law attenuations of fast-varying inputs. At slow and intermediate input modulations, the response gain decreased for increasing output mean firing rates. The gain also decreased significantly for increasing intensities of background synaptic activity, thus generalising earlier studies on F-I curves. We also found a direct correlation between the action potentials' onset rapidness and the neuronal bandwidth. Our novel results extend previous investigations of dynamical response properties to non-stationary and conductance-driven conditions, and provide computational neuroscientists with a novel set of observations that models must capture when aiming to replicate cortical cellular excitability.

Linear-nonlinear-time-warp-poisson models of neural activity.

Prominent models of spike trains assume only one source of variability - stochastic (Poisson) spiking - when stimuli and behavior are fixed. However, spike trains may also reflect variability due to internal processes such as planning. For example, we can plan a movement at one point in time and execute it at some arbitrary later time. Neurons involved in planning may thus share an underlying time course that is not precisely locked to the actual movement. Here we combine the standard Linear-Nonlinear-Poisson (LNP) model with Dynamic Time Warping (DTW) to account for shared temporal variability. When applied to recordings from macaque premotor cortex, we find that time warping considerably improves predictions of neural activity. We suggest that such temporal variability is a widespread phenomenon in the brain which should be modeled.

A common goodness-of-fit framework for neural population models using marked point process time-rescaling.

A critical component of any statistical modeling procedure is the ability to assess the goodness-of-fit between a model and observed data. For spike train models of individual neurons, many goodness-of-fit measures rely on the time-rescaling theorem and assess model quality using rescaled spike times. Recently, there has been increasing interest in statistical models that describe the simultaneous spiking activity of neuron populations, either in a single brain region or across brain regions. Classically, such models have used spike sorted data to describe relationships between the identified neurons, but more recently clusterless modeling methods have been used to describe population activity using a single model. Here we develop a generalization of the time-rescaling theorem that enables comprehensive goodness-of-fit analysis for either of these classes of population models. We use the theory of marked point processes to model population spiking activity, and show that under the correct model, each spike can be rescaled individually to generate a uniformly distributed set of events in time and the space of spike marks. After rescaling, multiple well-established goodness-of-fit procedures and statistical tests are available. We demonstrate the application of these methods both to simulated data and real population spiking in rat hippocampus. We have made the MATLAB and Python code used for the analyses in this paper publicly available through our Github repository at https://github.com/Eden-Kramer-Lab/popTRT .

Synchronous Spike Patterns in Macaque Motor Cortex during an Instructed-Delay Reach-to-Grasp Task.

UNLABELLED: The computational role of spike time synchronization at millisecond precision among neurons in the cerebral cortex is hotly debated. Studies performed on data of limited size provided experimental evidence that low-order correlations occur in relation to behavior. Advances in electrophysiological technology to record from hundreds of neurons simultaneously provide the opportunity to observe coordinated spiking activity of larger populations of cells. We recently published a method that combines data mining and statistical evaluation to search for significant patterns of synchronous spikes in massively parallel spike trains (Torre et al., 2013). The method solves the computational and multiple testing problems raised by the high dimensionality of the data. In the current study, we used our method on simultaneous recordings from two macaque monkeys engaged in an instructed-delay reach-to-grasp task to determine the emergence of spike synchronization in relation to behavior. We found a multitude of synchronous spike patterns aligned in both monkeys along a preferential mediolateral orientation in brain space. The occurrence of the patterns is highly specific to behavior, indicating that different behaviors are associated with the synchronization of different groups of neurons ("cell assemblies"). However, pooled patterns that overlap in neuronal composition exhibit no specificity, suggesting that exclusive cell assemblies become active during different behaviors, but can recruit partly identical neurons. These findings are consistent across multiple recording sessions analyzed across the two monkeys. SIGNIFICANCE STATEMENT: Neurons in the brain communicate via electrical impulses called spikes. How spikes are coordinated to process information is still largely unknown. Synchronous spikes are effective in triggering a spike emission in receiving neurons and have been shown to occur in relation to behavior in a number of studies on simultaneous recordings of few neurons. We recently published a method to extend this type of investigation to larger data. Here, we apply it to simultaneous recordings of hundreds of neurons from the motor cortex of macaque monkeys performing a motor task. Our analysis reveals groups of neurons selectively synchronizing their activity in relation to behavior, which sheds new light on the role of synchrony in information processing in the cerebral cortex.

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.

Analysis of spontaneous activity of superficial dorsal horn neurons in vitro: neuropathy-induced changes.

The superficial dorsal horn contains large numbers of interneurons which process afferent and descending information to generate the spinal nociceptive message. Here, we set out to evaluate whether adjustments in patterns and/or temporal correlation of spontaneous discharges of these neurons are involved in the generation of central sensitization caused by peripheral nerve damage. Multielectrode arrays were used to record from discrete groups of such neurons in slices from control or nerve damaged mice. Whole-cell recordings of individual neurons were also obtained. A large proportion of neurons recorded extracellularly showed well-defined patterns of spontaneous firing. Clock-like neurons (CL) showed regular discharges at ∼6 Hz and represented 9 % of the sample in control animals. They showed a tonic-firing pattern to direct current injection and depolarized membrane potentials. Irregular fast-burst neurons (IFB) produced short-lasting high-frequency bursts (2-5 spikes at ∼100 Hz) at irregular intervals and represented 25 % of the sample. They showed bursting behavior upon direct current injection. Of the pairs of neurons recorded, 10 % showed correlated firing. Correlated pairs always included an IFB neuron. After nerve damage, the mean spontaneous firing frequency was unchanged, but the proportion of CL increased significantly (18 %) and many of these neurons appeared to acquire a novel low-threshold A-fiber input. Similarly, the percentage of IFB neurons was unaltered, but synchronous firing was increased to 22 % of the pairs studied. These changes may contribute to transform spinal processing of nociceptive inputs following peripheral nerve damage. The specific roles that these neurons may play are discussed.

Joint analysis of spikes and local field potentials using copula.

Recent technological advances, which allow for simultaneous recording of spikes and local field potentials (LFPs) at multiple sites in a given cortical area or across different areas, have greatly increased our understanding of signal processing in brain circuits. Joint analysis of simultaneously collected spike and LFP signals is an important step to explicate how the brain orchestrates information processing. In this contribution, we present a novel statistical framework based on Gaussian copula to jointly model spikes and LFP. In our approach, we use copula to link separate, marginal regression models to construct a joint regression model, in which the binary-valued spike train data are modeled using generalized linear model (GLM) and the continuous-valued LFP data are modeled using linear regression. Model parameters can be efficiently estimated via maximum-likelihood. In particular, we show that our model offers a means to statistically detect directional influence between spikes and LFP, akin to Granger causality measure, and that we are able to assess its statistical significance by conducting a Wald test. Through extensive simulations, we also show that our method is able to reliably recover the true model used to generate the data. To demonstrate the effectiveness of our approach in real setting, we further apply the method to a mixed neural dataset, consisting of spikes and LFP simultaneously recorded from the visual cortex of a monkey performing a contour detection task.

A comparison of computational methods for detecting bursts in neuronal spike trains and their application to human stem cell-derived neuronal networks.

Accurate identification of bursting activity is an essential element in the characterization of neuronal network activity. Despite this, no one technique for identifying bursts in spike trains has been widely adopted. Instead, many methods have been developed for the analysis of bursting activity, often on an ad hoc basis. Here we provide an unbiased assessment of the effectiveness of eight of these methods at detecting bursts in a range of spike trains. We suggest a list of features that an ideal burst detection technique should possess and use synthetic data to assess each method in regard to these properties. We further employ each of the methods to reanalyze microelectrode array (MEA) recordings from mouse retinal ganglion cells and examine their coherence with bursts detected by a human observer. We show that several common burst detection techniques perform poorly at analyzing spike trains with a variety of properties. We identify four promising burst detection techniques, which are then applied to MEA recordings of networks of human induced pluripotent stem cell-derived neurons and used to describe the ontogeny of bursting activity in these networks over several months of development. We conclude that no current method can provide "perfect" burst detection results across a range of spike trains; however, two burst detection techniques, the MaxInterval and logISI methods, outperform compared with others. We provide recommendations for the robust analysis of bursting activity in experimental recordings using current techniques.

Beta oscillations in the parkinsonian primate: Similar oscillations across different populations.

Parkinson's disease (PD) is characterized by excessive beta band oscillations (BBO) in neuronal spiking activity across basal ganglia (BG) nuclei. High frequency stimulation of the subthalamic nucleus, an effective treatment for PD, suppresses these oscillations. There is still a heated debate on the origin and propagation of BBO and their association to clinical symptoms. The key prerequisite in addressing these issues is to obtain an accurate estimation of the subpopulation of oscillatory neurons and the magnitude of their oscillations. Studies have shown that neurons in different BG nuclei vary dramatically in the magnitude of their oscillations. However, the stochastic nature of neuronal activity subsamples the oscillatory neuronal rate functions, thus causing standard spectral analysis methods to be dramatically biased by biological and experimental factors such as variations in the neuronal firing rate across BG nuclei. In order to overcome these biases, and directly analyze the expression of BBO within BG nuclei, we used a novel objective method, the modulation index. This method reveals that unlike previous spectral results, individual neurons in the different nuclei display similar magnitudes of oscillations, whereas only the size of the oscillatory subpopulation varies between nuclei. During stimulation, the magnitude of the BBO does not change but the fraction of oscillatory neurons decreases in the globus pallidus internus, leading to a significant change in BG output. This non-biased oscillation quantification thus enables the reconstruction of oscillations at the single neuron and nuclei population levels, and calls for a reassessment of the role of BBO during PD.

Skellam process with resetting: a neural spike train model.

This paper introduces the Skellam process with resetting. Resetting is a modification that accommodates the modeling of neural spike trains. We show this as a biologically plausible model, which codes the information content of neural spike trains with three, potentially, time-varying functions. We show that the interspike interval distribution under this model follows a mixture of gamma distributions, a flexible class covering a wide range of commonly used models. Through simulation studies and the analyses of connected retinal ganglion and lateral geniculate nucleus cells, we evaluate the performance of this model. Copyright © 2016 John Wiley & Sons, Ltd.

A frequency-resolved mutual information rate and its application to neural systems.

The encoding and processing of time-dependent signals into sequences of action potentials of sensory neurons is still a challenging theoretical problem. Although, with some effort, it is possible to quantify the flow of information in the model-free framework of Shannon's information theory, this yields just a single number, the mutual information rate. This rate does not indicate which aspects of the stimulus are encoded. Several studies have identified mechanisms at the cellular and network level leading to low- or high-pass filtering of information, i.e., the selective coding of slow or fast stimulus components. However, these findings rely on an approximation, specifically, on the qualitative behavior of the coherence function, an approximate frequency-resolved measure of information flow, whose quality is generally unknown. Here, we develop an assumption-free method to measure a frequency-resolved information rate about a time-dependent Gaussian stimulus. We demonstrate its application for three paradigmatic descriptions of neural firing: an inhomogeneous Poisson process that carries a signal in its instantaneous firing rate; an integrator neuron (stochastic integrate-and-fire model) driven by a time-dependent stimulus; and the synchronous spikes fired by two commonly driven integrator neurons. In agreement with previous coherence-based estimates, we find that Poisson and integrate-and-fire neurons are broadband and low-pass filters of information, respectively. The band-pass information filtering observed in the coherence of synchronous spikes is confirmed by our frequency-resolved information measure in some but not all parameter configurations. Our results also explicitly show how the response-response coherence can fail as an upper bound on the information rate.

Quantification and classification of neuronal responses in kernel-smoothed peristimulus time histograms.

Peristimulus time histograms are a widespread form of visualizing neuronal responses. Kernel convolution methods transform these histograms into a smooth, continuous probability density function. This provides an improved estimate of a neuron's actual response envelope. We here develop a classifier, called the h-coefficient, to determine whether time-locked fluctuations in the firing rate of a neuron should be classified as a response or as random noise. Unlike previous approaches, the h-coefficient takes advantage of the more precise response envelope estimation provided by the kernel convolution method. The h-coefficient quantizes the smoothed response envelope and calculates the probability of a response of a given shape to occur by chance. We tested the efficacy of the h-coefficient in a large data set of Monte Carlo simulated smoothed peristimulus time histograms with varying response amplitudes, response durations, trial numbers, and baseline firing rates. Across all these conditions, the h-coefficient significantly outperformed more classical classifiers, with a mean false alarm rate of 0.004 and a mean hit rate of 0.494. We also tested the h-coefficient's performance in a set of neuronal responses recorded in humans. The algorithm behind the h-coefficient provides various opportunities for further adaptation and the flexibility to target specific parameters in a given data set. Our findings confirm that the h-coefficient can provide a conservative and powerful tool for the analysis of peristimulus time histograms with great potential for future development.

Hemodynamic and electrophysiological spontaneous low-frequency oscillations in the cortex: directional influences revealed by Granger causality.

We used a combined electrophysiological/hemodynamic system to examine low-frequency oscillations (LFOs) in spontaneous neuronal activities (spike trains and local field potentials) and hemodynamic signals (cerebral blood flow) recorded from the anesthetized rat somatosensory and visual cortices. The laser Doppler flowmetry (LDF) probe was tilted slightly to approach the area in which a microelectrode array (MEA) was implanted for simultaneous recordings. Spike trains (STs) were converted into continuous-time rate functions (CRFs) using the ST instantaneous firing rates. LFOs were detected for all three of the components using the multi-taper method (MTM). The frequencies of these LFOs ranged from 0.052 to 0.167 Hz (mean±SD, 0.10±0.026 Hz) for cerebral blood flow (CBF), from 0.027 to 0.26 Hz (mean±SD, 0.12±0.041 Hz) for the CRFs of the STs and from 0.04 to 0.19 Hz (mean±SD, 0.11±0.035 Hz) for local field potentials (LFPs). We evaluated the Granger causal relationships of spontaneous LFOs among CBF, LFPs and CRFs using Granger causality (GC) analysis. Significant Granger causal relationships were observed from LFPs to CBF, from STs to CBF and from LFPs to STs at approximately 0.1 Hz. The present results indicate that spontaneous LFOs exist not only in hemodynamic components but also in neuronal activities of the rat cortex. To the best of our knowledge, the present study is the first to identify Granger causal influences among CBF, LFPs and STs and show that spontaneous LFOs carry important Granger causal influences from neural activities to hemodynamic signals.