Glucose Metabolic Changes in the Brain and Muscles of Patients with Nonspecific Neck Pain Treated by Spinal Manipulation Therapy: A [18F]FDG PET Study

Abstract

Objective. The aim of this study was to investigate changes in brain and muscle glucose metabolism that are not yet known, using positron emission tomography with [18F]fluorodeoxyglucose ([18F]FDG PET). Methods. Twenty-one male volunteers were recruited for the present study. [18F]FDG PET scanning was performed twice on each subject: once after the spinal manipulation therapy (SMT) intervention (treatment condition) and once after resting (control condition). We performed the SMT intervention using an adjustment device. Glucose metabolism of the brain and skeletal muscles was measured and compared between the two conditions. In addition, we measured salivary amylase level as an index of autonomic nervous system (ANS) activity, as well as muscle tension and subjective pain intensity in each subject. Results. Changes in brain activity after SMT included activation of the dorsal anterior cingulate cortex, cerebellar vermis, and somatosensory association cortex and deactivation of the prefrontal cortex and temporal sites. Glucose uptake in skeletal muscles showed a trend toward decreased metabolism after SMT, although the difference was not significant. Other measurements indicated relaxation of cervical muscle tension, decrease in salivary amylase level (suppression of sympathetic nerve activity), and pain relief after SMT. Conclusion. Brain processing after SMT may lead to physiological relaxation via a decrease in sympathetic nerve activity.

Akie Inami 1 and Takeshi Ogura 1,2 and Shoichi Watanuki 1 and Md. Mehedi Masud 1,3 and Katsuhiko Shibuya 4 and Masayasu Miyake 1 and Rin Matsuda 1 and Kotaro Hiraoka 1 and Masatoshi Itoh 4 and Arlan W. Fuhr 5 and Kazuhiko Yanai 1,6 and Manabu Tashiro 1

Academic Editor:Jeng-Ren Duann

1, Division of Cyclotron Nuclear Medicine, Cyclotron and Radioisotope Center, Tohoku University, Sendai, Japan

2, Japan Chiropractic Doctor College, Sendai, Japan

3, Department of Nuclear Medicine, United Hospital, Dhaka, Bangladesh

4, Sendai Medical Imaging Clinic, Sendai, Japan

5, Activator Methods International, Ltd., Phoenix, AZ, USA

6, Department of Pharmacology, Graduate School of Medicine, Tohoku University, Sendai, Japan

Received 2 August 2016; Revised 24 November 2016; Accepted 15 December 2016; 12 January 2017

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1. Introduction

Spinal manipulation therapy (SMT), which is performed by healthcare practitioners such as chiropractors, osteopathic physicians, and physiotherapists, has been applied mainly to musculoskeletal problems such as neck pain or low back pain. Many investigators have performed various experiments to elucidate the mechanism underlying the clinical effects of SMT. The earliest studies analyzed the magnitude of the force applied to the vertebrae and the movement of vertebrae during SMT, with the forces exerted by the practitioner quantified using a flexible pressure mat [1, 2]. The rotation and relative movement of vertebrae after treatment using a mechanical adjusting device (activator adjusting instrument, AAI) [3] have also been reported [4-8]. Previous studies have suggested that SMT has beneficial clinical effects, including pain relief [9] and reduction of blood pressure [10, 11]. It is thought that biomechanical input by SMT could generate a physical response or reflex [12]; however, the mechanism of these clinical effects induced by SMT is still unknown.

In recent years, the findings of brain activation studies using imaging modalities such as functional magnetic resonance imaging, near-infrared spectroscopy, and positron emission tomography (PET) have contributed to advances in brain science [13]. Such studies are able to enhance our understanding of the neurophysiological effects of specific physical/psychological tasks by detecting associated activation and deactivation of brain regions. PET has also been used to elucidate the metabolic changes that alternative therapies induce in living tissues; for example, an ¹⁸ F-labeled glucose analog has been used to study cerebral metabolic changes after acupuncture [14] and aromatherapy [15], and [¹⁵ O]H₂ O has been used to assess changes in cerebral blood flow during massage [16]. Lystad and Pollard have noted the usefulness of neuroimaging techniques for gaining a better understanding of the neurophysiological effects of SMT [17].

The [¹⁸ F]-labeled glucose analog, fluorodeoxyglucose ([¹⁸ F]FDG), is thought to be able to visualize the energy metabolism of various tissues such as brain and muscles in vivo. One of important advantages of [¹⁸ F]FDG PET is that we can simultaneously measure the energy metabolism of brain and muscles. We have hypothesized that clinical effects such as muscle tension relaxation and pain relief are mediated by the regional brain activity induced by SMT. We have also expected any changes in muscular energy consumption in treatment condition. In addition, we have also expected significant correlation between the regional brain activity and the muscular energy consumption.

Previously, we reported preliminary results of our study to examine our initial hypothesis that a [¹⁸ F]FDG PET can visualize brain metabolic changes induced by SMT [18]. In this paper, we report our conclusive results on this issue. In addition, the present work tries to examine our additional hypothesis that the [¹⁸ F]FDG PET can visualize muscular metabolic changes induced by SMT, as well as their association with the regional brain metabolism.

2. Subjects and Methods

We recruited 21 male subjects (mean age ± SD: 26.4±5.9 years) with cervical pain and shoulder stiffness but without abnormalities in neck-to-shoulder MR images and without history of any treatments prior to the present study. We performed MRI examination of the neck to the shoulder area on all participants; the resulting MR images were used as a reference for the anatomical locations of cervical muscles in PET images. Female subjects were not included in this study because of physiological fluctuations due to the menstrual cycle [19-23].

SMT was applied using an AAI, Activator II (Activator Methods International, Ltd., Phoenix, USA), in accordance with the Activator Methods (AM) basic scan protocol [24, 25]. We utilized the AAI to apply impulses to specific vertebrae or joints (Figure 1). SMT was performed on the subject in a prone position without movements such as cervical rotation, lateral flexion, and extension, in order to prevent the muscular [¹⁸ F]FDG uptake due to muscle contractions during the therapeutic procedure. SMT was carried out on the whole spine, the scapulae, the ilium, and the sacrum, as necessary for each subject. The mean number of SMT-adjusted sites was 8 per subject.

Figure 1: The location of cervical adjustment by the activator adjusting instrument (AAI) in the treatment condition. Spinal manipulation therapy (SMT) by AAI was performed by contact on the joints and did not include muscle massage. SMT on all subjects was performed by the same chiropractor, who was an advanced practitioner of Activator Methods.

[figure omitted; refer to PDF]

PET scanning was performed twice on each subject according to the protocol given in Figure 2. The interval between the two scans ranged from 7 to 70 days (mean interval ± SD: 23±15 days). The subjects received SMT intervention which lasted for approximately 20 minutes including a diagnostic procedure. Soon after the treatment, [¹⁸ F]FDG-containing saline solution was injected to the subject through the left cubital vein (mean ± SD: 47.0±8.9 MBq) in quiet room with a dim light next door to the treatment room in the same building. On the other scan day, [¹⁸ F]FDG was injected to the subject after a 20-minute-long resting phase instead of SMT intervention. The subjects in both conditions were asked to sit in a relaxed position with their eyes closed for 30 minutes before PET scanning. The brain scan and the scan on the neck to the shoulder area of the subject were initiated after 30 minutes of [¹⁸ F]FDG injection, utilizing a PET scanner, SET2400W (Shimadzu, Inc., Kyoto, Japan). The PET scanning covered the entire brain in one scan, taking 10 minutes for the emission scan and another 5 minutes for the transmission scan for tissue attenuation correction. Scanning from the neck to the shoulder area in one scan took 5 minutes for the emission scan and another 5 minutes for the transmission scan for tissue attenuation correction. Images were acquired with a 128 × 128 matrix and reconstructed using Fourier rebinning and Ordered Subset Expectation Maximization Algorithm [18]. The intensity of subjective pain was evaluated using a 0-10 visual analog scale (VAS) before and after SMT in the treatment condition. VAS evaluation was not done in the resting condition except for in 9 subjects. Cervical muscle tension was measured bilaterally at the superior part of the trapezius muscle using a tissue hardness meter (Muscle Meter PEK-1, Imoto Inc., Kyoto, Japan); the mean value of three measurements was recorded. Additionally, salivary amylase levels were measured for each subject using an amylase monitor (Nipro Inc., Osaka, Japan) to evaluate changes in autonomic nervous system (ANS) function. Figure 2 shows the measurement points. Further details of the study protocol are described in our previous report [18]. The whole protocol was approved by the Ethics Committee of Tohoku University Graduate School of Medicine, Sendai, Japan (number 2008-115).

Figure 2: Diagram of the study protocol. Half of the subjects were randomly assigned to be scanned first in the resting condition and then in the treatment condition; the other half were scanned in the reverse order. The intensity of subjective pain was evaluated before and after spinal manipulation therapy only in the treatment condition. Muscle stiffness and salivary amylase were measured before and after the treatment or resting period. MT = muscle tension, VAS = visual analog scale.

[figure omitted; refer to PDF]

For data analysis, differences in values before and after treatment and between the control and treatment conditions were compared using paired t-tests for all measurements except for PET data. Brain PET images were analyzed using the voxel-wise statistical analysis software package Statistical Parametric Mapping 8 (SPM8; Functional Imaging Laboratory, London, UK) in order to identify regional glucose metabolic changes [26, 27]. An FDG brain template distributed by Montreal Neurological Institute (McGill University, Montreal, Canada) [28] was used for anatomical standardization (spatial normalization) of the PET images. Each voxel had dimensions of 2 × 2 × 2 mm in the normalized image. The normalized data were smoothed using an isotropic Gaussian kernel of 12 mm (for the x-, y-, and z-axes) to increase the signal-to-noise ratio by suppressing high-frequency noise. The threshold for the statistical significance of the voxel value height in the present study was set at p<0.05 with correction for multiple comparisons (family-wise error correction), while our previous study applied the compromised threshold of p<0.001 without correction. Voxel values of the PET images were compared between the resting and treatment conditions using a paired t-test.

PET images of the neck and shoulder regions were coregistered to the MR images of the same subject; regions of interest (ROIs) for cervical muscles were then manually drawn on the PET images using Dr. View software (Version 2.0, AJS, Tokyo, Japan), using the MR images as references. ROIs were drawn on the trapezius muscle at C7-T1 levels bilaterally, the splenius muscles, the semispinalis muscles, the elevator scapular muscles, and the trapezius muscles at C6-C7 levels bilaterally. The standardized uptake value (SUV) for each muscle was calculated using the following formula: [figure omitted; refer to PDF] SUVs were statistically examined using a paired t-test to compare muscle metabolism after SMT and after resting. In addition, the authors searched for specific brain regions associated with muscular energy consumption by applying linear correlation analysis using SUV of each muscle studied here as a parameter in SPM8.

3. Results

The SMT-associated regional brain metabolic changes (activation and deactivation) detected by the SPM8 analysis are shown in Table 1. Statistically significant areas were overlaid on the standard MRI brain template images (Figure 3). PET analysis of the cervical muscles showed a trend toward reduced metabolism (SUV) after SMT compared with the control condition; however, these changes were not statistically significant (Figure 4). In addition, no meaningful correlation was detected between muscular SUV and regional brain activity.

Table 1: Brain metabolic changes associated with the spinal manipulation therapy intervention.

Anatomical region	Coordinates x, y, and z (mm)	Brodmann area	Cluster equiv.	Voxel Z score
Activation
Broca's area	-34, 6, 28	44	110	5.84
ACC	2, 8, 40	32	228	5.72
SSAC	16, -26, 48	5	114	5.65
Wernicke's area	46, -48, 20	22	30	5.52
VAC	-6, -88, 38	19	31	5.32
CV	6, -62, -4	--	14	5.28
VC (V2)	24, -80, 6	18	17	5.24
Deactivation
IPL	-40, -40, 34	39/40	286	6.74
FP	-2, 68, -8	10	160	6.66
IFG PT	40, 28, 16	45	117	6.53
PSMA	30, 14, 44	6	348	6.34
PMC (M1)	-28, -18, 56	4	45	6.12
FEF/dl-PFC	-30, 12, 44	8/9	157	5.93
dl-PFC	-38, 26, 24	46	309	5.70
AG/FG	-40, -60, 8	39/37	177	5.58
ITG	-70, -22, -20	20	20	5.66
TP	40, 22, -44	38	12	5.64
CV (V1)	-14, -66, 24	17	46	5.29

Brain metabolic changes detected by SPM8 are presented (voxel height threshold p<0.05 with corrections for multiple comparisons, extent threshold 10 voxels minimum). The statistical significance of regional metabolic changes is given as Z scores [(Mea_ntreatment -Mea_ncontrol )/S_Dcontrol ].

ACC, anterior cingulate cortex; SSAC, somatosensory association cortex; VAC, visual association cortex; CV, cerebellar vermis; VC, visual cortex; IPL, inferior parietal lobule; FP, frontal pole; IFG, inferior frontal gyrus; PT, pars triangularis; PSMA, premotor area/supplementary motor area; PMC, primary motor cortex; FEF, frontal eye field; dl-PFC, dorsolateral prefrontal cortex; AG, angular gyrus; FC, fusiform gyrus; ITG, inferior temporal gyrus; TP, temporal pole.

Figure 3: Regional activation (left) and deactivation (right) after spinal manipulation therapy (SMT) using an activator adjusting instrument. Glucose metabolism is increased in regions including the anterior cingulate cortex and cerebellar vermis but is relatively reduced in many sites, including the prefrontal cortex, after SMT. The voxel height threshold is p<0.05, corrected for multiple comparisons; the extent threshold is 10 voxels minimum.

[figure omitted; refer to PDF]

Figure 4: Results of positron emission tomography analysis of cervical muscles (paired t-test). The error bars represent standard deviations. The results indicate a trend toward reduction of mean standardized uptake value (SUV) after SMT; however, the difference is not statistically significant. Trap., trapezius muscle; Splenius, splenius muscles; Lev. Scap., levator scapulae; Semispinalis, semispinalis muscles; C7-T1, between the seventh cervical spine and the first thoracic spine; C6-C7, between the sixth cervical spine and the seventh cervical spine.

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

In contrast, cervical muscle tension was significantly reduced bilaterally after SMT (p<0.0001 for both sides, Figure 5). Salivary amylase level decreased significantly after SMT (p=0.022) but increased significantly in the control condition (p=0.011, Figure 6). Comparisons of VAS pain scores in the treatment condition revealed a significant decrease after SMT (p<0.0001, Figure 7), while the comparison in the control condition showed nonsignificant difference (n=9).

Figure 5: Muscle tension is significantly reduced after spinal manipulation therapy. p[low *]<0.0001.

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

Figure 6: Changes in salivary amylase level. Salivary amylase level is reduced after spinal manipulation therapy but increased in the control condition. p[low *][low *]<0.05.

[figure omitted; refer to PDF]

Figure 7: Changes in subjective pain in the treatment condition. The pain scale score is significantly decreased after spinal manipulation therapy. p[low *]<0.0001.

[figure omitted; refer to PDF]

4. Discussion

The findings of the present study demonstrate how stimuli to the mechanoreceptors of the joints and skin during SMT are processed in the brain. Injected [¹⁸ F]FDG was absorbed into activated brain regions and visualized by PET. We observed multiple changes in brain activity after SMT.

SPM8 analyzes approximately half a million voxels of brain volume data simultaneously. Correction for multiple comparisons is therefore indispensable, making the statistical threshold extremely high. Since many studies have failed to detect significant differences in voxels after correction for multiple comparisons, the SPM8 development team proposed the use of a compromised threshold for voxel height (p<0.001) combined with a voxel extent threshold for the size of each voxel cluster (e.g., 10 voxels minimum), as used in our previous report [26, 27]. This technique has been useful for practical purposes but is prone to Type-1 errors. The significant voxel clusters that survived correction for multiple comparisons in the present work may therefore indicate more robust and reliable findings than those in our preliminary report [18]; we believe these results are worthy of reporting as conclusive findings, despite the fact that the intensity of the observed brain responses to SMT intervention was initially estimated to be much weaker.

As for regional brain metabolic changes after the SMT intervention, activation (increased metabolism) was detected in the dACC (Brodmann area [BA] 32), cerebellar vermis (CV), and somatosensory association cortex, and regional deactivation (decreased metabolism) was detected in regions including the prefrontal cortex (PFC) and temporal sites. Involvement of the ACC in cognitive and emotional processes has been recognized since Papez mentioned the idea in 1937 [29-31], and this area is also involved in placebo and opioid analgesia [32]. Some brain activation studies have also demonstrated activation of the ACC in response to pleasant or unpleasant stimuli such as massages or the olfactory stimulus of isovaleric acid, respectively [33-35]. The ACC is part of network that carries out cognitive processing based on individual factors such as experiences or emotions while making contact with the other regions in the network, such as the limbic system and cortex [32, 36]. Recently, specific features of the dACC functions and their connectivity to other brain regions are not fully elucidated yet, though they have been revealed little by little. For example, dACC is involved in cognitive functions, motivation, and reward-based decision-making as a part of the network including the CV and prefrontal cortex [37]. The neuronal activity changes in dACC and CV were detected in the present study.

The CV receives somatic sensory information from the spinal cord and via the vestibulospinal tract or reticular nuclei of the brainstem through the spinal cord, connecting indirectly or directly with motor cells on the ventral horn. These systems control involuntary muscular tension and reflexes. Our results suggest that stimulation of joints during SMT induced relaxation of reflexive muscle tension. The cerebellum is also thought to have a functional role as an integrator of multiple effector systems, including affective processing, pain modulation, and sensorimotor processing [38]. Recently, many studies have reported roles for the cerebellum in nonmotor functions [39, 40]. A study by Sacchetti et al. showed that the CV is activated during mental recall of emotional personal episodes in humans [39]. Lou et al. also reported that the CV, ACC, and some regions of the PFC are activated during relaxation mediation in yoga [41]. However, it is important to note that the PFC was deactivated after SMT in our study. Interestingly, Critchley et al. found that the ACC/dACC (BA24/32) and CV are specifically activated during biofeedback therapy [42], a technique for controlling one's tension to generate a state of relaxation. The activated areas in the present study are similar to those activated during biofeedback relaxation, indicating that the state of the brain after SMT may be similar to that induced by biofeedback therapy. Furthermore, our assessment of body responses in this study showed relaxation of muscle tension and decreased salivary amylase levels--phenomena that are associated with reduced sympathetic nerve activity. Salivary α -amylase levels correspond to plasma norepinephrine levels and are utilized as an accessible measure of sympathetic nervous reactivity in stress research, with lower levels indicating lower activity [43-45].

Ouchi et al. suggested that the comfortable sensation generated by back massage may be related to increased regional cerebral blood flow in the posterior brain--specifically, in the precuneus [16]. SMT stimuli to the joints may be processed differently from those of muscle massage, resulting in decreased sympathetic nerve activity. On the other hand, certain cervical muscles showed a tendency toward decreased glucose metabolism after SMT, although the difference was not statistically significant (Figure 4). Increasing the number of study subjects may increase the significance of this finding. The underlying mechanism of reduced muscular glucose uptake is not yet understood; however, an animal study suggests the involvement of sympathetic nerve activity [46]. Glucose uptake in skeletal muscles may thus be influenced directly or indirectly by the ANS.

Although the mechanism of muscle relaxation is still unknown, we hypothesize the involvement of (a) autonomic nervous activity and (b) improvement of the range of joint movement [47]. In the present study, we observed that SMT stimulus induced physical responses such as muscle tension relaxation, pain relief, and reduced amylase secretion. These changes may be associated with neural processing in the dACC and CV. Neural inputs evoked by SMT stimuli via various receptors in muscles, tendons, and joints may ascend to the somatosensory areas of the brain through the medial lemniscal system. The response signal descending from the brain may then adjust the postural muscles accordingly.

Muscular SUV, an index of energy consumption of skeletal muscles, did not show significant difference between the control and treatment conditions, while muscle tension did show significant difference. In addition, we could not detect meaningful correlation between the SUV of muscular glucose consumption and the regional brain glucose consumption.

As is the case with every study, our methodology is not without limitations. Ideally, the present subjects would have undergone PET scanning before and after the SMT intervention in both the control and treatment conditions. However, this protocol would have assigned 4 PET scans to each subject, resulting in unreasonably high radiation exposure for this kind of study. Furthermore, it is hard to perform serial PET scans within 1 or 2 hours of each other because of the relatively long physical half-life of the [¹⁸ F] nuclide (110 minutes). In addition, the number of subjects (n=21) was still relatively small for a clinical study; however, we define the present results as conclusive ones in order to minimize radiation exposure of the subjects, who were healthy other than their neck and shoulder symptoms and who are considered to be part of the general public. Previously, we reported preliminary findings with a sample size of 12; however, the study results were only at threshold level [18], and these results seemed to be prone to Type-1 errors because they were based on a compromised statistical examination without correction for multiple comparisons. On the other hand, no voxels survived with a standard statistical examination including correction for multiple comparisons, suggesting that these negative results may be prone to Type-2 errors. By raising the sample size to 21, now we are able to obtain robust results that survived even after correction for multiple comparisons. Thus, the authors believe the present results are more reliable as an evidence for further discussion on clinical effects of SMT interventions while our preliminary report was useful as a "proof of concept" study.

5. Conclusion

In summary, we observed metabolic changes in the brain and skeletal muscles, as well as reductions in subjective pain, muscle tension, and salivary amylase, after SMT intervention. These results may be associated with reduced sympathetic nerve activity, suggesting that SMT induces a kind of relaxation similar to that achieved by biofeedback. The brain response to SMT may reflect the psychophysiological relaxation that accompanies reduced sympathetic nerve activity.

Acknowledgments

The authors thank the staff of the Division of Cyclotron and Radioisotope Center, Tohoku University, for their technical support. They thank Japan Chiropractic Doctor College for funding this study. This work was also supported in part by a Grant-in-Aid for Scientific Research (no. 21590754 to Manabu Tashiro) from the Japan Society of Promotion of Sciences (JSPS) and by a grant from the Japan Society of Technology (JST) for research and education in molecular imaging.

References

[1] B. W. Hessell, W. Herzog, P. J. W. Conway, M. C. McEwen, "Experimental measurement of the force exerted during spinal manipulation using the Thompson technique,", Journal of Manipulative and Physiological Therapeutics , vol. 13, no. 8, pp. 448-453, 1990.

[2] W. Herzog, P. J. Conway, G. N. Kawchuk, E. M. Hasler, Y. Zhang, "Forces exerted during spinal manipulative therapy,", Spine , vol. 18, no. 9, pp. 1206-1212, 1993.

[3] A. W. Fuhr, J. M. Menke, "Status of activator methods chiropractic technique, theory, and practice,", Journal of Manipulative and Physiological Therapeutics , vol. 28, no. 2, pp. e1-e20, 2005.

[4] M. Nathan, T. S. Keller, "Measurement and analysis of the in vivo posteroanterior impulse response of the human thoracolumbar spine: a feasibility study,", Journal of Manipulative and Physiological Therapeutics , vol. 17, no. 7, pp. 431-441, 1994.

[5] D. B. Smith, A. W. Fuhr, B. P. Davis, "Skin accelerometer displacement and relative bone movement of adjacent vertebrae in response to chiropractic percussion thrusts,", Journal of Manipulative and Physiological Therapeutics , vol. 12, no. 1, pp. 26-37, 1989.

[6] A. W. Fuhr, D. B. Smith, "Accuracy of piezoelectric accelerometers measuring displacement of a spinal adjusting instrument,", Journal of Manipulative and Physiological Therapeutics , vol. 9, no. 1, pp. 15-21, 1986.

[7] T. S. Keller, C. J. Colloca, A. W. Fuhr, "Validation of the force and frequency characteristics of the Activator Adjusting Instrument: effectiveness as a mechanical impedance measurement tool,", Journal of Manipulative and Physiological Therapeutics , vol. 22, no. 2, pp. 75-86, 1999.

[8] G. N. Kawchuk, N. G. Prasad, R. C. McLeod, T. Liddle, T. Li, Q. Zhu, "Variability of force magnitude and force duration in manual and instrument-based manipulation techniques,", Journal of Manipulative and Physiological Therapeutics , vol. 29, no. 8, pp. 611-618, 2006.

[9] H. T. Vernon, M. S. Dhami, T. P. Howley, R. Annett, "Spinal manipulation and beta-endorphin: a controlled study of the effect of a spinal manipulation on plasma beta-endorphin levels in normal males,", Journal of Manipulative and Physiological Therapeutics , vol. 9, no. 2, pp. 115-123, 1986.

[10] R. G. Yates, D. L. Lamping, N. L. Abram, C. Wright, "Effects of chiropractic treatment on blood pressure and anxiety: a randomized, controlled trial,", Journal of Manipulative and Physiological Therapeutics , vol. 11, no. 6, pp. 484-488, 1988.

[11] G. Bakris, M. Dickholtz, P. M. Meyer, G. Kravitz, E. Avery, M. Miller, J. Brown, C. Woodfield, B. Bell, "Atlas vertebra realignment and achievement of arterial pressure goal in hypertensive patients: a pilot study,", Journal of Human Hypertension , vol. 21, no. 5, pp. 347-352, 2007.

[12] J. G. Pickar, "Neurophysiological effects of spinal manipulation,", The Spine Journal , vol. 2, no. 5, pp. 357-371, 2002.

[13] E. Salmon, R. S. J. Frackowiak, "Functional metabolic neuroimaging by positron-emission tomography in man,", Revue Neurologique , vol. 146, no. 8-9, pp. 459-477, 1990.

[14] M.-S. Park, Y.-Y. Sunwoo, K.-S. Jang, Y.-M. Han, M.-W. Kim, L.-S. Maeng, Y.-P. Hong, J.-H. Joo-Hyun, Y.-A. Chung, "Changes in brain FDG metabolism induced by acupuncture in healthy volunteers,", Acta Radiologica , vol. 51, no. 8, pp. 947-952, 2010.

[15] X. Duan, M. Tashiro, D. Wu, T. Yambe, Q. Wang, T. Sasaki, K. Kumagai, Y. Luo, S.-I. Nitta, M. Itoh, "Autonomic nervous function and localization of cerebral activity during lavender aromatic immersion,", Technology and Health Care , vol. 15, no. 2, pp. 69-78, 2007.

[16] Y. Ouchi, T. Kanno, H. Okada, E. Yoshikawa, T. Shinke, S. Nagasawa, K. Minoda, H. Doi, "Changes in cerebral blood flow under the prone condition with and without massage,", Neuroscience Letters , vol. 407, no. 2, pp. 131-135, 2006.

[17] R. P. Lystad, H. Pollard, "Functional neuroimaging: a brief overview and feasibility for use in chiropractic research,", Journal of the Canadian Chiropractic Association , vol. 53, no. 1, pp. 59-72, 2009.

[18] T. Ogura, M. Tashiro, M. Masud, S. Watanuki, K. Shibuya, K. Yamaguchi, M. Itoh, H. Fukuda, K. Yanai, "Cerebral metabolic changes in men after chiropractic spinal manipulation for neck pain,", Alternative Therapies in Health and Medicine , vol. 17, no. 6, pp. 12-17, 2011.

[19] L. Ossewaarde, G. A. van wingen, S. C. Kooijman, T. Bäckström, G. Fernández, E. J. Hermans, "Changes in functioning of mesolimbic incentive processing circuits during the premenstrual phase,", Social Cognitive and Affective Neuroscience , vol. 6, no. 5, pp. 612-620, 2011.

[20] J. M. Andreano, L. Cahill, "Menstrual cycle modulation of medial temporal activity evoked by negative emotion,", NeuroImage , vol. 53, no. 4, pp. 1286-1293, 2010.

[21] J.-C. Dreher, P. J. Schmidt, P. Kohn, D. Furman, D. Rubinow, K. F. Berman, "Menstrual cycle phase modulates reward-related neural function in women,", Proceedings of the National Academy of Sciences of the United States of America , vol. 104, no. 7, pp. 2465-2470, 2007.

[22] R.-J. Hwang, C.-H. Wu, L.-F. Chen, T.-C. Yeh, J.-C. Hsieh, "Female menstrual phases modulate human prefrontal asymmetry: a magnetoencephalographic study,", Hormones and Behavior , vol. 55, no. 1, pp. 203-209, 2009.

[23] C.-H. Tu, D. M. Niddam, H.-T. Chao, R.-S. Liu, R.-J. Hwang, T.-C. Yeh, J.-C. Hsieh, "Abnormal cerebral metabolism during menstrual pain in primary dysmenorrhea,", NeuroImage , vol. 47, no. 1, pp. 28-35, 2009.

[24] A. W. Fuhr, J. M. Menke, "Activator methods chiropractic technique,", Topics in Clinical Chiropractic , vol. 9, no. 3, pp. 30-43, 2002.

[25] T. S. Keller, C. J. Colloca, A. W. Fuhr, "In vivo transient vibration assessment of the normal human thoracolumbar spine,", Journal of Manipulative and Physiological Therapeutics , vol. 23, no. 8, pp. 521-530, 2000.

[26] K. J. Friston, C. D. Frith, P. F. Liddle, R. S. J. Frackowiak, "Comparing functional (PET) images: the assessment of significant change,", Journal of Cerebral Blood Flow and Metabolism , vol. 11, no. 4, pp. 690-699, 1991.

[27] K. J. Friston, K. J. Worsley, R. S. J. Frackowiak, J. C. Mazziotta, A. C. Evans, "Assessing the significance of focal activations using their spatial extent,", Human Brain Mapping , vol. 1, no. 3, pp. 210-220, 1994.

[28] K. J. Friston, J. Ashburner, C. D. Frith, J. -B. Poline, J. D. Heather, R. S. J. Frackowiak, "Spatial registration and normalization of images,", Human Brain Mapping , vol. 3, no. 3, pp. 165-189, 1995.

[29] J. W. Papez, "A proposed mechanism of emotion,", Archives of Neurology And Psychiatry , vol. 38, no. 4, pp. 725-743, 1937.

[30] G. Bush, P. Luu, M. I. Posner, "Cognitive and emotional influences in anterior cingulate cortex,", Trends in Cognitive Sciences , vol. 4, no. 6, pp. 215-222, 2000.

[31] A. J. Shackman, T. V. Salomons, H. A. Slagter, A. S. Fox, J. J. Winter, R. J. Davidson, "The integration of negative affect, pain and cognitive control in the cingulate cortex,", Nature Reviews Neuroscience , vol. 12, no. 3, pp. 154-167, 2011.

[32] P. Petrovic, E. Kalso, K. M. Petersson, M. Ingvar, "Placebo and opioid analgesia--imaging a shared neuronal network,", Science , vol. 295, no. 5560, pp. 1737-1740, 2002.

[33] L. Lindgren, G. Westling, C. Brulin, S. Lehtipalo, M. Andersson, L. Nyberg, "Pleasant human touch is represented in pregenual anterior cingulate cortex,", NeuroImage , vol. 59, no. 4, pp. 3427-3432, 2012.

[34] I. E. De Araujo, E. T. Rolls, M. I. Velazco, C. Margot, I. Cayeux, "Cognitive modulation of olfactory processing,", Neuron , vol. 46, no. 4, pp. 671-679, 2005.

[35] F. Grabenhorst, E. T. Rolls, B. A. Parris, A. A. D'Souza, "How the brain represents the reward value of fat in the mouth,", Cerebral Cortex , vol. 20, no. 5, pp. 1082-1091, 2010.

[36] I. Tracey, "Getting the pain you expect: mechanisms of placebo, nocebo and reappraisal effects in humans,", Nature Medicine , vol. 16, no. 11, pp. 1277-1283, 2010.

[37] A. Shenhav, J. D. Cohen, M. M. Botvinick, "Dorsal anterior cingulate cortex and the value of control,", Nature Neuroscience , vol. 19, no. 10, pp. 1286-1291, 2016.

[38] E. A. Moulton, J. D. Schmahmann, L. Becerra, D. Borsook, "The cerebellum and pain: passive integrator or active participator?,", Brain Research Reviews , vol. 65, no. 1, pp. 14-27, 2010.

[39] B. Sacchetti, B. Scelfo, P. Strata, "Cerebellum and emotional behavior,", Neuroscience , vol. 162, no. 3, pp. 756-762, 2009.

[40] P. Hernáez-Goñi, J. Tirapu-Ustárroz, L. Iglesias-Fernández, P. Luna-Lario, "The role of the cerebellum in the regulation of affection, emotion and behaviour,", Revista de Neurologia , vol. 51, no. 10, pp. 597-609, 2010.

[41] H. C. Lou, T. W. Kjaer, L. Friberg, G. Wildschiodtz, S. Holm, M. Nowak, "A ¹⁵ O-H₂ O PET study of meditation and the resting state of normal consciousness,", Human Brain Mapping , vol. 7, no. 2, pp. 98-105, 1999.

[42] H. D. Critchley, R. N. Melmed, E. Featherstone, C. J. Mathias, R. J. Dolan, "Brain activity during biofeedback relaxation a functional neuroimaging investigation,", Brain , vol. 124, no. 5, pp. 1003-1012, 2001.

[43] U. M. Nater, N. Rohleder, "Salivary alpha-amylase as a non-invasive biomarker for the sympathetic nervous system: current state of research,", Psychoneuroendocrinology , vol. 34, no. 4, pp. 486-496, 2009.

[44] D. A. Granger, K. T. Kivlighan, M. El-Sheikh, E. B. Gordis, L. R. Stroud, "Salivary α -amylase in biobehavioral research: recent developments and applications,", Annals of the New York Academy of Sciences , vol. 1098, pp. 122-144, 2007.

[45] U. M. Nater, N. Rohleder, J. Gaab, S. Berger, A. Jud, C. Kirschbaum, U. Ehlert, "Human salivary alpha-amylase reactivity in a psychosocial stress paradigm,", International Journal of Psychophysiology , vol. 55, no. 3, pp. 333-342, 2005.

[46] T. Nakamura, K. Kawahara, M. Kusunoki, T. Hara, I. Minejima, "Involvement of the nervous system in the control of glucose uptake in contracting hindlimb muscle in the rat,", Hormone and Metabolic Research , vol. 29, no. 12, pp. 593-598, 1997.

[47] N. J. Goddard, P. T. Gosling, "Intra-articular fluid pressure and pain in osteoarthritis of the hip,", The Journal of Bone & Joint Surgery--British Volume , vol. 70, no. 1, pp. 52-55, 1988.