Neuronal assembly detection and cell membership specification by principal component analysis

Beta bursts in the parkinsonian cortico-basal ganglia network form spatially discrete ensembles

Defining spatial synchronisation of pathological beta oscillations is important, given that many theories linking them to parkinsonian symptoms propose a reduction in the dimensionality of the coding space within and/or across cortico-basal ganglia structures. Such spatial synchronisation could arise from a single process, with widespread entrainment of neurons to the same oscillation. Alternatively, the partially segregated structure of cortico-basal ganglia loops could provide a substrate for multiple ensembles that are independently synchronized at beta frequencies. Addressing this question requires an analytical approach that identifies groups of signals with a statistical tendency for beta synchronisation, which is unachievable using standard pairwise measures. Here, we utilized such an approach on multichannel recordings of background unit activity (BUA) in the external globus pallidus (GP) and subthalamic nucleus (STN) in parkinsonian rats. We employed an adapted version of a principle and independent component analysis-based method commonly used to define assemblies of single neurons (i.e., neurons that are synchronized over short timescales). This analysis enabled us to define whether changes in the power of beta oscillations in local ensembles of neurons (i.e., the BUA recorded from single contacts) consistently covaried over time, forming a "beta ensemble". Multiple beta ensembles were often present in single recordings and could span brain structures. Membership of a beta ensemble predicted significantly higher levels of short latency (<5 ms) synchrony in the raw BUA signal and phase synchronisation with cortical beta oscillations, suggesting that they comprised clusters of neurons that are functionally connected at multiple levels, despite sometimes being non-contiguous in space. Overall, these findings suggest that beta oscillations do not comprise of a single synchronisation process, but rather multiple independent activities that can bind both spatially contiguous and non-contiguous pools of neurons within and across structures. As previously proposed, such ensembles provide a substrate for beta oscillations to constrain the coding space of cortico-basal ganglia circuits.

Circuit dynamics of superficial and deep CA1 pyramidal cells and inhibitory cells in freely moving macaques

Diverse neuron classes in hippocampal CA1 have been identified through the heterogeneity of their cellular/molecular composition. How these classes relate to hippocampal function and the network dynamics that support cognition in primates remains unclear. Here, we report inhibitory functional cell groups in CA1 of freely moving macaques whose diverse response profiles to network states and each other suggest distinct and specific roles in the functional microcircuit of CA1. In addition, pyramidal cells that were grouped by their superficial or deep layer position differed in firing rate, burstiness, and sharp-wave ripple-associated firing. They also showed strata-specific spike-timing interactions with inhibitory cell groups, suggestive of segregated neural populations. Furthermore, ensemble recordings revealed that cell assemblies were preferentially organized according to these strata. These results suggest that hippocampal CA1 in freely moving macaques bears a sublayer-specific circuit organization that may shape its role in cognition.

Circuit dynamics of superficial and deep CA1 pyramidal cells and inhibitory cells in freely-moving macaques

Diverse neuron classes in hippocampal CA1 have been identified through the heterogeneity of their cellular/molecular composition. How these classes relate to hippocampal function and the network dynamics that support cognition in primates remains unclear. Here we report inhibitory functional cell groups in CA1 of freely-moving macaques whose diverse response profiles to network states and each other suggest distinct and specific roles in the functional microcircuit of CA1. In addition, pyramidal cells that were segregated into superficial and deep layers differed in firing rate, burstiness, and sharp-wave ripple-associated firing. They also showed strata-specific spike-timing interactions with inhibitory cell groups, suggestive of segregated neural populations. Furthermore, ensemble recordings revealed that cell assemblies were preferentially organized according to these strata. These results suggest sublayer-specific circuit organization in hippocampal CA1 of the freely-moving macaques that may underlie its role in cognition.

Stimulus encoding by specific inactivation of cortical neurons

Neuronal ensembles are groups of neurons with correlated activity associated with sensory, motor, and behavioral functions. To explore how ensembles encode information, we investigated responses of visual cortical neurons in awake mice using volumetric two-photon calcium imaging during visual stimulation. We identified neuronal ensembles employing an unsupervised model-free algorithm and, besides neurons activated by the visual stimulus (termed "onsemble"), we also find neurons that are specifically inactivated (termed "offsemble"). Offsemble neurons showed faster calcium decay during stimuli, suggesting selective inhibition. In response to visual stimuli, each ensemble (onsemble+offsemble) exhibited small trial-to-trial variability, high orientation selectivity, and superior predictive accuracy for visual stimulus orientation, surpassing the sum of individual neuron activity. Thus, the combined selective activation and inactivation of cortical neurons enhances visual encoding as an emergent and distributed neural code.

Generalization of generative model for neuronal ensemble inference method

Various brain functions that are necessary to maintain life activities materialize through the interaction of countless neurons. Therefore, it is important to analyze functional neuronal network. To elucidate the mechanism of brain function, many studies are being actively conducted on functional neuronal ensemble and hub, including all areas of neuroscience. In addition, recent study suggests that the existence of functional neuronal ensembles and hubs contributes to the efficiency of information processing. For these reasons, there is a demand for methods to infer functional neuronal ensembles from neuronal activity data, and methods based on Bayesian inference have been proposed. However, there is a problem in modeling the activity in Bayesian inference. The features of each neuron's activity have non-stationarity depending on physiological experimental conditions. As a result, the assumption of stationarity in Bayesian inference model impedes inference, which leads to destabilization of inference results and degradation of inference accuracy. In this study, we extend the range of the variable for expressing the neuronal state, and generalize the likelihood of the model for extended variables. By comparing with the previous study, our model can express the neuronal state in larger space. This generalization without restriction of the binary input enables us to perform soft clustering and apply the method to non-stationary neuroactivity data. In addition, for the effectiveness of the method, we apply the developed method to multiple synthetic fluorescence data generated from the electrical potential data in leaky integrated-and-fire model.

Controlling neuronal assemblies: a fundamental function of respiration-related brain oscillations in neuronal networks

Respiration exerts profound influence on cognition, which is presumed to rely on the generation of local respiration-coherent brain oscillations and the entrainment of cortical neurons. Here, we propose an addition to that view by emphasizing the role of respiration in pacing cortical assemblies (i.e., groups of synchronized, coactive neurons). We review recent findings of how respiration directly entrains identified assembly patterns and discuss how respiration-dependent pacing of assembly activations might be beneficial for cognitive functions.

Complexity changes in functional state dynamics suggest focal connectivity reductions

The past two decades have seen an explosion in the methods and directions of neuroscience research. Along with many others, complexity research has rapidly gained traction as both an independent research field and a valuable subdiscipline in computational neuroscience. In the past decade alone, several studies have suggested that psychiatric disorders affect the spatiotemporal complexity of both global and region-specific brain activity (Liu et al., 2013; Adhikari et al., 2017; Li et al., 2018). However, many of these studies have not accounted for the distributed nature of cognition in either the global or regional complexity estimates, which may lead to erroneous interpretations of both global and region-specific entropy estimates. To alleviate this concern, we propose a novel method for estimating complexity. This method relies upon projecting dynamic functional connectivity into a low-dimensional space which captures the distributed nature of brain activity. Dimension-specific entropy may be estimated within this space, which in turn allows for a rapid estimate of global signal complexity. Testing this method on a recently acquired obsessive-compulsive disorder dataset reveals substantial increases in the complexity of both global and dimension-specific activity versus healthy controls, suggesting that obsessive-compulsive patients may experience increased disorder in cognition. To probe the potential causes of this alteration, we estimate subject-level effective connectivity via a Hopf oscillator-based model dynamic model, the results of which suggest that obsessive-compulsive patients may experience abnormally high connectivity across a broad network in the cortex. These findings are broadly in line with results from previous studies, suggesting that this method is both robust and sensitive to group-level complexity alterations.

Scalable and accurate method for neuronal ensemble detection in spiking neural networks

We propose a novel, scalable, and accurate method for detecting neuronal ensembles from a population of spiking neurons. Our approach offers a simple yet powerful tool to study ensemble activity. It relies on clustering synchronous population activity (population vectors), allows the participation of neurons in different ensembles, has few parameters to tune and is computationally efficient. To validate the performance and generality of our method, we generated synthetic data, where we found that our method accurately detects neuronal ensembles for a wide range of simulation parameters. We found that our method outperforms current alternative methodologies. We used spike trains of retinal ganglion cells obtained from multi-electrode array recordings under a simple ON-OFF light stimulus to test our method. We found a consistent stimuli-evoked ensemble activity intermingled with spontaneously active ensembles and irregular activity. Our results suggest that the early visual system activity could be organized in distinguishable functional ensembles. We provide a Graphic User Interface, which facilitates the use of our method by the scientific community.

Spectral estimation for detecting low-dimensional structure in networks using arbitrary null models

Discovering low-dimensional structure in real-world networks requires a suitable null model that defines the absence of meaningful structure. Here we introduce a spectral approach for detecting a network’s low-dimensional structure, and the nodes that participate in it, using any null model. We use generative models to estimate the expected eigenvalue distribution under a specified null model, and then detect where the data network’s eigenspectra exceed the estimated bounds. On synthetic networks, this spectral estimation approach cleanly detects transitions between random and community structure, recovers the number and membership of communities, and removes noise nodes. On real networks spectral estimation finds either a significant fraction of noise nodes or no departure from a null model, in stark contrast to traditional community detection methods. Across all analyses, we find the choice of null model can strongly alter conclusions about the presence of network structure. Our spectral estimation approach is therefore a promising basis for detecting low-dimensional structure in real-world networks, or lack thereof.

Patterned activation of action potential patterns during offline states in the neocortex: replay and non-replay

Action potential generation (spiking) in the neocortex is organized into repeating non-random patterns during both awake experiential states and non-engaged states ranging from inattention to sleep to anaesthesia—and even occur in slice preparations. Repeating patterns in a given population of neurons between states may imply a common means by which cortical networks can be engaged despite brain state changes, but super-imposed on this common firing is a variability that is both specific to ongoing inputs and can be re-shaped by experience. This similarity with specifically induced variance may allow for a range of processes including perception, memory consolidation and network homeostasis. Here, we review how patterned activity in neocortical populations has been studied and what it may imply for a cortex that must be both static and plastic.

This article is part of the Theo Murphy meeting issue ‘Memory reactivation: replaying events past, present and future’.

On the methods for reactivation and replay analysis

A major task in the history of neurophysiology has been to relate patterns of neural activity to ongoing external stimuli. More recently, this approach has branched out to relating current neural activity patterns to external stimuli or experiences that occurred in the past or future. Here, we aim to review the large body of methodological approaches used towards this goal, and to assess the assumptions each makes with reference to the statistics of neural data that are commonly observed. These methods primarily fall into two categories, those that quantify zero-lag relationships without examining temporal evolution, termed reactivation, and those that quantify the temporal structure of changing activity patterns, termed replay. However, no two studies use the exact same approach, which prevents an unbiased comparison between findings. These observations should instead be validated by multiple and, if possible, previously established tests. This will help the community to speak a common language and will eventually provide tools to study, more generally, the organization of neuronal patterns in the brain. This article is part of the Theo Murphy meeting issue 'Memory reactivation: replaying events past, present and future'.

Endocannabinoid degradation inhibitors ameliorate neuronal and synaptic alterations following traumatic brain injury

Our previous work showed that lateral fluid percussion injury to the sensorimotor cortex (SMC) of anesthetized rats increased neuronal synaptic hyperexcitability in layer 5 (L5) neurons in ex vivo brain slices 10 days postinjury. Furthermore, endocannabinoid (EC) degradation inhibition via intraperitoneal JZL184 injection 30 min postinjury attenuated synaptic hyperexcitability. This study tested the hypothesis that traumatic brain injury (TBI) induces synaptic and intrinsic neuronal alterations of L5 SMC pyramidal neurons and that these alterations are significantly attenuated by in vivo post-TBI treatment with EC degradation inhibitors. We tested the effects of systemically administered EC degradation enzyme inhibitors (JZL184, MJN110, URB597, or JZL195) with differential selectivity for fatty acid amide hydrolase (FAAH) and monoacylglycerol lipase (MAGL) on electrophysiological parameters in SMC neurons of TBI- and sham-treated rats 10 days post-TBI. We recorded intrinsic neuronal properties, including resting membrane voltage, input resistance, spike threshold, spiking responses to current input, voltage "sag" (rebound response to hyperpolarization-activated inward current), and burst firing. We also measured the frequency and amplitude of spontaneous excitatory postsynaptic currents. We then used the aggregate parameter sets (intrinsic + synaptic properties) to apply a machine learning classification algorithm to quantitatively compare neural population responses from each experimental group. Collectively, our electrophysiological and computational results indicate that sham neurons are the most distinguishable from TBI neurons. Administration of EC degradation inhibitors post-TBI exerted varying degrees of rescue, approximating the neuronal phenotype of sham neurons, with neurons from TBI/JZL195 (a dual MAGL/FAAH inhibitor) being most similar to neurons from sham rats.NEW & NOTEWORTHY This study elucidates neuronal properties altered by traumatic brain injury (TBI) in layer 5 of sensorimotor cortex, which may be implicated in post-TBI circuit dysfunction. We compared effects of systemic administration of four different endocannabinoid degradation inhibitors within a clinically relevant window postinjury. Electrophysiological measures and using a machine learning classification algorithm collectively suggest that pharmacological inhibitors targeting both monoacylglycerol lipase and fatty acid amide hydrolase (e.g., JZL195) may be most efficacious in attenuating TBI-induced neuronal dysfunction at site of injury.

Regulation of gamma-frequency oscillation by feedforward inhibition: A computational modeling study

Throughout the brain, reciprocally connected excitatory and inhibitory neurons interact to produce gamma-frequency oscillations. The emergent gamma rhythm synchronizes local neural activity and helps to select which cells should fire in each cycle. We previously found that such excitation-inhibition microcircuits, however, have a potentially significant caveat: the frequency of the gamma oscillation and the level of selection (i.e., the percentage of cells that are allowed to fire) vary with the magnitude of the input signal. In networks with varying levels of brain activity, such a feature may produce undesirable instability on the time and spatial structure of the neural signal with a potential for adversely impacting important neural processing mechanisms. Here we propose that feedforward inhibition solves the latter instability problem of the excitation-inhibition microcircuit. Using computer simulations, we show that the feedforward inhibitory signal reduces the dependence of both the frequency of population oscillation and the level of selection on the magnitude of the input excitation. Such a mechanism can produce stable gamma oscillations with its frequency determined only by the properties of the feedforward network, as observed in the hippocampus. As feedforward and feedback inhibition motifs commonly appear together in the brain, we hypothesize that their interaction underlies a robust implementation of general computational principles of neural processing involved in several cognitive tasks, including the formation of cell assemblies and the routing of information between brain areas.

Bayesian inference of neuronal assemblies

In many areas of the brain, both spontaneous and stimulus-evoked activity can manifest as synchronous activation of neuronal assemblies. The characterization of assembly structure and dynamics provides important insights into how brain computations are distributed across neural networks. The proliferation of experimental techniques for recording the activity of neuronal assemblies calls for a comprehensive statistical method to describe, analyze and characterize these high dimensional datasets. The performance of existing methods for defining assemblies is sensitive to noise and stochasticity in neuronal firing patterns and assembly heterogeneity. To address these problems, we introduce a generative hierarchical model of synchronous activity to describe the organization of neurons into assemblies. Unlike existing methods, our analysis provides a simultaneous estimation of assembly composition, dynamics and within-assembly statistical features, such as the levels of activity, noise and assembly synchrony. We have used our method to characterize population activity throughout the tectum of larval zebrafish, allowing us to make statistical inference on the spatiotemporal organization of tectal assemblies, their composition and the logic of their interactions. We have also applied our method to functional imaging and neuropixels recordings from the mouse, allowing us to relate the activity of identified assemblies to specific behaviours such as running or changes in pupil diameter.

Unsupervised Detection of Cell-Assembly Sequences by Similarity-Based Clustering

Neurons which fire in a fixed temporal pattern (i.e., "cell assemblies") are hypothesized to be a fundamental unit of neural information processing. Several methods are available for the detection of cell assemblies without a time structure. However, the systematic detection of cell assemblies with time structure has been challenging, especially in large datasets, due to the lack of efficient methods for handling the time structure. Here, we show a method to detect a variety of cell-assembly activity patterns, recurring in noisy neural population activities at multiple timescales. The key innovation is the use of a computer science method to comparing strings ("edit similarity"), to group spikes into assemblies. We validated the method using artificial data and experimental data, which were previously recorded from the hippocampus of male Long-Evans rats and the prefrontal cortex of male Brown Norway/Fisher hybrid rats. From the hippocampus, we could simultaneously extract place-cell sequences occurring on different timescales during navigation and awake replay. From the prefrontal cortex, we could discover multiple spike sequences of neurons encoding different segments of a goal-directed task. Unlike conventional event-driven statistical approaches, our method detects cell assemblies without creating event-locked averages. Thus, the method offers a novel analytical tool for deciphering the neural code during arbitrary behavioral and mental processes.

Unsupervised identification of states from voltage recordings of neural networks

BACKGROUND

Modern techniques for multi-neuronal recording produce large amounts of data. There is no automatic procedure for the identification of states in recurrent voltage patterns.

NEW METHOD

We propose NetSAP (Network States And Pathways), a data-driven analysis method that is able to recognize multi-neuron voltage patterns (states). To capture the subtle differences between snapshots in voltage recordings, NetSAP infers the underlying functional neural network in a time-resolved manner with a sliding window approach. Then NetSAP identifies states from a reordering of the time series of inferred networks according to a user-defined metric. The procedure for unsupervised identification of states was developed originally for the analysis of molecular dynamics simulations of proteins.

RESULTS

We tested NetSAP on neural network simulations of GABAergic inhibitory interneurons. Most simulation parameters are chosen to reproduce literature observations, and we keep noise terms as control parameters to regulate the coherence of the simulated signals. NetSAP is able to identify multiple states even in the case of high internal noise and low signal coherence. We provide evidence that NetSAP is robust for networks with up to about 50% of the neurons spiking randomly. NetSAP is scalable and its code is open source.

COMPARISON WITH EXISTING METHODS

NetSAP outperforms common analysis techniques, such as PCA and k-means clustering, on a simulated recording of voltage traces of 50 neurons.

CONCLUSIONS

NetSAP analysis is an efficient tool to identify voltage patterns from neuronal recordings.

Hippocampal Reactivation Extends for Several Hours Following Novel Experience

New memories are believed to be consolidated over several hours of post-task sleep. The reactivation or "replay" of hippocampal cell assemblies has been proposed to provide a key mechanism for this process. However, previous studies have indicated that such replay is restricted to the first 10-30 min of post-task sleep, suggesting that it has a limited role in memory consolidation. We performed long-duration recordings in sleeping and behaving male rats and applied methods for evaluating the reactivation of neurons in pairs as well as in larger ensembles while controlling for the continued activation of ensembles already present during pre-task sleep ("preplay"). We found that cell assemblies reactivate for up to 10 h, with a half-maximum timescale of ∼6 h, in sleep following novel experience, even when corrected for preplay. We further confirmed similarly prolonged reactivation in post-task sleep of rats in other datasets that used behavior in novel environments. In contrast, we saw limited reactivation in sleep following behavior in familiar environments. Overall, our findings reconcile the duration of replay with the timescale attributed to cellular memory consolidation and provide strong support for an integral role of replay in memory.SIGNIFICANCE STATEMENT Neurons that are active during an experience reactivate again afterward during rest and sleep. This replay of ensembles of neurons has been proposed to help strengthen memories, but it has also been reported that replay occurs only in the first 10-30 min of sleep, suggesting a circumscribed role. We performed long-duration recordings in the hippocampus of rats and found that replay persists for several hours in sleep following novel experience, far beyond the limits found in previous reports based on shorter recordings. These findings reconcile the duration of replay with the hours-long timescales attributed to memory consolidation.

Detecting neural assemblies in calcium imaging data

BACKGROUND

Activity in populations of neurons often takes the form of assemblies, where specific groups of neurons tend to activate at the same time. However, in calcium imaging data, reliably identifying these assemblies is a challenging problem, and the relative performance of different assembly-detection algorithms is unknown.

RESULTS

To test the performance of several recently proposed assembly-detection algorithms, we first generated large surrogate datasets of calcium imaging data with predefined assembly structures and characterised the ability of the algorithms to recover known assemblies. The algorithms we tested are based on independent component analysis (ICA), principal component analysis (Promax), similarity analysis (CORE), singular value decomposition (SVD), graph theory (SGC), and frequent item set mining (FIM-X). When applied to the simulated data and tested against parameters such as array size, number of assemblies, assembly size and overlap, and signal strength, the SGC and ICA algorithms and a modified form of the Promax algorithm performed well, while PCA-Promax and FIM-X did less well, for instance, showing a strong dependence on the size of the neural array. Notably, we identified additional analyses that can improve their importance. Next, we applied the same algorithms to a dataset of activity in the zebrafish optic tectum evoked by simple visual stimuli, and found that the SGC algorithm recovered assemblies closest to the averaged responses.

CONCLUSIONS

Our findings suggest that the neural assemblies recovered from calcium imaging data can vary considerably with the choice of algorithm, but that some algorithms reliably perform better than others. This suggests that previous results using these algorithms may need to be reevaluated in this light.

Uncovering Neuronal Networks Defined by Consistent Between-Neuron Spike Timing from Neuronal Spike Recordings

It is widely assumed that distributed neuronal networks are fundamental to the functioning of the brain. Consistent spike timing between neurons is thought to be one of the key principles for the formation of these networks. This can involve synchronous spiking or spiking with time delays, forming spike sequences when the order of spiking is consistent. Finding networks defined by their sequence of time-shifted spikes, denoted here as spike timing networks, is a tremendous challenge. As neurons can participate in multiple spike sequences at multiple between-spike time delays, the possible complexity of networks is prohibitively large. We present a novel approach that is capable of (1) extracting spike timing networks regardless of their sequence complexity, and (2) that describes their spiking sequences with high temporal precision. We achieve this by decomposing frequency-transformed neuronal spiking into separate networks, characterizing each network’s spike sequence by a time delay per neuron, forming a spike sequence timeline. These networks provide a detailed template for an investigation of the experimental relevance of their spike sequences. Using simulated spike timing networks, we show network extraction is robust to spiking noise, spike timing jitter, and partial occurrences of the involved spike sequences. Using rat multineuron recordings, we demonstrate the approach is capable of revealing real spike timing networks with sub-millisecond temporal precision. By uncovering spike timing networks, the prevalence, structure, and function of complex spike sequences can be investigated in greater detail, allowing us to gain a better understanding of their role in neuronal functioning.

Neuronal Assemblies Evidence Distributed Interactions within a Tactile Discrimination Task in Rats

Accumulating evidence suggests that neural interactions are distributed and relate to animal behavior, but many open questions remain. The neural assembly hypothesis, formulated by Hebb, states that synchronously active single neurons may transiently organize into functional neural circuits-neuronal assemblies (NAs)-and that would constitute the fundamental unit of information processing in the brain. However, the formation, vanishing, and temporal evolution of NAs are not fully understood. In particular, characterizing NAs in multiple brain regions over the course of behavioral tasks is relevant to assess the highly distributed nature of brain processing. In the context of NA characterization, active tactile discrimination tasks with rats are elucidative because they engage several cortical areas in the processing of information that are otherwise masked in passive or anesthetized scenarios. In this work, we investigate the dynamic formation of NAs within and among four different cortical regions in long-range fronto-parieto-occipital networks (primary somatosensory, primary visual, prefrontal, and posterior parietal cortices), simultaneously recorded from seven rats engaged in an active tactile discrimination task. Our results first confirm that task-related neuronal firing rate dynamics in all four regions is significantly modulated. Notably, a support vector machine decoder reveals that neural populations contain more information about the tactile stimulus than the majority of single neurons alone. Then, over the course of the task, we identify the emergence and vanishing of NAs whose participating neurons are shown to contain more information about animal behavior than randomly chosen neurons. Taken together, our results further support the role of multiple and distributed neurons as the functional unit of information processing in the brain (NA hypothesis) and their link to active animal behavior.

Dynamical networks: Finding, measuring, and tracking neural population activity using network science

Systems neuroscience is in a headlong rush to record from as many neurons at the same time as possible. As the brain computes and codes using neuron populations, it is hoped these data will uncover the fundamentals of neural computation. But with hundreds, thousands, or more simultaneously recorded neurons come the inescapable problems of visualizing, describing, and quantifying their interactions. Here I argue that network science provides a set of scalable, analytical tools that already solve these problems. By treating neurons as nodes and their interactions as links, a single network can visualize and describe an arbitrarily large recording. I show that with this description we can quantify the effects of manipulating a neural circuit, track changes in population dynamics over time, and quantitatively define theoretical concepts of neural populations such as cell assemblies. Using network science as a core part of analyzing population recordings will thus provide both qualitative and quantitative advances to our understanding of neural computation.

Motion contrast in primary visual cortex: a direct comparison of single neuron and population encoding

Features from outside the classical receptive field (CRF) can modulate the stimulus-driven activity of single cells in the primary visual cortex. This modulation, mediated by horizontal and feedback networks, has been extensively described as a variation of firing rate and is considered the basis of processing features as, for example, motion contrast. However, surround influences have also been identified in pairwise spiking or local field coherence. Yet, evidence about co-existence and integration of different neural signatures remains elusive. To compare multiple signatures, we recorded spiking and LFP activity evoked by stimuli exhibiting a motion contrast in the CRFs surround in anesthetized cat primary visual cortex. We chose natural-like scenes over gratings to avoid predominance of simple visual features, which could be easily represented by a rate code. We analyzed firing rates and phase-locking to low-gamma frequency in single cells and neuronal assemblies. Motion contrast was reflected in all measures but in semi-independent populations. Whereas activation of assemblies accompanied single neuron rates, their phase relations were modulated differently. Interestingly, only assembly phase relations mirrored the direction of movement of the surround and were selectively affected by thermal deactivation of visual interhemispheric connections. We argue that motion contrast can be reflected in complementary and superimposed neuronal signatures that can represent different surround features in independent neuronal populations.

Spontaneous Activity in the Zebrafish Tectum Reorganizes over Development and Is Influenced by Visual Experience

Spontaneous patterns of activity in the developing visual system may play an important role in shaping the brain for function. During the period 4-9 dpf (days post-fertilization), larval zebrafish learn to hunt prey, a behavior that is critically dependent on the optic tectum. However, how spontaneous activity develops in the tectum over this period and the effect of visual experience are unknown. Here we performed two-photon calcium imaging of GCaMP6s zebrafish larvae at all days from 4 to 9 dpf. Using recently developed graph theoretic techniques, we found significant changes in both single-cell and population activity characteristics over development. In particular, we identified days 5-6 as a critical moment in the reorganization of the underlying functional network. Altering visual experience early in development altered the statistics of tectal activity, and dark rearing also caused a long-lasting deficit in the ability to capture prey. Thus, tectal development is shaped by both intrinsic factors and visual experience.

An integrated calcium imaging processing toolbox for the analysis of neuronal population dynamics

The development of new imaging and optogenetics techniques to study the dynamics of large neuronal circuits is generating datasets of unprecedented volume and complexity, demanding the development of appropriate analysis tools. We present a comprehensive computational workflow for the analysis of neuronal population calcium dynamics. The toolbox includes newly developed algorithms and interactive tools for image pre-processing and segmentation, estimation of significant single-neuron single-trial signals, mapping event-related neuronal responses, detection of activity-correlated neuronal clusters, exploration of population dynamics, and analysis of clusters' features against surrogate control datasets. The modules are integrated in a modular and versatile processing pipeline, adaptable to different needs. The clustering module is capable of detecting flexible, dynamically activated neuronal assemblies, consistent with the distributed population coding of the brain. We demonstrate the suitability of the toolbox for a variety of calcium imaging datasets. The toolbox open-source code, a step-by-step tutorial and a case study dataset are available at https://github.com/zebrain-lab/Toolbox-Romano-et-al.

Cell assemblies at multiple time scales with arbitrary lag constellations

Hebb's idea of a cell assembly as the fundamental unit of neural information processing has dominated neuroscience like no other theoretical concept within the past 60 years. A range of different physiological phenomena, from precisely synchronized spiking to broadly simultaneous rate increases, has been subsumed under this term. Yet progress in this area is hampered by the lack of statistical tools that would enable to extract assemblies with arbitrary constellations of time lags, and at multiple temporal scales, partly due to the severe computational burden. Here we present such a unifying methodological and conceptual framework which detects assembly structure at many different time scales, levels of precision, and with arbitrary internal organization. Applying this methodology to multiple single unit recordings from various cortical areas, we find that there is no universal cortical coding scheme, but that assembly structure and precision significantly depends on the brain area recorded and ongoing task demands.

DOI:http://dx.doi.org/10.7554/eLife.19428.001

Using Matrix and Tensor Factorizations for the Single-Trial Analysis of Population Spike Trains

Advances in neuronal recording techniques are leading to ever larger numbers of simultaneously monitored neurons. This poses the important analytical challenge of how to capture compactly all sensory information that neural population codes carry in their spatial dimension (differences in stimulus tuning across neurons at different locations), in their temporal dimension (temporal neural response variations), or in their combination (temporally coordinated neural population firing). Here we investigate the utility of tensor factorizations of population spike trains along space and time. These factorizations decompose a dataset of single-trial population spike trains into spatial firing patterns (combinations of neurons firing together), temporal firing patterns (temporal activation of these groups of neurons) and trial-dependent activation coefficients (strength of recruitment of such neural patterns on each trial). We validated various factorization methods on simulated data and on populations of ganglion cells simultaneously recorded in the salamander retina. We found that single-trial tensor space-by-time decompositions provided low-dimensional data-robust representations of spike trains that capture efficiently both their spatial and temporal information about sensory stimuli. Tensor decompositions with orthogonality constraints were the most efficient in extracting sensory information, whereas non-negative tensor decompositions worked well even on non-independent and overlapping spike patterns, and retrieved informative firing patterns expressed by the same population in response to novel stimuli. Our method showed that populations of retinal ganglion cells carried information in their spike timing on the ten-milliseconds-scale about spatial details of natural images. This information could not be recovered from the spike counts of these cells. First-spike latencies carried the majority of information provided by the whole spike train about fine-scale image features, and supplied almost as much information about coarse natural image features as firing rates. Together, these results highlight the importance of spike timing, and particularly of first-spike latencies, in retinal coding.

Neural assemblies revealed by inferred connectivity-based models of prefrontal cortex recordings

We present two graphical model-based approaches to analyse the distribution of neural activities in the prefrontal cortex of behaving rats. The first method aims at identifying cell assemblies, groups of synchronously activating neurons possibly representing the units of neural coding and memory. A graphical (Ising) model distribution of snapshots of the neural activities, with an effective connectivity matrix reproducing the correlation statistics, is inferred from multi-electrode recordings, and then simulated in the presence of a virtual external drive, favoring high activity (multi-neuron) configurations. As the drive increases groups of neurons may activate together, and reveal the existence of cell assemblies. The identified groups are then showed to strongly coactivate in the neural spiking data and to be highly specific of the inferred connectivity network, which offers a sparse representation of the correlation pattern across neural cells. The second method relies on the inference of a Generalized Linear Model, in which spiking events are integrated over time by neurons through an effective connectivity matrix. The functional connectivity matrices inferred with the two approaches are compared. Sampling of the inferred GLM distribution allows us to study the spatio-temporal patterns of activation of neurons within the identified cell assemblies, particularly their activation order: the prevalence of one order with respect to the others is weak and reflects the neuron average firing rates and the strength of the largest effective connections. Other properties of the identified cell assemblies (spatial distribution of coactivation events and firing rates of coactivating neurons) are discussed.

Recoding a cocaine-place memory engram to a neutral engram in the hippocampus

The hippocampus provides the brain's memory system with a subset of neurons holding a map-like representation of each environment experienced. We found in mice that optogenetic silencing those neurons active in an environment unmasked a subset of quiet neurons, enabling the emergence of an alternative map. When applied in a cocaine-paired environment, this intervention neutralized an otherwise long-lasting drug-place preference, showing that recoding a spatial memory engram can alleviate associated maladaptive behavior.

The brain dynamics of linguistic computation

Neural oscillations at distinct frequencies are increasingly being related to a number of basic and higher cognitive faculties. Oscillations enable the construction of coherently organized neuronal assemblies through establishing transitory temporal correlations. By exploring the elementary operations of the language faculty—labeling, concatenation, cyclic transfer—alongside neural dynamics, a new model of linguistic computation is proposed. It is argued that the universality of language, and the true biological source of Universal Grammar, is not to be found purely in the genome as has long been suggested, but more specifically within the extraordinarily preserved nature of mammalian brain rhythms employed in the computation of linguistic structures. Computational-representational theories are used as a guide in investigating the neurobiological foundations of the human “cognome”—the set of computations performed by the nervous system—and new directions are suggested for how the dynamics of the brain (the “dynome”) operate and execute linguistic operations. The extent to which brain rhythms are the suitable neuronal processes which can capture the computational properties of the human language faculty is considered against a backdrop of existing cartographic research into the localization of linguistic interpretation. Particular focus is placed on labeling, the operation elsewhere argued to be species-specific. A Basic Label model of the human cognome-dynome is proposed, leading to clear, causally-addressable empirical predictions, to be investigated by a suggested research program, Dynamic Cognomics. In addition, a distinction between minimal and maximal degrees of explanation is introduced to differentiate between the depth of analysis provided by cartographic, rhythmic, neurochemical, and other approaches to computation.

Sleep-Dependent Reactivation of Ensembles in Motor Cortex Promotes Skill Consolidation

Despite many prior studies demonstrating offline behavioral gains in motor skills after sleep, the underlying neural mechanisms remain poorly understood. To investigate the neurophysiological basis for offline gains, we performed single-unit recordings in motor cortex as rats learned a skilled upper-limb task. We found that sleep improved movement speed with preservation of accuracy. These offline improvements were linked to both replay of task-related ensembles during non-rapid eye movement (NREM) sleep and temporal shifts that more tightly bound motor cortical ensembles to movements; such offline gains and temporal shifts were not evident with sleep restriction. Interestingly, replay was linked to the coincidence of slow-wave events and bursts of spindle activity. Neurons that experienced the most consistent replay also underwent the most significant temporal shift and binding to the motor task. Significantly, replay and the associated performance gains after sleep only occurred when animals first learned the skill; continued practice during later stages of learning (i.e., after motor kinematics had stabilized) did not show evidence of replay. Our results highlight how replay of synchronous neural activity during sleep mediates large-scale neural plasticity and stabilizes kinematics during early motor learning.

During non-REM sleep in rats, consolidation and offline improvements of a recently learned motor skill are linked to synchronous reactivation of task-related neural ensembles.

Sleep has been shown to help in consolidating learned motor tasks. In other words, sleep can induce “offline” gains in a new motor skill even in the absence of further training. However, how sleep induces this change has not been clearly identified. One hypothesis is that consolidation of memories during sleep occurs by “reactivation” of neurons engaged during learning. In this study, we tested this hypothesis by recording populations of neurons in the motor cortex of rats while they learned a new motor skill and during sleep both before and after the training session. We found that subsets of task-relevant neurons formed highly synchronized ensembles during learning. Interestingly, these same neural ensembles were reactivated during subsequent sleep blocks, and the degree of reactivation was correlated with several metrics of motor memory consolidation. Specifically, after sleep, the speed at which animals performed the task while maintaining accuracy was increased, and the activity of the neuronal assembles were more tightly bound to motor action. Further analyses showed that reactivation events occurred episodically and in conjunction with spindle-oscillations—common bursts of brain activity seen during sleep. This observation is consistent with previous findings in humans that spindle-oscillations correlate with consolidation of learned tasks. Our study thus provides insight into the neuronal network mechanism supporting consolidation of motor memory during sleep and may lead to novel interventions that can enhance skill learning in both healthy and injured nervous systems.

Improved estimation and interpretation of correlations in neural circuits

Ambitious projects aim to record the activity of ever larger and denser neuronal populations in vivo. Correlations in neural activity measured in such recordings can reveal important aspects of neural circuit organization. However, estimating and interpreting large correlation matrices is statistically challenging. Estimation can be improved by regularization, i.e. by imposing a structure on the estimate. The amount of improvement depends on how closely the assumed structure represents dependencies in the data. Therefore, the selection of the most efficient correlation matrix estimator for a given neural circuit must be determined empirically. Importantly, the identity and structure of the most efficient estimator informs about the types of dominant dependencies governing the system. We sought statistically efficient estimators of neural correlation matrices in recordings from large, dense groups of cortical neurons. Using fast 3D random-access laser scanning microscopy of calcium signals, we recorded the activity of nearly every neuron in volumes 200 μm wide and 100 μm deep (150-350 cells) in mouse visual cortex. We hypothesized that in these densely sampled recordings, the correlation matrix should be best modeled as the combination of a sparse graph of pairwise partial correlations representing local interactions and a low-rank component representing common fluctuations and external inputs. Indeed, in cross-validation tests, the covariance matrix estimator with this structure consistently outperformed other regularized estimators. The sparse component of the estimate defined a graph of interactions. These interactions reflected the physical distances and orientation tuning properties of cells: The density of positive 'excitatory' interactions decreased rapidly with geometric distances and with differences in orientation preference whereas negative 'inhibitory' interactions were less selective. Because of its superior performance, this 'sparse+latent' estimator likely provides a more physiologically relevant representation of the functional connectivity in densely sampled recordings than the sample correlation matrix.

Spontaneous neuronal network dynamics reveal circuit's functional adaptations for behavior

Spontaneous neuronal activity is spatiotemporally structured, influencing brain computations. Nevertheless, the neuronal interactions underlying these spontaneous activity patterns, and their biological relevance, remain elusive. Here, we addressed these questions using two-photon calcium imaging of intact zebrafish larvae to monitor the neuron-to-neuron spontaneous activity fine structure in the tectum, a region involved in visual spatial detection. Spontaneous activity was organized in topographically compact assemblies, grouping functionally similar neurons rather than merely neighboring ones, reflecting the tectal retinotopic map despite being independent of retinal drive. Assemblies represent all-or-none-like sub-networks shaped by competitive dynamics, mechanisms advantageous for visual detection in noisy natural environments. Notably, assemblies were tuned to the same angular sizes and spatial positions as prey-detection performance in behavioral assays, and their spontaneous activation predicted directional tail movements. Therefore, structured spontaneous activity represents "preferred" network states, tuned to behaviorally relevant features, emerging from the circuit's intrinsic non-linear dynamics, adapted for its functional role.

Monitoring Spiking Activity of Many Individual Neurons in Invertebrate Ganglia

Optical recording with fast voltage sensitive dyes makes it possible, in suitable preparations, to simultaneously monitor the action potentials of large numbers of individual neurons. Here we describe methods for doing this, including considerations of different dyes and imaging systems, methods for correlating the optical signals with their source neurons, procedures for getting good signals, and the use of Independent Component Analysis for spike-sorting raw optical data into single neuron traces. These combined tools represent a powerful approach for large-scale recording of neural networks with high temporal and spatial resolution.

Recent developments in VSD imaging of small neuronal networks

Voltage-sensitive dye (VSD) imaging is a powerful technique that can provide, in single experiments, a large-scale view of network activity unobtainable with traditional sharp electrode recording methods. Here we review recent work using VSDs to study small networks and highlight several results from this approach. Topics covered include circuit mapping, network multifunctionality, the network basis of decision making, and the presence of variably participating neurons in networks. Analytical tools being developed and applied to large-scale VSD imaging data sets are discussed, and the future prospects for this exciting field are considered.

Revealing cell assemblies at multiple levels of granularity

BACKGROUND

Current neuronal monitoring techniques, such as calcium imaging and multi-electrode arrays, enable recordings of spiking activity from hundreds of neurons simultaneously. Of primary importance in systems neuroscience is the identification of cell assemblies: groups of neurons that cooperate in some form within the recorded population.

NEW METHOD

We introduce a simple, integrated framework for the detection of cell-assemblies from spiking data without a priori assumptions about the size or number of groups present. We define a biophysically-inspired measure to extract a directed functional connectivity matrix between both excitatory and inhibitory neurons based on their spiking history. The resulting network representation is analyzed using the Markov Stability framework, a graph theoretical method for community detection across scales, to reveal groups of neurons that are significantly related in the recorded time-series at different levels of granularity.

RESULTS AND COMPARISON WITH EXISTING METHODS

Using synthetic spike-trains, including simulated data from leaky-integrate-and-fire networks, our method is able to identify important patterns in the data such as hierarchical structure that are missed by other standard methods. We further apply the method to experimental data from retinal ganglion cells of mouse and salamander, in which we identify cell-groups that correspond to known functional types, and to hippocampal recordings from rats exploring a linear track, where we detect place cells with high fidelity.

CONCLUSIONS

We present a versatile method to detect neural assemblies in spiking data applicable across a spectrum of relevant scales that contributes to understanding spatio-temporal information gathered from systems neuroscience experiments.

Influence of external inputs and asymmetry of connections on information-geometric measures involving up to ten neuronal interactions

The investigation of neural interactions is crucial for understanding information processing in the brain. Recently an analysis method based on information geometry (IG) has gained increased attention, and the property of the pairwise IG measure has been studied extensively in relation to the two-neuron interaction. However, little is known about the property of IG measures involving more neuronal interactions. In this study, we systematically investigated the influence of external inputs and the asymmetry of connections on the IG measures in cases ranging from 1-neuron to 10-neuron interactions. First, the analytical relationship between the IG measures and external inputs was derived for a network of 10 neurons with uniform connections. Our results confirmed that the single and pairwise IG measures were good estimators of the mean background input and of the sum of the connection weights, respectively. For the IG measures involving 3 to 10 neuronal interactions, we found that the influence of external inputs was highly nonlinear. Second, by computer simulation, we extended our analytical results to asymmetric connections. For a network of 10 neurons, the simulation showed that the behavior of the IG measures in relation to external inputs was similar to the analytical solution obtained for a uniformly connected network. When the network size was increased to 1000 neurons, the influence of external inputs almost disappeared. This result suggests that all IG measures from 1-neuron to 10-neuron interactions are robust against the influence of external inputs. In addition, we investigated how the strength of asymmetry influenced the IG measures. Computer simulation of a 1000-neuron network showed that all the IG measures were robust against the modulation of the asymmetry of connections. Our results provide further support for an information-geometric approach and will provide useful insights when these IG measures are applied to real experimental spike data.

An investigation of Hebbian phase sequences as assembly graphs

Hebb proposed that synapses between neurons that fire synchronously are strengthened, forming cell assemblies and phase sequences. The former, on a shorter scale, are ensembles of synchronized cells that function transiently as a closed processing system; the latter, on a larger scale, correspond to the sequential activation of cell assemblies able to represent percepts and behaviors. Nowadays, the recording of large neuronal populations allows for the detection of multiple cell assemblies. Within Hebb's theory, the next logical step is the analysis of phase sequences. Here we detected phase sequences as consecutive assembly activation patterns, and then analyzed their graph attributes in relation to behavior. We investigated action potentials recorded from the adult rat hippocampus and neocortex before, during and after novel object exploration (experimental periods). Within assembly graphs, each assembly corresponded to a node, and each edge corresponded to the temporal sequence of consecutive node activations. The sum of all assembly activations was proportional to firing rates, but the activity of individual assemblies was not. Assembly repertoire was stable across experimental periods, suggesting that novel experience does not create new assemblies in the adult rat. Assembly graph attributes, on the other hand, varied significantly across behavioral states and experimental periods, and were separable enough to correctly classify experimental periods (Naïve Bayes classifier; maximum AUROCs ranging from 0.55 to 0.99) and behavioral states (waking, slow wave sleep, and rapid eye movement sleep; maximum AUROCs ranging from 0.64 to 0.98). Our findings agree with Hebb's view that assemblies correspond to primitive building blocks of representation, nearly unchanged in the adult, while phase sequences are labile across behavioral states and change after novel experience. The results are compatible with a role for phase sequences in behavior and cognition.

Information-geometric measures estimate neural interactions during oscillatory brain states

The characterization of functional network structures among multiple neurons is essential to understanding neural information processing. Information geometry (IG), a theory developed for investigating a space of probability distributions has recently been applied to spike-train analysis and has provided robust estimations of neural interactions. Although neural firing in the equilibrium state is often assumed in these studies, in reality, neural activity is non-stationary. The brain exhibits various oscillations depending on cognitive demands or when an animal is asleep. Therefore, the investigation of the IG measures during oscillatory network states is important for testing how the IG method can be applied to real neural data. Using model networks of binary neurons or more realistic spiking neurons, we studied how the single- and pairwise-IG measures were influenced by oscillatory neural activity. Two general oscillatory mechanisms, externally driven oscillations and internally induced oscillations, were considered. In both mechanisms, we found that the single-IG measure was linearly related to the magnitude of the external input, and that the pairwise-IG measure was linearly related to the sum of connection strengths between two neurons. We also observed that the pairwise-IG measure was not dependent on the oscillation frequency. These results are consistent with the previous findings that were obtained under the equilibrium conditions. Therefore, we demonstrate that the IG method provides useful insights into neural interactions under the oscillatory condition that can often be observed in the real brain.

Volitional enhancement of firing synchrony and oscillation by neuronal operant conditioning: interaction with neurorehabilitation and brain-machine interface

In this review, we focus on neuronal operant conditioning in which increments in neuronal activities are directly rewarded without behaviors. We discuss the potential of this approach to elucidate neuronal plasticity for enhancing specific brain functions and its interaction with the progress in neurorehabilitation and brain-machine interfaces. The key to-be-conditioned activities that this paper emphasizes are synchronous and oscillatory firings of multiple neurons that reflect activities of cell assemblies. First, we introduce certain well-known studies on neuronal operant conditioning in which conditioned enhancements of neuronal firing were reported in animals and humans. These studies demonstrated the feasibility of volitional control over neuronal activity. Second, we refer to the recent studies on operant conditioning of synchrony and oscillation of neuronal activities. In particular, we introduce a recent study showing volitional enhancement of oscillatory activity in monkey motor cortex and our study showing selective enhancement of firing synchrony of neighboring neurons in rat hippocampus. Third, we discuss the reasons for emphasizing firing synchrony and oscillation in neuronal operant conditioning, the main reason being that they reflect the activities of cell assemblies, which have been suggested to be basic neuronal codes representing information in the brain. Finally, we discuss the interaction of neuronal operant conditioning with neurorehabilitation and brain-machine interface (BMI). We argue that synchrony and oscillation of neuronal firing are the key activities required for developing both reliable neurorehabilitation and high-performance BMI. Further, we conclude that research of neuronal operant conditioning, neurorehabilitation, BMI, and system neuroscience will produce findings applicable to these interrelated fields, and neuronal synchrony and oscillation can be a common important bridge among all of them.

Prediction of rat behavior outcomes in memory tasks using functional connections among neurons

BACKGROUND

Analyzing the neuronal organizational structures and studying the changes in the behavior of the organism is key to understanding cognitive functions of the brain. Although some studies have indicated that spatiotemporal firing patterns of neuronal populations have a certain relationship with the behavioral responses, the issues of whether there are any relationships between the functional networks comprised of these cortical neurons and behavioral tasks and whether it is possible to take advantage of these networks to predict correct and incorrect outcomes of single trials of animals are still unresolved.

METHODOLOGY/PRINCIPAL FINDINGS

This paper presents a new method of analyzing the structures of whole-recorded neuronal functional networks (WNFNs) and local neuronal circuit groups (LNCGs). The activity of these neurons was recorded in several rats. The rats performed two different behavioral tasks, the Y-maze task and the U-maze task. Using the results of the assessment of the WNFNs and LNCGs, this paper describes a realization procedure for predicting the behavioral outcomes of single trials. The methodology consists of four main parts: construction of WNFNs from recorded neuronal spike trains, partitioning the WNFNs into the optimal LNCGs using social community analysis, unsupervised clustering of all trials from each dataset into two different clusters, and predicting the behavioral outcomes of single trials. The results show that WNFNs and LNCGs correlate with the behavior of the animal. The U-maze datasets show higher accuracy for unsupervised clustering results than those from the Y-maze task, and these datasets can be used to predict behavioral responses effectively.

CONCLUSIONS/SIGNIFICANCE

The results of the present study suggest that a methodology proposed in this paper is suitable for analysis of the characteristics of neuronal functional networks and the prediction of rat behavior. These types of structures in cortical ensemble activity may be critical to information representation during the execution of behavior.

Neuronal assembly dynamics in supervised and unsupervised learning scenarios

The dynamic formation of groups of neurons--neuronal assemblies--is believed to mediate cognitive phenomena at many levels, but their detailed operation and mechanisms of interaction are still to be uncovered. One hypothesis suggests that synchronized oscillations underpin their formation and functioning, with a focus on the temporal structure of neuronal signals. In this context, we investigate neuronal assembly dynamics in two complementary scenarios: the first, a supervised spike pattern classification task, in which noisy variations of a collection of spikes have to be correctly labeled; the second, an unsupervised, minimally cognitive evolutionary robotics tasks, in which an evolved agent has to cope with multiple, possibly conflicting, objectives. In both cases, the more traditional dynamical analysis of the system's variables is paired with information-theoretic techniques in order to get a broader picture of the ongoing interactions with and within the network. The neural network model is inspired by the Kuramoto model of coupled phase oscillators and allows one to fine-tune the network synchronization dynamics and assembly configuration. The experiments explore the computational power, redundancy, and generalization capability of neuronal circuits, demonstrating that performance depends nonlinearly on the number of assemblies and neurons in the network and showing that the framework can be exploited to generate minimally cognitive behaviors, with dynamic assembly formation accounting for varying degrees of stimuli modulation of the sensorimotor interactions.

Diverse synchrony of firing reflects diverse cell-assembly coding in the prefrontal cortex

In the present paper, we focus on the coding by cell assemblies in the prefrontal cortex (PFC) and discuss the diversity of the coding, which results in stable and dynamic representations and the processing of various information in that higher brain region. The key activity that reflects cell-assembly coding is the synchrony of the firing of multiple neurons when animals are performing cognitive and memory tasks. First, we introduce some studies that have shown task-related synchrony of neuronal firing in the monkey PFC. These studies have reported fixed and several types of dynamic synchronous firing during working memory, long-term visual memory, and goal selection. The results of these studies have indicated that cell assemblies in the PFC can contribute to both the stability and the dynamics of various types of information. Second, we refer to rat studies and introduce the findings of cellular interactions that contribute to synchrony in working memory, learning-induced changes in synchrony in spatial tasks, and interactions of the PFC and hippocampus in dynamic synchrony. These studies have proposed neuronal mechanisms of cell-assembly coding in the PFC and its critical role in the learning of task demands in problematic situations. Based on the monkey and rat studies, we conclude that cell-assembly coding in the PFC is diverse and has various facets, which allow multipotentiality in the higher brain region. Finally, we discuss the problem of the sizes of cell assembly, how diverse the sizes are in the PFC, and the technical problems in their investigation. We introduce a unique spike-sorting method that can detect small and local cell assemblies that consist of closely neighboring neurons. Then, we describe the findings of our study that showed that the monkey PFC has both small and large cell assemblies, which have different roles in information coding in the working brain.

Canonical microcircuits for predictive coding

This Perspective considers the influential notion of a canonical (cortical) microcircuit in light of recent theories about neuronal processing. Specifically, we conciliate quantitative studies of microcircuitry and the functional logic of neuronal computations. We revisit the established idea that message passing among hierarchical cortical areas implements a form of Bayesian inference-paying careful attention to the implications for intrinsic connections among neuronal populations. By deriving canonical forms for these computations, one can associate specific neuronal populations with specific computational roles. This analysis discloses a remarkable correspondence between the microcircuitry of the cortical column and the connectivity implied by predictive coding. Furthermore, it provides some intuitive insights into the functional asymmetries between feedforward and feedback connections and the characteristic frequencies over which they operate.

Population-wide distributions of neural activity during perceptual decision-making

Cortical activity involves large populations of neurons, even when it is limited to functionally coherent areas. Electrophysiological recordings, on the other hand, involve comparatively small neural ensembles, even when modern-day techniques are used. Here we review results which have started to fill the gap between these two scales of inquiry, by shedding light on the statistical distributions of activity in large populations of cells. We put our main focus on data recorded in awake animals that perform simple decision-making tasks and consider statistical distributions of activity throughout cortex, across sensory, associative, and motor areas. We transversally review the complexity of these distributions, from distributions of firing rates and metrics of spike-train structure, through distributions of tuning to stimuli or actions and of choice signals, and finally the dynamical evolution of neural population activity and the distributions of (pairwise) neural interactions. This approach reveals shared patterns of statistical organization across cortex, including: (i) long-tailed distributions of activity, where quasi-silence seems to be the rule for a majority of neurons; that are barely distinguishable between spontaneous and active states; (ii) distributions of tuning parameters for sensory (and motor) variables, which show an extensive extrapolation and fragmentation of their representations in the periphery; and (iii) population-wide dynamics that reveal rotations of internal representations over time, whose traces can be found both in stimulus-driven and internally generated activity. We discuss how these insights are leading us away from the notion of discrete classes of cells, and are acting as powerful constraints on theories and models of cortical organization and population coding.

Information-geometric measures for estimation of connection weight under correlated inputs

The brain processes information in a highly parallel manner. Determination of the relationship between neural spikes and synaptic connections plays a key role in the analysis of electrophysiological data. Information geometry (IG) has been proposed as a powerful analysis tool for multiple spike data, providing useful insights into the statistical interactions within a population of neurons. Previous work has demonstrated that IG measures can be used to infer the connection weight between two neurons in a neural network. This property is useful in neuroscience because it provides a way to estimate learning-induced changes in synaptic strengths from extracellular neuronal recordings. A previous study has shown, however, that this property would hold only when inputs to neurons are not correlated. Since neurons in the brain often receive common inputs, this would hinder the application of the IG method to real data. We investigated the two-neuron-IG measures in higher-order log-linear models to overcome this limitation. First, we mathematically showed that the estimation of uniformly connected synaptic weight can be improved by taking into account higher-order log-linear models. Second, we numerically showed that the estimation can be improved for more general asymmetrically connected networks. Considering the estimated number of the synaptic connections in the brain, we showed that the two-neuron IG measure calculated by the fourth- or fifth-order log-linear model would provide an accurate estimation of connection strength within approximately a 10% error. These studies suggest that the two-neuron IG measure with higher-order log-linear expansion is a robust estimator of connection weight even under correlated inputs, providing a useful analytical tool for real multineuronal spike data.

Sodium concentration coding gives way to evaluative coding in cortex and amygdala

Typically, stimulus batteries used to characterize sensory neural coding span physical parameter spaces (e.g., concentration: from low to high). For awake animals, however, psychological variables (e.g., pleasantness/palatability) with complicated relationships to the physical often dominate neural responses. Here we pit physical and psychological axes against one another, presenting awake rats with a stimulus set including 4 NaCl concentrations (0.01, 0.1, 0.3, and 1.0 m) plus palatable (0.3 m sucrose) and aversive (0.001 m quinine) benchmarks, while recording the activity of neurons in two sites vital for NaCl taste processing, gustatory cortex (GC) and central amygdala (CeA). Since NaCl palatability (i.e., preference) follows a non-monotonic, "inverted-U-shaped" curve while concentration increases monotonically, this stimulus battery allowed us to test whether GC and CeA responses better reflect external or internal variables. As predicted, GC single-neuron and population responses reflected both parameters in separate response epochs: sodium concentration-related information appeared with the earliest taste-specific responses, giving way to palatability-related information, in an overlapping subset of neurons, several hundred milliseconds later. CeA single-neuron and population responses, meanwhile, contained only a brief period of concentration specificity, occurring just before palatability-related information emerged (simultaneously with, or slightly later than, in GC). Thus, cortex and amygdala both prominently reflect NaCl palatability late in their responses; CeA neurons largely respond to either palatable or aversive stimuli, while GC responses tend to reflect the entire palatability spectrum in a graded fashion.

Sodium concentration coding gives way to evaluative coding in cortex and amygdala

Typically, stimulus batteries used to characterize sensory neural coding span physical parameter spaces (e.g., concentration: from low to high). For awake animals, however, psychological variables (e.g., pleasantness/palatability) with complicated relationships to the physical often dominate neural responses. Here we pit physical and psychological axes against one another, presenting awake rats with a stimulus set including 4 NaCl concentrations (0.01, 0.1, 0.3, and 1.0 m) plus palatable (0.3 m sucrose) and aversive (0.001 m quinine) benchmarks, while recording the activity of neurons in two sites vital for NaCl taste processing, gustatory cortex (GC) and central amygdala (CeA). Since NaCl palatability (i.e., preference) follows a non-monotonic, "inverted-U-shaped" curve while concentration increases monotonically, this stimulus battery allowed us to test whether GC and CeA responses better reflect external or internal variables. As predicted, GC single-neuron and population responses reflected both parameters in separate response epochs: sodium concentration-related information appeared with the earliest taste-specific responses, giving way to palatability-related information, in an overlapping subset of neurons, several hundred milliseconds later. CeA single-neuron and population responses, meanwhile, contained only a brief period of concentration specificity, occurring just before palatability-related information emerged (simultaneously with, or slightly later than, in GC). Thus, cortex and amygdala both prominently reflect NaCl palatability late in their responses; CeA neurons largely respond to either palatable or aversive stimuli, while GC responses tend to reflect the entire palatability spectrum in a graded fashion.

Modeling and analyzing higher-order correlations in non-Poissonian spike trains

Measuring pairwise and higher-order spike correlations is crucial for studying their potential impact on neuronal information processing. In order to avoid misinterpretation of results, the tools used for data analysis need to be carefully calibrated with respect to their sensitivity and robustness. This, in turn, requires surrogate data with statistical properties common to experimental spike trains. Here, we present a novel method to generate correlated non-Poissonian spike trains and study the impact of single-neuron spike statistics on the inference of higher-order correlations. Our method to mimic cooperative neuronal spike activity allows the realization of a large variety of renewal processes with controlled higher-order correlation structure. Based on surrogate data obtained by this procedure we investigate the robustness of the recently proposed method empirical de-Poissonization (Ehm et al., 2007). It assumes Poissonian spiking, which is common also for many other estimation techniques. We observe that some degree of deviation from this assumption can generally be tolerated, that the results are more reliable for small analysis bins, and that the degree of misestimation depends on the detailed spike statistics. As a consequence of these findings we finally propose a strategy to assess the reliability of results for experimental data.