Fikret E. Kapucu

Burst analysis tool for developing neuronal networks exhibiting highly varying action potential dynamics.

In this paper we propose a firing statistics based neuronal network burst detection algorithm for neuronal networks exhibiting highly variable action potential dynamics. Electrical activity of neuronal networks is generally analyzed by the occurrences of spikes and bursts both in time and space. Commonly accepted analysis tools employ burst detection algorithms based on predefined criteria. However, maturing neuronal networks, such as those originating from human embryonic stem cells (hESCs), exhibit highly variable network structure and time-varying dynamics. To explore the developing burst/spike activities of such networks, we propose a burst detection algorithm which utilizes the firing statistics based on interspike interval (ISI) histograms. Moreover, the algorithm calculates ISI thresholds for burst spikes as well as for pre-burst spikes and burst tails by evaluating the cumulative moving average (CMA) and skewness of the ISI histogram. Because of the adaptive nature of the proposed algorithm, its analysis power is not limited by the type of neuronal cell network at hand. We demonstrate the functionality of our algorithm with two different types of microelectrode array (MEA) data recorded from spontaneously active hESC-derived neuronal cell networks. The same data was also analyzed by two commonly employed burst detection algorithms and the differences in burst detection results are illustrated. The results demonstrate that our method is both adaptive to the firing statistics of the network and yields successful burst detection from the data. In conclusion, the proposed method is a potential tool for analyzing of hESC-derived neuronal cell networks and thus can be utilized in studies aiming to understand the development and functioning of human neuronal networks and as an analysis tool for in vitro drug screening and neurotoxicity assays.

Spectral Entropy Based Neuronal Network Synchronization Analysis Based on Microelectrode Array Measurements

Synchrony and asynchrony are essential aspects of the functioning of interconnected neuronal cells and networks. New information on neuronal synchronization can be expected to aid in understanding these systems. Synchronization provides insight in the functional connectivity and the spatial distribution of the information processing in the networks. Synchronization is generally studied with time domain analysis of neuronal events, or using direct frequency spectrum analysis, e.g., in specific frequency bands. However, these methods have their pitfalls. Thus, we have previously proposed a method to analyze temporal changes in the complexity of the frequency of signals originating from different network regions. The method is based on the correlation of time varying spectral entropies (SEs). SE assesses the regularity, or complexity, of a time series by quantifying the uniformity of the frequency spectrum distribution. It has been previously employed, e.g., in electroencephalogram analysis. Here, we revisit our correlated spectral entropy method (CorSE), providing evidence of its justification, usability, and benefits. Here, CorSE is assessed with simulations and in vitro microelectrode array (MEA) data. CorSE is first demonstrated with a specifically tailored toy simulation to illustrate how it can identify synchronized populations. To provide a form of validation, the method was tested with simulated data from integrate-and-fire model based computational neuronal networks. To demonstrate the analysis of real data, CorSE was applied on in vitro MEA data measured from rat cortical cell cultures, and the results were compared with three known event based synchronization measures. Finally, we show the usability by tracking the development of networks in dissociated mouse cortical cell cultures. The results show that temporal correlations in frequency spectrum distributions reflect the network relations of neuronal populations. In the simulated data, CorSE unraveled the synchronizations. With the real in vitro MEA data, CorSE produced biologically plausible results. Since CorSE analyses continuous data, it is not affected by possibly poor spike or other event detection quality. We conclude that CorSE can reveal neuronal network synchronization based on in vitro MEA field potential measurements. CorSE is expected to be equally applicable also in the analysis of corresponding in vivo and ex vivo data analysis.

Analyzing the feasibility of time correlated spectral entropy for the assessment of neuronal synchrony

In this paper, we study neuronal network analysis based on microelectrode measurements. We search for potential relations between time correlated changes in spectral distributions and synchrony for neuronal network activity. Spectral distribution is quantified by spectral entropy as a measure of uniformity/complexity and this measure is calculated as a function of time for the recorded neuronal signals, i.e., time variant spectral entropy. Time variant correlations in the spectral distributions between different parts of a neuronal network, i.e., of concurrent measurements via different microelectrodes, are calculated to express the relation with a single scalar. We demonstrate these relations with in vivo rat hippocampal recordings, and observe the time courses of the correlations between different regions of hippocampus in three sequential recordings. Additionally, we evaluate the results with a commonly employed causality analysis method to assess the possible correlated findings. Results show that time correlated spectral entropy reveals different levels of interrelations in neuronal networks, which can be interpreted as different levels of neuronal network synchrony.

Joint analysis of extracellular spike waveforms and neuronal network bursts

BACKGROUND Neuronal networks are routinely assessed based on extracellular electrophysiological microelectrode array (MEA) measurements by spike sorting, and spike and burst statistics. We propose to jointly analyze sorted spikes and detected bursts, and hypothesize that the obtained spike type compositions of the bursts can provide new information on the functional networks. NEW METHOD Spikes are detected and sorted to obtain spike types and bursts are detected. In the proposed joint analysis, each burst spike is associated with a spike type, and the spike type compositions of the bursts are assessed. RESULTS The proposed method was tested with simulations and MEA measurements of in vitro human stem cell derived neuronal networks under different pharmacological treatments. The results show that the treatments altered the spike type compositions of the bursts. For example, 6-cyano-7-nitroquinoxaline-2,3-dione almost completely abolished two types of spikes which had composed the bursts in the baseline, while bursts of spikes of two other types appeared more frequently. This phenomenon was not observable by spike sorting or burst analysis alone, but was revealed by the proposed joint analysis. COMPARISON WITH EXISTING METHODS The existing methods do not provide the information obtainable with the proposed method: for the first time, the spike type compositions of bursts are analyzed. CONCLUSIONS We showed that the proposed method provides useful and novel information, including the possible changes in the spike type compositions of the bursts due to external factors. Our method can be employed on any data exhibiting sortable action potential waveforms and detectable bursts.

Quantification and automatized adaptive detection of in vivo and in vitro neuronal bursts based on signal complexity.

In this paper, we propose employing entropy values to quantify action potential bursts in electrophysiological measurements from the brain and neuronal cultures. Conventionally in the electrophysiological signal analysis, bursts are quantified by means of conventional measures such as their durations, and number of spikes in bursts. Here our main aim is to device metrics for burst quantification to provide for enhanced burst characterization. Entropy is a widely employed measure to quantify regularity/complexity of time series. Specifically, we investigate the applicability and differences of spectral entropy and sample entropy in the quantification of bursts in in vivo rat hippocampal measurements and in in vitro dissociated rat cortical cell culture measurement done with microelectrode arrays. For the task, an automatized and adaptive burst detection method is also utilized. Whereas the employed metrics are known from other applications, they are rarely employed in the assessment of burst in electrophysiological field potential measurements. Our results show that the proposed metrics are potential for the task at hand.

On the threshold based neuronal spike detection, and an objective criterion for setting the threshold

In this paper, we investigate the workings of threshold (TH) based spike detection for neuronal extracellular field potential spikes. Thresholding is the most used spike detection method. In general, it is employed by setting the TH as per convention and without considering either the undetected or spurious spikes. In this paper, we provide insight in to the workings of thresholding, and proposed a new objective way to set the TH based on spike count histogram analysis. We illustrate the method with 2D and 3D simulations and analysis of measured data.

Network-Wide Adaptive Burst Detection Depicts Neuronal Activity with Improved Accuracy

Neuronal networks are often characterized by their spiking and bursting statistics. Previously, we introduced an adaptive burst analysis method which enhances the analysis power for neuronal networks with highly varying firing dynamics. The adaptation is based on single channels analyzing each element of a network separately. Such kind of analysis was adequate for the assessment of local behavior, where the analysis focuses on the neuronal activity in the vicinity of a single electrode. However, the assessment of the whole network may be hampered, if parts of the network are analyzed using different rules. Here, we test how using multiple channels and measurement time points affect adaptive burst detection. The main emphasis is, if network-wide adaptive burst detection can provide new insights into the assessment of network activity. Therefore, we propose a modification to the previously introduced inter-spike interval (ISI) histogram based cumulative moving average (CMA) algorithm to analyze multiple spike trains simultaneously. The network size can be freely defined, e.g., to include all the electrodes in a microelectrode array (MEA) recording. Additionally, the method can be applied on a series of measurements on the same network to pool the data for statistical analysis. Firstly, we apply both the original CMA-algorithm and our proposed network-wide CMA-algorithm on artificial spike trains to investigate how the modification changes the burst detection. Thereafter, we use the algorithms on MEA data of spontaneously active chemically manipulated in vitro rat cortical networks. Moreover, we compare the synchrony of the detected bursts introducing a new burst synchrony measure. Finally, we demonstrate how the bursting statistics can be used to classify networks by applying k-means clustering to the bursting statistics. The results show that the proposed network wide adaptive burst detection provides a method to unify the burst definition in the whole network and thus improves the assessment and classification of the neuronal activity, e.g., the effects of different pharmaceuticals. The results indicate that the novel method is adaptive enough to be usable on networks with different dynamics, and it is especially feasible when comparing the behavior of differently spiking networks, for example in developing networks.

Phase Coupling in EEG Burst Suppression during Propofol Anesthesia

This work analyzes the quadratic phase coupling between EEG rhythms during burst suppression period of propofol anesthesia. The main goal is to specify the phase relation between different EEG rhythms and search for possible systematic behavior. In order to achieve this, the presented work focuses on the detection of quadratic phase coupling (QPC) based on bispectral analysis. The study indicates that frequencies of the dominant components of the bursts vary, so the detection of QPC is practically not as simple as detecting synchrony between couples of well defined rhythms.

On electrophysiological signal complexity during biological neuronal network development and maturation

Developing neuronal populations are assumed to increase their synaptic interactions and generate synchronized activity, such as bursting, during maturation. These effects may arise from increasing interactions of neuronal populations and increasing simultaneous intra-population activity in developing networks. In this paper, we investigated the neuronal network activity and its complexity by means of self-similarity during neuronal network development.

Evaluation of the effective and functional connectivity estimators for microelectrode array recordings during in vitro neuronal network maturation

During maturation, neurons and neuronal ensembles interact and build connections. Changes in the network structure have effects on the overall electrophysiological activity, and consequently on the observable connectivity. In this paper, we assessed effective and functional connectivities during neuronal network development by means of directed connectivity and synchronization, respectively.