Open Access

The modulation effect of longitudinal acupuncture on resting state functional connectivity in knee osteoarthritis patients

  • Xiaoyan Chen1,
  • Rosa B. Spaeth1,
  • Sonya G. Freeman1,
  • Donna Moxley Scarborough2,
  • Javeria A. Hashmi1,
  • Hsiao-Ying Wey1,
  • Natalia Egorova1,
  • Mark Vangel3, 4,
  • Jianren Mao5,
  • Ajay D. Wasan7,
  • Robert R. Edwards6,
  • Randy L. Gollub1, 3 and
  • Jian Kong1, 3Email author
Contributed equally
Molecular Pain201511:67

https://doi.org/10.1186/s12990-015-0071-9

Received: 22 June 2015

Accepted: 15 October 2015

Published: 29 October 2015

Abstract

Recent advances in brain imaging have contributed to our understanding of the neural activity associated with acupuncture treatment. In this study, we investigated functional connectivity across longitudinal acupuncture treatments in older patients with knee osteoarthritis (OA). Over a period of 4 weeks (six treatments), we collected resting state functional magnetic resonance imaging (fMRI) scans from 30 patients before and after their first, third and sixth treatments. Clinical outcome showed a significantly greater pain subscore on the Knee Injury and Osteoarthritis Outcome Score (KOOS) (indicative of improvement) with verum acupuncture than with sham acupuncture. Independent component analysis (ICA) of the resting state fMRI data showed that the right frontoparietal network (rFPN) and the executive control network (ECN) showed enhanced functional connectivity (FC) with the rostral anterior cingulate cortex/medial prefrontal cortex, a key region in the descending pain modulatory system, in the verum groups as compared to the sham group after treatments. We also found that the rFPN connectivity with the left insula is (1) significantly associated with changes in KOOS pain score after treatments, and (2) significantly enhanced after verum acupuncture treatments as compared to sham treatment. Analysis of the acupuncture needle stimulation scan showed that compared with sham treatment, verum acupuncture activated the left operculum/insula, which also overlaps with findings observed in resting state analysis. Our results suggest that acupuncture may achieve its therapeutic effect on knee OA pain by modulating functional connectivity between the rFPN, ECN and the descending pain modulatory pathway.

Clinical trial number: NCT01079390

Keywords

Acupuncture Knee osteoarthritis Chronic pain Resting state fMRI Independent component analysis (ICA) Knee injury and Osteoarthritis Outcome Score (KOOS) Right frontoparietal network Executive control network

Background

Osteoarthritis (OA) is a major public health problem among the elderly and is associated with considerable disability [1]. The knee is one of the most common body locations for OA. A recent analysis of data from the National Health and Nutrition Examination Survey III indicated that about 35 % of women and men aged 60 years and above had radiographic knee OA [2].

OA is a complex chronic pain condition due in part to both nociceptive and neuropathic mechanisms [3]. Previous studies have demonstrated central sensitization in OA patients [46] and increased activity in the periaqueductal grey (PAG) is associated with stimulation of the skin in referred pain areas of patients [6], indicating a pathology change in the central nervous system of OA patients.

Despite the high prevalence rate of OA, the treatment of OA is far from satisfactory [7, 8]. Pharmacological treatment of knee OA is often ineffective with unwanted and dangerous side effects [8, 9].

Arguably, acupuncture may be a promising treatment option for knee OA due to the effectiveness of the pain relief it provides [8, 1014], and the rarity of adverse effects [1518]. Nevertheless, the mechanisms underlying the effects of acupuncture treatment in knee OA patients are still poorly understood. As a unique treatment modality, studies have shown that acupuncture may produce an analgesic effect through the endogenous descending pain modulatory system [1922].

Brain imaging studies have shown that acupuncture needle stimulation [2327] can evoke widespread brain activity changes, and modulate the functional connectivity (FC) of the pain processing network [2835], which opens a new window for understanding the central mechanism of acupuncture treatment. In addition, investigators have found abnormal or disrupted FC in patients with various chronic pain disorders [3641], and that FC can be modulated by analgesics in both healthy [42] and chronic pain patients [43]. These studies suggested that FC may be a useful tool in both acupuncture and chronic pain research [34].

In the present study, we investigated changes in FC in knee OA patients across longitudinal acupuncture treatments. One challenge of acupuncture research is variability in the prescription of acupoints for the same disorders, i.e., the number of acupuncture acupoints (needles) used during treatment. As acupuncture is an invasive treatment, fewer needles are always preferred in acupuncture practice. In this study, in addition to comparing the effect/mechanism of real and sham acupuncture treatment, we also attempted to explore the effects of the number of acupoints after controlling for the total stimulation time by dividing the verum acupuncture group into two subgroups: low dose (two local acupoints around the knee), and high dose (two local acupoints around the knee plus an additional four local and distant acupoints).

Knee OA patients were randomized to receive either high dose verum acupuncture, low dose acupuncture or sham acupuncture. Using a longitudinal treatment design [8], each patient received 6 acupuncture treatments over 4 weeks. We measured brain activity during acupuncture needle stimulation and resting state functional connectivity using functional magnetic resonance imaging (fMRI) during the first, third and sixth treatment (Additional file 1: Figure S1a). We also administered the Knee Injury and Osteoarthritis Outcome Score (KOOS) [44] before and after the six-session acupuncture treatment period.

Results

Forty-four acupuncture naïve patients (19 females) aged 43–70 with a diagnosis of chronic osteoarthritis in the right and/or left knee participated in the study. Of the 44 patients who enrolled, 30 (13 females) completed all study procedures. Fourteen subjects did not complete the study due to problems with scheduling (6), ineligibility at screening (3), disinterest (2), claustrophobia (1), or inability to adhere to study requirements in the scanner (2). Of the 14 subjects who dropped out of the study, four dropped out after randomization (2 from the low dose verum acupuncture group and 2 from the high dose verum acupuncture group).

Sensation evoked by acupuncture treatment, clinical outcomes and blindness

As a treatment quality control, we measured the deqi sensation evoked by acupuncture treatment using the MASS [45, 46]. The average total MASS score (sum of the intensities of each sensation) differed significantly across the acupuncture treatment groups [high (12.95 ± 7.55) vs. low (15.22 ± 11.90) vs. sham (4.82 ± 3.79)] [F(2,27) = 4.21, p = 0.026]. Post hoc analysis showed the following differences: high vs low p = 0.55, high vs sham p = 0.04; low vs sham p = 0.01. There was no significant difference between the high and low dose acupuncture groups. A previous study by our group reported more details regarding the deqi sensations evoked by different treatment modalities in this sample [45].

We measured clinical outcomes using the KOOS [44], which is comprised of 5 subscales: (1) pain, (2) other symptoms, (3) function in daily living, (4) function in sport and recreation, and (5) knee-related quality of life at the beginning and end of treatment (Table 1). Because this study focused on brain networks involved in the pain process, we only reported the comparisons of KOOS pain changes before and after treatments in different groups. Please see our previously published paper on the comparison of other outcomes [45]. We found that both high and low dose verum acupuncture groups showed significant improvement or a trend (post- minus pre-treatment) as compared with sham acupuncture [p = 0.07 (high dose vs sham), p = 0.04 (low dose vs sham)]. The improvement in the two real acupuncture treatment groups was similar (p = 0.90, low dose vs high dose). Pooling the data from two real acupuncture groups using repeated measurements showed that there was a significant interaction between acupuncture mode (real vs. sham) and time (baseline vs. endpoint) for pain [F(1,28) = 5.596, p = 0.025]. Since our primary outcome, the pain KOOS subscale, showed the most significant difference between verum and sham acupuncture groups, and KOOS pain is the most relevant score for central pain processing, we used the KOOS pain subscale for the following fMRI regression analysis.
Table 1

Demographics and characteristics at baseline, and clinical outcomes before and after longitudinal acupuncture treatments

Mean ± SD

All

High dose

Low dose

Sham

N (F)

30 (13)

10 (2)

10 (7)

10 (4)

Age

58 ± 8

60 ± 9

58 ± 8

54 ± 7

Duration (treated knee, years)*

11 ± 8

10 ± 7

6 ± 6

16 ± 8

Pain

 Baseline

56 ± 14

59 ± 13

53 ± 9

56 ± 19

 Post-treatment

64 ± 14

70 ± 16

66 ± 11

56 ± 12

Symptoms

 Baseline

53 ± 16

57 ± 19

48 ± 11

54 ± 18

 Post-treatment

58 ± 15

60 ± 20

58 ± 15

55 ± 11

Adjusted daily living

 Baseline

64 ± 15

66 ± 12

61 ± 14

64 ± 20

 Post-treatment

72 ± 16

74 ± 18

76 ± 12

65 ± 15

Function in sport

 Baseline

30 ± 23

30 ± 18

31 ± 19

29 ± 31

 Post-treatment

41 ± 26

49 ± 33

48 ± 18

28 ± 21

Quality of life

 Baseline

39 ± 15

42 ± 17

38 ± 14

36 ± 16

 Post-treatment

38 ± 18

44 ± 20

41 ± 16

29 ± 17

Results reported as mean ± SD

Please note that greater score on the KOOS indicates improvement

* Significant main effect of group (high vs. low vs. sham)

N = 9 (one subject missing KOOS sport subscale score)

At the end of the study, we also assessed how well patients’ blinding was maintained throughout the study. Ninety percent (n = 27) of the subjects believed that the needle was inserted into the skin in every session. The 3 subjects who believed that the needle was not inserted were in the real acupuncture (low dose) group.

Independent component analysis

ICA identified three resting state networks (i.e., the rFPN, ECN and somato-sensory network (SMN) as shown in Additional file 2: Figure S2). SMN includes the supplementary motor area, sensorimotor cortex, and secondary somatosensory cortex and insula. The ECN covers several lateral and medial-frontal areas, including the anterior cingulate and paracingulate. The rFPN covers several frontoparietal areas including the bilateral dorsal lateral prefrontal gyri and the inferior and superior parietal lobules [47]. These networks are consistent with findings from previous published results [48].

Functional connectivity: comparing verum to sham group

The comparison in FC changes from pre-acupuncture treatment 1 to pre-acupuncture treatment 6 between the verum and sham groups showed stronger rFPN connectivity with the bilateral rostral anterior cingulate gyrus (rACC)/medial prefrontal cortex (MPFC), parahippocampus, right insula, superior and inferior temporal gyrus, inferior frontal gyrus, left insula/putamen and thalamus in the verum group than in the sham group (Table 2). Figure 1a identifies the rACC/MPFC as a representative brain region that showed significant difference between the verum and sham acupuncture groups. Further analysis showed that Z connectivity values (3 mm sphere around the peak) of the rACC increased gradually across treatment sessions in the verum group but not the sham group. Repeated measures ANOVA showed that there was a highly significant interaction between acupuncture mode (verum vs. sham) and time (across all three treatment sessions) for rACC/MPFC Fisher Z values [F(2,27) = 13.403, p < 0.0001]. After including age and duration of pain as covariates of no interest, the main finding that the verum group showed stronger rFPN connectivity with rACC/MPFC than did the sham group comparing FC changes from treatment 1 to treatment 6 remained valid (Fig. 1b).
Table 2

Regions that show connectivity change with pre-acupuncture resting state network after longitudinal acupuncture treatment (pre-treatment 6—pre-treatment 1) comparing acupuncture to sham group

Network

Group

Area

Z score

Peak coordinate (x, y, z)

Right frontoparietal network

Acu > Sham

rACC/MPFC

4.45

12, 36, −12

Right parahippocampus

3.58

32, 0, −18

Right midinsula/precentralgyrus

3.98

52, −2, 2

Right sup temporal/right inf temporal gyrus

4.15

38, 4, −34

Right inferior frontal/mid frontal gyrus

3.25

50, 32, −4

Left midinsula

4

−24, 4, 6

Left thalamus

3.77

−22, −18, 10

Left claustrum/parahippocampus

3.4

−38, −12, −14

Left sup temporal gyrus

3.46

−50, −6, 0

Sham > Acu

None

Executive control network

Acu > Sham

rACC

3.95

4, 32, 2

Left MPFC

3.4

−16, 56, −8

Left insula

3.72

−46, 2, 2

Left inferior frontal gyrus

3.33

−42, 36, 6

Sham > Acu

None

Somatosensory network

Acu > Sham

None

Sham > Acu

Dorsal ACC

4.74

−8, 18, 38

Bilateral sup frontal/left middle frontal/medial frontal gyrus

3.6

−20, 48, 42

Fig. 1

Comparison of the pre-acupuncture scans in treatments 6 and 1. a Connectivity between the right frontoparietal network and the rACC increases more in the verum group than in the sham group. (X = 12) b Results in a controlled for age and duration of pain. c Executive control network showed stronger connectivity with rACC after treatment in the verum group than in the sham group. (X = 4) d Results in c controlled for age and duration of pain. e Sensory-motor network showed reduced connectivity with dACC after real acupuncture treatment compared with sham group. (X = −8) f Results in e controlled for age and duration of pain

The ECN showed stronger connectivity with rACC/MPFC after verum acupuncture compared to sham acupuncture (Fig. 1c). Further analysis found that in the verum group, but not in the sham group, Z connectivity values of the rACC/MPFC increased gradually across treatment sessions. Repeated measures of ANOVA on Z values (3 mm sphere around the peak) showed that there was a highly significant interaction between acupuncture mode (verum vs. sham) and time (across all three treatment sessions) for the rACC Z value [F(2,27) = 11.85, p < 0.0001]. Other regions that showed significantly greater connectivity with the ECN in the verum group included the left mid insula and inferior frontal gyrus (Table 2). After controlling for age and duration of pain, the results remained the same—stronger connectivity between the ECN and the rACC/MPFC after acupuncture treatments was observed in the verum group relative to the sham group (Fig. 1d).

The SMN showed reduced connectivity with the dorsal ACC (dACC) in the verum group as compared to the sham group (Fig. 1e). Extracted Z values (3 mm sphere around the peak) from the regions showed a gradual decrease in pre-treatment resting state connectivity values across treatment sessions for the verum group, but not for the sham group. Repeated measures ANOVA showed that there was a highly significant interaction between acupuncture mode (verum vs. sham) and time (across all three treatment sessions) for the dACC Z value [F(2,27) = 12.361, p < 0.0001]. Other regions that showed reduced connectivity with the SMN included the bilateral superior frontal, the left middle frontal, and the medial frontal gyrus (Table 2). Controlling for age and duration of pain did not change the results, namely, reduced connectivity between the SMN and the dACC in the verum group compared to the sham group (Fig. 1f).

Association between functional connectivity and clinical outcome

To assess the association between the changes in resting state FC and the corresponding changes in clinical outcome following acupuncture treatment, we applied a regression analysis between the changes of pre-treatment resting state FC maps between treatments 1 and 6 and the corresponding changes in KOOS pain scores (post–pre) for all patients who completed the study. Results showed that the connectivity between the rFPN and the left insula/putamen increased after acupuncture treatment when the change in KOOS pain scores (post–pre) increased (greater KOOS pain score indicates improvement) (r = 0.61, p < 0.001) (Fig. 2a shown in green). Interestingly, this region overlapped with the findings from the group comparison between verum and sham acupuncture, i.e., the FC between the rFPN and left insula significantly increased after treatments in the verum acupuncture group as compared to the sham group (Fig. 2a shown in red). This result further confirmed the important role of the rFPN network in acupuncture modulation.
Fig. 2

a Shown in green: after treatment, the increase in functional connectivity between the rFPN and left insula/putamen positively correlated with the change in KOOS pain score. Shown in red: the verum group showed significant increase in connectivity between the rFPN and left insula/putamen compared to sham group. (X = −30) b Results in a after adjusting for age and duration of pain. The comparison between verum and sham shown in red were at a less conservative threshold of voxel-wise Z > 1.96 and a corrected cluster significance threshold of P < 0.05

After adjusting for age and duration of pain, the rFPN and the left insula/putamen connectivity change after acupuncture treatment remained positively correlated with the change in KOOS pain scores (post–pre) (Fig. 2b shown in green). Moreover, the FC between the rFPN and left insula remained increased after treatments in the verum groups as compared to the sham group with a slightly less stringent threshold (voxel-wise cluster forming threshold of Z > 1.96 and a corrected cluster significance threshold of P < 0.05) (Fig. 2b shown in red).

Other regions that showed a treatment-related increase in connectivity with the rFPN coupled with an increase in KOOS pain score include the temporal gyrus, parahippocampus, and cerebellum.

The ECN showed stronger connectivity with the right middle frontal gyrus coupled with an increase in KOOS pain score. The SMN did not show any significant changes in connectivity that correlated with a change in KOOS pain score (Table 3).
Table 3

Regions in which treatment-related changes in connectivity (pre-treatment 6–pre-treatment 1) were positively correlated with change in clinical KOOS pain (post–pre)

Clinical outcome

ICA networks

Area

Z value

Peak coordinate

KOOS change in pain

Right frontoparietal network

Left anterior insula

4.06

−34, 30, 8

Left posterior insula/left posterior parietal operculum

2.9

−52, −32, 20

Left putamen

3.94

−28, −10, −8

Left superior temporal/left mid temporal gyrus

4.02

−50, −52, 20

Left parahippocampus

3.32

−34, −32, −18

Cerebellum

4.39

18, −48, −38

Executive control network

Right mid frontal gyrus

3.99

38, 20, 32

Sensory-motor network

None

None

None

We observed similar results using the same analysis with a lower threshold (voxel-wise cluster forming threshold of Z > 1.96 and a corrected cluster significance threshold of P < 0.05) for 20 patients who only received verum acupuncture treatments.

Result of task related fMRI data analysis

Verum acupuncture needle stimulation produced significant fMRI signal increases in the left posterior operculum/secondary somatosensory cortex (S2) (x = 66, y = −24, z = 40, 101 voxels) (Fig. 3 shown in green) and the cuneus (x = 3, y = −96, z = 0, 180 voxel) compared to sham treatment. Interestingly, we found that the increased fMRI signal in the parietal operculum overlapped with the result of functional connectivity analysis, i.e., there was a significant association between the KOOS pain score changes and the corresponding rFPN connectivity to the left operculum/S2 changes after longitudinal treatment (r = 0.48, p = 0.007) (Fig. 3 shown in blue).
Fig. 3

Shown in green: during acupuncture needle stimulation, the left posterior parietal operculum showed more activation in the verum group than in the sham group. Shown in blue: after treatment, the change (increased) in connectivity between the right frontoparietal network and left posterior parietal operculum correlates with the change in KOOS pain in patients (Y = −34)

To further explore the brain activity changes evoked by acupuncture needle stimulation, we lowered the threshold to P < 0.005 voxel wise uncorrected with 50 voxels. Our results showed that verum acupuncture needle stimulation also produced significant fMRI signal increase in the left insula and putamen compared to sham acupuncture. Using the same threshold, we also found a significant fMRI signal increase in the bilateral operculum/S2 in both high and low dose groups when we compared brain activity during needle stimulation with brain activity at baseline (in the absence of needle stimulation) (Additional file 3: Figure S3). This finding suggests that there was a similar brain activation pattern in the two verum acupuncture treatment groups.

Discussion

In this study, we investigated the effects of longitudinal acupuncture treatment on brain functional connectivity and clinical outcome measurements in knee OA patients. The results showed that acupuncture treatment can significantly modulate brain functional connectivity and pain as assessed by a KOOS pain subscale in OA patients.

The results from previous studies suggest that acupuncture treatment can significantly modulate FC in healthy subjects [2830, 32, 33]. However, whether acupuncture treatment can modulate the FC in patient populations, which are often characterized by abnormal brain network functioning and requiring longitudinal treatment, remains to be answered. The results from more recent studies suggest that FC changes after medical treatments are significantly associated with clinical outcome changes [4951], implying that the investigation of FC changes can be a useful tool in clinical research. In our study, we found that more FC changes between the rFPN and the left insula correlated with less pain after longitudinal treatment. The verum acupuncture group exhibited a stronger increase in connectivity than did the sham acupuncture group after longitudinal treatment.

Previous studies suggest that the insula/operculum/S2 is involved in the perception of pain [5255] and the frontoparietal network is involved in cognitive control and pain modulation [47, 56]. Downar et al. suggested that the temporal/parietal junction, inferior frontal gyrus, insula and the ACC make up a multimodal cortical network that detects and controls sensory input from the environment [57]. Later, Vincent et al. identified the lateral prefrontal cortex, ACC, insula, and inferior parietal lobule as a frontoparietal control system engaged in cognitive control and decision-making processes [58]. The right lateral frontal cortex, which is part of the rFPN, has exhibited the ability to inhibit pain sensation and modulate pain in previous imaging studies [59, 60].

Anatomically, the posterior parietal cortex, which is also part of the rFPN, is interconnected with the insula, ACC, orbital prefrontal cortex, and parahippocampal cortex and receives projections from the somatosensory cortices [61]. Studies have shown that in knee OA patients, the lateral prefrontal cortex, superior and inferior parietal lobule, rACC, insula, parietal operculum, and limbic cortical areas are involved in altered pain processing [62, 63]. This suggests that the rFPN, rACC, insula, and parietal operculum are the main areas involved in the OA pain circuit. We now show that verum acupuncture modulates activity in these regions more compared to sham.

The ECN located in the anterior of the frontal cortex (including dorsal and medial prefrontal cortex, and anterior cingulate cortex) is involved in integrating information from the external environment with stored internal representations [64], controlling top-down attention during conflict processing of alternative responses [65], and monitoring conflict with subsequent adjustment in performance [6669].

Following longitudinal treatment, both the rFPN and the ECN showed stronger connectivity with the rACC/MPFC in the verum group compared to the sham group. The rACC/MPFC are key regions involved in the descending pain modulatory system [47, 7073]. Studies also showed that the rACC is involved in self-regulation of pain such as placebo analgesia [54, 60, 74]. It also forms a core network with the PAG and the RVM (ACC-PAG-RVM network) for pain modulation even in the absence of a painful stimulus [75].

In a previous study, investigators found that the functional connectivity fluctuations and structural connectivity between the PAG and the rACC/MPFC predicted the mind wandering away from pain [76]. In another study [47], we found that pre-test resting state FC between the rFPN and rACC/MPFC can significantly predict the conditioning placebo cue effect, implying that the intrinsic linkage between the two regions is crucial for self-modulation of pain. Similarly, activity in the rACC/MPFC, insula, and PAG was implicated in pain controllability [77]. In knee OA patients, the rACC co-varies with clinical knee OA pain [62], which suggests that the region is the essential node in pain modulation for knee OA patients.

In a more recent study using the same dataset [78], we found that baseline KOOS pain and sport scores were associated with the PAG FC with MPFC and hippocampus. Specifically, lower connectivity between the PAG and MPFC was associated with worse sport function (lower KOOS sport score), while higher PAG/hippocampus connectivity was associated with worse pain (lower KOOS pain score). Verum acupuncture-induced improvement in pain and sport scores (compared to sham) was associated with the modulation of PAG-MPFC and PAG-hippocampus connectivity—increase and decrease, respectively, suggesting the role of PAG in acupuncture-related improvement. This result is also consistent with a another study [79] using the same data set, in which we found acupuncture can prevent the thinning of cortical thickness at the posterior MPFC (a region in ECN) in knee OA patients, and the resting state FC between the posterior MPFC and the key descending pain modulatory system (rACC and PAG) is significantly stronger in real acupuncture as compared with sham after longitudinal treatment.

In light of the above findings and previous studies on the mechanisms involved in acupuncture analgesia [1921] demonstrating the involvement of the pain descending modulatory system, we speculate that the enhanced connectivity between the rFPN and ECN networks and the rACC/MPFC may suggest that verum acupuncture decreases patients’ pain experience by enhancing the descending pain modulatory system.

In the present study, we found that following longitudinal treatment, FC between the sensorimotor network and the dorsal ACC exhibited a greater decline in the verum group than in the sham group. It is well known that the dorsal ACC is associated with the affective component of pain [55, 60, 8082] and can be modulated by acupuncture needle stimulation [27]. We speculate that the decrease in the connectivity between these two regions may suggest that acupuncture can modulate (relieve) affective pain experience [83, 84].

In recent years, many investigators have explored acupuncture needle stimulation-evoked brain activity changes using fMRI [23, 24]. A wide range of brain networks are involved in acupuncture stimulation and each region may respond differently to needle stimulation [85]. In this study, we found that during acupuncture needle stimulation, there is stronger activation in the left parietal operculum in verum acupuncture as compared to sham acupuncture, which is consistent with findings from previous studies [23, 24, 26, 86]. A study by Pariente and colleagues [87] compared verum, sham (Streitberger needle) and open placebo acupuncture and found that verum acupuncture activated the insula, DLPFC, rACC, and midbrain compared to open sham; however, only the insula was more activated by verum acupuncture when accounting for expectation (comparison with the Streiberger needle stimulation). Our results are consistent with above findings.

Additionally, we found a positive association in the functional connectivity changes between the rFPN and the left parietal operculum (post–pre) and the corresponding changes in clinical KOOS pain score. These results suggest that lower levels of pain correlate with stronger connectivity between the rFPN and the left posterior parietal operculum after longitudinal acupuncture treatment. Investigators have identified the parietal operculum as one of the main regions related to OA pain [63] and pain perception [88]. Our results suggest that the parietal operculum may be a key region in acupuncture treatment that serves to enhance connectivity between the rFPN and the insula and rACC/MPFC, resulting in increased communication between these regions and restoration of the descending modulatory pathway.

Summing up, it seems that the insula/operculum activity appeared in all analyses: (1) fMRI signal increase during needling stimulation (verum vs sham); (2) resting state FC group differences, and (3) association between the resting state FC changes and the corresponding KOOS pain changes in all subjects. A review of the insula function [89] suggested that it represents a hub for (1) bottom-up detection of salient stimuli, (2) switching between large-scale attention and memory networks, (3) interacting with the anterior and posterior insula to modulate autonomic reactivity, and 4) coupling with the anterior cingulate to facilitate rapid access to the motor system. We observed brain activity and connectivity patterns consistent with all of these functions during verum acupuncture compared to sham. Therefore, in line with this model, we speculate that verum acupuncture needle stimulation activates the operculum/insula, which triggers the rFPN and ECN. The two networks further activate the descending pain modulatory network through the rACC/MPFC to achieve the modulation effect in the descending pain network (please note that although no differences in the PAG were directly observed in the current study, the intrinsic connectivity of the PAG and rACC/MPFC has been demonstrated in our previous study [75]). In parallel, the coupling between the insula/operculum and dACC modulates the sensory-motor network, interfering with the sensory motor changes evoked by knee pain to relieve pain experience (Fig. 4). Future studies are needed to elucidate the details of the modulation process.
Fig. 4

Hypothetical schematic illustration of the resting state functional connectivity modulated by acupuncture and its relevance to pain relief. Note that this diagram does not exhaustively describe all networks and brain regions involved in acupuncture analgesia but only summarizes the mechanisms likely relevant for the current study. Red indicates increased activation/connectivity; blue indicates decreased activation/connectivity. Verum acupuncture needle stimulation activated the operculum/insula, as suggested by the functional MRI analysis. The insula/operculum processes information about both the sensory component of pain (posterior insula), as well as cognitive-emotional aspects of pain (anterior insula) [101]. In its capacity as a nociceptive salience detection, affective and pain decision-making hub [102], it might (1) increase functional connectivity between the attention, cognitive control and appraisal networks (rFPN, ECN) and the rACC/MPFC, which is a key brain region for attention and descending pain modulation with a direct connection to the PAG [75], further inhibiting noxious input; (2) decrease functional connectivity between the SMN and dACC, representing reduced interaction between the sensory and affective components of pain processing, providing further relief from pain

One attribute of this study is that the needles were repeatedly manipulated. This differs significantly from the previous clinical trials [13, 14], which found no significant difference between real and sham acupuncture (in these trials, the needle stimulation is almost not applied). Based on Traditional Chinese Acupuncture theory, deqi (the sensations evoked by acupuncture needle stimulation) is key for the clinical achievement of acupuncture treatment [46]; previous studies [90] also suggested that strong or moderate stimulation is crucial for acupuncture analgesia. This may explain the group differences in KOOS pain observed with a relatively small sample size. Nevertheless, we would like to emphasize that our result is very preliminary and further study with a larger sample size is needed to validate this finding.

In this study, we also attempted to explore the effect of acupoint number. With a small sample size, our preliminary result suggested that if the total stimulation amount remains the same, the treatment effect and brain activation pattern is also similar in knee OA patients. Further study is needed to validate the findings with large sample size in different patient populations.

Several limitations of this study should be mentioned. First, the sample size is relatively small. However, we made several steps to increase the power in the following ways. On the one hand, by combining the two acupuncture groups, we increased the power of the verum group, creating a 2:1 ratio between the verum and sham treatment conditions that provided additional information on patients receiving active treatment [91]. On the other hand, due to the longitudinal design of the study, we acquired fMRI scans at three time points, which may have also increased the power of this study. All reported main brain imaging results survived the stringent multiple comparison correction.

Another limitation of this study was the potentially confounding effect of factors such as patient age and duration of chronic pain experienced by each patient. To address this concern, we repeated the analysis including age and duration of pain as covariates of no interest and re-analyzed the data. After this analysis, most of the key findings remained significant. Thus, our main results do not change after adjusting for these confounding factors. Furthermore, when KOOS pain scores were correlated with age and duration using multiple regression, we did not find any significant correlation. Further studies with larger sample sizes are needed to validate the findings of the study.

Thirdly, patients with knee OA very commonly report pain in other regions, including the neck and low back, and the contralateral knee and hips, which may also have influenced brain FC. However, we can assume these conditions were controlled for after randomization. In addition, the comparisons in this study were applied between pre- and post-treatment differences. Thus, we expect that these potential confounding factors should not influence the conclusions of this study.

Fourth, there are two patients who dropped in each real acupuncture treatment group. Due to a small sample size, this could potentially have influenced the conclusion of our study. Nevertheless, as we pointed out above, of the four patients who withdrew, one subject was unable to adhere to study requirements in the scanner and three had scheduling conflicts (the fact that the fMRI treatment session can only take place in certain time intervals significantly limited the flexibility of this study). None of these reasons, to the best of our knowledge, was related to the clinical outcome of treatments. We therefore believe it did not influence the validity of the study.

Finally, all subjects agreed not to use their regular medication during the acupuncture treatment period, except for painkillers as needed. However, we cannot completely exclude the potential influence of other treatments.

Conclusion

Our findings suggest that verum acupuncture may enhance connectivity in the descending pain modulation pathway through several networks, including the rFPN and ECN. A better understanding of the association between FC and clinical outcomes and how treatment can modulate FC may ultimately lead to the acceptance of acupuncture in mainstream medicine, and facilitate the development of mechanism-based therapies for chronic pain.

Methods

We briefly describe the experimental procedures below. Please also see previously published studies [45, 78, 79] for more details on the experimental procedure. The data has been used in previous studies to investigate the reliability of deqi sensation evoked by acupuncture needle stimulation [45], cortical thickness changes after longitudinal real and sham acupuncture treatments [79], and resting state FC changes of the PAG [78]. In this study, we used ICA to investigate the resting state FC changes before and after longitudinal acupuncture treatments, as well as the brain response during acupuncture needle stimulation. These results have not been reported before.

Subjects

The Institutional Review Board at Massachusetts General Hospital approved all study procedures. Subjects who enrolled in this study provided written informed consent before beginning any study procedures. We debriefed all subjects at the end of the study.

Patient recruitment and inclusion criteria

We recruited acupuncture naïve patients aged 40–70 with a diagnosis of chronic osteoarthritis (OA) in the right and/or left knee from Massachusetts General Hospital (MGH) and Brigham and Women’s Hospital (BWH).

Subjects exhibited Grade 2 or Grade 3 OA as defined by the Kellgren-Lawrence Scale used for radiographically grading knee OA. We excluded subjects who: (1) had undergone any interventional procedures for knee pain within 6 months prior to enrolling in the study, (2) intended to undergo surgery during the time of involvement in the study, (3) experienced knee pain due to causes other than OA, such as inflammation or malignancy, (4) had received a diagnosis of rheumatoid arthritis or other leg-related pain disorders, (5) were taking opioids, benzodiazepines, or other medications that may influence brain imaging outcome, or (6) presented MRI contraindications. All enrolled OA patients had an endogenous average pain intensity rating of >2 on the Brief Pain Inventory (BPI) scale at the first visit.

Experimental design

All subjects, regardless of the presence of bilateral or unilateral pain, received treatment on the knee that presented them with the most severe pain. We stratified subjects by knee and randomized them into one of three groups: high dose real acupuncture (6 acupoints), low dose real acupuncture (2 acupoints), or sham acupuncture (6 non-acupoints with Streitberger placebo needles [92]) (Additional file 4: Figure S4). We used the permuted block randomization stratified by the acupuncture-treated knee. The randomization number was provided by a biostatistician and only the acupuncturist was not blinded to the study randomization before the end of the study. The randomization was applied after the QST measurement and psychometric assessment in session 1. This time was chosen because the different acupuncture exposure took place immediately afterwards.

In this study, we were interested in whether the number of needles (acu points) can influence the treatment effect. Therefore, we had two groups of real acupuncture treatments with different numbers of acupoints, but an identical stimulation paradigm.

After an initial screening session, each subject engaged in a total of 6 acupuncture treatment sessions in 1 month (twice per week for the first 2 weeks, once per week for the last 2 weeks). Treatments 1, 3 and 6 occurred during a scan session with the patient lying in a 3-Tesla MRI scanner while we acquired fMRI data. All other treatments occurred in our behavioral testing room.

Acupuncture administration

Each acupuncture treatment session for subjects in both the verum and sham acupuncture groups lasted about 25 min. The two acupoints selected for the low dose acupuncture group [ST35 and Xi yian (extra point)] were located near the knee and each acupoint was stimulated a total of twelve times. The high dose acupuncture group received treatment at four additional acupoints (GB34, SP9, GB39 and SP6) (see Additional file 1: Figure S1b) and each acupoint was stimulated 4 times. All these points are well documented for the treatment of knee pain [8, 10, 11].

For subjects receiving sham acupuncture, the acupuncturist applied placebo needles [54, 9294] at six non-acupoints using a paradigm identical to that of the high dose verum needle administration. Rather than penetrating the skin like a verum needle, the sham needle retracts up the handle shaft when the acupuncturist presses it against the skin. Sham point 1 was located 1.5 cun posterior and inferior to GB4, sham points 2–3 were located 1.5 and 3 cun inferior to sham point 1, sham point 4 was located 1 cun posterior to the midpoint of K9 and K10, and sham 5–6 were located 1.5 cun inferior and superior to the sham point 4, respectively. All sham points were located on the lower leg where no meridians pass (Additional file 1: Figure S1b).

For all treatments, the acupuncturist placed a small plastic ring over the point and secured the ring with a thin strip of sterile plastic tape. This ensured patient blindness to the actual site of needle insertion and thus blindness to whether the treatment was verum or sham. During acupuncture, the acupuncturist stimulated one point at a time in a predetermined order, each for 10 s with 30-s breaks between each acupoint. (Additional file 1: Figure S1a). The total treatment and stimulation time for all three groups was the same. We randomized the specific starting acupoint across patients, but held it constant throughout all sessions for each individual patient. For consistency, we kept leg position, acupoint location, and needling parameters (1–2 cm depth, approximately 120 rotations/min, and moderate deqi sensations on a 0–10 scale) constant across groups. In fMRI session, there were two acupuncture scans per treatment, with each lasted for 9 min, followed by Massachusetts General Hospital Acupuncture Sensation Scale (MASS) rating.

As a treatment quality control, we measured the deqi sensation evoked by two separate acupuncture treatments (one midway through the treatment and one at the end of the treatment) using the MASS [45, 46].

Clinical outcomes

Knee injury and Osteoarthritis Outcome Score (KOOS)

We measured clinical outcomes using the KOOS [44], which is comprised of 5 subscales: (1) pain, (2) other symptoms, (3) function in daily living (ADL), (4) function in sport and recreation, and (5) knee-related quality of life (QOL). Each subscale allows for calculation of a normalized score, with 0 denoting the most extreme symptoms/pain and 100 denoting no symptoms/pain [44]. Based on previous studies [8], we selected the KOOS pain subscale and function in daily living subscale as our primary outcome measures. We used all other subscale measurements as secondary outcomes. Trained research assistants, blinded to treatment mode, administered the KOOS to all patients at baseline and again after the final treatment.

fMRI data acquisition

Each subject participated in three identical fMRI scanning sessions. Scan 1 (treatment 1) and scan 2 (treatment 3) were separated by approximately 7 days; scan 2 and scan 3 (treatment 6) were separated by approximately 14 days. We used a 3-axis gradient head coil in a 3-Tesla Siemens MRI system equipped for echo planar imaging (EPI). Each fMRI session included an anisotropic magnetization prepared rapid gradient-echo (MPRAGE) structural sequence followed by a 6-min resting state fMRI. The MPRAGE scanning parameters included TR of 2200 ms, echo time of 9.8 ms, flip angle of 7, field-of-view of 230 mm2, slice thickness of 1.2 mm. For the resting state, the scan acquisition included 47 slices with a thickness of 3 mm, a TR of 3000 ms, a TE of 30 ms, flip angle of 85°, field of view of 216 mm2 and a 3 × 3-mm in-plane spatial resolution. Following the first resting state fMRI, we conducted two functional scans during the 25-min acupuncture administration and a second resting state scan following acupuncture administration (Additional file 1: Figure S1a).

Data analysis

Independent component analysis for resting state fMRI data

We analyzed pre- and post-acupuncture resting state data for treatments 1, 3 and 6 using independent component analysis in FMRIB Software Library (FSL) [95], following the same processing steps as described in a previous study [47, 48]. We first applied a band pass filter between 0.01 and 0.1 Hz to the functional time series corrected for motion using MCFLIRT and slice timing, skull stripped using the Brain Extraction Tool (BET), and smoothed (full width at half maximum = 5 mm). Then we registered the data to their respective skull stripped anatomical volume and further registered to the MNI152 template using linear affine transformations with 12° of freedom. We then concatenated the functional data into 4D data and performed a probabilistic independent component analysis using Multivariate Exploratory Linear Optimized Decomposition into Independent Components (MELODIC) [96] on the data set to identify 20 resting state networks. We used an algorithm to search for similarities between our group-level networks and the template networks derived from 1414 healthy subjects to identify the corresponding network for our results [48].

We chose the sensory-motor network (SMN), executive control network (ECN), and right frontoparietal network (rFPN) for further analyses. Previous ICA studies have reported that these networks have been reliably observed and are believed to be associated with pain process and modulation [47, 56] as well as development of chronic pain [97, 98].

To perform group-level analyses of the association between subjects’ responses to acupuncture and resting state networks derived from ICA, we used a dual-regression technique [48]. In brief, we used the three previously defined networks as spatial regressors in a general linear model (GLM) to extract temporal dynamics associated with each spatial map. The resulting time courses served as temporal regressors in a GLM to generate subject-specific-maps of the whole brain for each subject.

Finally, group analyses were performed using the whole-brain subject-specific network maps from the second GLM. The results represent the strength of FC for each voxel within each of the chosen networks.

To explore the association between the clinical outcomes and FC changes, we performed regression analyses using the network connectivity “change maps” between pre-acupuncture treatment 1 and pre-acupuncture treatment 6 resting state scans and change in clinical outcome (KOOS pain score post–pre). To remove any baseline effect, we input all changes in KOOS pain scores as percent change [(postpre/pre) × 100 %].

To investigate the modulation effect of acupuncture treatment, we compared the changes in pre-acupuncture treatment 6 and pre-acupuncture treatment 1 in the verum group to that of the sham group. We used pre-acupuncture resting state fMRI to avoid the residual effect of acupuncture needle stimulation [99, 100]. Due to the small sample size of each group (10 patients), we decided to combine the two verum acupuncture groups (low and high dose groups) because (1) the total treatment time, needle stimulation time, and reported sensations were the same in both verum groups; (2) we found no clinical outcome differences between the two verum acupuncture groups and both verum acupuncture groups showed a trend towards improvement when compared to sham treatment, and (3) high and low dose acupuncture needle stimulations evoked similar activation patterns (Additional file 3: Figure S3); (4) reported sensations as measured by MASS were similar in both real acupuncture treatment groups, and both showed significant differences from the sham acupuncture group; (5) acupuncture prescription for Knee OA varies across different studies; combining two protocols may better represent clinical acupuncture in practice; and (6) combining two groups gave us more power to test our hypothesis. We performed all resting state analyses with a voxel-wise cluster forming threshold of Z > 2.3 and a corrected cluster significance threshold of P < 0.05.

To test the effect of acupuncture treatment across all three sessions, we extracted the peak Z values of each significant network from previous steps for both verum and sham group and performed repeated measure ANOVA with each treatment session as one time point using SPSS 18.0 Software (SPSS Inc., Chicago, IL, USA).

To control for the potentially confounding factors including age and duration of knee OA, we also repeated the above analysis including age and duration of pain as covariates of no interest.

Task-related fMRI data analysis

To explore how acupuncture needle stimulation affected brain activity changes, we performed data analysis using Statistical Parametric Mapping (SPM8, Wellcome Trust Centre for Neuroimaging, London, UK). During the preprocessing, we realigned all functional volumes, spatially normalized, and smoothed using an 8-mm full-width at half-maximum Gaussian kernel.

We used a general Linear Model (GLM) for first level analysis. Each individual design matrix contained the regressor for the main effect of acupuncture, combining acupuncture sessions 1 and 2 for each treatment session. To evaluate the effect of verum greater than sham acupuncture, we performed second level analyses using a full factorial design with two factors. The first factor had two levels (high and low dose verum acupuncture, and sham acupuncture) and the second factor had three levels (treatment 1, treatment 3 and treatment 6). By pooling all 3 treatment sessions together for each group, we estimated the effect of verum acupuncture as compared with sham acupuncture in the second level analyses. To further explore the brain activation evoked by needle stimulation during high dose and low dose acupuncture, we also calculated the fMRI signal change for each group compared to baseline separately. A threshold of voxel wise p < 0.001 uncorrected and p < 0.05 corrected (family-wise error, FWE) at cluster level was used for data analysis.

Notes

Declarations

Authors’ contributions

Experimental design: JK, JRM, AJW, RRE, RLG; data collection: RBS, DMS, JK; data analysis: XYC, JAH, HYW, MV; manuscript preparation/revision: JK, XYC, AJW, RRE, RLG, NE. All authors read and approved the final manuscript.

Acknowledgements

This work was supported by R21AT004497 (NIH/NCCIH), R03AT218317 (NIH/NIDA), R01AT006364 (NIH/NCCIH), R01AT008563 (NIH/NCCIH) to Jian Kong, R01AT005280 (NCCIH) to Randy Gollub, and P01 AT006663 to Bruce Rosen.

Competing interests

All authors declare that they have no competing interests.

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Authors’ Affiliations

(1)
Department of Psychiatry, Massachusetts General Hospital
(2)
Mass General Orthopaedics Sports Performance Center, Massachusetts General Hospital
(3)
MGH/MIT/HMS Athinoula A. Martinos Center for Biomedical Imaging
(4)
Department of Radiology, Massachusetts General Hospital
(5)
Departments of Anesthesiology, Massachusetts General Hospital
(6)
Departments of Anesthesiology, Perioperative, and Pain Medicine, and Psychiatry, Brigham and Women’s Hospital and Harvard Medical School
(7)
Departments of Anesthesiology and Psychiatry, University of Pittsburgh School of Medicine

References

  1. Ho-Pham LT, Lai TQ, Mai LD, Doan MC, Pham HN, Nguyen TV. Prevalence of radiographic osteoarthritis of the knee and its relationship to self-reported pain. PLoS One. 2014;9(4):e94563. doi:https://doi.org/10.1371/journal.pone.0094563.PubMed CentralView ArticlePubMedGoogle Scholar
  2. Nguyen US, Zhang Y, Zhu Y, Niu J, Zhang B, Felson DT. Increasing prevalence of knee pain and symptomatic knee osteoarthritis: survey and cohort data. Ann Intern Med. 2011;155(11):725–32. doi:https://doi.org/10.7326/0003-4819-155-11-201112060-00004.PubMed CentralView ArticlePubMedGoogle Scholar
  3. Hochman JR, Davis AM, Elkayam J, Gagliese L, Hawker GA. Neuropathic pain symptoms on the modified painDETECT correlate with signs of central sensitization in knee osteoarthritis. Osteoarthr Cartil. 2013;21(9):1236–42. doi:https://doi.org/10.1016/j.joca.2013.06.023.View ArticlePubMedGoogle Scholar
  4. Arendt-Nielsen L, Nie H, Laursen MB, Laursen BS, Madeleine P, Simonsen OH, et al. Sensitization in patients with painful knee osteoarthritis. Pain. 2010;149(3):573–81.View ArticlePubMedGoogle Scholar
  5. Finan PH, Buenaver LF, Bounds SC, Hussain S, Park RJ, Haque UJ, et al. Discordance between pain and radiographic severity in knee osteoarthritis: findings from quantitative sensory testing of central sensitization. Arthritis Rheum. 2013;65(2):363–72. doi:https://doi.org/10.1002/art.34646.View ArticlePubMedGoogle Scholar
  6. Gwilym SE, Keltner JR, Warnaby CE, Carr AJ, Chizh B, Chessell I, et al. Psychophysical and functional imaging evidence supporting the presence of central sensitization in a cohort of osteoarthritis patients. Arthritis Rheum. 2009;61(9):1226–34.View ArticlePubMedGoogle Scholar
  7. Gay MC, Philippot P, Luminet O. Differential effectiveness of psychological interventions for reducing osteoarthritis pain: a comparison of Erikson [correction of Erickson] hypnosis and Jacobson relaxation. Eur J Pain. 2002;6(1):1–16. doi:https://doi.org/10.1053/eujp.2001.0263.View ArticlePubMedGoogle Scholar
  8. Berman BM, Lao L, Langenberg P, Lee WL, Gilpin AM, Hochberg MC. Effectiveness of acupuncture as adjunctive therapy in osteoarthritis of the knee: a randomized, controlled trial. Ann Intern Med. 2004;141(12):901–10.View ArticlePubMedGoogle Scholar
  9. Felson DT, Lawrence RC, Hochberg MC, McAlindon T, Dieppe PA, Minor MA, et al. Osteoarthritis: new insights. Part 2: treatment approaches. Ann Intern Med. 2000;133(9):726–37.View ArticlePubMedGoogle Scholar
  10. Witt C, Brinkhaus B, Jena S, Linde K, Streng A, Wagenpfeil S, et al. Acupuncture in patients with osteoarthritis of the knee: a randomised trial. Lancet. 2005;366(9480):136–43.View ArticlePubMedGoogle Scholar
  11. Scharf HP, Mansmann U, Streitberger K, Witte S, Kramer J, Maier C, et al. Acupuncture and knee osteoarthritis: a three-armed randomized trial. Ann Intern Med. 2006;145(1):12–20.View ArticlePubMedGoogle Scholar
  12. Mavrommatis CI, Argyra E, Vadalouka A, Vasilakos DG. Acupuncture as an adjunctive therapy to pharmacological treatment in patients with chronic pain due to osteoarthritis of the knee: a 3-armed, randomized, placebo-controlled trial. Pain. 2012;153(8):1720–6.View ArticlePubMedGoogle Scholar
  13. Hinman RS, McCrory P, Pirotta M, Relf I, Forbes A, Crossley K, et al. Acupuncture for chronic knee pain a randomized clinical trial. JAMA. 2014;312(13):1313–22.View ArticlePubMedGoogle Scholar
  14. Suarez-Almazor M, Looney C, Kiu Y, Cox V, Pietz K, Marcus D, et al. A randomized controlled trial of acupuncture for osteoarthritis of the knee: effects of patient-practitioner communication. Arthritis Care Res. 2010;62:1229–36.View ArticleGoogle Scholar
  15. MacPherson H, Thomas K, Walters S, Fitter M. The York acupuncture safety study: prospective survey of 34 000 treatments by traditional acupuncturists. BMJ. 2001;323(7311):486–7.PubMed CentralView ArticlePubMedGoogle Scholar
  16. White A, Hayhoe S, Hart A, Ernst E. Adverse events following acupuncture: prospective survey of 32 000 consultations with doctors and physiotherapists. BMJ. 2001;323(7311):485–6.PubMed CentralView ArticlePubMedGoogle Scholar
  17. Melchart D, Weidenhammer W, Streng A, Reitmayr S, Hoppe A, Ernst E, et al. Prospective investigation of adverse effects of acupuncture in 97 733 patients. Arch Intern Med. 2004;164(1):104–5. doi:https://doi.org/10.1001/archinte.164.1.104.View ArticlePubMedGoogle Scholar
  18. Ochi JW. Acupuncture instead of codeine for tonsillectomy pain in children. Int J Pediatr Otorhinolaryngol. 2013;77(12):2058–62. doi:https://doi.org/10.1016/j.ijporl.2013.10.008.View ArticlePubMedGoogle Scholar
  19. Han JS. Acupuncture analgesia: areas of consensus and controversy. Pain. 2011;152(3 Suppl):S41–8.View ArticlePubMedGoogle Scholar
  20. Hans JS. Physiology of acupuncture: review of thirty years of research. J Altern Complement Med. 1997;(Suppl 1):S101–8.Google Scholar
  21. Zhao ZQ. Neural mechanism underlying acupuncture analgesia. Prog Neurobiol. 2008;85(4):355–75.View ArticlePubMedGoogle Scholar
  22. Dougherty DD, Kong J, Webb M, Bonab AA, Fischman AJ, Gollub RL. A combined [11C]diprenorphine PET study and fMRI study of acupuncture analgesia. Behav Brain Res. 2008;193(1):63–8 (PMCID: PMC2538486).Google Scholar
  23. Chae Y, Chang DS, Lee SH, Jung WM, Lee IS, Jackson S, et al. Inserting needles into the body: a meta-analysis of brain activity associated with acupuncture needle stimulation. J Pain. 2013;14(3):215–22.View ArticlePubMedGoogle Scholar
  24. Huang W, Pach D, Napadow V, Park K, Long X, Neumann J, et al. Characterizing acupuncture stimuli using brain imaging with FMRI—a systematic review and meta-analysis of the literature. PLoS One. 2012;7(4):e32960.PubMed CentralView ArticlePubMedGoogle Scholar
  25. Kong J, Ma L, Gollub RL, Wei J, Yang X, Li D, et al. A pilot study of functional magnetic resonance imaging of the brain during manual and electroacupuncture stimulation of acupuncture point (LI-4 Hegu) in normal subjects reveals differential brain activation between methods. J Altern Complement Med. 2002;8(4):411–9.View ArticlePubMedGoogle Scholar
  26. Kong J, Gollub RL, Webb JM, Kong JT, Vangel MG, Kwong K. Test-retest study of fMRI signal change evoked by electroacupuncture stimulation. Neuroimage. 2007;34(3):1171–81 (PMCID: PMC1994822).Google Scholar
  27. Leung A, Zhao Y, Shukla S. The effect of acupuncture needle combination on central pain processing—an fMRI study. Mol Pain. 2014;10:23. doi:https://doi.org/10.1186/1744-8069-10-23.PubMed CentralView ArticlePubMedGoogle Scholar
  28. Zyloney CE, Jensen K, Polich G, Loiotile RE, Cheetham A, LaViolette PS, et al. Imaging the functional connectivity of the Periaqueductal Gray during genuine and sham electroacupuncture treatment. Mol Pain. 2010;6:80.PubMed CentralView ArticlePubMedGoogle Scholar
  29. Dhond RP, Yeh C, Park K, Kettner N, Napadow V. Acupuncture modulates resting state connectivity in default and sensorimotor brain networks. Pain. 2008;136(3):407–18.PubMed CentralView ArticlePubMedGoogle Scholar
  30. Liu P, Qin W, Zhang Y, Tian J, Bai L, Zhou G, et al. Combining spatial and temporal information to explore function-guide action of acupuncture using fMRI. J Magn Reson Imaging. 2009;30(1):41–6.View ArticlePubMedGoogle Scholar
  31. Bai L, Qin W, Tian J, Dong M, Pan X, Chen P, et al. Acupuncture modulates spontaneous activities in the anticorrelated resting brain networks. Brain Res. 2009;1279:37–49. doi:https://doi.org/10.1016/j.brainres.2009.04.056.View ArticlePubMedGoogle Scholar
  32. Qin W, Tian J, Bai L, Pan X, Yang L, Chen P, et al. FMRI connectivity analysis of acupuncture effects on an amygdala-associated brain network. Mol Pain. 2008;4:55.PubMed CentralView ArticlePubMedGoogle Scholar
  33. Liu B, Chen J, Wang J, Liu X, Duan X, Shang X, et al. Altered small-world efficiency of brain functional networks in acupuncture at ST36: a functional MRI study. PLoS One. 2012;7(6):e39342. doi:https://doi.org/10.1371/journal.pone.0039342.PubMed CentralView ArticlePubMedGoogle Scholar
  34. Sun R, Yang Y, Li Z, Li Y, Cheng S, Zeng F. Connectomics: a New Direction in Research to Understand the Mechanism of Acupuncture. Evid Based Complement Alternat Med. 2014;2014:568429. doi:https://doi.org/10.1155/2014/568429.PubMed CentralPubMedGoogle Scholar
  35. Zhang Y, Qin W, Liu P, Tian J, Liang J, von Deneen KM, et al. An fMRI study of acupuncture using independent component analysis. Neurosci Lett. 2009;449(1):6–9. doi:https://doi.org/10.1016/j.neulet.2008.10.071.View ArticlePubMedGoogle Scholar
  36. Apkarian AV, Baliki MN, Geha PY. Towards a theory of chronic pain. Prog Neurobiol. 2009;87(2):81–97.PubMed CentralView ArticlePubMedGoogle Scholar
  37. Bushnell MC. How neuroimaging studies have challenged us to rethink: is chronic pain a disease? J Pain. 2009;10(11):1113–20.View ArticlePubMedGoogle Scholar
  38. Davis KD, Moayedi M. Central mechanisms of pain revealed through functional and structural MRI. J Neuroimmune Pharmacol. 2013;8(3):518–34.View ArticlePubMedGoogle Scholar
  39. 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.doi:https://doi.org/10.1038/nn.3153.
  40. Loggia ML, Kim J, Gollub RL, Vangel MG, Kirsch I, Kong J, et al. Default mode network connectivity encodes clinical pain: an arterial spin labeling study. Pain. 2013;154:24–33.PubMed CentralView ArticlePubMedGoogle Scholar
  41. Kong J, Spaech R, Wey HY, Cheetham A, Cook A, Jensen K et al. S1 is associated with chronic low back pain: a functional and structural MRI study. Molecular pain. 2013;9(43):doi:https://doi.org/10.1186/744-8069-9-43.
  42. Niesters M, Khalili-Mahani N, Martini C, Aarts L, van Gerven J, van Buchem MA, et al. Effect of subanesthetic ketamine on intrinsic functional brain connectivity: a placebo-controlled functional magnetic resonance imaging study in healthy male volunteers. Anesthesiology. 2012;117(4):868–77. doi:https://doi.org/10.1097/ALN.0b013e31826a0db3.View ArticlePubMedGoogle Scholar
  43. Harris RE, Napadow V, Huggins JP, Pauer L, Kim J, Hampson J, et al. Pregabalin rectifies aberrant brain chemistry, connectivity, and functional response in chronic pain patients. Anesthesiology. 2013;119(6):1453–64. doi:https://doi.org/10.1097/ALN.0000000000000017.View ArticlePubMedGoogle Scholar
  44. Roos EM, Toksvig-Larsen S. Knee injury and Osteoarthritis Outcome Score (KOOS) – validation and comparison to the WOMAC in total knee replacement. Health Quality Life Outcomes. 2003;1(17):1–10.Google Scholar
  45. Spaeth RB, Camhi S, Hashmi JA, Vangel M, Wasan AD, Edwards RR, et al. A longitudinal study of the reliability of acupuncture deqi sensations in knee osteoarthritis. Evid Based Complement Alternat Med. 2013;2013:204259.PubMed CentralPubMedGoogle Scholar
  46. Kong J, Gollub R, Huang T, Polich G, Napadow V, Hui K, et al. Acupuncture de qi, from qualitative history to quantitative measurement. J Altern Complement Med. 2007;13(10):1059–70.View ArticlePubMedGoogle Scholar
  47. Kong J, Jensen K, Loiotile R, Cheetham A, Wey HY, Tan T, et al. Functional connectivity of frontoparietal network predicts cognitive modulation of pain. Pain. 2013;154(3):459–67.PubMed CentralView ArticlePubMedGoogle Scholar
  48. Biswal BB, Mennes M, Zuo XN, Gohel S, Kelly C, Smith SM, et al. Toward discovery science of human brain function. Proc Natl Acad Sci USA. 2010;107(10):4734–9.PubMed CentralView ArticlePubMedGoogle Scholar
  49. Goldin PR, Ziv M, Jazaieri H, Hahn K, Heimberg R, Gross JJ. Impact of cognitive behavioral therapy for social anxiety disorder on the neural dynamics of cognitive reappraisal of negative self-beliefs: randomized clinical trial. JAMA Psychiatry. 2013;70(10):1048–56. doi:https://doi.org/10.1001/jamapsychiatry.2013.234.PubMed CentralView ArticlePubMedGoogle Scholar
  50. Pavuluri MN, Ellis JA, Wegbreit E, Passarotti AM, Stevens MC. Pharmacotherapy impacts functional connectivity among affective circuits during response inhibition in pediatric mania. Behav Brain Res. 2012;226(2):493–503. doi:https://doi.org/10.1016/j.bbr.2011.10.003.PubMed CentralView ArticlePubMedGoogle Scholar
  51. Fang J, Rong P, Hong Y, Fan Y, Liu J, Wang H, et al. Transcutaneous vagus nerve stimulation modulates default mode network in major depressive disorder. Biol Psychiatry. 2015;. doi:https://doi.org/10.1016/j.biopsych.2015.03.025.Google Scholar
  52. Apkarian AV, Bushnell MC, Treede RD, Zubieta JK. Human brain mechanisms of pain perception and regulation in health and disease. Eur J Pain. 2005;9(4):463–84.View ArticlePubMedGoogle Scholar
  53. Tracey I, Mantyh PW. The cerebral signature for pain perception and its modulation. Neuron. 2007;55(3):377–91.View ArticlePubMedGoogle Scholar
  54. Kong J, Gollub RL, Rosman IS, Webb JM, Vangel MG, Kirsch I, et al. Brain activity associated with expectancy-enhanced placebo analgesia as measured by functional magnetic resonance imaging. J Neurosci. 2006;26(2):381–8.View ArticlePubMedGoogle Scholar
  55. Kong J, Loggia ML, Zyloney C, Tu P, Laviolette P, Gollub RL. Exploring the brain in pain: activations, deactivations and their relation. Pain. 2010;148:257–67.PubMed CentralView ArticlePubMedGoogle Scholar
  56. Smith SM, Fox PT, Miller KL, Glahn DC, Fox PM, Mackay CE, et al. Correspondence of the brain’s functional architecture during activation and rest. Proc Natl Acad Sci USA. 2009;106(31):13040–5.PubMed CentralView ArticlePubMedGoogle Scholar
  57. Downar J, Crawley AP, Mikulis DJ, Davis KD. A multimodal cortical network for the detection of changes in the sensory environment. Nat Neurosci. 2000;3(3):277–83. doi:https://doi.org/10.1038/72991.View ArticlePubMedGoogle Scholar
  58. Vincent JL, Kahn I, Snyder AZ, Raichle ME, Buckner RL. Evidence for a frontoparietal control system revealed by intrinsic functional connectivity. J Neurophysiol. 2008;100(6):3328–42.PubMed CentralView ArticlePubMedGoogle Scholar
  59. Mayer EA, Berman S, Suyenobu B, Labus J, Mandelkern MA, Naliboff BD, et al. Differences in brain responses to visceral pain between patients with irritable bowel syndrome and ulcerative colitis. Pain. 2005;115(3):398–409. doi:https://doi.org/10.1016/j.pain.2005.03.023.View ArticlePubMedGoogle Scholar
  60. Kong J, Kaptchuk TJ, Polich G, Kirsch I, Gollub RL. Placebo analgesia: findings from brain imaging studies and emerging hypotheses. Rev Neurosci. 2007;18(3–4):173–90.PubMedGoogle Scholar
  61. Goldman-Rakic PS. Topography of cognition: parallel distributed networks in primate association cortex. Annu Rev Neurosci. 1988;11:137–56. doi:https://doi.org/10.1146/annurev.ne.11.030188.001033.View ArticlePubMedGoogle Scholar
  62. Parks EL, Geha PY, Baliki MN, Katz J, Schnitzer TJ, Apkarian AV. Brain activity for chronic knee osteoarthritis: dissociating evoked pain from spontaneous pain. Eur J Pain. 2011;15(8):843 e1–14. doi:https://doi.org/10.1016/j.ejpain.2010.12.007.
  63. Howard MA, Sanders D, Krause K, O’Muircheartaigh J, Fotopoulou A, Zelaya F, et al. Alterations in resting-state regional cerebral blood flow demonstrate ongoing pain in osteoarthritis: an arterial spin-labeled magnetic resonance imaging study. Arthritis Rheum. 2012;64(12):3936–46. doi:https://doi.org/10.1002/art.37685.View ArticlePubMedGoogle Scholar
  64. Miller EK. The prefrontal cortex and cognitive control. Nat Rev Neurosci. 2000;1(1):59–65.View ArticlePubMedGoogle Scholar
  65. Corbetta M, Shulman GL. Control of goal-directed and stimulus-driven attention in the brain. Nat Rev Neurosci. 2002;3(3):201–15.View ArticlePubMedGoogle Scholar
  66. Ridderinkhof KR, Ullsperger M, Crone EA, Nieuwenhuis S. The role of the medial frontal cortex in cognitive control. Science. 2004;306(5695):443–7.View ArticlePubMedGoogle Scholar
  67. Matsumoto K, Tanaka K. Neuroscience. Conflict and cognitive control. Science. 2004;303(5660):969–70.View ArticlePubMedGoogle Scholar
  68. Carter CS, van Veen V. Anterior cingulate cortex and conflict detection: an update of theory and data. Cogn Affect Behav Neurosci. 2007;7(4):367–79.View ArticlePubMedGoogle Scholar
  69. Crone EA, Wendelken C, Donohue SE, Bunge SA. Neural evidence for dissociable components of task-switching. Cereb Cortex. 2006;16(4):475–86.View ArticlePubMedGoogle Scholar
  70. Fields H. State-dependent opioid control of pain. Nat Rev Neurosci. 2004;5(7):565–75.View ArticlePubMedGoogle Scholar
  71. Tracey I, Ploghaus A, Gati JS, Clare S, Smith S, Menon RS, et al. Imaging attentional modulation of pain in the periaqueductal gray in humans. J Neurosci. 2002;22(7):2748–52.PubMedGoogle Scholar
  72. Jensen KB, Loitoile R, Kosek E, Petzke F, Carville S, Fransson P, et al. Patients with fibromyalgia display less functional connectivity in the brain’s pain inhibitory network. Mol Pain. 2012;8(1):32.PubMed CentralView ArticlePubMedGoogle Scholar
  73. Yu R, Gollub R, Spaetha R, Napadowa V, Wasana A, Kong J. Disrupted functional connectivity of the periaqueductal gray in chronic low back pain. NeuroImage. 2014;6:100–8.Google Scholar
  74. Eippert F, Bingel U, Schoell ED, Yacubian J, Klinger R, Lorenz J, et al. Activation of the opioidergic descending pain control system underlies placebo analgesia. Neuron. 2009;63(4):533–43.View ArticlePubMedGoogle Scholar
  75. Kong J, Tu PC, Zyloney C, Su TP. Intrinsic functional connectivity of the periaqueductal gray, a resting fMRI study. Behav Brain Res. 2010;211(2):215–9.PubMed CentralView ArticlePubMedGoogle Scholar
  76. Kucyi A, Salomons TV, Davis KD. Mind wandering away from pain dynamically engages antinociceptive and default mode brain networks. Proc Natl Acad Sci USA. 2013;110(46):18692–7. doi:https://doi.org/10.1073/pnas.1312902110.PubMed CentralView ArticlePubMedGoogle Scholar
  77. Salomons TV, Johnstone T, Backonja MM, Davidson RJ. Perceived controllability modulates the neural response to pain. J Neurosci. 2004;24:7199–203.View ArticlePubMedGoogle Scholar
  78. Egorova N, Gollub R, J K. Repeated verum but not placebo acupuncture normalizes connectivity in brain regions dysregulated in chronic pain. NeuroImage. 2015. doi:https://doi.org/10.1016/j.nicl.2015.09.012.
  79. Chen XY, Spaeth RB, Retzepi K, Ott D, Kong J. Acupuncture modulates cortical thickness and functional connectivity in knee osteoarthritis patients. Sci Rep 2014;4(6482). doi:https://doi.org/10.1038/srep06482.
  80. Rainville P, Duncan GH, Price DD, Carrier B, Bushnell MC. Pain affect encoded in human anterior cingulate but not somatosensory cortex. Science. 1997;277:968–71.View ArticlePubMedGoogle Scholar
  81. Price DD. Psychological and neural mechanisms of the affective dimension of pain. Science. 2000;288:1769–72.View ArticlePubMedGoogle Scholar
  82. Kong J, Gollub RL, Polich G, Kirsch I, Laviolette P, Vangel M et al. A functional magnetic resonance imaging study on the neural mechanisms of hyperalgesic nocebo effect. J Neurosci. 2008;28(49):13354–62 (PMCID: PMC2649754).Google Scholar
  83. Zhang Y, Meng X, Li A, Xin J, Berman BM, Lao L, et al. Electroacupuncture alleviates affective pain in an inflammatory pain rat model. Eur J Pain. 2012;16(2):170–81. doi:https://doi.org/10.1016/j.ejpain.2011.07.002.PubMed CentralView ArticlePubMedGoogle Scholar
  84. Fang J, Jin Z, Wang Y, Li K, Kong J, Nixon EE, et al. The salient characteristics of the central effects of acupuncture needling: limbic-paralimbic-neocortical network modulation. Hum Brain Mapp. 2009;30(4):1196–206.View ArticlePubMedGoogle Scholar
  85. Qin W, Bai L, Dai J, Liu P, Dong M, Liu J, et al. The temporal-spatial encoding of acupuncture effects in the brain. Mol Pain. 2011;7:19. doi:https://doi.org/10.1186/1744-8069-7-19.PubMed CentralView ArticlePubMedGoogle Scholar
  86. Kong J, Ma L, Gollub RL, Wei J, Yang X, Li D, et al. A pilot study of functional magnetic resonance imaging of the brain during manual and electroacupuncture stimulation of acupuncture point (LI-4 Hegu) in normal subjects reveals differential brain activation between methods. J Altern Complement Med. 2002;8(4):411–9. doi:https://doi.org/10.1089/107555302760253603.View ArticlePubMedGoogle Scholar
  87. Pariente J, White P, Frackowiak RS, Lewith G. Expectancy and belief modulate the neuronal substrates of pain treated by acupuncture. Neuroimage. 2005;25(4):1161–7.View ArticlePubMedGoogle Scholar
  88. Eickhoff SB, Schleicher A, Zilles K, Amunts K. The human parietal operculum. I. Cytoarchitectonic mapping of subdivisions. Cereb Cortex. 2006;16(2):254–67. doi:https://doi.org/10.1093/cercor/bhi105.View ArticlePubMedGoogle Scholar
  89. Menon V, Uddin LQ. Saliency, switching, attention and control: a network model of insula function. Brain Struct Funct. 2010;214(5–6):655–67. doi:https://doi.org/10.1007/s00429-010-0262-0.PubMed CentralView ArticlePubMedGoogle Scholar
  90. Melzack R. Folk Medicine and the sensory modulation of pain. In: Wall PD, Melzack R, editors. Textbook of Pain. London: Churchill Livingstone; 1994. p. 1209–17.Google Scholar
  91. Dumville J, Hahn S, Miles J, Torgerson D. The use of unequal randomisation ratios in clinical trials: a review. Contemp Clin Trials. 2006;27:1–12.View ArticlePubMedGoogle Scholar
  92. Streitberger K, Kleinhenz J. Introducing a placebo needle into acupuncture research. Lancet. 1998;352:364–5.View ArticlePubMedGoogle Scholar
  93. Kong J, Fufa DT, Gerber AJ, Rosman IS, Vangel MG, Gracely RH, et al. Psychophysical outcomes from a randomized pilot study of manual, electro, and sham acupuncture treatment on experimentally induced thermal pain. J Pain. 2005;6(1):55–64. doi:https://doi.org/10.1016/j.jpain.2004.10.005.View ArticlePubMedGoogle Scholar
  94. Zhu D, Gao Y, Chang J, Kong J. Placebo acupuncture devices: considerations for acupuncture research. Evid Based Complement Alternat Med. 2013;2013:628907.PubMed CentralPubMedGoogle Scholar
  95. Smith SM, Jenkinson M, Woolrich MW, Beckmann CF, Behrens TE, Johansen-Berg H, et al. Advances in functional and structural MR image analysis and implementation as FSL. Neuroimage. 2004;23(Suppl 1):S208–19.View ArticlePubMedGoogle Scholar
  96. Beckmann CF, Smith SM. Probabilistic independent component analysis for functional magnetic resonance imaging. IEEE Trans Med Imaging. 2004;23(2):137–52.View ArticlePubMedGoogle Scholar
  97. Xue T, Yuan K, Zhao L, Yu D, Dong T, Cheng P, et al. Intrinsic brain network abnormalities in migraines without aura revealed in resting-state fMRI. PLoS One. 2012;7(12):e52927. doi:https://doi.org/10.1371/journal.pone.0052927.PubMed CentralView ArticlePubMedGoogle Scholar
  98. Russo A, Tessitore A, Giordano A, Corbo D, Marcuccio L, De Stefano M, et al. Executive resting-state network connectivity in migraine without aura. Cephalalgia. 2012;32(14):1041–8. doi:https://doi.org/10.1177/0333102412457089.View ArticlePubMedGoogle Scholar
  99. Bai L, Qin W, Tian J, Liu P, Li L, Chen P, et al. Time-varied characteristics of acupuncture effects in fMRI studies. Hum Brain Mapp. 2009;30(11):3445–60. doi:https://doi.org/10.1002/hbm.20769.View ArticlePubMedGoogle Scholar
  100. Bai L, Tian J, Zhong C, Xue T, You Y, Liu Z, et al. Acupuncture modulates temporal neural responses in wide brain networks: evidence from fMRI study. Mol Pain. 2010;6:73. doi:https://doi.org/10.1186/1744-8069-6-73.PubMed CentralView ArticlePubMedGoogle Scholar
  101. Wiech K, Jbabdi S, Lin CS, Andersson J, Tracey I. Differential structural and resting state connectivity between insular subdivisions and other pain-related brain regions. Pain. 2014;155(10):2047–55. doi:https://doi.org/10.1016/j.pain.2014.07.009.PubMed CentralView ArticlePubMedGoogle Scholar
  102. Wiech K, Lin CS, Brodersen KH, Bingel U, Ploner M, Tracey I. Anterior insula integrates information about salience into perceptual decisions about pain. J Neurosci. 2010;30(48):16324–31.View ArticlePubMedGoogle Scholar

Copyright

© Chen et al. 2015

Advertisement