|Home | About | Journals | Submit | Contact Us | Français|
Previous studies have defined low-frequency, spatially consistent networks in resting fMRI data which may reflect functional connectivity. We sought to explore how a complex somatosensory stimulation, acupuncture, influences intrinsic connectivity in two of these networks: the default mode network (DMN) and sensorimotor network (SMN). We analyzed resting fMRI data taken before and after verum and sham acupuncture. Electrocardiography data was used to infer autonomic modulation through measures of heart rate variability (HRV). Probabilistic independent component analysis was used to separate resting fMRI data into DMN and SMN components. Following verum, but not sham, acupuncture there was increased DMN connectivity with pain (anterior cingulate cortex (ACC), periaqueductal gray), affective (amygdala, ACC), and memory (hippocampal formation, middle temporal gyrus) related brain regions. Furthermore, increased DMN connectivity with the hippocampal formation, a region known to support memory and interconnected with autonomic brain regions, was negatively correlated with acupuncture-induced increase in a sympathetic related HRV metric (LFu), and positively correlated with a parasympathetic related metric (HFu). Following verum, but not sham, acupuncture there was also increased SMN connectivity with pain related brain regions (ACC, cerebellum). We attribute differences between verum and sham acupuncture to more varied and stronger sensations evoked by verum acupuncture. Our results demonstrate for the first time that acupuncture can enhance the post-stimulation spatial extent of resting brain networks to include anti-nociceptive, memory, and affective brain regions. This modulation and sympathovagal response may relate to acupuncture analgesia and other potential therapeutic effects.
Acupuncture is an ancient Eastern healing modality which may have therapeutic effects in the treatment of chronic pain [8, 50] and other illnesses. A review of the literature suggests that acupuncture influences a diverse set of bodily organs and functions including brain processing [25, 26] and autonomic nervous system (ANS) activity [1, 44].
It has long been known that the analgesic effects of acupuncture may actually peak long after cessation of active needle stimulation . Furthermore, clinical acupuncture typically involves an active needle manipulation phase followed by a longer period of rest. Thus, although many fMRI studies have investigated brain response to acute needle stimulation, none have evaluated acupuncture effects on functional connectivity that may extend beyond the active needle stimulation phase, i.e. sustained acupuncture effects.
To investigate the sustained effects of acupuncture into post-stimulus rest, we evaluated functional connectivity changes in “resting state networks,” RSNs. Previous fMRI studies have found that in a task-free state (i.e. rest), multiple distributed brain areas demonstrate temporal correlation or intrinsic “functional connectivity” in low frequency ranges [11, 22, 32, 46]. For example, studies have found correlation in resting fMRI signal from sensorimotor cortices of opposite hemispheres . This RSN has been referred to as the sensorimotor network, or SMN . Resting connectivity has also been described in the default mode network (DMN) [14, 30, 35], which involves brain regions putatively engaged in self-referential cognition that are “deactivated” during a variety of externally focused task conditions [for review see 38]. Pain is known to interact with the DMN. Both acute pain  and acupuncture  are known to induce deactivation in DMN regions, while chronic pain may be associated with less pronounced DMN deactivation in response to a visual attention task . Furthermore, perception of somatosensory stimuli near sensory threshold is facilitated by decreased DMN activity in a brief pre-event resting period .
Our study aimed to resolve the change in resting functional connectivity following a complex somatosensory stimulation, acupuncture, using probabilistic Independent Component Analysis (pICA), which can separate fMRI data into independent RSNs, including the DMN and SMN . Furthermore, recent reports have suggested that ANS response to acupuncture (measured by heart rate variability, HRV, ) may be linked to clinical response for different chronic pain populations . As variable peripheral ANS response to acupuncture may influence post-stimulation brain activity, we also sought to explore the relationship between acupuncture induced ANS outflow to the heart as estimated by HRV and acupuncture induced changes in RSN connectivity. To our knowledge, this study constitutes the first investigation of the sustained effects of somatosensory stimulation on resting state connectivity.
Data was collected for 15 healthy, right handed  adults (8 female), 21−33 years of age. Subjects were recruited via fliers/newsletters adhering to MGH guidelines for distribution at neighboring academic institutions and hospitals. Subjects were screened to assure their safety and compatibility with MRI recording.
The present study involved verum manual acupuncture (ACUP) and sham acupuncture (SHAM). Both ACUP and SHAM consisted of a total of 150 seconds of active stimulation in a blocked design paradigm. A 5.5 minute rest run was completed before, and after, each active stimulation run (Figure 1). During rest runs, subjects were asked to lie quietly and fixate on a centrally presented “+”. Past studies have demonstrated robust RSN activity for both eyes closed, and eyes open, with and without fixation . Fixation was chosen in order to minimize drowsiness and the potential for sleep.
The order of ACUP and SHAM runs was randomized across subjects to mitigate order effects. These two runs (and their respective rest runs) were separated by structural scans which lasted at least 15 minutes total in order to maximize washout of acute stimulation effects within the scan session. Data from the rest run immediately following acupuncture were compared to data from the rest run immediately preceding stimulation (i.e. after – before acupuncture).
All acupuncture was performed by the same licensed (and experienced) acupuncturist at left PC-6. This acupoint, which is on the volar aspect of the forearm, is 2 cun (~5cm) proximal to the transverse wrist crease, between the tendons of the palmaris longus and flexor carpi radialis muscles. For ACUP, acupuncture was performed by first inserting a non-magnetic (pure silver), 0.23mm diameter, 30mm long acupuncture needle (Asahi Industry, Inc., Kawaguchi, Japan) into PC-6 to a depth of approximately 1.5cm. During stimulation, the needle was manually twirled (~±180°) at 0.5 Hz.
For SHAM, needle insertion was first simulated by poking the skin over PC-6 with a 5.88 von Frey monofilament passed through a needle guide tube (similar to the one used for verum acupuncture). During the scan run, the acupuncturist performed non-insertive cutaneous stimulation over the acupoint (tapping at 0.5 Hz) using the same monofilament. Subjects were acupuncture-naive and were not informed of a sham acupuncture condition, only that there would be “different forms” of acupuncture during fMRI. Subjects lay supine in the scanner with their vision of distal body regions blocked by the MRI head coil, preventing them from viewing the intervention occurring at their periphery. The issue of appropriate controls in acupuncture research remains controversial . Each type of control addresses different aspects of the acupuncture intervention. For example, non-insertive stimulation at the acupoint controls for the cutaneous somatosensory (and potentially cognitive) aspect of the intervention. Thus, any differences found between verum and sham acupuncture are most likely attributable to needle penetration and stimulation of deeper, subcutaneous mechanoreceptors and/or nociceptors. In contrast, insertive needling at a nonacupoint controls for the insertive nature, and potentially cognitive aspects, of the intervention. Thus, any differences between verum and sham might be more strongly attributed to acupoint-specificity, which has been posited by traditional East Asian medical theory. In this experiment, we chose to emphasize the stimulation of deep receptors as the active component of acupuncture and used a non-insertive sham acupuncture control, delivered over the same cutaneous location on the body. Our sham stimulation aimed to control for both superficial, cutaneous somatosensory effects over the acupoint, as well as the cognitive processing induced by subjects expecting “acupuncture” stimulation.
FMRI Data were acquired using a 3T Siemens Tim Trio MRI System (Siemens Medical, Erlangen, Germany) equipped for echo planer imaging with a 12 channel head coil. During scanning, subjects remained in the supine position with their heads immobilized by cushioned supports. Subjects wore earplugs throughout the experiment to attenuate MRI gradient noise. Structural scans consisted of a T1-weighted anatomical (TR/TE = 2.73/3.19 ms, flip angle = 7°, FOV = 256×256 mm; slice thickness = 1.33 mm) and a multi-echo fieldmap scan (f.a.=55°, TR/TE1/TE2 = 500/3.38/5.84 ms). Blood oxygenation level-dependent (BOLD) functional imaging was performed using a gradient echo T2*-weighted pulse sequence (TR/TE=3000/30 ms, flip angle = 90°, FOV = 200 × 200 mm, 48 coronal slices, slice thickness = 3.0 mm with 0.6 mm interslice gap, matrix = 64 × 64, 110 time points for a total of 5.5minutes). Image collection was preceded by 5 dummy scans. Coronal slices were used to minimize susceptibility artifact (gradient distortion) in subcortical limbic regions .
Data analysis was performed using a combination of analysis packages including FSL (FMRIB's Software Library) and AFNI. Functional data were pre-processed to correct for magnetic field inhomogeneities caused by magnetic susceptibility differences in neighboring tissues within the head (FSL-FUGUE) and motion corrected to compensate for any head movements using a linear (affine) transformation procedure (FSL-FLIRT). Brain extraction was performed on functional data using FSL-BET while skull stripping of structural data for alignment utilized Freesurfer Software [20, 28]. Functional data were smoothed using a Gaussian kernel of FWHM 5 mm; and high-pass temporal filtering (f = 1/250 s) was also performed.
Independent component analysis (ICA) was performed on all rest runs using FSL-MELODIC . The use of probabilistic Independent Component Analysis (pICA) to derive RSNs has been successfully performed in the past by several groups [6, 34] in order to evaluate functional connectivity networks. The independent components (IC) maps are in the form of a Z-statistic for which activated/deactivated areas may be considered temporally coherent. While some groups have utilized seed-voxel approaches to investigate resting state functional connectivity, these methods require a priori selection of regions of interest which imposes a subjective element to data analysis. ICA also has the advantage of being able to better separate/remove low frequency components due to physiological noise which may more readily contaminate seed voxel analyses of functional connectivity .
In the present analysis, the number of output IC's was not restricted and subsequent component masking and selection procedures were similar to those applied in Greicius et al. (2007). Specifically, the ICA components were first screened for the presence of high frequencies (> 0.1Hz) by rejecting those components which had more than 25% of their total Fourier spectrum power in frequencies above 0.1Hz. Our investigation focused on acupuncture modulation of DMN and SMN activity, as regions within these networks have been found to respond in our block-design fMRI studies of acupuncture . Templates of both the DMN and SMN were then used to select the “best-fit” components for each rest run in all subjects. Typically, the DMN includes the inferior parietal lobule (IPL) (~BA 40, 39), the posterior cingulate (~BA 30, 23, 31) and precuneus (~BA 7), areas of the inferior, medial and superior frontal gyri (~BA 8, 9, 10, 47), the hippocampal formation, and the lateral temporal cortex (~BA 21) . The SMN includes bilateral primary and secondary somatosensory areas (SI and SII), primary motor cortex (MI), and the SMA [6, 11]. The templates were RSN group maps (thresholded at p<0.05) from healthy adults which were provided by Beckmann , who also used pICA for their derivation (see Figure 2).
The best-fit component was selected by calculating the average z-score of voxels both inside and outside of the template then selecting the component for which the ratio of values inside/outside was the greatest. Subsequent group analysis was performed on the selected component maps using a mixed effects model (FSL-FEAT). It is important to note that these component maps contained whole-brain data and were un-masked, thereby preserving any significant connectivity that regions within the DMN and SMN may have with regions outside the classically defined borders outlined in the templates. Paired t-tests were performed to determine differences in spatial extent between resting state networks (DMN and SMN) before vs. after the acupuncture runs, thresholded at a corrected cluster significance level of p=0.05. Unpaired t-tests versus nil were performed to derive group maps for SMN and DMN connectivity both before and after acupuncture, also thresholded at a corrected cluster significance level of p=0.05. Anatomical localization was assisted by use of a brain atlas . Cerebellar localization was aided by recent work in MRI-guided parcellation .
In order to infer cardiac sympathovagal modulation by acupuncture, subjects were outfitted with MRI-compatible chest ECG electrodes and a nasal cannula to measure end tidal CO2. Electrocardiogram (ECG) and respiratory data were collected at 200Hz with a system specially designed to mitigate MRI-induced artifact in the ECG trace (InVivo Magnitude CV, Invivo Research Inc., Orlando, Florida), on a laptop acquisition system (Labview 7.0, DAQCard 6024E, National Instruments, Austin, Texas) for further analysis.
ECG data were processed with the WFDB (WaveForm DataBase) Software Package  and MATLAB 6.5 (The Mathworks, Inc. Natick, MA). Data were automatically annotated with careful manual correction in order to form an RR-interval time series. A power spectrum of this time series was then calculated using the maximum entropy (all-poles) method . Low frequency (LF) and high frequency (HF) power were estimated by integrating within the 0.04−0.15Hz and 0.15−0.50Hz frequency bands, respectively. These metrics were calculated in normalized units by dividing the LF and HF values by the total power minus that in the very low frequency range (0−0.04Hz), and termed LFu and HFu . While some controversy in interpretation remains, LFu is more influenced by cardio-sympathetic activity, while HFu is influenced solely by cardiovagal activity . The respiratory data were used to ensure that respiratory rate did not drop below 0.15Hz or rise above 0.40Hz, a requisite for subdividing the LF and HF power bands and attributing sympathetic and parasympathetic interpretation to LFu, and HFu. For correlation analysis, normalized individual response to ACUP and SHAM was calculated by subtracting the value in the pre-stimulation rest from the value during active stimulation. To investigate the relation between these metrics and imaging results (post-pre change in mean z-statistic within ROIs significant from the paired t-test above), Pearson correlation coefficients were calculated, significant at p<0.05.
Following the ACUP and SHAM portions of the scan session, subjects rated the intensity of sensations they felt during the active stimulation run. Subjects were presented with a 10 point visual analog scale (VAS) and were asked to rate the intensity of sensations commonly associated with the experience of deqi¥, and listed on the MGH Acupuncture Sensation Scale (MASS). Subjects also assessed anxiety/relaxation on a scale ranging from very anxious (5) to fully relaxed (−5). Additionally, subjects were asked to assess the extent of “spreading” that may have occurred for any of the listed sensations. A modified version of this procedure has been successfully used by our group in the past to assess psychophysical response in conjunction with neuroimaging [41, 54]. In order to quantify the total amplitude of deqi experienced we used the MASS-Index, which is calculated by an exponentially decreasing weighted sum of all sensations . This index attempts to balance breadth and depth of sensations as well as the number of different sensations chosen by the subject. The MASS index was compared between stimulation groups using a paired t-test, significant at p<0.05 (SPSS 10.0.7, Chicago, Illinois). Furthermore, frequency counts of specific sensations were compared between stimulation type with a Pearson Chi-squared test, significant at p<0.05. An omnibus test (Fisher's combined probability test) was used to test if the commonality of sensations elicited by ACUP differed from those elicited by SHAM.
To investigate the relation between the intensity of elicited sensations and imaging results (post-pre change in mean z-statistic within ROIs significant from the paired t-test above), Pearson correlation coefficients were calculated. Sensations included in this analysis were those reported by greater than 60% of subjects, as too many nil values would adversely skew the distribution. Results were deemed significant at p<0.05, Bonferroni corrected for multiple comparisons due to the large number of sensations tested (9).
Finally, in order to investigate the persistence of sensations into the rest run following active stimulation, we also asked subjects to rate the intensity of sensation during the rest run before and after verum and sham stimulation. For these data, sensation intensity difference between verum and sham acupuncture was contrast using a paired t-test, significant at p<0.05.
Group maps of the best fitting independent components for the DMN and SMN demonstrated consistently robust spatial distribution relative to their respective templates, both before and after verum (Figure 2) and sham (data not shown) acupuncture. The best fitting component for the DMN demonstrated resting connectivity within the inferior parietal lobule, posterior cingulate, and medial areas of the inferior, middle and superior frontal gyri, as well as the precuneus. The best fitting component for the SMN involved bilateral primary somatosensory and motor cortices, secondary somatosensory cortex and the supplementary motor area (SMA).
When comparing DMN connectivity after ACUP with that before, there was increased connectivity of this network to a number of areas (Figure 3, Table 1) including the amygdala, hippocampal formation, periaqueductal gray (PAG), substantia nigra (SN), middle temporal gyrus (MTG, ~BA 21), supplementary motor area (SMA), and anterior cingulate (ACC, ~BA 24), posterior parietal (~BA 7), and primary visual (V1) cortices. There were no areas of decreased connectivity following ACUP. Conversely, after SHAM, the DMN demonstrated increased connectivity with only the temporo-occipital cortex (~BA 37/39) and decreased connectivity with areas of the middle (~BA 21) and inferior temporal (~BA 20) gyri (Figure 3, Table 1). This “decreased” connectivity resulted from a greater z-score before SHAM than after (2 rightmost columns, Table 1)
After ACUP, the SMN demonstrated increased connectivity with the anterior cingulate cortex (~BA 32/24), pre-SMA (~BA 8/6), and cerebellum (simplex, lateral hemispheric zone 1, VI l) (Figure 4, Table 1). Conversely, after SHAM, the SMN demonstrated decreased connectivity with the dorsolateral prefrontal cortex (~BA 8, greater before versus after, Figure 4, Table 1).
In order to explore the relationship between changes in brain functional connectivity and acupuncture-induced ANS modulation, we calculated the inter-subject Pearson correlation coefficient between the change in HRV metrics (LFu, HFu) during acupuncture minus the pre-stimulus rest and the change in z-statistic in ROIs taken from significant regions from the paired t-test above for both DMN and SMN. The metrics demonstrated consistent correlation with increasing hippocampal formation connectivity to the DMN following ACUP (LFu: r = −0.73, HFu: r = 0.68, p<0.01, Figure 5). This result suggests that increased parasympathetic and decreased sympathetic modulation during ACUP is associated with increased post-ACUP DMN connectivity with the hippocampal formation. In addition, the change in post-ACUP DMN connectivity with the MTG also demonstrated a negative correlation with LFu (r = −0.70, p<0.01). These data passed the Grubb's test for outliers (p<0.05). No other significant correlations were found for either ACUP or SHAM, though it should be remembered that correlations were tested only for regions significant in the fMRI group analysis (paired t-test). Thus, while regions such as the insula (particularly the right anterior insula) are associated with the central autonomic network and have demonstrated brain response correlated with HRV , the insula was not a region which demonstrated either increased or decreased connectivity to either the DMN or SMN following acupuncture and was not included in this analysis.
A statistical analysis found that MASS Index, a measure of deqi intensity, was greater for ACUP compared to SHAM stimulation (paired t-test, p<0.01). Furthermore, there was a greater intensity (paired t-test, p<0.05) of spreading for ACUP (3.6±2.4) compared to SHAM (1.7±1.5). Interestingly, SHAM was associated with greater relaxation (−1.2±2.4, μ±σ) compared to ACUP (0.2±2.3, μ±σ). There was no significant difference in the intensity of sharp pain (ACUP: 3.7±2.7; SHAM: 4.9±2.6) or throbbing (ACUP: 3.9±1.9; SHAM: 2.6±2.4). An omnibus test (Fisher's combined probability test) found that ACUP and SHAM also differed in the commonality of sensations, p(n=14) < 0.001. Differences also existed in the prevalence of specific deqi sensations elicited (Figure 6). The prevalence of “aching” (ACUP: 86.7% of subjects, SHAM: 40.0%, p<0.01), “fullness” (ACUP: 46.7%, SHAM: 13.3%, p<0.05), and “dull pain” (ACUP: 93.3%, SHAM: 40.0%, p<0.005) was found to be greater for ACUP (uncorrected for multiple comparisons due to significant omnibus test). It should also be noted that while sharp pain was reported for both conditions, when debriefed, most subjects who reported sharp pain described this sensation as only transient, occurring typically at the start of a stimulation block.
We also wanted to explore the relationship between changes in brain functional connectivity and the intensity of sensations evoked by acupuncture. We calculated the inter-subject Pearson correlation coefficient between the intensity of acupuncture-evoked sensations and the change in z-statistic in ROIs taken from the paired t-test above for both DMN and SMN. We found a trend for increasing ACUP evoked “soreness” correlating with increasing SMN connectivity with pre-SMA (r = 0.65, p=0.081, Bonferroni corrected). This result is consistent with the hypothesis that changes in resting SMN connectivity, a network containing brain regions subserving the sensory-discriminative aspects of pain, are related to the sensations evoked by acupuncture stimulation. No significant correlations were found for DMN ROIs, nor for any ROIs from the sham acupuncture paired t-test.
For the rest runs following ACUP and SHAM, sensation intensity scores were collected for 11 subjects. Of these, four subjects reported feeling any sensation at some point during the rest run immediately following ACUP (0.9±1.3, overall μ±σ) and two following SHAM (0.4±0.8). Ultimately, there was no significant difference between mean sensation intensity in the rest runs following ACUP versus those following SHAM (p>0.1).
Acupuncture has shown promise in treating chronic pain and other disorders . Acupuncture analgesia for chronic pain is typically delayed from the needle stimulation phase and may result from modulation of brain processing and/or autonomic response sustained into the post-stimulation period. We used pICA to evaluate acupuncture modulation of functional connectivity within RSNs consistent with the DMN and SMN. Following ACUP, but not SHAM, there was increased DMN connectivity with limbic anti-nociceptive (ACC, PAG), affective (amygdala, ACC), and memory (hippocampal formation, MTG) related brain regions. Furthermore, increased DMN connectivity with the hippocampal formation, a region known to support memory and interconnected with autonomic motor nuclei, was anti-correlated with acupuncture-induced increase in a sympathetic related HRV metric (LFu), and correlated with a parasympathetic related HRV metric (HFu). In addition, following ACUP, there was greater connectivity between SMN and pain processing and motor planning brain regions (ACC, pre-SMA, and cerebellum), changes not seen following SHAM. These results demonstrate for the first time that a somatosensory stimulation can have sustained effects on brain resting state connectivity.
Previous studies have defined the DMN as a set of regions demonstrating increased hemodynamic activity at baseline relative to a variety of active (i.e. goal oriented) task conditions . The precise function of the DMN is debated, but evidence suggests that it may support autobiographical memory and decreased externally focused attention . Importantly, the DMN is mainly comprised of cortical “association areas” and thus has extensive anatomical connections with not only primary sensory areas but also subcortical regions involved in pain processing, affect, and memory encoding/retrieval.
Our results found that ACUP enhanced resting DMN connectivity with the ACC, PAG, and SN - areas putatively involved in affective pain processing/evaluation, pain anticipation and endogenous anti-nociception. The dorsal subregion of the ACC (BA 24) is activated by acute painful stimuli . This limbic area is also activated in endogenous anti-nociception , affective pain processing , and attention to painful stimuli . Moreover, ACUP increased resting DMN connectivity with other anti-nociceptive regions such as the PAG. Animal studies of acupuncture suggest that opioidergic modulation of the PAG plays a critical role in acupuncture analgesia via inhibition of afferent pain signaling . In addition, acupuncture related opioidergic modulation of the PAG also supports sympatho-inhibition [37, 61]. It is likely that these functions are inter-related, as the sympathetic nervous system is known to up-regulate pain circuitry and vice versa . The SN is the primary source of dopaminergic (DA) input to the striatum. The DA system projects diffusely to many brain regions, and the SN has been implicated in pain and non-pain somatosensory perception, especially when stimuli have high behavioral context [for review see 16]. We propose that increased DMN connectivity with brain areas involved in anti-nociception and sympatho-inhibition forms the basis of an integrated network for top-down regulation of acupuncture analgesia.
Other areas demonstrating increased connectivity with DMN following ACUP are most often implicated in memory encoding/retrieval (hippocampal formation) and affective memory (amygdala) but are also modulated during experimental pain [4, 9, 13, 24]. Our previous block design fMRI studies demonstrated that acupuncture reduces activity within the amygdala and hippocampal formation , a potential region of the DMN [14, 36]. Our data suggests that this reduction may be sustained minutes after acupuncture stimulation and manifest as increased connectivity between these regions and the DMN. Interestingly, it has been suggested that chronic pain leads to ‘pain memories’ which can persist long after the peripheral pain source has resolved, leading to centrally driven pain [29, 43]. Non-pain somatosensory memories have also been noted . We propose that acupuncture sensations, which combine dull pain and deep, unusual somato-sensations, can form somatosensory memories by binding the DMN with pain-related, affective and memory related regions. If this sustained connectivity can be coupled with brain processing in anti-nociceptive brain regions, over multiple treatment sessions, the linkage may help erode centrally driven chronic pain memories. While speculative, this hypothesis is supported by clinical acupuncture practice, which typically demands multiple, regularly spaced treatment sessions.
Increased DMN connectivity in the above regions was not found following another complex somatosensory intervention, SHAM. Our psychophysical data demonstrated that compared to SHAM, ACUP induced more intense and more varied deqi sensations which spread more extensively away from the acupoint site. These differences in psychophysics were consistent with our past studies  and other data comparing verum and non-insertive shams (e.g. placebo-needles), which may better mimic needle insertion than needle manipulation . Furthermore, only a small number of subjects experienced mild deqi sensations during rest (following both ACUP and SHAM). Thus, our results are likely due to the sensations experienced during stimulation and not bottom-up somatosensory persistence for either ACUP or SHAM. It is possible that increased somatosensory saliency (greater deqi) for ACUP promotes connectivity of DMN with areas supporting sensory stimulus encoding seen in our data, such as the posterior parietal cortex and SMA. Interestingly, tasks of greater cognitive complexity  and pain intensity  have been shown to more readily deactivate DMN brain regions. Our results demonstrate that increased somatosensory saliency may also influence DMN connectivity in the resting state minutes after active stimulation.
We also sought to explore the relationship between acupuncture-induced autonomic modulation and subsequent change in RSN connectivity. Increased DMN connectivity with the hippocampal formation was anti-correlated with acupuncture-induced increase in sympathetic modulation, and correlated with parasympathetic modulation. Interestingly, the hippocampus has polysynaptic connectivity with autonomic brainstem nuclei , while vagal nerve stimulation enhances memory , an effect that may be supported by hippocampal activity. Our data contribute to these results by suggesting that when acupuncture induces a sympathovagal shift toward parasympathetic, the hippocampal formation has greater connectivity with the DMN, a brain network thought to subserve self-referential cognition and autobiographical memory . While pain is known to modulate autonomic outflow , the response is typically increased sympathetic tone. However, increased hippocampal/DMN connectivity following acupuncture was associated decreased sympathetic and increased parasympathetic modulation, suggesting that acupuncture modulation of DMN connectivity may not result from classic pain-like aspects of this complex intervention. Interestingly, other groups have found that acupuncture-induced decrease of LF/HF ratio, relating a shift in sympathovagal balance toward greater parasympathetic modulation, may relate to acupuncture therapeutic effects on chronic pain . Our correlative results connecting HRV metrics with enhanced DMN connectivity should to be replicated in a clinical population to investigate this possible connection with clinical efficacy.
Following ACUP, the SMN demonstrated increased connectivity with the ACC, pre-SMA and cerebellum. The subregion of the ACC demonstrating increased connectivity (BA 24/32) may support attention as well as pain/anti-nociceptive processing . The pre-SMA supports motor planning , while the cerebellar simplex is related to cognitive processing  and is modulated by acupuncture  and pain . Enhanced connectivity of these areas with the SMN may support linkage of lower order sensory/discriminative with higher-order affective/evaluative processing inherent to a complex sensory stimulus such as acupuncture.
In the current study we evaluated the sustained effects of acupuncture on brain activity during rest following stimulation. Unfortunately, the delay which constitutes an optimal washout period for acupuncture's neurobiological effects remains unknown. To avoid contamination of the control conditions (SHAM), post-stimulus rest runs were separated by structural scanning, a minimum of 15 minutes. Furthermore, the “before” and “after” rest runs immediately preceded and followed active stimulation, respectively. Thus, our results reflect the direct influence of either ACUP or SHAM. Finally, the order of the ACUP and SHAM portions of the scan session was pseudo-randomized to avoid any order effects.
The current study evaluated acupuncture modulation of resting functional connectivity within the brain's DMN and SMN. Following verum acupuncture there was increased DMN connectivity with pain (ACC, PAG), affective (amygdala, ACC), and memory (hippocampal formation, MTG) related brain regions. Following acupuncture there was also increased SMN connectivity with pain related brain regions (ACC, cerebellum). We attributed differences between verum and sham acupuncture to more varied and stronger sensations evoked by verum acupuncture. Furthermore, increased DMN connectivity with the hippocampal formation, a region known to support memory and interconnected with autonomic brain regions, was anti-correlated with acupuncture-induced increase in a sympathetic related HRV metric (LFu), and correlated with a parasympathetic related HRV metric (HFu). Our results demonstrate for the first time that a complex somatosensory stimulation, acupuncture, can modulate the post-stimulation spatial extent of resting state connectivity networks. Modulation of these networks and sympathovagal response to acupuncture may relate to acupuncture analgesia as well as other acupuncture associated therapeutic effects. Additionally, neuroimaging evaluation of chronic pain states has been notoriously difficult. Intrinsic functional connectivity, which we have demonstrated to be altered following our acute pain-like stimulus, may also be altered and/or dysfunctional in chronic pain patients – a hypothesis we plan to study in the future.
This research was supported by grants from NCCAM, NIH (P01-AT002048, K01-AT002166, F05-AT003770). We would like to thank Randy Buckner and Bruce Rosen for helpful comments during manuscript preparation. There are no conflicts of interest for any author.
Publisher's Disclaimer: This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final citable form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.
¥Deqi corresponds to a multitude of different pain-like and non-pain sensations experienced by a needled subject and is a correlate of effective treatment .