Abstract

Global normalization is often used as a preprocessing step for dispelling the “nuisance effects.” However, it has been shown in cognitive and emotion tasks that this preprocessing step might greatly distort statistical results when the orthogonality assumption of global normalization is violated. The present study examines this issue in fMRI acupuncture studies. Thirty healthy subjects were recruited to evaluate the impacts of the global normalization on the BOLD responses evoked by acupuncture stimulation during De-qi sensation and tactile stimulation during nonpainful sensations. To this end, we compared results by conducting global normalization (PSGS) and not conducting global normalization (NO PSGS) based on a proportional scaling model. The orthogonality assumption of global normalization was violated, and significant changes between BOLD responses for NO PSGS and PSGS were shown in most subjects. Extensive deactivations of acupuncture in fMRI were the non-specifically pernicious consequences of global normalization. The central responses of acupuncture during De-qi are non-specifically activation-dominant at the somatosensory-related brain network, whose statistical power is specifically enhanced by PSGS. In conclusion, PSGS should be unjustified for acupuncture studies in fMRI. The differences including the global normalization or not may partly contribute to conflicting results and interpretations in previous fMRI acupuncture studies.

1. Introduction

Blood oxygen level-dependent (BOLD) functional magnetic resonance imaging (fMRI) is increasingly used to investigate central mechanism of acupuncture [17]. However, few hypotheses have been recognized without debate in fMRI acupuncture studies. Various kinds of differences could bring about the disparate results in fMRI acupuncture studies [8], particularly the differences in data processing and statistical analysis methods [811]. As a questionable data processing step which is at risk of changing the interpretation of the statistical results both qualitatively (changing the sign of the results) and quantitatively [1214], the normalization of the global signal was adopted in many acupuncture fMRI studies [3, 5, 7, 1521]. This paper attempts to evaluate the validity and impact of the normalization of the global signal in acupuncture fMRI studies.

It is known that global normalization was usually employed to dispel global effects, which were considered as “nuisance effects” [12, 13, 22]. Its employment was restricted to the assumption that global effects are orthogonal to the variation of the task-induced effect [12, 13, 22]. In this case, the global signal (estimator of the global effects which is defined by the simple average of all voxels within the brain at each time point) was usually normalized in a proportional scaling or ANCOVA model [1214, 22, 23]. When the task-induced effect is strong enough to contaminate the global signal and to break the orthogonality assumption, the global normalization would be unjustified, often reduce activations, and introduce artificial deactivations [13, 14, 22, 24, 25]. The impact of global normalization in fMRI studies in the domain of emotion research [24], other cognitive tasks [12, 13], and on resting state correlations [25] has been evaluated. So far, however, none of the acupuncture studies in fMRI took the impacts of global normalization into account. Therefore, it raised questions on whether the orthogonality assumption of global normalization is satisfied in acupuncture fMRI studies or not. Furthermore, if violated, the impacts of global normalization on activations and deactivations need to be determined in acupuncture fMRI studies.

The current study set out to (1) assess the validity of using proportional scaling to normalize the global signal in acupuncture fMRI studies and (2) examine the impacts of global normalization on activations and deactivations in acupuncture fMRI studies. To this end, we compared results of the same dataset obtained during acupuncture stimulation by conducting and not conducting global normalization based on a proportional scaling model. Besides acupuncture, a superficial tactile stimulus (noninvasive) was employed to explore the acupuncture-specific results.

2. Materials and Methods

2.1. Subjects

Participants were recruited from a group of 30 college students (15 males and 15 females; ages years). All subjects were right-handed with normal or corrected-to-normal vision. Subjects were acupuncture-naïve, had no history of major medical illnesses, head trauma, neuropsychiatric disorders, or any prescription medications one month preceding the experiment, and did not have any contraindications to exposure to a high magnetic field. All subjects gave written and informed consent after the experimental procedures were fully explained. All research procedures were approved by the West China Hospital Subcommittee on Human Studies and were conducted in accordance with the Declaration of Helsinki.

2.2. Experimental Procedures

For each subject, functional scans of tactile stimulation were collected prior to anatomical scans, and functional scans of acupuncture stimulation were taken after the anatomical scans. The intervals between runs were 5–10 minutes. The order of tactile stimulation and acupuncture stimulation was not given to the subjects. During fMRI scanning, all subjects were instructed to keep their eyes closed to prevent them from actually observing the procedures. Subjects were told that acupuncture was to be performed with different techniques that would generate different sensations during needling. Because the subjects were acupuncture naïve, they would be unable to discriminate the tactile stimulation control from real acupuncture until they experienced real acupuncture.

Acupuncture stimulation was performed at acupoint ST36 on the right leg (Zusanli, located in the tibialis anterior muscle four fingerbreadths below the lower margin of the patella and one fingerbreadth lateral from the anterior crest of the tibia). The fMRI paradigm for the acupuncture stimulation run lasted for 8 minutes and consisted of three one-minute acupuncture manipulations (Figure 1(a)). The needle was inserted perpendicularly to a depth of 2-3 cm before the scan started. A one-minute baseline period was held preceding the first acupuncture stimulation. The interval between the first two acupuncture manipulations was two minutes, while the second and third acupuncture manipulations were separated by an interval of one minute. Scanning was then continued for another minute after the third manipulation. During the acupuncture procedure, the needle was rotated manually clockwise and counterclockwise for one minute at a rate of 60 times per minute. The stimulation was administered by a balanced “tonifying and reducing” technique using a sterile disposable 38 gauge stainless steel acupuncture needle (0.3 mm × 40 mm). After the scan ended, the needle was extracted. The tactile stimulation run was performed at acupoint ST36 on the right leg with a size 5.88 Von Frey monofilament. The paradigm of tactile stimulation was matched to that of the acupuncture stimulation run (Figure 1(b)). During the tactile manipulations, the monofilament tapped the skin gently at a rate of 60 times per minute. The acupuncture and tactile procedures were conducted by the same experienced and licensed acupuncturist.

In the end, the subjects were facilitated by the acupuncturist to quantify their sensations using a 10-point visual analogue scale (VAS) to rate their De-qi experience felt during the acupuncture run [11, 26, 27]. To quantitatively summarize the full multivariate breadth and depth of the De-qi sensations for each subject, the VAS index was calculated [11, 26, 27]. Because sharp pain was considered an inadvertent noxious stimulation and evoked different BOLD responses to the De-qi sensation [3], we excluded subjects from further analysis if they experienced sharp pain. Among the 30 participants, two reported sharp pain during both runs, three experienced sharp pain only in the acupuncture run, and one felt sharp pain only in the tactile run. Thus, 25 subjects in the acupuncture run and 27 subjects in the tactile run remained for data analysis.

2.3. fMRI Scanning Procedure

Imaging data were collected from a 3T Siemens scanner (Allegra; Siemens Medical System) at the Huaxi MR Research Center, West China Hospital of Sichuan University, Chengdu, China. A standard birdcage head coil was used, along with restraining foam pads to minimize head motion and to diminish scanner noise. Thirty axial slices (FOV = 240 mm × 240 mm, matrix = 64 × 64, thickness = 5 mm) parallel to the AC-PC plane covering the whole brain were obtained using a T2*-weighted single-shot, gradient-recalled echo planar imaging (EPI) sequence (TR = 2,000 ms, TE = 30 ms, and flip angle = 90°). The scan covered the entire brain including the cerebellum and brainstem. After the functional run, high-resolution structural information on each subject was acquired using 3D MRI sequences with a voxel size of 1 mm3 for anatomical localization (TR = 2.7 s, TE = 3.39 ms, matrix = 256 × 256, FOV = 256 mm × 256 mm, flip angle = 7°, in-plane resolution = 1 mm × 1 mm, slice thickness = 1 mm).

2.4. fMRI Data Analysis

Preprocessing and statistical analyses at both the individual level and group level were performed using the Statistical Parametric Mapping software (SPM5, http://www.fil.ion.ucl.ac.uk/spm/). Initially, the first 5 time points were discarded in order to avoid the instability of the initial MRI signal. The remaining images were realigned to the first volume. Three subjects in the acupuncture run and four subjects in the tactile run exceeded our rigorous motion threshold of less than 1 mm spatial displacement in any direction. Ultimately, 22 subjects (10 males) in the acupuncture run and 23 subjects (12 males) in the tactile run remained. Subsequently, the images were normalized to the standard EPI template and resampled to a voxel size of 3 × 3 × 3 mm3 and then smoothed spatially using a 6 mm full-width-at-half maximum (FWHM) isotropic Gaussian kernel to decrease spatial noise. The global normalization by proportional scaling was applied for a single processing flow (PSGS) and was not applied to another processing flow (NO PSGS). Then, the time series from each voxel were high-pass filtered (1/235-Hz cutoff) to remove low-frequency noise and signal drift. For each subject, the preprocessed fMRI data were then submitted for fixed-effects model analyses using the general linear model (GLM) performed at each voxel across the whole brain. After acquiring the contrast images, individual level analyses were accomplished, and statistical parametric maps for the statistics (spmT) were then generated for each contrast image. At the group level, the random-effects model analysis was performed based on inference images (i.e., -test for contrast images) from the individual level analysis.

In order to examine whether the orthogonality assumption was satisfied or not, the correlation coefficients ( -value) of the global signal and the reference vector (the experimental paradigm convolved the hemodynamic response function) (viz., “GS-RV -value”) were calculated for each subject. Besides, the averaged -value of the reference vector and time series of the BOLD signal across the whole brain voxels (viz., “AVER-WB-RV -value”) were calculated for each subject for NO PSGS and PSGS, respectively. The changes in the “AVER-WB-RV -value” (PSGS–NO PSGS) were also acquired. Then, the correlations of the “GS-RV -value” and “AVER-WB-RV -value” of NO PSGS, as well as the correlations of the “GS-RV -value” and the changes of the “AVER-WB-RV -value” were calculated.

To evaluate the impact of global normalization, firstly, the number of activations and deactivations for each subject and for the group results was calculated for PSGS and NO PSGS, respectively. At the individual level, two thresholds were selected, a conservative one: , uncorrected and a liberal one: , uncorrected. At the group level, three thresholds were selected, a conservative one: , uncorrected, a moderate one: , uncorrected, and a liberal one: , uncorrected. The percentages of the changes in the activations and deactivations numbers ((PSGS-NO PSGS)/NO PSGS*100%) for each subject and for the group results were calculated, respectively. The ratios of activation divided by deactivation for each subject and for the group results were also compared between NO PSGS and PSGS. The histogram of the whole brain’s -values for the group results was plotted for PSGS and NO PSGS, respectively. Secondly, a paired -test of the BOLD responses was carried out based on the contrast file of each subject between PSGS and NO PSGS. Two thresholds were selected, a moderate one: , uncorrected, and a liberal one: , uncorrected. The histogram of the whole brain’s -values of the paired -test results was also plotted for each group. Thirdly, the conjunction analyses were implemented based on group results map of the one-sample -test for PSGS and NO PSGS. A conservative threshold: , uncorrected, a moderate threshold: , uncorrected, and a liberal threshold: , uncorrected were selected. Four categories were defined for activations and deactivations, respectively, namely, “stronger,” “weaker,” “disappeared” and “arisen.” The “stronger” category was the regions which were significant for both NO PSGS and PSGS, and their -values were larger for PSGS than NO PSGS. The “weaker” category was the regions which were significant for both NO PSGS and PSGS and their -values were smaller for PSGS than NO PSGS. The “disappeared” category was the regions which were significant for NO PSGS but were nonsignificant for PSGS. The “arisen” category was the regions which were nonsignificant for NO PSGS but were significant for PSGS. Finally, the amplitude of BOLD signal changes (ABSC) for PSGS and NO PSGS was calculated, respectively, masked by the “stronger” and, “weaker” regions of activations in the conjunction analyses maps.

For exploring the authentic acupuncture-induced BOLD response, the group results of the one sample -test in the acupuncture run for NO PSGS were mapped and listed. The specific BOLD response evoked by acupuncture stimulation was explored based on a two-sample -test between the acupuncture run and tactile run for NO PSGS.

3. Results

3.1. Correlations among the Reference Vector, Global Signal, and Averaged -Value of the Whole Brain (Figures 2(a) and 2(b) and Figures S1 and S2)

For acupuncture stimulation, the “GS-RV -values” were significant for eight subjects ( , corrected, 6 positive correlations, and 2 negative correlations in Figure 2(a)). For tactile stimulation, seven subjects’ “GS-RV -values” were significant ( , corrected, all positive correlations in Figure 2(b)). Figures and (see Supplementary Material available online at http://dx.doi.org/10.1155/2012/467061) showed the “GS-RV -value” of each subject and plotted the RV and GS of representative subjects. The correlations of the “GS-RV -value” and the “AVER-WB-RV -value” (NO PSGS) are shown in Figures 2(a) and 2(b). Most linear relations were between these two groups of values for acupuncture ( ) and tactile stimulation ( ), indicating that the “GS-RV -value” which averaged the time series and the “AVER-WB-RV -value” which averaged the correlation coefficient were consistent.

3.2. Impact of PSGS at the Individual Level (Figures 2(c)2(f), Table 2)

The correlation of the “GS-RV -value” and the changes in the “AVER-WB-RV -value” (PSGS-NO PSGS) are shown in Figures 2(c) and 2(d). Most linear relations (inverse correlation) were between these two groups of values for acupuncture ( ) and tactile stimulation ( ), indicating that the higher the correlation between GS and RV, the greater the changes (decrease) in the BOLD responses produced by the normalization of global signal. The mean and standard deviation (SD) of the AVER-WB-RV -value across subjects for NO PSGS and PSGS are shown in Figure 2(e). The definite values and the differences between NO PSGS and PSGS are listed in Table S1. The mean values ranged from 0.018 to −0.005 for acupuncture and from 0.036 to −0.011 for tactile stimulation, which were both significantly changed ( for acupuncture and for tactile stimulation).

Figure 2(f) shows the mean and SD of the number of activations and deactivations ( , uncorrected) across subjects for NO PSGS and PSGS. The definite values and the ratio of changes between NO PSGS and PSGS are listed in Table 1. The activations of most subjects were reduced, while the deactivations were increased. The mean ratios of changes of the activations were −33% for both acupuncture and tactile stimulations. The deactivation increased 25% on average for acupuncture and 83% on average for tactile stimulation. Very similar results are manifested in Table which show the number of activations and deactivations at , uncorrected for NO PSGS and PSGS. The ratios of activation divided by deactivation for NO PSGS and PSGS are listed in Table S3. For NO PSGS, the ratio was 194% or 174% for acupuncture and 297% or 189% for tactile stimulation at , uncorrected or , uncorrected, indicating activations were greater in number than deactivations for both stimulations. However, for PSGS, the number of activations was equivalent to and even less than the number of deactivations in both stimulations (the ratio was 105% for acupuncture and 109% for tactile stimulation at , uncorrected, while it was 88% for acupuncture and 78% for tactile stimulation at , uncorrected).

3.3. Impact of PSGS at the Group Level (Figure 3, Table 1, and Table S3)

Histograms of the one sample -test group results for NO PSGS and PSGS are shown in Figure 3. The right part of Figure 3, Table 1, and Table shows the number of group-level activations and deactivations and their ratios of changes for NO PSGS and PSGS. For NO PSGS, the distribution of acupuncture was bell-shaped with a mean greater than zero. The distribution of tactile stimulation for NO PSGS was positively skewed. Group-level deactivations at , uncorrected, were few in both stimulations for NO PSGS (59 and 87 voxels). However, for PSGS, the distribution for both acupuncture and tactile stimulations was negatively skewed with the peak at about ( , uncorrected). Deactivations consistently increased at , uncorrected, , uncorrected and , uncorrected, compared to NO PSGS for both stimulations, increasing 315%, 475%, and 232% for acupuncture and 1608%, 234%, and 130% for tactile stimulation. Activations were reduced by 20% and 34% for acupuncture and 59% and 57% for tactile stimulation at , uncorrected, and , uncorrected. Interestingly, the variation trends for the number of activations at , uncorrected, were different between stimulations, that is, significantly increased (+59%) in acupuncture and significantly reduced (−60%) in tactile stimulation for PSGS. The ratios of activation divided by deactivation ranged from 1441%, 710%, and 598% for NO PSGS to 551%, 99%, and 68% for PSGS in acupuncture and from 2564%, 378%, and 245% for NO PSGS to 61%, 47%, and 46% for PSGS in tactile stimulation.

3.4. Paired -Test between NO PSGS versus PSGS (Figure 4, Table 2)

Figure 4(c) shows the histogram of the paired -test between PSGS and NO PSGS. Most -values were less than zero in both stimulations, indicating the individual-level’s -values of most voxels for PSGS were less than those for NO PSGS. The peak value of the distribution was less for tactile stimulation than for acupuncture. More voxels were shown in tactile stimulation at , uncorrected, or , uncorrected. Significant brain regions of the paired -test between PSGS and NO PSGS are mapped in Figures 4(a) and 4(b). For acupuncture, the BOLD responses of the ipsilateral inferior frontal gyrus, bilateral primary somatosensory cortex (SI), bilateral inferior parietal lobule (IPL), contralateral postcentral gyrus, contralateral insula, and contralateral superior temporal gyrus were significantly reduced in PSGS. For tactile stimulation, more regions’ BOLD responses were significantly reduced for PSGS, including the bilateral inferior frontal gyrus, bilateral middle frontal gyrus, bilateral precentral gyrus, SI, postcentral gyrus, bilateral IPL, bilateral thalamus, bilateral insula, bilateral claustrum, bilateral superior temporal gyrus, bilateral transverse temporal gyrus, and some regions in the cerebellum. Detailed results including coordinates, maximum -values and sizes are summarized in Table 2.

3.5. Changes of One-Sample-Based Group Results between NO PSGS and PSGS (Figures 5 and 6, Figure , and Tables 3 and 4)

The changes in activations and deactivations between PSGS and NO PSGS are mapped in Figures 5 and 6, respectively. The patterns of activations are similar between stimulations. For activations in acupuncture (Figure 5(a)), the dominant results were “stronger,” mainly including the brainstem, inferior frontal gyrus, precentral gyrus, cingulate gyrus, SI, postcentral gyrus and IPL, thalamus, insula, superior temporal gyrus, and declive of the cerebellum. The “weaker,” “disappeared,” and “arisen” regions were few and were located around the “stronger” regions. For activations in tactile stimulation (Figure 5(b)), the dominant results were “weaker” and “disappeared.” These regions mainly included the brainstem, inferior frontal gyrus, middle frontal gyrus, precentral gyrus, SI, postcentral gyrus and IPL, thalamus, insula, putamen, superior temporal gyrus and several cerebellar regions (declive, inferior semilunar lobule, pyramis, tuber, and uvula). All of these regions were bilateral. Very few regions were “stronger.” For deactivations, both stimulations were “arisen” dominant. Furthermore, the patterns of the “arisen” regions were quite similar between stimulations, including the medial frontal gyrus, middle frontal gyrus, superior frontal gyrus, precentral gyrus, SI, superior parietal lobule (SPL), posterior cingulate, precuneus, parahippocampal gyrus, cuneus, and middle occipital gyrus. For tactile stimulation, the culmen of the cerebellum was also “arisen.” Parts of the precentral gyrus, SI, SPL, precuneus, and cuneus were “stronger” in the center of “arisen” regions for tactile stimulation. Detailed results of brain regions in conjunctional maps including BA indices and sizes are summarized in Tables 3 and 4. Figure shows similar results at different thresholds ( , uncorrected, and , uncorrected, resp.).

3.6. Nonspecific and Specific BOLD Responses Evoked by Acupuncture Stimulation (Figure 7, Figure , and Table 5)

The group results for the one-sample -test for NO PSGS in acupuncture are shown in Figure 7(a). Significant activations were present in the bilateral brainstem, the ipsilateral inferior frontal gyrus, the bilateral precentral gyrus, the contralateral cingulate gyrus, the bilateral SI, the contralateral postcentral gyrus, the bilateral IPL, the bilateral thalamus, the bilateral insula, the contralateral claustrum, the bilateral superior temporal gyrus and the contralateral dentate, and the declive, pyramis, and uvula of the cerebellum. The two-sample -test results between acupuncture and tactile stimulation for NO PSGS are shown in Figure 7(b). Compared to tactile stimulation, the BOLD responses of the ipsilateral precentral gyrus, the bilateral SPL, and the bilateral precuneus were more activated in acupuncture stimulation. Detailed results including coordinates, maximum -values, and sizes are summarized in Table 5. Figure shows similar results at other thresholds ( , uncorrected for a one-sample -test and , uncorrected for a two-sample -test).

3.7. De-qi’s Influence on Our Results

Considering that De-qi is associated with a remedial mechanism of acupuncture [3, 28], the relationship between individual responses and De-qi sensations should be elucidated. No direct correlations were shown between the VAS index and the number of activations or deactivations, as well as between the VAS index and the “AVER-WR-RV -value” (all ), which were consistent with our previous study [11].

4. Discussion

The present study first examined the underlying validity of using PSGS in acupuncture fMRI studies. Our results demonstrated that the orthogonality assumption was violated and prominent changes between acupuncture-evoked BOLD responses for NO PSGS and PSGS were shown in most subjects. Secondly, the impact of PSGS on the BOLD responses evoked by acupuncture stimulation was evaluated. For most subjects and the group results, the positive-correlation-dominant BOLD responses were changed into negative-correlation-dominant by PSGS, which was the nonspecific impact of PSGS for both acupuncture and tactile stimulations. In particular, our results implicated that extensive deactivations found in acupuncture fMRI studies may be the effect of global normalization other than a plausible acupuncture-specific effect. In addition, the specific impact of PSGS on acupuncture was that the statistical power of several regions in the somatosensory-related brain network was increased for PSGS, which had the most significant activations for NO PSGS. Finally, our results indicated that the “on-off”-based central BOLD responses during De-qi evoked by acupuncture stimulation were nonspecifically activation-dominant and were mainly around the somatosensory-related brain network.

4.1. Global Normalization Is Unjustified for Exploring BOLD Responses Evoked by Acupuncture Stimulation

Many previous studies demonstrated that adverse consequences were generated by global normalization when the orthogonality assumption was violated [13, 14, 22, 24, 29]. Therefore, researchers suggested that it was necessary to examine the orthogonality assumption before the application of global normalization. First, global normalization should not be applied on the condition that the “GS-RV -value” is large [22]. Second, it should not be applied on the condition that global normalization significantly changes the results [12, 13, 24], even with a small “GS-RV -value” [14]. Our results for acupuncture showed that the larger the “GS-RV -value,” the stronger the relationship of the experimental design and the BOLD signal changes in the whole brain (Figure 2(a)) and the more serious the impact of global normalization (Figure 2(c)). These dependences had a linear relation. The “GS-RV -values” of many subjects were significant (Figure 2, Figure , and Table ). Significant changes between BOLD responses for NO PSGS and PSGS were demonstrated, mainly reversing the activation-dominant trend of BOLD responses to the deactivation-dominant ones (Figures 2 and 3, Table 1, and Tables and ). Therefore, we supposed that it is unjustified for most subjects in acupuncture studies in fMRI. In this sense, we suggest that one should take precautions when interpreting previous conclusions using global normalization.

When the orthogonality assumption was violated, the main consequences caused by PSGS were that (i) the level of activity from regions that were actually associated with the experimental paradigm would be reduced, potentially decreasing statistical power and (ii) the level of activity from regions that were nonassociative with the experimental paradigm might manifest an apparent negative association [12, 14, 2224, 30]. Our results indicated that for subjects whose global signals were positively correlated to the reference vector, the overall level of activations was reduced, whereas deactivations increased (Figures 3 and 4, Table 3, Figure , and Tables and ). On the contrary, for subjects whose global signals were negatively correlated with the reference vector, the overall level of activations was significantly increased, whereas deactivations were significantly reduced by PSGS. On the group level, the improper application of global normalization in our study manifested as reducing activations and introducing artificial deactivations (Figures 2, 3, and 4, Table 1, and Tables and ) according to previous findings [13, 14, 22, 24]. As an inference, if the proportion of negative correlation (positive correlation) in a group was changed, the impact of global normalization should also be changed. Therefore, we inferred that the differences for the application of global normalization, as well as the different proportions of negative correlations in a group when global normalization was applied, may partly explain the disparate results in previous fMRI acupuncture studies. Furthermore, the amplitude of activations was consistently decreased after PSGS (Table ), which was also in accordance with previous reports [22, 30]. We also adopted an additional tactile study to examine the robustness of the current findings. The results drew similar conclusions, suggesting that they were just the fatal influence of a data processing pipeline.

4.2. Extensive Deactivations Introduced by Global Normalization Are Nonspecific Spurious Effects

The role of BOLD deactivation induced by acupuncture is a pivotal but controversial topic [10, 31]. A series of papers believed that the extensive deactivations in the LPNN, highly overlapped with DMN, were specific BOLD responses during De-qi evoked by acupuncture stimulation [3, 5, 19, 21, 3133], which may reflect the neuroregulation mechanism of acupuncture [3, 5, 32]. However, more fMRI-based acupuncture studies failed to discover the significantly extensive deactivations [4, 7, 8, 10, 3442]. Our results demonstrated that the central BOLD responses during De-qi evoked by acupuncture stimulation are activation-dominant (Figures 3 and 7, Table 5, and Table ). For NO PSGS, no deactivations in grey matter were shown at , uncorrected (Figure 7, Table 5), and very few deactivations were shown at , uncorrected (Figure 6 and Figure ). For PSGS, few deactivations were significant at , uncorrected (Figure ), while a great amount of deactivations were significant at , uncorrected (Figure 6). These deactivations were mainly at the occipital cortex and around the DMN. Comparing these results between acupuncture and tactile stimulation, we found that the pattern of deactivations was similar (Figure 6, Figure , and Table 4). These results are in accordance with findings in previous studies [13, 14, 22, 24]. Taken together, we inferred that the significant deactivations in acupuncture were spurious effects due to incorrect global normalization.

The deactivations in DMN were found in many fMRI studies across various fields [4348]. Recent studies about the origins of the deactivations have demonstrated that deactivations in the DMN merely reflected the temporary suspension of activity of the semantic-related network by an externally passive task with reallocation of attentional resources, rather than a specific BOLD response pattern evoked by a task [49, 50]. Interestingly, global normalization often surrounded deactivations in the DMN in fMRI during the early stage [1214, 29, 30], because these deactivations could be enhanced by incorrect global normalization [22]. Reviewing studies about deactivation in acupuncture in fMRI, we found that most of them adopted global normalization, but none of them discussed its impact for introducing deactivations [3, 5, 15, 19, 20]. Consequently, we conjectured that most deactivations reported in previous acupuncture fMRI studies were probably nonspecific manifestations of the DMN under an externally passive task [10]. Moreover, since the BOLD responses evoked by acupuncture stimulation were activation-dominant, the deactivations in the DMN were enhanced after incorrect global normalization.

4.3. Heightened Effects of the Most Significant Activations Produced by Global Normalization Are Specific for Acupuncture

Although most impacts of PSGS are consistent and similar between acupuncture and tactile stimulation, remarkably specific results for acupuncture produced by PSGS were shown, that is, the “stronger” activations in several regions of the somatosensory-related brain network, which were the most significant activations in both acupuncture and tactile stimulation groups for NO PSGS (Figure 6(a) and Table 3). The “stronger” activations in these regions indicated that their group-level statistical power for the acupuncture group was enhanced by PSGS. However, these regions’ activations were “weaker” or “disappeared” in the tactile group (Figure 6(b) and Table 3). For tactile stimulation, all measures indicated that PSGS decreased the activations, which is similar to previous studies [13, 14, 22, 24]. In contrast, the heightened effect of the activations in these regions seemingly contradicted the reduction of the activations on the overall level and the amplitudes of BOLD signal changes for acupuncture. In fact, this heightened effect was also consistent with findings in previous studies [22, 51]. Zarahn et al. proposed that despite the reduced activations’ strength, the inclusion of the global regressor might still affect a reduction in otherwise unexplained variance large enough to cause a net increase in sensitivity [51]. We inferred that the processing of PSGS would more effectively reduce the experimental-paradigm-unrelated variations in the acupuncture group rather than in the tactile group. Furthermore, we supposed that the proportion of experimental-paradigm-unrelated variations of the BOLD signals during acupuncture stimulation might be significantly larger than those during tactile stimulation, which may reflect the complexity of the underlying characteristics of BOLD responses in acupuncture.

Based on the one-sample -test, our results indicate that the “on-off”-modal central BOLD responses during De-qi evoked by acupuncture stimulation were activation-dominant and nonspecific, mainly in the somatosensory, motor and sensory integration brain network (Figure 7, Figure , and Table 5). However, this does not mean that there is no specific effect of acupuncture on the central nervous system. We suggested that exploration of un-“on-off”-modal variations may capture the specific BOLD responses evoked by acupuncture stimulation. For this purpose, the novel experimental design including several conditions in a single run based on better quantification of the acupuncture stimulation and the subjective sensation should be adopted [2, 4]. Besides, the application of data-driven methods may lead to a more comprehensive understanding of the underlying mechanism of acupuncture [5256]. As recommended recently, pattern-based approaches and statistical inferences should be widely adopted in future fMRI-based acupuncture studies [5759].

4.4. Limitations

This paper still has space for improvement in upcoming studies. Firstly, due to the limitation of length, only the method of PSGS was detailed, which was more popular in previous acupuncture studies in fMRI. We also checked the results of global normalization as an additive term in the GLM model in partial datasets and acquired similar results (figures not shown). Secondly, maps and tables under the corrected threshold are not shown due to the unequally absolute -values among groups. The conservative threshold in this paper is , uncorrected, which is parallel to the -values of , FWE corrected. Actually, most results in this paper were gradually changed and became independent of the threshold. Finally, all of the results, conclusions, and inferences limited to the “on-off” experimental modal detections using other paradigms should also be covered.

5. Conclusions

Taken together, the present study demonstrated that the risks and detrimental consequences regarding the standard use of global normalization in fMRI are particularly relevant to acupuncture research. This preprocessing step might be unjustified for acupuncture studies in fMRI. Particularly, extensive deactivations induced by acupuncture might be the nonspecific pernicious consequences of global normalization, which may contribute to conflicting results and interpretations in acupuncture fMRI studies. The central responses evoked by acupuncture stimulation associated with De-qi were nonspecifically activation-dominant and were mainly around the somatosensory-related brain network, whose statistical power was specifically enhanced by PSGS. The specific BOLD responses evoked by acupuncture stimulation may be hidden in the experimental-paradigm-unrelated variations and should be explored by data-driven methods.

Authors’ Contribution

J. Sun and W. Qin contributed equally to this work.

Acknowledgments

This study was supported by the Project for the National Key Basic Research and Development Program (973) under Grants no. 2012CB518501and 2011CB707702, the National Natural Science Foundation of China under Grants no. 30930112, 30970774, 81000640, 81000641, 81101036, 81101108, 31150110171, 30901900, 81271644, and 31200837, the Fundamental Research Funds for the Central Universities, and fund no. 81030027.

Supplementary Materials

Supplementary Material include Supplemental Results, four figures and four tables

  1. Supplementary Material