# Interaction between different cells in olfactory bulb and synchronous kinematic analysis.

1. IntroductionOlfactory system is an important part in sensory nervous system. Since there are many olfactory sensory neurons in our nose, we can sense the odour [1] and distinguish thousands of odours successfully. Each olfactory sensory neuron can express a specific odorant receptor. These receptors pass the odour information to olfactory bulb. Then the information will be the information will be transmitted to the entorhinal cortex of the brain by olfactory bulb, resulting in the sense of smell. As the first transfer station of the olfactory system, olfactory bulb plays an important role in processing olfactory information. When nerve signals coming from the depths of nasal get to the olfactory bulb, it will generate corresponding spatial and temporal coding [2]. The process of encoding information, in fact, is the process of further processing and handling information. Many scholars are discussing how the information processed by olfactory bulb is transferred synchronously to the entorhinal cortex.

In the olfactory bulb, there is an interaction between excitatory MCs and inhibitory GCs. The synaptic connections have several functions: (1) the lateral inhibition from GCs makes the receptive field of MCs smaller; (2) the interaction plays a role in oscillations formation in the olfactory bulb and also acts to synchronize the spikes of the MCs; (3) these synapses may be able to store the olfactory memory. Local field potential can be produced by the interaction between MCs and GCs [3]. The strength of synaptic connections determines the degree of the synchronization of the spikes, which is conducive to the integration of information. In this way, it can be more reliable to activate neurons of following olfactory cortex [4]. Thus there is a guarantee of transmitting signal and it lays a foundation for processing different olfactory information. Bazhenov et al. [5] and some other scholars put forward a network structural model of locust antenna, and then they verified the phenomenon of synchronous oscillation. On the concept of many-are-equal (MAE), Brody and Hopfield [6] further checked that olfactory system can identify odour by the synchronization of spikes of MCs. When the odorousness is in a certain range, the synchronization of group cells is unchanged. Some scholars like Horcholle-Bossavit et al. [7] studied the encoding mechanism of locust olfactory system through space-time coding model. The simulation results indicated that excitatory neurons were phase-locked in different oscillatory period and showed synchronization.

In conclusion, the synchronization of the spikes of MCs plays an important role in information encoding. The interaction between MCs and GCs in olfactory bulb is particularly critical. So it becomes necessary to analyse these two kinds of cells in detail. In this paper, we set up the dynamical model of MC and GC in the olfactory bulb. The simulation figures show the firing patterns of single MC and single GC, as well as these two kinds of cells having a coupling relationship. It shows that MCs can be synchronized by inhibition from GC, changing from irregular spiking to synchronous spiking. We wonder if the synchronization still exists when MCs are connected to different neuronal network of GCs and whether the degree of synchronization is different. To answer this question, we set up simple olfactory network models that consist of two MCs and nine GCs, discussing the influence of ring-like and grid-like neuronal network of GCs on the synchronization of two MCs.

2. Dynamics Model of Olfactory System

In this paper, MCs and GCs are implemented as conduction base models following Hodgkin-Huxley kinetics [8]. MCs are modeled using a single compartment while GCs are modeled using two compartments, representing a soma and dendrites. The membrane capacitance for all cells is 1 [mu]F/[cm.sup.2] and therefore not mentioned specifically in the voltage equations. In addition, the following units were used: conductance per unit area in ms/[cm.sup.2], current per unit area in [micro]A/[cm.sup.2], and voltage in mV and time (t) in ms. The membrane potentials are calculated with the following kinetic equation:

d[v.sub.k]/dt = -[I.sub.L] - [summation over i] [I.sub.i] - [summation over s][I.sub.s] + [summation over e][I.sub.e] + [summation over j][g.sub.k,j]([v.sub.j] - [v.sub.k]), (1)

where [I.sub.L], [I.sub.i], [I.sub.s], and [I.sub.e] are the leak, ionic, synaptic, and external currents, respectively. [g.sub.k,j], [v.sub.j], and [v.sub.k] represent the electronic conductance and voltages between adjacent compartments in a given cell.

There are two sodium currents [I.sub.Na] and [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] and three potassium currents [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], in MCs. The GC soma has a sodium current [I.sub.Na] and three potassium currents [I.sub.K], [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. The GC dendrites have one sodium current [I.sub.Na] and one potassium current [I.sub.K]. All the intrinsic currents are described by the following equation:

[I.sub.i] = [g.sub.i][m.sup.M][h.sup.H] (v - [E.sub.i]), (2)

where [g.sub.i] is the maximal conductance and [E.sub.i] the reversal potential. The activation and inactivation variables m and h raised to the power M and H, respectively; follow the kinetic equation as follows:

dm/dt = ([m.sub.[infinity]] - m)/[[tau].sub.m], dh/dt = ([h.sub.[infinity]] - h)/[[tau].sub.t], (3)

The parameters [m.sub.[infinity]] and [h.sub.[infinity]] and the time constants [tau] are not the same for different currents. Values for some parameters of MCs and GCs are given in Table 1 [8, 9].

For the current [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] in the MC and the currents [I.sub.Na] and [I.sub.K] in the GC dendrites and soma, we used values as shown in Figure 1 [9].

Synaptic connections between neurons are in the form of channel connectivity [10]. When MCs excite GCs, GCs will receive excitatory effect coming from receptors like APMA and NMDA. Accordingly, MCs will receive inhibitory effect from receptor GABA when inhibited by GCs. The strength of synaptic connections between cells is expressed by excitatory and inhibitory strength. The synaptic equations and relevant parameters are described in Table 2 [8].

3. Simulation Results and Analysis

3.1. Simulation of the Action Potential. Figures 2 and 3 are the spike trains of a single MC and GC under different external currents, respectively. Figure 2(a) is the membrane potential of MC without any external stimulus. When the stimulus gradually increases to 2 nA and 4 nA, the firing frequency increases at the same time. For GC, the spiking frequency also increases when stimulus varies from 0 nA to 10 nA and 20 nA. It shows that the response of cell changes with different stimulus, consistent with the sensitivity to outside stimuli of the olfactory system.

We can find that the firing frequency of MC has a great change when the stimulus only increases by 2 nA. While for GC, the spiking frequency just changes a little with stimulus adding 10 nA. By comparing these two group figures, maybe we can conclude that the firing patterns of MC change more obviously than GC with different stimulus.

3.2. Interaction between MC and GC. In the olfactory bulb, MC has an excitatory effect on GC, while GC has an inhibitory effect on MC. Cells coupled by synaptic connections and the firing patterns change according to different connection strength. We choose to change GABA conductance as an example to investigate the firing patterns of MC under different inhibition strength. In addition, we also discussed how the spike trains of MC varied when coupled with a single GC.

As is shown in Figures 4 and 5, the firing frequency of GC increases when it gets an excitatory effect from MC. And when MC gets an inhibitory effect from GC, the spiking frequency decreases. This phenomenon becomes more obvious as the inhibitory strength increases, as shown in Figures 6, 7, 8, and 9.

(1) MC Having Excitatory on GC. See Figures 4 and 5.

(2) GC Having an Inhibitory Effect on MC. See Figures 6,7, 8, and 9.

(3) MCs Spike Synchronously When Coupled with GC. The response of five MCs to stimulus without synaptic connections is shown in Figure 10(a). MC population spike asynchronously, for each MC getting different stimulus. When MCs got an inhibitory effect from GC, the firing rate of MCs decreases and they spike synchronously, as shown in Figure 11(a).

Figures 10(b) and 11(b) are the corresponding LFP plots of Figures 10(a) and 11(a). LFP reflects the excitement of local network, aiming to analyse the collaborative behavior of MC population. LFP focuses on the collaborative function between cells, hiding the individual role of each neuron. It is of great importance to discuss synchronous oscillation in the olfactory network. The LFP is defined as th mean value of the membrane potential of all MCs. When MC population spikes asynchronously, the LFP shows high frequency and low amplitude. While MC population fires synchronously, the LFP shows low frequency and high amplitude.

With the increase in the number of MCs, the curve of different membrane potential will be confused. Therefore, when the number reaches 20, we take the raster plot to characterize the spiking activity of cell population. Figure 12(a) is the raster plot of MCs firing irregularly with no synaptic connections. With an inhibition from GC, MCs spike regularly with a slow speed, as shown in Figure 13(a). Figures 12(b) and 13(b) are the corresponding LFP plot of Figures 12(a) and 13(a).

3.3. Network Model of Olfactory Bulb. Based on the previous work, we found that MCs spiked irregularly before; now they can spike synchronously with inhibition from a single GC. We wonder if the synchronization still exists when MCs connected to different neuronal network of GCs and whether the degree of synchronization is different.

To figure out this question, we set up simple olfactory network models consisting of two MCs, nine GCs, and several synaptic connections. We bring in GABA receptor in GC-MC and GC-GC synapses and also bring in AMPA and NMDA receptor in MC-GC synapses. We discussed how the ring-like and grid-like neuronal network of GCs influenced the synchronization of two MCs, respectively [11, 12].

In this model, the connection between MCs and two different topologies of GCs are shown in Figure 14. The ring-like neuronal network of GC only considers nearest-neighbour neurons while the grid-like neuronal network is with all nearest neighbour couplings. We have already known that coupled MCs could achieve synchronization with the variation of coupling strength without inhibition from GCs [13]. In order to exclude

the interference, we take the coupling strength of MCs as 0.01. In this way, there is weak coupling relationship between two MCs.

In the present paper, the synchronization degree of two spike trains is quantified by the ISI-distance [14], which is suitable for both two and multiple neurons. It uses the interspike interval (ISI) instead of spike frequency as the basic element of comparison. The method is a simple complementary approach that extracts information from the interspike intervals by evaluating the ratio of the instantaneous firing rates. It is parameter-free, time scale independent, and easy to visualize. Firstly, we take the instantaneous ISI-ratio between two interspike intervals, [x.sup.1.sub.ISI] and [x.sup.2.sub.ISI], and then normalize it. When two spike trains are the same, the quantity of ISI-ratio 2(t) becomes zero. While if the first spike train is much higher or lower than the second spike train, it approaches -1 and 1, respectively. In order to derive a measure of spike train distance, the absolute ISI-distance D is integrated over time. The more synchronous the coupled neurons are, the less the ISI-distance D is, and the complete synchronization state of the coupled neurons is achieved when D is equal to 0. The calculation formula is as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (4)

The membrane potential curve of two MCs without inhibition from GCs networks is shown in Figure 15. When coupled with ring-like and grid-like neuronal networks of GCs, the membrane potential plots of MCs have some change, as illustrated in Figures 16 and 17, respectively.

As shown in Figure 15, the solid line and dotted line are the spike trains of two MCs without connection with GCs network, respectively. These two MCs spike in different moment, and D is equal to 3.6575 x [10.sup.-3]. Compared with Figure 15, the spiking frequency of MCs decreases when coupled with ring-like network of GCs. At the same time, we can find that there is an apparent synchronization between two spike trains and D is equal to 0.7041 x [10.sup.-3]. When coupled with grid-like network, the response is similar to Figure 16. Two MCs spike synchronously in a low speed, and D is equal to 0.5811 x [10.sup.-3] at the time. Comparing Figure 15 with Figures 16 and 17, it shows that MCs spike irregularly without connection with GCs, turning to regular spiking slowly when coupled with GCs networks.

Comparing Figures 16 and 17, we can find similar synchronous spike trains of MCs when coupled with different topology network of GCs. MCs fire synchronously, but the synchronization degree is a little different. According to ISI-distance D, it shows that MCs become more synchronous when coupled with grid-like neuronal network of GCs.

4. Conclusion

In this paper, we have done some simulation work based on the dynamical model of MC and GC in the olfactory bulb. The simulation results show the firing patterns of a single MC and a single GC, as well as these two kinds of cells having a coupling relationship. In addition, we set up simple olfactory network models, discussing the influence of ring-like and grid-like neuronal networks of GCs on the synchronization of two MCs. ISI-distance method has been used to measure the synchronization degree of MCs. In summary, we have got four outcomes.

(1) The spike trains of MC and GC change with different external currents. When there is a larger external current, the spiking frequency of neuron increases accordingly. By comparison, we may conclude that the firing patterns of MC change more obviously than GC with different stimulus.

(2) MC has an excitatory effect on GC. While GC has an inhibitory effect on MC. The firing frequency of GC increases when it gets an excitatory effect from MC. And when MC gets an inhibitory effect from GC, the spiking frequency decreases. Cells coupled by synaptic connections and the firing patterns change accordingly with different connection strength.

(3) MCs spike asynchronously when there is no connection with GC, and LFP shows high frequency and low amplitude. When MCs got an inhibitory effect from GC, the firing rate of MCs decreases and they spike synchronously at the same time. The LFP shows low frequency and high amplitude.

(4) In the simple olfactory network models, we discussed how the ring-like and grid-like neuronal networks of GCs influenced the synchronization of two MCs, respectively. Through numerical analysis, it shows that MCs spike irregularly without connection with GCs, turning to regular spiking slowly when coupled with GCs networks. Compared with ring-like network, MCs become more synchronous when coupled with grid-like neuronal network of GCs.

Since the synchronization of the spikes of MCs plays an important role in information encoding and the interaction between MCs and GCs in olfactory bulb is particularly critical, we have made a detailed analysis on excited MC and inhibitory GC. These results are important both for understanding the sense of smell and for understanding the mechanisms of neural computation. The study of the neural basis of olfactory system lays a foundation for our later work about analysing how the neurons in olfactory bulb play a role in information coding comprehensively and quantificationally, starting from the overall structure of olfactory bulb. It can also help us study the odour pattern recognition based on olfactory bulb physiological structure. In this paper, our research about olfactory cognitive function is below the entorhinal cortex. In the future work, we will discuss the formation mechanism of smell more roundly, combining the information processing process of entorhinal cortex.

http://dx.doi.org/10.1155/2014/808792

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China under (Grant nos. 11232005, 11002055, and 11172086) and The Ministry of Education Doctoral Foundation (20120074110020).

References

[1] J. Tan and M. Luo, "Olfactory information processing by the olfactory bulb," Journal of Biological Physics, vol. 26, no. 3, pp. 194-208, 2010.

[2] S. Wolf, "Distributed processing and temporal codes in neuronal networks," Cognitive Neurodynamics, vol. 3, no. 3, pp. 189-196, 2009.

[3] Q. Zhen, "New development of olfactory recognition model," Journal of Wenzhou Medical College, vol. 39, no. 1, 2009.

[4] R. Wang, Z. Zhang, J. Qu, and J. Cao, "Phase synchronization motion and neural coding in dynamic transmission of neural information," IEEE Transactions on Neural Networks, vol. 22, no. 7, pp. 1097-1106, 2011.

[5] M. Bazhenov, M. Stopfer, M. Rabinovich et al., "Model of transient oscillatory synchronization in the locust antennal lobe," Neuron, vol. 30, no. 2, pp. 553-567, 2001.

[6] C. D. Brody and J. J. Hopfield, "Simple networks for spike-timing-based computation, with application to olfactory processing," Neuron, vol. 37, no. 5, pp. 843-852, 2003.

[7] G. Horcholle-Bossavit, B. Quenet, and O. Foucart, "Oscillation and coding in a formal neural network considered as a guide for plausible simulations of the insect olfactory system," BioSystems, vol. 89, no. 1-3, pp. 244-256, 2007

[8] J. N. Brea, L. M. Kay, and N. J. Kopell, "Biophysical model for gamma rhythms in the olfactory bulb via subthreshold oscillations," Proceedings of the National Academy of Sciences of the United States of America, vol. 106, no. 51, pp. 21954-21959, 2009.

[9] U. S. Bhalla and J. M. Bower, "Exploring parameter space in detailed single neuron models: simulations of the mitral and granule cells of the olfactory bulb," Journal of Neurophysiology, vol. 69, no. 6, pp. 1948-1965, 1993.

[10] D. Friedman and B. W. Strowbridge, "Both electrical and chemical synapses mediate fast network oscillations in the olfactory bulb," Journal of Neurophysiology, vol. 89, no. 5, pp. 2601-2610, 2003.

[11] J. Qu, R. Wang, Y. Du, and J. Cao, "Synchronization study in ring-like and grid-like neuronal networks," Cognitive Neurodynamics, vol. 6, no. 1, pp. 21-31, 2012.

[12] Y. Liu, R. Wang, Z. Zhang, and X. Jiao, "Analysis of stability of neural network with inhibitory neurons," Cognitive Neurodynamics, vol. 4, no. 1, pp. 61-68, 2010.

[13] Y. Du, R. Wang, F. Han, Q. Lu, and J. Qu, "Firing pattern and synchronization property analysis in a network model of the olfactory bulb," Cognitive Neurodynamics, vol. 6, no. 2, pp. 203-209, 2012.

[14] T. Kreuz, J. S. Haas, A. Morelli, H. D. I. Abarbanel, and A. Politi, "Measuring spike train synchrony," Journal of Neuroscience Methods, vol. 165, no. 1, pp. 151-161, 2007.

Ailei Xu (1,2), Ying Du (1,2), Rubin Wang (1,2), and Jinde Cao (3,4)

(1) School of Information Science and Engineering, East China University of Science and Technology, Shanghai 200237, China

(2) School of Science, East China University of Science and Technology, Shanghai 200237, China

(3) Department of Mathematics, Southeast University, Nanjing 210096, China

(4) Department of Mathematics, Faculty of Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia

Correspondence should be addressed to Rubin Wang; rbwang@163.com

Received 12 February 2014; Accepted 9 March 2014; Published 14 May 2014

Academic Editor: Guanghui Wen

TABLE 1: Some parameters in the dynamical model. Mitral cells Na ion M = 3, H = 1, [E.sub.Na] = 45, [a.sub.m] = 0.32 (v + channel 50) /{1 /exp [/(v + 50) /4]}, [b.sub.m] = 0.28(v + 23)/{exp [(v + 23)/5] /1}, [a.sub.h] = 0.128/exp [(v + 46)/18], = 4/{1 + exp [/(v + 23)/5]}, [[m.sub.[infinity]] = [a.sub.[m]/([a.sub.m] + [b.sub.m]), [h.sub.[infinity]] = [a.sub.h]/ ([a.sub.h] + [b.sub.h]), [[tau].sub.m] = 1/ ([a.sub.m] + [b.sub.m]), [[tau].sub.h] = 1 ([a.sub.h] + [b.sub.h]) [K.sub.s] [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] ion channel [K.sub.A] [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] ion channel [Na.sub.p] [I.sub.Nap] = g x [sub.Nap] (1/(exp(-(v + 51) ion channel /5) + 1))(v - 45) Granule cell soma [K.sub.A] [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] ion channel [K.sub.M] [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] ion channel TABLE 2: Equations and some parameters of synaptic currents. Synaptic currents GABA channel [i.sub.GABA] = [g.sub.GABA] (exp(-t/[[tau].sub.2]) - exp(-t/[[tau].sub.1]))(v-[e.sub.GABA]), [[tau].sub.1] = 1, [[tau].sub.2] = 200, [e.sub.GABA] = -80 APMA channel [i.sub.AMPA] = [g.sub.AMPA] exp(- t/[tau]) (v - [e.sub.AMPA]), [tau] = 3, [e.sub.AMPA] = 0 NMDA channel [i.sub.NMDA] = mg(Ron + Roff)(v - [e.sub.NMDA]), Ron' = -Ron/2.6, mg = 1/(1 + exp(-0.062v)/3.57), Roff' = -0.035 Roff, [e.sub.NMDA] = 0

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Xu, Ailei; Du, Ying; Wang, Rubin; Cao, Jinde |

Publication: | Discrete Dynamics in Nature and Society |

Article Type: | Report |

Date: | Jan 1, 2014 |

Words: | 3561 |

Previous Article: | The evolution of the strategies of innovation cooperation in scale-free network. |

Next Article: | The effects of the tractor and semitrailer routing problem on mitigation of carbon dioxide emissions. |

Topics: |