Spike-Centered Jitter Can Mistake Temporal Structure

The olivary input to the cerebellum dissociates sensory events from movement plans

Neurons in the inferior olive are thought to anatomically organize the Purkinje cells (P-cells) of the cerebellum into computational modules, but what is computed by each module? Here, we designed a saccade task in marmosets that dissociated sensory events from motor events and then recorded the complex and simple spikes of hundreds of P-cells. We found that when a visual target was presented at a random location, the olive reported the direction of that sensory event to one group of P-cells, but not to a second group. However, just before movement onset, it reported the direction of the planned movement to both groups, even if that movement was not toward the target. At the end of the movement if the subject experienced an error but chose to withhold the corrective movement, only the first group received information about the sensory prediction error. We organized the P-cells based on the information content of their olivary input and found that in the group that received sensory information, the simple spikes were suppressed during fixation, then produced a burst before saccade onset in a direction consistent with assisting the movement. In the second group, the simple spikes were not suppressed during fixation but burst near saccade deceleration in a direction consistent with stopping the movement. Thus, the olive differentiated the P-cells based on whether they would receive sensory or motor information, and this defined their contributions to control of movements as well as holding still.

Uncovering the organization of neural circuits with Generalized Phase Locking Analysis

Despite the considerable progress of in vivo neural recording techniques, inferring the biophysical mechanisms underlying large scale coordination of brain activity from neural data remains challenging. One obstacle is the difficulty to link high dimensional functional connectivity measures to mechanistic models of network activity. We address this issue by investigating spike-field coupling (SFC) measurements, which quantify the synchronization between, on the one hand, the action potentials produced by neurons, and on the other hand mesoscopic "field" signals, reflecting subthreshold activities at possibly multiple recording sites. As the number of recording sites gets large, the amount of pairwise SFC measurements becomes overwhelmingly challenging to interpret. We develop Generalized Phase Locking Analysis (GPLA) as an interpretable dimensionality reduction of this multivariate SFC. GPLA describes the dominant coupling between field activity and neural ensembles across space and frequencies. We show that GPLA features are biophysically interpretable when used in conjunction with appropriate network models, such that we can identify the influence of underlying circuit properties on these features. We demonstrate the statistical benefits and interpretability of this approach in various computational models and Utah array recordings. The results suggest that GPLA, used jointly with biophysical modeling, can help uncover the contribution of recurrent microcircuits to the spatio-temporal dynamics observed in multi-channel experimental recordings.

Error correction and improved precision of spike timing in converging cortical networks

The brain propagates neuronal signals accurately and rapidly. Nevertheless, whether and how a pool of cortical neurons transmits an undistorted message to a target remains unclear. We apply optogenetic white noise signals to small assemblies of cortical pyramidal cells (PYRs) in freely moving mice. The directly activated PYRs exhibit a spike timing precision of several milliseconds. Instead of losing precision, interneurons driven via synaptic activation exhibit higher precision with respect to the white noise signal. Compared with directly activated PYRs, postsynaptic interneuron spike trains allow better signal reconstruction, demonstrating error correction. Data-driven modeling shows that nonlinear amplification of coincident spikes can generate error correction and improved precision. Over multiple applications of the same signal, postsynaptic interneuron spiking is most reliable at timescales ten times shorter than those of the presynaptic PYR, exhibiting temporal coding. Similar results are observed in hippocampal region CA1. Coincidence detection of convergent inputs enables messages to be precisely propagated between cortical PYRs and interneurons.

•

PYR-to-interneuron spike transmission exhibits error correction and improved precision

•

Interneuron precision is higher when a larger pool of presynaptic PYRs is recruited

•

Error correction and improved precision are consistent with coincidence detection

•

Interneurons activated by synaptic transmission act as temporal coders

Deconvolution improves the detection and quantification of spike transmission gain from spike trains

Accurate detection and quantification of spike transmission between neurons is essential for determining neural network mechanisms that govern cognitive functions. Using point process and conductance-based simulations, we found that existing methods for determining neuronal connectivity from spike times are highly affected by burst spiking activity, resulting in over- or underestimation of spike transmission. To improve performance, we developed a mathematical framework for decomposing the cross-correlation between two spike trains. We then devised a deconvolution-based algorithm for removing effects of second-order spike train statistics. Deconvolution removed the effect of burst spiking, improving the estimation of neuronal connectivity yielded by state-of-the-art methods. Application of deconvolution to neuronal data recorded from hippocampal region CA1 of freely-moving mice produced higher estimates of spike transmission, in particular when spike trains exhibited bursts. Deconvolution facilitates the precise construction of complex connectivity maps, opening the door to enhanced understanding of the neural mechanisms underlying brain function.

Monosynaptic inference via finely-timed spikes

Observations of finely-timed spike relationships in population recordings have been used to support partial reconstruction of neural microcircuit diagrams. In this approach, fine-timescale components of paired spike train interactions are isolated and subsequently attributed to synaptic parameters. Recent perturbation studies strengthen the case for such an inference, yet the complete set of measurements needed to calibrate statistical models is unavailable. To address this gap, we study features of pairwise spiking in a large-scale in vivo dataset where presynaptic neurons were explicitly decoupled from network activity by juxtacellular stimulation. We then construct biophysical models of paired spike trains to reproduce the observed phenomenology of in vivo monosynaptic interactions, including both fine-timescale spike-spike correlations and firing irregularity. A key characteristic of these models is that the paired neurons are coupled by rapidly-fluctuating background inputs. We quantify a monosynapse's causal effect by comparing the postsynaptic train with its counterfactual, when the monosynapse is removed. Subsequently, we develop statistical techniques for estimating this causal effect from the pre- and post-synaptic spike trains. A particular focus is the justification and application of a nonparametric separation of timescale principle to implement synaptic inference. Using simulated data generated from the biophysical models, we characterize the regimes in which the estimators accurately identify the monosynaptic effect. A secondary goal is to initiate a critical exploration of neurostatistical assumptions in terms of biophysical mechanisms, particularly with regards to the challenging but arguably fundamental issue of fast, unobservable nonstationarities in background dynamics.

Decoupling of timescales reveals sparse convergent CPG network in the adult spinal cord

During the generation of rhythmic movements, most spinal neurons receive an oscillatory synaptic drive. The neuronal architecture underlying this drive is unknown, and the corresponding network size and sparseness have not yet been addressed. If the input originates from a small central pattern generator (CPG) with dense divergent connectivity, it will induce correlated input to all receiving neurons, while sparse convergent wiring will induce a weak correlation, if any. Here, we use pairwise recordings of spinal neurons to measure synaptic correlations and thus infer the wiring architecture qualitatively. A strong correlation on a slow timescale implies functional relatedness and a common source, which will also cause correlation on fast timescale due to shared synaptic connections. However, we consistently find marginal coupling between slow and fast correlations regardless of neuronal identity. This suggests either sparse convergent connectivity or a CPG network with recurrent inhibition that actively decorrelates common input.

Methods for identification of spike patterns in massively parallel spike trains

Temporally, precise correlations between simultaneously recorded neurons have been interpreted as signatures of cell assemblies, i.e., groups of neurons that form processing units. Evidence for this hypothesis was found on the level of pairwise correlations in simultaneous recordings of few neurons. Increasing the number of simultaneously recorded neurons increases the chances to detect cell assembly activity due to the larger sample size. Recent technological advances have enabled the recording of 100 or more neurons in parallel. However, these massively parallel spike train data require novel statistical tools to be analyzed for correlations, because they raise considerable combinatorial and multiple testing issues. Recently, various of such methods have started to develop. First approaches were based on population or pairwise measures of synchronization, and later led to methods for the detection of various types of higher-order synchronization and of spatio-temporal patterns. The latest techniques combine data mining with analysis of statistical significance. Here, we give a comparative overview of these methods, of their assumptions and of the types of correlations they can detect.

Computational Neuroscience: Mathematical and Statistical Perspectives

Mathematical and statistical models have played important roles in neuroscience, especially by describing the electrical activity of neurons recorded individually, or collectively across large networks. As the field moves forward rapidly, new challenges are emerging. For maximal effectiveness, those working to advance computational neuroscience will need to appreciate and exploit the complementary strengths of mechanistic theory and the statistical paradigm.

Pyramidal Cell-Interneuron Circuit Architecture and Dynamics in Hippocampal Networks

Excitatory control of inhibitory neurons is poorly understood due to the difficulty of studying synaptic connectivity in vivo. We inferred such connectivity through analysis of spike timing and validated this inference using juxtacellular and optogenetic control of presynaptic spikes in behaving mice. We observed that neighboring CA1 neurons had stronger connections and that superficial pyramidal cells projected more to deep interneurons. Connection probability and strength were skewed, with a minority of highly connected hubs. Divergent presynaptic connections led to synchrony between interneurons. Synchrony of convergent presynaptic inputs boosted postsynaptic drive. Presynaptic firing frequency was read out by postsynaptic neurons through short-term depression and facilitation, with individual pyramidal cells and interneurons displaying a diversity of spike transmission filters. Additionally, spike transmission was strongly modulated by prior spike timing of the postsynaptic cell. These results bridge anatomical structure with physiological function.