Evidence-Based Complementary and Alternative Medicine

Volume 2013 (2013), Article ID 160357, 11 pages

http://dx.doi.org/10.1155/2013/160357

Research Article

Modulation of Brain Electroencephalography Oscillations by Electroacupuncture in a Rat Model of Postincisional Pain

Jing Wang,¹ Jing Wang,² Xuezhu Li,¹ Duan Li,³ Xiao-Li Li,² Ji-Sheng Han,1,4 and You Wan1,4

¹Neuroscience Research Institute, Peking University, Beijing 100191, China

²State Key Lab of Cognitive Neuroscience and Learning, Beijing Normal University, Beijing 100875, China

³Institute of Electrical Engineering, Yanshan University, Qinhuangdao 066004, China

⁴Key Laboratory for Neuroscience, Ministry of Education/Ministry of Health, Beijing 100191, China

Received 21 January 2013; Revised 1 March 2013; Accepted 21 March 2013

Academic Editor: Lijun Bai

Copyright © 2013 Jing Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The present study aimed to investigate how ongoing brain rhythmical oscillations changed during the postoperative pain and whether electroacupuncture (EA) regulated these brain oscillations when it relieved pain. We established a postincisional pain model of rats with plantar incision to mimic the clinical pathological pain state, tested the analgesic effects of EA, and recorded electroencephalography (EEG) activities before and after the EA application. By analysis of power spectrum and bicoherence of EEG, we found that in rats with postincisional pain, ongoing activities at the delta-frequency band decreased, while activities at theta-, alpha-, and beta-frequency bands increased. EA treatment on these postincisional pain rats decreased the power at high-frequency bands especially at the beta-frequency band and reversed the enhancement of the cross-frequency coupling strength between the beta band and low-frequency bands. After searching for the PubMed, our study is the first time to describe that brain oscillations are correlated with the processing of spontaneous pain information in postincisional pain model of rats, and EA could regulate these brain rhythmical frequency oscillations, including the power and cross-frequency couplings.

1. Introduction

Brain rhythmical oscillations in the low (delta, theta, and alpha) and high (beta and gamma) frequencies of electroencephalography (EEG) have been demonstrated to be linked to broad varieties of perceptual, sensorimotor, and cognitive operations [1]. Interaction of oscillations at different frequencies, for example, cross-frequency phase synchronization between alpha, beta, and gamma oscillations, could be observed during working memory, perception, and consciousness [2].

Pain, as a perception, is subserved by an extended network of brain areas [3], or different brain networks are involved in the perception of pain [4]. Previous studies on acute pain disclosed that painful stimulation altered the activities of different frequency oscillations [5, 6], including their power and phase couplings.

It is also noticed that most EEG neurophysiological studies on acute pain were based on phasic pain models with experimental pain induced by short-lasting, noninvasive painful stimuli (e.g., laser noxious heat stimulation). When we consider clinical situations, it is obvious that tonic pathological pain models could better mimic clinical pain than phasic pain models [7].

Postoperative incisional pain is common in clinic. In the rat model of incisional pain, persistent pain existed for several days after the hind paw incision with peripheral and central sensitization [8, 9]. On human subjects after the incision, an imaging study observed increased brain activities in the anterior cingulate cortex (ACC), the insular cortex, the thalamus, the frontal cortex, and the somatosensory cortex [10], and it would be useful to know how the brain EEG oscillation changes in tonic pathological pain states like postoperative pain (postincisional pain in the rat model).

Acupuncture has been widely used in clinical settings. Acupuncture or electroacupuncture (EA) has therapeutic effects in various painful conditions, and these effects could last for a long period of time even hours after acupuncture application being terminated [11–13]. Neuroimaging studies revealed that acupuncture or EA application elicited widespread changes in cerebrocerebellar brain regions [14, 15], largely overlapped with the neural networks for both pain transmission and perception. Acupuncture could directly affect EEG activities on healthy volunteers as well as on animals [16, 17]. Experimental and clinical evidence indicated that pain could affect cognitive processes [18], default-mode network dynamics [19] and even decreased the grey matter volume of brain regions [20]. To examine the neural consequences of acupuncture or EA treatment on pain, it would be useful to determine how acupuncture or EA treatment modulates cortical activities under tonic pathological pain conditions.

We established a rat model of plantar incision to mimic the clinic pain and observed the analgesic effects of EA on this model. On this basis, with EEG study, we further investigated changes of spontaneous brain oscillations in the incisional pain and the EA modulation on EEG oscillations.

2. Materials and Methods

All experimental procedures were in accordance with the guideline of the International Association for the Study of Pain [21] and were approved by the Animal Care and Use Committee of our university. The behavioral experimenters were kept blind.

2.1. Animals and Housing

Adult Sprague-Dawley male rats were provided by the Department of Experimental Animal Sciences of our university (30 rats for the behavior test, weighing 180–230 g; and 16 rats for EEG recordings, weighing 300–350 g at the beginning of the experiment). They were housed individually in cage with free access to food and water. The temperature was maintained at about 22°C under natural light/dark cycles. Rats were habituated to the environment and handled daily for one week before the experiment.

2.2. The Plantar Incisional Pain Model of Rats

The rat model of incisional pain employed a 1 cm longitudinal incision with muscle involvement [22]. Briefly, the animal was placed in a sealed glass container with 5% isoflurane mixed with air to induce anesthesia and delivered 1.5–2% isoflurane via a nose cone to maintain the anesthesia during the following surgical operation. The left hind paw of the rat was sterilized with 10% povidone-iodine, and a sterile no. 11 scalpel blade was used to make a 1 cm long incision through the skin and fascia of the plantar hind paw including the underlying muscles, beginning 0.5 cm from the heel. The wound was closed with two mattress sutures of 5–0 nylon, covering with the antibiotic ointment.

2.3. EA Application

Before the incisional surgery, the rat was loosely immobilized by the bandage on a metal mesh floor with the head, hind legs, and the tail protruding and habituated for 3-4 days, at least 20 min a day, in order to minimize the discomfort and the tension during the EA operation and application.

Two hours after the plantar incision, rats were divided into three groups: the restriction (incision) group, which received plantar incision in the left hind paw and was loosely immobilized by the bandage; the sham EA group, which received needling at “acupoints” without electrical stimulation after the incision; and the EA group, which received EA application bilaterally after the incision. EA was applied according to the routine procedure as in our previous reports [23–26]. Stainless-steel needles (0.4 mm in diameter, 4 mm in length) were inserted into the “acupoints” on each hind leg. Two commonly used acupoints, “Zusanli” (ST36, 4 mm lateral to the anterior tubercle of the tibia, which is marked by a notch) and “Sanyinjiao” (SP6, 3 mm proximal to the medial malleolus, at the posterior border of the tibia), were stimulated with square waves of 0.2 ms in pulse width and 2 Hz in frequency from a Han’s Acupoint Nerve Stimulator (HANS, LH series, manufactured in our university). The EA intensities were increased in a stepwise manner at 1-2-3 mA, with each intensity lasting for 10 min.

2.4. Assessment of Mechanical Allodynia

Mechanical allodynia was assessed by measuring the 50% paw withdrawal threshold (PWT) as described in our previous reports [26, 27]. The 50% PWT in response to a series of von Frey filaments (Semmes-Weinstein Monofilaments, North Coast Medial Inc., San Jose, CA, USA) was determined by the up-down method [28]. The rat was restricted on a metal mesh floor covered with an inverted clear plastic cage (18 × 8 × 8 cm) and allowed a 15–20 min period for habituation. Eight von Frey filaments with approximately equal logarithmic incremental (0.224) bending forces were chosen (0.41, 0.70, 1.20, 2.00, 3.63, 5.50, 8.50, and 15.10 g). Each trial started with a von Frey force of 2.00 g delivered perpendicularly to the plantar surface of the left hind paw adjacent to the wound near the medial heel. An abrupt withdrawal of the foot during the stimulation or immediately after the removal of the filament was recorded as a positive response. Once a positive or a negative response was evoked, the next weaker or stronger filament was applied, respectively. This procedure was terminated until 6 stimuli after the first change in response occurred. The 50% PWT was calculated using the following formula:

where is the value of the last von Frey filament used (in log unit), is a value measured from the pattern of positive/negative responses, and which is the average interval (in log unit) between the von Frey filaments [29]. The mechanical allodynia of the left hind paw was tested at 6 time points, that is, before the incision, after the incision but before the EA application, 10, 20, and 30 min during the EA application, and 30 min after the EA application. The data were analyzed with one-way ANOVA followed by Tukey post hoc test, and was chosen as statistically significant level.

2.5. Electrode Implantation for EEG Recording

Sixteen male rats were anesthetized with sodium pentobarbital (50 mg/kg, i.p.). After removing the scalp of the rat and exposing the skull, 14 stainless steel screws (tip diameter 1 mm, impedance 300–350 Ω) with sockets were implanted bilaterally as epidural electrodes into the skull to record the cortical EEGs. The locations of these electrodes were determined by the method from Shaw et al. [30]: anterior frontal (FL1, FR1), anterior (A) +4.5 mm, lateral (L) ±1.5 mm; centrofrontal (FL2, FR2, PR1, PL1), A ±1.5 mm, L ±4.5 mm; lateral frontal (PL2, PR2), A 0.0 mm, L ±4.5 mm; frontooccipital (LFL, RFR, LPL, RPR), A −4.5 mm, L ±1.5 mm for LFL, RFR, and A −3.0 mm, L ±4.5 mm for LPL, RPR. The reference and the ground electrodes were positioned 2 mm and 4 mm caudal to the lambda, respectively (Figure S1 of the Supplementary Material available at http://dx.doi.org/10.1155/2013/160357). The electrodes were fixed to the skull with dental cement and had no any connections with muscles. Penicillin (6 × 10⁴ U, i.m.) was administrated for 3 consecutive days to prevent possible infection.

2.6. EEG Signals Collection and Analysis

After one week recovery from the surgery, rats were habituated to the restriction for 3-4 days as the above mentioned, then the EEG rats were randomly divided into EA group or control group (restriction only) ( for each group). Rats were loosely immobilized on a metal mesh floor for the convenience of EA application during EEG recordings. EEG signals were collected in awake rats during three sessions, that is, before the plantar incision, 1.5–2 h after the plantar incision, and after the EA application. Each session lasted for 25–30 min.

The ASA-Lab EEG/ERP recording system (ANT Inc., The Netherlands) was used. Data were analyzed offline with Matlab (The Mathworks, Natick, MA, USA) EEGLAB software. The movement artifact or the baseline drift was removed from all channels, and the channels with impedance values above 25 k Ohms were also discarded. The EEG signals were digitized at a sampling rate of 256 Hz, rereferenced to an average of residual channels, and filtered through a 1–45 Hz band pass to avoid the interference of 50 Hz signals.

The wavelet power spectrum was used to obtain the power of the on going EEG activities [31]. The Morlet wavelet transform was employed with the wavelet central angle frequency of 6 (). Five spectral bands were examined: 1–4, 4–8, 8–13, 13–30, 30–45 Hz, corresponding to delta, theta, alpha, beta, and gamma bands, respectively [32, 33], with a step of 0.5 Hz.

The bispectral analysis, including the amplitude and the phase information, is used to quantify the degree of quadratic phase coupling (QPC) among different frequency components of a signal [34]. Bicoherence method is the normalized form of the bispectral analysis; it is independent of the amplitude of the signal, therefore; it can be used as an indicator of phase coupling in nonlinear signals. In this study, general harmonic wavelet bicoherence was employed to measure the comodulation of oscillations between two frequency bands [35]. Signals were divided into a series of 2-second epochs, with an overlap of 75%. For each epoch, bicoherence values were computed in all pairs of frequencies from 1 to 45 Hz, with a step of 1 Hz and a bandwidth of 2 Hz. The same epoch as the power analysis in the above was used to calculate the filtered wavelet bicoherence value (FIWBIC).

2.7. Statistical Analysis of EEG Data

For power spectral data, a paired t-test was used to analyze at which frequency ranges the change of the EEG power was significant. In the EA and the restriction groups, the change of power was expressed as the percentage relative to the power value during the session before the incision, which was defined as (After − Before)/Before × 100%, and paired t-tests were used in the interior-group. Unpaired t-tests were further used between the two groups.

For the wavelet bicoherence data, the statistical analysis was focused on their characteristics to determine if a significant difference existed between different recording sessions; so prior to the comparison, the total bicoherence value at the frequency bands   and   was extracted, which was defined as

where is the bicoherence value (FIWBIC). This value is a measure of the degree of QPC between frequency bands and can be used to measure the phase coupling strength between different waves [35]. Then, a Wilcoxon rank-sum test was conducted to determine significant difference, with value less than 0.05 as a statistically significant standard.

3. Results

3.1. EA Treatment Attenuates Mechanical Allodynia in the Incisional Pain Model of Rats

We established the plantar incisional pain model of rats and applied 2 Hz EA to “Zusanli” and “Sanyinjiao” acupoints in bilateral sides of the hind paw. A 50% of PWT to calibrated von Frey filaments was measured over time. Mechanical allodynia was attenuated significantly during EA application, and this antinociceptive effect could last for at least 30 min after EA application compared to that in the sham EA group or in the restriction group (Figure 1).

Figure 1: The effect of EA treatment on mechanical allodynia in the postincisional pain model of rats. Mechanical allodynia of the left hind paw was tested before incision, before EA treatment, 10, 20, and 30 min during EA, and 30 min after EA application. A 50% of paw withdrawal threshold (PWT) significantly increased after EA. , compared with the restriction group; compared with the sham-treatment group. Data were presented as means ± SEM, .

3.2. Changes of EEG Oscillations in the Incisional Pain Model of Rats

To examine the ongoing brain activities induced by the plantar incision, we calculated the power of different oscillations for each recording session (Figure 2). Compared with the session before the incision, the relative spontaneous EEG power of delta-frequency oscillation decreased, while the power in theta, alpha, and beta bands increased in rats after plantar incision. In addition, no significant change in gamma power was observed.

(a)

(b)

Figure 2: Changes of EEG power spectra in the postincisional pain model of rats. (a) The average of absolute EEG power before the incision (black line) and after the incision (red line). The -axis represented power value, and -axis represented frequency bands. (b) The relative power (normalized to the overall power) at 5 frequency bands. A significant increase of relative power at theta, alpha, and beta frequencies as well as a decrease at delta frequency was observed after the incision (red columns), compared with the power before the plantar surgery (black columns). , compared with before incision. All data were expressed as means ± SEM.

Furthermore, we compared the power of each channel at different frequency bands between different sessions. The corresponding topographies were shown in Figure 3(a). The locations where power changed significantly at different frequency bands during the sessions after the plantar incision and after EA application were summarized in Figure 3(b), respectively. We found that the changes mainly located at bilateral fronto-parietal lobes in low and medium frequency bands after the plantar incision.

(a)

(b)

Figure 3: Locations of electrodes with significant EEG power change in the postincisional pain model of rats. (a) Topographic mapping of EEG during two sessions before (upper) and after (lower) the plantar incision. Averaged EEG power densities ranged from delta to gamma bands. Values were color-coded and plotted at the corresponding position on the planar projection of the epidural surface and interpolated between electrodes (dots). (b) Locations of electrodes showing statistically significant difference () of power changes at different bands.

Synchronization of networks is often reflected by cross-frequency interaction, and the power spectrum could not reflect the phase information. In order to measure the degree of phase couplings at different frequency bands, we computed the filtered wavelet bicoherence (FIWBIC) values with general harmonic wavelet bicoherence [35]. As shown in Figure 4, the FIWBIC values were divided into bands from delta to gamma. A statistical analysis of estimated average of synchronization of FIWBIC values at different bifrequency bands was conducted on the 16 samples from sessions before the incision and compared to that from sessions after the incision. The Wilcoxon rank-sum test at was performed. The boxplots showed that the mean of the synchronization values among the beta, the alpha, and the theta increased, that is, the cross-frequency coupling strength between delta and alpha, and between theta and theta/alpha/beta, between alpha and alpha/beta during the session after the incision was significantly strengthened. Meanwhile, the coupling in the delta band before the incision was higher compared to that after the incision.

Figure 4: Phase couplings at different frequency bands: statistical analysis of the filtered wavelet bicoherence (FIWBIC) values in the postincisional pain model of rats. Boxplots stand for the phase coupling at local frequency bands, and the -axis represents the mean of synchronization. The synchronization values increased between delta and alpha, between theta and theta/alpha/beta, between alpha and alpha/beta, whereas the value decreased at delta band after the plantar incision. Wilcoxon rank-sum test was used. , , , .

3.3. EA Treatment-Induced EEG Power Changes in the Incisional Pain Model of Rats

In view of behavior results described above, both the sham EA treatment group and the restriction (without EA) group showed similar pain behaviors, and we chose the restriction group as control. The 16 rats were randomized into two groups: one group received EA application for 30 min from 2 h to 2.5 h after the plantar incision, and the other remained in restriction as control. We compared the power change at different oscillation bands in two groups (Figure 5). Figures 5(a) and 5(b) showed the change rate of averaged absolute power after EA treatment or restriction, separately. Figures 5(c) and 5(d) presented the change rate of relative power during the sessions before and after the EA application relative to the session before the incision in two groups, respectively. Only after EA treatment, the negative change rate relative to the baseline (before the incision) showed significant reduction in high-frequency bands, especially in the beta band. A tendency of reduction was also found in the gamma band. Additionally, the change of power values during the two sessions before EA or restriction was not different between two groups (Figure S2).

(a)

(b)

(c)

(d)

Figure 5: Changes of EEG power spectra after EA treatment in the postincisional pain model of rats. (a) The averaged absolute EEG power after EA treatment (red curve) compared with that before EA treatment (pink curve). The -axis represented power value, and -axis represented frequency bands. (b) The averaged absolute EEG power after restriction (dark blue color) compared with that before restriction (light blue curve). The -axis represents power value, and the -axis represents the frequency bands. (c) Change of the relative power in the EA group. The relative power value in EA group was normalized as the percentage relative to that before the incision. There was a significant decrease at the beta band after EA (red bar) compared with that after the incision (pink bar). (d) Change of the relative power in the restriction group. No significant changes were observed (light blue bar versus dark blue bar). (paired t-test).

Furthermore, after EA application, the change at the beta band mainly located over the electrodes of FR1, FR2, PR2, PL2, RFR, and LFL (Figure 6(b)), corresponding to the right fronto-parietal lobe, the left posterior parietal lobe, and the bilateral temporal lobes; meanwhile, the change at the gamma band located in the bilateral frontal lobes and the left posterior parietal lobe according to the location of electrodes (Figure 6(b)).

(a)

(b)

Figure 6: Locations of electrodes with significant EEG power change after EA treatment in incisional pain model of rats. (a) Topographic mapping of EEG power in beta band (left) and gamma band (right) after EA application compared with that before EA application. Averaged EEG power density values were color-coded and plotted at the corresponding position on the planar projection of the epidural surface and interpolated between electrodes (dots). (b) Electrodes with statistically significant difference () of power change in the beta and gamma band.

Similarly, we computed FIWBIC values in EA treatment group and the restriction group, respectively. The statistical results of Wilcoxon rank-sum tests at at the different frequency bands in two groups are shown in Figure 7 and Figure S3, respectively. From Figure 7, it is shown that after EA treatment, the strength of bifrequency coupling between beta and delta/theta/alpha attenuated significantly, whereas the coupling of beta-gamma bands became greater. In contrast, no significant difference was found in the strength of phase coupling between bands after restriction compared with that after the incision (Figure S3).

Figure 7: Phase couplings at different frequency bands: statistical analysis of the filtered wavelet bicoherence (FIWBIC) values after EA treatment in the postincisional pain model of rats. Boxplots stand for the phase coupling at local frequency bands, and the -axis represents the mean of synchronization. After EA application, the synchronization value between beta and delta/theta/alpha decreased, that is, synchronization became weaker; whereas the value between beta and gamma increased, that is, synchronization became stronger. Wilcoxon rank-sum test was used. , , , .

4. Discussion

In the incisional pain model of rats, by analysis of power spectrum and bicoherence of EEG, we found that the ongoing activities at the delta band decreased, while the activities at theta, alpha, and beta bands increased in the plantar-incisional pain rats; EA treatment decreased the power in high-frequency bands, especially at the beta band, and reversed the enhancement of the cross-frequency coupling strength between the beta and low-frequency bands.

4.1. Brain EEG Oscillations in Postincisional Pain

With rat EEG, we observed significant changes of the ongoing power spectra in different frequency oscillations ranging from the delta band to the beta band except the gamma band we also noticed the obvious changes of bi-frequency coherence during the postincisional pain. It is well known that pain perception is a multidimensional experience with sensory-discriminative, affective-emotional, and cognitive-evaluative components [3], and different brain areas are involved in different dimensions. Moreover, different cognitive systems are related to neuronal networks of different sizes and distribution, networks of different sizes oscillate at different frequencies, and mutual interactions of cross-frequency oscillations could be well positioned to regulate the multinetwork integration [36]. Thus, it might be conceivable that extensive range of neuronal oscillations participated in pathological pain cortical processing; any oscillation itself is not sufficient for integrating all the distributed information required for pain perception.

It has been accepted that delta-band oscillation is associated with compromised neuronal function [37]; theta-band is linked with emotional arousal [38]; beta-band may relate to the maintenance of the current sensorimotor or cognitive state, and its pathological enhancement may result in an abnormal persistence of the status quo and a deterioration of flexible behavior and cognitive control [39]. Thus, in our results, the reduced delta activity, as well as the increased theta and beta activities, might reflect the cortical overactivation induced by different pain dimensions in the resting state. Alpha oscillation is well known to mainly serve as a top-down controlled inhibitory mechanism [38]. One possible explanation for the increased alpha oscillation is that the evaluative component of pain, represented by paw lifting, was weakened by the activation of the descending pain inhibition network. Meanwhile, no significant change in the gamma oscillation was found during the postincisional pain. Owning that the gamma band is well recognized to be linked with feature binding, working memory, attention, or sensory selection [1], whereas our current study focused on the sustained, spontaneous pain in acute stage, it may be speculated that sensory selection or pain memory was not engaged in the pain processing. More importantly, the findings on alterations of phase coupling strength across distinct oscillations might reflect the synchronization of neuronal networks involved in different pain dimensions and the integration of pain inhibition and facilitation networks. It is also in accordance with the idea of neuronal processing with various simultaneous oscillations [36].

Topography findings showed the change of power in different frequency bands mainly located over bilateral frontal and parietal cortices (Figure 3). These are locations of the primary somatosensory cortex and the anterior cingulate cortex in rat. In this study, we selected the recording session from 1.5 to 2 h after the plantar incision. Given anatomical and physiological connections between right and left hemispheres, it was reasonable that the change of power occurred in bilateral cortices via communication and integration of persistent painful information in the whole brain.

4.2. EA Modulation on Incisional Pain-Related Brain Oscillations

In the present study, we investigated the modulation of EA treatment on the incisional pain-related brain oscillations in rats. As a basis, we firstly investigated the analgesic effects of EA treatment in the incisional pain rats. We found that EA treatment on “Zusanli” and “Sanyinjiao” acupoints on bilateral sides in hind paws relieved the mechanical allodynia in rats after the plantar incision. These results, in line with other reports [11–13], suggested that EA was effective in relieving plantar incisional pain and this antinociceptive effect could maintain at least 30 min.

We further explored the neural oscillation mechanism of EA analgesia. From our results, brain oscillations induced by EA application and those in nociceptive processing were not identical; the change of power after EA occurred at the beta band. It provides evidence that neuronal networks participated in EA treatment are different from those in pain perception. The networks oscillate at different frequencies, although brain areas activated by both EA treatment and incisional pain itself have large overlap from imaging studies [14, 15, 40], which is also in line with our results of topographic mapping.

As we mentioned before, the enhancement of the beta oscillation might be linked to the deterioration of flexible behavior and cognitive control [39]. Accumulating evidence suggests that EA facilitates the descending pain inhibitory pathway by increasing the release of opioid peptides in the central nervous system. Therefore, in our results, the decrease of activities at the beta band is reasonable, because EA treatment facilitates the descending inhibitory system of pain.

Prior studies indicated that inhibitory GABAergic interneurons network played a key role in the modulation of beta and gamma oscillations. Elevated endogenous GABA levels could cause the elevation of beta power [41, 42]. In the EA antinociceptive effects, GABA agonists showed reverse interaction with opioid receptor agonists [43, 44]. In addition, EA could decrease intracerebral GABA content in the cortex on the lesioned side in the rat model of Parkinson’s disease [45]. Taken together, it is conceivable that the decreased power of high-frequency oscillations in our results might reflect the inhibition of GABAergic interneurons induced by EA application in pain situation.

It has been accepted that low-frequency oscillations might be involved in the integration across widely spatially distributed neural assemblies and high-frequency oscillations (beta and gamma bands) distributed over a more limited topographic area. The integration of different local high-frequency oscillations is mediated by the large scale interactions of low-frequency oscillations. By analysis of the cross-frequency couplings, we found that EA treatment reversed the enhancement of the couplings between low-frequency bands and beta band in incisional pain and also strengthened the couplings between beta and gamma band. It may be speculated that EA exerted antinociceptive effect by modulation on power and cross-frequency coupling strength, which disturbed the cortex excitability and the multinetwork integration of nociceptive information in incisional pain.

In conclusion, the present study suggests that broad frequency oscillations ranging from delta-to beta-frequency bands are correlated with the cortical processing of the nociceptive information in the plantar incisional pain rats. EA reverses the increased beta power and the cross-frequency couplings between the beta and low-frequency bands induced by postincisional pain, suggesting that EA could regulate the neuronal networks involved in the central processing and the integration of spontaneous nociceptive information. These results can deepen our understanding in the central neuromodulatory mechanisms of EA analgesia.

Conflict of Interests

The authors have declared that no conflict of interests exists.

Acknowledgments

This research was supported by Grants from the National Basic Research Program of the Ministry of Science and Technology of China (2013CB531905), Key Project of Ministry of Education of China (109003), and the “111” Project of the Ministry of Education of China (B07001). Authors would like to thank Zifang Zhao for his help in optimization in Matlab computation, Dr. Miguel A. L. Nicolelis, Dr. Miguel Santos Pairs Vieira, Dr. Eric Thomson, and Gary Lehew in Duke University, USA, for their discussion and kind help in editing.

References

C. E. Schroeder and P. Lakatos, “Low-frequency neuronal oscillations as instruments of sensory selection,” Trends in Neurosciences, vol. 32, no. 1, pp. 9–18, 2009.

S. Palva and J. M. Palva, “New vistas for α-frequency band oscillations,” Trends in Neurosciences, vol. 30, no. 4, pp. 150–158, 2007.

X. Moisset and D. Bouhassira, “Brain imaging of neuropathic pain,” NeuroImage, vol. 37, supplement 1, pp. S80–S88, 2007.

A. V. Apkarian, M. C. Bushnell, R. D. Treede, and J. K. Zubieta, “Human brain mechanisms of pain perception and regulation in health and disease,” European Journal of Pain, vol. 9, no. 4, pp. 463–484, 2005.

C. Babiloni, F. Babiloni, F. Carducci et al., “Human brain oscillatory activity phase-locked to painful electrical stimulations: a multi-channel EEG study,” Human Brain Mapping, vol. 15, no. 2, pp. 112–123, 2002.

J. Wang, D. Li, X. Li et al., “Phase-amplitude coupling between theta and gamma oscillations during nociception in rat electroencephalography,” Neuroscience Letters, vol. 499, no. 2, pp. 84–87, 2011.

A. Mouraux and G. D. Iannetti, “Nociceptive laser-evoked brain potentials do not reflect nociceptive-specific neural activity,” Journal of Neurophysiology, vol. 101, no. 6, pp. 3258–3269, 2009.

M. M. Hämäläinen, G. F. Gebhart, and T. J. Brennan, “Acute effect of an incision on mechanosensitive afferents in the plantar rat hindpaw,” Journal of Neurophysiology, vol. 87, no. 2, pp. 712–720, 2002.

C. F. Villarreal, V. A. V. Kina, and W. A. Prado, “Participation of brainstem nuclei in the pronociceptive effect of lesion or neural block of the anterior pretectal nucleus in a rat model of incisional pain,” Neuropharmacology, vol. 47, no. 1, pp. 117–127, 2004.

E. M. Pogatzki-Zahn, C. Wagner, A. Meinhardt-Renner et al., “Coding of incisional pain in the brain: a functional magnetic resonance imaging study in human volunteers,” Anesthesiology, vol. 112, no. 2, pp. 406–417, 2010.

S. J. Wang, H. Y. Yang, and G. S. Xu, “Acupuncture alleviates colorectal hypersensitivity and correlates with the regulatory mechanism of TrpV1 and p-ERK,” Evidence-Based Complementary and Alternative Medicine, vol. 2012, Article ID 483123, 10 pages, 2012.

K. H. Chang, R. Won, I. Shim, H. Lee, and B. H. Lee, “Effects of electroacupuncture at BL60 on formalin-induced pain in rats,” Evidence-Based Complementary and Alternative Medicine, vol. 2012, Article ID 324039, 7 pages, 2012.

A. Li, Y. Zhang, L. Lao, et al., “Serotonin receptor 2A/C is involved in electroacupuncture inhibition of pain in an osteoarthritis rat model,” Evidence-based Complementary and Alternative Medicine, vol. 2011, Article ID 619650, 6 pages, 2011.

G. T. Lewith, P. J. White, and J. Pariente, “Investigating acupuncture using brain imaging techniques: the current state of play,” Evidence-Based Complementary and Alternative Medicine, vol. 2, no. 3, pp. 315–319, 2005.

J. Fang, X. Wang, H. Liu, et al., “The limbic-prefrontal network modulated by electroacupuncture at CV4 and CV12,” Evidence-Based Complementary and Alternative Medicine, vol. 2012, Article ID 515893, 11 pages, 2012.

S. F. Hsu, C. Y. Chen, M. D. Ke, C. H. Huang, Y. T. Sun, and J. G. Lin, “Variations of brain activities of acupuncture to TE5 of left hand in normal subjects,” American Journal of Chinese Medicine, vol. 39, no. 4, pp. 673–686, 2011.

G. Litscher, “Ten years evidence-based high-tech acupuncture part 3: a short review of animal experiments,” Evidence-Based Complementary and Alternative Medicine, vol. 7, no. 2, pp. 151–155, 2010.

D. A. Seminowicz and K. D. Davis, “Pain enhances functional connectivity of a brain network evoked by performance of a cognitive task,” Journal of Neurophysiology, vol. 97, no. 5, pp. 3651–3659, 2007.

M. N. Baliki, P. Y. Geha, A. V. Apkarian, and D. R. Chialvo, “Beyond feeling: chronic pain hurts the brain, disrupting the default-mode network dynamics,” Journal of Neuroscience, vol. 28, no. 6, pp. 1398–1403, 2008.

R. Ruscheweyh, M. Deppe, H. Lohmann et al., “Pain is associated with regional grey matter reduction in the general population,” Pain, vol. 152, no. 4, pp. 904–911, 2011.

M. Zimmermann, “Ethical guidelines for investigations of experimental pain in conscious animals,” Pain, vol. 16, no. 2, pp. 109–110, 1983.

T. J. Brennan, E. P. Vandermeulen, and G. F. Gebhart, “Characterization of a rat model of incisional pain,” Pain, vol. 64, no. 3, pp. 493–501, 1996.

Y. Wan, S. G. Wilson, J. S. Han, and J. S. Mogil, “The effect of genotype on sensitivity to electroacupuncture analgesia,” Pain, vol. 91, no. 1-2, pp. 5–13, 2001.

C. Huang, Z. Q. Huang, Z. P. Hu et al., “Electroacupuncture effects in a rat model of complete Freund's adjuvant-induced inflammatory pain: antinociceptive effects enhanced and tolerance development accelerated,” Neurochemical Research, vol. 33, no. 10, pp. 2107–2111, 2008.

C. Huang, H. Long, Y. S. Shi, J. S. Han, and Y. Wan, “Ketamine enhances the efficacy to and delays the development of tolerance to electroacupuncture-induced antinociception in rats,” Neuroscience Letters, vol. 375, no. 2, pp. 138–142, 2005.

G. G. Xing, F. Y. Liu, X. X. Qu, J. S. Han, and Y. Wan, “Long-term synaptic plasticity in the spinal dorsal horn and its modulation by electroacupuncture in rats with neuropathic pain,” Experimental Neurology, vol. 208, no. 2, pp. 323–332, 2007.

F. Y. Liu, Y. N. Sun, F. T. Wang, et al., “Activation of satellite glial cells in lumbar dorsal root ganglia contributes to neuropathic pain after spinal nerve ligation,” Brain Research, vol. 1427, pp. 65–77, 2012.

S. R. Chaplan, F. W. Bach, J. W. Pogrel, J. M. Chung, and T. L. Yaksh, “Quantitative assessment of tactile allodynia in the rat paw,” Journal of Neuroscience Methods, vol. 53, no. 1, pp. 55–63, 1994.

W. J. Dixon, “Efficient analysis of experimental observations,” Annual Review of Pharmacology and Toxicology, vol. 20, pp. 441–462, 1980.

F. Z. Shaw, R. F. Chen, H. W. Tsao, and C. T. Yen, “A multichannel system for recording and analysis of cortical field potentials in freely moving rats,” Journal of Neuroscience Methods, vol. 88, no. 1, pp. 33–43, 1999.

X. Li, X. Yao, J. Fox, and J. G. Jefferys, “Interaction dynamics of neuronal oscillations analysed using wavelet transforms,” Journal of Neuroscience Methods, vol. 160, no. 1, pp. 178–185, 2007.

D. Mantini, M. G. Perrucci, C. Del Gratta, G. L. Romani, and M. Corbetta, “Electrophysiological signatures of resting state networks in the human brain,” Proceedings of the National Academy of Sciences of the United States of America, vol. 104, no. 32, pp. 13170–13175, 2007.

Y. F. Zuo, J. Y. Wang, J. H. Chen et al., “A comparison between spontaneous electroencephalographic activities induced by morphine and morphine-related environment in rats,” Brain Research, vol. 1136, no. 1, pp. 88–101, 2007.

A. A. Koronovskii and A. E. Khramov, “Wavelet bicoherence analysis as a method for investigating coherent structures in an electron beam with an overcritical current,” Plasma Physics Reports, vol. 28, no. 8, pp. 666–681, 2002.

X. Li, D. Li, L. J. Voss, and J. W. Sleigh, “The comodulation measure of neuronal oscillations with general harmonic wavelet bicoherence and application to sleep analysis,” NeuroImage, vol. 48, no. 3, pp. 501–514, 2009.

J. M. Palva, S. Monto, S. Kulashekhar, and S. Palva, “Neuronal synchrony reveals working memory networks and predicts individual memory capacity,” Proceedings of the National Academy of Sciences of the United States of America, vol. 107, no. 16, pp. 7580–7585, 2010.

M. Steriade, “Grouping of brain rhythms in corticothalamic systems,” Neuroscience, vol. 137, no. 4, pp. 1087–1106, 2006.

G. G. Knyazev, “Motivation, emotion, and their inhibitory control mirrored in brain oscillations,” Neuroscience and Biobehavioral Reviews, vol. 31, no. 3, pp. 377–395, 2007.

A. K. Engel and P. Fries, “Beta-band oscillations-signalling the status quo?” Current Opinion in Neurobiology, vol. 20, no. 2, pp. 156–165, 2010.

K. K. Hui, J. Liu, N. Makris, et al., “Acupuncture modulates the limbic system and subcortical gray structures of the human brain: evidence from fMRI studies in normal subjects,” Human Brain Mapping, vol. 9, no. 1, pp. 13–25, 2000.

S. D. Hall, I. M. Stanford, N. Yamawaki et al., “The role of GABAergic modulation in motor function related neuronal network activity,” NeuroImage, vol. 56, no. 3, pp. 1506–1510, 2011.

S. D. Muthukumaraswamy, J. F. Myers, S. J. Wilson, et al., “The effects of elevated endogenous GABA levels on movement-related network oscillations,” Neuroimage C, vol. 66, pp. 36–41, 2012.

G. B. Kaplan, K. A. Leite-Morris, M. Joshi, M. H. Shoeb, and R. J. Carey, “Baclofen inhibits opiate-induced conditioned place preference and associated induction of Fos in cortical and limbic regions,” Brain Research, vol. 987, no. 1, pp. 122–125, 2003.

T. C. Chen, Y. Y. Cheng, W. Z. Sun, and B. C. Shyu, “Differential regulation of morphine antinociceptive effects by endogenous enkephalinergic system in the forebrain of mice,” Molecular Pain, vol. 4, article 41, 2008.

J. Du, Z. L. Sun, J. Jia, X. Wang, and X. M. Wang, “High-frequency electro-acupuncture stimulation modulates intracerebral gamma-aminobutyric acid content in rat model of Parkinson's disease,” Sheng Li Xue Bao, vol. 63, no. 4, pp. 305–310, 2011.