Neural Plasticity
Volume 2018 (2018), Article ID 6109723, 15 pages
https://doi.org/10.1155/2018/6109723

Research Article

Electroacupuncture Treatment Alleviates the Remifentanil-Induced Hyperalgesia by Regulating the Activities of the Ventral Posterior Lateral Nucleus of the Thalamus Neurons in Rats

Hong-Yan Zhao,¹ Ling-Yu Liu,¹ Jie Cai,¹ Yan-Jun Cui,² and Guo-Gang Xing1,3,4

¹Neuroscience Research Institute and Department of Neurobiology, School of Basic Medical Sciences, Peking University, Beijing 100083, China
²Department of Internal Medicine, Peking University Hospital, Beijing 100871, China
³The Second Affiliated Hospital of Xinxiang Medical University, Xinxiang, Henan 453002, China
⁴Key Laboratory for Neuroscience, Ministry of Education/National Health and Family Planning Commission, Beijing 100083, China

Correspondence should be addressed to Yan-Jun Cui; cyj@pku.edu.cn and Guo-Gang Xing; ggxing@bjmu.edu.cn

Received 2 May 2018; Accepted 4 September 2018; Published 11 November 2018

Academic Editor: Xavier Navarro

Copyright © 2018 Hong-Yan Zhao 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.

Mechanisms underlying remifentanil- (RF-) induced hyperalgesia, a phenomenon that is generally named as opioid-induced hyperalgesia (OIH), still remain elusive. The ventral posterior lateral nucleus (VPL) of the thalamus, a key relay station for the transmission of nociceptive information to the cerebral cortex, is activated by RF infusion. Electroacupuncture (EA) is an effective method for the treatment of pain. This study aimed to explore the role of VPL in the development of OIH and the effect of EA treatment on OIH in rats. RF was administered to rats via the tail vein for OIH induction. Paw withdrawal threshold (PWT) in response to mechanical stimuli and paw withdrawal latency (PWL) to thermal stimulation were tested in rats for the assessment of mechanical allodynia and thermal hyperalgesia, respectively. Spontaneous neuronal activity and local field potential (LFP) in VPL were recorded in freely moving rats using the in vivo multichannel recording technique. EA at 2 Hz frequency (pulse width 0.6 ms, 1–3 mA) was applied to the bilateral acupoints “Zusanli” (ST.36) and “Sanyinjiao” (SP.6) in rats. The results showed that both the PWT and PWL were significantly decreased after RF infusion to rats. Meanwhile, both the spontaneous neuronal firing rate and the theta band oscillation in VPL LFP were increased on day 3 post-RF infusion, indicating that the VPL may promote the development of RF-induced hyperalgesia by regulating the pain-related cortical activity. Moreover, 2 Hz-EA reversed the RF-induced decrease both in PWT and PWL of rats and also abrogated the RF-induced augmentation of the spontaneous neuronal activity and the power spectral density (PSD) of the theta band oscillation in VPL LFP. These results suggested that 2 Hz-EA attenuates the remifentanil-induced hyperalgesia via reducing the excitability of VPL neurons and the low-frequency (theta band) oscillation in VPL LFP.

1. Introduction

As an ideal ultrashort-acting opioid, remifentanil (RF) is widely used in clinical anesthesia. But after long-time or high-dose infusion, patients tend to appear hyperalgesia after RF withdrawal, a phenomenon that is generally named as opioid-induced hyperalgesia (OIH). The clinical symptoms of OIH are manifested as an increased sensitivity to noxious stimuli, and preclinical experiments have shown that mechanical and thermal hyperalgesia appeared after RF withdrawal [1–3]. How to prevent and treat the remifentanil-induced hyperalgesia and subsequently to investigate its underlying mechanism has become a research hotspot in recent years.

The ventral posterior lateral nucleus (VPL) of the thalamus is an important relay station for the transmission of nociceptive information to the cerebral cortex and sends up a sensory aspect of nociceptive information to the primary somatosensory cortex (SI) and other cortical areas for further processing [4]. In patients with neuropathic pain, the development of pain and pain modulation are related to the power of oscillations (e.g., theta, alpha, and beta oscillations) in local field potentials (LFPs) [5]. In animal models of neuropathic pain, increased spontaneous discharges, enlarged pain receptive fields, and rhythmic oscillatory firing in VPL neurons are observed [6–8]. In rats with central pain after spinal cord injury, neurons in the VPL firing spontaneously at a high rate, with increased after-discharges and evoked responses [9]. Both preclinical and clinical studies have found that VPL and other brain regions can be activated by RF [10–12]. But there has been less evidence for explaining the role of VPL in OIH. We thus first investigated the role of VPL in the development of OIH by in vivo multichannel recording in this study.

Electroacupuncture (EA), a type of acupuncture with electrical stimulation, is a safe and effective method for the treatment of pain and is widely applied in clinical settings [13–15]. Neuroimaging studies have shown that acupuncture or EA applications elicit extensive changes in multiple brain regions [16, 17], which overlapped with most of the neural networks of pain transmission and perception, including the spinal cord, the thalamus, and the cerebral cortex [18]. In addition, EA treatment also inhibits the activity of nociceptive neurons in the VPL [19], and 2 Hz-EA relieves the remifentanil-induced hyperalgesia in rats [20]. We therefore speculated that 2 Hz-EA might alleviate OIH by altering the activity of VPL neurons in rats.

In the present study, we aimed to explore the role and the underlying mechanisms of VPL in the development of RF-induced hyperalgesia in rats by using the in vivo multichannel recording technique. Furthermore, we investigated whether EA treatment could alleviate the RF-induced hyperalgesia by regulating the activity of VPL neurons in rats. We found that both the excitability of VPL neurons and the activity of theta band power spectral density (PSD) in VPL LFP are increased after RF infusion, indicating that the VPL may promote the development of RF-induced hyperalgesia by regulating the pain-related cortical activity. Moreover, we demonstrated that 2 Hz-EA attenuates the remifentanil-induced hyperalgesia via reducing the excitability of VPL neurons and the activity of low-frequency (theta band) oscillation in VPL LFP.

2. Materials and Methods

2.1. Animals

Male Sprague-Dawley rats weighing 180~220 g at the beginning of the experiment were provided by the Department of Experimental Animal Sciences, Peking University Health Science Center. The rats were kept in separated cages, free to get food and water, and kept in the animal feeding box controlled by temperature (24 ± 1°C), humidity (50~55%), and irradiation (12 : 12 h light : dark cycle). All animal experimental procedures were carried out in accordance with the guidelines of the International Association for the Study of Pain [21] and were approved by the Animal Care and Use Committee of Peking University (No. LA2016060, 2016). The behavioural experimenters were kept blind from the groupings of the rats.

2.2. Drug Administration

A 24-gauge disposable catheter (24 G, flow: 23 ml/min, L = 20 mm, Shanghai Puyi Medical Instruments Co. Ltd.) was flushed with heparinized saline and inserted into the caudal vein. The caudal vein catheter was used to infuse normal saline (Datong Huida Pharmaceutical Co. Ltd.) and Ruijie® (remifentanil hydrochloride for injection, H20030197, Yichang Humanwell Pharmaceutical Co. Ltd., China) dissolved in normal saline. Rats were anesthetized by pentobarbital sodium (50 mg·kg⁻¹, i.p.). The rats were assigned randomly to a treatment group that received remifentanil (RF) intravenously or to a control group that received normal saline (vehicle) via the tail vein. The treatment protocol was performed as follows: rats in RF group received remifentanil infusion (1.5 μg·kg⁻¹·min⁻¹), while rats in vehicle group rats received 0.9% normal saline infusion (0.1 ml·kg⁻¹·min⁻¹). The infusion lasted for 2 hours, and spontaneous breaths and eyelash reflexes were monitored to ensure the rats’ wellbeing.

3. Behavioural Tests

3.1. Assessment of Mechanical Allodynia

Mechanical allodynia was assessed by measuring 50% paw withdrawal threshold (PWT) in rats as described in our previous reports [22–25]. In brief, the 50% PWT in response to a series of von Frey filaments (Stoelting, Wood Dale, IL) was determined by the up-and-down method [26]. Rats were placed on a metal mesh floor covered with an inverted clear plastic cage (18 × 8 × 8 cm) and allowed for a habituation period of 20 minutes. Eight von Frey filaments (0.4, 0.7, 1.2, 2.0, 3.6, 5.5, 8.5, and 15.1 g) with approximately equal logarithmic increments (0.224) of bending force were selected. Each test began with a delivery of a 2.0 g von Frey filament force perpendicular to the plantar surface of the hind paw for approximately 2~3 s. An abrupt withdrawal of the foot immediately after stimulation or immediately after removal of the filament was recorded as a positive response. Whenever there was a positive or negative response, the next weaker or stronger filament was applied, respectively. This procedure was completed until the six stimuli after the first change in response were observed. The 50% PWT was calculated using

where is the value of the final von Frey filament used (in log units), is a value measured from the pattern of positive/negative responses, and δ = 0.224, which is the average interval (in log units) between the von Frey filaments [27]. If an animal responded to the lowest filament, a value of 0.25 g was assigned. If an animal did not respond to the highest filament, the value was recorded as 15.0 g. The allodynic rats were defined as 50% PWT which is less than 4.0 g (i.e., withdrawal in response to innoxious tactile stimulation). Testing sessions were performed on day 1 prior to RF infusion as well as on days 1, 3, 5, 7, 9, and 11 post-RF infusion.

3.2. Assessment of Thermal Hyperalgesia

Thermal hyperalgesia of the hind paws was tested according to our previous studies [23, 25, 28, 29]. Briefly, rats were adapted to the acrylic enclosures on a clear glass plate maintained at 30°C for at least 30 minutes. A radiant heat source was focused onto the plantar surface of the hind paw. Measurements of paw withdrawal latency (PWL) were taken by a timer that was started by the activation of the heat source and stopped when withdrawal of the paw was detected with a photo detector. A maximal cut-off time of 30 seconds was applied to avoid unnecessary tissue damage. Three measurements of PWL to each hind paw were taken and averaged as the result of each test session. The hind paw was alternately tested between consecutive tests with intervals greater than 5 minutes. Testing sessions were performed on day 1 prior to RF infusion as well as on days 1, 3, 5, 7, 9, and 11 post-RF infusion.

3.3. Electroacupuncture (EA) Application

Rats were restrained in rodent holders with their hind legs and tails protruding [30]. Briefly, a pair of stainless steel needles (0.2 mm in diameter, 5 mm in length) was inserted into acupoints “Zusanli” (ST.36, 4 mm lateral to the anterior tubercle of the tibia, which is marked by a notch) and “Sanyinjiao” (SP.6, 3 mm proximal to the medial malleolus, at the posterior border of the tibia). Electrical stimulation of a constant current generated from Han’s Acupoint Nerve Stimulator (HANS, LH202H, Beijing Astronautics and Aeronautics Aviation University, Beijing, China) was delivered to bilateral hind limbs simultaneously. This electrical stimulation was set as a square wave of 0.6 ms in pulse width and 2 Hz in frequency. The electrical pulses were delivered from 1 to 3 mA in a step of 1 mA increment. Each step of the stimulus lasted for 10 min. The EA stimulation was given for 30 min on day 1 prior to RF infusion, continued during the 2-hour RF infusion, and three times (once per day) on days 1, 2, and 3 after RF infusion. Animals allocated into the mock EA groups were then placed in the same apparatus and underwent needle insertion in the same acupoints, but no electrical current was applied to them. EA procedures were always carried out by the same experimenter.

4. In Vivo Multichannel Recording

4.1. Surgery

Initial anesthesia was performed on the rats by sodium pentobarbital injection (50 mg/kg, i.p.) prior to implantation of the microelectrode array. A supplemental dose (1/3 of the original dose) of pentobarbital sodium was administered to the animals as needed to maintain proper anesthetic depth during surgery. Rats were mounted on a stereotaxic device (Reward Life Science Technology Co. Ltd., Shenzhen), and an array of eight nickel chromium alloy wire PEG 2000-insulated microwires (14 μm in diameter, arranged in a 4 × 2 configuration, 250 μm spacing between each microwire, STABLOHM 675, California Fine Wire Company, USA) was slowly lowered into the right VPL (−3.0 mm AP, −3.0 mm ML, and −6.0 mm DV) according to the rat atlas of Paxinos and Watson [31]. Six stainless steel screws were driven into the skull to serve as anchors for holding the microwires in place after implantation. Prior to surgery, rats received penicillin injection (16,000 IU, i.m.) to prevent infection, and then the animals were allowed to recover for one week before the start of recording (see Figure 1).

(a)

(b)

Figure 1: In vivo multichannel recording. (a) Surgery for the electrode implantation and the in vivo multichannel recording in freely moving rats. (b) Experimental procedure of in vivo multichannel recording. Data was collected on day 1 before and day 3 after remifentanil infusion.

4.2. Electrophysiological Recording

In vivo multichannel recording was performed as described elsewhere [32]. Briefly, rats were placed in a transparent plastic chamber (44 × 44 × 44 cm) in a quiet room at 22 ± 1°C and allowed to move freely throughout the recording period. Electrophysiological signals were recorded after the animals adapted to the experimental environment. Data were collected at baseline (before infusion) and on day 3 after RF infusion (see Figure 1(b)). Multiunit neuronal activities were collected through the implanted microwire assemblies that were connected to a preamplifier via a headstage plug and a lightweight cable, and the ground wire was used as a reference. The outputs of the preamplifier were filtered (0.5 and 5 kHz, 6 dB cut-off) and sent to a multichannel spike-sorting device (Plexon Inc.) for online signal processing. Local field potentials (LFPs) were recorded from the right VPL by microwire arrays with a multichannel data acquisition system (Plexon Inc.). The LFP signals were transmitted from the headstage assemblies to the preamplifier via a lightweight cable. LFPs were collected at a sampling frequency of 10 kHz, amplified (300x), and band-pass filtered (0.3–500 Hz). LFP signals were filtered into four frequency bands: theta (4–8 Hz), alpha (9–12 Hz), beta (13–30 Hz), and gamma (31–100 Hz), in which the theta, alpha, and beta bands belong to the low-frequency band (4–30 Hz), while the gamma band belongs to the high-frequency band (31–100 Hz). In a transparent plastic chamber, LFP oscillations and the video recording for the animals’ behaviours were simultaneously recorded during a 10-minute resting state. Spike train activity was analysed by NeuroExplorer (Plexon Inc.). Waveform capture and frequency distribution histogram were used to verify the classification of a single unit. Different waveforms were individually distinguished by setting multiple time-voltage windows using Offline Sorter software (Plexon Inc.). The time stamps of these waveforms were then stored on a personal computer for off-line analysis.

4.3. Data Analysis

According to a previous report [32], the neuronal firing rate was quantified for each neuron to construct rate histograms with a time range from the beginning to 5 min after recording in different groups. In brief, the bin size used to calculate the histograms was 2 seconds. Bin counts for each trial were calculated using the analysis program NeuroExplorer (Plexon Inc.), and the results were exported to Matlab (MathWorks Inc.) in spreadsheet form. The firing rates for all neurons were normalized and ranked into a spreadsheet for further statistical analysis. Power spectral density (PSD) was analysed using the NeuroExplorer (Plexon Inc.) and Matlab toolbox Chronux (http://chronux.org). Local field potential (LFPs) data were resampled at 1 kHz and analysed separately for each 5-minute recording period. The 50 Hz frequency was removed from the digitized signal prior to any analysis with a frequency resolution of ~1 Hz. The calculated parameter values were set as follows: Params. fpass = [0–100] (frequency domain for analysis between 0 and 100 Hz); Params. fpass = [4–8] (frequency domain for analysis between 4 and 8 Hz). The time-varying power spectra was then calculated by fast Fourier transformation (FFT). Spectral units were normalized by raw power spectral densities. Power spectral analysis was carried out to calculate PSD.

4.4. Histology

After the termination of the experiment, the rats were deeply anesthetized with pentobarbital sodium and the tip positions of the electrodes were identified by a 2 mA, 10 s DC current (anode current) through the electrodes to produce a thermal lesion of the near area tissue. The animals were then sacrificed and perfused with 0.9% saline followed by 4% paraformaldehyde. After fixing in paraformaldehyde at 4°C overnight, the brains were transferred to a 20% sucrose solution in saline for cryoprotection. Coronal sections of 30 μm were cut on a microtome, mounted on charged slides, and stained with neutral red dye for 10–15 min. Recording sites were determined under a light microscope (see Figure 2). Data of those sites deflecting from the target area were excluded from analysis.

(a)

(b)

Figure 2: Histological verification of the recording electrodes sites in the right ventral posterior lateral nucleus (VPL) of the thalamus region. (a) Photomicrograph of a histological section stained with neutral red to illustrate the position of the recording electrode sites made by an electrolytic lesion through the respective electrodes (as arrow shows). (b) Representations of three coronal sections through the rat VPL are shown in sequence from anterior to posterior, respectively. The numbers in the right margin indicate millimeters posterior to the bregma. The filled circles in the right hemisphere show the approximate positions of the recording electrodes corresponding to some representative rats in the different groups.

4.5. Statistical Analysis

As described in previous reports [24, 33], statistical analyses were performed with GraphPad Prism 7.0 for Windows (GraphPad Software Inc., La Jolla, CA). All data were expressed as mean ± standard error of the mean. Two-way analysis of variance (ANOVA) followed by the Bonferroni post hoc test was used for multiple comparisons. Differences with were considered statistically significant.

5. Results

5.1. Remifentanil Induces Mechanical Allodynia and Thermal Hyperalgesia in Rats

Compared to the vehicle group, infusion of RF (1.5 μg·kg⁻¹·min⁻¹ for 120 minutes) to rats induced a significant decrease both in the paw withdrawal threshold (PWT) in response to mechanical stimuli (on day 3 postinfusion, 13.3 ± 0.6 g vehicle vs. 3.1 ± 0.5 g RF, , two-way ANOVA, vehicle and 13 RF) (see Figure 3(a)) and the paw withdrawal latency (PWL) to thermal stimulation (on day 3 postinfusion, 16.1 ± 0.8 sec vehicle vs. 11.9 ± 0.3 sec RF, , two-way ANOVA, vehicle and 13 RF) (see Figure 3(b)). The results suggested that RF infusion resulted in a stable OIH in rats. Moreover, inclined-plate test was performed to assess the locomotor function of experimental animals. Expectedly, no significant motor dysfunction was observed in rats that received RF infusion (data not shown). In addition, we also performed an elevated plus-maze (EPM) test to determine whether RF infusion could produce an anxiety-like behaviour in rats. The results revealed that there were no significant differences in the time spent in the open arms as well as the number of entries in the open arms between the vehicle- and the RF-treated rats on day 3 after RF infusion (data not shown). These results suggested that RF infusion could not induce anxiety-like behaviour in rats.

(a)

(b)

Figure 3: Effects of RF infusion on the paw withdrawal threshold (PWT) in response to mechanical stimuli and the paw withdrawal latency (PWL) to thermal stimulation in rats. (a) PWT. (b) PWL. Note that RF infusion results in a significant mechanical allodynia and thermal hyperalgesia as measured by a decreased PWT and PWL in rats. , , or versus the corresponding vehicle group, two-way ANOVA followed by the Bonferroni post hoc test, vehicle and 13 RF.

5.2. Remifentanil Enhances the Spontaneous Neuronal Firing Rate and the Activity of Theta Band Oscillation in VPL LFP

The numbers of simultaneously recorded neurons were depicted as follows: 26 neurons pre- and 28 neurons postinfusion in the vehicle group and 30 neurons pre- and 30 neurons postinfusion in the RF group.

To investigate the ongoing neuronal activity in the VPL associated with RF-induced hyperalgesia, we examined the changes in mean spontaneous firing rate of VPL neurons by in vivo multichannel recording in free-moving rats. Compared with the vehicle group, RF infusion enhanced the mean spontaneous firing rate of VPL neurons on day 3 postinfusion (spikes/s, 0.2 ± 0.02 vehicle vs. 0.5 ± 0.06 RF, , two-way ANOVA, to 30 cells (from four to six rats) per group; see Figures 4(a), 4(c)). Similar results were observed on the spontaneous firing density as depicted in the cluster plot within 400 seconds pre- and post-RF infusion in the vehicle- and RF-infusion groups, respectively (see Figure 4(b)). However, as compared with baseline, the vehicle group rats did not show significant changes in the mean spontaneous firing rate of VPL neurons on day 3 post-RF infusion (, two-way ANOVA, to 30 cells (from four to six rats) per group; see Figure 4(a), 4(c)).

(a)

(b)

(c)

Figure 4: Effects of RF infusion on spontaneous neuronal activity in VPL neurons. (a) Representative spontaneous activity of VPL neurons before and after RF infusion in the vehicle- and RF-treated rats. Upper: representative traces of the neuronal discharges. Scale bar: 400 μV, 10 sec. Middle: scatter plot of the spikes. Lower: histograms show the number of action potentials (spikes) per second (bin width, 2 seconds). (b) Cluster plots depicting activity of VPL neurons (rows) before and after RF infusion in the two groups are shown. Warm and cool colors indicate the increase and decrease, respectively, in the spikes of multichannel recorded neurons. (c) Averaged spontaneous activity of VPL neurons (spikes per second) before and after RF infusion in vehicle- and RF-treated rats. versus the corresponding vehicle group, two-way ANOVA followed by the Bonferroni post hoc test, to 30 cells (four to six rats) per group.

To further investigate the time-frequency relationship of theta, alpha, beta, and gamma oscillations of LFP in the VPL, we calculated the time-varying power spectra by FFT (see Figure 5(a)). The spectrum units were normalized by the raw power spectral density (PSD). Then, we performed power spectral analysis to compare the PSD of theta, alpha, beta, and gamma oscillations in the VPL LFP on day 3 post-RF infusion. The results showed that RF infusion induced a significant increase in the low-frequency theta (4–10 Hz) band activity of VPL LFP (W/Hz, 0.08 ± 0.01 vehicle vs. 0.13 ± 0.02 RF, , two-way ANOVA, to 6 rats per group; see Figure 5(b)).

(a)

(b)

Figure 5: Effects of RF infusion on the activity of VPL LFP in rats. (a) An example of the LFP and the time-frequency representation in the vehicle- and RF-treated rats on day 3 post-RF infusion. Representatives of 4~100 Hz (left) and 4~8 Hz (right) frequency bands of LFP from the vehicle and RF groups are shown. The power spectra are normalized. Warm and cool colors indicate the increase and decrease, respectively, in field potential power of VPL LFP. (b) Summary of the power spectral density of theta, alpha, beta, and gamma bands on day 3 after RF infusion. Note that RF infusion induces a significant increase in activity of low-frequency theta (4–10 Hz) band but not alpha, beta, and gamma bands in VPL LFP. versus the corresponding vehicle group, two-way ANOVA followed by the Bonferroni post hoc test, to 6 rats per group.

5.3. EA Treatment Attenuates the Remifentanil-Induced Mechanical Allodynia and Thermal Hyperalgesia in Rats

Furthermore, we investigated the therapeutic effect of EA on OIH in a rat model of RF-induced hyperalgesia. The EA stimulation was performed by using Han’s Acupoint Nerve Stimulator (HANS, LH202H) applied to the bilateral hind limbs of rats at the acupoints Zusanli (ST.36) and Sanyinjiao (SP.6) simultaneously (see Figures 6(a)–6(b)). The EA stimulation was given for 30 min per day before RF infusion, continued in the 2-hour infusion, and three times (once per day) on days 1, 2, and 3 post-RF infusion (see Figure 6(c)). The results showed that compared with the mock EA group, EA at 2 Hz frequency reversed the RF-induced decrease both in the PWT (on day 1 post-RF infusion, 10.5 ± 1.4 g EA vs. 3.9 ± 1.1 g mock EA, ; see Figure 6(d)) and in the PWL (on day 1 post-RF infusion, 13.8 ± 0.7 sec EA vs. 11.7 ± 0.4 sec RF, ; see Figure 6(e)) in rats (two-way ANOVA, to 16 rats per group). These results suggested that the 2 Hz-EA treatment could alleviate the RF-induced hyperalgesia in rats.

(a)

(b)

(c)

(d)

(e)

Figure 6: Effects of 2 Hz-EA treatment on the RF-induced hyperalgesia in rats. (a–b) Equipment and acupoints of EA stimulation ((a) Han’s Acupoint Nerve Stimulator, LH202H; (b) a schematic diagram indicates the location of EA stimulation in rats). (c) Experimental procedure of the EA treatment. The EA stimulation was given for 30 min on the day preinfusion, continued in the 2-hour infusion, and three times (once per day) on days 1, 2, and 3 post-RF infusion. (d–e) Effects of EA treatment on the paw withdrawal threshold (PWT) in response to mechanical stimuli and the paw withdrawal latency (PWL) to thermal stimulation in rats that received RF infusion. Note that 2 Hz-EA treatment reverses the RF-induced decreases both in the PWT and PWL in rats. or versus the corresponding mock EA group, two-way ANOVA followed by the Bonferroni post hoc test, to 16 rats per group.

5.4. EA Treatment Abrogates the RF-Induced Augmentation of the Spontaneous Neuronal Activity and the PSD of Theta Band Oscillation in VPL LFP

The numbers of simultaneously recorded neurons were depicted as follows: 24 neurons pre- and 17 neurons postinfusion in the mock EA group and 20 neurons pre- and 16 neurons postinfusion in the EA group.

Compared with the mock EA group, the 2 Hz-EA treatment could rescue the RF-induced increase in the spontaneous neuronal activity of VPL neurons on day 3 post-RF infusion (spikes/s, 0.4 ± 0.04 EA vs. 0.6 ± 0.04 mock EA , two-way ANOVA, to 24 cells (from four to five rats) per group; see Figures 7(a), 7(c)). Similar results were observed on the spontaneous firing density as depicted in the cluster plot within 400 seconds pre-RF and post-RF infusion in the EA and mock EA groups (see Figure 7(b)).

(a)

(b)

(c)

Figure 7: Effects of 2 Hz-EA treatment on the spontaneous neuronal activity of VPL neurons in rats that received RF infusion. (a) Representative spontaneous activity of VPL neurons before and after RF infusion in the EA- and mock EA-treated rats. Upper: representative traces of the neuronal discharges. Scale bar: 400 μV, 10 sec. Middle: scatter plot of the spikes. Lower: histograms show the number of action potentials (spikes) per second (bin width, 2 seconds). (b) Cluster plot depicting activity of VPL neurons (rows) before and after RF infusion in the EA- and mock EA-treated rats are shown. Warm and cool colors indicate the increase and decrease, respectively, in the spikes of multichannel recorded neurons. (c) Averaged spontaneous activity of VPL neurons (spikes per second) before and after RF infusion in EA- and mock EA-treated rats. versus the corresponding mock EA group, two-way ANOVA followed by the Bonferroni post hoc test, to 24 cells (four to five rats) per group.

Through the power spectrum analysis of the time-frequency spectrogram in the LFP, we found that the 2 Hz-EA treatment also abrogated the RF-induced increase in the power spectral densities (PSD) of theta band oscillation (4–10 Hz) in VPL LFP (W/Hz, 0.06 ± 0.01 EA vs. 0.1 ± 0.04 mock EA, , two-way ANOVA, to 5 rats per group; see Figure 8).

(a)

(b)

Figure 8: Effects of 2 Hz-EA treatment on the activity of VPL LFP in rats that received RF infusion. (a) An example of the LFP and the time-frequency representation in the EA- and mock EA-treated rats on day 3 post-RF infusion. Representatives of 4~100 Hz (left) and 4~8 Hz (right) frequency bands of LFP from the EA- and mock EA-treated rats are shown. The power spectra are normalized. Warm and cool colors indicate the increase and decrease, respectively, in field potential power of VPL LFP. (b) Summary of the power spectral density of theta, alpha, beta, and gamma bands on day 3 after RF infusion in the EA and mock EA groups. Note that 2 Hz-EA treatment abrogates the RF-induced increase in the power spectral density (PSD) of theta band oscillation (4–10 Hz) in rats on day 3 post-RF infusion. versus the corresponding mock EA group, two-way ANOVA followed by the Bonferroni post hoc test, to 5 rats per group.

6. Discussion

6.1. Remifentanil-Induced Hyperalgesia in Rats

In order to study the underlying mechanisms of OIH, we first developed an animal model of RF-induced hyperalgesia in rats. In line with previous findings [10], we found that the RF-induced hyperalgesia was associated with the dose and duration of remifentanil infusion. Clinically, remifentanil is usually administered intravenously to patients at a dose of 0.1~0.3 μg·kg⁻¹·min⁻¹. Since the equivalent dose in rats is 6.259 times of that in humans, the corresponding dose in rats is approximately 0.625~1.875 μg·kg⁻¹·min⁻¹. A study of mice showed that the half effective doses (ED₅₀) of RF-induced thermal and mechanical hyperalgesia were 1.7 μg·kg⁻¹·min⁻¹ (95% confidence interval: 1.3~2.1) and 1.26 μg·kg⁻¹·min⁻¹ (1.0~1.6), respectively [34]. Since a larger dose of RF can induce a more significant OIH phenomenon [35], we thus chose a dose of 1.5 μg·kg⁻¹·min⁻¹ to develop a rat model of RF-induced hyperalgesia. A previous study has shown that the RF-induced hyperalgesia requires sufficient drug exposure time, and the duration of RF exposure is the main factor influencing the hyperalgesia induction [10]. Other studies also demonstrated that after 90 min of continuous infusion of RF, the area of mechanical hyperalgesia on the skin of rats is significantly enlarged, and more than 2 hour duration of infusion increases the incidence of RF-induced hyperalgesia [36, 37]. Based on our preliminary data, a long-term infusion (3 hours) of RF could result in severe respiratory depression, we thus developed a rat model of RF-induced hyperalgesia by infusing RF at 1.5 μg·kg⁻¹·min⁻¹ for 2 hours in this study. In line with previous findings [1, 20], our present study revealed that after a 2-hour continuous infusion of RF, the rats showed significant mechanical allodynia and thermal hyperalgesia that lasted for almost one week. Moreover, we also observed that the postoperative hyperalgesia caused by RF infusion in a rat model of plantar incision could last for at least two weeks (data not shown), which far exceeded the recovery period of incisional pain only (within a week). Apart from the pain sensitivity, we also have tested the locomotor function and the anxiety-like behaviour in rats after RF infusion. We found that no significant motor dysfunction and anxiety-like behaviour appeared in rats on day 3 post-RF infusion. We speculated that it is possibly due to the time of our observation was relatively earlier after RF infusion; although the hyperalgesia was most severe on day 3 post-RF infusion, it may have not caused the pain-related negative emotions yet.

6.2. Remifentanil Enhances the Excitability of VPL Neurons in Rats

The thalamus is a gateway for the transmission of sensory information from the spinal cord to the cerebral cortex, where the VPL is projected to the primary somatosensory cortex (SI), secondary somatosensory cortex (SII), and other cortical areas for further processing. After the systematic infusion of RF in human volunteers, the functional magnetic resonance imaging (fMRI) results showed that the cerebral blood flow in the brain regions such as the anterior cingulate cortex (ACC) and the VPL was increased significantly, suggesting that these brain regions can be activated by RF infusion [38, 39]. In addition, systematic administration of morphine could increase the spontaneous firing rate of the VPL neurons [40]. Since the SI receive dense projections from the VPL, a major target for ascending sensory tracts in the spinal cord and brain stem, a large number of neurons in the SI are activated in a rat model of RF-induced hyperalgesia [41]. Our present results showed that the pain-related spontaneous firing rate of VPL neurons was increased in RF-treated rats. These data indicated that the increased excitability of VPL neurons may lead to a hyperexcitability in the pain-related brain regions such as SI and SII, thereby promoting the development of RF-induced hyperalgesia. In fact, evidence from human studies has revealed that the thalamocortical dysrhythmia and altered thalamic burst firing play a vital role in neuropathic pain [42, 43]. In patients with neuropathic pain, the thalamic neurons fire low-threshold calcium spike (LTS) bursts, which are thought to be mediated by T-type calcium channels [44]. These abnormal bursts are critical for the generation of theta rhythms during sleep [45, 46] as well as delta and theta oscillations in thalamic LFPs [44]. Hence, the LTS bursts of thalamic relay neurons exert a rhythmic influence on thalamocortical modules in the theta frequency band, which may underlie the pathophysiology of RF-induced hyperalgesia or chronic neuropathic pain [47].

6.3. Remifentanil Increases the Activity of Theta frequency Oscillation in VPL LFP

Brain rhythmical oscillations in the low (delta, theta, and alpha) and high (beta and gamma) frequencies of electroencephalography are linked to broad varieties of perceptual, sensorimotor, and cognitive operations [48]. For example, oscillatory brain activities in the high-frequency gamma band play an important role in selecting and integrating sensory-relevant information into a coherent perception [49]. Accumulative evidence has demonstrated that gamma band oscillations also participate in pain perception [50–52] and the abnormal working memory of chronic pain [53]. An increased gamma band oscillatory activity was also observed in the frontal-central region in subjects suffering from tonic muscle pain [52] or tonic heat pain [54]. The beta band oscillation may be involved in the maintenance of the current sensorimotor or cognitive state, and the extraordinary enhancement of the beta oscillation may result in an abnormal persistence of the current situation and a deterioration of flexible behaviours and cognitive controls [55]. With respect to the low-frequency oscillations, the delta band oscillation is associated with compromised neuronal function [56], the alpha band oscillation mainly serves as a top-down controlled inhibitory mechanism [57], and the theta band oscillation promotes neural plasticity [58] and participates in the process of information processing for touch and pain [59, 60]. Theta oscillation also plays a crucial network-level role in hippocampal learning and memory; thus, suppression of the theta rhythm impairs learning and memory [61]. These findings support our understanding for the role of low-frequency oscillations in pain processing [62].

Moreover, changes in wide-range band oscillations have been seen in several acute and chronic pain models. Under noxious stimulation, decreased power of alpha and beta bands and increased power of the gamma band have been reported [63, 64]. In a separate study, the authors found that the power of delta, gamma, and epsilon bands was generally increased whereas the power of alpha and beta bands was decreased in rats that received noxious laser heat stimulation [65]. In neuropathic rats, acute mechanical pain stimulation leads to increased alpha and decreased beta and gamma oscillations in the medial prefrontal cortex (mPFC) LFP, while the delta oscillation is decreased and the gamma oscillation is dynamically decreased during chronic neuropathic pain condition [66]. Additionally, in freely moving rats with carrageenan-induced inflammation, application of von Frey mechanical stimulation to the ipsilateral paw results in a significant increase in delta, theta, and alpha bands of LFP in the ACC, suggesting that there are significant changes in the low-frequency oscillations of LFP occurred during inflammatory pain [67]. So, an increased theta band oscillation of VPL LFP in RF-treated rats may suggest the role of low-frequency oscillation, particularly in the VPL, in RF-induced hyperalgesia.

In addition, cross-frequency coupling has functional significance in cortical information processing, thereby contributing to the integration of relevant regions in the brain [68]. In human beings, painful laser stimulation enhance gamma oscillatory responses in the pain-related amygdala and hippocampal regions, and these gamma responses are significantly coupled with the phases of theta and alpha rhythms during pain processing [69]. Painful cutaneous laser stimuli also induce event-related gamma band activity in the lateral thalamus in humans, and thalamic cross-frequency coupling analysis indicates that the phase of the lower frequency activity (theta to beta) modulates the amplitude of the higher frequency activity (low and high gamma) more strongly during the cutaneous laser stimuli [70]. In rats with visceral pain, an increased theta band power is observed in the ACC, and the neuronal spike activity in the ACC becomes synchronized with ongoing theta oscillations of LFP [71]. Moreover, an enhanced synchronization of thalamo-ACC theta band LFPs was consistent with an increased neuronal communication between the two regions [71]. These results reveal that theta band oscillations and theta frequency phase-locking may serve as prominent features of neuronal activity in the ACC and a potential neural mechanism underlying acute visceral pain. Consistently, our present data showed that the activity of the theta band oscillation in VPL LFP was significantly increased on day 3 post-RF infusion. These data suggest that the VPL may promote the development of RF-induced hyperalgesia by regulating the pain-related cortical areas.

It is accepted that thalamocortical oscillations are critical for sensory perception [72], and pain is known to disrupt synchrony in thalamocortical oscillations. Thalamocortical coherence has been reported to be increased in patients with neurogenic pain [42] but decreased in a rat model of central pain [73]. In a separate study, the authors also reported that the thalamocortical coherence between VPL and SI is decreased while the theta power amplitude in the thalamus and cortex is increased in rat models of acute and chronic pain [74]. It has been suggested that the thalamus regulates and coordinates with the cerebral cortex through direct anatomical connections rather than distal synapses and that the low-frequency theta rhythm in the thalamus can directly regulate high-frequency cortical activity in the beta or higher range [75, 76]. Also, the clinical studies have found that the power of the theta band in the VPL was negatively correlated with pain relief; that is, the power of the theta band was increased during the onset of pain [8]. In agreement with these findings, our present results showed that the power spectral density (PSD) in the theta band of VPL LFP was significantly increased after RF infusion, which support our hypothesis that the VPL may promote the development of RF-induced hyperalgesia by coupling with high-frequency oscillatory activity in other pain-related cortical areas.

6.4. The Analgesic Effects of EA on the Remifentanil-Induced Hyperalgesia in Rats

The frequency is one of the important parameters for EA treatment. Low-frequency (2 Hz) EA was effective for the treatment of inflammatory pain and neuropathic pain [13, 77], while high-frequency (100 Hz) EA was better for the spinal injury-induced muscle crampy pain [78]. In this study, we validated these findings by showing that 2 Hz-EA indeed has a therapeutic effect on RF-induced hyperalgesia in rats. Additionally, we found that the current intensity is another key parameter for EA treatment. Several studies have shown that the low-intensity (1 mA) EA is superior to the high-intensity (5 mA) EA for the analgesic effect of 2 Hz-EA on the CCI-induced neuropathic pain. For patients with postoperative pain, transcutaneous electrical acupoint stimulation (TEAS) with strong current intensity is more effective [79]. We have tested the analgesic effect of 2 Hz-EA on RF-induced hyperalgesia with low-intensity (0.5–1.5 mA) and high-intensity (1–3 mA) EA in our preliminary experiments and found that the high-intensity EA could effectively alleviate the RF-induced hyperalgesia. Moreover, acupoints selection is also an important factor influencing the efficacy of EA treatment. Zusanli (ST.36) and Sanyinjiao (SP.6) are usually applied to produce the analgesic effect [14]. In this study, we found that 2 Hz-EA had a significant analgesic effect on the RF-induced hyperalgesia in rats by stimulating bilateral ST.36 and SP.6 acupoints.

6.5. Mechanisms Underlying the Analgesic Effects of EA on the Remifentanil-Induced Hyperalgesia

It is reported that EA treatment inhibits the activity of nociceptive neurons in the VPL [19]. Studies have found that 94% of the VPL neurons can be activated by innoxious and noxious stimuli in peripheral receptive fields, while 6% of the VPL neurons only respond to noxious stimuli, and there is no neuron in the VPL that only respond to innoxious stimuli [80]. In line with these findings, we discovered that 2 Hz-EA can significantly inhibit the increased spontaneous firing rate of VPL neurons in rats with RF-induced hyperalgesia. We thus assumed that these inhibitory effects of 2 Hz-EA on spontaneous activity of VPL neurons should be, at least in part, responsible for the analgesic effects of EA on RF-induced hyperalgesia.

Analysis of LFP in the VPL showed that 2 Hz-EA also abrogates the increased power spectral density (PSD) of the theta band LFP in the VPL in RF-treated rats. Our results are consistent with previous findings showing that the low-frequency (2–10 Hz) stimulation of the thalamus reduces the theta band oscillation of the cortex and the nociceptive behaviours [81–83]. As mentioned above, thalamocortical oscillations are critical for sensory perception [72] and that the low-frequency theta rhythm in the thalamus can directly regulate high-frequency cortical activity in the beta, gamma, or higher ranges [75, 76]. Therefore, coupling between theta and gamma is involved in various acute and chronic pain processing [66, 71]. On the other hand, some clinical and preclinical studies have demonstrated that the power of the theta band oscillations in the VPL is negatively correlated with pain relief and EA at 2 Hz frequency reverses the increased beta power and the cross-frequency couplings between beta and low-frequency theta band oscillations induced by postincisional pain [8], implying that EA regulates the neuronal networks involved in the central processing and in the integration of spontaneous nociceptive information. Together these findings with our present results, we speculated that EA may serve as an “antagonist” for theta band oscillation of VPL LFP and be able to disrupt the cross-frequency couplings between low-frequency theta band and high-frequency cortical activities. Thus, the analgesic effects of 2 Hz-EA on the RF-induced hyperalgesia is likely mediated by suppressing the low-frequency (theta band) oscillation of LFP in the VPL.

7. Conclusions

In summary, our present results demonstrated that a 2-hour remifentanil infusion induces a substantial pain hypersensitivity in rats and also enhances the spontaneous neuronal firing rate and the activity of the theta band LFP in the VPL. 2 Hz-EA attenuates the remifentanil-induced hyperalgesia via inhibiting the excitability of VPL neurons and the activity of low-frequency theta band oscillation in VPL LFP. The analgesic effects of EA on RF-induced hyperalgesia through its "antagonistic" role in theta band oscillation of VPL LFP also provides an extra evidence for the contribution of neural oscillation to the pathogenesis of OIH.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors have no conflicts of interest to declare.

Acknowledgments

The present work was supported by grants from the National Natural Science Foundation of China (81671085, 81371237, and 61527815), the National Basic Research Program of China (973 Program, 2013CB531905), and the open project for the basic research program from the Institute of Acupuncture and Moxibustion, China Academy of Chinese Medical Sciences (ZZ04004).

References

Y. Z. Li, X. H. Tang, C. Y. Wang et al., “Glycogen synthase kinase-3β inhibition prevents remifentanil-induced postoperative hyperalgesia via regulating the expression and function of AMPA receptors,” Anesthesia and Analgesia, vol. 119, no. 4, pp. 978–987, 2014.

S. H. Kim, N. Stoicea, S. Soghomonyan, and S. D. Bergese, “Intraoperative use of remifentanil and opioid induced hyperalgesia/acute opioid tolerance: systematic review,” Frontiers in Pharmacology, vol. 5, p. 108, 2014.

C. Pasero and M. McCaffery, “Opioid-induced hyperalgesia,” Journal of Perianesthesia Nursing, vol. 27, no. 1, pp. 46–50, 2012.

C. T. Yen and P. L. Lu, “Thalamus and pain,” Acta Anaesthesiologica Taiwanica, vol. 51, no. 2, pp. 73–80, 2013.

H. Luo, Y. Huang, X. du et al., “Dynamic neural state identification in deep brain local field potentials of neuropathic pain,” Frontiers in Neuroscience, vol. 12, p. 237, 2018.

B. C. Hains, C. Y. Saab, and S. G. Waxman, “Alterations in burst firing of thalamic VPL neurons and reversal by Na (v)1.3 antisense after spinal cord injury,” Journal of Neurophysiology, vol. 95, no. 6, pp. 3343–3352, 2006.

B. C. Hains, C. Y. Saab, and S. G. Waxman, “Changes in electrophysiological properties and sodium channel Nav 1.3 expression in thalamic neurons after spinal cord injury,” Brain, vol. 128, no. 10, pp. 2359–2371, 2005.

Y. Huang, H. Luo, A. L. Green, T. Z. Aziz, and S. Wang, “Characteristics of local field potentials correlate with pain relief by deep brain stimulation,” Clinical Neurophysiology, vol. 127, no. 7, pp. 2573–2580, 2016.

M. B. Gerke, A. W. Duggan, L. Xu, and P. J. Siddall, “Thalamic neuronal activity in rats with mechanical allodynia following contusive spinal cord injury,” Neuroscience, vol. 117, no. 3, pp. 715–722, 2003.

R. Ishida, T. Nikai, T. Hashimoto, T. Tsumori, and Y. Saito, “Intravenous infusion of remifentanil induces transient withdrawal hyperalgesia depending on administration duration in rats,” Anesthesia and Analgesia, vol. 114, no. 1, pp. 224–229, 2012.

C. H. Liu, D. N. Greve, G. Dai, J. J. A. Marota, and J. B. Mandeville, “Remifentanil administration reveals biphasic phMRI temporal responses in rat consistent with dynamic receptor regulation,” NeuroImage, vol. 34, no. 3, pp. 1042–1053, 2007.

V. Fodale, D. Schifilliti, C. Praticò, and L. B. Santamaria, “Remifentanil and the brain,” Acta Anaesthesiologica Scandinavica, vol. 52, no. 3, pp. 319–326, 2008.

R. Zhang, L. Lao, K. Ren, and B. M. Berman, “Mechanisms of acupuncture-electroacupuncture on persistent pain,” Anesthesiology, vol. 120, no. 2, pp. 482–503, 2014.

Z. Q. Zhao, “Neural mechanism underlying acupuncture analgesia,” Progress in Neurobiology, vol. 85, no. 4, pp. 355–375, 2008.

Z. Y. Ju, K. Wang, H. S. Cui et al., “Acupuncture for neuropathic pain in adults,” The Cochrane Database of Systematic Reviews, vol. 12, 2017.

R. P. Dhond, N. Kettner, and V. Napadow, “Neuroimaging acupuncture effects in the human brain,” The Journal of Alternative and Complementary Medicine, vol. 13, no. 6, pp. 603–616, 2007.

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, 319 pages, 2005.

G. Biella, M. L. Sotgiu, G. Pellegata, E. Paulesu, I. Castiglioni, and F. Fazio, “Acupuncture produces central activations in pain regions,” NeuroImage, vol. 14, no. 1, pp. 60–66, 2001.

P. J. Rong, J. J. Zhao, L. L. Yu et al., “Function of nucleus ventralis posterior lateralis thalami in acupoint sensitization phenomena,” Evidence-Based Complementary and Alternative Medicine, vol. 2015, Article ID 516851, 6 pages, 2015.

C. Shi, Y. Liu, W. Zhang et al., “Intraoperative electroacupuncture relieves remifentanil-induced postoperative hyperalgesia via inhibiting spinal glial activation in rats,” Molecular Pain, vol. 13, pp. 174480691772563–174480691772229, 2017.

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

M. J. Li, L. Y. Liu, L. Chen, J. Cai, Y. Wan, and G. G. Xing, “Chronic stress exacerbates neuropathic pain via the integration of stress-affect-related information with nociceptive information in the central nucleus of the amygdala,” Pain, vol. 158, no. 4, pp. 717–739, 2017.

D. Fang, L. Y. Kong, J. Cai et al., “Interleukin-6-mediated functional up-regulation of TRPV1 receptors in dorsal root ganglion neurons through the activation of JAK/PI3K signaling pathway: roles in the development of bone cancer pain in a rat model,” Pain, vol. 156, no. 6, pp. 1–1144, 2015.

H. Jiang, D. Fang, L. Y. Kong et al., “Sensitization of neurons in the central nucleus of the amygdala via the decreased GABAergic inhibition contributes to the development of neuropathic pain-related anxiety-like behaviors in rats,” Molecular Brain, vol. 7, no. 1, p. 72, 2014.

Q. Zheng, D. Fang, J. Cai, Y. Wan, J. S. Han, and G. G. Xing, “Enhanced excitability of small dorsal root ganglion neurons in rats with bone cancer pain,” Molecular Pain, vol. 8, p. 1744-8069-8-24, 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, no. 1, pp. 441–462, 1980.

M. Liu, H. Yang, D. Fang et al., “Upregulation of P2X3 receptors by neuronal calcium sensor protein VILIP-1 in dorsal root ganglions contributes to the bone cancer pain in rats,” Pain, vol. 154, no. 9, pp. 1551–1568, 2013.

Q. Zheng, D. Fang, M. Liu et al., “Suppression of KCNQ/M (Kv7) potassium channels in dorsal root ganglion neurons contributes to the development of bone cancer pain in a rat model,” Pain, vol. 154, no. 3, pp. 434–448, 2013.

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, pp. 5–13, 2001.

G. Paxinos and C. Watson, The Rat Brain in Stereotaxic Coordinates, Elsevier Academic Press, San Diego, 3th edition, 2005.

N. Wang, Y. Zhang, J. Y. Wang, G. Gao, and F. Luo, “Effects of pentobarbital anesthesia on nociceptive processing in the medial and lateral pain pathways in rats,” Neuroscience Bulletin, vol. 26, no. 3, pp. 188–196, 2010.

L. Chen, S. Li, J. Cai et al., “Activation of CRF/CRFR1 signaling in the basolateral nucleus of the amygdala contributes to chronic forced swim-induced depressive-like behaviors in rats,” Behavioural Brain Research, vol. 338, pp. 134–142, 2018.

W. Cui, Y. Li, S. Li et al., “Systemic lidocaine inhibits remifentanil-induced hyperalgesia via the inhibition of cPKCgamma membrane translocation in spinal dorsal horn of rats,” Journal of Neurosurgical Anesthesiology, vol. 21, no. 4, pp. 318–325, 2009.

C. J. Hayhurst and M. E. Durieux, “Differential opioid tolerance and opioid-induced hyperalgesia: a clinical reality,” Anesthesiology, vol. 124, no. 2, pp. 483–488, 2016.

J. P. Laulin, P. Maurette, J. B. Corcuff, C. Rivat, M. Chauvin, and G. Simonnet, “The role of ketamine in preventing fentanyl-induced hyperalgesia and subsequent acute morphine tolerance,” Anesthesia & Analgesia, vol. 94, no. 5, pp. 1263–1269, 2002.

Z. D. Cooper, Y. N. T. Truong, Y. G. Shi, and J. H. Woods, “Morphine deprivation increases self-administration of the fast- and short-acting mu-opioid receptor agonist remifentanil in the rat,” The Journal of Pharmacology and Experimental Therapeutics, vol. 326, no. 3, pp. 920–929, 2008.

Y. L. Su, J. Huang, N. Wang, J. Y. Wang, and F. Luo, “The effects of morphine on basal neuronal activities in the lateral and medial pain pathways,” Neuroscience Letters, vol. 525, no. 2, pp. 173–178, 2012.

M. Leppä, A. Korvenoja, S. Carlson et al., “Acute opioid effects on human brain as revealed by functional magnetic resonance imaging,” NeuroImage, vol. 31, no. 2, pp. 661–669, 2006.

A. A. Abdul Aziz, D. P. Finn, R. Mason, and V. Chapman, “Comparison of responses of ventral posterolateral and posterior complex thalamic neurons in naive rats and rats with hindpaw inflammation: mu-opioid receptor mediated inhibitions,” Neuropharmacology, vol. 48, no. 4, pp. 607–616, 2005.

Q. Wang, X. Zhao, S. Li, S. Han, Z. Peng, and J. Li, “Phosphorylated CaMKII levels increase in rat central nervous system after large-dose intravenous remifentanil,” Medical Science Monitor Basic Research, vol. 19, pp. 118–125, 2013.

J. Sarnthein and D. Jeanmonod, “High thalamocortical theta coherence in patients with neurogenic pain,” NeuroImage, vol. 39, no. 4, pp. 1910–1917, 2008.

E. G. Jones, “Thalamocortical dysrhythmia and chronic pain,” Pain, vol. 150, no. 1, pp. 4-5, 2010.

D. Jeanmonod, M. Magnin, and A. Morel, “Low-threshold calcium spike bursts in the human thalamus. Common physiopathology for sensory, motor and limbic positive symptoms,” Brain, vol. 119, no. 2, pp. 363–375, 1996.

M. Deschenes, M. Paradis, J. P. Roy, and M. Steriade, “Electrophysiology of neurons of lateral thalamic nuclei in cat: resting properties and burst discharges,” Journal of Neurophysiology, vol. 51, no. 6, pp. 1196–1219, 1984.

V. Crunelli, D. Cope, and S. Hughes, “Thalamic T-type Ca²⁺ channels and NREM sleep,” Cell Calcium, vol. 40, no. 2, pp. 175–190, 2006.

Y. Huang, A. L. Green, J. Hyam, J. Fitzgerald, T. Z. Aziz, and S. Wang, “Oscillatory neural representations in the sensory thalamus predict neuropathic pain relief by deep brain stimulation,” Neurobiology of Disease, vol. 109, Part A, pp. 117–126, 2018.

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.

C. S. Herrmann, M. H. J. Munk, and A. K. Engel, “Cognitive functions of gamma-band activity: memory match and utilization,” Trends in Cognitive Sciences, vol. 8, no. 8, pp. 347–355, 2004.

Z. G. Zhang, L. Hu, Y. S. Hung, A. Mouraux, and G. D. Iannetti, “Gamma-band oscillations in the primary somatosensory cortex—a direct and obligatory correlate of subjective pain intensity,” The Journal of Neuroscience, vol. 32, no. 22, pp. 7429–7438, 2012.

E. Schulz, E. S. May, M. Postorino et al., “Prefrontal gamma oscillations encode tonic pain in humans,” Cerebral Cortex, vol. 25, no. 11, pp. 4407–4414, 2015.

L. Li, X. Liu, C. Cai et al., “Changes of gamma-band oscillatory activity to tonic muscle pain,” Neuroscience Letters, vol. 627, pp. 126–131, 2016.

H. Cardoso-Cruz, M. Sousa, J. B. Vieira, D. Lima, and V. Galhardo, “Prefrontal cortex and mediodorsal thalamus reduced connectivity is associated with spatial working memory impairment in rats with inflammatory pain,” Pain, vol. 154, no. 11, pp. 2397–2406, 2013.

W. Peng, L. Hu, Z. Zhang, and Y. Hu, “Changes of spontaneous oscillatory activity to tonic heat pain,” PLoS One, vol. 9, no. 3, article e91052, 2014.

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.

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.

T. Masquelier, E. Hugues, G. Deco, and S. J. Thorpe, “Oscillations, phase-of-firing coding, and spike timing-dependent plasticity: an efficient learning scheme,” The Journal of Neuroscience, vol. 29, no. 43, pp. 13484–13493, 2009.

G. Michail, C. Dresel, V. Witkovský, A. Stankewitz, and E. Schulz, “Neuronal oscillations in various frequency bands differ between pain and touch,” Frontiers in Human Neuroscience, vol. 10, p. 182, 2016.

E. Schulz, L. Tiemann, T. Schuster, J. Gross, and M. Ploner, “Neurophysiological coding of traits and states in the perception of pain,” Cerebral Cortex, vol. 21, no. 10, pp. 2408–2414, 2011.

N. McNaughton, M. Ruan, and M. A. Woodnorth, “Restoring theta-like rhythmicity in rats restores initial learning in the Morris water maze,” Hippocampus, vol. 16, no. 12, pp. 1102–1110, 2006.

P. Fries, J. H. Schröder, P. R. Roelfsema, W. Singer, and A. K. Engel, “Oscillatory neuronal synchronization in primary visual cortex as a correlate of stimulus selection,” The Journal of Neuroscience: The Official Journal of the Society for Neuroscience, vol. 22, no. 9, pp. 3739–3754, 2002.

M. Ploner, J. Gross, L. Timmermann, B. Pollok, and A. Schnitzler, “Pain suppresses spontaneous brain rhythms,” Cerebral Cortex, vol. 16, no. 4, pp. 537–540, 2006.

M. Hauck, C. Domnick, J. Lorenz, C. Gerloff, and A. K. Engel, “Top-down and bottom-up modulation of pain-induced oscillations,” Frontiers in Human Neuroscience, vol. 9, p. 375, 2015.

X. Li, Z. Zhao, J. Ma et al., “Extracting neural oscillation signatures of laser-induced nociception in pain-related regions in rats,” Frontiers in Neural Circuits, vol. 11, p. 71, 2017.

B. Fu, S. N. Wen, B. Wang, K. Wang, J. Y. Zhang, and S. J. Liu, “Acute and chronic pain affects local field potential of the medial prefrontal cortex in different band neural oscillations,” Molecular Pain, vol. 14, p. 174480691878568, 2018.

A. L. Harris-Bozer and Y. B. Peng, “Inflammatory pain by carrageenan recruits low-frequency local field potential changes in the anterior cingulate cortex,” Neuroscience Letters, vol. 632, pp. 8–14, 2016.

R. T. Canolty and R. T. Knight, “The functional role of cross-frequency coupling,” Trends in Cognitive Sciences, vol. 14, no. 11, pp. 506–515, 2010.

C. C. Liu, J. H. Chien, J. H. Kim et al., “Cross-frequency coupling in deep brain structures upon processing the painful sensory inputs,” Neuroscience, vol. 303, pp. 412–421, 2015.

J. H. Kim, J. H. Chien, C. C. Liu, and F. A. Lenz, “Painful cutaneous laser stimuli induce event-related gamma-band activity in the lateral thalamus of humans,” Journal of Neurophysiology, vol. 113, no. 5, pp. 1564–1573, 2015.

J. Wang, B. Cao, T. R. Yu et al., “Theta-frequency phase-locking of single anterior cingulate cortex neurons and synchronization with the medial thalamus are modulated by visceral noxious stimulation in rats,” Neuroscience, vol. 298, pp. 200–210, 2015.

M. Steriade, “Central core modulation of spontaneous oscillations and sensory transmission in thalamocortical systems,” Current Opinion in Neurobiology, vol. 3, no. 4, pp. 619–625, 1993.

D. A. Seminowicz, L. Jiang, Y. Ji, S. Xu, R. P. Gullapalli, and R. Masri, “Thalamocortical asynchrony in conditions of spinal cord injury pain in rats,” The Journal of Neuroscience, vol. 32, no. 45, pp. 15843–15848, 2012.

B. W. LeBlanc, T. R. Lii, A. E. Silverman, R. T. Alleyne, and C. Y. Saab, “Cortical theta is increased while thalamocortical coherence is decreased in rat models of acute and chronic pain,” Pain, vol. 155, no. 4, pp. 773–782, 2014.

Y. B. Saalmann, M. A. Pinsk, L. Wang, X. Li, and S. Kastner, “The pulvinar regulates information transmission between cortical areas based on attention demands,” Science (New York, N.Y.), vol. 337, no. 6095, pp. 753–756, 2012.

M. Malekmohammadi, W. J. Elias, and N. Pouratian, “Human thalamus regulates cortical activity via spatially specific and structurally constrained phase-amplitude coupling,” Cerebral Cortex, vol. 25, no. 6, pp. 1618–1628, 2015.

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.

H. W. Dong, L. H. Wang, M. Zhang, and J. S. Han, “Decreased dynorphin A (1-17) in the spinal cord of spastic rats after the compressive injury,” Brain Research Bulletin, vol. 67, no. 3, pp. 189–195, 2005.

B. Wang, J. Tang, P. F. White et al., “Effect of the intensity of transcutaneous acupoint electrical stimulation on the postoperative analgesic requirement,” Anesthesia and Analgesia, vol. 85, no. 2, pp. 406–413, 1997.

S. W. Yang, K. A. Follett, J. G. Piper, and T. J. Ness, “The effect of morphine on responses of nucleus ventroposterolateralis neurons to colorectal distension in the rat,” Brain Research Bulletin, vol. 48, no. 6, pp. 609–614, 1999.

J. Liu, H. J. Lee, A. J. Weitz et al., “Frequency-selective control of cortical and subcortical networks by central thalamus,” Elife, vol. 4, article e09215, 2015.

C. Y. Saab and L. F. Barrett, “Thalamic bursts and the epic pain model,” Frontiers in Computational Neuroscience, vol. 10, 2017.

B. W. LeBlanc, B. Cross, K. A. Smith et al., “Thalamic bursts down-regulate cortical theta and nociceptive behavior,” Scientific Reports, vol. 7, no. 1, p. 2482, 2017.