- Open Access
Plasticity changes in forebrain activity and functional connectivity during neuropathic pain development in rats with sciatic spared nerve injury
Molecular Brainvolume 11, Article number: 55 (2018)
Neuropathic pain is a major worldwide health problem. Although central sensitization has been reported in well-established neuropathic conditions, information on the acute brain activation patterns in response to peripheral nerve injury is lacking. This study first mapped the brain activity in rats immediately following spared nerve injury (SNI) of the sciatic nerve. Using blood-oxygenation-level-dependent functional magnetic resonance imaging (BOLD-fMRI), we observed sustained activation in the bilateral insular cortices (ICs), primary somatosensory cortex (S1), and cingulate cortex. Second, this study sought to link this sustained activation pattern with brain sensitization. Using manganese-enhanced magnetic resonance imaging (MEMRI), we observed enhanced activity in the ipsilateral anterior IC (AIC) in free-moving SNI rats on Days 1 and 8 post-SNI. Furthermore, enhanced functional connectivity between the ipsilateral AIC, bilateral rostral AIC, and S1 was observed on Day 8 post-SNI. Chronic electrophysiological recording experiments were conducted to confirm the tonic neuronal activation in selected brain regions. Our data provide evidence of tonic activation-dependent brain sensitization during neuropathic pain development and offer evidence that the plasticity changes in the IC and S1 may contribute to neuropathic pain development.
Chronic pain is a major health problem affecting up to 20% of the global population. Although acute pain can be properly managed, most people with chronic pain do not have access to adequate pain relief . Among the most difficult cases are those with neuropathic pain initiated by primary lesions or dysfunction in the somatosensory nervous system.
Many studies have investigated brain sensitization in well-established chronic neuropathic conditions. Higher tonic activity or sensitization to peripheral and central stimuli have been detected in many brain regions of humans with such conditions [2,3,4] and in animal models [5, 6]. Central sensitization at the spinal level was reported as activation-dependent and possibly induced by repeated stimulation [7,8,9,10]. In addition, some brain areas have been reported as having similar activation-dependent sensitization properties [11, 12]. By investigating the responses to nerve injury in the peripheral and central brain areas, previous studies have detected the early onset of ectopic discharge in injured nerve fibers [13,14,15]. One study of trigeminal neuropathic pain in rats showed that nerve injury induced the long-term enhancement of activity in trigeminal ganglion neurons and thalamic neurons . Although strong and continuous brain activation is a possible cause of brain sensitization during neuropathic pain development, experimental data are lacking.
Practically, the study of brain reactions to severe peripheral nerve injury in human subjects is difficult. Therefore, human subjects for such studies are usually selected from people who have received a diagnosis of chronic pain [2, 16, 17]. By contrast, using animal models enables researchers to study brain reactions immediately following peripheral nerve injury [18, 19].
To observe the dynamic changes of brain activity during nerve injury and the possible sustained activation, we designed a nerve cutting device for inducing spared nerve injury (SNI) in the functional magnetic resonance imaging (fMRI) chamber. We used a recoverable anesthesia protocol during blood-oxygenation-level-dependent (BOLD)-fMRI [20,21,22,23], and thus could observe neuropathic pain development in the same rat after the BOLD-fMRI experiment.
To correlate the sustained activation with long-term plasticity changes in the brain, we used manganese-enhanced magnetic resonance imaging (MEMRI) to detect plasticity changes in brain activity on Days 1 and 8 after SNI. Recent studies have used MEMRI to assess brain activity in free-moving animals [24,25,26]. Mn2+ can accumulate in excitable cells via voltage-gated Ca2+ channels and ionotropic glutamate receptors [27,28,29], thereby enhancing the signal in T1 contrast images and providing high spatial resolution mapping of brain activations.
Finally, we confirmed the tonic neuronal activations observed in the MR with direct neuronal recording in the insula cortices (ICs) and primary somatosensory cortex (S1). Sustained brain activations induced by SNI were observed by resting fMRI and electrophysiological recording, and these changes correlated with long-term brain plasticity in rats with neuropathic pain.
This study used 55 male and 15 female Sprague Dawley rats (National Laboratory Animal Center, Taipei, Taiwan) aged 8–10 weeks and weighing 250–350 g. Vendor health reports indicated that the rats were free of known viruses, bacteria, and parasites. All the rats were housed pairwise in type 3H cages under a 12-h dark/light cycle with an environmental temperature of 22 °C. The cages were filled with C-grade Sani-Chips and food and water were available ad libitum. All experimental procedures were approved by the Institutional Animal Care and Use Committee of National Taiwan University. This study adhered to the guidelines established by the Council of Agriculture of Taiwan for the experimental use of animals.
Seven rats (four males and three females) were used in the BOLD-fMRI study of dynamic brain activity changes during nerve injury. In the MEMRI study of long-lasting changes of brain activity in free-moving rats, the sham (n = 10) and SNI (n = 10) groups were observed on Day 1 and additional independent sham (n = 10) and SNI (n = 8) groups were observed on Day 8 after SNI surgery. At both time points, the naïve (n = 10) group was observed. Three additional rats were used to validate the adverse behavioral effects of Mn2+ injection. All the rats in the MEMRI study were male. The electrophysiological study observed the following four groups of rats: ipsilateral rostral anterior insular cortex (RAIC), contralateral RAIC, contralateral anterior cingulate cortex (ACC), and contralateral S1. Each group consisted of three rats. Because skull growth in male rats is relatively rapid, which likely affects electrophysiological recording quality, we used female rats in our electrophysiological study. The BOLD-fMRI data (mixed-gender) revealed similar result when comparing with the electrophysiological data (female only) in acute response to the nerve injury. The basic design of this study is diagrammatically shown in the Fig. 1.
SNI neuropathic pain model
The SNI neuropathic pain model used followed the one established by Decosterd and Woolf . Briefly, rats were anesthetized using ketamine hydrochloride (75 mg/kg, i.p.) and xylazine (15 mg/kg i.p.). An incision was made on the left lateral thigh to expose the sciatic nerve at the level of its trifurcation into the sural, tibial, and common peroneal nerves. The left tibial and common peroneal nerves were tightly ligated using 6.0 silk. Two to 4 mm of the distal nerve was removed and the sural nerve was left intact. After nerve transection, we sutured the muscle by using 6.0 silk and closed the incision with 4.0 silk. Subsequently, the rats were recovered from anesthesia and their allodynia behavior was observed 8 days later. The rats received lincomycin (Lita Pharmacy Co. Ltd., Taiwan; 30 mg/kg/day, i.m. for 3 days) after the surgery to prevent infection.
To determine whether the rats had developed neuropathic pain after SNI surgery, we assessed the mechanical allodynia through von Frey filament (North Coast Medical, Inc., Morgan Hill, USA) testing. All behavioral tests were conducted between12:00 and 18:00. Each rat was tested in a 21 × 12 × 14-cm transparent cage after 5 min of habituation. For the mechanical sensitivity test, each rat was placed in a cage with an open wire mesh base and given 15 min for habituation. A set of eight von Frey filaments (0.4, 0.6, 1, 2, 4, 6, 8, and 15 g of bending force) were used to apply increasing amounts of bending force to the lateral plantar surface (the lateral plantar surface was innervated by the sural nerve). If the rat briskly withdrew or one or both of its hind paws flinched, the final lightest von Frey hair force was applied in the subsequent test. Six stimulations were applied in each session and the withdrawal patterns were recorded to determine the 50% withdrawal threshold according to the formula established by Chaplan .
Rat preparation for BOLD-fMRI study
All fMRI experiments were conducted between 12:00 and 18:00. To identify the tonic brain activation after nerve injury, SNI surgery was performed inside the MRI bore during fMRI acquisition. A specially designed nerve transection device consisting of a guidance head and needle with fork tips was implanted before MRI scanning (Fig. 2). Each rat was initially anesthetized using 5% isoflurane and orotracheally intubated for mechanical ventilation. Under isoflurane anesthesia (3.5% mixed with a 30% oxygen and 70% nitrogen mixture), a 1-cm-long incision was made on the lateral surface of the left thigh and the biceps femoris muscle was dissected to expose the sural, common peroneal, and tibial sciatic nerve branches. The nerve transection device was inserted at the caudal surface of the left thigh, through the muscle layers to the sciatic nerve branches. After the guidance head had been removed, the fork tips were exposed. A stainless steel wire was threaded through the needle and the common peroneal and tibial nerves were fixed to the notches of the fork tips by using a slack noose (Fig. 2). Finally, lidocaine hydrochloride jelly USP 2% (Akorn, Inc., Lake Forest, IL, USA) was applied on the wound carefully without contacting sciatic nerve branches, and then the incision was temporally sutured. This setup enabled the cutting of the common peroneal and tibial nerves during fMRI scanning by pulling on the stainless steel wire. In addition, a set of bipolar stainless steel electrodes each was inserted between the second and third digits of both the left forepaw and left hind paw, respectively. We used forepaw stimulation to induce S1 BOLD responses to confirm the suitability of the conditions for fMRI testing and hind paw stimulation to induce S1 BOLD responses in order to confirm that our settings did not damage the sciatic nerve before nerve transection.
Imaging data were acquired using a 7-Tesla scanner with a 30-cm diameter bore (Bruker Biospec 7030 USR, Ettlingen, Germany). The system was equipped with a 670-mT/m (175-μs rise time) actively shielded gradient system (Bruker, BGA12-S) with an inner diameter of 116 mm. A receive-only four-element phased array coil was used to receive radio frequency signals and a linear volume coil was used to transmit radio frequency pulses.
After the rats were secured on the MRI holder, a bolus of dexmedetomidine (0.5 mL; 0.025 mg/kg; Dexdormitor, Orion, Espoo, Finland) was injected into the tail vein. Fifteen minutes after the bolus injection, continuous intravenous infusion of dexmedetomidine (1 mL/h; 0.05 mg/kg/h) and rocuronium bromide (9 mg/kg/h; Sigma–Aldrich, St. Louis, USA) was initiated and the isoflurane concentration was adjusted to 0.5–1% for the entire scanning period . During the MRI scan, rectal temperature was measured using a thermocouple (Model 1025, SA Instruments, Inc., New York, USA) and maintained at 36.5–37.5 °C by using a circulated hot water bed. The end-tidal CO2 level was continuously monitored and adjusted to between 2.5–3.0%, a range previously calibrated for invasive blood gas sampling under identical conditions to ensure normal physiological conditions in the rat.
The anatomical images were acquired using a rapid acquisition with relaxation enhancement (RARE) sequence (10 coronal slices, thickness = 1 mm, repetition time (TR) = 2500 ms, echo time (TE) = 33 ms, matrix size = 160 × 160, field-of-view (FOV) = 25 × 25 mm, average = 2). To improve the magnetic field homogeneity of the acquisition site, local shimming of the brain area was performed before BOLD-fMRI data acquisition (Mapshim; Bruker BioSpin). BOLD-fMRI data were acquired using single-shot gradient-echo echo planar imaging (EPI) (10 coronal slices, thickness = 1 mm, TR = 2000 ms, TE = 22 ms, matrix size = 80 × 80, FOV = 25 × 25 mm, bandwidth = 200 kHz).
For standard forepaw and hind paw stimulation fMRI, an isolated stimulator (S48 Square Pulse Stimulator with Stimulus Isolation Unit, Grass Technologies, West Warwick, USA) was used to deliver constant current pulses to the stainless steel electrodes for forepaw and hind paw stimulation. These constant current pulses consisted of monophasic square wave electrical stimulation (0.5 ms, 2 mA, 9 Hz) divided into five blocks of 20-s on/off cycles. Ten dummy scans and ten additional baseline images were acquired, yielding 120 images in total. After the nerve was confirmed to be healthy, a 10-min fMRI session was conducted (yielding 300 images) with nerve transection at the 151st scan. After fMRI, we sutured the muscle using 6.0 silk and closed the incision with 4.0 silk. Subsequently, the rat was recovered from anesthesia and paralysis by receiving atipamezole hydrochloride (3 mg/kg, i.v.; ANTISEDAN, Orion, Espoo, Finland), for the reversal of the sedative and analgesic effects of dexmedetomidine, and sugammadex sodium (4–8 mg/kg, i.v.; Merck Sharp & Dohme Corp., Kenilworth, NJ, USA), for the reversal of the paralytic effect of rocuronium, and its allodynia behavior was tested 8 days later. The rat was administered lincomycin (Lita Pharmacy Co. Ltd., Taichung, Taiwan; 30 mg/kg/day, i.m. for 3 days) after scanning to prevent infection.
Electrode implantation and recordings
Each rat was anesthetized with sodium pentobarbital (50 mg/kg, i.p.), of which supplemental doses (16 mg/kg, i.p.) were administered when necessary. Craniotomies were performed to expose the brain surface vertical to the recording sites. The coordinates for implantation were as follows: (1) mediodorsal thalamic nucleus, medial part (MDM): right or left: 0.4 mm, posterior: 3 mm, depth: 4.5–5.5 mm; (2) rostral anterior insular cortex (RAIC): right or left: 3–5 mm, anterior: 1–4 mm, depth: 5–6 mm; (3) mediodorsal thalamic nucleus, lateral part (MDL): right: 1.2 mm, posterior: 2.5 mm, depth: 4.5–5.5 mm; (4) anterior cingulate cortex (ACC): right: 1.2 mm, anterior: 1–4 mm, depth: 2.5–3 mm; (5) ventral posterolateral thalamic nucleus (VP): right: 2.8–3.3 mm, posterior: 2.8 mm, depth: 5.5–6.5 mm; (6) S1: right: 1–3 mm, posterior: 1 mm, depth: 600–800 μm. For thalamus recording, we used a bundled microarray electrode consisting of seven tungsten microwires with diameters of 35 μm bare and 50 μm insulated (#100211; California Fine Wire) in a 29-G guide tube . The electrode set for the cortex recording consisted of eight stainless steel microwires arranged in a 2–3-mm-wide array. Each rat was implanted with a matching thalamic and cortical target set (MDM-RAIC, MDL-ACC and VP-SI) of electrodes. Only those rats with good evoked cortical responses under low intensity (below 10 μA) thalamic stimulation were used to ensure the accuracy of the cortical implantations. Four stainless steel screws were set in each rat’s skull to serve as anchors for the electrode sets. To ground the array electrodes, a copper wire was fixed around the anchoring screw positioned at the occipital bone. When all the electrodes were in place, the surface of the skull was covered with dental cement and the wound was sutured. Lincomycin hydrochloride (30 mg/kg, i.m.) was administrated to prevent infection.
After recovering for 1 week, single unit activities were recorded while SNI surgery was conducted. Multiple-channel cortical unit activities were transmitted to a multichannel acquisition processor system (MAP, Plexon, Dallas, USA) through a connecting cable. To record single unit activities while performing SNI surgery, the anesthesia procedure and surgery preparation were identical to those of the fMRI experiment. Single unit recording was initiated 5 min before sciatic nerve transection and continued for 25 min after transection, resulting in 30 min of continuous recording. Spike signals were amplified 7000–32,000-fold, bandpass-filtered at 250 Hz–13 kHz, and digitized at 40 kHz. Well isolated single unit activities in one area were linearly added for a quantitative estimation of the activity change in that cortical area.
Animal preparation for MEMRI study
Two groups of rats underwent MEMRI testing. The first group was tested to determine the cumulative brain activity during the first 24 h of neuropathic pain initiation, whereas the second group was tested to determine the brain activity on Day 8 following neuropathic pain initiation.
Each rat in the first group was initially anesthetized using 5% isoflurane and maintained in a state of anesthesia by using 3.5% isoflurane. SNI surgery was performed as previously described. Following SNI or sham surgery, a solution of 120 mM MnCl2 in saline (75 mg/kg; 2.25 mL/h; MnCl2-4H2O, Sigma–Aldrich, St. Louis, USA) was injected into the tail vein. The MnCl2 injection protocol followed that of a previous study . During surgery, each rat’s body temperature was maintained at 36.5–37.5 °C by using a feedback-controlled heating pad. After surgery, each rat was administered lincomycin (30 mg/kg/day, i.m. for 3 days) to prevent infection. MEMRI scanning was performed 24 h after MnCl2 infusion and each rat’s allodynia behavior was observed on Days 3 and 8 after scanning.
Following our previous study, each rat in the second group underwent initial SNI or sham surgery under ketamine hydrochloride anesthesia (75 mg/kg, i.p.) and xylazine (15 mg/kg i.p.) . Lincomycin (30 mg/kg/day, i.m. for 3 days) was administered after surgery and all rats were tested for allodynia behavior 1 week later. As previously described, rats in which neuropathic pain was observed were administered MnCl2 infusion under isoflurane anesthesia and MEMRI scanning was performed 24 h later.
Imaging data were acquired using the 7-Tesla scanner. The receive-only four-element-phased array coil and linear volume coil were used to receive radio frequency signals and transmit radio frequency pulses, respectively. Each rat was initially anesthetized by 5% isoflurane mixed with oxygen and maintained in a state of anesthesia under 2% isoflurane during scanning. The T1-weighted images were acquired using a RARE sequence (30 coronal slices, thickness = 0.7 mm, TR = 870 ms, TE = 8.4 ms, matrix size = 160 × 160, FOV = 25 × 25 mm, average = 16). During scanning, each rat’s respiration rate was monitored using a pressure sensor placed below the abdomen. A reading within the range of 40–60 breath/min signified stability. Rectal temperature was measured using a thermocouple and maintained at 36.5–37.5 °C by using a circulated hot water bed.
Image processing and analysis
Statistical parametric maps were generated using the SPM version 8 (www.fil.ion.ucl.ac.uk/spm). All raw images were enlarged by a factor of ten to correlate the image dimensions to human data, thereby facilitating the use of SPM, which was developed for use on humans. For the fMRI data, EPI images obtained from a single session were realigned and resliced to their averaged image to minimize movement artifacts. The realigned images were further coregistered to the average image between subjects. Subsequently, these images were smoothed using a Gaussian kernel with a full width at half maximum (FWHM) of 8 mm (in the space of the enlarged images) to reduce white noise and blur the images. Individual statistical maps were plotted using a general linear model with a hemodynamic response function. A 0.0078-Hz high-pass filter was applied to remove slow signal drift. Group analysis of individual statistical maps was conducted using a single sample t test. The significant BOLD response was determined by an individual voxel threshold of P < .02 with a cluster size threshold of 23 continuous voxels (volume = 2.25 mm3). According to the AlphaSim procedure , this threshold combination provides a probability of false-positive clusters within the brain area of P < .05.
In the MEMRI data, the images of each subject were coregistered together and the average image from the coregistered images served as the template image. Subsequently, all MEMRI images were coregistered to the template image again and these fine-aligned images were smoothed using a Gaussian kernel with a FWHM of 3 mm (in the space of the enlarged images). The extra brain area was then removed. Statistical parametric maps were generated and one-way analysis of variance (ANOVA) was performed to compare the groups for differences. Grand mean scaling and global normalization scaling were used and an absolute threshold masking of 0 was applied to exclude the extra brain area from the analysis. Significance was determined by an individual voxel threshold of P < .05 with a cluster size threshold of 60 continuous voxels (volume = 1.025 mm3). According to the AlphaSim procedure , this threshold combination provides a probability of false positive clusters within the brain area of P < .03. The colors of significant areas in statistical maps were encoded from statistical t values. Warm and cold colors indicate increases and decreases in BOLD (or MEMRI) signals, respectively.
Resting-State fMRI Data Analysis Toolkit V1.6 (http://restfmri.net/forum/index.php) was used to conduct functional connectivity analysis of the fMRI data. For fMRI data acquired before and after SNI, EPI images acquired within the first 4 min (first 120 images) and final 4 min (final 120 images) of the 10-min fMRI session were used to generate correlation maps of the pre- and post-SNI groups, respectively. All preprocessed image series were detrended and bandpass-filtered (0.01–0.08 Hz) before the functional connectivity analysis. The correlation maps were transformed into z-statistical parametric maps using REST, and the pre- and post-SNI groups were compared by conducting a paired t test with SPM. The multiple regression analysis function in SPM was used for functional connectivity analysis of the MEMRI images. To determine the brain areas with functional connections to specific regions of interest (ROIs), the inter-subject variability of neural activities in selected ROIs were used as the regressor. These activities were quantified and normalized based on the following equation: average signal intensity within an ROI/average signal intensity of the entire brain.
All data are expressed as the mean ± standard error. A P value of <.05 was considered statistically significant. Brain activity in the sham and SNI groups was compared by conducting an independent two-sample t test. The results of the hind paw withdrawal threshold for both groups were compared through one-way repeated measures ANOVA and Tukey’s post hoc multiple comparison.
Immediate functional brain changes after spared nerve injury surgery
To analyze the acute brain responses to nerve transection and determine the sustained activated brain areas possibly related to neuropathic pain development, we developed a nerve transection device for conducting SNI inside the MRI machine (Fig. 2). Before the formal fMRI experiment of the brain responses to SNI, forepaw and hind paw stimulation–evoked BOLD responses were obtained to confirm the rat’s physiological condition and the healthiness of the sciatic nerve, with the nerve transection devise implanted, respectively. The robust and identical S1 BOLD responses to forepaw stimulation indicated that the rats’ physiological conditions were suitable for the fMRI experiment. The S1 BOLD responses to hind paw stimulation indicated that each rat’s sciatic nerve was functionally intact before SNI (Fig. 3).
Transection of the left tibial and common peroneal nerves produced significant positive BOLD responses in the medial thalamus, hypothalamus, and S1 hind limb area (S1HL) contralateral to the injured nerve, and cingulate cortex (CC) and RAICs bilaterally. Significant negative BOLD responses were also found in the bilateral caudate putamen (CPu) (Fig. 4A). Moreover, examining each ROI, the increases in the BOLD signal in the bilateral ICs and contralateral S1HL area were continuously enhanced throughout the acquisition period after nerve transection (more than 5 min); the CC exhibited a sustained low level of activation; whereas the CPu showed a robust transient negative BOLD response. The auditory cortex (Aud), selected as a negative control, didn’t show significant response during the nerve transection (Fig. 5).
Multi-channel, multiple single unit recording experiments were performed in a separate series of experiment outside the MR chamber to independently examine the neuronal responses in selected forebrain regions to nerve transection. These brain regions included the bilateral RAIC, S1HL and ACC. The locations of the recording sites are shown in Fig. 6. The multiunit recording of the responses to nerve transection in the bilateral RAIC and S1HL area confirmed the sustained enhancement of neural activity (more than 20 min after transection) in both brain regions (Fig. 7).
After fMRI, the rats recovered from anesthesia. The 50% withdrawal threshold of the ipsilateral hind paw was tested on day 3 and day 8 after the fMRI experiment. The thresholds were significantly lower on day 3 and day 8 when comparing with pre-experiment baseline (Fig. 8). This indicated our newly designed nerve cutting device did successfully produced SNI rats.
The significantly enhanced activity immediately after SNI in the RAIC raised the question whether a change in the RAIC functional connectivity network occurred. The EPI images acquired within the first 4 min and the final 4 min of the fMRI experiments were used to generate correlation maps of the pre- and post-SNI groups, respectively. The correlation maps seeded in the contralateral RAIC (the RAIC ROI is shown in Fig. 9a) were transformed into z-statistical parametric maps for comparison between the pre- and post-SNI groups. We found that the functional connections between the contralateral RAIC, contralateral S1, and ipsilateral AIC were enhanced immediately after SNI (Fig. 9b). Upon using the ipsilateral AIC cluster as the ROI to observe the signal time course during the fMRI experiment, we found that the ipsilateral AIC exhibited sustained activation immediately after SNI (Fig. 9c). These data imply that the sustained activities in the ipsilateral AIC observed during the MEMRI experiment may have started at the time point immediately after SNI.
Subacute changes of functional brain activities and connectivity probed with MEMRI 1 day and 8 day after SNI
In addition to using BOLD-fMRI to study the acute brain responses to nerve transection, we used MEMRI to study the 24-h cumulative brain activity on Day 1 after SNI surgery (Fig. 1b). We found reduced activity in the contralateral S1HL area and enhanced activity in the ipsilateral AIC in the SNI surgery group compared to the naïve and the sham surgery groups (Fig. 4b). To determine whether these brain activity changes would be long-term, we conducted MEMRI to study the brain activity on Day 8 after SNI surgery. We found ipsilateral AIC and RAIC activities in the SNI group were higher than their counter parts of the sham group (Fig. 4c) on the Day 8. The MEMRI results showed that only in ipsilateral AIC and RAIC, there were long-term changes. Activities in the S1HL area and CC were not significantly changed on day 1 and day 8 after SNI surgery (Fig. 4b and c); in contrast with the fMRI results (Fig. 4a) and unit recording results (Fig. 7) immediately after SNI surgery. These data showed that SNI surgery induced sustained activation in the IC, and the enhanced activity lasted for at least 8 days.
The observation of a continuous increase in ipsilateral AIC activity after SNI surgery implies that long-term plasticity changes may be triggered within the ipsilateral AIC-related networks. To determine how the ipsilateral AIC interacted with the other brain areas in the SNI neuropathic pain rats, we used the MEMRI images acquired on Day 8 after surgery to extract the ipsilateral AIC activity of various subjects and used these quantified activities as regressors to search for other brain areas with synchronized activity fluctuations. The ROI of the ipsilateral AIC was defined as the area of significant activation in the statistical map in Fig. 4C (SNI > sham and naïve, Day 8) that overlapped within the anatomical AIC region. The green area in Fig. 10A represents the seeding ROI. We found significant cooperation between the ipsilateral AIC versus bilateral RAIC and S1HL area contralateral to the side with nerve injury in the SNI rats (Fig. 10b and c), consistent with our fMRI data (Fig. 9b). This functional correlation was not observed among the sham rats (Fig. 10c); however, although the ipsilateral AIC activity was significantly higher in the SNI rats than the sham rats (P = .0038), the contralateral S1HL activity was similar in both groups (Fig. 10d), indicating that under similar somatosensory input conditions, the ipsilateral AIC exhibited an increase in sensory-related activity only in the SNI rats.
We performed a regression analysis using the signal fluctuations extracted from significant clusters in the ipsilateral RAIC area in Fig. 4c (SNI > sham and naïve, D8). The yellow area in Fig. 11a denotes the seeding ROI. We observed a positive correlation between the ipsilateral RAIC and superficial layer of the ACC and prelimbic areas (Fig. 11b). Notably, the ipsilateral RAIC correlated strongly with the bilateral locus coeruleus (LC) (Fig. 11b) and enhanced activity in the peri-coeruleus regions such as the parabrachial nucleus and nucleus raphe magnus (Fig. 11c). To identify the relative location between the negatively correlated LC and activated peri-coeruleus regions, we overlapped Fig. 11b and c (Fig. 11d). The strong negatively correlated areas are denoted in green (P < .003) and the positively correlated areas are shown in magenta (Fig. 11d). Because the LC contributes to descending pain modulations, these results imply that the descending pain modulation pathway may be involved in the development of neuropathic pain in the SNI model.
This study combined functional brain imaging and electrophysiological recording methods to map sustained brain activation regions during early SNI neuropathic pain development in rats. We found nerve injury induced immediate and sustained activity increase in the bilateral ICs and contralateral S1. The elevated activity in the ipsilateral IC was long-lasting and could be observed on Days 1 and 8 after SNI. Plasticity changes in functional connectivity among the ipsilateral AIC, contralateral S1HL, and bilateral RAIC were observed during the first few days of nerve injury. In addition, the enhanced functional connectivity between the RAIC and LC that spread to the brainstem implied that the descending pain modulation system may be involved in the development of neuropathic pain.
Many functional brain imaging studies of neuropathic pain have been conducted, most of which involved humans with well-established chronic pain. Such studies have focused on the brain responses to various stimulations [36,37,38,39,40] and various brain network patterns [41,42,43]. However, how nerve injury causes neuropathic pain chronification remains unresolved. Although central sensitization at the spinal level caused by peripheral nerve injury has been reported as activation-dependent [44,45,46,47], how nerve injury initially affects brain activity has not been mapped using functional brain imaging methods. Because nerve transection directly activated the peripheral nerve, our findings of brain activation pattern through SNI are similar to those reported in previous studies that use direct peripheral nerve stimulation on rats . The novelty of the current study is the finding that activity in the bilateral ICs, CC, and S1 increased tonically after nerve transection. Such sustained activity may provide a condition for brain sensitization [11, 12].
Some researchers postulate that chronic pain is the result of nociceptive memory in the brain [2, 49, 50]. The center part of this hypothesis is that, similar to memory consolidation, neuropathic pain chronifies inside the brain through a period of sustained high frequency inputs. We tested this hypothesis using resting fMRI, MEMRI and electrophysiological recording experiments. BOLD-fMRI enables the detection of evoked neural activity on the basis of the neurovascular coupling [51,52,53,54], however, it is difficult using BOLD-fMRI to quantify the basal neural activity. Although a resting-state fMRI design can be employed to observe large-scale networks inside the brain , determining whether activity in a specific brain area is high or low over long periods is difficult. By contrast, MEMRI enables the mapping of cumulative brain activity in free-moving animals over relatively long periods of 3–24 h [24,25,26, 33, 56, 57]. Furthermore, a recent study used inter-subject variability of brain activity to analyze functional networks during memory consolidation . In the present study, a similar strategy was employed to investigate neuropathic pain chronification. Rather than using the average signals within the atlas-defined ROIs, we used the average signals within the activation cluster as regressors to assess voxel-wise functional connections. This is because we assumed that only parts of the anatomically defined brain regions exhibited abnormal functional connections; for example, we observed that only part of the contralateral S1HL exhibited plasticity change in functional connectivity with the AIC, possibly because a portion of the right S1HL has been deprived of its input from the tibial and common peroneal nerves.
In many studies, the IC has been consistently activated in various pain paradigms. Some studies have suggested that plasticity in the IC is crucial to the maintenance of neuropathic pain [58,59,60]. However, the IC is involved in various sensory modalities, as well as motor and emotional functions with complex connections with various brain areas . Craig suggested that the IC is a multimodal homeostatic or interoceptive integration area . Its multimodal input property enables it to serve as a multimodal magnitude estimator or salience detector [61, 63]. Some studies have found that the IC encodes pain intensity [64, 65]. However, parcellation of the IC into a pain magnitude estimator or salience detector remains controversial [66, 67]. The IC is currently considered a multidimensional integration site for pain . Recent studies have shown that people with IC lesions can rate the magnitude of evoked acute pain. Because the S1 has been reported to encode pain intensity [69, 70], it may serve as a redundancy system for magnitude estimation . In the present study, we found the ipsilateral AIC exhibited consistently enhanced activity after SNI surgery; and the functional connection between the AIC with the contralateral S1HLalso enhanced under the neuropathic pain condition. These data thereby imply that a plasticity change of the pain magnitude estimation standard may have occurred in the SNI neuropathic pain rats.
Peripheral nerve injury induces S1 reorganization and causes the input-deprived cortex to be occupied by neighborhood expansion . Our data showed that activity within the S1HL decreased on Day 1 after SNI surgery but recovered to a similar level as that of the sham group by Day 8. This recovery may be the result of cortical reorganization and the expanded representation of the neighboring nerves. The association between S1 reorganization and neuropathic pain development is not yet understood. Furthermore, we found the S1HL area began to cooperate with the AIC under the neuropathic pain condition; however, the overall S1HL activity did not differ between the SNI and sham groups. In our previous study, we observed no variations in S1HL responses to von Frey hair stimulation between SNI and naïve rats . Whether and how the S1 is involved in allodynia is the subject of debate. Some studies on people with neuropathic pain have reported increase activity in the S1 with allodynic stimulations [40, 72, 73], whereas others with similar experimental designs have not observed any variations in S1 response [37, 74]. Nevertheless, our data shows the enhanced functional connection between the reorganized S1HL area and AIC in SNI neuropathic pain rats, thereby shedding new light on the association between the S1 and neuropathic pain development.
In this study, the RAIC exhibited sustained activity in response to SNI and strongly enhanced activity on Day 8 after SNI. A previous anatomical study revealed that the RAIC received considerable input from the centrolateral thalamus (CL) and mediodorsal thalamus (MD) , which showed strong activity in response to SNI in the fMRI result. The CL and MD are major sites of termination for the spinothalamic tract [76,77,78,79] and are implicated in the sensorimotor integration of nociceptive processing . Our data showed that after SNI, RAIC activity was strongly and negatively correlated with the LC and several brainstem areas, whereas the peri-coeruleus region exhibited enhanced neural activity. According to a previous anatomical study, the RAIC projects to the peri-coeruleus region, rostroventral medulla, lateral hypothalamus, parabrachial area, dorsal raphe, and periaqueductal gray , all of which have been directly linked to brainstem-mediated descending modulation in spinal nociceptive neurons [81,82,83,84,85,86,87]. However, a previous study demonstrated that the activation of RAIC projection neurons can cause hyperalgesia in a top-down manner . Two groups of RAIC projection neurons respectively act upon two independent subcortical nuclei to modulate the nociceptive threshold. The first group projects to the gamma-aminobutyric acid (GABA) interneurons in the peri-coerulear zone, and these GABAergic interneurons inhibit the LC neurons and affect the noradrenergic bulbospinal projections. The second group projects to the amygdala and causes hyperalgesia when activated . Therefore, our data may imply a plasticity change in the pain threshold through a top-down modulation from the RAIC after SNI.
The limitation of the present study is the relatively small number of rats used in the BOLD fMRI and its supplementary electrophysiological studies. The classical spared nerve injury model transected the 2 larger branch of the sciatic nerve and leaves only the minor sural nerve intact. This is a major trauma to the rat. SNI is a very robust neuropathic pain model with minimal variations. This has been described by Decosterd & Woolf in the abstract of their original paper: “The spared nerve injury model results in early (<24 h), prolonged (>6 months), robust (all animals are responders) behavioral modifications.” . The robustness of the SNI model has been repeated in our lab in several previously published papers [34, 89] and can be clearly observed in the Fig. 8 of the present study. Nevertheless, the small number of animal used prevents a meaningful analysis of possible difference between male and female rats.
In summary, this study constitutes the first experimental observation of the forebrain regions with tonically enhanced activity induced by peripheral nerve injury and the first study to correlate this tonic brain activity with long-term brain plasticity. The SNI-induced plasticity changes in the ipsilateral AIC, bilateral RAICs, and S1 may contributed to the neuropathic pain development through changes in the sensorimotor integration of nociceptive information. Moreover, the negative functional correlation between the RAIC and LC suggested a change in the descending pain modulation system under the condition of neuropathic pain development.
anterior cingulate cortex
anterior insular cortex
blood-oxygenation-level-dependent functional magnetic resonance imaging
echo planar imaging
full width at half maximum
mediodorsal thalamic nucleus, lateral part
mediodorsal thalamic nucleus, medial part
manganese-enhanced magnetic resonance imaging
rostral anterior insular cortex
regions of interest
primary somatosensory cortex
S1 hind limb
spared nerve injury
ventral posterolateral thalamic nucleus
Goldberg DS, McGee SJ. Pain as a global public health priority. BMC Public Health 2011;11(1):1.
Apkarian AV, Baliki MN, Geha PY. Towards a theory of chronic pain. Prog Neurobiol. 2009;87(2):81–97.
Li X-Y, Ko H-G, Chen T, Descalzi G, Koga K, Wang H, et al. Alleviating neuropathic pain hypersensitivity by inhibiting PKMzeta in the anterior cingulate cortex. Science. 2010;330:1400–4.
Tan W, Yao W-L, Hu R, Lv Y-Y, Wan L, Zhang C-H, et al. Alleviating neuropathic pain mechanical allodynia by increasing Cdh1 in the anterior cingulate cortex. Mol Pain. 2015;11:56.
Han J, Kwon M, Cha M, Tanioka M, Hong S-K, Bai SJ, et al. Plasticity-related PKMζ signaling in the insular cortex is involved in the modulation of neuropathic pain after nerve injury. Neural Plast. 2015:2015.
Li X-Y, Ko H-G, Chen T, Descalzi G, Koga K, Wang H, et al. Alleviating neuropathic pain hypersensitivity by inhibiting PKMζ in the anterior cingulate cortex. Science. 2010;330(6009):1400–4.
Woolf CJ. Evidence for a central component of post-injury pain hypersensitivity. Nature. 1983.
Woolf CJ, Wall PD. Relative effectiveness of C primary afferent fibers of different origins in evoking a prolonged facilitation of the flexor reflex in the rat. J Neurosci. 1986;6(5):1433–42.
Cook AJ, Woolf CJ, Wall PD, McMahon SB. Dynamic receptive field plasticity in rat spinal cord dorsal horn following C-primary afferent input. Nature. 1987;325(6100):151–3.
Ji R-R, Kohno T, Moore KA, Woolf CJ. Central sensitization and LTP: do pain and memory share similar mechanisms? Trends Neurosci. 2003;26(12):696–705.
Collingridge GL, Isaac JT, Wang YT. Receptor trafficking and synaptic plasticity. Nat Rev Neurosci. 2004;5(12):952–62.
Morimoto K, Fahnestock M, Racine RJ. Kindling and status epilepticus models of epilepsy: rewiring the brain. Prog Neurobiol. 2004;73(1):1–60.
Ali Z, Ringkamp M, Hartke TV, Chien HF. Flavahan Na, Campbell JN, et al. uninjured C-fiber nociceptors develop spontaneous activity and alpha-adrenergic sensitivity following L6 spinal nerve ligation in monkey. J Neurophysiol. 1999;81:81.
Liu C-N, Wall PD, Ben-Dor E, Michaelis M, Amir R, Devor M. Tactile allodynia in the absence of C-fiber activation: altered firing properties of DRG neurons following spinal nerve injury. Pain. 2000;85(3):503–21.
Tseng W-T, Tsai M-L, Iwata K, Yen C-T. Long-term changes in trigeminal ganglionic and thalamic neuronal activities following inferior alveolar nerve transection in behaving rats. J Neurosci. 2012;32:16051–63.
Baliki MN, Petre B, Torbey S, Herrmann KM, Huang L, Schnitzer TJ, et al. Corticostriatal functional connectivity predicts transition to chronic back pain. Nat Neurosci. 2012;15:1117–9.
Burgmer M, Pfleiderer B, Maihöfner C, Gaubitz M, Wessolleck E, Heuft G, et al. Cerebral mechanisms of experimental hyperalgesia in fibromyalgia. Eur J Pain. 2012;16(5):636–47.
Kim S, Nabekura J. Rapid synaptic remodeling in the adult somatosensory cortex following peripheral nerve injury and its association with neuropathic pain. J Neurosci. 2011;31:5477–82.
Li R, Hettinger PC, Machol JA, Liu X, Stephenson J, Pawela CP, et al. Cortical plasticity induced by different degrees of peripheral nerve injuries: a rat functional magnetic resonance imaging study under 9.4 tesla. Journal of brachial plexus and peripheral nerve injury. 2013;8(1):4.
Chao THH, Chen JH, Yen CT. Repeated BOLD-fMRI imaging of deep brain stimulation responses in rats. PLoS One. 2014;9.
Weber R, Ramos-Cabrer P, Wiedermann D, Van Camp N, Hoehn M. A fully noninvasive and robust experimental protocol for longitudinal fMRI studies in the rat. NeuroImage 2006;29:1303–1310.
Niskanen J-P, Airaksinen AM, Sierra A, Huttunen JK, Nissinen J, Karjalainen PA, et al. Monitoring Functional Impairment and Recovery After Traumatic Brain Injury in Rats by fMRI. J Neurotrauma. 2012(ja).
Weber R, Ramos-Cabrer P, Justicia C, Wiedermann D, Strecker C, Sprenger C, et al. Early prediction of functional recovery after experimental stroke: functional magnetic resonance imaging, electrophysiology, and behavioral testing in rats. J Neurosci. 2008;28(5):1022–9.
Eschenko O, Canals S, Simanova I, Beyerlein M, Murayama Y, Logothetis NK. Mapping of functional brain activity in freely behaving rats during voluntary running using manganese-enhanced MRI: implication for longitudinal studies. NeuroImage. 2010;49:2544–55.
Hattori S, Hagihara H, Ohira K, Aoki I, Saga T, Suhara T, et al. In vivo evaluation of cellular activity in αCaMKII heterozygous knockout mice using manganese-enhanced magnetic resonance imaging (MEMRI). Front Integr Neurosci. 2013;7:76.
Chen KH, Chen DY, Liang KC. Functional connectivity changes during consolidation of inhibitory avoidance memory in rats: a manganese-enhanced MRI study. The Chinese journal of physiology. 2013;56:269–81.
Itoh K, Sakata M, Watanabe M, Aikawa Y, Fujii H. The entry of manganese ions into the brain is accelerated by the activation of N-methyl-d-aspartate receptors. Neuroscience. 2008;154:732–40.
Hankir MK, Parkinson JR, Bloom SR, Bell JD. The effects of glutamate receptor agonists and antagonists on mouse hypothalamic and hippocampal neuronal activity shown through manganese enhanced MRI. NeuroImage. 2012;59:968–78.
Silva AC, Bock NA. Manganese-enhanced MRI: An exceptional tool in translational neuroimaging. Schizophr Bull 2008;34:595–604.
Decosterd I, Woolf CJ. Spared nerve injury: an animal model of persistent peripheral neuropathic pain. Pain. 2000;87(2):149–58.
Chaplan SR, Bach FW, Pogrel JW, Chung JM, Yaksh TL. Quantitative assessment of tactile allodynia in the rat paw. J Neurosci Methods. 1994;53:55–63.
Tseng W-T, Yen C-T, Tsai M-L. A bundled microwire array for long-term chronic single-unit recording in deep brain regions of behaving rats. J Neurosci Methods. 2011;201(2):368–76.
Chuang KH, Koretsky AP, Sotak CH. Temporal changes in the T1 and T2 relaxation rates (??R1 and ??R2) in the rat brain are consistent with the tissue-clearance rates of elemental manganese. Magn Reson Med. 2009;61:1528–32.
Lin H-C, Huang Y-H, Chao T-HH, Lin W-Y, Sun W-Z, Yen C-T. Gabapentin reverses central hypersensitivity and suppresses medial prefrontal cortical glucose metabolism in rats with neuropathic pain. Mol Pain. 2014;10:63.
Ward BD. Simultaneous Inference for fMRI Data. 2000:1–16.
Apkarian AV, Hashmi JA, Baliki MN. Pain and the brain: specificity and plasticity of the brain in clinical chronic pain. Pain. 2011;152:S49–64.
Witting N, Kupers RC, Svensson P, Jensen TS. A PET activation study of brush-evoked allodynia in patients with nerve injury pain. Pain. 2006;120:145–54.
Schweinhardt P, Glynn C, Brooks J, McQuay H, Jack T, Chessell I, et al. An fMRI study of cerebral processing of brush-evoked allodynia in neuropathic pain patients. NeuroImage. 2006;32:256–65.
Samuelsson M, Leffler AS, Hansson P. Dynamic mechanical allodynia: on the relationship between temporo-spatial stimulus parameters and evoked pain in patients with peripheral neuropathy. Pain. 2005;115:264–72.
Peyron R, Schneider F, Faillenot I, Convers P, Barral F-G, Garcia-Larrea L, et al. An fMRI study of cortical representation of mechanical allodynia in patients with neuropathic pain. Neurology. 2004;63:1838–46.
Baliki MN, Geha PY, Apkarian AV, Chialvo DR. Beyond feeling: chronic pain hurts the brain, disrupting the default-mode network dynamics. J Neurosci. 2008;28:1398–403.
Welsh RC, Chen AC, Taylor SF. Low-frequency BOLD fluctuations demonstrate altered thalamocortical connectivity in schizophrenia. Schizophr Bull. 2010;36:713–22.
Cardoso-Cruz H, Lima D, Galhardo V. Impaired spatial memory performance in a rat model of neuropathic pain is associated with reduced hippocampus-prefrontal cortex connectivity. J Neurosci. 2013;33:2465–80.
Sandkühler J, Liu X. Induction of long-term potentiation at spinal synapses by noxious stimulation or nerve injury. Eur J Neurosci. 1998;10:2476–80.
Gruber-Schoffnegger D, Drdla-Schutting R, Hönigsperger C, Wunderbaldinger G, Gassner M, Sandkühler J. Induction of thermal hyperalgesia and synaptic long-term potentiation in the spinal cord lamina I by TNF-α and IL-1β is mediated by glial cells. J Neurosci. 2013;33:6540–51.
Zhang H-M, Zhou L-J, Hu X-D, Hu N-W, Zhang T, Liu X-G. Acute nerve injury induces long-term potentiation of C-fiber evoked field potentials in spinal dorsal horn of intact rat. Sheng li xue bao : [Acta physiologica Sinica]. 2004;56:591–6.
Liu X-G, Zhou L-J. Long-term potentiation at spinal C-fiber synapses: a target for pathological pain. Curr Pharm Des. 2015;21:895–905.
Cho YR, Jones SR, Pawela CP, Li R, Kao DS, Schulte ML, et al. Cortical brain mapping of peripheral nerves using functional magnetic resonance imaging in a rodent model. J Reconstr Microsurg. 2008;24:551–7.
Yi M, Zhang H. Nociceptive memory in the brain: cortical mechanisms of chronic pain. J Neurosci. 2011;31:13343–5.
Mansour AR, Farmer MA, Baliki MN, Apkarian AV. Chronic pain: the role of learning and brain plasticity. Restor Neurol Neurosci. 2014;32:129–39.
Pasley B, Freeman R. Neurovascular coupling. Scholarpedia. 2008;3:5340.
Harder DR, Alkayed NJ. Lange aR, Gebremedhin D, Roman RJ. Functional hyperemia in the brain: hypothesis for astrocyte-derived vasodilator metabolites. Stroke. 1998;29:229–34.
Attwell D, Iadecola C. The neural basis of functional brain imaging signals. Trends Neurosci. 2002;25(12):621–5.
Logothetis NK. The neural basis of the blood-oxygen-level-dependent functional magnetic resonance imaging signal. Philos Trans R Soc Lond Ser B Biol Sci. 2002;357:1003–37.
Biswal B, Yetkin FZ, Haughton VM, Hyde JS. Functional connectivity in the motor cortex of resting human brain using echo-planar MRI. Magnetic resonance in medicine : official journal of the Society of Magnetic Resonance in Medicine / Society of Magnetic Resonance in Medicine. 1995;34:537–41.
Silva AC, Lee JH, Aoki I, Koretsky AP. Manganese-enhanced magnetic resonance imaging (MEMRI): methodological and practical considerations. NMR Biomed. 2004;17:532–43.
Jeong K-Y, Lee C, Cho J-H, Kang J-H, Na H-S. New method of manganese-enhanced magnetic resonance imaging (MEMRI) for rat brain research. Experimental animals / Japanese Association for Laboratory Animal Science. 2012;61:157–64.
Ferrier J, Bayet-Robert M, Dalmann R. El Guerrab a, Aissouni Y, Graveron-Demilly D, et al. cholinergic neurotransmission in the posterior insular cortex is altered in preclinical models of neuropathic pain: key role of muscarinic M2 receptors in donepezil-induced Antinociception. J Neurosci. 2015;35:16418–30.
Qiu S, Zhang M, Liu Y, Guo Y, Zhao H, Song Q, et al. GluA1 phosphorylation contributes to postsynaptic amplification of neuropathic pain in the insular cortex. J Neurosci. 2014;34:13505–15.
Han J, Kwon M, Cha M, Tanioka M. Hong S-k, Bai SJ, et al. Plasticity-Related PKM ζ Signaling in the Insular Cortex Is Involved in the Modulation of Neuropathic Pain after Nerve Injury Neural Plasticity. 2015;2015:1–10.
Circuitry and functional aspects of the insular lobe in primates including humans, (1996).
Craig AD. Interoception: the sense of the physiological condition of the body. Curr Opin Neurobiol. 2003;13:500–5.
Downar J, Crawley AP, Mikulis DJ, Davis KD. A cortical network for the detection of novel events across multiple sensory modalities. NeuroImage. 2001;13:S310-S.
Coghill RC, Sang CN, Maisog JM, Iadarola MJ. Pain intensity processing within the human brain: a bilateral, distributed mechanism. J Neurophysiol. 1999;82:1934–43.
Baliki MN, Geha PY. Apkarian aV. parsing pain perception between nociceptive representation and magnitude estimation. J Neurophysiol. 2009;101:875–87.
Moayedi M, Weissman-Fogel I. Is the insula the "how much" intensity coder? J Neurophysiol. 2009;102:1345–7.
Downar J, Mikulis DJ, Davis KD. Neural correlates of the prolonged salience of painful stimulation. NeuroImage. 2003;20:1540–51.
Brooks JCW, Tracey I. The insula: A multidimensional integration site for pain. Pain2007. p. 1–2.
Ea M, Keaser ML, Gullapalli RP, Greenspan JD. Regional intensive and temporal patterns of functional MRI activation distinguishing noxious and innocuous contact heat. J Neurophysiol. 2005;93:2183–93.
Rainville P, Carrier B, Hofbauer RK, Bushnell MC, Duncan GH. Dissociation of sensory and affective dimensions of pain using hypnotic modulation. Pain. 1999;82:159–71.
Kaas JH, Merzenich MM, Killackey HP. The reorganization of somatosensory cortex following peripheral nerve damage in adult and developing mammals. Annu Rev Neurosci. 1983;6:325–56.
Petrovic P, Ingvar M, Stone-Elander S, Petersson KM, Hansson P. A PET activation study of dynamic mechanical allodynia in patients with mononeuropathy. Pain. 1999;83:459–70.
Peyron R, Faillenot I, Pomares FB, Le Bars D, Garcia-Larrea L, Laurent B. Mechanical allodynia in neuropathic pain. Where are the brain representations located? A positron emission tomography (PET) study. Eur J Pain 2013;17:1327–1337.
Kim CE, Kim YK, Chung G, Im HJ, Lee DS, Kim J, et al. Identifying neuropathic pain using 18F-FDG micro-PET: a multivariate pattern analysis. NeuroImage. 2014;86:311–6.
Jasmin L, Granato A, Ohara PT. Rostral agranular insular cortex and pain areas of the central nervous system: a tract-tracing study in the rat. J Comp Neurol. 2004;468(3):425–40.
Cliffer K, Burstein R, Giesler G. Distributions of spinothalamic, spinohypothalamic, and spinotelencephalic fibers revealed by anterograde transport of PHA-L in rats. J Neurosci. 1991;11(3):852–68.
Ma W, Peschanski M, Ralston HJ. Fine structure of the spinothalamic projections to the central lateral nucleus of the rat thalamus. Brain Res. 1987;414(1):187–91.
Ganchrow D. Intratrigeminal and thalamic projections of nucleus caudalis in the squirrel monkey (Saimiri sciureus): a degeneration and autoradiographic study. J Comp Neurol. 1978;178(2):281–311.
Craig A. Distribution of trigeminothalamic and spinothalamic lamina I terminations in the macaque monkey. J Comp Neurol. 2004;477(2):119–48.
Itoh K, Mizuno N. Topographical arrangement of thalamocortical neurons in the centrolateral nucleus (CL) of the cat, with special reference to a spino-thalamo-motor cortical path through the CL. Exp Brain Res. 1977;30(4):471–80.
Aimone L, Bauer C, Gebhart G. Brain-stem relays mediating stimulation-produced antinociception from the lateral hypothalamus in the rat. J Neurosci. 1988;8(7):2652–63.
Basbaum AI, Fields HL. Endogenous pain control systems: brainstem spinal pathways and endorphin circuitry. Annu Rev Neurosci. 1984;7(1):309–38.
Carstens E, Fraunhoffer M, Suberg S. Inhibition of spinal dorsal horn neuronal responses to noxious skin heating by lateral hypothalamic stimulation in the cat. J Neurophysiol. 1983;50(1):192–204.
Chiang C, Hu JW, Sessle BJ. Parabrachial area and nucleus raphe magnus-induced modulation of nociceptive and nonnociceptive trigeminal subnucleus caudalis neurons activated by cutaneous or deep inputs. J Neurophysiol. 1994;71(6):2430–45.
Chiang CY, Sessle BJ, Hu JW. Parabrachial area and nucleus raphe magnus-induced modulation of electrically evoked trigeminal subnucleus caudalis neuronal responses to cutaneous or deep A-fiber and C-fiber inputs in rats. Pain. 1995;62(1):61–8.
Millan MJ. Descending control of pain. Prog Neurobiol. 2002;66(6):355–474.
Prado WA, Faganello FA. The anterior pretectal nucleus participates as a relay station in the glutamate-, but not morphine-induced antinociception from the dorsal raphe nucleus in rats. Pain. 2000;88(2):169–76.
Jasmin L, Rabkin SD, Granato A, Boudah A, Ohara PT. Analgesia and hyperalgesia from GABA-mediated modulation of the cerebral cortex. Nature. 2003;424(6946):316–20.
Lee K-S, Huang Y-H, Yen C-T. Periaqueductal gray stimulation suppresses spontaneous pain behavior in rats. Neurosci Lett. 2012;514(1):42–5.
The authors thank 7 T animal MRI Core Lab of the Neurobiology and Cognitive Science Center, National Taiwan University for technical and facility supports.
Availability of data and material
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
This study was supported by a grant from Ministry of Science and Technology, Taiwan, MOST105–2311-B-002 -029.
Ethics approval and consent to participate
All experimental procedures were approved by the Institutional Animal Care and Use Committee of National Taiwan University. This study adhered to the guidelines established by the Council of Agriculture of Taiwan for the experimental use of animals.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.