Abstract

The linkage between brain response to acupuncture and subsequent analgesia remains poorly understood. Our aim was to evaluate this linkage in chronic pain patients with carpal tunnel syndrome (CTS). Brain response to electroacupuncture (EA) was evaluated with functional MRI. Subjects were randomized to 3 groups: (1) EA applied at local acupoints on the affected wrist (PC-7 to TW-5), (2) EA at distal acupoints (contralateral ankle, SP-6 to LV-4), and (3) sham EA at nonacupoint locations on the affected wrist. Symptom ratings were evaluated prior to and following the scan. Subjects in the local and distal groups reported reduced pain. Verum EA produced greater reduction of paresthesia compared to sham. Compared to sham EA, local EA produced greater activation in insula and S2 and greater deactivation in ipsilateral S1, while distal EA produced greater activation in S2 and deactivation in posterior cingulate cortex. Brain response to distal EA in prefrontal cortex (PFC) and brain response to verum EA in S1, SMA, and PFC were correlated with pain reduction following stimulation. Thus, while greater activation to verum acupuncture in these regions may predict subsequent analgesia, PFC activation may specifically mediate reduced pain when stimulating distal acupoints.

1. Introduction

Acupuncture, a component of traditional Chinese medicine, has been commonly applied to alleviate symptoms of patients with chronic pain [1]. Carpal tunnel syndrome (CTS) is mainly driven by partial deafferentation secondary to compression of the median nerve within the carpal tunnel [2]. CTS clinically manifests as slowing of median nerve conduction velocity, local pain, and paresthesia. Recent randomized controlled trials (RCT) for CTS have shown that acupuncture produced significant improvement in symptoms, with effects similar to steroid treatment [3] and night splinting [4]. A recent RCT demonstrated that acupuncture also reduced CTS symptoms significantly greater than placebo [5].

Noninvasive brain imaging techniques, such as functional MRI (fMRI), have offered an unprecedented window into how the human brain responds to acupuncture needle stimulation [6, 7]. However, very few of these studies have been performed in chronic patient populations, and even fewer have evaluated how brain response to acupuncture relates to the alleviation of clinical symptoms or even evoked-pain ratings. Zhang et al. demonstrated that fMRI response to transcutaneous electrical acupoint stimulation in brain areas such as secondary somatosensory cortex (S2), insula, and primary motor cortex (M1) was associated with reduced heat pain ratings in healthy adults [8]. More recently, Yang et al. found that acupuncture increased metabolism in regions including the prefrontal cortex (PFC), insula, and cingulate, while also decreasing pain levels in acute migraine patients in a positron emission tomography (PET) study [9]. Harris et al. used PET with 11C-carfentanil in fibromyalgia patients and found that long-term increases in resting mu-opioid receptor binding in regions including insula, cingulate, and basal ganglia following 4 weeks of acupuncture correlated with decreased clinical pain levels over this same time period [10]. In CTS patients, compared to healthy adults, manual acupuncture needling has been found to produce more robust fMRI response in several brain areas including insula, cingulate, S1, and PFC, when controlling for the effects of sham (noninserted cutaneous tactile) acupuncture [11]. Enhanced processing in S1 and PFC was particularly interesting, given the fact that these patients demonstrated altered somatosensory-stimulation-evoked brain response in these areas [12] and that S1 activity was specifically modulated by five weeks of acupuncture treatment [13]. Unfortunately, changes in clinical pain were not evaluated in this study. Thus, the association between the brain circuitry processing acupuncture stimulation and postacupuncture clinical outcomes such as pain reduction is currently unknown.

In this cross-sectional study, CTS subjects were randomized to 3 groups: (1) EA applied at local acupoints on the affected wrist (PC-7 to TW-5), (2) EA at distal acupoints (contralateral ankle, SP-6 to LV-4), and (3) sham EA at nonacupoint locations on the affected wrist. In addition to fMRI data acquired during EA, we also evaluated changes in clinical symptoms following acupuncture and correlated fMRI response to EA with changes in symptoms. We hypothesized that both local and distal EA would produce greater symptom reduction compared to sham acupuncture. Furthermore, we hypothesized that the magnitude of symptom reduction would correlate with activation in brain regions previously associated with somatosensory, affective, and cognitive processing of pain and paresthesia in CTS subjects.

2. Methods

2.1. Subjects

Subjects, aged 20 to 60, with a 3-month or greater history of pain and/or paresthesia in median-nerve-innervated areas were enrolled. All subjects were examined for eligibility by a physiatrist at Spaulding Rehabilitation Hospital, which included a physical exam for Phalen’s maneuver [14] and Durkan’s sign [15] and testing of median and ulnar sensory nerve conduction (NCS: Cadwell Sierra EMG/NCS Device, Kennewick, WA). NCS inclusion criteria consisted of median nerve sensory latency greater than 3.7 ms or median nerve sensory latency greater than 0.5 ms compared to ulnar nerve. Exclusion criteria consisted of contraindications to MRI, history of diabetes mellitus, cardiovascular, respiratory, or neurological illnesses, rheumatoid arthritis, wrist fracture with direct trauma to median nerve, current usage of prescriptive opioid medication, thenar atrophy, previous acupuncture treatment (manual, EA, and TENS) for CTS, nerve entrapment other than median nerve, cervical radiculopathy or myelopathy, generalized peripheral neuropathy, blood dyscrasia or coagulopathy or current use of anticoagulation therapy. History of axis I psychiatric diagnosis (substance use disorder, psychotic disorder, or bipolar disorder), and use of psychotropic medications were also exclusions for this study. Chronic symptomatology for all eligible subjects was evaluated using the Boston Carpal Tunnel Syndrome Questionnaire (BCTSQ) [16].

A total of 59 CTS subjects ( years old, mean ± SD, 49 Female) were enrolled in this study. For the subjects who had diagnosed bilateral CTS, the more affected hand was determined as the test hand. CTS subjects were randomized to one of the three study arms: (1) local verum EA ( , 17F, 14 right hand affected), (2) distal verum EA ( , 14F, 13R), and (3) sham EA ( , 18F, 11R). All study protocols were approved by the Massachusetts General Hospital and Partners Human Research Committee. Written informed consent was obtained from all subjects.

2.2. Acupuncture Procedure

For local verum EA, MRI-compatible titanium needles (0.2 mm in diameter, 35–50 mm in length, DongBang Acupuncture Inc. Boryeong, Korea) were inserted and deqi sensation elicited at acupoints PC7 (pericardium 7, 1st wrist crease) and TW5 (triple-warmer 5, dorsal aspect of forearm), local to the more affected hand (Figure 1(a)). PC-7 was chosen because it is close to the CTS lesion, and this set of points was found to reduce pain and paresthesia in our previous study [13].

For distal verum EA, MRI-compatible titanium needles were inserted and deqi sensation elicited at acupoints SP6 (spleen 6, medial aspect of lower leg) and LV4 (liver 4, anterior aspect of the ankle) at the ankle on the contralateral side to the more affected hand (Figure 1(b)). Distal acupoints LV4 and SP6 were chosen based on mirror point methods common in acupuncture practice, where acupoints on the leg/ankle can be used to treat symptoms on the opposite arm/wrist [17].

For sham EA, MRI-compatible blunt-tipped acupuncture needles were placed with a single tap but not inserted percutaneously, over sham points, SH1 (2 cun, or roughly 2-3 cm, distal and slightly volar to acupoint SI-7, which is 5 cun proximal to ulnar edge of the transverse wrist crease, ulnar forearm) and SH2 (1 cun distal and slightly volar to SI-7), on the more affected hand (Figure 1(a)).

For all 3 groups, needles were connected to a constant current electroacupuncture (EA) device (HANS LH202H, Neuroscience Research Center, Peking University, Beijing, China). A licensed acupuncturist trained to place and stimulate acupuncture needles in the scanner performed these procedures. For verum EA, current stimulation frequency was set to 2 Hz, while the current intensity was set just prior to the functional scan. The acupuncturist gradually increased the current intensity until the subject felt a moderately strong but not painful sensation. For sham EA, electrodes were attached to the needles, but no electrical current was passed. Subjects were instructed that the current intensity was set to a predetermined level and that they may or may not feel any sensation at the needle sites. All subjects were instructed to close their eyes and focus on the acupuncture stimulation while remaining as still as possible.

2.3. Data Acquisition

All imaging data were acquired on a 3.0T Siemens Trio (Siemens Medical, Erlangen, Germany) equipped with 32-channel head coil. Structural imaging data were acquired with a multiecho MPRAGE T1-weighted pulse sequence (TR = 2530 ms, TE1/TE2 = 1.64/30.0 ms, TI = 1200 ms, flip angle = 7°, FOV = , slices = 176, sagittal acquisition, spatial resolution =  mm3).

Functional imaging (fMRI) data were acquired using a gradient echo BOLD T2*-weighted pulse sequence (TR/TE = 2000/30 ms, FOV =  mm, 32 axial slices parallel to anterior/posterior commissural plane, voxel size = 3.125 × 3.125 × 3.6 mm, flip angle = 90°). Subjects lay supine in the scanner with earplugs to attenuate acoustic gradient switching noise. For verum EA, stimulation was performed using an event-related design similar to our previously published approach for manual acupuncture [18] (2-second stimulation events with randomized ISI, 6–12 seconds, and total scan time 5 minutes and 6 seconds, Figure 1(c)). For sham EA, procedures were identical, but no electricity was passed through the needles.

Symptoms were assessed using a 0–10 VAS scale for both pain and paresthesia (tingling) at the hand/wrist. The scale ranged from none (0) to unbearable (10) and was administered at the beginning and the end of the MRI session. In addition, subjects were asked to rate the intensity of acupuncture-evoked sensations after the scan session using the MGH Acupuncture Sensation Scale (MASS) instrument [19].

2.4. Data Analysis

Statistical analyses for behavioral and clinical data were performed with SPSS (SPSS version 10.0.7, Chicago, Il). Median sensory nerve velocities and motor nerve latencies were compared with those of the ulnar nerve in all CTS subjects using a Student’s -test, significant at .

Changes in VAS scores for pain and paresthesia were compared between groups using a mixed model ANOVA with interaction of Group (local, distal, and sham) × Time (pre and post). Student’s -test was used to compare change in VAS scores between verum (combined local and distal) and sham groups. EA current intensities were compared between local and distal acupoint groups using a Student’s -test significant at . A one-way ANOVA was performed in order to compare evoked acupuncture sensations, as well as the MASS Index (a composite metric of deqi sensation) between groups. Post hoc testing was performed with the Tukey test.

FMRI data were preprocessed using the FMRIB software Library (FSL v.4.1), Freesurfer (v.5.1.), and AFNI (v.2.). FMRI data were coregistered with each subject’s structural MRI data using boundary-based registration (BB registration, Freesurfer [20]). Preprocessing included slice timing correction, motion correction, high pass filtering with a cut-off period of 50 sec, and spatial smoothing with Gaussian kernel at FWHM = 5 mm (Feat, FSL). Preprocessed fMRI data were then analyzed using a general linear model for all subjects (Feat, FSL), with the explanatory variable set by the event-related design. The resultant parameter estimates and variances from all subjects were transformed to standard MNI space (FNIRT, FSL) in order to perform nonflipped (compare with below) group analyses. Registration was ensured by visualization (AFNI, afni).

In order to better assess brain response for structures known to be lateralized relative to somatosensory input (i.e., S1 and M1, and thalamus), parameter estimates and variances of subjects whose more affected hand was the left hand (and thus experienced local EA on the left hand or distal EA on the right ankle) also had their fMRI parameter estimates flipped across the midsagittal plane before passing them up to the flipped group analysis. Therefore, a total of 21 fMRI datasets (local: 8, distal: 5, sham: 8) were analyzed by flipping across the midsagittal plane.

Group maps were calculated using a mixed effects statistical model (FLAME, Feat, FSL). Difference maps were calculated with an ANOVA using a mixed effect model (FLAME, Feat, FSL). As we found no differences in brain response between local and distal groups, a combined group map for “verum EA” was also calculated in order to increase statistical power and for future testing. Whole brain regression analysis was performed for local, distal, and sham groups, as well as the combined verum group in order to identify brain regions associated with symptom reduction. For the regression analysis with symptom reduction, changes in VAS pain and paresthesia scores (post and pre) were demeaned and added as explanatory variable to the model. All statistical maps were thresholded with cluster forming threshold at (voxel wise threshold ), and cluster corrected for multiple comparisons at [21].

3. Results

3.1. Demographic and Clinical Features

Bilateral CTS was diagnosed in 44/59 (75%) subjects, while unilateral CTS was diagnosed in 15/59 (25%) subjects. The more affected hand was the right hand in 38/59 (64.4%) subjects and the left hand in 21/59 (35.6%) subjects. Pain was the more severe symptom in 11/59 (18.3%) subjects; paresthesia was more severe in 37/59 (63.3%) subjects, and pain and paresthesia were observed with equal severity in 11/59 (18.3%) subjects. BCTSQ assessment of symptoms on the scale 1 to 5 demonstrated that pain and paresthesia were moderate ( , , mean ± S.D., resp.). Pain and paresthesia ratings were positively correlated ( , ); that is, those subjects with greater pain also reported greater paresthesia. Self-reported symptom duration was years (mean ± S.D.) and was positively correlated with subjects’ age ( , ). A significant correlation was found between VAS pain and paresthesia scores before acupuncture (pre; , ). In addition, VAS pain scores at baseline were significantly correlated with the BCTSQ pain score ( , ). However, VAS paresthesia scores were not significantly correlated with BCTSQ paresthesia scores ( , ).

Phalen’s test was positive on the right hand in 39/59 (66.1%) subjects and on the left hand in 34/59 (57.6%) subjects. Durkan’s test was positive on the right hand in 26/59 (44.0%) subjects and on the left hand in 24/59 (40.7%) subjects. For the study hand, Phalen’s test was positive in 51/59 (86.4%) subjects, while Durkan’s test was positive in 38/59 (64.4%) subjects. Median nerve sensory velocities were significantly slower compared to ulnar nerve sensory velocities (median:  m/s, ulnar: , mean ± SD, ). Furthermore, median nerve motor latencies were significantly longer compared to ulnar nerve motor latencies (median:  ms, ulnar: , mean ± S.D., ).

The regular usage of night splints was reported in 33/59 (56.7%) subjects. Subjects’ occupational status could be described as “full-time work” in 38/59 (64.4%) subjects and part-time in 8/59 (13.6%) subjects. Mean body mass index (BMI) was (mean ± SD).

3.2. Symptom Change following EA

We found a significant main effect of Time (pre versus post) for VAS score in pain and paresthesia ( , , , , resp., Figure 2). There was neither a significant main effect of Group ( , , resp.) nor interaction of Time × Group for either pain or paresthesia ( , , resp.). Post hoc testing revealed that VAS scores for pain showed significant reductions for local and distal groups but not for the sham group (local: , ; distal: , ; sham: , ; mean ± SD). VAS scores for paresthesia showed significant reductions for the local group but not for the distal or sham groups (local: , ; distal: , ; sham: , ; mean ± SD). Also, the combined verum (local and distal) group produced a greater reduction in VAS scores for paresthesia compared to the sham group (verum: , sham: , mean ± S.D., ). A statistically significant reduction for verum compared to sham acupuncture was not seen in VAS scores for pain (verum: , sham: , mean ± SD, ).

We also found that change in VAS pain scores correlated with change in VAS paresthesia scores for both the local ( , ) and sham groups ( , ) but not for the distal group ( , ).

3.3. Electroacupuncture Current Intensity and EA-Evoked Sensations

EA current intensity did not differ between local and distal groups (local:  mA, distal: , mean ± S.D., ). For acupuncture-evoked sensation, an ANOVA demonstrated a significant main effect of Group for Mass index (MI, , , Figure 3). Post hoc testing revealed that the local group produced greater MI scores compared to the sham group ( , local: , distal: , sham: , mean ± S.D.). In regard to individual MASS sensations, a significant main effect of Group was also detected for soreness ( , ), aching ( , ), deep pressure ( , ), tingling ( , ), dull pain ( , ), sharp pain ( , ), and unpleasantness ( , ) (Figure 3). Post hoc testing revealed that the local group produced greater sensation compared to the sham group for soreness, aching, deep pressure, tingling, sharp pain, spreading, and unpleasantness ( ). Also, the distal group produced greater sensation compared to sham acupuncture for soreness and sharp pain ( ), while the local group produced greater sensation compared to the distal group for dull pain and unpleasantness ( , Figure 3). EA current intensity was correlated with numbness for the local group ( , ). EA current intensity was correlated with soreness ( , ) and cold ( , ) for the distal group. In addition, Mass index was not correlated with pain reduction (local: , distal: , sham: ) or paresthesia reduction (local: , distal: , sham: ) in any group ( ). EA current intensity was not correlated with pain reduction (local: , distal: ) or paresthesia reduction (local: , distal: ) in either the local or distal group ( ).

3.4. fMRI during Acupuncture

Local EA produced activation in bilateral insulae, secondary somatosensory cortex (S2), superior temporal gyrus (STG), and contralateral postcentral gyrus (S1), with deactivation in cuneus and ipsilateral S1 (Table 1, Figure 4). Distal EA produced activation in bilateral insulae, S2, and premotor cortex (PMC), while no deactivation was found (Figure 4). Sham EA showed no significant fMRI response (Table 1). There were no significant differences in brain response between local and distal groups (Table 1); hence we also calculated a combined verum (local and distal) EA group map, which showed activation in bilateral S2, premotor (PMC), supramarginal (SMG), and superior temporal (STG) gyri, as well as anterior and posterior insulae (a.Ins, p.Ins) and thalamus. Activation was also noted in right presupplementary motor area (pSMA) and middle and inferior frontal gyri (MFG, IFG) (Table 1, Figure 4). Verum EA also produced deactivation in bilateral occipital gyrus, right medial prefrontal cortex (PFC), inferior temporal gyrus (ITG), and S1 (Table 1, Figure 4).

In comparison to fMRI response for the sham group, local EA produced greater activation in bilateral S2 and right frontal insular cortex (FIC), while greater deactivation was found in right inferior temporal gyrus (ITG) (Table 2, Figure 4). Distal EA produced greater activation in bilateral S2 and deactivation in posterior cingulate cortex (PCC). Combined (local and distal) verum EA produced greater activation in bilateral S2 and insula and deactivation in PCC, precuneus, and right S1.

Whole brain regression analysis for change in VAS pain score in the distal group revealed a significant correlation in right prefrontal cortex (PFC). Thus, greater activation in PFC was associated with greater pain reduction following distal EA (Table 3, Figure 5). No significant correlations were found for either local or sham groups. For the combined verum EA group, a significant negative correlation was found between changes in VAS pain scores and brain activity in right S1 and bilateral supplemental motor area (SMA) and prefrontal cortex (PFC) (Table 3, Figure 6). Thus, greater activation in S1, SMA, and PFC was associated with greater pain reduction following verum EA.

4. Discussion

This study investigated how brain response to EA was associated with symptom reduction following stimulation in CTS subjects. Our main finding was that pain was reduced following verum (both local and distal) EA and that greater brain activation in PFC, SMA, and S1 in response to verum EA was associated with more pronounced pain reduction. Thus, greater activation in these regions may be a biomarker for immediate analgesia following EA.

Brain response to verum EA produced activation in several regions including dorsolateral prefrontal cortex (dlPFC), pre-SMA (pSMA), S1, bilateral S2, and insula, with the latter two regions also showing greater activation compared to sham EA. Deactivation was noted in ipsilateral S1 and default mode network areas such as the medial prefrontal cortex and lateral temporal cortex, also with evidence of greater deactivation compared to sham EA. These findings are consistent with multiple previous acupuncture fMRI studies [6, 7], which support the veracity of the event-related fMRI data used to correlate with clinical outcomes in this study.

Similar to our findings, linking dorsolateral PFC (dlPFC) activation with pain reduction following EA, Yang et al. used FDG PET and found that acupuncture produced increased brain metabolism in the middle frontal gyrus (a subregion of the dlPFC) and concomitant pain reduction in acute migraine patients [9]. The prefrontal cortex is known to modulate pain [22]. Chronic pain patients demonstrated reduced dlPFC gray matter volume [23], while transcranial magnetic stimulation to this region reduced placebo analgesia [24]. Moreover, recent studies have suggested that PFC dopamine levels mediate pain sensitivity [25]. As EA has been shown to modulate fMRI activity in dopaminergic source regions (i.e., substantia nigra) in a time-dependent manner [26], analgesia related with PFC response to EA in CTS subjects may prove to be mediated by this neuromodulatory catecholamine.

PFC response may also differentiate the underlying analgesic mechanisms for local versus distal EA. Both local and distal EA reduced pain, and we found no differences in brain response between these groups. However, this result does not necessarily suggest that there is no acupoint specificity for brain-linked mediators of acupuncture analgesia in CTS. For instance, greater PFC activation to distal, but not local EA, was associated with greater pain reduction. The tissue at local acupoints (i.e., PC-7) is adjacent to the carpal tunnel, which is known to be perturbed by the increased pressure, fibrosis, swelling, and various biochemical changes within the carpal tunnel [27], while the tissue at distal acupoints is not. Thus, the mechanism by which local EA reduces pain may involve changes in local signaling from the wrist lesion (bottom, up), while the mechanism by which distal EA reduces pain may instead involve changes in brain processing (top, down), particularly in the PFC. Interestingly, multiple previous studies have found that PFC activity supports placebo [24, 28] and expectation-mediated analgesia [29]. Future studies should further explore how PFC activity relates to analgesia for different forms of acupuncture and how this is similar to or differs from contextually-mediated analgesic phenomena.

Brain response to verum EA in right S1 was also correlated with pain reduction. We found that greater activation in this region was associated with reduced pain, while more pronounced deactivation was associated with worsening pain following verum EA. The somatotopy of this cluster is consistent with the hand area, which would be ipsilateral for the local group and contralateral (but well outside the leg area) for the distal group. In general, S1 regions outside of the contralateral somatotopic representation for somatosensory stimulated body areas are deactivated [30], and this is clear in the scatterplot for our data as well deactivation is noted in this right S1 cluster for most subjects. Interhemispheric communication modulates subcortical relay of sensory information [31]. Therefore, the balance of left and right S1 activity may modulate the amount of sensory input, which might serve to diminish the spontaneous afferent signal (e.g., pain, paresthesia) from the affected hands in our study.

SMA activation was also correlated with pain reduction following verum EA. SMA is a cortical region that modulates communication between the somatosensory and motor systems and has been shown to be activated by painful stimulation in fMRI [32] and involved in pain control [33]. Greater activation following verum EA may reflect greater transfer of EA-induced somatosensory inputs to the motor system [34], fostering a more normalized sensorimotor communication, compared to the sporadic afference coming from diffuse paresthesia. Thus, EA (whether local or distal) may also reduce pain by supplying regulated somatosensory input to the brain.

While verum acupuncture was found to reduce paresthesia significantly more than sham acupuncture, the same was not true for pain. Firstly, these results suggest that compared to pain, paresthesia may be less susceptible to modulation by sham acupuncture. In fact, paresthesia is a hallmark of CTS and other peripheral neuropathic pain disorders stemming from compression of the nerve trunk, and may be more dependent on peripheral factors such as hand positioning and temperature. Our data suggest that reduction of paresthesia is dependent on therapies with significant somatosensory afference, such as verum EA, regardless of whether this afference comes from the site of the lesion. Perception of paresthesia in CTS may be less centralized compared to pain and may thus be less amenable to placebo effects. However, we should also note that controversy regarding the use of sham acupuncture as a control for verum acupuncture exists, as sham procedures are not physiologically inert [35, 36]. In fact, sham acupuncture has been shown to reduce aversive symptom more readily than a placebo pill [37]. Interestingly, while analgesic outcomes may be similar between verum and sham EA, the brain mechanisms supporting this analgesia may differ substantially. In our study, brain response was more profound for verum compared to sham acupuncture, and activity in specific regions in response to verum EA was correlated with pain reductions. These results further support the growing evidence [10] that brain response may serve as an objective marker that differentiates verum and sham acupuncture more readily than subjective pain report does.

Several limitations to our study should be noted. Although brain response to verum or distal EA showed significant correlation with short-term pain reduction, we did not have enough statistical power to show significant correlations between brain response to local or sham EA alone and short-term pain reduction. In addition, our study examined brain response to EA stimulation associated with short-term clinical outcomes, and future studies should extend our analysis to longer-term outcomes. Future studies also should apply alternative fMRI approaches, such as functional connectivity [38], to evaluate how brain connectivity response to both verum and sham acupuncture mediates analgesia in CTS.

In conclusion, EA stimulation at both local and distal acupoints reduced pain, and greater brain activation in PFC, SMA, and S1 was associated with more pronounced pain reduction. Thus, greater PFC, SMA, and S1 activation may support acupuncture analgesia. Brain response in PFC, SMA, and S1 could serve as predictive biomarkers to identify patients more likely to benefit from acupuncture therapy.

Acknowledgment

This work was supported by NCCAM, National Institutes of Health (R01-AT004714, R01-AT004714-02S1, P01-AT002048).