Citation: Shun Sakuma, Yuko Mizuno-Matsumoto, Yoshi Nishitani, Shinichi Tamura. Simulation of Spike Wave Propagation and Two-to-one Communication with Dynamic Time Warping[J]. AIMS Neuroscience, 2016, 3(4): 474-486. doi: 10.3934/Neuroscience.2016.4.474
[1] | Shinichi Tamura, Yoshi Nishitani, Chie Hosokawa . Feasibility of Multiplex Communication in a 2D Mesh Asynchronous Neural Network with Fluctuations. AIMS Neuroscience, 2016, 3(4): 385-397. doi: 10.3934/Neuroscience.2016.4.385 |
[2] | Yoshi Nishitani, Chie Hosokawa, Yuko Mizuno-Matsumoto, Tomomitsu Miyoshi, Shinichi Tamura . Classification of Spike Wave Propagations in a Cultured Neuronal Network: Investigating a Brain Communication Mechanism. AIMS Neuroscience, 2017, 4(1): 1-13. doi: 10.3934/Neuroscience.2017.1.1 |
[3] | Yoshi Nishitani, Chie Hosokawa, Yuko Mizuno-Matsumoto, Tomomitsu Miyoshi, Shinichi Tamura . Learning process for identifying different types of communication via repetitive stimulation: feasibility study in a cultured neuronal network. AIMS Neuroscience, 2019, 6(4): 240-249. doi: 10.3934/Neuroscience.2019.4.240 |
[4] | Yoshi Nishitani, Chie Hosokawa, Yuko Mizuno-Matsumoto, Tomomitsu Miyoshi, Shinichi Tamura . Effect of correlating adjacent neurons for identifying communications: Feasibility experiment in a cultured neuronal network. AIMS Neuroscience, 2018, 5(1): 18-31. doi: 10.3934/Neuroscience.2018.1.18 |
[5] | Eduardo Mercado III . Learning-Related Synaptic Reconfiguration in Hippocampal Networks: Memory Storage or Waveguide Tuning?. AIMS Neuroscience, 2015, 2(1): 28-34. doi: 10.3934/Neuroscience.2015.1.28 |
[6] | Eduardo Mercado III . Relating Cortical Wave Dynamics to Learning and Remembering. AIMS Neuroscience, 2014, 1(3): 185-209. doi: 10.3934/Neuroscience.2014.3.185 |
[7] | Shun Sakuma, Yuko Mizuno-Matsumoto, Yoshi Nishitani, Shinichi Tamura . Learning Times Required to Identify the Stimulated Position and Shortening of Propagation Path by Hebb’s Rule in Neural Network. AIMS Neuroscience, 2017, 4(4): 238-253. doi: 10.3934/Neuroscience.2017.4.238 |
[8] | Dimitris Pinotsis, Karl Friston . Gamma Oscillations and Neural Field DCMs Can Reveal Cortical Excitability and Microstructure. AIMS Neuroscience, 2014, 1(1): 18-38. doi: 10.3934/Neuroscience.2014.1.18 |
[9] | Jhunlyn Lorenzo, Stéphane Binczak, Sabir Jacquir . A multilayer-multiplexer network processing scheme based on the dendritic integration in a single neuron. AIMS Neuroscience, 2022, 9(1): 76-113. doi: 10.3934/Neuroscience.2022006 |
[10] | Damien Moore, Paul D Loprinzi . The association of self-reported physical activity on human sensory long-term potentiation. AIMS Neuroscience, 2021, 8(3): 435-447. doi: 10.3934/Neuroscience.2021023 |
The brain is a large network system that transmits information through spikes, which are short electrical signals in neurons. An action potential is transmitted to neurons with a time delay compared with the action potential in the preceding neuron. Over 0.5 msec passes before change occurs in the potential in the neuron receiving the input, and this is called synaptic propagation delay [1]. When a neuron fires, it is not able to fire for a certain length of time, and this is called a refractory period [2]. Synaptic propagation delays and refractory periods vary in different neurons. In addition, these times can fluctuate and are considered as a kind of noise. Information processing in brain activity might be conducted with spike propagation. In actuality information processing in the brain is rather stable despite the possible variabilities. Spike propagation and information transmission mechanisms are important research targets.
The first theory of information architecture was the cell-assembly theory proposed by Hebb in 1949 [3,4]. This theory was based on previous findings in associative memory and cell assembly. Okada et al. examined the relationship between associative memory and sparse coding [5]. The second theory involves the synfire chain model proposed by Abel [6]. This theory states that neuronal groups fire in a synchronous time pattern. However, the transfer of information between neurons and the function of the communication have not been fully elucidated. A question that remains is how neural activity propagates through cortical networks that are connected through synapses. Tanaka et al. showed that cortical networks use recurrent circuits in which a large number of neurons are bound to each other [7]. We developed a time-shift diagramming method that can be used to visualize the propagation of brain waves that communicate information [8,9,10].
The questions of how information communication is controlled, what constructs the information, and how the controlled information communication is constructed are unanswered. The answers to these questions are essential in investigations of the mechanisms underlying information communication. In order to resolve this question, we need to decode the sequence pattern of spike activity (analyses of the time-series patterns of firing) rather than examining the rate of spikes or action potential waveforms. In this study, we conducted a physiological experiment in cultured neurons from the rat hippocampus and recorded the spike trains. A significantly increased percentage of M-sequences spikes were recorded in the spike trains [11]. The raster plot analyses showed that linear-feedback shift-register circuits generated pseudo-random sequences including M-sequences [12]. However, we were unable to elucidate the meaning of the codes. In addition, we have shown that spikes propagate in neural networks because spike waves can be observed as code flow [13]. We then showed that cultured neuronal networks can be used in simulations with mesh-type two-dimensional neural network models composed of neurons that are modeled by integration and firing without leakage. In order to understand information processing in the brain, we simulated the firing activity of a large number of neurons in a neural network [14]. In our study, we found that stimulations at different sites were associated with the detection of different waves. In our previous study, we observed the generation of a number of spatiotemporal forms of spike wave propagations by various stimulated neurons in a cultural neuronal network [15]. We stimulated one of the two transmission neuron groups, and used the Dynamic Time Warping (DTW) method [16] to determine whether remote neurons receiving inputs can be used to identify the transmission neuronal group that was stimulated. We then confirmed the existence of two types of neurons: one that can identify the stimulated neuron and one that cannot. However, in that study, we were unable to elucidate the mechanisms underlying the results. With culture, parameters such as connection weights and synaptic delays are fixed and cannot be manipulated. In contrast, we can change the parameters in simulations. By changing the parameters and analyzing the different results obtained with different conditions, we can elucidate the mechanisms underlying the identification of the stimulated group of neurons.
In this study, we simulated spike responses to stimulations with various synaptic propagation delays and refractory periods. The simulation was conducted under the condition that the weights of the synapses were all fixed because we were only interested in examining the effects of fluctuations in the synaptic propagation delays and refractory periods without any confounding influence of variations in the synaptic weights. In order to focus on communication, we analyzed the information-flow of the network.
A 25 × 25 two-dimensional neural network was implemented (Figure 2). We used an integrate-and-fire model without any leakage as the neuronal model [17,18,19]. As well known, neurons have the property of all-or-none. Therefore, we implement an accepting period which is more stringent than the leak in the model. Accepting periods randomly accept input spikes or ignore. We will verify that neurons can communicate under such condition. The leak was ignored for simplicity. Each neuron had connection weights to and from eight neighboring neurons. We generated random weights of the synapses in the beginning of the experiment, and these weights were fixed through the experiment because we were only interested in examining the effects of fluctuations in the synaptic propagation delays and refractory periods. Therefore, any potential effects of synaptic weight variations were excluded. Three neurons were simultaneously stimulated at 0.1 msec, as shown in Figure 2, because we found in a number of preliminary simulations that the stimulation of three or more neurons stabilizes information propagation. We characterized the three neurons as being in the stimulated neuron group or transmitting neuron group. Spike waves were propagated from the stimulated neuron group to the other neurons. In our previous wet-lab experiments, we applied a time sampling rate of 0.1 msec, which was considered a bin. Thus, the time unit of a bin was 0.1 msec. The instantaneous variances in the synaptic propagation delays and refractory periods were both set to 0.167, 0.333, 0.500, 0.667, 1.000, or 2.000 [bin2] (bin = 0.1 msec) [17,18]. The stimulations (T, L, and D) were applied, as shown in Figure 2.
We calculated the spike-interval sequences of one neuron between 0-120 msec for 10 trials (Figure 3).
We obtained the spike-interval sequences for 625 neurons (total, 6, 250 sequences) with stimulation T (Figure 4) and stimulation L.
We calculated DTW in a combination of 10 trials (total, 45 sets) to calculate the spike-interval time differences between the trial sets of stimulation T (Figure 5). DTW was calculated among 10 trials for each neuron. In addition, we conducted DTW in a combination of 10 trials (total, 45 sets) to calculate the spike-interval time differences between the trial sets of stimulation L. DTW is an algorithm that is used to measure the differences between two signal sequences with different time scales or expansions. When the difference is large, the value approaches 1, while the value approaches 0 when the difference is small. We measured the DTW values in a combination of 10 trials with the same stimulated groups (T or L) and calculated the average of the 90 sets (two trials of 45 sets). We refer to these calculations as Local DTW (See Figure 5).
We calculated DTW for the combination of the 10 trials (total, 100 sets) in stimulations L and T in order to calculate the spike-interval time differences between the trial sets (Figure 6). We calculated the average of the 100 sets. We refer to these calculations as Inter DTW.
We used two-sided t-tests at the 5% significance level to compare the Local DTW and Inter DTW results. We determined the number of neurons in which the Inter DTW was significantly larger than the Local DTW. The neurons with significantly larger test statistics were considered able to identify stimulations T or L.
Figure 7 shows the results of a simulation of spike propagation in 10 msec. The transmitting neurons were stimulated at 0.0 msec. In all of the stimulations (T, L, and D), the spike waves spread in all directions. The propagation route was changed in each trial.
Figure 8 shows the results of the positions of the firing neurons with different variances. We set the variances of the synaptic delay and refractory period as 0.167 (yellow in Figure 8) and 2.0 (red in Figure 8), respectively. The blue in Figure 8 indicates the results with both 0.167 and 2.0 variances. Figure 8 shows the positions 0.0-6.0 msec after the stimulation. The propagation speed for a variance of 2.0 (red) was faster than the speed for a variance of 0.167 (yellow). Thus, when the synaptic delays and refractory periods were changed, the propagation speeds increased as the variances increased.
In order to ascertain whether the difference between the mean value of Local DTW (90 set of trial combination) and the mean value of Inter DTW (100 set of trial combination) is statistically significant or not, a two-tailed t-test was conducted with a significance level of 5% at each cell. Among them, we picked up and showed the value of Local DTW and Inter DTW (mean ± SD) at cell number 98, 132, 313, 504 and 590 in Table 1. As results, the mean value of Inter DTW was significantly greater than that of Local DTW (p = 0.03) at cell number 590, while no significant difference was observed (p = 0.79) at cell number 313.
Cell No. | Inter DTW | Local DTW | Significance |
98 | 0.029 ± 0.026 | 0.020 ± 0.020 | S |
132 | 0.037 ± 0.023 | 0.039 ± 0.023 | NS |
313 | 0.022 ± 0.018 | 0.023 ± 0.018 | NS |
504 | 0.029 ± 0.025 | 0.021 ± 0.016 | S |
590 | 0.020 ± 0.018 | 0.014 ± 0.014 | S |
As described in Section 2.5, the neurons with significantly larger test statistics were considered able to identify stimulations T or L (Identifiable neurons).
Figure 9 shows the percentages of neurons that could identify the stimulated neurons for neurons with significantly larger test statistics vs. neurons with variances in the synaptic delay and refractory period. The vertical axis shows the percentage of identifying neurons among the 625. The horizontal axis shows the variances of the synaptic delay and refractory period in bin2. When the variances of the synaptic delay and refractory period increased from 0.167 to 0.667, the percentage of identifiable neurons with different firing time profiles for stimulations T and L increased. However, when the variance was over 0.667, the percentage of neurons that could identify the stimulated neurons decreased. In the case of the variances of the synaptic delay and refractory being 0, the t-test can’t be performed because the variance is 0. Therefore, the result is not described. Since when variance is around 0.6, there was no special difference in the spike waves, we didn’t show them in figure 8.
We examined whether it is possible to discriminate between stimulus T and stimulus L individually for each neuron. We should have multiply compared at the stage of raw data. However, for the sake of simplicity, we only t-tested the difference between the mean values of Local DTW and Inter DTW for stimulus T and stimulus L, respectively. In order to integrate them, we evaluated the rate of Identifiable neurons as Figure 9.
The spike propagation results showed that the propagation route and velocity changed in response to alterations in the variances of the synaptic propagation delay and refractory period. This occurred because the synaptic delays and refractory periods with larger variances resulted in faster propagation of the spikes. As a result, the stimulus arrived earlier. The DTW results showed many signaling routes that were spatially different compared with the route for the transmitting neurons to the receiving neurons. They were transmitted in parallel (Figure 10).
The spikes that reach a neuron first result in firing of the neuron. Thus, the approximate flow of information is through this spatially representative route. When the variances in the synaptic propagation delays and refractory periods were small, the information passed through the thick route in each trial of Figure 10.
The routes to each transmission neuron from the receiving neuron often overlapped from the half-way point. In that case, the receiving neuron was difficult to use to identify the transmitting neurons (Figure 11). If there is large instantaneous variation in the synaptic propagation delays and refractory periods, the information is thought to pass through temporally various representative routes (Figure 11).
Thus, if the variations in the synaptic propagation delays and refractory periods increase, the spike wave will pass through the various representative routes. This then increases the probability of identification. When the variance is too large with a large bin2, the receiving neuron will receive spike profiles with too much disturbance, which makes identification of the spike wave itself difficult. These results suggested that changes in the propagation routes of the firing of neurons required a little variance in the synaptic propagation delays and refractory periods. Thus, variation may stabilize multiplex communication. In this paper, we simulated the case where the variances in synaptic propagation delays and refractory periods are the same. In that case, we explained a reason that the curve of figure 9 has a peak, which seems having high possibility. However, it is necessary to evaluate it by changing the variances in synaptic propagation delays and refractory periods, separately. We like to explain that in another paper.
Many previous studies have examined neural networks based on their firing patterns from a macro point of view, while we studied neural networks from the point of view of communication [3,4,5,6,7]. In order to understand information processing in the brain, we conducted simulations that assumed interactions of the firing activities of a large number of neural networks in the present study. We used an integrate-and-fire neuronal model without leakage and a 25 × 25 two-dimensional neural network. We showed how stimulation of the neurons was transmitted to the required neurons in the simulation as spike propagation occurred in 10 msec after the transmitting neurons were stimulated at 0.1 msec in the simulation. The two neuron groups were also stimulated with synaptic propagation delays and refractory periods with different variances. When the fluctuations in the synaptic delays and refractory periods were changed, the propagation speed increased as the variance increased. These results suggested that to change the propagation route of the firing of the neurons, some variance in the synaptic propagation delays and refractory periods is required. We examined the percentages of neurons with significantly larger test statistics vs. those with variances in the synaptic delays and refractory periods. The results suggested that variations in the synaptic delays and refractory periods improved the stability of multiplex communication.
In future studies, we need to examine whether multi-directional spike waves can be identified in more than two directions and the three-dimensional structure of the network. The present culture experiment examined not just one stimulated neurons and receiving neuron, but plural. Thus, several neurons were simultaneously recorded because an external electrode was used in the culture experiment. The simulation in the present study was used to examine neuron identification. It is possible to analyze the simulation with more detail. Thus, future studies need to examine the raw characteristics of each neuron.
This study was supported in part by the Grant-in-Aid for Scientific Research of Exploratory Research JP21656100, JP25630176, JP16K12524 and Scientific Research (A) JP22246054 of Japan Society for the Promotion of Science.
The authors declare that there is no conflict of interest regarding the publication of this paper.
[1] |
Eugene M. Izhikevich (2006) Polychronization: Computation with spikes. Neural Compu 18: 245-282. doi: 10.1162/089976606775093882
![]() |
[2] | Hiroyoshi Miyagawa, Masashi Inoue (2013) Biophysics of neurons. Maruzen Publishing. |
[3] | Donald Olding Hebb (1972) Textbook of Psychology, Philadelphia: Saunders, Pa, USA, 3rd edition. |
[4] | Anders Lansner (2009) Associative memory models: from the cell-assembly theory to biophysically detailed cortex simulations, Trends Neurosci 32: 178-186. |
[5] | Masato Okada (1996) Notions of Associative Memory and Sparse Coding, Neural Networks 9: 1429-1458. |
[6] |
Yuval Aviell, David Horn, Moshe Abeles (2004) Synfire waves in small balanced networks. Neurcomputing 58-60: 123-127. doi: 10.1016/j.neucom.2004.01.032
![]() |
[7] |
Takuma Tanaka, Takeshi Kaneko, Toshio Aoyagi (2009) Recurrent infomax generates cell assemblies, neuronal avalanches, and simple cell-like selectivity. Neural Computation 21: 1038-1067. doi: 10.1162/neco.2008.03-08-727
![]() |
[8] |
Yuko Mizuno-Matsumoto, Masatsugu Ishijima, Kazuhiro Shinosaki, et al. (2001) Transient Global Amnesia (TGA) in an MEG Study. Brain Topography 13: 269-274. doi: 10.1023/A:1011176612377
![]() |
[9] | Yuko Mizuno-Matsumoto, Toshiki Yoshimine, Yasuo Nii, et al. (2001) Landau-Kleffner Syndrome: Localization of Epileptogenic lesion Using Wavelet-Crosscorrelation Analysis, Epilepsy Behavior 2: 288-294. |
[10] |
Yuko Mizuno-Matsumoto, Satoshi Ukai, Ryosuke Ishii, et al. (2005) Wavelet-crosscorrelation analysis: Non-stationary analysis of neurophysiological signals. Brain Topogr 17: 237-252. doi: 10.1007/s10548-005-6032-2
![]() |
[11] | Shinichi Tamura, Shigenori Nakano, Kozo Okazaki (1985) Optical code-multiplex transmission by Gold-sequences, IEEE/OSA J. Lightwave Tech 1: 21-127. |
[12] | Yoshi Nishitani, Chie Hosokawa, Yuko Mizuno-Matsumoto, et al. (2012) Detection of M-sequences from spike sequence in neuronal networks. Computational Intelligence and Neuroscience 2012, Article ID 862579, 9 pages. doi:10.1155/2012/862579. |
[13] | Shinichi Tamura, Yoshi Nishitani, Chie Hosokawa, et al. (2016) Spike code flow in cultured neuronal networks. Computational Intelligence and Neuroscience 2016, Article ID 7267691, 11 pages. http://dx.doi.org/10.1155/2016/7267691. |
[14] |
Shinichi Tamura, Yoshi Nishitani, Yakuya Kamimura, et al. (2013) Multiplexed Spatiotemporal Communication Model in Artificial Neural Networks. Auto Control Intelligent System 1: 121-130. doi: 10.11648/j.acis.20130106.11. doi: 10.11648/j.acis.20130106.11
![]() |
[15] | Yoshi Nishitani, Chie Hosokawa, Yuko Mizuno-Matsumoto, et al. (2016) Variance of spatiotemporal spiking patterns by different stimulated neurons in cultured neuronal networks. Int J Academic Res Reflect 4: 11-19. |
[16] | Müller Meinard (2007) Information Retrieval for Music and Motion, Springer 69-84. |
[17] | Wulfram Gerstner, Werner M. Kistler (2002) Spiking Neuron Models. Single Neurons, Populations, Plasticity, Cambridge University Press. |
[18] | Shinichi Tamura, Yoshi Nishitani, Chie Hosokawa, et al. (2016) Simulation of code spectrum and code flow of cultured neuronal networks. Comput Intell Neurosci 2016, Article ID 7186092, 12 pages. |
[19] |
Shinichi Tamura Yoshi Nishitani, Chie Hosokawa (2016) Feasibility of multiplex communication in a 2D mesh asynchronous neural network with fluctuations. AIMS Neurosci 3: 385-397. DOI: 10.3934/Neuroscience.2016.4.385. doi: 10.3934/Neuroscience.2016.4.385
![]() |
1. | Shun Sakuma, Yuko Mizuno-Matsumoto, Yoshi Nishitani, Shinichi Tamura, Learning Times Required to Identify the Stimulated Position and Shortening of Propagation Path by Hebb’s Rule in Neural Network, 2017, 4, 2373-7972, 238, 10.3934/Neuroscience.2017.4.238 | |
2. | Jewel YunHsuan Lee, Michael Stiber, Dong Si, 2018, Machine Learning of Spatiotemporal Bursting Behavior in Developing Neural Networks, 978-1-5386-3646-6, 348, 10.1109/EMBC.2018.8512358 |
Cell No. | Inter DTW | Local DTW | Significance |
98 | 0.029 ± 0.026 | 0.020 ± 0.020 | S |
132 | 0.037 ± 0.023 | 0.039 ± 0.023 | NS |
313 | 0.022 ± 0.018 | 0.023 ± 0.018 | NS |
504 | 0.029 ± 0.025 | 0.021 ± 0.016 | S |
590 | 0.020 ± 0.018 | 0.014 ± 0.014 | S |