Email updates

Keep up to date with the latest news and content from Molecular Pain and BioMed Central.

Open Access Highly Accessed Research

The analgesic effect of electroacupuncture on acute thermal pain perception-a central neural correlate study with fMRI

Shivshil Shukla1, Artour Torossian2, Jeng-Ren Duann3 and Albert Leung4*

Author Affiliations

1 Anesthesia Service, Veteran Administrations San Diego Healthcare System; 3350 La Jolla Village Drive, MC 125, San Diego, CA92161, USA

2 Vanderbilt University, School of Medicine; D-3300 Medical Center North, Nashville, TN 37232-2104, USA

3 Swartz Center for Computational Neuroscience, Institute for Neural Computation, The University of California, San Diego; 9500 Gilman Drive # 0559, La Jolla, CA 92093-0559, USA

4 Department of Anesthesiology, The University of California, San Diego, School of Medicine; Veteran Administrations San Diego Healthcare System; 9300 Campus Point Drive, A-113, MC 7651, La Jolla, CA 92037-1300, USA

For all author emails, please log on.

Molecular Pain 2011, 7:45  doi:10.1186/1744-8069-7-45


The electronic version of this article is the complete one and can be found online at: http://www.molecularpain.com/content/7/1/45


Received:8 December 2010
Accepted:7 June 2011
Published:7 June 2011

© 2011 Shukla et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Electrical acupuncture (EA) has been utilized in acute pain management. However, the neuronal mechanisms that lead to the analgesic effect are still not well defined. The current study assessed the intensity [optimal EA (OI-EA) vs. minimal EA (MI-EA)] effect of non-noxious EA on supraspinal regions related to noxious heat pain (HP) stimulation utilizing an EA treatment protocol for acute pain and functional magnetic resonance imaging (fMRI) with correlation in behavioral changes. Subjects underwent five fMRI scanning paradigms: one with heat pain (HP), two with OI-EA and MI-EA, and two with OI-EA and HP, and MI-EA and HP.

Results

While HP resulted in activations (excitatory effect) in supraspinal areas known for pain processing and perception, EA paradigms primarily resulted in deactivations (suppressive effect) in most of these corresponding areas. In addition, OI-EA resulted in a more robust supraspinal sedative effect in comparison to MI-EA. As a result, OI-EA is more effective than MI-EA in suppressing the excitatory effect of HP in supraspinal areas related to both pain processing and perception.

Conclusion

Intensities of EA plays an important role in modulating central pain perception.

Background

Recent studies which explored the potential use of acupuncture in preemptive and acute pain analgesia indicated that acupuncture could reduce intraoperative opioid requirement and improve postoperative pain control [1-4]. However, the neuronal events leading to the analgesic effect of acupuncture on acute pain perception are largely unknown. In the context of acute pain treatment, a particular acupuncture treatment system called the tendinomuscular meridian has been described in acupuncture textbooks as an effective means to treat acute pain located in extremities [5-7]. A previous study adopting this acute pain acupuncture treatment paradigm and a thermal human experimental pain model with peripheral sensory testing have demonstrated that a short duration of unilateral electrical acupuncture (EA) at one lower extremity could transiently reduce heat pain (HP) perception with corresponding thermal threshold elevations in bilateral lower extremities, suggesting this acupuncture treatment protocol may have a direct effect on central neuromodulation of pain [8]. A follow-up study using the same treatment protocol with an extended duration of stimulation resulted in a prolonged period of analgesic benefit [9]. However, aside from the stimulation duration effect, how the intensity of the stimulation may affect the analgesic and corresponding central neuromodulatory effect is still unknown. Clinically, the adequate dosing of acupuncture is commonly assessed in terms of the intensity of deqi sensation. Therefore, it is mechanistically important to understand how acupuncture may affect cortical neuromodulation in the acute thermal pain state based on the intensity of the stimulation. Here we hypothesize that the neuromodulatory effect of acupuncture on acute pain is stimulation intensity dependent and different intensities of stimulation will result in different degrees/patterns of cortical activation which subsequently result in different degrees of alteration in pain perception when stimulated with the same intensity of noxious thermal simulation. To test our hypothesis, we conducted the current study utilizing a published acute pain treatment acupuncture paradigm, thermal noxious stimulation and functional magnetic resonance imaging (fMRI), to: 1) assess the baseline supraspinal blood oxygen level dependent (BOLD) response related to HP stimulation and two (minimal and optimal) intensities of EA with behavioral correlation in pain and deqi (tingling sensation) VAS scores respectively; and 2) assess the direct effect of the different intensities of EA on supraspinal response to subsequent HP stimulation. The authors are aware of different sensations associated with deqi sensation and chose "tingling sensation" as the primary acupuncture intensity assessment based on a previous study which established this particular sensation as the most prominent sensation associated with EA [10].

Results

10 right-handed normal subjects (6 females and 4 males) were enrolled for the study.

Thermal threshold and HP VAS scores

The average threshold (°C ± SD) for cold, warm, cold pain and heat pain for the subjects (n = 10) were 27.8 ± 1.8, 35.9 ± 2.1, 11.2 ± 10.2 and 46.5 ± 1.1 respectively. The average HP VAS scores (± SD) in the HP only, HP with minimal intensity EA (MI-EA) and HP with optimal intensity EA (OI-EA) paradigms were 41.0 ± 8.0, 37.5 ± 8.9 and 13.8 ± 4.6 respectively. Paired samples T-test showed significant (P < 0.05) differences in HP VAS scores between HP only and HP with OI-EA paradigms, and also between HP with MI-EA and HP with OI-EA paradigms. However no statistically significant difference was observed between HP only and HP after MI-EA VAS scores.

EA intensities and Deqi sensation

For the stimulation intensity assessment, the average deqi VAS scores (± SD) for MI-EA and OI-EA were 13.1 ± 8.1 and 59.9 ± 8.8 respectively. For the EA scanning paradigms, the deqi VAS scores were 12.5 ± 7.8 and 61 ± 8.6 for the MI-EA and OI-EA respectively. The average numerical intensities for MI-EA and OI-EA were 3.0 ± 0.5 and 8.0 ± 0.5 respectively on a 0-10 dialing scale. The difference between the MI-EA and OI-EA deqi VAS scores was statistically significant (P < 0.001).

fMRI findings

An inflated cortical representation of identified brain areas of activation (positively correlated BOLD signal) and deactivation (negatively correlated BOLD signals) in all five paradigms are shown in Figure 1. The result of within-group and between-group analyses is summarized as follows:

thumbnailFigure 1. The inflated group (n = 10) cortical representation of brain activations and deactivations (P < 0.01) in all five paradigms. HP: heat pain; MI-EA: Minimal intensity electrical acupuncture; OI-EA: Optimal intensity electrical acupuncture; A: Anterior; P: Posterior; LH: Left Hemisphere.

HP only

The within-group random effects analysis of the HP only paradigms (Table 1) showed significant (P < 0.01) activations at the right thalamus, right secondary somatosensory cortex (SII) and right parietal lobe (BA7,39), right insular cortex (BA13), right premotor cortex(BA6&44), bilateral frontal association cortex (BA9), right dorsal anterior cingulate cortex(BA 32) and left parietal lobe (BA39), and deactivations at the left medial frontal cortex(BA11), left temporal limbic cortex (BA38), left lateral temporal lobe(BA21) and left occipital lobe (BA19).

Table 1. Brain regions of activities related to Heat Pain (Rt: Right; Lt: Left; SII: Secondary Somatosensory cortex)

Electrical acupuncture (OI-EA and MI-EA)

The within-group random effects analysis of the MI-EA paradigm (Table 2) showed significant (P < 0.01) activations at the right caudate nucleus, right dorsolateral prefrontal cortex (BA46) and right inferior frontal cortex (BA44) and the left frontal association cortex (BA45), and significant (P < 0.01) deactivation at the right SI (BA1), bilateral SII (BA5 &7), bilateral frontal cortices (BA10), and bilateral dorsal posterior cingulate cortex (BA31). On the other hand, the within-group random effects analysis of the OI-EA paradigm showed only significant (P < 0.01) deactivations at the right SII (BA7), left dorsal anterior cingulate cortex (BA32), left premotor cortex(BA6&8), left Raphe nucleus, left parietal lobe (BA39) and bilateral amygdale.

Table 2. Brain Region of activities related to Minimal(MI-EA) and Optimal(OI-EA) intensities of EA (Rt: right; Lt: Left; SI: Primary somatosensory cortex; SII: Secondary Somatosensory cortex)

Subsequently, the between-group random effects analysis (OI-EA > MI-EA) showed OI-EA achieved a significant (p < 0.01) degree of deactivation than MI-EA (Figure 2) in the left amygdala, right Raphe nucleus within the reticular formation, bilateral frontal lobes (BA6, 8, 10 11) and right inferior frontal lobe (BA 44).

thumbnailFigure 2. Brain regions with significant (P < 0.01) deactivation with OI-EA > MI-EA group comparison.

Heat pain with electrical acupuncture

In the HP with MI-EA within-group random effects analysis (Table 3), significant activations were found (P < 0.01) at the right thalamus, right caudate nucleus, right frontal lobe(BA47), right prefrontal association cortex(BA10), right limbic association cortex (BA29), bilateral visual cortices (19) and left dorsal lateral prefrontal cortex. No significant areas of deactivation were extracted from the analysis. In the HP with OI-EA within-group random effects analysis, significant (P < 0.01) activations were extracted at the right premotor cortex (BA44), right parietal association cortex (BA39) and right visual association cortex (BA37). For areas of deactivation, the analysis showed a significant effect (P < 0.01) at the left primary somatosensory cortex (SI/BA3), left SII (BA7), and left premotor cortex (BA6 &8). To assess both the overall treatment and the intensity effect of EA on HP, a second level 2-factor (EA &intensity) ANOVA analysis was then performed to assess the specific areas of activation identified with HP stimulation for the two HP with EA paradigms. Except for bilateral frontal association cortex (BA 39), all the other areas including SII, insular, right premotor cortex, anterior cingulate cortex showed a significant (< 0.01) EA × intensity effect. The average t-values of these heat pain related brain regions after the two different EA paradigms were shown in Figure 3. In addition, to assess the effect of EA on the spinothalamic pain pathway, a side-by-side comparison(Figure 4) was also conducted to assess the effect of MI-EA and OI-EA on the thalamic activation induced by HP and showed that OI-EA completely suppressed the HP-induced thalamic activation, whereas, MI-EA moderately reduced the thalamic activation.

Table 3. Brain regions of activities related to Heat Pain given under Minimal and Optimal intensities of EA(MI-EA: minimal intensity electrical acupuncture; OI-EA: optimal intensity electrical acupuncture; Rt: Right; Lt: Left)

thumbnailFigure 3. Comparison of the average T-value in the heat pain (HP) related regions (average voxel x, y, z coordinates) after the two different modes of EA: minimal intensity electrical acupuncture (MI-EA) and optimal intensity electrical acupuncture (OI-EA). * indicates regions with significant difference in EA × Intensity interaction between the two paradigms.

thumbnailFigure 4. The relative degree of right thalamic activation observed with heat pain (HP), HP with minimal intensity electrical acupuncture (MI-EA) and HP with optimal intensity electrical acupuncture (OI-EA).

Discussion

Recent studies with human experimental pain and functional imaging techniques have provided insightful information regarding areas of the central nervous systems that may play a role in encoding noxious stimulation [11-15]. These pain-related regions-of-activities (ROA) include SI and SII, thalamus, insular cortex, cerebellum, amygdale, prefrontal cortex, premotor and motor cortex. It is also known that collateral pain ascending pathways exist and individual afferent neurons often project in more than one of these pathways [16-18]. These pathways include a spinohypothalamic pathway, a spinopontoamygdaloid pathway, and a component of the spinothalamic pathway that projects to specific midline thalamic nuclei which further projects to limbic cortical areas such as the anterior cingulate cortex (ACC) and insular cortex (IC). Another component of the spinothalamic pathway projects to somatosensory relay nuclei of the thalamus that further relays nociceptive information to SI and SII, which are anatomically interconnected with a ventrally directed cortico-limbic somatosensory pathway that integrates somatosensory input with other sensory modalities such as vision and audition, and learning and memory. In addition, from SI and SII the pathway projects into the posterior parietal cortical areas and the insular cortex, and from the insular cortex to amygdala, perirhinal cortex (medial temporal lobe, BA35 &36), and hippocampus [19]. The latter supraspinal components also have a dual convergence from the ascending spinopontoamygdaloid pathway. These pain related ROA can be further functionally divided into the lateral and the medial systems [20-24]. The lateral system, which is thought to be responsible for the initial noxious signal encoding, carries a somatosensory/discriminatory sensory function, and consists of SI, SII, premotor cortex, motor cortex, and cerebellum, whereas, the medial system, composed of the amygdala, cingulate gyrus and insular cortex, is thought to underlie the affective component of the pain experience. In addition, other cortical areas located in the prefrontal cortical region are thought to play a role in pain modulation [25-28]. In the current study, the areas of supraspinal activities detected from HP stimuli delivered via the block imaging design is consistent with the previous studies.

In the area of acupuncture research, functional imaging studies with various settings of non-noxious acupuncture demonstrated activities in the somatosensory components of pain and deactivations in the affective areas of pain processing such as amygdala, hippocampus and cingulate cortex [29-32]. Although these previous functional imaging studies have distinguished the difference in supraspinal activities between non-painful and painful deqi sensations, how the intensity of non-noxious EA may affect subsequent supraspinal activities in HP modulation is unknown. To the authors' best knowledge, this is the first fMRI study to assess the preemptive effect of acupuncture on central acute thermal pain processing. In the current study, HP alone predominantly induced supraspinal response in the somatosensory (thalamus, contralateral SII), emotional (ACC, insular) and modulatory (bilateral prefrontal cortices, contralateral premotor cortex) components of pain perception. These findings are consistent with results from previous functional imaging studies as discussed above.

In the two EA paradigms with different intensities, significant differences in the supraspinal response between the OI-EA and MI-EA paradigms were noted. First, OI-EA appeared to induce a predominantly deactivation response at the ipsilateral cortical regions that were related to HP processing, whereas, MI-EA induced a mixture of activations and deactivations in both ipsilateral and contralateral hemispheres. Second, while both EA paradigms resulted in a similar degree of deactivation in the SII and ACC, the result of the between-group comparison (OI-EA > MI-EA) demonstrated that the OI-EA generated a much more significant level of deactivation in the contralateral Raphe nucleus, inferior frontal cortex and amygdala, suggesting a more effective role of OI-EA in suppressing the limbic (medial pain pathway) system in comparison to the MI-EA. Moreover, in the HP with EA paradigms, HP with OI-EA resulted in deactivations of the ipsilateral somatosensory cortex(SI, SII) and insular cortex, and no activation in the medial limbic system, whereas, HP with MI-EA resulted in activations in the right limbic system, left dorsal lateral prefrontal cortex, right caudate nucleus, right prefrontal and frontal cortex. In addition to the observed medial pain pathway suppression, this contrast in analgesic effect between MI-EA and OI-EA can be further attributed by the observation that OI-EA can more effectively suppress HP induced contralateral spinothalamic activation as shown in Figure 4. Therefore the overall result from the current study suggests that OI-EA preemptively induces a general sedative effect in supraspinal areas that are essential for HP perception and this effect does not appear to be mediated via any known supraspinal pain modulatory components as none of the prefrontal cortical areas have been activated under OI-EA. In addition, this intensity dependent sedative/suppressing effect is equally effective in both the medial and the lateral pain pathways as the ANOVA showed a significant EA × Intensity interaction effect in several HP related supraspinal areas including the ACC, right SII, insular cortex and premotor cortex.

Previous studies demonstrated that the neuronal analgesic mechanisms of acupuncture were primarily mediated via A-delta afferent fibers which were known to modulate C-fiber mediated HP transmission [33-37]. However, the precise level of neuronal interaction is largely unknown. The result of the current study suggests that EA if given at the optimal non-noxious intensity induces a general supraspinal suppression state which in turn counteracts or minimizes the subsequent excitatory state caused by HP induced C-fiber mediated activation via the spinothalamic and spinopontoamygdaloid pain pathways. This current observation can be added to the overall understanding in the central neuromodulatory mechanisms of EA mediated analgesia.

Conclusion

In short, by adopting an acupuncture treatment paradigm for acute thermal pain and an experimental HP paradigm, we demonstrated the intensity dependent analgesic effect of EA and the corresponding supraspinal activities. EA at optimal non-noxious intensity can generate a general sedative state in supraspinal areas related to both affective and discriminatory aspects of pain perception and this preemptive supraspinal response is essential in preventing the excitatory state associated with subsequent noxious stimuli.

Methods

With the Institution Human Subject Protection Committee approval, healthy volunteers were enrolled for the studies based on the following inclusion and exclusion criteria:

Inclusion criteria:

Age 18 to 80;

Male and female;

No analgesics for the past 2 weeks;

Absence of neuropathic pain states;

Exclusion criteria:

History of psychological illness;

History of claustrophobia;

Lack of ability to understand the experimental protocol and to adequately communicate in English;

Pregnancy;

Pending litigation;

History of head trauma, history of trauma or surgery to lower extremities or low back;

History of any metallic implant in the body as listed in the Institute fMRI Center screening list;

Pre-scanning neurosensory threshold assessment

To be consistent in the study, the location of the thermal thresholds measurement and stimulation was marked at the medial aspect of the left calf between the 6th and 7th marking of an elastic band which consisted of a total of 13 increments, extending from the medial malleolus to the medial tibial plateau. Non-noxious and noxious thermal thresholds including cold and warm, cold and hot pain thresholds were measured by using a Thermal Sensory Analyzer (Medoc Advanced Medical Systems, Minneapolis). This device consisted of a thermode measuring 46 × 29 mm. The temperature of the thermode could either rise or fall (at a rate of 1.2 degrees Celsius/sec for cold and warm sensations, and 3 degrees Celsius/sec for cold and hot pain), depending on the sensations that were being tested. The subject signaled the onset of feeling the tested sensation by pressing a switch, which in turn reversed the temperature change and returned the temperature of the thermode to the 32 degree Celsius baseline. The computer then recorded the temperature of the thermode when the switch was pressed. The average value of testing result was automatically calculated by the computer and displayed on the screen. This method of peripheral sensory testing has been well established in literature and has been used extensively in pain-related studies [8,9,38-40].

FMRI scanning

Subjects were placed comfortably in a supine position in a scanner with their eyes covered by an eye shield. A Facial-cervical Collar Restraint (FCCR) Device was applied to minimize head movement [41]. The following 5 scanning paradigms (Figure 5) were conducted in a random order:

thumbnailFigure 5. The study paradigms with heat pain (HP), two intensities of electrical acupuncture (MI-EA &OI-EA) with and without HP.

1) HP stimuli with pre-measured subject-specific HP thresholds were delivered to the subjects' premarked locations at the left medial calves for 15 seconds with a 60-second of baseline temperature (32°C). The stimulus was repeated four times in a box-cart fashion. The heat pain was provided in an oscillating pattern at the target heat pain threshold ± 0.5°C per oscillation with a total of 5 oscillations per stimulation block to avoid pain wind-up via a fMRI-compatible thermode. At the end of the scanning session, the subjects were asked to rate the overall heat pain score on a VAS score.

2) For the two EA stimulation paradigms, one-inch-long, 36G fMRI-compatible gold plated sterile needles were inserted at the LR1 and SP1 points. The location and method of needle placement used in the study were described in the previously published studies [8,9]. Both the needles and grounding electrodes were linked to a 6-volt ES-160 (Electro-Therapeutic Devices Inc., Markham, Ontario, Canada) clinical acupuncture stimulation device, which consisted of a digital display of the stimulation paradigm. Electrical stimulation was provided at a constant frequency of 5 Hz with a pulse width of 300 microseconds. A stimulation intensity assessment session was conducted prior to the scanning. For the MI-EA stimulation threshold determination, the amplitude of the stimulation was gradually increased from 0 to the level when the subject first noticed the tingling sensation. The process was repeated twice. The average stimulation intensity in a 0-10 range and the patients' average VAS rating of deqi sensation (degree of tingling) were recorded. For the optimal intensity EA (OI-EA) stimulation threshold determination, the amplitude was increased from 0 to a maximally tolerated intensity without any sharp pain or discomfort.

After the sensation of EA has completely subsided, the scanning sessions were then conducted. Two separate scanning runs were conducted in a random order for the two different intensities. EA was given for 15 seconds with 60 seconds of resting period in between stimulations. The stimulation was repeated 4 times to complete the scanning session. At the end of each scanning session, the subjects were asked to rate the intensity of tingling sensation felt during the study.

3) Two separate runs of 15 seconds of EA with two different intensities followed by 15 seconds of HP stimulation at subjects' specific heat pain thresholds and 60 seconds of baseline temperature at 32°C were conducted with 4 repetitions to complete the scanning sequence. The subjects were asked to rate the intensity of overall HP at the end of scanning sessions.

In between each scanning paradigm, a minimal of 15 minutes of washout period was provided to ensure either the HP sensation or the EA related deqi sensation had completely subsided.

FMRI Images were obtained via a 3T GE scanner with T2*- weighted EPI-sequence (TE = 30 ms, TR = 2.0 s, α = 90°, TH = 4 mm, 32 slices, FOV = 220 × 220 mm2, MA = 64 × 64). Two T1-weighted images were acquired: one for spatially normalizing the functional images and the other for anatomical details.

Behavioral data analysis

A paired sample t-test was used to compare the VAS scores of hot pain and de qi sensations.

fMRI data analysis

Each individual subject's functional and anatomical data sets were processed, aligned and prepared in Brain Voyager for within - and between-group random effects analyses based on steps described by Goebel et al.[42].

Preprocessing of functional data

Raw functional data (dicom format) was loaded and converted into Brain Voyager's internal "FMR" data format. Standard sequence of preprocessing steps including slice scan time correction, head motion correction, drift removal and spatial smoothing with Gaussian filter (FWHM = 5 mm) were conducted for each paradigm data set of each subject.

Preprocessing of anatomical data

The anatomical data (dicom format) of each subject was loaded and converted into Brain Voyager's internal "VMR" data format. Intensity inhomogeneities correction was applied and the data was then resampled to 1-mm resolution, and transformed into AC-PC and Talairach standard space. The three spatial transformations were combined and applied backward in one step to avoid quality loss due to successive data sampling. The two affine transformations, iso-voxel scaling and AC-PC transformation, were concatenated to form a single 4 × 4 transformation matrix. For each voxel coordinates in the target (Talairach) space a piece affine "Un-Talairah" step was performed, followed by application of the inverted spatial transformation matrix. The computed coordinates were used to sample the data points in the original 3-D space using sinc interpolation.

Brain segmentation

For 3-D visualization, the brain was segmented from surrounding head tissue using an automatic "brain peeling" tool. The tool analyzes the local intensity histogram in small volumes (20 × 20 × 20 voxels) to define thresholds for an adaptive region-growing technique. This step resulted in the automatic labeling of voxels containing the white and gray matter of the brain, but also other high-intensity head tissue. The next step consisted of a sequence of morphological erosions to remove tissue at the border of the segmented data. By "shrinking" the segmented data, this step separated subparts, which were connected by relatively thin "bridges" with each other. By determining the largest connected component after the erosion step, the brain was separated from other head tissue. Finally, the sequence of erosions was reversed but restricted to voxels in the neighborhood of the largest connected component.

Cortex segmentation

In order to perform a cortex-based data analysis, the gray/white matter boundary was segmented using largely automatic segmentation routines [43]. Following the correction of inhomogeneities of signal intensity across space as described above, the white/gray matter border was segmented with a region-growing method using an analysis of intensity histograms. Morphological operations were used to smooth the borders of the segmented data and to separate the left from the right hemisphere. Each segmented hemisphere was finally submitted to a "bridge removal" algorithm, which ensures the creation of topologically correct mesh representations [43]. The borders of the two resulting segmented subvolumes were tessellated to produce a surface reconstruction of the left and right hemisphere. For better visualization of the areas of activities including those in the sulcus, the resulting meshes were transformed into inflated cortical representations by performing repeated small morphing steps until the central sulcus are visible. The inflated cortical meshes were used as the reference meshes for functional data (maps and time courses) projection. For subsequent cortex-based analysis, the inflated cortical meshes were used to sample the functional data at each vertex (node), resulting in a mesh time course ("MTC") dataset for each run of each subject.

Normalization of functional data

To transform the functional data into Talairach space, the functional time series data of each subject was first coregistered with the subject's 3-D anatomical dataset, followed by the application of the same transformation steps as performed for the 3-D anatomical dataset (see above). This step results in normalized 4-D volume time course ("VTC") data. In order to avoid quality loss due to successive data sampling, normalization was performed in a single step combining a functional-anatomical affine transformation matrix, a rigid-body AC-PC transformation matrix, and a piecewise affine Talairach grid scaling step. As described for the anatomical normalization procedure, these steps were performed backward, starting with a voxel in Talairach space and sampling the corresponding data in the original functional space. In the context of the functional-anatomical alignment, some manual adjustment was necessary to reduce as much as possible the geometrical distortions of the echo-planar images, which exhibited linear scaling in the phase-encoding direction. The necessary scaling adjustment was done interactively using appropriate transformation and visualization tools of Brain Voyager QX.

GLM analysis

For each run of each subject's block, a protocol file (PRT) was derived representing the onset and duration of the events for the different stimulation conditions. In order to account for hemodynamic delay and dispersion, each of the predictors was derived by convolution of an appropriate box-car waveform with a double-gamma hemodynamic response function [44] to extracted brain regions with both positively and negatively correlated blood oxygen level dependent (BOLD) responses. Within group random effect analysis was conducted for each paradigm and areas of activation (positively correlated BOLD) and deactivation (negatively correlated BOLD) were recorded. Between-group random effect analyses were also performed between MI-EA and OI-EA paradigms and a second level 2-factor ANOVA (EA and intensity) was also performed to assess treatment and intensity interaction effect at regions of interest related to HP stimulation.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

AL carried out the experiment, supervised the data analysis and prepared the manuscript. SS carried the data analysis, and graphic and manuscript preparation. AT and JR assisted in conducting the experiment and data analysis. All authors read and approved the final manuscript.

Acknowledgements

The authors thank the National Institute of Health for funding support (NIH NCCAM (1K08-1695-01) and Mr. Martin Bales, a licensed acupuncturist trained in Traditional Chinese Medicine for reviewing this paper.

References

  1. Wu S, Sapru A, Stewart MA, Milet MJ, Hudes M, Livermore LF, Flori HR: Using acupuncture for acute pain in hospitalized children.

    Pediatr Crit Care Med 2009, 10:291-296. PubMed Abstract | Publisher Full Text OpenURL

  2. Lin YC, Tassone RF, Jahng S, Rahbar R, Holzman RS, Zurakowski D, Sethna NF: Acupuncture management of pain and emergence agitation in children after bilateral myringotomy and tympanostomy tube insertion.

    Paediatr Anaesth 2009, 19:1096-1101. PubMed Abstract | Publisher Full Text OpenURL

  3. Sahmeddini MA, Farbood A, Ghafaripuor S: Electro-acupuncture for pain relief after nasal septoplasty: a randomized controlled study.

    J Altern Complement Med 2010, 16:53-57. PubMed Abstract | Publisher Full Text OpenURL

  4. Sim CK, Xu PC, Pua HL, Zhang G, Lee TL: Effects of electroacupuncture on intraoperative and postoperative analgesic requirement.

    Acupunct Med 2002, 20:56-65. PubMed Abstract | Publisher Full Text OpenURL

  5. Helms J: The Basic, Clinical, and Speculative Science of Acupuncture. In Acupuncture Energetics A Clinical Approach for Physicians. Edited by Helms J. Berkeley: Medical Acupuncture Publishers; 1995:19-43. OpenURL

  6. Helms J: The Science of Acupuncture. In Acupuncture Energetics A Clinical Approach for Physicians. Edited by Helms J. Berkeley: Medical Acupuncture Publishers; 1995:19-70. OpenURL

  7. Helms J: The Tendinomuscular Meridian Subsystem. In Acupuncture Energetics A Clinical Approach for Physicians. 1st edition. Edited by Helms J. Berkeley: Medical Acupuncture Publishers; 1995:103-130. OpenURL

  8. Leung A, Khadivi B, Duann JR, Cho ZH, Yaksh T: The effect of Ting point (tendinomuscular meridians) electroacupuncture on thermal pain: a model for studying the neuronal mechanism of acupuncture analgesia.

    J Altern Complement Med 2005, 11:653-661. PubMed Abstract | Publisher Full Text OpenURL

  9. Leung AY, Kim SJ, Schulteis G, Yaksh T: The effect of acupuncture duration on analgesia and peripheral sensory thresholds.

    BMC Complement Altern Med 2008, 8:18. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  10. Leung AY, Park J, Schulteis G, Duann JR, Yaksh T: The electrophysiology of de qi sensations.

    J Altern Complement Med 2006, 12:743-750. PubMed Abstract | Publisher Full Text OpenURL

  11. Brooks JC, Zambreanu L, Godinez A, Craig AD, Tracey I: Somatotopic organisation of the human insula to painful heat studied with high resolution functional imaging.

    Neuroimage 2005, 27:201-209. PubMed Abstract | Publisher Full Text OpenURL

  12. Bowsher D: Representation of somatosensory modalities in pathways ascending from the spinal anterolateral funiculus to the thalamus demonstrated by lesions in man.

    Eur Neurol 2005, 54:14-22. PubMed Abstract | Publisher Full Text OpenURL

  13. Brooks JC, Nurmikko TJ, Bimson WE, Singh KD, Roberts N: fMRI of thermal pain: effects of stimulus laterality and attention.

    Neuroimage 2002, 15:293-301. PubMed Abstract | Publisher Full Text OpenURL

  14. Tracey I, Becerra L, Chang I, Breiter H, Jenkins L, Borsook D, Gonzalez RG: Noxious hot and cold stimulation produce common patterns of brain activation in humans: a functional magnetic resonance imaging study.

    Neurosci Lett 2000, 288:159-162. PubMed Abstract | Publisher Full Text OpenURL

  15. Becerra LR, Breiter HC, Stojanovic M, Fishman S, Edwards A, Comite AR, Gonzalez RG, Borsook D: Human brain activation under controlled thermal stimulation and habituation to noxious heat: an fMRI study.

    Magn Reson Med 1999, 41:1044-1057. PubMed Abstract | Publisher Full Text OpenURL

  16. Price DD: Psychological and neural mechanisms of the affective dimension of pain.

    Science 2000, 288:1769-1772. PubMed Abstract | Publisher Full Text OpenURL

  17. Bernard JF, Besson JM: The spino(trigemino)pontoamygdaloid pathway: electrophysiological evidence for an involvement in pain processes.

    J Neurophysiol 1990, 63:473-490. PubMed Abstract | Publisher Full Text OpenURL

  18. Burstein R, Cliffer KD, Giesler GJ Jr: Direct somatosensory projections from the spinal cord to the hypothalamus and telencephalon.

    J Neurosci 1987, 7:4159-4164. PubMed Abstract | Publisher Full Text OpenURL

  19. Friedman DP, Murray EA, O'Neill JB, Mishkin M: Cortical connections of the somatosensory fields of the lateral sulcus of macaques: evidence for a corticolimbic pathway for touch.

    Journal of Comparative Neurology 1986, 252:323-347. PubMed Abstract | Publisher Full Text OpenURL

  20. Neugebauer V, Li W: Processing of nociceptive mechanical and thermal information in central amygdala neurons with knee-joint input.

    J Neurophysiol 2002, 87:103-112. PubMed Abstract | Publisher Full Text OpenURL

  21. Neugebauer V, Li W: Differential sensitization of amygdala neurons to afferent inputs in a model of arthritic pain.

    J Neurophysiol 2003, 89:716-727. PubMed Abstract | Publisher Full Text OpenURL

  22. Tanimoto S, Nakagawa T, Yamauchi Y, Minami M, Satoh M: Differential contributions of the basolateral and central nuclei of the amygdala in the negative affective component of chemical somatic and visceral pains in rats.

    Eur J Neurosci 2003, 18:2343-2350. PubMed Abstract | Publisher Full Text OpenURL

  23. LeDoux J: The emotional brain, fear, and the amygdala.

    Cell Mol Neurobiol 2003, 23:727-738. PubMed Abstract | Publisher Full Text OpenURL

  24. Bornhovd K, Quante M, Glauche V, Bromm B, Weiller C, Buchel C: Painful stimuli evoke different stimulus-response functions in the amygdala, prefrontal, insula and somatosensory cortex: a single-trial fMRI study.

    Brain 2002, 125:1326-1336. PubMed Abstract | Publisher Full Text OpenURL

  25. Graff-Guerrero A, Gonzalez-Olvera J, Fresan A, Gomez-Martin D, Mendez-Nunez JC, Pellicer F: Repetitive transcranial magnetic stimulation of dorsolateral prefrontal cortex increases tolerance to human experimental pain.

    Brain Res Cogn Brain Res 2005, 25:153-160. PubMed Abstract | Publisher Full Text OpenURL

  26. Hardy SG, Haigler HJ: Prefrontal influences upon the midbrain: a possible route for pain modulation.

    Brain Res 1985, 339:285-293. PubMed Abstract | Publisher Full Text OpenURL

  27. Lorenz J, Minoshima S, Casey KL: Keeping pain out of mind: the role of the dorsolateral prefrontal cortex in pain modulation.

    Brain 2003, 126:1079-1091. PubMed Abstract | Publisher Full Text OpenURL

  28. Wiech K, Seymour B, Kalisch R, Stephan KE, Koltzenburg M, Driver J, Dolan RJ: Modulation of pain processing in hyperalgesia by cognitive demand.

    Neuroimage 2005, 27:59-69. PubMed Abstract | Publisher Full Text OpenURL

  29. Hui KK, Marina O, Liu J, Rosen BR, Kwong KK: Acupuncture, the limbic system, and the anticorrelated networks of the brain.

    Auton Neurosci 2010. OpenURL

  30. Asghar AU, Green G, Lythgoe MF, Lewith G, MacPherson H: Acupuncture needling sensation: the neural correlates of deqi using fMRI.

    Brain Res 2010, 1315:111-118. PubMed Abstract | Publisher Full Text OpenURL

  31. Hui KK, Marina O, Claunch JD, Nixon EE, Fang J, Liu J, Li M, Napadow V, Vangel M, Makris N, et al.: Acupuncture mobilizes the brain's default mode and its anti-correlated network in healthy subjects.

    Brain Res 2009, 1287:84-103. PubMed Abstract | Publisher Full Text OpenURL

  32. Fang J, Jin Z, Wang Y, Li K, Kong J, Nixon EE, Zeng Y, Ren Y, Tong H, Wang P, Hui KK: The salient characteristics of the central effects of acupuncture needling: limbic-paralimbic-neocortical network modulation.

    Hum Brain Mapp 2009, 30:1196-1206. PubMed Abstract | Publisher Full Text OpenURL

  33. Chung JM, Lee KH, Hori Y, Endo K, Willis WD: Factors influencing peripheral nerve stimulation produced inhibition of primate spinothalamic tract cells.

    Pain 1984, 19:277-293. PubMed Abstract | Publisher Full Text OpenURL

  34. Chung JM, Fang ZR, Hori Y, Lee KH, Willis WD: Prolonged inhibition of primate spinothalamic tract cells by peripheral nerve stimulation.

    Pain 1984, 19:259-275. PubMed Abstract | Publisher Full Text OpenURL

  35. Yarnitsky D, Ochoa JL: Warm and cold specific somatosensory systems. Psychophysical thresholds, reaction times and peripheral conduction velocities.

    Brain 1991, 114:1819-1826. PubMed Abstract | Publisher Full Text OpenURL

  36. Verdugo R, Ochoa JL: Quantitative somatosensory thermotest. A key method for functional evaluation of small calibre afferent channels.

    Brain 1992, 115:893-913. PubMed Abstract | Publisher Full Text OpenURL

  37. Wahren LK, Torebjork E, Jorum E: Central suppression of cold-induced C fibre pain by myelinated fibre input.

    Pain 1989, 38:313-319. PubMed Abstract | Publisher Full Text OpenURL

  38. Leung A, Li E, Fallah A, Schulteis G, Novak E, Duann JR, Yaksh TL: The effect of needle combination on the analgesic efficacy of the tendinomuscular meridians(TMM) systems.

    Medical Acupuncture 2007, 19(4):191-199. Publisher Full Text OpenURL

  39. Leung A, Wallace MS, Ridgeway B, Yaksh T: Concentration-effect relationship of intravenous alfentanil and ketamine on peripheral neurosensory thresholds, allodynia and hyperalgesia of neuropathic pain.

    Pain 2001, 91:177-187. PubMed Abstract | Publisher Full Text OpenURL

  40. Leung AY, Wallace MS, Schulteis G, Yaksh TL: Qualitative and quantitative characterization of the thermal grill.

    Pain 2005, 116:26-32. PubMed Abstract | Publisher Full Text OpenURL

  41. Leung A, Duann JR, Davis M, Li E, Fallah A, Yaksh TL: Facial-cervical Collar Restraint (FCCR) Device in reducing Head Motion during a Noxious Stimulus Study.

    Neuroimage 2006, 31:504. OpenURL

  42. Goebel R, Esposito F, Formisano E: Analysis of functional image analysis contest (FIAC) data with brainvoyager QX: From single-subject to cortically aligned group general linear model analysis and self-organizing group independent component analysis.

    Hum Brain Mapp 2006, 27:392-401. PubMed Abstract | Publisher Full Text OpenURL

  43. Kriegeskorte N, Goebel R: An efficient algorithm for topologically correct segmentation of the cortical sheet in anatomical mr volumes.

    Neuroimage 2001, 14:329-346. PubMed Abstract | Publisher Full Text OpenURL

  44. Friston KJ, Fletcher P, Josephs O, Holmes A, Rugg MD, Turner R: Event-related fMRI: characterizing differential responses.

    Neuroimage 1998, 7:30-40. PubMed Abstract | Publisher Full Text OpenURL