Compression at Myofascial Trigger Point on Chronic Neck Pain Provides Pain Relief through the Prefrontal Cortex and Autonomic Nervous System: A Pilot Study

Abstract

Compression at myofascial trigger points (MTrPs), known as ‘ischemic compression’, has been reported to provide immediate relief of musculoskeletal pain and reduce the sympathetic activity that exacerbates chronic pain. We conducted a pilot study to investigate the possible involvement of the prefrontal cortex in pain relief obtained by MTrP compression in the present study, we analyzed the relationships among prefrontal hemodynamic activity, activity of the autonomic nervous system, and subjective pain in patients with chronic neck pain, with and without MTrP compression. Twenty-one female subjects with chronic neck pain were randomly assigned to two groups: MTrP compression (n = 11) or Non-MTrP compression (n = 10). Compression for 30 s was conducted 4 times. During the experiment, prefrontal hemodynamic activity [changes in Oxy-hemoglobin (Hb), Deoxy-Hb, and Total-Hb concentrations] and autonomic activity based on heart rate variability were monitored by using near infrared spectroscopy (NIRS) and electrocardiography (ECG), respectively. The results indicated that MTrP compression significantly reduced subjective pain compared with Non-MTrP compression. The spectral frequency-domain analyses of heart rate variability (HRV) indicated that a low frequency (LF) component of HRV was decreased, and a high frequency (HF) component of HRV was increased during MTrP compression, while LF/HF ratio was decreased during MTrP compression. In addition, prefrontal hemodynamic activity was significantly decreased during MTrP compression compared with Non-MTrP compression. Furthermore, changes in autonomic activity were significantly correlated with changes in subjective pain and prefrontal hemodynamic activity. Along with previous studies indicating a role for sympathetic activity in the exacerbation of chronic pain, the present results suggest that MTrP compression in the neck region alters the activity of the autonomic nervous system via the prefrontal cortex to reduce subjective pain.

Introduction

The prevalence of neck pain in the general population has been increasing in recent years (Hoy et al., 2010), and women complain more frequently of neck pain than men (Guez et al., 2002; Côté et al., 2004; Fejer et al., 2006). Musculoskeletal neck pain is reportedly caused by myofascial trigger points (MTrPs) in the neck and shoulder muscles (Simons et al., 1999). MTrPs are hypersensitive spots in palpable taut bands of skeletal muscle fibers (Simons et al., 1999), and recent clinical studies have reported that patients with chronic neck pain have a larger number of MTrPs in the upper trapezius muscle than healthy subjects (Fernández-de-las-Peñas et al., 2007; Muñoz-Muñoz et al., 2012). These studies suggest that MTrPs are responsible for chronic neck pain.

Abnormalities of the autonomic nervous system are thought to be involved in the maintenance of chronic musculoskeletal pain (Passatore and Roatta, 2006; Martinez-Lavin, 2007). Excessive activity of the sympathetic nervous system and reduced activity of the parasympathetic nervous system have been reported in patients with chronic pain such as neck-shoulder pain and fibromyalgia (Gockel et al., 1995; Furlan et al., 2005; Hallman et al., 2011). Altered sympathetic activity with sweating, vasoconstriction, vasodilation, and piloerection has also been reported in MTrP regions (Simons et al., 1999). Increased sympathetic activity is thought to exacerbate spontaneous local pain at MTrPs in patients with chronic neck and shoulder pain (Ge et al., 2006). Thus, abnormalities of the autonomic nervous system seem to be involved in the chronic pain associated with MTrPs.

Recent human neuroimaging studies have shown that several brain regions including the prefrontal cortex are involved in the processing of pain information as well as autonomic regulation (Hallman and Lyskov, 2012). In those with chronic pain, anatomical and functional abnormalities have been found in the medial prefrontal cortex (mPFC) and higher ratings of spontaneous pain were associated with increased activity in the mPFC of patients with chronic pain (Baliki et al., 2006). Furthermore, a study using functional magnetic resonance imaging (fMRI) has suggested that the mPFC is involved in the generation of autonomic responses (Critchley et al., 2000). These results imply that the prefrontal cortex is involved in chronic pain through abnormal activation of the autonomic nervous system.

Compression at MTrPs is an effective massage technique for musculoskeletal pain, and can provide immediate relief from chronic neck pain (Hou et al., 2002), with increased parasympathetic activity and decreased sympathetic activity (Delaney et al., 2002; Takamoto et al., 2009). MTrP compression is considered to exert its therapeutic effects through peripheral, spinal, supraspinal, and autonomic pathways (Bialosky et al., 2009; Takamoto et al., 2009). However, its effects on pain perception and autonomic control in the prefrontal cortex remain unclear.

In the present study, we hypothesized that MTrP compression might affect the mPFC involved in autonomic regulation, which in turn induces pain relief. To investigate possible involvement of the mPFC in pain relief by MTrP compression, we analyzed the relationships among prefrontal hemodynamic activity measured using near infrared spectroscopy (NIRS), activity of the autonomic nervous system based on heart rate variability (HRV), and subjective pain in patients with chronic neck pain, with and without treatment by MTrP compression. NIRS is a non-invasive neuroimaging technique that can measure changes in oxy-hemoglobin (Oxy-Hb), deoxy-hemoglobin (Deoxy-Hb), and total hemoglobin (Total-Hb) in the cortical surface elicited by local cortical neuronal activities (Ferrari and Quaresima, 2012). It is noted that NIRS can be applied with less body and head restriction in a limited space, compared with the other non-invasive methods including fMRI and positron emission tomography (PET). Thus, NIRS allows to measure brain activity in a similar condition to that in actual clinical environments.

Materials and Methods

Subjects

The participants comprised 21 female patients (aged 20–31 years; 23.4 ± 0.9 years, mean age ± standard error [SE]) who had complained of neck pain lasting more than 3 months from the time of onset. All patients were diagnosed with myofascial-pain syndrome by a well-trained clinical practitioner with a national license to practice as an acupuncturist. Inclusion criteria were (1) a palpable band or hardened nodules in the upper trapezius muscle, and (2) myofascial pain emanating from a well-localized area in the palpable band. These patients were randomly assigned to two groups (MTrP and Non-MTrP groups) (Figure 1), based on a random allocation software. All patients were treated in strict compliance with the Declaration of Helsinki and the U.S. code of Federal Regulations for the Protection of Human Subjects. The experiments were conducted with the understanding and informed written consent of each subject, and approved by the Clinical Research and Ethics Committee at the University of Toyama.

FIGURE 1

Figure 1. Participant flow diagram.

Muscle Compression

Ischemic compression consisted of constant deep pressure by the thumb, which was applied to MTrPs or Non-MTrPs in the trapezius muscle. The locations of the MTrPs in the trapezius muscle were identified by a practitioner with a national license to practice acupuncture in Japan who had more than 6 years' experience in MTrPs diagnosis and treatment. The identification of MTrPs was based on the minimal diagnostic criteria recommended by Simons et al. (1999) and Gerwin et al. (1997), requiring (1) the presence of a hypersensitive tender point in a taut band in the muscle, (2) recognition by the patient of the pain as “familiar” when compression is applied to the point, and (3) induction of pain when the muscle including the MTrP is stretched. For the controls, Non-MTrPs were defined as points 2 cm away and proximal to a given MTrP, where a taut band was not detected and local and referred pain was not induced when compression was applied to the point.

The intensity of compression was set at a point midway between the pressure pain threshold (PPT) and maximally tolerable pain (MTP), which corresponded to a “moderately painful average” intensity in each patient (Hou et al., 2002). PPTs and MTPs were measured using a digital algometer (Atkins et al., 1992; Bendtsen et al., 1994). The device consisted of a capacitance pressure sensor with a circular tip (diameter, 6 mm) (HV-2-556: SHOWA SOKKI, Japan) and data collection hardware connected to the sensor (SME-101A: Kyowa Dengyou, Japan). The sensor was attached to the pad of the thumb. Pressure intensity was expressed in Arbitrary Unit (AU); 550 AU is equivalent to about 1 kg. Pressure intensity was displayed on a device monitor that was not visible to the palpating practitioner or the patients. Two experimenters worked together during data collection. Experimenter 1 (palpating practitioner), who was blinded to the monitor recordings, determined the PPT and MTP based on the behavioral responses of the patient, while Experimenter 2 recorded the pressure intensity values from the device monitor. Experimenter 1 applied steady, gradually increasing pressure to the identified MTrPs or Non-MTrPs using the thumb with the pressure sensor. The patients were instructed to press a button when they initially felt pain and also when they felt severe pain, where 0 = no pain, 4 = pain threshold (PPT), and 10 = severe pain. Compression ceased when patients felt severe pain, and Experimenter 2 recorded the pressure intensity values. The PPT and MTP were recorded 3 times for each stimulation site, and the intensity of compression in the later experiment was set at an intensity midway between the mean PPT and mean MTP for each subject.

Protocol

The patients were fitted with a NIRS head pad on the frontal region (Figure 2A) and ECG electrodes on the chest. The patients lay in a supine position on a bed with the cervical spine in a neutral position. The patients were instructed to close their eyes and relax as much as possible. Muscle compression was started 1 min after the commencement of NIRS and ECG recordings. Compression was applied to MTrPs or Non-MTrPs for 30 s, followed by a rest for 120 s. This cycle was repeated 4 times (4 cycles). The practitioner maintained constant pressure (i.e., a pressure intensity midway between the PPT and MTP) by monitoring the pressure intensity of the thumb.

FIGURE 2

Figure 2. Location of NIRS probes and channels. (A) Location of the NIRS head pad on the frontal region. (B) The arrangement of probes (sources and detectors) and channels. (C) Spatial estimation of NIRS channels. The coordinates of each NIRS channel were normalized to MNI (Montreal Neurological Institute) space using virtual registration. The red dots indicate the coordinates of the channels: a, right dorsolateral prefrontal cortex (rDLPFC) (Brodmann Areas 45, 46); b, right dorsomedial prefrontal cortex (rDMPFC) (Brodmann Area 10); c, central dorsomedial prefrontal cortex (cDMPFC) (Area 10); d, left dorsomedial prefrontal cortex (lDMPFC) (Area 10); e, left dorsolateral prefrontal cortex (lDLPFC) (Areas 45, 46).

Psychophysical Data

To subjectively evaluate the compression treatment, the subjects were instructed to report pain levels in the neck before and after the experiment using a 100-grade visual analog scale (VAS) that is measured in mm. In subjective pain scoring, score 0 indicted “no pain,” and score 100, “strongest pain perception that the subjects have ever had.” At the end of the experiment, the subjects were also asked to report the pleasantness of the compression in each cycle on a scale of –10 to 10, where the extremes represented the most unpleasant and pleasant experiences, respectively (the comfort/discomfort score). The comfort/discomfort scores from 4 cycles were averaged. Furthermore, the subjects reported the pain intensity during compression in each cycle using a 10-grade scale, where a score of 0 indicated “no pain,” a score of 4 indicated “pain sensation at a threshold level,” and score 10 indicated “strongest pain perception” (pain intensity score). The pain intensity scores from 4 cycles were averaged.

Measurement of the Activity of the Autonomic Nervous System

Autonomic activity was measured using a system for measuring HRV (Makin AD2; Suwa Trust, Japan). The V5 chest leads of the ECG were used. RR (beat-to-beat) intervals were analyzed by the maximum-entropy method to calculate the HRV spectrum. The low-frequency (LF) (0.04–0.15 Hz) and high-frequency (HF) (0.15–0.40 Hz) components of heart variability (HRV) were calculated. In addition, the LF/HF ratio was obtained. The LF/HF ratio has been used as an index of a sympatho-vagal balance, the LF component as an index of sympathetic modulation, and the HF component as an index of parasympathetic (vagal) modulation (Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology, 1996; Pagani et al., 1997; Thayer et al., 2010). However, recent studies indicate that meanings of these indices are more complex (see Limitations) (Billman, 2011, 2013). Furthermore, these indices are also modulated by several non-autonomic nervous factors such as respiration and heart rate themselves (Billman, 2013; Gasior et al., 2016). Therefore, respiratory rates were estimated from the R-R intervals of ECG data using a free software (Kubios HRV ver.2.2; http://www.kubios.com; University of Eastern Finland, Kuopio, Finland; Tarvainen et al., 2014).

Measurement of Cerebral Hemodynamic Activity in the Prefrontal Cortex

A NIRS instrument (Spectratech OEG-16; Spectratech Inc., Yokohama, Japan) was used to measure cerebral hemodynamic activity. Six optical sources and 6 detectors (resulting in a total of 16 recording channels) were fixed on the frontal area (Figures 2A,B). The bottom horizontal line of the probes was placed on the Fp1–Fp2 line in the 10–20 EEG system. The sources and detectors were positioned so that distance between any 2 probes was set at 3 cm. The midpoints between a source and a detector were called “channels,” and hemodynamic activity was detected by pairs of source and detector using two different wavelengths (770 and 840 nm). Changes in the hemoglobin (Hb) concentration (ΔOxy-Hb, ΔDeoxy-Hb, and ΔTotal-Hb [ΔOxy-Hb + ΔDeoxy-Hb]) from the baseline were estimated using a modified Lambert-Beer law (Seiyama et al., 1988; Wray et al., 1988). After the recording, the 3-D locations of the NIRS probes were measured using a digitizer (Shimadzu Co. Ltd., Japan) with reference to the nasion and bilateral external auditory meatus.

Data Analysis

The data are presented as the mean ± standard error (SE). Normality of the data distribution was assessed by Shapiro-Wilk test. The homogeneity of variance in all variables with normality was assessed by the Levene's test. Student's t-test (or Mann-Whitney U test) and ANOVA were used to compare the data in analyzed parameters between the MTrP and Non-MTrP groups. When the data during compression were compared between the MTrP and Non-MTrP groups, the data (VAS, heart rate, respiratory rate, and HRV parameters) during compression were corrected for the baselines.

The HF and LF components were reported as percentages (HF%, LF%, respectively), and normalized by dividing them by the sum of all components (LF + HF). This normalization of LF and HF is recommended since it tends to minimize the effect on the values of LF and HF components of the changes in total power (Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology, 1996; Pagani et al., 1997). Furthermore, the low values of HF, LF, and LF/HF were also normalized by natural logarithm transformation (lnLF, lnHF, lnLF/HF, respectively). Changes in HRV after the treatments were assessed by subtracting the mean value of each parameter in the 30-s rest period before the application of compression from the mean value for the 30-s period during compression. All autonomic data (changes in HR, respiratory rate, and HRV parameters) were averaged across 4 cycles, since the all autonomic parameters except respiratory rate showed no significant main effect of “cycle” (P > 0.05) by repeated measures one-way ANOVAs with a factor of “cycle” in each treatment group while respiratory data showed no significant effect of “cycle” (P > 0.05) by Friedman's test. Student's t-test was used to compare changes in autonomic activity between the MTrP and Non-MTrP groups. Findings of P < 0.05 were considered significant.

To identity the anatomical locations of NIRS channels in each subject, the 3D locations of the NIRS probes and channels in each subject were spatially normalized to a standard coordinate system using NIRS SPM software (statistical parametric mapping: http://bisp.kaist.ac.kr/NIRS-SPM, version 3.1; Ye et al., 2009); the coordinates for each NIRS channel were normalized to MNI (Montreal Neurological Institute) space using virtual registration (Tsuzuki et al., 2007). We then identified the Brodmann areas corresponding to the NIRS channels of each subject using MRIcro software (www.mricro.com, version 1.4). We divided the channels into 5 regions: the dorsolateral prefrontal cortex (DLPFC; Brodmann areas 45 and 46) in each hemisphere corresponding to Ca (rDLPFC) and Ce (lDLPFC) in Figures 2, 3 regions in the dorsomedial prefrontal cortex (DMPFC; Brodmann area 10) corresponding to Cb (rDMPFC), Cc (cDMPFC), and Cd (lDMPFC) in Figure 2. In two subjects, the locations of NIRS channels were determined by stereotaxic superimposition on the surface of the 3-D MRI reconstructed brain of each subject. For 3-D MRI, thin-slice 3-D sagittal T1-weighted gradient echo MR images were obtained at 1.5 T using a specific protocol tailored for reconstruction (Takeuchi et al., 2009). These two subjects had the following protocol: (TR/TE/NSA) 25/5/1, flip angle 10, FOV 87.5 cm, matrix 256 9 256, 1.0 mm contiguous slices, obtained in a plane parallel to the brain stem.

FIGURE 3

Figure 3. Comparison of changes in subjective pain scores after compression between the MTrPs and Non-MTrP groups. Ischemic compression at MTrPs significantly decreased subjective pain scores compared with Non-MTrP compression. Error bars indicate SEM. VAS, visual analog scale; **P < 0.01.

Cerebral hemodynamic changes during compression at MTrPs and Non-MTrPs were converted to effect sizes. Effect sizes can adjust for the influence of different path-length factors among different subjects and cortical regions (Schroeter et al., 2003; Suzuki et al., 2008). The effect sizes of hemodynamic responses were calculated according to the following formula: effect size = [(mean Oxy-Hb levels during compression for 30 s) − (mean Oxy-Hb levels during the rest period for 30 s before start of compression)]/[standard deviation of Oxy-Hb levels during the rest period of 30 s before the start of compression]. For each channel, the effect sizes from 4 cycles were averaged. The data from the channels in each brain region were then averaged to give the mean hemodynamic responses in each area for each patient. These data for compression at MTrPs and Non-MTrPs were compared using a repeated-measures two-way ANOVA (treatment × brain region). Findings of P < 0.05 were considered significant.

We also analyzed correlations among the effect sizes of hemodynamic responses in the DMPFC, changes in possible autonomic activity with compression (HF%, LF%, LF/HF, lnHF, lnLF, and lnLF/HF), and changes in subjective pain scores with compression (VAS scores) using simple regression analysis.

All data analyses were performed using SPSS 19.0 (IBM Inc., New York, USA). A p < 0.05 was considered statistically significant.

Sample Size

Sample size for comparison of two independent samples (two-tailed t-test) was estimated using a free sample size calculator (https://www.stat.ubc.ca/~rollin/stats/ssize/n2.html by Dr. Rollin Brant, University of British Columbia) as n = 10 based on the following condition; changes in VAS in the MTrP group = −35, changes in VAS in the Non-MTrP group = −16, standard deviation (SD) = 15, level of significance = 0.05, statistical power = 0.8. Our previous preliminary data were used for this sample size estimation.

Results

Baseline Characteristics and Sensations Evoked by Compression in the Two Groups

Table 1 shows the baseline clinical characteristics of the patients (age, VAS score for neck pain, PPT, and intensity of compression, heart rate, and respiratory rate). There were no significant differences between the two groups in the baseline characteristics of these parameters (Student's t-test, P > 0.05). Furthermore, there were no significant differences in the sensations evoked by compression between the two groups (pain intensity scores, comfort/discomfort scores) (Student's t-test, P > 0.05).

TABLE 1

Table 1. Baseline characteristics and sensations evoked by compression in the two groups.

Changes in Subjective Pain Scores in the Neck

Figure 3 compares the changes in subjective pain scores resulting from compression in the MTrP and Non-MTrP groups. Compression at MTrPs significantly ameliorated subjective pain scores compared with compression at Non-MTrPs (Student's t-test, P < 0.01).

Changes in Heart Rate Variability from the Pre-treatment Baseline

There was no significant difference in changes in heart rates during compression between the MTrP and Non-MTrP groups (Student's t-test, P > 0.05) (Supplementary Figure 1A), nor significant difference in changes in respiratory rates during compression between the MTrP and Non-MTrP groups (Mann-Whitney U test, P > 0.05) (Supplementary Figure 1B). Figure 4 compares changes in autonomic responses during compression between the MTrP and Non-MTrP groups. The HF components (HF%) were significantly greater in the MTrP group compared with the Non-MTrP group (Student's t-test, P < 0.01) (A), while the LF components (LF%) were significantly greater in the MTrP group compared with the Non-MTrP group (Student's t-test, P < 0.01) (B). The LF/HF ratio was significantly less in the MTrP group compared with the Non-MTrP group (Student's t-test, P < 0.01) (C). These results suggest that compression at MTrPs suppressed sympathetic activity.

FIGURE 4

Figure 4. Comparison of changes in the activity of the autonomic nervous system during ischemic compression between the MTrP and Non-MTrP groups. (A) High frequency (HF) components of HR variability (HF%) were significantly increased during ischemic compression at MTrPs compared with Non-MTrPs. (B) Low frequency (LF) components of HR variability (LF%) was significantly decreased during ischemic compression at MTrPs compared with Non-MTrPs. (C) LF/HF ratio was significantly decreased during ischemic compression at MTrPs compared with Non-MTrPs. Error bars indicate SE; **P < 0.01.

We also analyzed logarithmically transformed HRV parameters (Supplementary Figures 1C–E). The trend of the results was essentially consistent with that in Figure 4. Changes in lnLF during compression were significantly less in the MTrP group compared with the Non-MTrP group (Student's t-test, P < 0.01) (C), while there were no significant differences in changes in lnHF during compression between the MTrPs and Non-MTrPs (Student's t-test, P > 0.05) (D). Furthermore, changes in lnLF/HF during compression were significantly less in the MTrP group compared with the Non-MTrP group (Student's t-test, P < 0.01) (E).

Prefrontal Hemodynamic Responses

Figure 5 shows typical examples of Oxy-Hb maps 25 s after starting compression at Non-MTrPs (A) and MTrPs (B). The NIRS data were stereotaxically superimposed on the 3D-MRIs of the brain of the subject to construct the Oxy-Hb concentration maps. The Oxy-Hb concentration in the DMPFC gradually increased during compression at Non-MTrPs (A), while it decreased during compression at MTrPs (B).

FIGURE 5

Figure 5. Time course and topographical maps of Oxy-Hb changes during Non-MTrP (A) and MTrP (B) compression. Topographical maps indicate Oxy-Hb maps 25 s after starting ischemic compression. NIRS image data were superimposed on individual 3D-MRIs. Yellow dots on the 3D-MRI indicate NIRS channels. Red lines indicate changes in Oxy-Hb concentrations from baseline activity. Note that the Oxy-Hb concentration in the prefrontal cortex was decreased during MTrP compression.

Figure 6 compares the mean effect sizes of hemodynamic responses in the 5 prefrontal regions of the MTrP and Non-MTrP groups. Statistical analysis of the data using a repeated-measures 2-way ANOVA with “treatment” (MTrP vs. Non-MTrP) and “brain region” as factors indicated that there was a significant main effect of treatment [F(1, 19) = 6.624, P < 0.05], but not brain region [F(2.697, 51.236) = 1.672, P > 0.05], and no significant interaction between treatment and brain region [F(2.697, 51.236) = 1.179, P > 0.05]. The results indicate that compression at MTrPs significantly decreased hemodynamic activity in the prefrontal cortex. In contrast, there was no significant differences in Deoxy-Hb between the MTrP and Non-MTrP group in the same statistical analysis using a repeated measures 2-way ANOVA; there was no significant main effect of treatment [F(1, 19) = 0.019, P > 0.05], nor significant interaction between treatment and brain region [F(2.8, 53.73) = 0.664, P > 0.05].

FIGURE 6

Figure 6. Comparison of hemodynamic responses in the prefrontal regions during ischemic compression between MTrPs and Non-MTrP groups. Hemodynamic responses are shown as effect sizes. Statistical analysis by repeated-measures 2-way ANOVA indicated that MTrP compression significantly decreased the effect size compared with Non-MTrP compression. Error bars indicate SE. rDLPFC, right dorsolateral prefrontal cortex; rDMPFC, right dorsomedial prefrontal cortex; cDMPFC, central dorsomedial prefrontal cortex; lDMPFC, left dorsomedial prefrontal cortex; lDLPFC, left dorsolateral prefrontal cortex.

Relationships among Subjective Pain, Autonomic Activity, and Prefrontal Hemodynamic Responses

Figures 7A–C show the correlations between changes in autonomic activity with compression and changes in subjective pain scores. Changes in HF% were significantly and negatively correlated with changes in subjective pain scores, r² = 0.272, F(1, 20) = 7.092, P < 0.05 (A). In contrast, changes in LF% were significantly and positively correlated with changes in subjective pain scores, r² = 0.272, F(1, 20) = 7.092, P < 0.05 (B). Furthermore, changes in LF/HF were significantly and positively correlated with changes in subjective pain scores, r² = 0.285, F(1, 20) = 7.573, P < 0.05 (C).

FIGURE 7

Figure 7. Correlations between changes in autonomic activity and changes in subjective pain scores (A–C), and between changes in autonomic activity and changes in hemodynamic responses in the DMPFC (D–F). (A) A negative correlation was observed between changes in HF% and changes in subjective pain scores. VAS, visual analog scale. (B) A positive correlation was observed between changes in LF% and changes in subjective pain scores. (C) A positive correlation was observed between changes in LF/HF ratios and changes in subjective pain scores. (D) A negative correlation was observed between changes in HF% and changes in the hemodynamic response in the DMPFC. (E) A positive correlation was observed between changes in LF% and changes in the hemodynamic response in the DMPFC. (F) A positive correlation was observed between changes in LF/HF ratios and changes in the hemodynamic response in the DMPFC.

Figures 7D–F show correlations between changes in autonomic activity and the mean effect sizes of cerebral hemodynamic responses in the cDMPFC during compression. Changes in HF% were significantly and negatively correlated with the effect sizes of hemodynamic responses in the DMPFC, r² = 0.235, F(1, 20) = 5.830, P < 0.05 (D). In contrast, changes in the LF% were significantly and positively correlated with the effect size of hemodynamic responses in the DMPFC, r² = 0.235, F(1, 20) = 5.830, P < 0.05 (E). Furthermore, changes in the LF/HF ratio were significantly and positively correlated with the effect size of hemodynamic responses in the DMPFC, r² = 0.192, F(1, 20) = 4.514, P < 0.05 (F). However, there was no significant correlation between changes in subjective pain scores and the effect size of hemodynamic responses in the DMPFC (data not shown).

We also analyzed logarithmically transformed parameters (lnLF, lnHF, and lnLF/HF) (Supplementary Figure 2). The results were essentially consistent with those in Figure 7. Changes in lnLF were significantly and positively correlated with changes in subjective pain scores, r² = 0.299, F(1, 20) = 8.121, P < 0.05 (A). Furthermore, changes in lnLF/HF were significantly and positively correlated with changes in subjective pain scores, r² = 0.332, F(1, 20) = 9.439, P < 0.05 (B). However, there was no significant correlation between changes in lnHF and changes in subjective pain scores (data not shown). On the other hand, changes in lnHF were significantly and negatively correlated with the effect size of hemodynamic responses in the DMPFC, r² = 0.205, F(1, 20) = 4.890, P < 0.05 (C). Furthermore, changes in lnLF/HF were significantly and positively correlated with the effect size of hemodynamic responses in the DMPFC, r² = 0.247, F(1, 20) = 6.223, P < 0.05 (D). However, there was no significant correlation between changes in lnLF and the effect size of hemodynamic responses in the DMPFC (data not shown).

Discussion

Effects of MTrP Compression on Pain Perception and the Activity of the Autonomic Nervous System

In the present study, compression at MTrPs significantly ameliorated subjective pain compared with compression at Non-MTrPs (Figure 3). It also significantly increased HRV parameters, which are believed to reflect parasympathetic activity, and inhibited HRV parameters, which are believed to reflect sympathetic activity (Figure 4, Supplementary Figure 1). Furthermore, these changes in autonomic activity were significantly correlated with changes in subjective pain scores (Figure 7, Supplementary Figure 2). Consistent with the present study, previous physiological studies have reported that compression at MTrPs increases the activity of the parasympathetic nervous system and decreases sympathetic activity (Delaney et al., 2002; Takamoto et al., 2009). Furthermore, a reduction of sympathetic activity during thermotherapy for the neck has been associated with decreased neck stiffness and fatigue (Yasui et al., 2010). These findings suggest that altered sympathetic activity might exacerbate neck pain, and that compression at MTrPs might reduce neck pain by suppressing sympathetic activity.

Acute or chronic muscle overload induces hyperactivation of the neuromuscular junction, which leads to an excessive release of acetylcholine (Gerwin et al., 2004; Bron and Dommerholt, 2012). This excess of acetylcholine at the motor endplates leads to the formation of knots in muscle fibers that are characterized by continuous contraction. Formation of contraction knots then leads to the development of local ischemia and hypoxia. The loss of energy and O2 supply in contracting muscle causes the release of sensitizing noxious substances that lead to increased hypersensitivity and pain. Increased activity of the sympathetic nervous system might exacerbate these processes and cause MTrP formation. The sympathetic nervous system controls extrafusal and intrafusal muscle fibers through collateral branches (Selkowitz, 1992; Bombardi et al., 2006), and enhances the release of acetylcholine from motor nerve terminals that is mediated by α- and β-adrenoceptors on motor nerve terminals (Gerwin et al., 2004).

By contrast, a reduction in sympathetic activity by MTrP compression might suppress MTrP processes and the associated pain. Compression at MTrPs in the upper trapezius muscle has been shown to ameliorate the perception of pain and reduce spontaneous electrical activity (SEA) in motor end-plates at MTrPs (Kostopoulos et al., 2008). SEA is thought to be induced by excessive acetylcholine release (Gerwin et al., 2004). Pain at MTrPs and SEA is also decreased by the injection of sympatholytic agents, such as α1-adrenergic antagonists (Hubbard and Berkoff, 1993; McNulty et al., 1994; Hong and Simons, 1998) and relaxation techniques to reduce sympathetic activity (Banks et al., 1998). These findings suggest that suppression of sympathetic activity by MTrP compression might reduce acetylcholine release and decrease muscle contraction.

Furthermore, reduced local blood flow has been observed in trapezius myalgia and was correlated with pain intensity (Larsson et al., 1990), which might be ascribed to vasoconstriction due to increased sympathetic activity. Compression at MTrPs can induce reactive hyperemia in the MTrP region (Simons et al., 1999). Taken together, these results suggest that compression at MTrPs induces pain relief through inhibition of sympathetic activity, which (1) might increase the peripheral blood flow and subsequent removal of noxious substances, and (2) might block the excessive release of acetylcholine. However, other alternative hypotheses for generation of MTrPs such as the neurogenic hypothesis (Srbely, 2010) and the neurophysiologic hypothesis (Partanen et al., 2010) have been also proposed. Further studies are required to fully understand the role of sympathetic activity in myofascial pain syndrome.

Effects of MTrP Compression on Prefrontal Hemodynamic Responses

In the present study, the Oxy-Hb concentration was significantly decreased in the PFC including the DMPFC during compression at MTrPs compared with Non-MTrPs, while there was no significant difference in Deoxy-Hb between the MTrP and Non-MTrP groups. Thus, the present results did not show typical hemodynamic changes (i.e., the response patterns of Oxy-Hb were not opposite to those of Deoxy-Hb concentration in the present study). However, changes in Deoxy-Hb were reported to be not consistent across individuals and across tasks (Hoshi et al., 2001; Toichi et al., 2004; Sato et al., 2005), while the Oxy-Hb concentration was stronger correlated with fMRI BOLD signals than Deoxy-Hb concentration (Strangman et al., 2002; Yamamoto and Kato, 2002). These findings suggest that Oxy-Hb concentration may be the most consistent parameter for cortical activity (Okamoto et al., 2006).

Decreases in fMRI signals were associated with local decreases in neuronal activity (Shmuel et al., 2006). Concentration of inhibitory neurotransmitter GABA was associated with decreases in signal intensity in fMRI (Northoff et al., 2007; Muthukumaraswamy et al., 2009). In NIRS studies, a decrease in Oxy-Hb concentration might correspond to inhibition of brain activity (Seitz and Roland, 1992; Shmuel et al., 2002; Stefanovic et al., 2004). Thus, the decrease in Oxy-Hb suggests that prefrontal activity was suppressed by compression at MTrPs. However, it is noted that compression at Non-MTrPs increased the Oxy-Hb concentration in the prefrontal cortex. The opposite hemodynamic responses in MTrP and Non-MTrP compression cannot be ascribed to differences in the psychophysical characteristics of compression between MTrPs and Non-MTrPs; there were no significant differences in compression intensity, pain intensity scores, and comfort/discomfort scores between MTrP and Non-MTrP compression. Further studies are required to determine the physiological factors contributing to differences in prefrontal hemodynamic responses between compression at MTrPs and Non-MTrPs.

Furthermore, changes in Oxy-Hb concentration in the DMPFC were significantly correlated with changes in HRV parameters (HF%, LF%, LF/HF; lnHF, lnLF, lnLF/HF) during compression (Figure 7, Supplementary Figure 2). A similar correlation has been reported for thermotherapy on the neck region (Yasui et al., 2010). Previous neuroanatomical and non-invasive imaging studies have reported that the PFC including the DMPFC has direct connections with the hypothalamus and brainstem, which are involved in autonomic and behavioral responses to pain (Ongür et al., 1998; Hadjipavlou et al., 2006). Increased activity in the DMPFC during the experience of pain has been correlated with decreased skin blood flow and increased skin conductance responses, suggesting increased sympathetic activity (Seifert et al., 2013). Furthermore, EEG gamma-band oscillation in the DMPFC and autonomic functions coherently changed in response to mental stress, and an increase in gamma-band oscillation went ahead of the autonomic fluctuation (Umeno et al., 2003). Chronic pain states such as chronic low back pain and sympathetically mediated chronic pain were found to be associated with hyperactivity in the mPFC including the DMPFC (Apkarian et al., 2001; Baliki et al., 2008). On the other hand, the dorsolateral PFC is implicated in inhibition of pain perception (Lorenz et al., 2003; Brighina et al., 2004), and gray matter atrophy in this brain region was reported in patients with chronic neck pain with MTrP and chronic low back pain (Fritz et al., 2016; Niddam et al., 2017). These findings suggest that unbalanced activity within the PFC might be associated with chronic pain. Thus, the available findings suggest that hyperactivity in the mPFC including the DMPFC might be involved in altered autonomic activity in the chronic pain state, and that compression at MTrPs might suppress sympathetic activity via the mPFC.

Possible Physiological Mechanisms for the Effects of MTrP Compression

Different cerebral hemodynamic and autonomic nervous system responses were observed for compression at MTrPs and Non-MTrPs in the present study. The spatial locations of 70% of MTrPs have been reported to overlap with the locations of traditional acupuncture points (Melzack et al., 1997). Needling stimulation at both acupuncture points and MTrPs has been shown to significantly induce a stronger specific sensation known as “deqi” compared with needling stimulation at Non-MTrPs and non-acupuncture points (Roth et al., 1997; Takamoto et al., 2010). Deqi sensation is a composite of unique sensations described as aching, soreness, pressure, heaviness, fullness, warmth, cooling, numbness, tingling, and dull pain (Kong et al., 2007). Deqi sensations have also been induced by pressure stimulation at acupuncture points (Yip and Tse, 2004, 2006; Li et al., 2007). Furthermore, sustained pressure stimulation at MTrPs has also elicited deqi-like sensations including aching and dull pain in the distal or proximal areas from the point of stimulation (Simons, 1995; Clark, 2008; Delany, 2014). Thus, stimulation at both MTrPs and acupoints induces similar deqi sensations. MTrPs and acupuncture points are considered to be locations where the sensitization of polymodal-type receptors occurs (Kawakita and Itoh, 2002). Polymodal-type receptors are present on afferent fibers with low conduction velocity (A-delta and C-fibers) (Almeida et al., 2004), which might be involved in deqi sensations evoked by needling stimulation at acupuncture points (Lu, 1983; Wang et al., 1985). The physiological effects induced by compression at MTrPs in the present study might be attributed to deqi sensations.

Needling at acupuncture points that induced deqi sensations has been shown to increase local blood flow, whereas the stimulation of non-acupuncture points changes the local blood flow only slightly (Kuo et al., 2004). Furthermore, the number of episodes of deqi sensations induced by needling at acupuncture points was correlated with a decrease in sympathetic activity and an increase in parasympathetic activity (Sakai et al., 2007), and a significant negative correlation was observed between the intensity of the specific acupuncture sensations and heart rate responses (Beissner et al., 2012). A neuroimaging study reported that decreased activity in the prefrontal cortex was associated with both increased intensity of deqi sensation and heart rate responses during needling at acupuncture points (Beissner et al., 2012). Furthermore, heart rate responses were correlated with fMRI BOLD signals in the mPFC during needling at acupuncture points, and a greater decrease in mPFC activity was associated with a greater decrease in heart rates (Napadow et al., 2013). A previous NIRS study also showed that needling at MTrPs, which induced deqi sensations, decreased Oxy-Hb concentrations in the DMPFC and supplementary motor cortex (Takamoto et al., 2010). Taken together, different therapeutic methods employing either compression or needling at MTrPs and/or acupuncture points, might provide pain relief through common afferent nerve fibers that induce deqi sensations. Further studies are required to investigate the physiological mechanisms that induce cerebral hemodynamic and autonomic nervous system responses by compression at MTrPs.

Conclusions

We conducted a pilot study to investigate the effects of compression at MTrPs in the neck region of patients with chronic neck pain on subjective pain perception, prefrontal hemodynamic activity, and possible autonomic activity using NIRS and HRV analyses. Compression at MTrPs significantly improved subjective pain scores compared with compression at Non-MTrPs. HRV parameters that are believed to reflect parasympathetic activity were significantly increased during compression at MTrPs compared with that induced by compression at Non-MTrP, and HRV parameters that are believed to reflect sympathetic activity were decreased. Furthermore, compression at MTrPs significantly decreased the Oxy-Hb concentration in the DMPFC compared with that induced by Non-MTrP compression. Changes in HRV parameters that are believed to reflect sympathetic activity were positively correlated with changes in Oxy-Hb concentrations in the DMPFC, and were positively correlated with changes in subjective pain scores during ischemic compression. The present results, along with those of previous studies, suggest that the effects of compression at MTrPs might be mediated through inhibitory effects on DMPFC activity, which might be beneficial for treating chronic pain in which hyperactivity of the sympathetic nervous system is involved.

Limitations

First of all, we analyzed non-invasively HRV, but not by direct microneurography of sympathetic and parasympathetic nervous activity. This non-invasive method is widely used in clinical and experimental researches (Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology, 1996; Lombardi and Stein, 2011), and normalized HRV parameters (HF%, LH%, LH/HF) were effective to evaluate muscle sympathetic nervous activity (Pagani et al., 1997). However, recent reviews argue that meaning of these HRV indices is more complex in that both LF and HF components may be modulated not only by sympathetic and parasympathetic neural activity but also by other non-autonomic nervous factors (Lombardi and Stein, 2011; Billman, 2013; Heathers, 2014; Gasior et al., 2016), which should be taken into account.

Second, it is noted that Oxy-Hb signals includes both cerebral (brain) and extra cerebral (scalp, skull, cerebrospinal fluid) components (Fukui et al., 2003). Our previous study using a multi-distance NIRS probe arrangement indicated that mean contribution ratio of the brain to whole hemodynamic signals at 3 cm of probe distance was 66.6% in the DMPFC (Ishikuro et al., 2014). Further studies were required to determine contribution ratio of the brain during compression at MTrPs to dissociate extra cerebral and cerebral components of the hemodynamic responses.

Finally, some of the correlations among autonomic nervous activity, neck pain and PFC hemodynamic activity were not high with a substantial variance in the data in the present study. Furthermore, causal relationships among these parameters remain unknown. Further studies with larger sample size, in which inhibition of the DMPFC using transcranial direct current stimulation (tDCS) should be also tested, are required.

Author Contributions

Designing work: HNishijo and KT; Data acquisition: YM and KT; Data analysis and interpretation: All of the authors.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This work was supported partly by a research grant from Japan Judo Therapist Association and Grants-in-Aid for Scientific Research (C) (16K01499).

Supplementary Material

The Supplementary Material for this article can be found online at: https://[see PDF]/article/10.3389/fnins.2017.00186/full#supplementary-material

References

Almeida, T. F., Roizenblatt, S., and Tufik, S. (2004). Afferent pain pathways: a neuroanatomical review. Brain. Res. 1000, 40–56. doi: 10.1016/j.brainres.2003.10.073

Apkarian, A. V., Thomas, P. S., Krauss, B. R., and Szeverenyi, N. M. (2001). Prefrontal cortical hyperactivity in patients with sympathetically mediated chronic pain. Neurosci. Lett. 311, 193–197. doi: 10.1016/S0304-3940(01)02122-X

Atkins, C. J., Zielinski, A., Klinkhoff, A. V., Chalmers, A., Wade, J., Williams, D., et al. (1992). An electronic method for measuring joint tenderness in rheumatoid arthritis. Arthritis. Rheum. 35, 407–410. doi: 10.1002/art.1780350408

Baliki, M. N., Chialvo, D. R., Geha, P. Y., Levy, R. M., Harden, R. N., Parrish, T. B., et al. (2006). Chronic pain and the emotional brain: specific brain activity associated with spontaneous fluctuations of intensity of chronic back pain. J. Neurosci. 26, 12165–12173. doi: 10.1523/JNEUROSCI.3576-06.2006

Baliki, M. N., Geha, P. Y., Jabakhanji, R., Harden, N., Schnitzer, T. J., and Apkarian, A. V. (2008). A preliminary fMRI study of analgesic treatment in chronic back pain and knee osteoarthritis. Molecular. Pain. 4:47. doi: 10.1186/1744-8069-4-47

Banks, S. L., Jacobs, D. W., Gevirtz, R., and Hubbard, D. R. (1998). Effects of autogenic relaxation training on electromyographic activity in active myofascial trigger points. J. Musculoskelet. Pain. 6, 23–32. doi: 10.1300/J094v06n04_03

Beissner, F., Deichmann, R., Henke, C., and Bär, K. J. (2012). Acupuncture–deep pain with an autonomic dimension? Neuroimage 60, 653–660. doi: 10.1016/j.neuroimage.2011.12.045

Bendtsen, L., Jensen, R., Jensen, N. K., and Olesen, J. (1994). Muscle palpation with controlled finger pressure: new equipment for the study of tender myofascial tissues. Pain 59, 235–239. doi: 10.1016/0304-3959(94)90076-0

Bialosky, J. E., Bishop, M. D., Price, D. D., Robinson, M. E., and George, S. Z. (2009). The mechanisms of manual therapy in the treatment of musculoskeletal pain: a comprehensive model. Man. Ther. 14, 531–538. doi: 10.1016/j.math.2008.09.001

Billman, G. E. (2011). Heart rate variability – a historical perspective. Front. Physiol. 2:86. doi: 10.3389/fphys.2011.00086

Billman, G. E. (2013). The LF/HF ratio does not accurately measure cardiac sympatho-vagal balance. Front. Physiol. 4:26. doi: 10.3389/fphys.2013.00026

Bombardi, C., Grandis, A., Chiocchetti, R., Bortolami, R., Johansson, H., and Lucchi, M. L. (2006). Immunohistochemical localization of alpha(1a)-adrenoreceptors in muscle spindles of rabbit masseter muscle. Tissue. Cell 38, 121–125. doi: 10.1016/j.tice.2005.12.003

Brighina, F., Piazza, A., Vitello, G., Aloisio, A., Palermo, A., Daniele, O., et al. (2004). rTMS of the prefrontal cortex in the treatment of chronic migraine: a pilot study. J. Neurol. Sci. 227, 67–71. doi: 10.1016/j.jns.2004.08.008

Bron, C., and Dommerholt, J. D. (2012). Etiology of myofascial trigger points. Curr. Pain Headache Rep. 16, 439–444. doi: 10.1007/s11916-012-0289-4

Clark, G. T. (2008). Classification, causation and treatment of masticatory myogenous pain and dysfunction. Oral. Maxillofac. Surg. Clin. North. Am. 20, 145–157. doi: 10.1016/j.coms.2007.12.003

Côté, P., Cassidy, J. D., Carroll, L. J., and Kristman, V. (2004). The annual incidence and course of neck pain in the general population: a population-based cohort study. Pain 112, 267–273. doi: 10.1016/j.pain.2004.09.004

Critchley, H. D., Elliott, R., Mathias, C. J., and Dolan, R. J. (2000). Neural activity relating to generation and representation of galvanic skin conductance responses: a functional magnetic resonance imaging study. J. Neurosci. 20, 3033–3040.

Delaney, J. P., Leong, K. S., Watkins, A., and Brodie, D. (2002). The short-term effects of myofascial trigger point massage therapy on cardiac autonomic tone in healthy subjects. J. Adv. Nurs. 37, 364–371. doi: 10.1046/j.1365-2648.2002.02103.x

Delany, J. (2014). “Massage, bodywork, and touch therapies,” in Fundamentals of Complementary and Alternative Medicine, 5th Edn., ed M. S. Micozzi (London: Elsevier Health Sciences), 247–274.

Fejer, R., Kyvik, K. O., and Hartvigsen, J. (2006). The prevalence of neck pain in the world population: a systematic critical review of the literature. Eur. Spine. J. 15, 834–848. doi: 10.1007/s00586-004-0864-4

Fernández-de-las-Peñas, C., Alonso-Blanco, C., and Miangolarra, J. C. (2007). Myofascial trigger points in subjects presenting with mechanical neck pain: a blinded, controlled study. Man. Ther. 12, 29–33. doi: 10.1016/j.math.2006.02.002

Ferrari, M., and Quaresima, V. (2012). A brief review on the history of human functional near-infrared spectroscopy (fNIRS) development and fields of application. Neuroimage 63, 921–935. doi: 10.1016/j.neuroimage.2012.03.049

Fritz, H. C., McAuley, J. H., Wittfeld, K., Hegenscheid, K., Schmidt, C. O., Langner, S., et al. (2016). Chronic back pain is associated with decreased prefrontal and anterior insular gray matter: results from a population-based cohort study. J. Pain 17, 111–118. doi: 10.1016/j.jpain.2015.10.003

Fukui, Y., Ajichi, Y., and Okada, E. (2003). Monte Carlo prediction of near-infrared light propagation in realistic adult and neonatal head models. Appl. Opt. 42, 2881–2887. doi: 10.1364/AO.42.002881

Furlan, R., Colombo, S., Perego, F., Atzeni, F., Diana, A., Barbic, F., et al. (2005). Abnormalities of cardiovascular neural control and reduced orthostatic tolerance in patients with primary fibromyalgia. J. Rheumatol. 32, 1787–1793.

Gasior, J. S., Sacha, J., Jeleń, P. J., Zieliński, J., and Przybylski, J. (2016). Heart rate and respiratory rate influence on heart rate variability repeatability: effects of the correction for the prevailing heart rate. Front. Physiol. 7:356. doi: 10.3389/fphys.2016.00356

Ge, H. Y., Fernández-de-las-Peñas, C., and Arendt-Nielsen, L. (2006). Sympathetic facilitation of hyperalgesia evoked from myofascial tender and trigger points in patients with unilateral shoulder pain. Clin. Neurophysiol. 117, 1545–1550. doi: 10.1016/j.clinph.2006.03.026

Gerwin, R. D., Dommerholt, J., and Shah, J. P. (2004). An expansion of Simons' integrated hypothesis of trigger point formation. Curr. Pain. Headache. Rep. 8, 468–475. doi: 10.1007/s11916-004-0069-x

Gerwin, R. D., Shannon, S., Hong, C. Z., Hubbard, D., and Gevirtz, R. (1997). Interrater reliability in myofascial trigger point examination. Pain 69, 65–73. doi: 10.1016/S0304-3959(96)03248-4

Gockel, M., Lindholm, H., Alaranta, H., Viljanen, A., Lindquist, A., and Lindholm, T. (1995). Cardiovascular functional disorder and stress among patients having neck–shoulder symptoms. Ann. Rheum. Dis. 54, 494–497. doi: 10.1136/ard.54.6.494

Guez, M., Hildingsson, C., and Toolanen, G. (2002). The prevalence of neck pain: a population-based study from northern Sweden. Acta. Orthop. Scand. 73, 455–459. doi: 10.1080/00016470216329

Hadjipavlou, G., Dunckley, P., Behrens, T. E., and Tracey, I. (2006). Determining anatomical connectivities between cortical and brainstem pain processing regions in humans: a diffusion tensor imaging study in healthy controls. Pain 123, 169–178. doi: 10.1016/j.pain.2006.02.027

Hallman, D., and Lyskov, E. (2012). “Autonomic regulation in musculoskeletal pain,” in Pain in Perspective, ed S. Ghosh (Rijeka: InTech), 35–62.

Hallman, D. M., Lindberg, L. G., Arnetz, B. B., and Lyskov, E. (2011). Effects of static contraction and cold stimulation on cardiovascular autonomic indices, trapezius blood flow and muscle activity in chronic neck–shoulder pain. Eur. J. Appl. Physiol. 111, 1725–1735. doi: 10.1007/s00421-010-1813-z

Heathers, J. A. (2014). Everything Hertz: methodological issues in short-term frequency-domain HRV. Front. Physiol. 5:177. doi: 10.3389/fphys.2014.00177

Hong, C. Z., and Simons, D. G. (1998). Pathophysiologic and electrophysiologic mechanisms of myofascial trigger points. Arch. Phys. Med. Rehabil. 79, 863–872. doi: 10.1016/S0003-9993(98)90371-9

Hoshi, Y., Kobayashi, N., and Tamura, M. (2001). Interpretation of near-infrared spectroscopy signals: a study with a newly developed perfused rat brain model. J. Appl. Physiol. (1985). 90, 1657–1662.

Hou, C. R., Tsai, L. C., Cheng, K. F., Chung, K. C., and Hong, C. Z. (2002). Immediate effects of various physical therapeutic modalities on cervical myofascial pain and trigger-point sensitivity. Arch. Phys. Med. Rehabili. 83, 1406–1414. doi: 10.1053/apmr.2002.34834

Hoy, D. G., Protani, M., De, R., and Buchbinder, R. (2010). The epidemiology of neck pain. Best. Pract. Res. Clin. Rheumatol. 24, 783–792. doi: 10.1016/j.berh.2011.01.019

Hubbard, D. R., and Berkoff, G. M. (1993). Myofascial trigger points show spontaneous needle EMG activity. Spine (Phila Pa 1976). 18, 1803–1807. doi: 10.1097/00007632-199310000-00015

Ishikuro, K., Urakawa, S., Takamoto, K., Ishikawa, A., Ono, T., and Nishijo, H. (2014). Cerebral functional imaging using near-infrared spectroscopy during repeated performances of motor rehabilitation tasks tested on healthy subjects. Front. Hum. Neurosci. 8:292. doi: 10.3389/fnhum.2014.00292

Kawakita, K., and Itoh, K. (2002). The polymodal receptor hypothesis of acupuncture and moxibustion, and its rational explanation of acupuncture points. International. Congress. Series. 1238, 63–68. doi: 10.1016/S0531-5131(02)00411-9

Kong, J., Gollub, R., Huang, T., Polich, G., Napadow, V., Hui, K., et al. (2007). Acupuncture de qi, from qualitative history to quantitative measurement. J. Altern. Complement. Med. 13, 1059–1070. doi: 10.1089/acm.2007.0524

Kostopoulos, D., Nelson, A. J., Ingber, R. S., and Larkin, R. W. (2008). Reduction of spontaneous electrical activity and pain perception of trigger points in the upper trapezius muscle through trigger point compression and passive stretching. J. Musculoskelet. Pain 16, 267–279. doi: 10.1080/10582450802479594

Kuo, T. C., Chen, Z. S., Chen, C. H., Ho, F. M., Lin, C. W., and Chen, Y. J. (2004). The physiological effect of de qi during acupuncture. J. Health. Sci. 50, 336–342. doi: 10.1248/jhs.50.336

Larsson, S. E., Bodegard, L., Henriksson, K. G., and Oberg, P. A. (1990). Chronic trapezius myalgia. Morphology and blood flow studied in 17 patients. Acta Orthop. Scand. 61, 394–398. doi: 10.3109/17453679008993548

Li, X., Hirokawa, M., Inoue, Y., Sugano, N., Qian, S., and Iwai, T. (2007). Effects of acupressure on lower limb blood flow for the treatment of peripheral arterial occlusive diseases. Surg. Today 37, 103–108. doi: 10.1007/s00595-006-3347-x

Lombardi, F., and Stein, P. K. (2011). Origin of heart rate variability and turbulence: an appraisal of autonomic modulation of cardiovascular function. Front. Physiol. 2:95. doi: 10.3389/fphys.2011.00095

Lorenz, J., Minoshima, S., and Casey, K. L. (2003). Keeping pain out of mind: the role of the dorsolateral prefrontal cortex in pain modulation. Brain 126, 1079–1091. doi: 10.1093/brain/awg102

Lu, G. W. (1983). Characteristics of afferent fiber innervation on acupuncture points zusanli. Am. J. Physiol. 245, R606–R612.

Martinez-Lavin, M. (2007). Biology and therapy of fibromyalgia. Stress, the stress response system, and fibromyalgia. Arthritis. Res. Ther. 9:216. doi: 10.1186/ar2146

McNulty, W. H., Gevirtz, R. N., Hubbard, D. R., and Berkoff, G. M. (1994). Needle electromyographic evaluation of trigger point response to a psychological stressor. Psychophysiology 31, 313–316. doi: 10.1111/j.1469-8986.1994.tb02220.x

Melzack, R., Stillwell, D. M., and Fox, E. J. (1997). Trigger points and acupuncture points for pain: correlations and implications. Pain 3, 3–23. doi: 10.1016/0304-3959(77)90032-X

Muñoz-Muñoz, S., Muñoz-García, M. T., Alburquerque-Sendín, F., Arroyo-Morales, M., and Fernández-de-las-Peñas, C. (2012). Myofascial trigger points, pain, disability, and sleep quality in individuals with mechanical neck pain. J. Manipulative. Physiol. Ther. 35, 608–613. doi: 10.1016/j.jmpt.2012.09.003

Muthukumaraswamy, S. D., Edden, R. A., Jones, D. K., Swettenham, J. B., and Singh, K. D. (2009). Resting GABA concentration predicts peak gamma frequency and fMRI amplitude in response to visual stimulation in humans. Proc. Natl. Acad. Sci. U.S.A. 106, 8356–8361. doi: 10.1073/pnas.0900728106

Napadow, V., Lee, J., Kim, J., Cina, S., Maeda, Y., Barbieri, R., et al. (2013). Brain correlates of phasic autonomic response to acupuncture stimulation: an event-related fMRI study. Hum. Brain. Mapp. 34, 2592–2606. doi: 10.1002/hbm.22091

Niddam, D. M., Lee, S.-H., Su, Y.-T., and Chan, R.-C. (2017). Brain structural changes in patients with chronic myofascial pain. Eur. J. Pain. 21, 148–158. doi: 10.1002/ejp.911

Northoff, G., Walter, M., Schulte, R. F., Beck, J., Dydak, U., Henning, A., et al. (2007). GABA concentrations in the human anterior cingulate cortex predict negative BOLD responses in fMRI. Nat. Neurosci. 10, 1515–1517. doi: 10.1038/nn2001

Okamoto, M., Matsunami, M., Dan, H., Kohata, T., Kohyama, K., and Dan, I. (2006). Prefrontal activity during taste encoding: an fNIRS study. Neuroimage 31, 796–806. doi: 10.1016/j.neuroimage.2005.12.021

Ongür, D., An, X., and Price, J. L. (1998). Prefrontal cortical projections to the hypothalamus in macaque monkeys. J. Comp. Neurol. 401, 480–505.

Pagani, M., Montano, N., Porta, A., Malliani, A., Abboud, F. M., Birkett, C., et al. (1997). Relationship between spectral components of cardiovascular variabilities and direct measures of muscle sympathetic nerve activity in humans. Circulation 95, 1441–1448. doi: 10.1161/01.CIR.95.6.1441

Partanen, J. V., Ojala, T. A., and Arokoski, J. P. A. (2010). Myofascial syndrome and pain: a neurophysiological approach. Pathophysiology 17, 19–28. doi: 10.1016/j.pathophys.2009.05.001

Passatore, M., and Roatta, S. (2006). Influence of sympathetic nervous system on sensorimotor function: whiplash associated disorders (WAD) as a model. Eur. J. Appl. Physiol. 98, 423–449. doi: 10.1007/s00421-006-0312-8

Roth, L. U., Maret-Maric, A., Adler, R. H., and Neuenschwander, B. E. (1997). Acupuncture points have subjective (needling sensation) and objective (serum cortisol increase) specificity. Acupunct. Med. 15, 2–5. doi: 10.1136/aim.15.1.2

Sakai, S., Hori, E., Umeno, K., Kitabayashi, N., Ono, T., and Nishijo, H. (2007). Specific acupuncture sensation correlates with EEGs and autonomic changes in human subjects. Autonom. Neurosci. 133, 158–169. doi: 10.1016/j.autneu.2007.01.001

Sato, H., Fuchino, Y., Kiguchi, M., Katura, T., Maki, A., Yoro, T., et al. (2005). Intersubject variability of near-infrared spectroscopy signals during sensorimotor cortex activation. J. Biomed. Opt. 10, 44001. doi: 10.1117/1.1960907

Schroeter, M. L., Zysset, S., Kruggel, F., and von Cramon, D. Y. (2003). Age dependency of the hemodynamic response as measured by functional near-infrared spectroscopy. Neuroimage. 19, 555–564. doi: 10.1016/S1053-8119(03)00155-1

Seifert, F., Schuberth, N., De Col, R., Peltz, E., Nickel, F. T., and Maihöfner, C. (2013). Brain activity during sympathetic response in anticipation and experience of pain. Hum. Brain Mapp. 34, 1768–1782. doi: 10.1002/hbm.22035

Seitz, R. J., and Roland, P. E. (1992). Vibratory stimulation increases and decreases the regional cerebral blood flow and oxidative metabolism: a positron emission tomography (PET) study. Acta. Neurol. Scand. 86, 60–67. doi: 10.1111/j.1600-0404.1992.tb08055.x

Seiyama, A., Hazeki, O., and Tamuram, M. (1988). Noninvasive quantitative analysis of blood oxygenation in rat skeletal muscle. J. Biochem. 103, 419–424. doi: 10.1093/oxfordjournals.jbchem.a122285

Selkowitz, D. M. (1992). The sympathetic nervous system in neuromotor function and dysfunction and pain. A brief review and discussion. Funct. Neurol. 7, 89–95.

Shmuel, A., Augath, M., Oeltermann, A., and Logothetis, N. K. (2006). Negative functional MRI response correlates with decreases in neuronal activity in monkey visual area V1. Nat. Neurosci. 9, 569–577. doi: 10.1038/nn1675

Shmuel, A., Yacoub, E., Pfeuffer, J., Van de Moortele, P. F., Adriany, G., Hu, X., et al. (2002). Sustained negative BOLD, blood flow and oxygen consumption response and its coupling to the positive response in the human brain. Neuron 36, 1195–1210. doi: 10.1016/S0896-6273(02)01061-9

Simons, D. G. (1995). The nature of myofascial trigger points. Clin. J. Pain. 11, 83–84.

Simons, D. G., Travell, J. G., and Simons, L. S. (1999). Myofascial Pain and Dysfunction, The Trigger Point Manual, The Upper Half of Body, Vol. 1, 2nd Edn. Baltimore, MD: Williams and Wilkins Press.

Srbely, J. Z. (2010). New trends in the treatment and management of myofascial pain syndrome. Curr. Pain. Headache. Rep. 14, 346e352. doi: 10.1007/s11916-010-0128-4

Stefanovic, B., Warnking, J. M., and Pike, G. B. (2004). Hemodynamic and metabolic responses to neuronal inhibitation. Neuroimage 22, 771–778. doi: 10.1016/j.neuroimage.2004.01.036

Strangman, G., Culver, J. P., Thompson, J. H., and Boas, D. A. (2002). A quantitative comparison of simultaneous BOLD fMRI and NIRS recordings during functional brain activation. Neuroimage 17, 719–731. doi: 10.1006/nimg.2002.1227

Suzuki, M., Miyai, I., Ono, T., and Kubota, K. (2008). Activities in the frontal cortex and gait performance are modulated by preparation. An fNIRS study. Neuroimage 39, 600–607. doi: 10.1016/j.neuroimage.2007.08.044

Takamoto, K., Hori, E., Urakawa, S., Sakai, S., Ishikawa, A., Kohno, S., et al. (2010). Cerebral hemodynamic responses induced by specific acupuncture sensations during needling at trigger points: a near-infrared spectroscopic study. Brain. Topogr. 23, 279–291. doi: 10.1007/s10548-010-0148-8

Takamoto, K., Sakai, S., Hori, E., Urakawa, S., Umeno, K., Ono, T., et al. (2009). Compression on trigger points in the leg muscle increases parasympathetic nervous activity based on heart rate variability. J. Physiol. Sci. 59, 191–197. doi: 10.1007/s12576-009-0025-y

Takeuchi, M., Hori, E., Takamoto, K., Tran, A. H., Satoru, K., Ishikawa, A., et al. (2009). Brain cortical mapping by simultaneous recording of functional near infrared spectroscopy and electroencephalograms from the whole brain during right median nerve stimulation. Brain. Topogr. 22, 197–214. doi: 10.1007/s10548-009-0109-2

Tarvainen, M. P., Niskanen, J. P., Lipponen, J. A., Ranta-Aho, P. O., and Karjalainen, P. A. (2014). Kubios HRV – heart rate variability analysis software. Comput. Methods. Prog. Biomed. 113, 210–220. doi: 10.1016/j.cmpb.2013.07.024

Task Force of the European Society of Cardiology the North American Society of Pacing Electrophysiology (1996). Heart rate variability: standards of measurement, physiological interpretation and clinical use. Eur. Heart. J. 17, 354–381. doi: 10.1093/oxfordjournals.eurheartj.a014868

Thayer, J. F., Yamamoto, S. S., and Brosschot, J. F. (2010). The relationship of autonomic imbalance, heart rate variability and cardiovascular disease risk factors. Int. J. Cardiol. 141, 122–131. doi: 10.1016/j.ijcard.2009.09.543

Toichi, M., Findling, R. L., Kubota, Y., Calabrese, J. R., Wiznitzer, M., McNamara, N. K., et al. (2004). Hemodynamic differences in the activation of the prefrontal cortex: attention vs. higher cognitive processing. Neuropsychologia 42, 698–706. doi: 10.1016/j.neuropsychologia.2003.08.012

Tsuzuki, D., Jurcak, V., Singh, A. K., Okamoto, M., Watanabe, E., and Dan, I. (2007). Virtual spatial registration of stand-alone fNIRS data to MNI space. Neuroimage 34, 1506–1518. doi: 10.1016/j.neuroimage.2006.10.043

Umeno, K., Hori, E., Tabuchi, E., Takakura, H., Miyamoto, K., Ono, T., et al. (2003). Gamma-band EEGs predict autonomic responses during mental arithmetic. Neuroreport 14, 477–480. doi: 10.1097/00001756-200303030-00036

Wang, K. M., Yao, S. M., Xian, Y. L., and Hou, Z. L. (1985). A study on the receptive field of acupoints and the relationship between characteristics of needling sensation and groups of afferent fibres. Sci. Sin. B 28, 963–971.

Wray, S., Cope, M., Delpy, D. T., Wyatt, J. S., and Reynolds, E. O. (1988). Characterization of the near infrared absorption spectra of cytochrome aa3 and haemoglobin for the non-invasive monitoring of cerebral oxygenation. Biochim. Biophys. Acta 933, 184–192. doi: 10.1016/0005-2728(88)90069-2

Yamamoto, T., and Kato, T. (2002). Paradoxical correlation between signal in functional magnetic resonance imaging and deoxygenated hemoglobin content in capillaries: a new theoretical explanation. Phys. Med. Biol. 47, 1121–1141. doi: 10.1088/0031-9155/47/7/309

Yasui, H., Takamoto, K., Hori, E., Urakawa, S., Nagashima, Y., Yada, Y., et al. (2010). Significant correlation between autonomic nervous activity and cerebral hemodynamics during thermotherapy on the neck. Autonom. Neurosci. 156, 96–103. doi: 10.1016/j.autneu.2010.03.011

Ye, J. C., Tak, S., Jang, K. E., Jung, J., and Jang, J. (2009). NIRS-SPM: statistical parametric mapping for near-infrared spectroscopy. Neurimage 44, 428–447. doi: 10.1016/j.neuroimage.2008.08.036

Yip, Y. B., and Tse, S. H. (2004). The effectiveness of relaxation acupoint stimulation and acupressure with aromatic lavender essential oil for non-specific low back pain in Hong Kong: a randomised controlled trial. Complement. Ther. Med. 12, 28–37. doi: 10.1016/j.ctim.2003.12.003

Yip, Y. B., and Tse, S. H. (2006). An experimental study on the effectiveness of acupressure with aromatic lavender essential oil for sub-acute, non-specific neck pain in Hong Kong. Complement. Ther. Clin. Pract. 12, 18–26. doi: 10.1016/j.ctcp.2005.09.005