The effects of pooling on spike train correlations

Between-area communication through the lens of within-area neuronal dynamics

A core problem in systems and circuits neuroscience is deciphering the origin of shared dynamics in neuronal activity: Do they emerge through local network interactions, or are they inherited from external sources? We explore this question with large-scale networks of spatially ordered spiking neuron models where a downstream network receives input from an upstream sender network. We show that linear measures of the communication between the sender and receiver networks can discriminate between emergent or inherited population dynamics. A match in the dimensionality of the sender and receiver population activities promotes faithful communication. In contrast, a nonlinear mapping between the sender to receiver activity, for example, through downstream emergent population-wide fluctuations, can impair linear communication. Our work exposes the benefits and limitations of linear measures when analyzing between-area communication in circuits with rich population-wide neuronal dynamics.

Reconciliation of weak pairwise spike-train correlations and highly coherent local field potentials across space

Multi-electrode arrays covering several square millimeters of neural tissue provide simultaneous access to population signals such as extracellular potentials and spiking activity of one hundred or more individual neurons. The interpretation of the recorded data calls for multiscale computational models with corresponding spatial dimensions and signal predictions. Multi-layer spiking neuron network models of local cortical circuits covering about $1\,{\text{mm}^{2}}$ have been developed, integrating experimentally obtained neuron-type-specific connectivity data and reproducing features of observed in-vivo spiking statistics. Local field potentials can be computed from the simulated spiking activity. We here extend a local network and local field potential model to an area of $4\times 4\,{\text{mm}^{2}}$, preserving the neuron density and introducing distance-dependent connection probabilities and conduction delays. We find that the upscaling procedure preserves the overall spiking statistics of the original model and reproduces asynchronous irregular spiking across populations and weak pairwise spike-train correlations in agreement with experimental recordings from sensory cortex. Also compatible with experimental observations, the correlation of local field potential signals is strong and decays over a distance of several hundred micrometers. Enhanced spatial coherence in the low-gamma band around $50\,\text{Hz}$ may explain the recent report of an apparent band-pass filter effect in the spatial reach of the local field potential.

Topographic Organization of Correlation Along the Longitudinal and Transverse Axes in Rat Hippocampal CA3 Due to Excitatory Afferents

The topographic organization of afferents to the hippocampal CA3 subfield are well-studied, but their role in influencing the spatiotemporal dynamics of population activity is not understood. Using a large-scale, computational neuronal network model of the entorhinal-dentate-CA3 system, the effects of the perforant path, mossy fibers, and associational system on the propagation and transformation of network spiking patterns were investigated. A correlation map was constructed to characterize the spatial structure and temporal evolution of pairwise correlations which underlie the emergent patterns found in the population activity. The topographic organization of the associational system gave rise to changes in the spatial correlation structure along the longitudinal and transverse axes of the CA3. The resulting gradients may provide a basis for the known functional organization observed in hippocampus.

Interneuronal correlations at longer time scales predict decision signals for bistable structure-from-motion perception

Perceptual decisions are thought to depend on the activation of task-relevant neurons, whose activity is often correlated in time. Here, we examined how the temporal structure of shared variability in neuronal firing relates to perceptual choices. We recorded stimulus-selective neurons from visual area V5/MT while two monkeys (Macaca mulatta) made perceptual decisions about the rotation direction of structure-from-motion cylinders. Interneuronal correlations for a perceptually ambiguous cylinder stimulus were significantly higher than those for unambiguous cylinders or for random 2D motion during passive viewing. Much of the difference arose from correlations at relatively long timescales (hundreds of milliseconds). Choice-related neural activity (quantified as choice probability; CP) for ambiguous cylinders was positively correlated with interneuronal correlations and was specifically associated with their long timescale component. Furthermore, the slope of the long timescale - but not the instantaneous - component of the correlation predicted higher CPs towards the end of the trial i.e. close to the decision. Our results suggest that the perceptual stability of structure-from-motion cylinders may be controlled by enhanced interneuronal correlations on longer timescales. We propose this as a potential signature of top-down influences onto V5/MT processing that shape and stabilize the appearance of 3D-motion percepts.

Correlated states in balanced neuronal networks

Understanding the magnitude and structure of interneuronal correlations and their relationship to synaptic connectivity structure is an important and difficult problem in computational neuroscience. Early studies show that neuronal network models with excitatory-inhibitory balance naturally create very weak spike train correlations, defining the "asynchronous state." Later work showed that, under some connectivity structures, balanced networks can produce larger correlations between some neuron pairs, even when the average correlation is very small. All of these previous studies assume that the local network receives feedforward synaptic input from a population of uncorrelated spike trains. We show that when spike trains providing feedforward input are correlated, the downstream recurrent network produces much larger correlations. We provide an in-depth analysis of the resulting "correlated state" in balanced networks and show that, unlike the asynchronous state, it produces a tight excitatory-inhibitory balance consistent with in vivo cortical recordings.

Variations in Static Force Control and Motor Unit Behavior with Error Amplification Feedback in the Elderly

Error amplification (EA) feedback is a promising approach to advance visuomotor skill. As error detection and visuomotor processing at short time scales decline with age, this study examined whether older adults could benefit from EA feedback that included higher-frequency information to guide a force-tracking task. Fourteen young and 14 older adults performed low-level static isometric force-tracking with visual guidance of typical visual feedback and EA feedback containing augmented high-frequency errors. Stabilogram diffusion analysis was used to characterize force fluctuation dynamics. Also, the discharge behaviors of motor units and pooled motor unit coherence were assessed following the decomposition of multi-channel surface electromyography (EMG). EA produced different behavioral and neurophysiological impacts on young and older adults. Older adults exhibited inferior task accuracy with EA feedback than with typical visual feedback, but not young adults. Although stabilogram diffusion analysis revealed that EA led to a significant decrease in critical time points for both groups, EA potentiated the critical point of force fluctuations [Formula: see text], short-term effective diffusion coefficients (Ds), and short-term exponent scaling only for the older adults. Moreover, in older adults, EA added to the size of discharge variability of motor units and discharge regularity of cumulative discharge rate, but suppressed the pooled motor unit coherence in the 13-35 Hz band. Virtual EA alters the strategic balance between open-loop and closed-loop controls for force-tracking. Contrary to expectations, the prevailing use of closed-loop control with EA that contained high-frequency error information enhanced the motor unit discharge variability and undermined the force steadiness in the older group, concerning declines in physiological complexity in the neurobehavioral system and the common drive to the motoneuronal pool against force destabilization.

Altered Motor Unit Discharge Coherence in Paretic Muscles of Stroke Survivors

After a cerebral stroke, a series of changes at the supraspinal and spinal nervous system can alter the control of muscle activation, leading to persistent motor impairment. However, the relative contribution of these different levels of the nervous system to impaired muscle activation is not well understood. The coherence of motor unit (MU) spike trains is considered to partly reflect activities of higher level control, with different frequency band representing different levels of control. Accordingly, the objective of this study was to quantify the different sources of contribution to altered muscle activation. We examined the coherence of MU spike trains decomposed from surface electromyogram (sEMG) of the first dorsal interosseous muscle on both paretic and contralateral sides of 14 hemispheric stroke survivors. sEMG was obtained over a range of force contraction levels at 40, 50, and 60% of maximum voluntary contraction. Our results showed that MU coherence increased significantly in delta (1–4 Hz), alpha (8–12 Hz), and beta (15–30 Hz) bands on the affected side compared with the contralateral side, but was maintained at the same level in the gamma (30–60 Hz) band. In addition, no significant alteration was observed across medium–high force levels (40–60%). These results indicated that the common synaptic input to motor neurons increased on the paretic side, and the increased common input can originate from changes at multiple levels, including spinal and supraspinal levels following a stroke. All these changes can contribute to impaired activation of affected muscles in stroke survivors. Our findings also provide evidence regarding the different origins of impaired muscle activation poststroke.

Alterations in Neural Control of Constant Isometric Contraction with the Size of Error Feedback

Discharge patterns from a population of motor units (MUs) were estimated with multi-channel surface electromyogram and signal processing techniques to investigate parametric differences in low-frequency force fluctuations, MU discharges, and force-discharge relation during static force-tracking with varying sizes of execution error presented via visual feedback. Fourteen healthy adults produced isometric force at 10% of maximal voluntary contraction through index abduction under three visual conditions that scaled execution errors with different amplification factors. Error-augmentation feedback that used a high amplification factor (HAF) to potentiate visualized error size resulted in higher sample entropy, mean frequency, ratio of high-frequency components, and spectral dispersion of force fluctuations than those of error-reducing feedback using a low amplification factor (LAF). In the HAF condition, MUs with relatively high recruitment thresholds in the dorsal interosseous muscle exhibited a larger coefficient of variation for inter-spike intervals and a greater spectral peak of the pooled MU coherence at 13-35 Hz than did those in the LAF condition. Manipulation of the size of error feedback altered the force-discharge relation, which was characterized with non-linear approaches such as mutual information and cross sample entropy. The association of force fluctuations and global discharge trace decreased with increasing error amplification factor. Our findings provide direct neurophysiological evidence that favors motor training using error-augmentation feedback. Amplification of the visualized error size of visual feedback could enrich force gradation strategies during static force-tracking, pertaining to selective increases in the discharge variability of higher-threshold MUs that receive greater common oscillatory inputs in the β-band.

The mechanics of state-dependent neural correlations

Simultaneous recordings from large neural populations are becoming increasingly common. An important feature of population activity is the trial-to-trial correlated fluctuation of spike train outputs from recorded neuron pairs. Similar to the firing rate of single neurons, correlated activity can be modulated by a number of factors, from changes in arousal and attentional state to learning and task engagement. However, the physiological mechanisms that underlie these changes are not fully understood. We review recent theoretical results that identify three separate mechanisms that modulate spike train correlations: changes in input correlations, internal fluctuations and the transfer function of single neurons. We first examine these mechanisms in feedforward pathways and then show how the same approach can explain the modulation of correlations in recurrent networks. Such mechanistic constraints on the modulation of population activity will be important in statistical analyses of high-dimensional neural data.

The proportion of common synaptic input to motor neurons increases with an increase in net excitatory input

α-Motor neurons receive synaptic inputs from spinal and supraspinal centers that comprise components either common to the motor neuron pool or independent. The input shared by motor neurons--common input--determines force control. The aim of the study was to investigate the changes in the strength of common synaptic input delivered to motor neurons with changes in force and with fatigue, two conditions that underlie an increase in the net excitatory drive to the motor neurons. High-density surface electromyogram (EMG) signals were recorded from the tibialis anterior muscle during contractions at 20, 50, and 75% of the maximal voluntary contraction force (in 3 sessions separated by at least 2 days), all sustained until task failure. EMG signal decomposition identified the activity of a total of 1,245 motor units. The coherence values between cumulative motor unit spike trains increased with increasing force, especially for low frequencies. This increase in coherence was not observed when comparing two subsets of motor units having different recruitment thresholds, but detected at the same force level. Moreover, the coherence values for frequencies <5 Hz increased at task failure with respect to the beginning of the contractions for all force levels. In conclusion, the results indicated that the relative strength of common synaptic input to motor neurons increases with respect to independent input when the net excitatory drive to motor neurons increases as a consequence of a change in force and fatigue.

Role of input correlations in shaping the variability and noise correlations of evoked activity in the neocortex

Recent analysis of evoked activity recorded across different brain regions and tasks revealed a marked decrease in noise correlations and trial-by-trial variability. Given the importance of correlations and variability for information processing within the rate coding paradigm, several mechanisms have been proposed to explain the reduction in these quantities despite an increase in firing rates. These models suggest that anatomical clusters and/or tightly balanced excitation-inhibition can generate intrinsic network dynamics that may exhibit a reduction in noise correlations and trial-by-trial variability when perturbed by an external input. Such mechanisms based on the recurrent feedback crucially ignore the contribution of feedforward input to the statistics of the evoked activity. Therefore, we investigated how statistical properties of the feedforward input shape the statistics of the evoked activity. Specifically, we focused on the effect of input correlation structure on the noise correlations and trial-by-trial variability. We show that the ability of neurons to transfer the input firing rate, correlation, and variability to the output depends on the correlations within the presynaptic pool of a neuron, and that an input with even weak within-correlations can be sufficient to reduce noise correlations and trial-by-trial variability, without requiring any specific recurrent connectivity structure. In general, depending on the ongoing activity state, feedforward input could either increase or decrease noise correlation and trial-by-trial variability. Thus, we propose that evoked activity statistics are jointly determined by the feedforward and feedback inputs.

Oscillation-induced signal transmission and gating in neural circuits

Reliable signal transmission constitutes a key requirement for neural circuit function. The propagation of synchronous pulse packets through recurrent circuits is hypothesized to be one robust form of signal transmission and has been extensively studied in computational and theoretical works. Yet, although external or internally generated oscillations are ubiquitous across neural systems, their influence on such signal propagation is unclear. Here we systematically investigate the impact of oscillations on propagating synchrony. We find that for standard, additive couplings and a net excitatory effect of oscillations, robust propagation of synchrony is enabled in less prominent feed-forward structures than in systems without oscillations. In the presence of non-additive coupling (as mediated by fast dendritic spikes), even balanced oscillatory inputs may enable robust propagation. Here, emerging resonances create complex locking patterns between oscillations and spike synchrony. Interestingly, these resonances make the circuits capable of selecting specific pathways for signal transmission. Oscillations may thus promote reliable transmission and, in co-action with dendritic nonlinearities, provide a mechanism for information processing by selectively gating and routing of signals. Our results are of particular interest for the interpretation of sharp wave/ripple complexes in the hippocampus, where previously learned spike patterns are replayed in conjunction with global high-frequency oscillations. We suggest that the oscillations may serve to stabilize the replay.

Rhythmic activity in the brain is ubiquitous, its functions are debated. Here we show that it may contribute to the reliable transmission of information within brain areas. We find that its effect is particularly strong if we take nonlinear coupling into account. This experimentally found neuronal property implies that inputs which arrive nearly simultaneously can have a much stronger impact than expected from the sum of their individuals strengths. In such systems, rhythmic activity supports information transmission even if its positive and negative part exactly cancels all the time. Further, the information transmission can adapt to the oscillation frequency to optimally benefit from it. Finally, we show that rhythms with different frequencies may enable or disable communication channels, and are thus suitable for the steering of information flow.

Physiology and Impact of Horizontal Connections in Rat Neocortex

Cortical information processing at the cellular level has predominantly been studied in local networks, which are dominated by strong vertical connectivity between layers. However, recent studies suggest that the bulk of axons targeting pyramidal neurons most likely originate from outside this local range, emphasizing the importance of horizontal connections. We mapped a subset of these connections to L5B pyramidal neurons in rat somatosensory cortex with photostimulation, identifying intact projections up to a lateral distance of 2 mm. Our estimates of the spatial distribution of cells presynaptic to L5B pyramids support the idea that the majority is located outside the local volume. The synaptic physiology of horizontal connections does not differ markedly from that of local connections, whereas the layer and cell-type-dependent pattern of innervation does. Apart from L2/3, L6A provides a strong source of horizontal connections. Implementing our data into a spiking neuronal network model shows that more horizontal connections promote robust asynchronous ongoing activity states and reduce noise correlations in stimulus-induced activity.

Propagating synchrony in feed-forward networks

Coordinated patterns of precisely timed action potentials (spikes) emerge in a variety of neural circuits but their dynamical origin is still not well understood. One hypothesis states that synchronous activity propagating through feed-forward chains of groups of neurons (synfire chains) may dynamically generate such spike patterns. Additionally, synfire chains offer the possibility to enable reliable signal transmission. So far, mostly densely connected chains, often with all-to-all connectivity between groups, have been theoretically and computationally studied. Yet, such prominent feed-forward structures have not been observed experimentally. Here we analytically and numerically investigate under which conditions diluted feed-forward chains may exhibit synchrony propagation. In addition to conventional linear input summation, we study the impact of non-linear, non-additive summation accounting for the effect of fast dendritic spikes. The non-linearities promote synchronous inputs to generate precisely timed spikes. We identify how non-additive coupling relaxes the conditions on connectivity such that it enables synchrony propagation at connectivities substantially lower than required for linearly coupled chains. Although the analytical treatment is based on a simple leaky integrate-and-fire neuron model, we show how to generalize our methods to biologically more detailed neuron models and verify our results by numerical simulations with, e.g., Hodgkin Huxley type neurons.

Correlations in background activity control persistent state stability and allow execution of working memory tasks

Working memory (WM) requires selective information gating, active information maintenance, and rapid active updating. Hence performing a WM task needs rapid and controlled transitions between neural persistent activity and the resting state. We propose that changes in correlations in neural activity provides a mechanism for the required WM operations. As a proof of principle, we implement sustained activity and WM in recurrently coupled spiking networks with neurons receiving excitatory random background activity where background correlations are induced by a common noise source. We first characterize how the level of background correlations controls the stability of the persistent state. With sufficiently high correlations, the sustained state becomes practically unstable, so it cannot be initiated by a transient stimulus. We exploit this in WM models implementing the delay match to sample task by modulating flexibly in time the correlation level at different phases of the task. The modulation sets the network in different working regimes: more prompt to gate in a signal or clear the memory. We examine how the correlations affect the ability of the network to perform the task when distractors are present. We show that in a winner-take-all version of the model, where two populations cross-inhibit, correlations make the distractor blocking robust. In a version of the mode where no cross inhibition is present, we show that appropriate modulation of correlation levels is sufficient to also block the distractor access while leaving the relevant memory trace in tact. The findings presented in this manuscript can form the basis for a new paradigm about how correlations are flexibly controlled by the cortical circuits to execute WM operations.

A generative spike train model with time-structured higher order correlations

Emerging technologies are revealing the spiking activity in ever larger neural ensembles. Frequently, this spiking is far from independent, with correlations in the spike times of different cells. Understanding how such correlations impact the dynamics and function of neural ensembles remains an important open problem. Here we describe a new, generative model for correlated spike trains that can exhibit many of the features observed in data. Extending prior work in mathematical finance, this generalized thinning and shift (GTaS) model creates marginally Poisson spike trains with diverse temporal correlation structures. We give several examples which highlight the model's flexibility and utility. For instance, we use it to examine how a neural network responds to highly structured patterns of inputs. We then show that the GTaS model is analytically tractable, and derive cumulant densities of all orders in terms of model parameters. The GTaS framework can therefore be an important tool in the experimental and theoretical exploration of neural dynamics.

Short-term synaptic depression and stochastic vesicle dynamics reduce and shape neuronal correlations

Correlated neuronal activity is an important feature in many neural codes, a neural correlate of a variety of cognitive states, as well as a signature of several disease states in the nervous system. The cellular and circuit mechanics of neural correlations is a vibrant area of research. Synapses throughout the cortex exhibit a form of short-term depression where increased presynaptic firing rates deplete neurotransmitter vesicles, which transiently reduces synaptic efficacy. The release and recovery of these vesicles are inherently stochastic, and this stochasticity introduces variability into the conductance elicited by depressing synapses. The impact of spiking and subthreshold membrane dynamics on the transfer of neuronal correlations has been studied intensively, but an investigation of the impact of short-term synaptic depression and stochastic vesicle dynamics on correlation transfer is lacking. We find that short-term synaptic depression and stochastic vesicle dynamics can substantially reduce correlations, shape the timescale over which these correlations occur, and alter the dependence of spiking correlations on firing rate. Our results show that short-term depression and stochastic vesicle dynamics need to be taken into account when modeling correlations in neuronal populations.

Factors influencing the estimates of correlation between motor unit activities in humans

Background

Alpha motoneurons receive common synaptic inputs from spinal and supraspinal pathways. As a result, a certain degree of correlation can be observed between motoneuron spike trains during voluntary contractions. This has been studied by using correlation measures in the time and frequency domains. These measures are interpreted as reflecting different types of connectivity in the spinal networks, although the relation between the degree of correlation of the output motoneuron spike trains and of their synaptic inputs is unclear.

Methodology/Principal Findings

In this study, we analyze theoretically this relation and we complete this analysis by simulations and experimental data on the abductor digiti minimi muscle. The results demonstrate that correlation measures between motoneuron output spike trains are inherently influenced by the discharge rate and that this influence cannot be compensated by normalization. Because of the influence of discharge rate, frequency domain measures of correlation (coherence) do not identify the full frequency content of the common input signal when computed from pairs of motoneurons. Rather, an increase in sampling rate is needed by using cumulative spike trains of several motoneurons. Moreover, the application of averaging filters to the spike trains influences the magnitude of the estimated correlation levels calculated in the time, but not in the frequency domain (coherence).

Conclusions

It is concluded that the analysis of coherence in different frequency bands between cumulative spike trains of a sufficient number of motoneurons provides information on the spectrum of the common synaptic input. Nonetheless, the absolute values of coherent peaks cannot be compared across conditions with different cumulative discharge rates.

Correlated neural variability in persistent state networks

Neural activity that persists long after stimulus presentation is a biological correlate of short-term memory. Variability in spiking activity causes persistent states to drift over time, ultimately degrading memory. Models of short-term memory often assume that the input fluctuations to neural populations are independent across cells, a feature that attenuates population-level variability and stabilizes persistent activity. However, this assumption is at odds with experimental recordings from pairs of cortical neurons showing that both the input currents and output spike trains are correlated. It remains unclear how correlated variability affects the stability of persistent activity and the performance of cognitive tasks that it supports. We consider the stochastic long-timescale attractor dynamics of pairs of mutually inhibitory populations of spiking neurons. In these networks, persistent activity was less variable when correlated variability was globally distributed across both populations compared with the case when correlations were locally distributed only within each population. Using a reduced firing rate model with a continuum of persistent states, we show that, when input fluctuations are correlated across both populations, they drive firing rate fluctuations orthogonal to the persistent state attractor, thereby causing minimal stochastic drift. Using these insights, we establish that distributing correlated fluctuations globally as opposed to locally improves network's performance on a two-interval, delayed response discrimination task. Our work shows that the correlation structure of input fluctuations to a network is an important factor when determining long-timescale, persistent population spiking activity.

Sensitivity of noisy neurons to coincident inputs

How do neurons compute? Two main theories compete: neurons could temporally integrate noisy inputs (rate-based theories) or they could detect coincident input spikes (spike timing-based theories). Correlations at fine timescales have been observed in many areas of the nervous system, but they might have a minor impact. To address this issue, we used a probabilistic approach to quantify the impact of coincidences on neuronal response in the presence of fluctuating synaptic activity. We found that when excitation and inhibition are balanced, as in the sensory cortex in vivo, synchrony in a very small proportion of inputs results in dramatic increases in output firing rate. Our theory was experimentally validated with in vitro recordings of cortical neurons of mice. We conclude that not only are noisy neurons well equipped to detect coincidences, but they are so sensitive to fine correlations that a rate-based description of neural computation is unlikely to be accurate in general.

Membrane potential and spike train statistics depend distinctly on input statistics

A description of how the activity of a population of neurons reflects the structure of its inputs is essential for understanding neural coding. Many studies have examined how inputs determine spiking statistics, while comparatively little is known about membrane potentials. We examine how membrane potential statistics are related to input and spiking statistics. Surprisingly, firing rates and membrane potentials are sensitive to input current modulations in distinct regimes. Additionally, the correlation between the membrane potentials of two uncoupled cells and the correlation between their spike trains reflect input correlations in distinct regimes. Our predictions are experimentally testable, provide insight into the filtering properties of neurons, and indicate that care needs to be taken when interpreting neuronal recordings that reflect a combination of subthreshold and spiking activity.