Spontaneous activity in cortical neurons is stereotyped and non-Poisson.

Neurons fire even in the absence of sensory stimulation or task demands. Numerous theoretical studies have modeled this spontaneous activity as a Poisson process with uncorrelated intervals between successive spikes and a variance in firing rate equal to the mean. Experimental tests of this hypothesis have yielded variable results, though most have concluded that firing is not Poisson. However, these tests say little about the ways firing might deviate from randomness. Nor are they definitive because many different distributions can have equal means and variances. Here, we characterized spontaneous spiking patterns in extracellular recordings from monkey, cat, and mouse cerebral cortex neurons using rate-normalized spike train autocorrelation functions (ACFs) and a logarithmic timescale. If activity was Poisson, this function should be flat. This was almost never the case. Instead, ACFs had diverse shapes, often with characteristic peaks in the 1-700 ms range. Shapes were stable over time, up to the longest recording periods used (51 min). They did not fall into obvious clusters. ACFs were often unaffected by visual stimulation, though some abruptly changed during brain state shifts. These behaviors may have their origin in the intrinsic biophysics and dendritic anatomy of the cells or in the inputs they receive.

Mesoscale cortex-wide neural dynamics predict self-initiated actions in mice several seconds prior to movement.

Volition - the sense of control or agency over one's voluntary actions - is widely recognized as the basis of both human subjective experience and natural behavior in nonhuman animals. Several human studies have found peaks in neural activity preceding voluntary actions, for example the readiness potential (RP), and some have shown upcoming actions could be decoded even before awareness. Others propose that random processes underlie and explain pre-movement neural activity. Here, we seek to address these issues by evaluating whether pre-movement neural activity in mice contains structure beyond that present in random neural activity. Implementing a self-initiated water-rewarded lever-pull paradigm in mice while recording widefield [Ca++] neural activity we find that cortical activity changes in variance seconds prior to movement and that upcoming lever pulls could be predicted between 3 and 5 s (or more in some cases) prior to movement. We found inhibition of motor cortex starting at approximately 5 s prior to lever pulls and activation of motor cortex starting at approximately 2 s prior to a random unrewarded left limb movement. We show that mice, like humans, are biased toward commencing self-initiated actions during specific phases of neural activity but that the pre-movement neural code changes over time in some mice and is widely distributed as behavior prediction improved when using all vs. single cortical areas. These findings support the presence of structured multi-second neural dynamics preceding self-initiated action beyond that expected from random processes. Our results also suggest that neural mechanisms underlying self-initiated action could be preserved between mice and humans.

Uncovering the effect of different brain regions on behavioral classification using recurrent neural networks.

As our ability to record neural activity from a larger number of brain areas increases, we need to develop tools to understand how this activity is related to ongoing behavior. Recurrent neural networks (RNNs) have been shown to perform successful classification for sequence data. However, they are black box models: once trained, it is difficult to uncover the mechanisms that they are using to classify. In this study, we analyze the effect of RNNs on classifying behavior using a simulated dataset and a widefield neural activity dataset as mice perform a self-initiated behavior. We show that RNNs are comparable to, or outperform, traditional classification methods such as Support Vector Machine (SVM), and can also lead to accurate prediction of behavior. Using dimensionality reduction, we visualize the activity of the RNNs to better understand the classification mechanisms of the RNNs. Finally, we are able to accurately pinpoint the effect of different regions on behavioral classification. This study highlights the utility and interpretability of RNNs while classifying behavior using neural activity from different regions.

Voltage distributions in extracellular brain recordings.

Extracellular recordings of brain voltage signals have many uses, including the identification of spikes and the characterization of brain states via analysis of local field potential (LFP) or EEG recordings. Though the factors underlying the generation of these signals are time varying and complex, their analysis may be facilitated by an understanding of their statistical properties. To this end, we analyzed the voltage distributions of high-pass extracellular recordings from a variety of structures, including cortex, thalamus, and hippocampus, in monkeys, cats, and rodents. We additionally investigated LFP signals in these recordings as well as human EEG signals obtained during different sleep stages. In all cases, the distributions were accurately described by a Gaussian within ±1.5 standard deviations from zero. Outside these limits, voltages tended to be distributed exponentially, that is, they fell off linearly on log-linear frequency plots, with variable heights and slopes. A possible explanation for this is that sporadically and independently occurring events with individual Gaussian size distributions can sum to produce approximately exponential distributions. For the high-pass recordings, a second explanation results from a model of the noisy behavior of ion channels that produce action potentials via Hodgkin-Huxley kinetics. The distributions produced by this model, relative to the averaged potential, were also Gaussian with approximately exponential flanks. The model also predicted time-varying noise distributions during action potentials, which were observed in the extracellular spike signals. These findings suggest a principled method for detecting spikes in high-pass recordings and transient events in LFP and EEG signals.NEW & NOTEWORTHY We show that the voltage distributions in brain recordings, including high-pass extracellular recordings, the LFP, and human EEG, are accurately described by a Gaussian within ±1.5 standard deviations from zero, with heavy, exponential tails outside these limits. This offers a principled way of setting event detection thresholds in high-pass recordings. It also offers a means for identifying event-like, transient signals in LFP and EEG recordings which may correlate with other neural phenomena.

LFP clustering in cortex reveals a taxonomy of Up states and near-millisecond, ordered phase-locking in cortical neurons.

During slow-wave sleep and anesthesia, mammalian cortex exhibits a synchronized state during which neurons shift from a largely nonfiring to a firing state, known as an Up-state transition. Up-state transitions may constitute the default activity pattern of the entire cortex (Neske GT. Front Neural Circuits 9: 88, 2016) and could be critical to understanding cortical function, yet the genesis of such transitions and their interaction with single neurons is not well understood. It was recently shown that neurons firing at rates >2 Hz fire spikes in a stereotyped order during Up-state transitions (Luczak A, McNaughton BL, Harris KD. Nat Rev Neurosci 16: 745-755, 2015), yet it is still unknown if Up states are homogeneous and whether spiking order is present in neurons with rates <2 Hz (the majority). Using extracellular recordings from anesthetized cats and mice and from naturally sleeping rats, we show for the first time that Up-state transitions can be classified into several types based on the shape of the local field potential (LFP) during each transition. Individual LFP events could be localized in time to within 1-4 ms, more than an order of magnitude less than in previous studies. The majority of recorded neurons synchronized their firing to within ±5-15 ms relative to each Up-state transition. Simultaneous electrophysiology and wide-field imaging in mouse confirmed that LFP event clusters are cortex-wide phenomena. Our findings show that Up states are of different types and point to the potential importance of temporal order and millisecond-scale signaling by cortical neurons.NEW & NOTEWORTHY During cortical Up-state transitions in sleep and anesthesia, neurons undergo brief periods of increased firing in an order similar to that occurring in awake states. We show that these transitions can be classified into distinct types based on the shape of the local field potential. Transition times can be defined to <5 ms. Most neurons synchronize their firing to within ±5-15 ms of the transitions and fire in a consistent order.

A Visual Guide to Sorting Electrophysiological Recordings Using 'SpikeSorter'.

Few stand-alone software applications are available for sorting spikes from recordings made with multi-electrode arrays. Ideally, an application should be user friendly with a graphical user interface, able to read data files in a variety of formats, and provide users with a flexible set of tools giving them the ability to detect and sort extracellular voltage waveforms from different units with some degree of reliability. Previously published spike sorting methods are now available in a software program, SpikeSorter, intended to provide electrophysiologists with a complete set of tools for sorting, starting from raw recorded data file and ending with the export of sorted spikes times. Procedures are automated to the extent this is currently possible. The article explains and illustrates the use of the program. A representative data file is opened, extracellular traces are filtered, events are detected and then clustered. A number of problems that commonly occur during sorting are illustrated, including the artefactual over-splitting of units due to the tendency of some units to fire spikes in pairs where the second spike is significantly smaller than the first, and over-splitting caused by slow variation in spike height over time encountered in some units. The accuracy of SpikeSorter's performance has been tested with surrogate ground truth data and found to be comparable to that of other algorithms in current development.

Mapping cortical mesoscopic networks of single spiking cortical or sub-cortical neurons.

Understanding the basis of brain function requires knowledge of cortical operations over wide-spatial scales, but also within the context of single neurons. In vivo, wide-field GCaMP imaging and sub-cortical/cortical cellular electrophysiology were used in mice to investigate relationships between spontaneous single neuron spiking and mesoscopic cortical activity. We make use of a rich set of cortical activity motifs that are present in spontaneous activity in anesthetized and awake animals. A mesoscale spike-triggered averaging procedure allowed the identification of motifs that are preferentially linked to individual spiking neurons by employing genetically targeted indicators of neuronal activity. Thalamic neurons predicted and reported specific cycles of wide-scale cortical inhibition/excitation. In contrast, spike-triggered maps derived from single cortical neurons yielded spatio-temporal maps expected for regional cortical consensus function. This approach can define network relationships between any point source of neuronal spiking and mesoscale cortical maps.

Verification of multichannel electrode array integrity by use of cross-channel correlations.

BACKGROUND: The use of multichannel electrode arrays (MEAs) presents a number of practical challenges to experimenters including correctly labelling different recording channel locations and identifying sites that may be non-functional or short-circuited. These challenges are likely to increase as the number of sites used in recording increases. NEW METHOD: This paper presents a simple method for assessing MEA integrity based on the observation that physiologically induced signal correlations between nearby channels fall off with distance. Channels that violate this relationship are flagged as being potentially problematic. RESULTS: The method is able to present to the user a list of potentially faulty channels for further inspection. Underlying problems include non-functional, shorted and mislocalised channels and channels carrying spurious noisy signals unrelated to those on other channels. COMPARISON WITH EXISTING METHODS: Computational methods which automatically screen MEAs for faulty electrode channels do not appear to exist in the literature. Currently a user would have to examine single channels, or channel pairs, individually, which would be very time-consuming. CONCLUSIONS: Shorted or mislocalised channels may be more prevalent in MEA recordings than users suspect. The paper presents a simple screening method for identifying such channels prior to carrying out spike-sorting.

Spike detection methods for polytrodes and high density microelectrode arrays.

This paper compares the ability of different methods to detect and resolve spikes recorded extracellularly with polytrode and high-density microelectrode arrays (MEAs). Detecting spikes on such arrays is more complex than with single electrodes or tetrodes since a single spike from a neuron may cause threshold crossings on several adjacent channels, giving rise to multiple events. These initial events have to be recognized as belonging to a single spike. Combining them is, in essence, a clustering problem. A conflicting need is to be able to resolve spike waveforms that occur close together in space and time. We first evaluated three different detection methods, using simulated data in which spike shape waveforms obtained from real recordings were added to noise with an amplitude and temporal structure similar to that found in real recordings. Performance was assessed by calculating the percentage of correctly identified spikes vs. the false positive rate. Using the best of these detection methods, two different methods for avoiding multiple detections per spike were tested: one based on windowing and the other based on clustering. Using parameters that avoided spatial and temporal duplication, the spatiotemporal resolution of the two methods was next evaluated. The method based on clustering gave slightly better results. Both methods could resolve spikes occurring 1 ms or more apart, regardless of their spatial separation. There was no restriction on the temporal resolution of spike pairs for units more than 200 µm apart.

Modeling development in retinal afferents: retinotopy, segregation, and ephrinA/EphA mutants.

During neural development, neurons extend axons to target areas of the brain. Through processes of growth, branching and retraction these axons establish stereotypic patterns of connectivity. In the visual system, these patterns include retinotopic organization and the segregation of individual axons onto different subsets of target neurons based on the eye of origin (ocular dominance) or receptive field type (ON or OFF). Characteristic disruptions to these patterns occur when neural activity or guidance molecule expression is perturbed. In this paper we present a model that explains how these developmental patterns might emerge as a result of the coordinated growth and retraction of individual axons and synapses responding to position-specific markers, trophic factors and spontaneous neural activity. This model derives from one presented earlier (Godfrey et al., 2009) but which is here extended to account for a wider range of phenomena than previously described. These include ocular dominance and ON-OFF segregation and the results of altered ephrinA and EphA guidance molecule expression. The model takes into account molecular guidance factors, realistic patterns of spontaneous retinal wave activity, trophic molecules, homeostatic mechanisms, axon branching and retraction rules and intra-axonal signaling mechanisms that contribute to the survival of nearby synapses on an axon. We show that, collectively, these mechanisms can account for a wider range of phenomena than previous models of retino-tectal development.

Spike sorting for polytrodes: a divide and conquer approach.

In order to determine patterns of neural activity, spike signals recorded by extracellular electrodes have to be clustered (sorted) with the aim of ensuring that each cluster represents all the spikes generated by an individual neuron. Many methods for spike sorting have been proposed but few are easily applicable to recordings from polytrodes which may have 16 or more recording sites. As with tetrodes, these are spaced sufficiently closely that signals from single neurons will usually be recorded on several adjacent sites. Although this offers a better chance of distinguishing neurons with similarly shaped spikes, sorting is difficult in such cases because of the high dimensionality of the space in which the signals must be classified. This report details a method for spike sorting based on a divide and conquer approach. Clusters are initially formed by assigning each event to the channel on which it is largest. Each channel-based cluster is then sub-divided into as many distinct clusters as possible. These are then recombined on the basis of pairwise tests into a final set of clusters. Pairwise tests are also performed to establish how distinct each cluster is from the others. A modified gradient ascent clustering (GAC) algorithm is used to do the clustering. The method can sort spikes with minimal user input in times comparable to real time for recordings lasting up to 45 min. Our results illustrate some of the difficulties inherent in spike sorting, including changes in spike shape over time. We show that some physiologically distinct units may have very similar spike shapes. We show that RMS measures of spike shape similarity are not sensitive enough to discriminate clusters that can otherwise be separated by principal components analysis (PCA). Hence spike sorting based on least-squares matching to templates may be unreliable. Our methods should be applicable to tetrodes and scalable to larger multi-electrode arrays (MEAs).

A multi-component model of the developing retinocollicular pathway incorporating axonal and synaptic growth.

During development, neurons extend axons to different brain areas and produce stereotypical patterns of connections. The mechanisms underlying this process have been intensively studied in the visual system, where retinal neurons form retinotopic maps in the thalamus and superior colliculus. The mechanisms active in map formation include molecular guidance cues, trophic factor release, spontaneous neural activity, spike-timing dependent plasticity (STDP), synapse creation and retraction, and axon growth, branching and retraction. To investigate how these mechanisms interact, a multi-component model of the developing retinocollicular pathway was produced based on phenomenological approximations of each of these mechanisms. Core assumptions of the model were that the probabilities of axonal branching and synaptic growth are highest where the combined influences of chemoaffinity and trophic factor cues are highest, and that activity-dependent release of trophic factors acts to stabilize synapses. Based on these behaviors, model axons produced morphologically realistic growth patterns and projected to retinotopically correct locations in the colliculus. Findings of the model include that STDP, gradient detection by axonal growth cones and lateral connectivity among collicular neurons were not necessary for refinement, and that the instructive cues for axonal growth appear to be mediated first by molecular guidance and then by neural activity. Although complex, the model appears to be insensitive to variations in how the component developmental mechanisms are implemented. Activity, molecular guidance and the growth and retraction of axons and synapses are common features of neural development, and the findings of this study may have relevance beyond organization in the retinocollicular pathway.

Feedback decoding of spatially structured population activity in cortical maps.

A mechanism is proposed by which feedback pathways model spatial patterns of feedforward activity in cortical maps. The mechanism can be viewed equivalently as readout of a content-addressable memory or as decoding of a population code. The model is based on the evidence that cortical receptive fields can often be described as a separable product of functions along several dimensions, each represented in a spatially ordered map. Given this, it is shown that for an N-dimensional map, accurate modeling and decoding of x(N) feedforward activity patterns can be done with Nx fibers, N of which must be active at any one time. The proposed mechanism explains several known properties of the cortex and pyramidal neurons: (1) the integration of signals by dendrites with a narrow tangential distribution, that is, apical dendrites; (2) the presence of fast-conducting feedback projections with broad tangential distributions; (3) the multiplicative effects of attention on receptive field profiles; and (4) the existence of multiplicative interactions between subthreshold feedforward inputs to basal dendrites and inputs to apical dendrites.

Python for large-scale electrophysiology.

Electrophysiology is increasingly moving towards highly parallel recording techniques which generate large data sets. We record extracellularly in vivo in cat and rat visual cortex with 54-channel silicon polytrodes, under time-locked visual stimulation, from localized neuronal populations within a cortical column. To help deal with the complexity of generating and analysing these data, we used the Python programming language to develop three software projects: one for temporally precise visual stimulus generation ("dimstim"); one for electrophysiological waveform visualization and spike sorting ("spyke"); and one for spike train and stimulus analysis ("neuropy"). All three are open source and available for download (http://swindale.ecc.ubc.ca/code). The requirements and solutions for these projects differed greatly, yet we found Python to be well suited for all three. Here we present our software as a showcase of the extensive capabilities of Python in neuroscience.

Retinal wave behavior through activity-dependent refractory periods.

In the developing mammalian visual system, spontaneous retinal ganglion cell (RGC) activity contributes to and drives several aspects of visual system organization. This spontaneous activity takes the form of spreading patches of synchronized bursting that slowly advance across portions of the retina. These patches are non-repeating and tile the retina in minutes. Several transmitter systems are known to be involved, but the basic mechanism underlying wave production is still not well-understood. We present a model for retinal waves that focuses on acetylcholine mediated waves but whose principles are adaptable to other developmental stages. Its assumptions are that a) spontaneous depolarizations of amacrine cells drive wave activity; b) amacrine cells are locally connected, and c) cells receiving more input during their depolarization are subsequently less responsive and have longer periods between spontaneous depolarizations. The resulting model produces waves with non-repeating borders and randomly distributed initiation points. The wave generation mechanism appears to be chaotic and does not require neural noise to produce this wave behavior. Variations in parameter settings allow the model to produce waves that are similar in size, frequency, and velocity to those observed in several species. Our results suggest that retinal wave behavior results from activity-dependent refractory periods and that the average velocity of retinal waves depends on the duration a cell is excitatory: longer periods of excitation result in slower waves. In contrast to previous studies, we find that a single layer of cells is sufficient for wave generation. The principles described here are very general and may be adaptable to the description of spontaneous wave activity in other areas of the nervous system.

Visual cortex: more wiggle room for the brain.

Experiments in which one eye of a ferret is removed at birth show subtle effects on the development of visual cortex maps that are in agreement with those predicted by theory.

A model for the thick, thin and pale stripe organization of primate V2.

Models based on the idea of dimension reduction have been successful in describing the patterns of ocular dominance, spatial frequency and orientation preference found in primate V1. It is shown here that this approach can be extended to describe the organization of thick, thin and pale cytochrome oxidase stripes of primate V2 given an appropriately constructed stimulus space which includes a 3-valued variable which co-varies with color, orientation and disparity. The model successfully describes several aspects of V2 organization, including the fact that there are two pale stripes for each thick and thin stripe and the strong tendency for stripes to run perpendicular to the V1 border. In addition it predicts the presence of reversals in the direction of mapping of retinal eccentricity which should be more common in the pale stripes than elsewhere.

Cerebral cortex: the singular precision of visual cortex maps.

A remarkable new technique, two-photon confocal fluorescence microscopy, has revealed an extraordinarily precise organization in the visual cortex. The methodology seems set to become the tool of choice for studying cortical maps.