The effects of chiropractic spinal manipulation on central processing of tonic pain - a pilot study using standardized low-resolution brain electromagnetic tomography (sLORETA).

The objectives of the study were to investigate changes in pain perception and neural activity during tonic pain due to altered sensory input from the spine following chiropractic spinal adjustments. Fifteen participants with subclinical pain (recurrent spinal dysfunction such as mild pain, ache or stiffness but with no pain on the day of the experiment) participated in this randomized cross-over study involving a chiropractic spinal adjustment and a sham session, separated by 4.0 ± 4.2 days. Before and after each intervention, 61-channel electroencephalography (EEG) was recorded at rest and during 80 seconds of tonic pain evoked by the cold-pressor test (left hand immersed in 2 °C water). Participants rated the pain and unpleasantness to the cold-pressor test on two separate numerical rating scales. To study brain sources, sLORETA was performed on four EEG frequency bands: delta (1-4 Hz), theta (4-8 Hz), alpha (8-12 Hz) and beta (12-32 Hz). The pain scores decreased by 9% after the sham intervention (p < 0.05), whereas the unpleasantness scores decreased by 7% after both interventions (p < 0.05). sLORETA showed decreased brain activity following tonic pain in all frequency bands after the sham intervention, whereas no change in activity was seen after the chiropractic spinal adjustment session. This study showed habituation to pain following the sham intervention, with no habituation occurring following the chiropractic intervention. This suggests that the chiropractic spinal adjustments may alter central processing of pain and unpleasantness.

RCT Entities:   Population Interventions Outcomes

Introduction

Changes in the way the human brain processes pain, as well as the capacity to modulate the pain experience, underlies the pathogenesis of most chronic pain conditions[1]. Typically, the perception of pain is induced by a potential damaging stimulus in the periphery, which activates peripheral nerves. The pain signal is transmitted to the spinal cord and further to deep centers within the brain. However, the pain system is not hard-wired, but rather a complex dynamic process with advanced modulatory properties, such as activation of descending pathways that can modulate pain perception, and that is also responsible for the recurrence and chronicity of the pain experience[2-4]. During the last two decades of research, there has been a growing interest to understand the underlying mechanisms of the descending inhibition of pain. Pain modulating processes in humans are present at multiple levels of the human nervous system; at the spinal level in the dorsal horn[5], in the brainstem, including neurons in the rostral ventromedial medulla and periaqueductal grey (RVM-PAG)[6], and at the cortical level[7]. Cognitive and affective processes within the limbic system are also known to influence the RVM-PAG system, and are also involved in the top-down cortico-limbic-brainstem inhibition[8,9]. In the literature, these kinds of central nervous system modulations of the pain experience are collectively known as conditioned pain modulation (CPM).

Chiropractic is one option that is utilized clinically to help people who suffer with musculoskeletal pain, both acute and chronic. Multiple clinical and basic science studies have shown that chiropractic care can reduce pain[10-12], however the mechanisms involved are not well understood[13-16]. Chiropractic care is a holistic approach to health with a particular focus on the relationship between the spine and nervous system[17]. Traditionally, the main focus of chiropractic care has been the location, analysis and correction of vertebral subluxations[17]. Vertebral subluxations have been defined as ‘a self- perpetuating, central segmental motor control problem that involves a joint, such as a vertebral motion segment, that is not moving appropriately, resulting in ongoing maladaptive neural plastic changes that interfere with the central nervous system’s ability to self-regulate, self-organize, adapt, repair and heal’[18]. Chiropractors use a combination of indicators of spinal dysfunction to identify and characterize vertebral subluxations[19]. They then use a variety of manual techniques, the most common being specific high-velocity, low amplitude adjustments, to correct the vertebral subluxation[20,21]. This is sometimes referred to as spinal manipulation in the research literature, however, spinal manipulation is also used in the literature for other forms of non-specific thrusts delivered to spinal regions as opposed to specific dysfunctional vertebrae[22]. Research has shown that chiropractic spinal adjustments alter the afferent input from the spine which leads to changes in central nervous system (CNS) function[23]. It has been postulated that spinal adjustments stimulates the dysfunctional paraspinal tissues, which modulates central processing and potentially sensitization, and therefore alters the perception of pain cortically, in a ‘top-down’ manner or both[13,23-25]. The majority of the studies which have shown neuroplastic effects following a chiropractic intervention have focused on changes in somatosensory processing, motor control and functional performance[26-29]. These neural plastic changes have been assessed with both spinal and cortical measures[30,31], however, knowledge about the effects of spinal adjustments on pain and neural markers for altered CNS function due to this alteration in sensory input from the spine is still lacking[14,32].

As it is becoming clear that chiropractors impact brain function consistently, it is very likely that chiropractic care influences not only the biomechanical movement patterns of the spine and improves proprioceptive processing of the spine[33,34] but also directly impacts the so called ‘pain matrix’ in the brain and thus has a CPM effect on a person’s perception of pain. Chiropractic care may alter the way a person ‘feels’ pain and may therefore help chronic musculoskeletal pain sufferers by improving the top-down cortico-limbic-brainstem inhibition of pain. Alternatively, chiropractic spinal adjustments may also alter pain perception via the neuroendocrine system, or segmental inhibition at the spinal cord level, or even motivational-affective modulation[10,13,23,24,32]. Further research is therefore required to elucidate the basic neurophysiological mechanisms of spinal adjustments and their effect on pain perception[13,14,32].

Tonic pain models can be used in experiments to study pain, as experimental tonic pain produces a similar sensory experience to that of a clinical setting[35-37]. The cold-pressor (CP) test has been shown to be a robust stimulus to induce tonic pain[38] and mimics clinical pain due to its high level of unpleasantness[39]. Using the CP test together with electroencephalography (EEG), Gram et al. showed that EEG spectral indices change in parallel with pain intensity[38], further validating this method to assess central processing of tonic pain. Shao et al., using EEG and standardized low-resolution brain electromagnetic tomography (sLORETA), showed that tonic cold pain induced significant changes of source power across different frequency bands in many brain regions, including the prefrontal, primary and secondary somatosensory (S1 and S2), insular and cingulate cortices[40]. These same regions have been consistently shown to be activated or deactivated during pain perception by various neuroimaging studies using functional magnetic resonance imaging (fMRI), magnetoencephalography (MEG) and positron emission tomography (PET) (see e.g.[7,41-43]). Specifically, Shao et al. found that the oscillatory activities that significantly correlated with subjective pain ratings were found in the prefrontal and cingulate regions[40]. The authors recommended the testing of pain treatment methods, especially to explore whether they could alter the neural oscillations in the particular frequency bands in the specific brain regions they had identified in their study. These specific brain regions and frequency bands identified during the CP test were the 4–8 Hz (theta) oscillatory activity in the prefrontal cortex, the 8–12 Hz (alpha) oscillatory activity in the anterior cingulate cortex (ACC) and the 12–18 Hz (beta) oscillatory activity in the posterior cingulate cortex (PCC).

In the present study, sLORETA was therefore used for source localization of EEG recorded pre and post chiropractic spinal adjustments and a sham intervention. The sLORETA is a linear inverse algorithm that estimates the 3D distribution of the cortical generators of the EEG and provides the lowest localization error compared to several other linear inverse algorithms[44]. The sLORETA was done on four frequency bands: (i) delta (1–4 Hz), (ii) theta (4–8 Hz), (iii) alpha (8–12 Hz), and (iv) beta (12–32 Hz). These frequency bands were chosen to cover the frequency changes noted by Shao et al.[40] and because they have been the most commonly reported with LORETA analysis in frequency domain[45-47].

The hypothesis of the current cross-over study was that a single session of chiropractic spinal adjustments would alter the neural response to tonic pain, in comparison to a sham spinal adjustment session. Therefore, the aims of this pilot study were to:

Assess whether there are differences in subjective pain perception and unpleasantness levels during the CP test following sham and chiropractic interventions, and

Compare the brain activity (sLORETA solution) underlying different frequency bands of the CP-EEG data before and after each intervention.

Methods

The study was conducted according to the Declaration of Helsinki. The North Denmark Region Committee on Health Research Ethics approved the study (N-20150033). The study was one-way blinded; therefore, the participants did not know which type of intervention they received. This study was registered retrospectively on 24th August 2018 with the Australian New Zealand Clinical Trial Registry (trial registration number ACTRN12618001420235). The study design is given in Supplementary Figure S1.

Subjects

Fifteen subclinical pain subjects (10 males, (mean ± SD) age = 32.1 ± 7.2 years, BMI = 24.02 ± 5.94 kg/m2) participated in the study. Subclinical pain refers to recurrent spinal dysfunction such as mild pain, ache or stiffness for which treatment is not yet sought, and most importantly no pain on the day of experimental assessment, to avoid the confounding effects of altered resting pain levels. All subjects gave their written informed consent to participate in the study.

Before entering the study, the subjects were introduced to the lab environment and were assessed by a chiropractor with 15 years of experience to assure that the subjects pass the eligibility criteria to participate in the study. Subjects were included if they were aged between 18 and 50 years, had a history of recurring spinal dysfunction such as mild pain, ache or stiffness without a history of known trauma. Subjects were ineligible to participate if they exhibited actual pain on the day(s) of the experiment, had no evidence of spinal dysfunction, had absolute contraindications to chiropractic spinal adjustments, had experienced previous significant adverse reactions to chiropractic care or spinal manipulation, or if they had sought treatment for their pain symptoms. The subjects were also required to be fluent in the understanding of written and spoken English to participate in the study.

Experimental protocol

The subjects participated in two experimental sessions; sham and chiropractic, separated by 4.0 ± 4.2 days at Aalborg University Hospital. Each session consisted of a 1-minute resting state EEG recording, followed by a CP test during which EEG was also recorded. The resting EEG was recorded in order to observe the central processing of tonic pain itself (i.e. difference between pre-intervention CP-EEG and pre-intervention resting EEG) and this would subsequently aid in interpreting the changes in the central processing of tonic pain following the interventions. Figure 1 shows the overview of the experiment.

Figure 1

Methodology overview. (A) Shows the order and duration of the resting state and cold- pressor test during which EEG was also recorded. (B) Shows the length and position of EEG used for analysis after it was truncated. The circle in (C) shows the time point when the participant rated the pain and unpleasantness to cold-pressor. (D) Shows the pairs used for comparisons in sLORETA. The intervention lasted approximately 10 minutes. Abbreviations: sLORETA = standardized low-resolution brain electromagnetic tomography.

Interventions

The chiropractic spinal adjustment intervention and sham intervention were similar to those used in previous studies that have investigated the neurophysiological effects of chiropractic spinal adjustments[23,30,31]. The two researchers who carried out the interventions made a concerted effort to maintain the same levels of language and professional dealing with the participants. We advertised about the recruitment on a Facebook page. All the participants (upon questioning) were completely naïve to chiropractic care, so had no previous experience with a chiropractor, and therefore did not know what to expect. Participants were blinded to nature of intervention (sham or adjustment intervention). At the end of the second session, before telling them which session a real chiropractic intervention was, and which one was a sham intervention, the participants were asked whether they thought one session could have been a sham. Out of the 15 subjects, three subjects said that one session felt like a sham. Out of these three, two were correct and one was wrong, as this third person thought the sham was real and the real session was a sham.

Chiropractic spinal adjustment intervention

The chiropractic spinal adjustments carried out were high-velocity, low-amplitude thrusts to the spine or pelvic joints, which is a standard adjustment technique used by chiropractors and is also referred to as spinal manipulation. The sites chosen for the spinal adjustments were based on the clinical indicators of spinal and pelvic joint dysfunction[19], which were: tenderness to palpation of the relevant joints; manual palpation for restricted intersegmental range of motion; palpable asymmetric intervertebral muscle tension, and any abnormal or blocked joint play and end-feel of the joints. These indicators have been shown to be reliable when used in combination to identify level of subluxation in the spine[48]. Multiple levels of the spine were adjusted in each participant by a registered chiropractor during the chiropractic adjustment session.

Sham intervention

The sham intervention acted as a physiological control. During the sham intervention, the chiropractor simulated a spinal adjustment session, which included passive and active movements of the participants head, spine, and body, similar to what is done during the chiropractic adjustment session. Care was taken during the sham session to not take any joints to their end range of motion or to cause a cavitation in the spine. This was done to limit the afferent barrage from large diameter afferents in the paraspinal muscles to the central nervous system, while controlling for the body movements, touch and vestibular input associated with setting up to provide chiropractic adjustments.

Cold-pressor test

The CP test was performed using a circulating water bath (Grant, Fischer Scientific, Slangerup, Denmark). The water was cooled to 2 °C and the subjects immersed their left hand in the water up to the wrist for 80 seconds (Fig. 1A).

Pain and unpleasantness scores

The participants rated their pain and unpleasantness on two separate numerical rating scales after their hand had been in the water for 60 seconds (Fig. 1C). The two scales ranged from 0 (no pain/unpleasantness) to 10 (maximum pain/unpleasantness), and the experimenter noted down the scores.

EEG

EEG was recorded at a sampling rate of 1000 Hz in a dimly lit room using a 61-channel cap (MEQNordic A/S, Jyllinge, Denmark) and Synamp system (Neuroscan Compumedics, El Paso, TX, USA). The reference electrode was just above AFz. During the recordings, the participants lay in a supine position and were instructed to relax, keep their focus on a point and reduce eye blinking.

The EEG was preprocessed offline. The following preprocessing steps were performed using Neuroscan 4.3.1 (Neuroscan, El Paso, TX, USA): (1) noisy channels were interpolated using their neighboring channels; (2) afterward, a 50 Hz notch filter was applied, (3) followed by a band-pass filter of 1 to 70 Hz. For the remaining preprocessing steps, MATLAB 2015b (The MathWorks, Inc., Natick, MA, USA.) was used: (4) EEG was truncated for further analysis: the resting EEG from 2 to 58 s was taken, and for the CP part, 72 s of EEG were taken, starting at 8 s from the onset of the stimulus (Fig. 1B). The first 8 seconds of CP-EEG data were removed to avoid artifacts including muscle contractions caused by the immediate unpleasantness after immersing the hand into the cold water[38]; (5) to reduce the computational load in the sLORETA matrix calculations, the filtered EEG was downsampled by a factor of 4; and finally, (6) for obtaining smooth power spectral density, the EEG was divided into epochs with length of 8 s to facilitate the averaging procedure in sLORETA.

sLORETA

The underlying sources of the EEG were estimated using the sLORETA software package, version 20151222[44] (available at http://www.uzh.ch/keyinst/loreta). The sLORETA was done in the frequency domain to localize neural oscillators on the average referenced EEG. The EEG was average referenced in the sLORETA software. Cross-spectral matrices for each subject were computed in sLORETA software for four frequency bands: (i) delta (1–4 Hz), (ii) theta (4–8 Hz), (iii) alpha (8–12 Hz), and (iv) beta (12–32 Hz). The cross-spectral matrices for each participant were then averaged as the input for sLORETA source analysis. The sLORETA software was used to estimate the statistical differences in brain activity (in the four EEG frequency bands) between:

The baselines of both experimental sessions for (i) resting state EEG and (ii) CP-EEG to make sure there were no unexpected differences,

The baseline CP-EEG and baseline resting state EEG to observe how the brain processes tonic pain (Fig. 1D),

The post-intervention and baseline EEG for (i) resting state and (ii) CP to find the effect of each intervention on neural activity (Fig. 1D).

The cortical gray matter was divided into 6239 voxels with a resolution of 5 mm3. The head model and electrode coordinates according to Montreal Neurological Institute average MRI brain map (MNI-152)[49] were used.

Spectral analysis

The EEG power spectral analysis was performed using FieldTrip toolbox[50] to assist the sLORETA results. To estimate the differences in the power spectrum of the four EEG frequency bands between the conditions mentioned in the section above, the power spectra between 1 and 32 Hz of the EEG were calculated using Fourier basis with an Hanning window of 1 s, which was followed by computation of the average power of each frequency band.

Statistics

The data are presented as mean ± SD unless otherwise indicated. The statistical significance threshold was p < 0.05.

Two-way repeated measures analysis of variance (ANOVA) was performed to identify the changes in the pain and unpleasantness scores with time (before and after) and intervention (sham and chiropractic) as the two factors. If the overall significance in the ANOVA test was found, all pairwise multiple comparisons procedures (Student-Newman-Keuls Method) were performed in order to assess where the differences were. The software package SigmaStat version 3.0 (SPSS Inc. Chicago, IL, USA) was used for the above statistical analysis.

The statistical analysis for source localization was done using the sLORETA software’s built-in statistics tool using statistical non-parametric mapping[51] which adjusted for multiple comparisons by utilizing Fisher’s random permutation test with 5000 randomizations. To compare current sources in different frequency bands, paired two-tailed Student’s t-test was used to compare the baselines of chiropractic and sham sessions; the baseline CP and baseline resting state; and the post-sessions and baselines.

Non-parametric cluster-based permutation test[52] was used to identify the differences in the EEG power spectrum between the baselines of chiropractic and sham sessions; the baseline CP and baseline resting state; and the post-sessions and baselines. The clusters were defined as two or more continuous channel-power pairs each with p < 0.05 from the paired two-tailed t-test with respect to the conditions. The t-values within each cluster were added to get the cluster-level statistics and the maximum of cluster-level statistics was used as the test statistic. A cluster was considered significant if its Monte Carlo probability for each tail exceeded the threshold of 0.025 compared to the reference distribution approximated by Monte Carlo method with 5000 permutations.

Results

All fifteen of the enrolled subjects successfully completed the experiment and data from all subjects were used for the analysis.

The pain and unpleasantness scores of the CP test are summarized in Table 1. The time effects on the pain scores (F = 6.7, p < 0.05) and unpleasantness scores (F = 9.6, p < 0.05) were significant. The posthoc test revealed that the pain scores decreased after the sham intervention (p < 0.05), whereas the unpleasantness scores decreased after both interventions (both p < 0.05). There were no significant interactive effects present between the time and intervention.

Table 1

Pain and unpleasant scores.

Scores	Sham Intervention		Spinal Adjustment Intervention
	Before	After	Before	After
Pain	8.33 ± 1.54	7.53 ± 1.73*	8.27 ± 1.39	7.77 ± 1.88
Unpleasantness	8.73 ± 1.62	8.07 ± 2.19*	8.53 ± 2.07	7.90 ± 2.45*

*Significant differences (p < 0.05).

Effects on source location

The sLORETA analysis showed no differences in both the resting state EEG baselines and the CP-EEG baselines in all frequency bands (all p > 0.05).

The comparison between baseline CP-EEG and baseline resting state EEG showed a widespread increase in cortical activity during pain compared to resting state in all frequency bands. The results from the sLORETA analysis are summarized in Table 2 and Fig. 2.

Table 2

sLORETA localized EEG cortical sources with significant differences between the baseline cold-pressor and baseline resting state.

Brain structure (Brodmann Area)	Delta (1–4 Hz)	Theta (4–8 Hz)	Alpha (8–12 Hz)	Beta (12–32 Hz)
Angular Gyrus (39)	24	19	—	5
Anterior Cingulate (24, 25, 32, 33)	114	—	61	57
Cingulate Gyrus (23, 24, 31, 32)	192	—	—	38
Cuneus (17, 18, 19, 23, 30)	161	—	—	49
Extra-Nuclear (13)	—	—	10	—
Fusiform Gyrus (18, 19, 20, 21, 32, 37)	119	70	148	87
Inferior Frontal Gyrus (9, 10, 11, 13, 44, 45, 46, 47)	185	25	295	83
Inferior Occipital Gyrus (18, 19)	18	18	28	5
Inferior Parietal Gyrus (40)	—	—	—	25
Inferior Parietal Lobule (39, 40)	183	54	15	—
Inferior Temporal Gyrus (20, 37)	42	36	93	19
Insula (13)	68	22	64	35
Lingual Gyrus (17, 18, 19)	143	—	121	70
Medial Frontal Gyrus (6, 8, 9, 10, 11, 25, 32, 46)	97	—	92	39
Middle Frontal Gyrus (6, 8, 9, 10, 11, 46, 47)	276	95	121	13
Middle Occipital Gyrus (18, 19, 37)	68	51	65	—
Middle Temporal Gyrus (12, 21, 22, 37, 38, 39)	103	84	175	44
Orbital Gyrus (11, 47)	23	—	30	12
Parahippocampal Gyrus (19, 27, 28, 30, 34, 35, 36, 37)	172	56	101	139
Postcentral Gyrus (1, 2, 3, 40, 43)	116	32	6	—
Posterior Cingulate (23, 29, 30, 31)	85	—	13	60
Precentral Gyrus (4, 6, 9, 43, 44)	115	56	48	—
Precuneus (7, 9, 31)	55	—	—	69
Rectal Gyrus (11)	39	—	44	36
Subcallosal Gyrus (25, 34)	18	—	18	17
Sub-Gyral (6, 20, 21)	6	—	17	—
Superior Frontal Gyrus (6, 8, 9, 10, 11)	171	10	—	10
Superior Occiptial Gyrus (17, 19)	7	6	—	—
Superior Parietal Lobe (7)	—	6	—	—
Superior Parietal Lobule (7)	8	—	—	23
Superior Temporal Gyrus (19, 22, 38, 39, 41, 42)	146	78	229	86
Supramarginal Gyrus (40)	45	24	21	—
Transverse Temporal Gyrus (41, 42)	10	—	5	5
Uncus (20, 28, 34, 36, 38)	64	19	64	57

The number of voxels with significant power changes (p < 0.05) is listed.

Figure 2

sLORETA CP vs Rest. Slice views of source locations with the changes in activity (t-values) between baseline cold-pressor and resting EEG for each frequency band. A significant increase (p < 0.05) in activity (yellow colored) in all frequency bands can be seen during the pain condition. Abbreviations: CP = cold-pressor; sLORETA = standardized low-resolution brain electromagnetic tomography.

Neither of the two interventions changed the resting state EEG (all p > 0.05).

The brain activity underlying the CP test decreased following the sham intervention (Fig. 3 and Table 3) but showed no differences following the spinal manipulation. Following the sham session, the most significant decreases (p < 0.05) in activity due to CP were seen in the delta (Brodmann area 32, cingulate gyrus, limbic lobe) and alpha (Brodmann area 42, transverse temporal gyrus, temporal lobe) bands, whereas marginally significant decrease (p = 0.05) were seen in the theta (Brodmann area 9, medial frontal gyrus, frontal lobe) and beta (Brodmann area 22, superior temporal gyrus, temporal lobe) bands. The complete list of brain regions showing changes in activity can be seen in Table 3.

Figure 3

sLORETA CP(Treatment) vs CP (Baseline). Slice views of source locations with the changes in activity (t-values) during the cold-pressor test for each frequency band after sham and spinal adjustment sessions compared to respective baseline activity. The decrease in activity in all frequency bands can be seen after the sham session with a significant decrease (p < 0.05) in activity (magenta colored) in delta and alpha bands; and marginally significant decreases (p = 0.05) in the theta and beta bands. Abbreviations: CP = cold-pressor; sLORETA = standardized low-resolution brain electromagnetic tomography.

Table 3

sLORETA localized EEG cortical sources during the cold-pressor test with significant differences after sham session compared to baseline activity.

Brain structure (Brodmann Area)	Delta (1–4 Hz)	Theta (4–8 Hz)	Alpha (8–12 Hz)	Beta (12–32 Hz)
Cingulate Gyrus (24, 32)	5	—	—	—
Fusiform Gyrus (42)	—	—	1	—
Inferior Parietal Lobule (20)	—	—	1	—
Inferior Temporal Gyrus (20, 21, 22, 40)	—	—	4	—
Medial Frontal Gyrus (6, 9, 32)	4	1	—	—
Middle Frontal Gyrus (6)	7	—	—	—
Middle Temporal Gyrus (20, 21, 22, 40, 42)	—	—	9	—
Precuneus (31)	—	—	—	1
Sub-Gyral (6, 31)	1	—	—	1
Superior Frontal Gyrus (6)	10	—	—	—
Superior Temporal Gyrus (20, 22, 40)	—	—	7	1
Supramarginal Gyrus (21, 22)	—	—	6	—
Transverse Temporal Gyrus (21)	—	—	1	—

The number of voxels with significant power changes (p < 0.05) is listed for delta and alpha bands; and with marginal significant changes (p = 0.05) for the theta and beta bands.

Effects on power spectrum

Similar to the sLORETA results, the spectral analysis showed no differences in both the resting state EEG baselines (Fig. 4A) and the CP-EEG baselines (Fig. 4B) in all frequency bands (all p > 0.05).

Figure 4

Spectral analysis. Differences in power (μV2/Hz) of the four frequency bands. No differences in (A) resting state and (B) CP-EEG baselines. (C) A significant increase (p < 0.05) in power (asterisks) in all frequency bands can be seen during the pain condition. No effect of both sham (D) and the chiropractic (E) interventions on the resting state. (F) The decrease in power of CP-EEG in all frequency bands can be seen after the sham session with a significant decrease (p < 0.05) in beta band. (G) There was slight increase in the overall power of CP-EEG in all frequency bands except alpha band after the chiropractic session. Abbreviations: CP = cold-pressor; BL = baseline; TR = treatment; Chiro = chiropractic.

The comparison between baseline CP-EEG and baseline resting state EEG showed a widespread increase in EEG power during pain as compared to resting state in all frequency bands (Fig. 4C).

Neither of the two interventions changed the resting state EEG power (Fig. 4D,E) (all p > 0.05).

The EEG power underlying the CP test decreased following the sham intervention (Fig. 4F), with significant decrease in the beta band. Although non-significant for the delta, theta, alpha band, the trend of decreased brain activity was similar to the sLORETA results. The EEG power underlying the CP test showed a slight (non-significant) increase following the spinal manipulation (Fig. 4G) in all frequency bands except the alpha band.

Discussion

In this study, we investigated the effects of altering sensory input from dysfunctional areas of the spine (by spinal adjustments) on pain perception and the neural activity (in terms of source localization of EEG) obtained during the cold-pressor test. The pain scores due to the cold-pressor test decreased by 9% after the sham session, whereas the unpleasantness scores decreased by 7% after both interventions. The brain activity associated with the cold-pressor test following the sham session decreased significantly in the delta and alpha bands; and showed a marginally significant decrease in theta and beta bands, whereas there were no changes in EEG source localization following the chiropractic spinal adjustment session.

EEG activity during tonic pain

This study found increased neural activity in all frequency bands when CP was compared with resting state. The increase in cortical activity due to tonic pain is supported by other studies using EEG[53-55], PET[56] and fMRI[57,58]. In this study, the brain regions with a change in activity included the insula and anterior cingulate cortex, which are among the most often reported active regions during pain perception in many neuroimaging studies using fMRI and PET[59,60]. Hence, the insula and anterior cingulate cortex likely have an important role in the processing of pain. We found increased delta, theta and beta activities, which is consistent with findings of studies utilizing tonic pain in healthy volunteers[40,53-55,61]. The increase of the delta, theta and beta activities in the anterior cingulate and insula cortices, among other regions in this study, likely implies negative feelings to the pain induced by CP, as it has been reported by many studies[58,60,62,63] that these regions are associated with emotional aspects of pain processing.

There was an increase in alpha activity during the CP test. Alpha oscillations have been found to be increased over frontal or parieto-occipital regions during tonic pain using the CP test[53] and hypertonic saline injection[61]. On the other hand, alpha power and underlying neural activity have also been shown to decrease following tonic pain[37,40,64]. The type or intensities of stimulus used in these studies may be the reason for inconsistency with the current study. For example, Shao et al. used a 10 °C CP test[40], Babiloni et al. used CO2-laser stimulation[64] and heat stimulation was used by Nir et al.[37]. The changes in the alpha EEG band are associated with attention processes[65] and anticipation to pain[64]. Therefore, the increase of source activity underlying the alpha band in this study is likely due to the subjects’ attention to pain.

After the sham intervention there was a decrease in pain scores during the tonic pain. This is not an unexpected result, since habituation to pain is a normal reaction in humans and animals to continuous or repetitive painful stimuli, which decreases the perceived pain and pain-related responses[66-68]. There was also a decrease in the neural activity underlying EEG during tonic pain following the sham intervention, indicating the pain scores related to the underlying activity in these specific brain regions. The regions which showed the most significant decrease were: cingulate gyrus, limbic lobe (delta) and transverse temporal gyrus, temporal lobe (alpha). Marginally significant decrease was seen in the medial frontal gyrus, frontal lobe (theta) and superior temporal gyrus, temporal lobe (beta). The areas that showed a decrease due to the CP test are a subgroup of brain areas, which were also activated due to the CP test itself. Hence, the decrease of activity in these brain areas likely has a role in the inhibition of pain, as seen in the decreased pain scores. In this study, there were about 10 minutes between the first and second EEG recording during tonic pain. The after-effects of tonic pain can last up to 30 minutes in humans[69,70] and therefore, the decreased cortical activity, in combination with the decreased pain/unpleasantness perception of the CP stimulus after the sham intervention, can most likely be attributed to central habituation to the CP-induced pain. It is possible that there is some form of placebo effect occurring as well, as most subjects in this study were novices to chiropractic and therefore did not know what to expect.

Although habituation to pain is a normal reaction in humans and animals to continuous or repetitive pain stimuli[66-68], it is a multifactorial event that is not well understood[71]. It has been hypothesized that habituation actually involves dual competing processes of depression (habituation) and facilitation (sensitization), that combine to give a behavioral or perceptual outcome[72]. Habituation to pain may lead to maladaptive plastic changes in the neural system. Maladaptive neuroplasticity has been reported to be induced by chronic pain[73]. These maladaptive neuroplastic changes cause individual sufferers to experience symptoms and functional disturbance, rather than the pain itself[74-76].

Spinal adjustment intervention

In the present study, there were no differences in neural activity due to CP-induced pain following the spinal adjustment session. This may be due to the spinal adjustments having an impact on pain habituation. It is possible that the altered afferent input from the spine following spinal adjustments modulates the interplay between the dual competing processes of depression (habituation) and facilitation (sensitization). Previous studies that have investigated the hypoalgesic effects of spinal manipulation on temperature-induced pain have produced conflicting results[14]. Millan et al. suggested that this may be related to the type of pain fibers (C-fiber or A-delta fibers) that are stimulated during different testing protocols, with C-fiber mediated pain being more influenced by a spinal manipulation intervention than A-delta fiber mediated pain[14]. Previous studies have also suggested that spinal manipulation may reduce the central sensitization of pain[16,25]. This seems somewhat paradoxical based on the findings of the present study because a lessening of sensitization would be expected to result in an increase in habituation, as opposed to a decrease[72]. It is possible that instead of simply decreasing central sensitization, spinal adjustments may in some cases ‘reset’ the facilitatory and inhibitory processes associated with habituation. This may be due to spinal adjustments resulting in altered afferent paraspinal tissue input that affects the manner in which the somatosensory cortex integrates subsequent afferent information, as has been previously hypothesized[10,23,25,27]. For example, it has been shown that spinal adjustments alter processing in the prefrontal cortex[31]. The prefrontal cortex has been shown to impact the degree to which the insular cortex is activated during repeated cold pressor stimulation, thus can alter the way in which a person’s brain habituates to a cold-pain stimulus[77]. Gaining a greater understanding of these potential mechanisms and processes may be important when considering the effects of chiropractic care on both acute and chronic pain.

Study considerations

The sham and experimental interventions were separated by up to 10 days. This may induce some time effects in the data. Furthermore, this was a pilot study with 15 subjects where the neural response to tonic pain was assessed after a single session of spinal adjustments. We recognize that the number of subjects was not high. However, due to the little variation in response and the high sensitivity of the experimental models, the numbers are within the normal range for such explorative studies[10,78-82]. To further validate these findings, future studies should look at the effects of chiropractic care on central processing of tonic pain in a larger population and over a longer period of chiropractic care.

It would also be worth considering to include a non-intervention session, as mind-set has been shown to impact habituation to cold-pain stimuli[83]. From the current study design, it is therefore hard to make firm conclusions about the sham only effects. It is impossible to be sure whether our results following the sham were due to the habituation, placebo effect, or both. Having a separate no-intervention control session would have helped elucidate this.

Finally, it would also have been interesting to analyze the gamma band for a study such as this. However, to do that, longer EEG recordings are required to facilitate the removal of the artifacts, especially those related to muscles as the spectrum of EMG overlaps that of the gamma band. The independent component analysis (ICA) can be used for this purpose but it requires approximately 20 to 30 times squared number of channels amount of data points, and it is not easy from a practical (pain) point of view to record this amount of EEG during the CP-test as most of the subjects cannot tolerate pain this long.

Conclusion

This study showed a habituation to pain response following the sham intervention, with no changes in the neural processing of tonic pain following the chiropractic spinal adjustment session. Changes to spinal function with chiropractic spinal adjustments appears to affect the way in which the central nervous system responds to repeated pain stimuli. However, this needs to be further explored before concrete conclusions can be made. Future studies should investigate the long-term effects of chiropractic care on pain processing in sub-clinical pain populations as well as in patients suffering from acute and chronic pain.