|
1.INTRODUCTIONAcupuncture and moxibustion is an important part of Chinese medicine, and more than two thousand years of clinical experience have fully demonstrated its effectiveness[1]. Modern clinical and experimental studies have shown that acupuncture is highly effective in postoperative analgesia and neurological disorders[2,3]. Acupuncture treatment has long lacked normative standards of diagnostic and treatment and technical indicators, and the learning of acupuncture manipulation was often based on the method of apprenticeship. The efficacy of acupuncture is closely related to acupuncture frequency. Different acupuncture frequencies evoke different neuronal responding activity and have different regulatory mechanisms for human body functions[4]. Therefore, this paper introduces statistical models to analyze the neuronal spiking activity evoked by acupuncture, explore the responding mechanism of acupuncture with different frequencies, and provide a scientific basis for the selection of acupuncture frequency in clinical treatment. Spike is the most fundamental unit of communication during neural response. With the rapid development of multi-electrode acquisition techniques, experimenters have been able to acquire neuroelectric signals from hundreds of channels simultaneously. Therefore, the premise of all research work is to accurately obtain the spiking trains of single neuron from the collected signals of electrodes. Most classification algorithms were based on the discharge waveform. Defining one or several time windows for the spiking waveforms of specific neurons, and then sliding and comparing them on the time axis, and the spiking waveform with the same shape is the spike of that neuron. However, in this classification algorithm, both superposition of spiking waveforms and background noise can produce variable spiking waveform[5,6]. At this time, the classification effect will be significantly worse, the variable spiking waveform will be eliminated, and a large amount of spiking information will be omitted. In order to solve this problem, this paper takes the results of classification algorithm based on waveform as a priori knowledge, and uses Bayesian statistical method to optimize the classification results. The optimized classification results can effectively identify the variable spiking waveforms, greatly reduce the missed detection rate of spike, and improve the accuracy of classification results[7]. Neuron model analysis is an important method to study neural information encoding. In order to better analyze various phenomena of neurons, scholars have established a variety of different neuron models based on the real biological neural system. For example, Hodgkin-Huxley (H-H) neuron model[8], Hindmarsh-Rose (HR) neuron model[9], Morris-Lecar (M-L) neuron model[10], FitzHugh-Nagumo (FHN) neuron model[11], Leaky Integrate-and-Fire (LIF) neuron model[12]. In addition, mathematical models with abstract concepts such as the cascade model can be used to describe a variety of dynamic processes within neurons or in neural networks. When selecting a coding model, it is important to be able to interpret the model’s parameters and predicted results not only from a statistical perspective, but also from a biophysical perspective. In the real process of neuronal spike, the spiking activity of neurons is affected by external stimuli, the neuron’s own spike-history and the spike-history of other neurons[13]. Therefore, we used a probabilistic statistical model to encode the spiking activity of neurons, and add the coupling of stimulus and spike-history to the probabilistic statistical model. Then the Maximum Likelihood Estimation (MLE) method was used to estimate the relevant parameters and explain the responding mechanism of acupuncture of different frequencies[14,15]. 2.METHOD2.1Spiking sorting by Bayesian statistical modelIn this section, the results of algorithm for traditional waveform classification are taken as prior knowledge. and Bayesian statistical model is used to optimize the classification results, identify the variable waveforms that are missed or falsely detected, and accurately obtain the spiking trains of single neurons. According to the results of the algorithms for traditional waveform classification, the prior probability function of the spiking trains dj is: where dj(t) is a binary variable. When neuron j evokes spikes at time t, dj(t) = 1, otherwise, dj(t) = 0. Assuming that the membrane voltage recorded by the electrode at the time t is a linear superposition of all the spiking waveforms in the neuron cluster by different forms at time t. The generative model of the membrane voltage at time t is established as follow[7]: where, the column vector represents the spikes recorded by all electrodes at time t, nc is the number of muti-electrodes, and represents the spikes recorded by i th electrode at time t. All spiking waveforms are discretized into nτ parts on average, and the matrix denotes the τ th segment waveform of neuron j recorded by all electrodes. is the background white noise with zero mean. According to (1), the generative model for spikes throughout the entire observation period is obtained: where, W * D represents the convolution operation in (2), and η is a multivariable Gaussian white noise. Therefore, V − W * D follows a multivariate Gaussian distribution. Consequently, under the conditions of spike times and spiking waveforms, the conditional probability function of the spikes generated is: where Λ is the covariance matrix of the multidimensional Gaussian variable η, which not only represents the covariance of noise over time but also the covariance of noise caused by the spatial positions of the electrodes in the collection process. Based on the prior function, the conditional probability function, and Bayesian theory, the posterior probability function of spikes evoked by the spiking trains and spiking waveforms is: The spiking trains and the spiking waveforms are estimated by making the posterior probability obtain its maximum value[16-18]. 2.2Information coding by probabilistic statistical modelThe spiking activity of neurons is affected by many factors at the same time, and the study of neural activity usually involves three different covariates. First, neural spiking activity is often associated with external biological and behavioral covariates, such as sensory stimuli and behaviors or specific motor outputs[13]. Second, the neuron’s current spiking activity is also related to its own spike-history[19]. Finally, the advancement of modern technology allows simultaneous record spiking activity of multi-neuron, and show that the spiking activity of a particular neuron is related to the spike-history of other neurons[15]. In this paper, a probabilistic statistical model is used to study the encoding and decoding process of acupuncture neurons. Firstly, the encoding model of acupuncture neuronal response was established using the framework of linear-nonlinear-Poisson process. In linear filtering, the external covariates, spike-history of self and the spike-history of other neurons are analyzed simultaneously. The three components undergo linear filtering and are then subjected to a nonlinear exponential filtering to generate the intensity function of the Poisson process[20,21]: where ki is the coupling parameter of external stimulation to neuron i, hii is the coupling parameter of own spike-history of the neuron i, and hij is the coupling parameter of other neurons. Let θi = (ki, hii, hij), is the parameter set of the encoding model. Assuming that the number of spiking events occurring in the time window dt is nt, then there is nt ~ Poiss[λdt]. Thus, the logarithmic probability function generated by neuronal spiking trains is as follows: where c is a constant that does not depend on the parameter θ. Then, the MLE was used to fit the spiking trains evoked by acupuncture to the model and estimate the model parameters. The logarithmic probability function log P(D | θ) is the logarithmic likelihood function of the model, and we can obtain the MLE of the model parameters by setting it to its maximum value[22,23]. 3.RESULT3.1Classification of acupuncture electrical signalsUsing the spiking activity of neurons in the dorsal horn of the spinal cord from an acupuncture experiment on the Zusanli (ST36) acupoint in rats as the raw data. We perform both traditional classification algorithms and Bayesian statistical model based spiking classification, and compare their classification effects. In the acupuncture experiment, healthy adult rats were selected as experimental subjects, the acupuncture frequencies were 30 times/min, 60 times/min, 120 times/min and 180 times/min respectively. Taking the stimulation frequency of 60 times/minute as an example, we selected the classification results of the two algorithms in 0.3 seconds to compare, as shown in Figure 1. Fig. 1 (a) shows the spiking trains of two neurons detected by the traditional classification algorithm. Fig. 1 (b) shows the spiking trains of two neurons detected in the Bayesian statistical model. Fig. 1 (c) shows the raw data collected by the electrodes. In the figure, the blue dots are the spikes extracted by both algorithms, the red dots are the spikes extracted by the Bayesian statistical model, and the green dots are the spikes omitted by the Bayesian statistical model. From the Fig. 1, it can be seen that the Bayesian statistical model has effectively solved the problem of superposition of multiple spiking waveforms. In particular, in the red dotted box, the spikes recorded by the electrode are concentrated. It shows that the responding activity of neurons is intense, and the Bayesian statistical model can detect all the spikes. In addition, when dealing with noise, such as the spiking waveforms corresponding to the spike in the green dotted line box is within the background noise range, such spikes are eliminated. It can be seen that the Bayesian statistical model reduces the missing rate of traditional classification algorithms and greatly improves the accuracy of classification results. Under the conditions of acupuncture at four frequencies, the number of spiking events of two neurons estimated by the two classification algorithms is shown in Table 1. It can be concluded that the number of spikes detected by these two neurons in the Bayesian statistical model is significantly increased, and some even twice the number of spiking events in the traditional classification algorithm. Table 1.The number of spiking events in two classification algorithms.
In addition, we found that after the stimulation was increased to 120 times/min, although the frequency of the stimulation continued to increase, the production of spikes did not increase anymore. This is due to the limited number and availability of ion channels on neuron cell membranes, as well as the saturation of neurotransmitter release. When the frequency of stimulation is too high, the release of ion channels and neurotransmitters reaches saturation, resulting in no more spikes being generated. This phenomenon helps to maintain the stability and accuracy of the nervous system, ensuring the reliability of nerve signals during transmission. In medical research and practice, the stimulation frequency of neurons can be adjusted according to this principle to achieve the purpose of treatment and research. In clinical treatment, the commonly used frequency of stimulation is 30 times/minute to 120 times/minute. Therefore, in the follow-up study, we mainly analyzed the neuronal spiking activity under the three stimulus frequencies of 30 times/minute, 60 times/minute and 120 times/minute. 3.2Encoding and decoding of acupuncture spiking activityWe verify the accuracy of the MLE through a set of simulation data. First, the intensity function was used to generate simulated spiking trains, and the model parameters were set to the real values in Table 2. Then, the MLE is used to fit the model and estimate the parameters of the model. The estimated results are shown in Table 2. By comparing the real value and the estimated value in Table 2, the relative error distribution of the estimated model parameters ranges from 0.62% to 4.47%. It is concluded that the MLE can accurately fit the model and estimate the model parameters. Table 2.Real and estimated values of model parameters.
In order to further study the difference of neuronal responding under the condition of acupuncture manipulation with different frequencies, a probabilistic statistical model was used to encode and decode the neuronal responding activity evoked by acupuncture of different frequencies, and the coupling parameters were estimated by MLE. The estimated results are shown in Table 3. We found that the coupling parameters of stimulus k were significantly lower than the coupling parameters of spike-history h under all three frequencies of acupuncture. With the increase of frequency, the more historical spikes of neurons, the smaller the coupling parameters of stimulus. It is concluded that the spike-history is the main factor to evoke the neuronal response in the experiment of neuronal spiking activity of acupuncture. Table 3.The coupling parameters of different acupuncture frequencies.
4.CONCLUSIONSIn this paper, statistical models were used to analyze the neuronal responding activity evoked by acupuncture with different frequencies, and to reveal the mechanism of acupuncture with different frequencies. Firstly, in order to obtain accurate neuron spiking trains, Bayesian statistical model is used to optimize the results of algorithm for traditional waveform classification. Then, the neuronal spiking events evoked by different frequencies of acupuncture were statistically analyzed, and it was found that the number of neuronal spiking events increased gradually with the increase of acupuncture frequency. However, after the acupuncture frequency reached 120/min, the increase of stimulation frequency did not produce more neuronal spiking events due to the frequency adaptation of neurons to the saturation mechanism. Finally, a probabilistic statistical model was used to encode and decode the spiking activity of acupuncture neurons. We found that the coupling parameters of stimulus were significantly lower than the coupling parameters of spike-history under different acupuncture frequencies, indicating that spike-history is the main factor that induces the responding activity of acupuncture neurons. According to the study of the relationship between acupuncture frequency and neuronal spiking activity, this may provide a scientific basis for future clinical treatment. ACKNOWLEDGEMENTSThis paper is supported by the Fund of Scientific Research Project of Tianjin Education Commission (Grant Nos. 2022ZD038, 2022KJ114, 2022KJ106) and Research Initiation Project of Tianjin University of Technology and Education (KYQD202358). REFERENCESLiu, B., Chen, B., Guo, Y., and Tian, L.,
“Acupuncture – a national heritage of China to the world: international clinical research advances from the past decade,”
Acupuncture and Herbal Medicine, 1
(2), 65
–73
(2021). https://doi.org/10.1097/HM9.0000000000000017 Google Scholar
Smith, C. A., Collins, C. T., Levett, K. M., Armour, M., Dahlen, H. G., Tan, A. L., and Mesgarpour, B.,
“Acupuncture or acupressure for pain management during labour,”
Cochrane Library,
(2020). Google Scholar
Fan, Z., Dou, B., Wang, J., Wu, Y., Du, S., Li, J., Yao, K., Li, Y., Wang, S., et al.,
“Effects and mechanisms of acupuncture analgesia mediated by afferent nerves in acupoint microenvironments,”
Frontiers in Neuroscience, 17
(2024). https://doi.org/10.3389/fnins.2023.1239839 Google Scholar
Li, N., Guo, Y., Gong, Y., Zhang, Y., Fan, W., Yao, K., Chen, Z., Dou, B., Lin, X., et al.,
“The Anti-Inflammatory Actions and Mechanisms of Acupuncture from Acupoint to Target Organs via Neuro-Immune Regulation,”
Journal of Inflammation Research Volume, 14 7191
–7224
(2021). https://doi.org/10.2147/JIR.S341581 Google Scholar
Yang, N., Ma, S., Yang, J., Li, T., and Liu, C.,
“Standardizing therapeutic parameters of acupuncture in vascular dementia rats,”
Brain and Behavior, 10
(10),
(2020). https://doi.org/10.1002/brb3.v10.10 Google Scholar
Wu, S.-C., and Swindlehurst, A. L.,
“Direct feature extraction from multi-electrode recordings for spike sorting,”
Digital Signal Processing, 75 222
–231
(2018). https://doi.org/10.1016/j.dsp.2018.01.016 Google Scholar
Pillow, J. W., Shlens, J., Chichilnisky, E. J., and Simoncelli, E. P.,
“A Model-Based spike sorting algorithm for removing correlation artifacts in Multi-Neuron recordings,”
PloS One, 8
(5), e62123
(2013). https://doi.org/10.1371/journal.pone.0062123 Google Scholar
Hodgkin, A. L., and Huxley, A. F.,
“A quantitative description of membrane current and its application to conduction and excitation in nerve,”
Journal of Physiology, 117
(4), 500
–544
(1952). https://doi.org/10.1113/jphysiol.1952.sp004764 Google Scholar
Hindmarsh, J. L., and Rose, R. M.,
“A model of the nerve impulse using two first-order differential equations,”
Nature, 296
(5853), 162
–164
(1982). https://doi.org/10.1038/296162a0 Google Scholar
Morris, C., and Lecar, H.,
“Voltage oscillations in the barnacle giant muscle fiber,”
Biophysical Journal, 35
(1), 193
–213
(1981). https://doi.org/10.1016/S0006-3495(81)84782-0 Google Scholar
FitzHugh, R.,
“Impulses and physiological states in theoretical models of nerve membrane,”
Biophysical Journal, 1
(6), 445
–466
(1961). https://doi.org/10.1016/S0006-3495(61)86902-6 Google Scholar
Richardson, M. J. E., Brunel, N., and Hakim, V.,
“From subthreshold to Firing-Rate resonance,”
Journal of Neurophysiology, 89
(5), 2538
–2554
(2003). https://doi.org/10.1152/jn.00955.2002 Google Scholar
Truccolo, W., Eden, U. T., Fellows, M. R., Donoghue, J. P., and Brown, E. N.,
“A point process framework for relating neural spiking activity to spiking history, neural ensemble, and extrinsic covariate effects,”
Journal of Neurophysiology, 93
(2), 1074
–1089
(2005). https://doi.org/10.1152/jn.00697.2004 Google Scholar
Pillow, J. W., Paninski, L., Uzzell, V. J., Simoncelli, E. P., and Chichilnisky, E. J.,
“Prediction and Decoding of Retinal Ganglion Cell Responses with a Probabilistic Spiking Model,”
the Journal of Neuroscience, 25
(47), 11003
–11013
(2005). https://doi.org/10.1523/JNEUROSCI.3305-05.2005 Google Scholar
Pillow, J. W., Shlens, J., Paninski, L., Sher, A., Litke, A. M., Chichilnisky, E. J., and Simoncelli, E. P.,
“Spatiotemporal correlations and visual signalling in a complete neuronal population,”
Nature, 454
(7207), 995
–999
(2008). https://doi.org/10.1038/nature07140 Google Scholar
Wood, F., and Black, M. J.,
“A nonparametric Bayesian alternative to spike sorting,”
Journal of Neuroscience Methods, 173
(1), 1
–12
(2008). https://doi.org/10.1016/j.jneumeth.2008.04.030 Google Scholar
Yger, P., Spampinato, G. L., Esposito, E., Lefebvre, B., Deny, S., Gardella, C., Stimberg, M., Jetter, F., Zeck, G.,
“A spike sorting toolbox for up to thousands of electrodes validated with ground truth recordings in vitro and in vivo,”
eLife, 7
(2018). https://doi.org/10.7554/eLife.34518 Google Scholar
Prentice, J.S., Homann, J., Simmons, K.D., Tkačik, G., Balasubramanian, V., and Nelson, P.C.,
“Fast, scalable, Bayesian spike identification for Multi-Electrode arrays,”
PloS One, 6
(7),
(2011). https://doi.org/10.1371/journal.pone.0019884 Google Scholar
Keat, J., Reinagel, P., Reid, R. C., and Meister, M.,
“Predicting every spike,”
Neuron, 30
(3), 803
–817
(2001). https://doi.org/10.1016/S0896-6273(01)00322-1 Google Scholar
Paninski, L., Shoham, S., Fellows, M. R., Hatsopoulos, N.G., and Donoghue, J.P.,
“Superlinear population encoding of dynamic hand trajectory in primary motor cortex,”
the Journal of Neuroscience, 24
(39), 8551
–8561
(2004). https://doi.org/10.1523/JNEUROSCI.0919-04.2004 Google Scholar
Pillow, J. W., and Latham, P. E.,
“Neural characterization in partially observed populations of spiking neurons,”
Neural Information Processing Systems, 20 1161
–1168
(2007). Google Scholar
Kulkarni, J. E., and Paninski, L.,
“Common-input models for multiple neural spike-train data,”
Network, 18
(4), 375
–407
(2007). https://doi.org/10.1080/09548980701625173 Google Scholar
Smith, A. C., and Brown, E. N.,
“Estimating a State-Space Model from Point Process Observations,”
Neural Computation, 15
(5), 965
–991
(2003). https://doi.org/10.1162/089976603765202622 Google Scholar
|