Separating burst from background spikes in multichannel neuronal recordings using return map analysis

Fluctuation Scaling of Neuronal Firing and Bursting in Spontaneously Active Brain Circuits

We employed high-density microelectrode arrays to investigate spontaneous firing patterns of neurons in brain circuits of the primary somatosensory cortex (S1) in mice. We recorded from over 150 neurons for 10min in each of eight different experiments, identified their location in S1, sorted their action potentials (spikes), and computed their power spectra and inter-spike interval (ISI) statistics. Of all persistently active neurons, 92% fired with a single dominant frequency - regularly firing neurons (RNs) - from 1 to 8Hz while 8% fired in burst with two dominant frequencies - bursting neurons (BNs) - corresponding to the inter-burst (2-6Hz) and intra-burst intervals (20-160Hz). RNs were predominantly located in layers 2/3 and 5/6 while BNs localized to layers 4 and 5. Across neurons, the standard deviation of ISI was a power law of its mean, a property known as fluctuation scaling, with a power law exponent of 1 for RNs and 1.25 for BNs. The power law implies that firing and bursting patterns are scale invariant: the firing pattern of a given RN or BN resembles that of another RN or BN, respectively, after a time contraction or dilation. An explanation for this scale invariance is discussed in the context of previous computational studies as well as its potential role in information processing.

Sensitivity to Stimulus Irregularity Is Inherent in Neural Networks

Behavior is controlled by complex neural networks in which neurons process thousands of inputs. However, even short spike trains evoked in a single cortical neuron were demonstrated to be sufficient to influence behavior in vivo. Specifically, irregular sequences of interspike intervals (ISIs) had a more reliable influence on behavior despite their resemblance to stochastic activity. Similarly, irregular tactile stimulation led to higher rates of behavioral responses. In this study, we identify the mechanisms enabling this sensitivity to stimulus irregularity (SSI) on the neuronal and network levels using simulated spiking neural networks. Matching in vivo experiments, we find that irregular stimulation elicits more detectable network events (bursts) than regular stimulation. Dissecting the stimuli, we identify short ISIs-occurring more frequently in irregular stimulations-as the main drivers of SSI rather than complex irregularity per se. In addition, we find that short-term plasticity modulates SSI. We subsequently eliminate the different mechanisms in turn to assess their role in generating SSI. Removing inhibitory interneurons, we find that SSI is retained, suggesting that SSI is not dependent on inhibition. Removing recurrency, we find that SSI is retained due to the ability of individual neurons to integrate activity over short timescales ("cell memory"). Removing single-neuron dynamics, we find that SSI is retained based on the short-term retention of activity within the recurrent network structure ("network memory"). Finally, using a further simplified probabilistic model, we find that local network structure is not required for SSI. Hence, SSI is identified as a general property that we hypothesize to be ubiquitous in neural networks with different structures and biophysical properties. Irregular sequences contain shorter ISIs, which are the main drivers underlying SSI. The experimentally observed SSI should thus generalize to other systems, suggesting a functional role for irregular activity in cortex.

Burst Detection Methods

'Bursting', defined as periods of high-frequency firing of a neuron separated by periods of quiescence, has been observed in various neuronal systems, both in vitro and in vivo. It has been associated with a range of neuronal processes, including efficient information transfer and the formation of functional networks during development, and has been shown to be sensitive to genetic and pharmacological manipulations. Accurate detection of periods of bursting activity is thus an important aspect of characterising both spontaneous and evoked neuronal network activity. A wide variety of computational methods have been developed to detect periods of bursting in spike trains recorded from neuronal networks. In this chapter, we review several of the most popular and successful of these methods.

Neural Coding With Bursts-Current State and Future Perspectives

Neuronal action potentials or spikes provide a long-range, noise-resistant means of communication between neurons. As point processes single spikes contain little information in themselves, i.e., outside the context of spikes from other neurons. Moreover, they may fail to cross a synapse. A burst, which consists of a short, high frequency train of spikes, will more reliably cross a synapse, increasing the likelihood of eliciting a postsynaptic spike, depending on the specific short-term plasticity at that synapse. Both the number and the temporal pattern of spikes in a burst provide a coding space that lies within the temporal integration realm of single neurons. Bursts have been observed in many species, including the non-mammalian, and in brain regions that range from subcortical to cortical. Despite their widespread presence and potential relevance, the uncertainties of how to classify bursts seems to have limited the research into the coding possibilities for bursts. The present series of research articles provides new insights into the relevance and interpretation of bursts across different neural circuits, and new methods for their analysis. Here, we provide a succinct introduction to the history of burst coding and an overview of recent work on this topic.

Anti-correlations in the degree distribution increase stimulus detection performance in noisy spiking neural networks

Neuronal circuits in the rodent barrel cortex are characterized by stable low firing rates. However, recent experiments show that short spike trains elicited by electrical stimulation in single neurons can induce behavioral responses. Hence, the underlying neural networks provide stability against internal fluctuations in the firing rate, while simultaneously making the circuits sensitive to small external perturbations. Here we studied whether stability and sensitivity are affected by the connectivity structure in recurrently connected spiking networks. We found that anti-correlation between the number of afferent (in-degree) and efferent (out-degree) synaptic connections of neurons increases stability against pathological bursting, relative to networks where the degrees were either positively correlated or uncorrelated. In the stable network state, stimulation of a few cells could lead to a detectable change in the firing rate. To quantify the ability of networks to detect the stimulation, we used a receiver operating characteristic (ROC) analysis. For a given level of background noise, networks with anti-correlated degrees displayed the lowest false positive rates, and consequently had the highest stimulus detection performance. We propose that anti-correlation in the degree distribution may be a computational strategy employed by sensory cortices to increase the detectability of external stimuli. We show that networks with anti-correlated degrees can in principle be formed by applying learning rules comprised of a combination of spike-timing dependent plasticity, homeostatic plasticity and pruning to networks with uncorrelated degrees. To test our prediction we suggest a novel experimental method to estimate correlations in the degree distribution.

Euchromatin histone methyltransferase 1 regulates cortical neuronal network development

Heterozygous mutations or deletions in the human Euchromatin histone methyltransferase 1 (EHMT1) gene cause Kleefstra syndrome, a neurodevelopmental disorder that is characterized by autistic-like features and severe intellectual disability (ID). Neurodevelopmental disorders including ID and autism may be related to deficits in activity-dependent wiring of brain circuits during development. Although Kleefstra syndrome has been associated with dendritic and synaptic defects in mice and Drosophila, little is known about the role of EHMT1 in the development of cortical neuronal networks. Here we used micro-electrode arrays and whole-cell patch-clamp recordings to investigate the impact of EHMT1 deficiency at the network and single cell level. We show that EHMT1 deficiency impaired neural network activity during the transition from uncorrelated background action potential firing to synchronized network bursting. Spontaneous bursting and excitatory synaptic currents were transiently reduced, whereas miniature excitatory postsynaptic currents were not affected. Finally, we show that loss of function of EHMT1 ultimately resulted in less regular network bursting patterns later in development. These data suggest that the developmental impairments observed in EHMT1-deficient networks may result in a temporal misalignment between activity-dependent developmental processes thereby contributing to the pathophysiology of Kleefstra syndrome.

MUSIC-Expected maximization gaussian mixture methodology for clustering and detection of task-related neuronal firing rates

Researchers often rely on simple methods to identify involvement of neurons in a particular motor task. The historical approach has been to inspect large groups of neurons and subjectively separate neurons into groups based on the expertise of the investigator. In cases where neuron populations are small it is reasonable to inspect these neuronal recordings and their firing rates carefully to avoid data omissions. In this paper, a new methodology is presented for automatic objective classification of neurons recorded in association with behavioral tasks into groups. By identifying characteristics of neurons in a particular group, the investigator can then identify functional classes of neurons based on their relationship to the task. The methodology is based on integration of a multiple signal classification (MUSIC) algorithm to extract relevant features from the firing rate and an expectation-maximization Gaussian mixture algorithm (EM-GMM) to cluster the extracted features. The methodology is capable of identifying and clustering similar firing rate profiles automatically based on specific signal features. An empirical wavelet transform (EWT) was used to validate the features found in the MUSIC pseudospectrum and the resulting signal features captured by the methodology. Additionally, this methodology was used to inspect behavioral elements of neurons to physiologically validate the model. This methodology was tested using a set of data collected from awake behaving non-human primates.

A comparison of computational methods for detecting bursts in neuronal spike trains and their application to human stem cell-derived neuronal networks

We provide an unbiased quantitative assessment of eight existing methods for identifying bursts in neuronal spike trains. We reveal limitations in a number of commonly used burst detection techniques and provide recommendations for the best practice for accurate identification of bursts using existing techniques. An analysis of the ontogeny of bursting activity in a novel data set of recordings from human induced pluripotent stem cell-derived neuronal networks, using the highest-performing burst detectors from our study, is also presented.

Accurate identification of bursting activity is an essential element in the characterization of neuronal network activity. Despite this, no one technique for identifying bursts in spike trains has been widely adopted. Instead, many methods have been developed for the analysis of bursting activity, often on an ad hoc basis. Here we provide an unbiased assessment of the effectiveness of eight of these methods at detecting bursts in a range of spike trains. We suggest a list of features that an ideal burst detection technique should possess and use synthetic data to assess each method in regard to these properties. We further employ each of the methods to reanalyze microelectrode array (MEA) recordings from mouse retinal ganglion cells and examine their coherence with bursts detected by a human observer. We show that several common burst detection techniques perform poorly at analyzing spike trains with a variety of properties. We identify four promising burst detection techniques, which are then applied to MEA recordings of networks of human induced pluripotent stem cell-derived neurons and used to describe the ontogeny of bursting activity in these networks over several months of development. We conclude that no current method can provide “perfect” burst detection results across a range of spike trains; however, two burst detection techniques, the MaxInterval and logISI methods, outperform compared with others. We provide recommendations for the robust analysis of bursting activity in experimental recordings using current techniques.

Assembling A Multi-Feature EEG Classifier for Left–Right Motor Imagery Data Using Wavelet-Based Fuzzy Approximate Entropy for Improved Accuracy

An EEG classifier is proposed for application in the analysis of motor imagery (MI) EEG data from a brain–computer interface (BCI) competition in this study. Applying subject-action-related brainwave data acquired from the sensorimotor cortices, the system primarily consists of artifact and background removal, feature extraction, feature selection and classification. In addition to background noise, the electrooculographic (EOG) artifacts are also automatically removed to further improve the analysis of EEG signals. Several potential features, including amplitude modulation, spectral power and asymmetry ratio, adaptive autoregressive model, and wavelet fuzzy approximate entropy (wfApEn) that can measure and quantify the complexity or irregularity of EEG signals, are then extracted for subsequent classification. Finally, the significant sub-features are selected from feature combination by quantum-behaved particle swarm optimization and then classified by support vector machine (SVM). Compared with feature extraction without wfApEn on MI data from two data sets for nine subjects, the results indicate that the proposed system including wfApEn obtains better performance in average classification accuracy of 88.2% and average number of commands per minute of 12.1, which is promising in the BCI work applications.

Quantitative Analysis of Rat Dorsal Root Ganglion Neurons Cultured on Microelectrode Arrays Based on Fluorescence Microscopy Image Processing

Microelectrode Arrays (MEA) are devices for long term electrophysiological recording of extracellular spontaneous or evocated activities on in vitro neuron culture. This work proposes and develops a framework for quantitative and morphological analysis of neuron cultures on MEAs, by processing their corresponding images, acquired by fluorescence microscopy. The neurons are segmented from the fluorescence channel images using a combination of segmentation by thresholding, watershed transform, and object classification. The positioning of microelectrodes is obtained from the transmitted light channel images using the circular Hough transform. The proposed method was applied to images of dissociated culture of rat dorsal root ganglion (DRG) neuronal cells. The morphological and topological quantitative analysis carried out produced information regarding the state of culture, such as population count, neuron-to-neuron and neuron-to-microelectrode distances, soma morphologies, neuron sizes, neuron and microelectrode spatial distributions. Most of the analysis of microscopy images taken from neuronal cultures on MEA only consider simple qualitative analysis. Also, the proposed framework aims to standardize the image processing and to compute quantitative useful measures for integrated image-signal studies and further computational simulations. As results show, the implemented microelectrode identification method is robust and so are the implemented neuron segmentation and classification one (with a correct segmentation rate up to 84%). The quantitative information retrieved by the method is highly relevant to assist the integrated signal-image study of recorded electrophysiological signals as well as the physical aspects of the neuron culture on MEA. Although the experiments deal with DRG cell images, cortical and hippocampal cell images could also be processed with small adjustments in the image processing parameter estimation.