Movement is governed by rotational neural dynamics in spinal motor networks.

Although the generation of movements is a fundamental function of the nervous system, the underlying neural principles remain unclear. As flexor and extensor muscle activities alternate during rhythmic movements such as walking, it is often assumed that the responsible neural circuitry is similarly exhibiting alternating activity1. Here we present ensemble recordings of neurons in the lumbar spinal cord that indicate that, rather than alternating, the population is performing a low-dimensional 'rotation' in neural space, in which the neural activity is cycling through all phases continuously during the rhythmic behaviour. The radius of rotation correlates with the intended muscle force, and a perturbation of the low-dimensional trajectory can modify the motor behaviour. As existing models of spinal motor control do not offer an adequate explanation of rotation1,2, we propose a theory of neural generation of movements from which this and other unresolved issues, such as speed regulation, force control and multifunctionalism, are readily explained.

Emerging principles of spacetime in brains: Meeting report on spatial neurodynamics.

How do neurons and networks of neurons interact spatially? Here, we overview recent discoveries revealing how spatial dynamics of spiking and postsynaptic activity efficiently expose and explain fundamental brain and brainstem mechanisms behind detection, perception, learning, and behavior.

Why Firing Rate Distributions Are Important for Understanding Spinal Central Pattern Generators.

Networks in the spinal cord, which are responsible for the generation of rhythmic movements, commonly known as central pattern generators (CPGs), have remained elusive for decades. Although it is well-known that many spinal neurons are rhythmically active, little attention has been given to the distribution of firing rates across the population. Here, we argue that firing rate distributions can provide an important clue to the organization of the CPGs. The data that can be gleaned from the sparse literature indicate a firing rate distribution, which is skewed toward zero with a long tail, akin to a normal distribution on a log-scale, i.e., a "log-normal" distribution. Importantly, such a shape is difficult to unite with the widespread assumption of modules composed of recurrently connected excitatory neurons. Spinal modules with recurrent excitation has the propensity to quickly escalate their firing rate and reach the maximum, hence equalizing the spiking activity across the population. The population distribution of firing rates hence would consist of a narrow peak near the maximum. This is incompatible with experiments, that show wide distributions and a peak close to zero. A way to resolve this puzzle is to include recurrent inhibition internally in each CPG modules. Hence, we investigate the impact of recurrent inhibition in a model and find that the firing rate distributions are closer to the experimentally observed. We therefore propose that recurrent inhibition is a crucial element in motor circuits, and suggest that future models of motor circuits should include recurrent inhibition as a mandatory element.

A Late Critical Period for Frequency Modulated Sweeps in the Mouse Auditory System.

Neuronal circuits are shaped by experience during time windows of increased plasticity in postnatal development. In the auditory system, the critical period for the simplest sounds-pure frequency tones-is well defined. Critical periods for more complex sounds remain to be elucidated. We used in vivo electrophysiological recordings in the mouse auditory cortex to demonstrate that passive exposure to frequency modulated sweeps (FMS) from postnatal day 31 to 38 leads to long-term changes in the temporal representation of sweep directions. Immunohistochemical analysis revealed a decreased percentage of layer 4 parvalbumin-positive (PV+) cells during this critical period, paralleled with a transient increase in responses to FMS, but not to pure tones. Preventing the PV+ cell decrease with continuous white noise exposure delayed the critical period onset, suggesting a reduction in inhibition as a mechanism for this plasticity. Our findings shed new light on the dependence of plastic windows on stimulus complexity that persistently sculpt the functional organization of the auditory cortex.

White Noise Background Improves Tone Discrimination by Suppressing Cortical Tuning Curves.

The brain faces the difficult task of maintaining a stable representation of key features of the outside world in noisy sensory surroundings. How does the sensory representation change with noise, and how does the brain make sense of it? We investigated the effect of background white noise (WN) on tuning properties of neurons in mouse A1 and its impact on discrimination performance in a go/no-go task. We find that WN suppresses the activity of A1 neurons, which surprisingly increases the discriminability of tones spectrally close to each other. To confirm the involvement of A1, we optogenetically excited parvalbumin-positive (PV+) neurons in A1, which have similar effects as WN on both tuning properties and frequency discrimination. A population model suggests that the suppression of A1 tuning curves increases frequency selectivity and thereby improves discrimination. Our findings demonstrate that the cortical representation of pure tones adapts during noise to improve sensory acuity.

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.

Hybrid Scheme for Modeling Local Field Potentials from Point-Neuron Networks.

With rapidly advancing multi-electrode recording technology, the local field potential (LFP) has again become a popular measure of neuronal activity in both research and clinical applications. Proper understanding of the LFP requires detailed mathematical modeling incorporating the anatomical and electrophysiological features of neurons near the recording electrode, as well as synaptic inputs from the entire network. Here we propose a hybrid modeling scheme combining efficient point-neuron network models with biophysical principles underlying LFP generation by real neurons. The LFP predictions rely on populations of network-equivalent multicompartment neuron models with layer-specific synaptic connectivity, can be used with an arbitrary number of point-neuron network populations, and allows for a full separation of simulated network dynamics and LFPs. We apply the scheme to a full-scale cortical network model for a ∼1 mm2 patch of primary visual cortex, predict laminar LFPs for different network states, assess the relative LFP contribution from different laminar populations, and investigate effects of input correlations and neuron density on the LFP. The generic nature of the hybrid scheme and its public implementation in hybridLFPy form the basis for LFP predictions from other and larger point-neuron network models, as well as extensions of the current application with additional biological detail.

Spike-Based Bayesian-Hebbian Learning of Temporal Sequences.

Many cognitive and motor functions are enabled by the temporal representation and processing of stimuli, but it remains an open issue how neocortical microcircuits can reliably encode and replay such sequences of information. To better understand this, a modular attractor memory network is proposed in which meta-stable sequential attractor transitions are learned through changes to synaptic weights and intrinsic excitabilities via the spike-based Bayesian Confidence Propagation Neural Network (BCPNN) learning rule. We find that the formation of distributed memories, embodied by increased periods of firing in pools of excitatory neurons, together with asymmetrical associations between these distinct network states, can be acquired through plasticity. The model's feasibility is demonstrated using simulations of adaptive exponential integrate-and-fire model neurons (AdEx). We show that the learning and speed of sequence replay depends on a confluence of biophysically relevant parameters including stimulus duration, level of background noise, ratio of synaptic currents, and strengths of short-term depression and adaptation. Moreover, sequence elements are shown to flexibly participate multiple times in the sequence, suggesting that spiking attractor networks of this type can support an efficient combinatorial code. The model provides a principled approach towards understanding how multiple interacting plasticity mechanisms can coordinate hetero-associative learning in unison.

Computing the Local Field Potential (LFP) from Integrate-and-Fire Network Models.

Leaky integrate-and-fire (LIF) network models are commonly used to study how the spiking dynamics of neural networks changes with stimuli, tasks or dynamic network states. However, neurophysiological studies in vivo often rather measure the mass activity of neuronal microcircuits with the local field potential (LFP). Given that LFPs are generated by spatially separated currents across the neuronal membrane, they cannot be computed directly from quantities defined in models of point-like LIF neurons. Here, we explore the best approximation for predicting the LFP based on standard output from point-neuron LIF networks. To search for this best "LFP proxy", we compared LFP predictions from candidate proxies based on LIF network output (e.g, firing rates, membrane potentials, synaptic currents) with "ground-truth" LFP obtained when the LIF network synaptic input currents were injected into an analogous three-dimensional (3D) network model of multi-compartmental neurons with realistic morphology, spatial distributions of somata and synapses. We found that a specific fixed linear combination of the LIF synaptic currents provided an accurate LFP proxy, accounting for most of the variance of the LFP time course observed in the 3D network for all recording locations. This proxy performed well over a broad set of conditions, including substantial variations of the neuronal morphologies. Our results provide a simple formula for estimating the time course of the LFP from LIF network simulations in cases where a single pyramidal population dominates the LFP generation, and thereby facilitate quantitative comparison between computational models and experimental LFP recordings in vivo.

Power laws from linear neuronal cable theory: power spectral densities of the soma potential, soma membrane current and single-neuron contribution to the EEG.

Power laws, that is, power spectral densities (PSDs) exhibiting 1/f(α) behavior for large frequencies f, have been observed both in microscopic (neural membrane potentials and currents) and macroscopic (electroencephalography; EEG) recordings. While complex network behavior has been suggested to be at the root of this phenomenon, we here demonstrate a possible origin of such power laws in the biophysical properties of single neurons described by the standard cable equation. Taking advantage of the analytical tractability of the so called ball and stick neuron model, we derive general expressions for the PSD transfer functions for a set of measures of neuronal activity: the soma membrane current, the current-dipole moment (corresponding to the single-neuron EEG contribution), and the soma membrane potential. These PSD transfer functions relate the PSDs of the respective measurements to the PSDs of the noisy input currents. With homogeneously distributed input currents across the neuronal membrane we find that all PSD transfer functions express asymptotic high-frequency 1/f(α) power laws with power-law exponents analytically identified as α∞(I) = 1/2 for the soma membrane current, α∞(p) = 3/2 for the current-dipole moment, and α∞(V) = 2 for the soma membrane potential. Comparison with available data suggests that the apparent power laws observed in the high-frequency end of the PSD spectra may stem from uncorrelated current sources which are homogeneously distributed across the neural membranes and themselves exhibit pink (1/f) noise distributions. While the PSD noise spectra at low frequencies may be dominated by synaptic noise, our findings suggest that the high-frequency power laws may originate in noise from intrinsic ion channels. The significance of this finding goes beyond neuroscience as it demonstrates how 1/f(α) power laws with a wide range of values for the power-law exponent α may arise from a simple, linear partial differential equation.

Frequency dependence of signal power and spatial reach of the local field potential.

Despite its century-old use, the interpretation of local field potentials (LFPs), the low-frequency part of electrical signals recorded in the brain, is still debated. In cortex the LFP appears to mainly stem from transmembrane neuronal currents following synaptic input, and obvious questions regarding the 'locality' of the LFP are: What is the size of the signal-generating region, i.e., the spatial reach, around a recording contact? How far does the LFP signal extend outside a synaptically activated neuronal population? And how do the answers depend on the temporal frequency of the LFP signal? Experimental inquiries have given conflicting results, and we here pursue a modeling approach based on a well-established biophysical forward-modeling scheme incorporating detailed reconstructed neuronal morphologies in precise calculations of population LFPs including thousands of neurons. The two key factors determining the frequency dependence of LFP are the spatial decay of the single-neuron LFP contribution and the conversion of synaptic input correlations into correlations between single-neuron LFP contributions. Both factors are seen to give low-pass filtering of the LFP signal power. For uncorrelated input only the first factor is relevant, and here a modest reduction (<50%) in the spatial reach is observed for higher frequencies (>100 Hz) compared to the near-DC ([Formula: see text]) value of about [Formula: see text]. Much larger frequency-dependent effects are seen when populations of pyramidal neurons receive correlated and spatially asymmetric inputs: the low-frequency ([Formula: see text]) LFP power can here be an order of magnitude or more larger than at 60 Hz. Moreover, the low-frequency LFP components have larger spatial reach and extend further outside the active population than high-frequency components. Further, the spatial LFP profiles for such populations typically span the full vertical extent of the dendrites of neurons in the population. Our numerical findings are backed up by an intuitive simplified model for the generation of population LFP.

LFPy: a tool for biophysical simulation of extracellular potentials generated by detailed model neurons.

Electrical extracellular recordings, i.e., recordings of the electrical potentials in the extracellular medium between cells, have been a main work-horse in electrophysiology for almost a century. The high-frequency part of the signal (≳500 Hz), i.e., the multi-unit activity (MUA), contains information about the firing of action potentials in surrounding neurons, while the low-frequency part, the local field potential (LFP), contains information about how these neurons integrate synaptic inputs. As the recorded extracellular signals arise from multiple neural processes, their interpretation is typically ambiguous and difficult. Fortunately, a precise biophysical modeling scheme linking activity at the cellular level and the recorded signal has been established: the extracellular potential can be calculated as a weighted sum of all transmembrane currents in all cells located in the vicinity of the electrode. This computational scheme can considerably aid the modeling and analysis of MUA and LFP signals. Here, we describe LFPy, an open source Python package for numerical simulations of extracellular potentials. LFPy consists of a set of easy-to-use classes for defining cells, synapses and recording electrodes as Python objects, implementing this biophysical modeling scheme. It runs on top of the widely used NEURON simulation environment, which allows for flexible usage of both new and existing cell models. Further, calculation of extracellular potentials using the line-source-method is efficiently implemented. We describe the theoretical framework underlying the extracellular potential calculations and illustrate by examples how LFPy can be used both for simulating LFPs, i.e., synaptic contributions from single cells as well a populations of cells, and MUAs, i.e., extracellular signatures of action potentials.

Modeling the spatial reach of the LFP.

The local field potential (LFP) reflects activity of many neurons in the vicinity of the recording electrode and is therefore useful for studying local network dynamics. Much of the nature of the LFP is, however, still unknown. There are, for instance, contradicting reports on the spatial extent of the region generating the LFP. Here, we use a detailed biophysical modeling approach to investigate the size of the contributing region by simulating the LFP from a large number of neurons around the electrode. We find that the size of the generating region depends on the neuron morphology, the synapse distribution, and the correlation in synaptic activity. For uncorrelated activity, the LFP represents cells in a small region (within a radius of a few hundred micrometers). If the LFP contributions from different cells are correlated, the size of the generating region is determined by the spatial extent of the correlated activity.

The local field potential reflects surplus spike synchrony.

While oscillations of the local field potential (LFP) are commonly attributed to the synchronization of neuronal firing rate on the same time scale, their relationship to coincident spiking in the millisecond range is unknown. Here, we present experimental evidence to reconcile the notions of synchrony at the level of spiking and at the mesoscopic scale. We demonstrate that only in time intervals of significant spike synchrony that cannot be explained on the basis of firing rates, coincident spikes are better phase locked to the LFP than predicted by the locking of the individual spikes. This effect is enhanced in periods of large LFP amplitudes. A quantitative model explains the LFP dynamics by the orchestrated spiking activity in neuronal groups that contribute the observed surplus synchrony. From the correlation analysis, we infer that neurons participate in different constellations but contribute only a fraction of their spikes to temporally precise spike configurations. This finding provides direct evidence for the hypothesized relation that precise spike synchrony constitutes a major temporally and spatially organized component of the LFP.

Electroencephalographic characterization of seizure activity in the synapsin I/II double knockout mouse.

We present a detailed comparison of the behavioral and electrophysiological development of seizure activity in mice genetically depleted of synapsin I and synapsin II (SynDKO mice), based on combined video and surface EEG recordings. SynDKO mice develop handling-induced epileptic seizures at the age of 2months. The seizures show a very regular behavioral pattern, where activity is initially dominated by truncal muscle contractions followed by various myoclonic elements. Whereas seizure behavior goes through clearly defined transitions, cortical activity as reflected by EEG recordings shows a more gradual development with respect to the emergence of different EEG components and the frequency of these components. No EEG pattern was seen to define a particular seizure behavior. However, myoclonic activity was characterized by more regular patterns of combined sharp waves and spikes. Where countable, the number of myoclonic jerks was significantly correlated to the number of such EEG complexes. Furthermore, some EEG recordings revealed epileptic regular discharges without clear behavioral seizure correlates. Our findings suggest that seizure behavior in SynDKO mice is not solely determined by cortical activity but rather reflects interplay between cortical activity and activity in other brain regions.

An automated online positioning system and simulation environment for multi-electrodes in extracellular recordings.

Extracellular recordings are a key tool to record the activity of neurons in vivo. Especially in the case of experiments with behaving animals, however, the tedious procedure of electrode placement can take a considerable amount of expensive and restricted experimental time. Furthermore, due to tissue drifts and other sources of variability in the recording setup, the position of the electrodes with respect to the recorded neurons can change causing low recording quality. The contributions of this work are threefold. We introduce a quality measure for the recording position of the electrode which should be maximized during recordings and is especially suitable for the use of multi-electrodes. An automated positioning system based on this quality measure is proposed. The system is able to find favorable recording positions and adapts the electrode position smoothly to changes of the neuron positions. Finally, we evaluate the system using a new simulator for extracellular recordings based on realistically reconstructed 3D neurons.

Intrinsic dendritic filtering gives low-pass power spectra of local field potentials.

The local field potential (LFP) is among the most important experimental measures when probing neural population activity, but a proper understanding of the link between the underlying neural activity and the LFP signal is still missing. Here we investigate this link by mathematical modeling of contributions to the LFP from a single layer-5 pyramidal neuron and a single layer-4 stellate neuron receiving synaptic input. An intrinsic dendritic low-pass filtering effect of the LFP signal, previously demonstrated for extracellular signatures of action potentials, is seen to strongly affect the LFP power spectra, even for frequencies as low as 10 Hz for the example pyramidal neuron. Further, the LFP signal is found to depend sensitively on both the recording position and the position of the synaptic input: the LFP power spectra recorded close to the active synapse are typically found to be less low-pass filtered than spectra recorded further away. Some recording positions display striking band-pass characteristics of the LFP. The frequency dependence of the properties of the current dipole moment set up by the synaptic input current is found to qualitatively account for several salient features of the observed LFP. Two approximate schemes for calculating the LFP, the dipole approximation and the two-monopole approximation, are tested and found to be potentially useful for translating results from large-scale neural network models into predictions for results from electroencephalographic (EEG) or electrocorticographic (ECoG) recordings.