Spike Synchronization Research Articles

Synchrony changes on different time-scales during in vitro neuronal network development

Synchrony changes on different time-scales during in vitro neuronal network development

Phenotyping of Gap-Junctional Coupling in the Mouse Retina.

In the mammalian retina, gap junctions, made of connexin proteins, are found in all neuronal cell types and are important for the transmission of rod photoreceptor signals, spike synchronization, noise reduction, and signal averaging. There are several methods available to assess gap junctional coupling in the retina: simultaneous electrical recordings from two adjacent cells, cut-loading, and intracellular injection of gap junction-permeable tracers. Here, we focus on the latter as it allows precise targeting of the cell of interest and is suitable to assess tracer coupling in a wide variety of retinal cell types, e.g., horizontal cells, amacrine cells, and ganglion cells. Tracer coupling experiments are usually performed in the intact retina and can provide information on the extent of coupling, the identity of synaptic partners, and (when combined with immunohistochemistry or pharmacology) the underlying connexin or the regulation of gap junctions.

Effects of topologies on signal propagation in feedforward networks

We systematically investigate the effects of topologies on signal propagation in feedforward networks (FFNs) based on the FitzHugh-Nagumo neuron model. FFNs with different topological structures are constructed with same number of both in-degrees and out-degrees in each layer and given the same input signal. The propagation of firing patterns and firing rates are found to be affected by the distribution of neuron connections in the FFNs. Synchronous firing patterns emerge in the later layers of FFNs with identical, uniform, and exponential degree distributions, but the number of synchronous spike trains in the output layers of the three topologies obviously differs from one another. The firing rates in the output layers of the three FFNs can be ordered from high to low according to their topological structures as exponential, uniform, and identical distributions, respectively. Interestingly, the sequence of spiking regularity in the output layers of the three FFNs is consistent with the firing rates, but their firing synchronization is in the opposite order. In summary, the node degree is an important factor that can dramatically influence the neuronal network activity.

Firing rate equations require a spike synchrony mechanism to correctly describe fast oscillations in inhibitory networks.

Recurrently coupled networks of inhibitory neurons robustly generate oscillations in the gamma band. Nonetheless, the corresponding Wilson-Cowan type firing rate equation for such an inhibitory population does not generate such oscillations without an explicit time delay. We show that this discrepancy is due to a voltage-dependent spike-synchronization mechanism inherent in networks of spiking neurons which is not captured by standard firing rate equations. Here we investigate an exact low-dimensional description for a network of heterogeneous canonical Class 1 inhibitory neurons which includes the sub-threshold dynamics crucial for generating synchronous states. In the limit of slow synaptic kinetics the spike-synchrony mechanism is suppressed and the standard Wilson-Cowan equations are formally recovered as long as external inputs are also slow. However, even in this limit synchronous spiking can be elicited by inputs which fluctuate on a time-scale of the membrane time-constant of the neurons. Our meanfield equations therefore represent an extension of the standard Wilson-Cowan equations in which spike synchrony is also correctly described.

Synchrony-induced modes of oscillation of a neural field model.

We investigate the modes of oscillation of heterogeneous ring networks of quadratic integrate-and-fire (QIF) neurons with nonlocal, space-dependent coupling. Perturbations of the equilibrium state with a particular wave number produce transient standing waves with a specific temporal frequency, analogously to those in a tense string. In the neuronal network, the equilibrium corresponds to a spatially homogeneous, asynchronous state. Perturbations of this state excite the network's oscillatory modes, which reflect the interplay of episodes of synchronous spiking with the excitatory-inhibitory spatial interactions. In the thermodynamic limit, an exact low-dimensional neural field model describing the macroscopic dynamics of the network is derived. This allows us to obtain formulas for the Turing eigenvalues of the spatially homogeneous state and hence to obtain its stability boundary. We find that the frequency of each Turing mode depends on the corresponding Fourier coefficient of the synaptic pattern of connectivity. The decay rate instead is identical for all oscillation modes as a consequence of the heterogeneity-induced desynchronization of the neurons. Finally, we numerically compute the spectrum of spatially inhomogeneous solutions branching from the Turing bifurcation, showing that similar oscillatory modes operate in neural bump states and are maintained away from onset.

Neuronal correlates of ketamine and walking induced gamma oscillations in the medial prefrontal cortex and mediodorsal thalamus.

Alterations in the function of the medial prefrontal cortex (mPFC) and its major thalamic source of innervation, the mediodorsal (MD) thalamus, have been hypothesized to contribute to the symptoms of schizophrenia. The NMDAR antagonist ketamine, used to model schizophrenia, elicits a brain state resembling early stage schizophrenia characterized by cognitive deficits and increases in cortical low gamma (40–70 Hz) power. Here we sought to determine how ketamine differentially affects spiking and gamma local field potential (LFP) activity in the rat mPFC and MD thalamus. Additionally, we investigated the ability of drugs targeting the dopamine D4 receptor (D4R) to modify the effects of ketamine on gamma activity as a measure of potential cognitive therapeutic efficacy. Rats were trained to walk on a treadmill to reduce confounds related to hyperactivity after ketamine administration (10 mg/kg s.c.) while recordings were obtained from electrodes chronically implanted in the mPFC and MD thalamus. Ketamine increased gamma LFP power in mPFC and MD thalamus in a similar frequency range, yet did not increase thalamocortical synchronization. Ketamine also increased firing rates and spike synchronization to gamma oscillations in the mPFC but decreased both measures in MD thalamus. Conversely, walking alone increased both firing rates and spike-gamma LFP correlations in both mPFC and MD thalamus. The D4R antagonist alone (L-745,870) had no effect on gamma LFP power during treadmill walking, although it reversed increases induced by the D4R agonist (A-412997) in both mPFC and MD thalamus. Neither drug altered ketamine-induced changes in gamma power or firing rates in the mPFC. However, in MD thalamus, the D4R agonist increased ketamine-induced gamma power and prevented ketamine’s inhibitory effect on firing rates. Results provide new evidence that ketamine differentially modulates spiking and gamma power in MD thalamus and mPFC, supporting a potential role for both areas in contributing to ketamine-induced schizophrenia-like symptoms.

Synchronicity and Rhythmicity of Purkinje Cell Firing during Generalized Spike-and-Wave Discharges in a Natural Mouse Model of Absence Epilepsy.

Absence epilepsy is characterized by the occurrence of generalized spike and wave discharges (GSWDs) in electrocorticographical (ECoG) recordings representing oscillatory activity in thalamocortical networks. The oscillatory nature of GSWDs has been shown to be reflected in the simple spike activity of cerebellar Purkinje cells and in the activity of their target neurons in the cerebellar nuclei, but it is unclear to what extent complex spike activity is implicated in generalized epilepsy. Purkinje cell complex spike firing is elicited by climbing fiber activation and reflects action potential firing in the inferior olive. Here, we investigated to what extent modulation of complex spike firing is reflected in the temporal patterns of seizures. Extracellular single-unit recordings in awake, head-restrained homozygous tottering mice, which suffer from a mutation in the voltage-gated CaV2.1 calcium channel, revealed that a substantial proportion of Purkinje cells (26%) showed increased complex spike activity and rhythmicity during GSWDs. Moreover, Purkinje cells, recorded either electrophysiologically or by using Ca2+-imaging, showed a significant increase in complex spike synchronicity for both adjacent and remote Purkinje cells during ictal events. These seizure-related changes in firing frequency, rhythmicity and synchronicity were most prominent in the lateral cerebellum, a region known to receive cerebral input via the inferior olive. These data indicate profound and widespread changes in olivary firing that are most likely induced by seizure-related activity changes in the thalamocortical network, thereby highlighting the possibility that olivary neurons can compensate for pathological brain-state changes by dampening oscillations.

Stochastic spike synchronization in a small-world neural network with spike-timing-dependent plasticity

We consider the Watts–Strogatz small-world network (SWN) consisting of subthreshold neurons which exhibit noise-induced spikings. This neuronal network has adaptive dynamic synaptic strengths governed by the spike-timing-dependent plasticity (STDP). In previous works without STDP, stochastic spike synchronization (SSS) between noise-induced spikings of subthreshold neurons was found to occur in a range of intermediate noise intensities. Here, we investigate the effect of additive STDP on the SSS by varying the noise intensity. Occurrence of a “Matthew” effect in synaptic plasticity is found due to a positive feedback process. As a result, good synchronization gets better via long-term potentiation of synaptic strengths, while bad synchronization gets worse via long-term depression. Emergences of long-term potentiation and long-term depression of synaptic strengths are intensively investigated via microscopic studies based on the pair-correlations between the pre- and the post-synaptic IISRs (instantaneous individual spike rates) as well as the distributions of time delays between the pre- and the post-synaptic spike times. Furthermore, the effects of multiplicative STDP (which depends on states) on the SSS are studied and discussed in comparison with the case of additive STDP (independent of states). These effects of STDP on the SSS in the SWN are also compared with those in the regular lattice and the random graph.

Synaptic Plasticity and Spike Synchronisation in Neuronal Networks

Brain plasticity, also known as neuroplasticity, is a fundamental mechanism of neuronal adaptation in response to changes in the environment or due to brain injury. In this review, we show our results about the effects of synaptic plasticity on neuronal networks composed by Hodgkin-Huxley neurons. We show that the final topology of the evolved network depends crucially on the ratio between the strengths of the inhibitory and excitatory synapses. Excitation of the same order of inhibition revels an evolved network that presents the rich-club phenomenon, well known to exist in the brain. For initial networks with considerably larger inhibitory strengths, we observe the emergence of a complex evolved topology, where neurons sparsely connected to other neurons, also a typical topology of the brain. The presence of noise enhances the strength of both types of synapses, but if the initial network has synapses of both natures with similar strengths. Finally, we show how the synchronous behaviour of the evolved network will reflect its evolved topology.

Synaptic convergence regulates synchronization-dependent spike transfer in feedforward neural networks

Correlated neural activities such as synchronizations can significantly alter the characteristics of spike transfer between neural layers. However, it is not clear how this synchronization-dependent spike transfer can be affected by the structure of convergent feedforward wiring. To address this question, we implemented computer simulations of model neural networks: a source and a target layer connected with different types of convergent wiring rules. In the Gaussian-Gaussian (GG) model, both the connection probability and the strength are given as Gaussian distribution as a function of spatial distance. In the Uniform-Constant (UC) and Uniform-Exponential (UE) models, the connection probability density is a uniform constant within a certain range, but the connection strength is set as a constant value or an exponentially decaying function, respectively. Then we examined how the spike transfer function is modulated under these conditions, while static or synchronized input patterns were introduced to simulate different levels of feedforward spike synchronization. We observed that the synchronization-dependent modulation of the transfer function appeared noticeably different for each convergence condition. The modulation of the spike transfer function was largest in the UC model, and smallest in the UE model. Our analysis showed that this difference was induced by the different spike weight distributions that was generated from convergent synapses in each model. Our results suggest that, the structure of the feedforward convergence is a crucial factor for correlation-dependent spike control, thus must be considered important to understand the mechanism of information transfer in the brain.

Does Layer 4 in the Barrel Cortex Function as a Balanced Circuit when Responding to Whisker Movements?

Neurons in one barrel in layer 4 (L4) in the mouse vibrissa somatosensory cortex are innervated mostly by neurons from the VPM nucleus and by other neurons within the same barrel. During quiet wakefulness or whisking in air, thalamic inputs vary slowly in time, and excitatory neurons rarely fire. A barrel in L4 contains a modest amount of neurons; the synaptic conductances are not very strong and connections are not sparse. Are the dynamical properties of the L4 circuit similar to those expected from fluctuation-dominated, balanced networks observed for large, strongly coupled and sparse cortical circuits? To resolve this question, we analyze a network of 150 inhibitory parvalbumin-expressing fast-spiking inhibitory interneurons innervated by the VPM thalamus with random connectivity, without or with 1600 low-firing excitatory neurons. Above threshold, the population-average firing rate of inhibitory cortical neurons increases linearly with the thalamic firing rate. The coefficient of variation CV is somewhat less than 1. Moderate levels of synchrony are induced by converging VPM inputs and by inhibitory interaction among neurons. The strengths of excitatory and inhibitory currents during whisking are about three times larger than threshold. We identify values of numbers of presynaptic neurons, synaptic delays between inhibitory neurons, and electrical coupling within the experimentally plausible ranges for which spike synchrony levels are low. Heterogeneity in in-degrees increases the width of the firing rate distribution to the experimentally observed value. We conclude that an L4 circuit in the low-synchrony regime exhibits qualitative dynamical properties similar to those of balanced networks.

The Electroencephalogram in Learning Disability: Special Emphasis on Interictal Epileptiform Discharges

This study dealt with the electroencephalograms (EEGs) of 187 children with learningdisability. Special attention was given to interictal epileptiform discharges. The patients who had interictal epileptiform discharge in their EEG was labeled as spike group, in contrast with those without such activity as the control group. In the whole group, 45.9% were completely normal and an additional 10.6% had positive spikes as the only finding. Epileptiform activity was seen in 32%, mainly focal (mostly centrotemporal or occipital) less often generalized, with bilaterally synchronous spike and waves complexes seen in 13 children. Slow wave abnormalities (mainly frontal or temporal), nearly always mild in degree, were seen mainly in the spike group. These different findings suggest that learning disability is a condition often with organic changes in the form of EEG abnormality, at times with epileptiform activity that could contribute to a deficit in learning.

Enhancement of synchronized activity between hippocampal CA1 neurons during initial storage of associative fear memory.

Learning and memory storage requires neuronal plasticity induced in the hippocampus and other related brain areas, and this process is thought to rely on synchronized activity in neural networks. We used paired whole-cell recording in vivo to examine the synchronized activity that was induced in hippocampal CA1 neurons by associative fear learning. We found that both membrane potential synchronization and spike synchronization of CA1 neurons could be transiently enhanced after task learning, as observed on day 1 but not day 5. On day 1 after learning, CA1 neurons showed a decrease in firing threshold and rise times of suprathreshold membrane potential changes as well as an increase in spontaneous firing rates, possibly contributing to the enhancement of spike synchronization. The transient enhancement of CA1 neuronal synchronization may play important roles in the induction of neuronal plasticity for initial storage and consolidation of associative memory. The hippocampus is critical for memory acquisition and consolidation. This function requires activity- and experience-induced neuronal plasticity. It is known that neuronal plasticity is largely dependent on synchronized activity. As has been well characterized, repetitive correlated activity of presynaptic and postsynaptic neurons can lead to long-term modifications at their synapses. Studies on network activity have also suggested that memory processing in the hippocampus may involve learning-induced changes of neuronal synchronization, as observed in vivo between hippocampal CA3 and CA1 networks as well as between the rhinal cortex and the hippocampus. However, further investigation of learning-induced synchronized activity in the hippocampus is needed for a full understanding of hippocampal memory processing. In this study, by performing paired whole-cell recording in vivo on CA1 pyramidal cells (PCs) in anaesthetized adult rats, we examined CA1 neuronal synchronization before and after associative fear learning. We first found in naive animals that there was a low level of membrane potential (MP) synchronization and spike synchronization of CA1 PCs. In conditioned animals, we found a significant enhancement of both MP synchronization and spike synchronization, as observed on day 1 after learning, and this enhancement was transient and not observed on day 5. Accompanying learning-induced synchronized activity was a decreased firing threshold and rise time of suprathreshold MP changes as well as an increased spontaneous firing rate, possibly contributing to the enhanced spike synchronization. The transiently enhanced CA1 neuronal synchronization may have important roles in generating neuronal plasticity for hippocampal storage and consolidation of associative memory traces.

Locking of correlated neural activity to ongoing oscillations.

Population-wide oscillations are ubiquitously observed in mesoscopic signals of cortical activity. In these network states a global oscillatory cycle modulates the propensity of neurons to fire. Synchronous activation of neurons has been hypothesized to be a separate channel of signal processing information in the brain. A salient question is therefore if and how oscillations interact with spike synchrony and in how far these channels can be considered separate. Experiments indeed showed that correlated spiking co-modulates with the static firing rate and is also tightly locked to the phase of beta-oscillations. While the dependence of correlations on the mean rate is well understood in feed-forward networks, it remains unclear why and by which mechanisms correlations tightly lock to an oscillatory cycle. We here demonstrate that such correlated activation of pairs of neurons is qualitatively explained by periodically-driven random networks. We identify the mechanisms by which covariances depend on a driving periodic stimulus. Mean-field theory combined with linear response theory yields closed-form expressions for the cyclostationary mean activities and pairwise zero-time-lag covariances of binary recurrent random networks. Two distinct mechanisms cause time-dependent covariances: the modulation of the susceptibility of single neurons (via the external input and network feedback) and the time-varying variances of single unit activities. For some parameters, the effectively inhibitory recurrent feedback leads to resonant covariances even if mean activities show non-resonant behavior. Our analytical results open the question of time-modulated synchronous activity to a quantitative analysis.

Cortical Spike Synchrony as a Measure of Input Familiarity.

Spike synchrony, which occurs in various cortical areas in response to specific perception, action, and memory tasks, has sparked a long-standing debate on the nature of temporal organization in cortex. One prominent view is that this type of synchrony facilitates the binding or grouping of separate stimulus components. We argue instead for a more general function: a measure of the prior probability of incoming stimuli, implemented by long-range, horizontal, intracortical connections. We show that networks of this kind-pulse-coupled excitatory spiking networks in a noisy environment-can provide a sufficient substrate for stimulus-dependent spike synchrony. This allows for a quick (few spikes) estimate of the match between inputs and the input history as encoded in the network structure. Given the ubiquity of small, strongly excitatory subnetworks in cortex, we thus propose that many experimental observations of spike synchrony can be viewed as signs of input patterns that resemble long-term experience-that is, of patterns with high prior probability.

Synchronization of Pulse-Coupled Phase Oscillators over Multi-FPGA Communication Links

This study proposes an implementation of pulse-coupled phase oscillators over multiple field-programmable-gate-array (FPGA) communication links. Two FPGAs are connected by a gigabit transceiver and a First-In First-Out interface. To verify the effect of communication delay between FPGAs on oscillator synchronization, we implement four oscillators on the multi-FPGA platform. We have successfully observed synchronization over two FPGAs correctly, despite of a 0.1 μs communication delay. The measurement results show that first spike synchronization requires 12.47 μs with a 3.2 Gbps communication throughput.

Detection and Evaluation of Spatio-Temporal Spike Patterns in Massively Parallel Spike Train Data with SPADE

Repeated, precise sequences of spikes are largely considered a signature of activation of cell assemblies. These repeated sequences are commonly known under the name of spatio-temporal patterns (STPs). STPs are hypothesized to play a role in the communication of information in the computational process operated by the cerebral cortex. A variety of statistical methods for the detection of STPs have been developed and applied to electrophysiological recordings, but such methods scale poorly with the current size of available parallel spike train recordings (more than 100 neurons). In this work, we introduce a novel method capable of overcoming the computational and statistical limits of existing analysis techniques in detecting repeating STPs within massively parallel spike trains (MPST). We employ advanced data mining techniques to efficiently extract repeating sequences of spikes from the data. Then, we introduce and compare two alternative approaches to distinguish statistically significant patterns from chance sequences. The first approach uses a measure known as conceptual stability, of which we investigate a computationally cheap approximation for applications to such large data sets. The second approach is based on the evaluation of pattern statistical significance. In particular, we provide an extension to STPs of a method we recently introduced for the evaluation of statistical significance of synchronous spike patterns. The performance of the two approaches is evaluated in terms of computational load and statistical power on a variety of artificial data sets that replicate specific features of experimental data. Both methods provide an effective and robust procedure for detection of STPs in MPST data. The method based on significance evaluation shows the best overall performance, although at a higher computational cost. We name the novel procedure the spatio-temporal Spike PAttern Detection and Evaluation (SPADE) analysis.

Detecting joint pausiness in parallel spike trains

BackgroundTransient periods with reduced neuronal discharge – called ‘pauses’ – have recently gained increasing attention. In dopamine neurons, pauses are considered important teaching signals, encoding negative reward prediction errors. Particularly simultaneous pauses are likely to have increased impact on information processing. Comparison with existing methodsAvailable methods for detecting joint pausing analyze temporal overlap of pauses across spike trains. Such techniques are threshold dependent and can fail to identify joint pauses that are easily detectable by eye, particularly in spike trains with different firing rates. New methodWe introduce a new statistic called pausiness that measures the degree of synchronous pausing in spike train pairs and avoids threshold-dependent identification of specific pauses. A new graphic termed the cross-pauseogram compares the joint pausiness of two spike trains with its time shifted analogue, such that a (pausiness) peak indicates joint pausing. When assessing significance of pausiness peaks, we use a stochastic model with synchronous spikes to disentangle joint pausiness arising from synchronous spikes from additional ‘joint excess pausiness’ (JEP). Parameter estimates are obtained from auto- and cross-correlograms, and statistical significance is assessed by comparison to simulated cross-pauseograms. ResultsOur new method was applied to dopamine neuron pairs recorded in the ventral tegmental area of awake behaving mice. Significant JEP was detected in about 20% of the pairs. ConclusionGiven the neurophysiological importance of pauses and the fact that neurons integrate multiple inputs, our findings suggest that the analysis of JEP can reveal interesting aspects in the activity of simultaneously recorded neurons.

Synchronization of bursting neurons with a slowly varying d. c. current

Bursting of neuronal firing is an interesting dynamical consequences depending on fast/slow dynamics. Certain cells in different brain regions produce spike-burst activity. We study such firing activity and its transitions to synchronization using identical as well as non-identical coupled bursting Morris-Lecar (M-L) neurons. Synchronization of different firing activity is a multi-time-scale phenomenon and burst synchronization presents the precursor to spike synchronization. Chemical synapses are one of the dynamical means of information processing between neurons. Electrical synapses play a major role for synchronous activity in a certain network of neurons. Synaptically coupled neural cells exhibit different types of synchronization such as in phase or anti-phase depending on the nature and strength of coupling functions and the synchronization regimes are analyzed by similarity functions. The sequential transitions to synchronization regime are examined by the maximum transverse Lyapunov exponents. Synchronization of voltage traces of two types of planar bursting mechanisms is explored for both kind of synapses under realistic conditions. The noisy influence effects on the transmission of signals and strongly acts to the firing activity (such as periodic firing and bursting) and integration of signals for a network. It has been examined using the mean interspike interval analysis. The transition to synchronization states of coupled and a network of bursting neurons may be useful for further research in information processing and even the origins of certain neurological disorders.

Hierarchical estimation of neural activity through explicit identification of temporally synchronous spikes

Extracellular recording from living neurons employing microelectrode arrays has attracted paramount attention in recent years as a way to investigate the functionality and disorders of the brain. To decipher useful information from the recorded signals, accurate and efficient neural spike activity detection and sorting becomes an essential prerequisite. Traditional approaches rely on thresholding to detect individual spikes and clustering to identify subset groups; however, these methods fail to identify temporally synchronous spikes due to neuronal synchrony. To address this challenge, we introduce a novel spike sorting algorithm incorporating both quantitative and probabilistic techniques to better approximate the ground truth information of the spike activity. A novel pre-clustering method for identifying key features that can form natural clusters and a dimension reduction technique for identifying the spiking activity are introduced. To address the temporal neuronal synchrony phenomenon leading to detection of multineural overlapped spikes, a procedure for template spike shape estimation and iterative recognition is developed employing the cross correlation methodology tailored to individual neuron's spike rate. A performance comparison between the proposed method and existing techniques in terms of the number of spikes identified and efficiency of sorting the spikes is presented. The outcome shows the effectiveness of the proposed method in identifying temporally synchronous spikes.