Next Article in Journal
Site-Specific Effects of Online rTMS during a Working Memory Task in Healthy Older Adults
Next Article in Special Issue
Differential Changes in Early Somatosensory Evoked Potentials between the Dominant and Non-Dominant Hand, Following a Novel Motor Tracing Task
Previous Article in Journal
Building, Hosting and Recruiting: A Brief Introduction to Running Behavioral Experiments Online
Previous Article in Special Issue
Adolescent Awkwardness: Alterations in Temporal Control Characteristics of Posture with Maturation and the Relation to Movement Exploration
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigating the Effects of Chiropractic Spinal Manipulation on EEG in Stroke Patients

1
Mech-Sense, Department of Gastroenterology and Hepatology, Aalborg University Hospital, 9000 Aalborg, Denmark
2
Department of Clinical Medicine, Aalborg University, 9000 Aalborg, Denmark
3
Centre for Chiropractic Research, New Zealand College of Chiropractic, Auckland 1060, New Zealand
4
Faculty of Health & Environmental Sciences, Health & Rehabilitation Research Institute, AUT University, Auckland 0627, New Zealand
5
Centre for Sensory-Motor Interactions, Department of Health Science and Technology, Aalborg University, 9220 Aalborg, Denmark
6
Riphah College of Rehabilitation Sciences, Riphah International University, Islamabad 46000, Pakistan
*
Author to whom correspondence should be addressed.
Submission received: 2 April 2020 / Revised: 14 April 2020 / Accepted: 26 April 2020 / Published: 27 April 2020
(This article belongs to the Special Issue Brain Plasticity and Motor Control—Series II)

Abstract

:
Objective: The purpose of this study was to evaluate the impact of chiropractic spinal manipulation on the early somatosensory evoked potentials (SEPs) and resting-state electroencephalography (EEG) recorded from chronic stroke patients. Methods: Seventeen male patients (53 ± 12 years old) participated in this randomized cross-over study. The patients received chiropractic spinal manipulation and control intervention, in random order, separated by at least 24 hours. EEG was recorded before and after each intervention during rest and stimulation of the non-paretic median nerve. For resting-state EEG, the delta-alpha ratio, brain-symmetry index, and power-spectra were calculated. For SEPs, the amplitudes and latencies of N20 and N30 peaks were assessed. Source localization was performed on the power-spectra of resting-state EEG and the N30 SEP peak. Results: Following spinal manipulation, the N30 amplitude increased by 39%, which was a significant increase compared to the control intervention (p < 0.01). The latency and changes to the strength of the cortical sources underlying the N30 peak were not significant. The N20 peak, the resting-state power-spectra, delta-alpha ratio, brain-symmetry index, and resting-state source localization showed no significant changes after either intervention. Conclusion: A single session of chiropractic spinal manipulation increased the amplitude of the N30 SEP peak in a group of chronic stroke patients, which may reflect changes to early sensorimotor function. More research is required to investigate the long-term effects of chiropractic spinal manipulation, to better understand what impact it may have on the neurological function of stroke survivors.

1. Introduction

Stroke is the second-most common cause of death globally, preceded only by ischemic heart disease [1]. It has a high prevalence, affecting approximately 200 persons per 100,000 [2] and requires extensive rehabilitation, with high economic and social costs (~€21,000 per patient in 2010) [3]. Upper limb function is often affected by stroke, and its improvement has been identified as one of the top ten research priorities by stroke survivors, caregivers, and clinicians [4].
The predictors known to identify inefficient upper limb sensorimotor recovery are increased stroke severity, more severe early somatosensory and motor impairments, and the presence of visuospatial neglect [5]. Somatosensory evoked potentials (SEPs) are frequently used as a measure of somatosensory processing after stroke, and abnormalities in the median nerve SEPs often predict functional outcomes for stroke survivors [6,7,8,9,10,11,12]. Compared to healthy controls or the same subjects’ unaffected hemisphere, typical median nerve SEP abnormalities are found in affected hemispheres of stroke survivors. These are abnormally long interpeak intervals (such as N9-N13, N13-N20), absent or severely reduced SEP peak amplitudes (such as N20 and P22 amplitudes), and prolonged latency of some SEP peaks (such as 2-5 standard deviations higher latency of N20) [7,10,11,13,14]. A recent systematic review suggested that acute, inter-hemispheric imbalances involving higher alpha event-related synchronization in the affected hemisphere could be more pronounced in those with moderate upper limb motor impairments [11].
The cortical and subcortical areas associated with the movement are interconnected [15], and the active networks can be seen during the assessment of the electrical resting-state activity of the brain. After a stroke, the reorganization of these electrophysiological networks can result in changes, for example, in the functional connectivity [16], and the power spectrum [17]. The changes in resting activity have been associated with motor dysfunction during movements [18]. Increased activity of the slower oscillations, i.e., increased power of the delta and theta frequency bands have been found to be related to brain damage due to, for example, stroke, brain hemorrhage, tumors and traumatic injury [19,20,21,22]. Finnigan et al. [23] suggested that the delta-alpha ratio (DAR) and pairwise-derived brain symmetry index (BSI) calculated from the resting-state electroencephalography (EEG) can be important predictors of neurological function in stroke survivors. Increased DAR and BSI were found in the (sub-) acute stage after stroke [23,24,25]. However, only BSI in lower frequency bands was increased in the chronic stroke survivors [17].
The possible mechanisms involved in post-stroke motor recovery consist of facilitation and modulation of neural plastic changes in the brain [26]. In the past two decades, many studies have shown that chiropractic spinal manipulation has a neural plastic effect on the central nervous system (CNS) [27,28,29]. Numerous studies have demonstrated altered central processing, such as those associated with somatosensory processing, sensorimotor integration, motor control, and pain after (usually a single session of) chiropractic spinal manipulation, which suggests that the chiropractic intervention can rapidly affect central neural function in a variety of ways [27,28,30,31,32,33,34,35,36,37,38]. Several of these studies have shown that chiropractic care can alter the amplitudes of several SEP peaks, in particular, the N20 and N30 peaks [32,36,39,40]. The most consistent change following chiropractic care is a reduction in amplitude of the N30 SEP peak [32,36,39,40], which has been shown to occur primarily in the prefrontal cortex [39]. The prefrontal cortex is highly involved in sensorimotor integration and motor control [32,34,35,38]. Thus, the changes in the prefrontal cortex might be the reason why a single session of chiropractic care could impact and improve force development, as has recently been shown in a chronic stroke population [41]. However, the effects of chiropractic care on post-stroke sensory and motor recovery have not yet been adequately studied.
Considering that N30 SEP peak amplitudes have been shown to decrease due to reduced cerebral blood flow [42], and previous studies have consistently shown a decrease in N30 SEP peak amplitude after chiropractic care [32,36,39,40] albeit, in a different population, we considered it to be worth investigating what effect chiropractic care has on N30 SEP peaks, as well as other measures of brain activity, in stroke survivors. Therefore, this study aimed to evaluate the effects of a single session of chiropractic spinal manipulation on the resting-state EEG, early somatosensory evoked potentials, and the strength of its corresponding brain sources in stroke survivors.

2. Methods

The study used a randomized controlled cross-over design and was conducted at Railway General Hospital in Rawalpindi, Pakistan. The Riphah International University Research Ethics Committee, Pakistan, approved the study (ref # Riphah/RCRS/REC/000118). The study was also approved by the New Zealand College of Chiropractic Research Review Committee. The study was conducted in accordance with the Declaration of Helsinki.

2.1. Subjects

Nineteen stroke patients (all males, 53 ± 12 years old) participated in this study. By using a purposive sampling technique, the subjects were recruited from the outpatient facility of the rehabilitation department, where they came for their conventional physical rehabilitation. The subjects gave their written informed consent to participate in the study. The patient details are given in Table 1. Although it was not part of the inclusion criteria, all subjects recruited were males and naïve to chiropractic care.
Before enrolling in the study, the subjects were introduced to the lab environment. Subjects were eligible to participate if they had suffered from a stroke at least 12 weeks before their involvement in the study and had a Fugl-Meyer Assessment score (combined upper and lower limb) of less than or equal to 85 (i.e., they had significant motor impairment [43]). The subjects were ineligible to participate if they showed no evidence of spinal dysfunction (i.e., presence of vertebral subluxation indicators identified by a chiropractor), had absolute contraindications to spinal adjustments (including spinal fracture, atlantoaxial instability, spinal infection, spinal tumor, or cauda equina syndrome), or previously had a significant adverse response to chiropractic care.

2.2. Experimental Protocol

The subjects participated in two sessions, chiropractic and control, in random order, separated by at least 24 hours. The balanced randomization scheme was generated using Minimizer (Microsoft Corp. Redmond, WA, USA). The subjects were not informed beforehand that one of the interventions would be a control intervention.
Each session consisted of recording resting-state EEG, followed by SEPs evoked by electrical stimulation of the median nerve of the non-paretic limb. Both recordings were done before and after the intervention. During each session, the subjects were seated comfortably in a chair, in front of a screen and were asked to keep their eyes open and be relaxed to reduce the contamination of EEG signals.
The study was single-blinded; therefore, the subjects did not know which intervention they received. The data analysis was performed by blinding the analyst through assigning random numbers to each dataset (recording) of every subject.

2.3. Interventions

The chiropractic spinal manipulation and control interventions were similar to those used in previous studies [27,28,37,39,41] that have investigated the neurophysiological effects of chiropractic spinal adjustments. The same chiropractor performed the actual and control adjustments. At the end of the second session, the subjects were asked if they felt that they had undergone active treatment in each session (‘yes’ or ‘no’).

2.3.1. Chiropractic Manipulation

The standard adjustment techniques used by the chiropractors, also known as spinal manipulation, were used in the chiropractic spinal adjustment session. The chiropractor performed manual high-velocity low-amplitude adjustments or instrument-assisted adjustments to the spine or pelvic joints [44]. The chiropractor used standard clinical indicators of spinal and pelvic dysfunction to decide where to adjust. [45] These indicators included tenderness to palpation, restricted intersegmental motion, muscle asymmetry, and blocked jopint play or end-feel. Chiropractic adjustments were applied to multiple spinal segments if required.

2.3.2. Control Manipulation

The control intervention was performed by the same chiropractor who provided the chiropractic intervention. In the control session, the chiropractor interacted with the patient in a similar way to the active session, including assessing the spine and pelvis for dysfunction and then moving and setting up the patient as if they were going to apply an adjustive thrust. However, during the adjustment set up, the chiropractor took care not to provide an adjustive thrust or to take a vertebral segment that was deemed to be subluxated to tension. The control session was designed to control for the interaction and time taken during the chiropractic intervention and to control for the mechanoreceptive input associated with the chiropractor assessing the patient’s spine, while ensuring the afferent input associated with the adjustive thrust was minimized.

2.4. Median Nerve Stimulation

The median nerve was stimulated using electrical pulses delivered by the electrical stimulator (Digitimer DS7AH, Hertfordshire, UK) to evoke SEPs. The stimulation electrodes (Neuroline 700, AMBU A/S, Ballerup, Denmark) were placed at the wrist. The motor threshold was defined as the lowest current intensity, which elicited a visible twitch of the thumb. Before and after each intervention, a total of 1000 electrical pulses were given to the median nerve of the non-paretic limb. The stimulation pulse was monophasic, with a width of 0.2 ms and a frequency of 2.3 Hz.

2.5. EEG

The EEG was recorded at a sampling rate of 2048 Hz from 62 channels using a REFA amplifier (TMSi, Twente, The Netherlands) according to the 10-20 electrode system. The ground electrode was placed at AFz. The impedance of the electrodes was kept below 10 kΩ. The subjects were asked to focus on a white fixation cross with black background displayed in the center of a computer screen while minimizing eye blinks, eye movements, and facial movements.
The EEG was analyzed offline using EEGLAB version 14.1.1 [46], ERPLAB version 6.1.4 [47], and FieldTrip version 20180912 [48] running on MATLAB 2015b (The MathWorks, Inc., Natick, MA, USA). We developed custom scripts in MATLAB utilizing EEGLAB, ERPLAB, FieldTrip, and MATLAB functions to perform the analysis. The analysis pipeline has been previously described in detail [49].
The raw EEG was truncated to keep up to an additional 30 s of data in the beginning and end of the recordings to reduce filter artifacts. The PREP pipeline [50] was used to identify bad channels, remove line noise, and obtain the average referenced data.

2.5.1. Resting-State EEG

The resting-state EEG was recorded for 5 minutes. For artifact detection and correction, the PREPed EEG was high-pass filtered using a Kaiser windowed FIR filter (β = 5.653) with a cutoff frequency of 1 Hz and an order of 4948 equivalent to the transition bandwidth of 1.5 Hz. After cleaning the data, the channels of the subjects with the right affected hemisphere were mirrored in the sagittal plane to model all subjects having the lesion(s) in the left hemisphere. The data from subjects with at least 2 min of clean EEG per session were used for further analysis. The analyses were performed on the data downsampled to 256 Hz to reduce the computational time and remove redundant frequency spectrum not required in the analysis. The EEG was segmented into 2 s long epochs to encompass two cycles of the lowest frequency of interest (1 Hz). For any session with clean data longer than 2 min, 60 epochs were randomly selected to equalize the length of data across sessions and subjects.

Spatio-Spectral Power

The power spectra were calculated between 1 and 80 Hz using the Fourier basis with a 2 s wide Hanning window. Afterwards, the average power of each classical frequency band (delta (1–4 Hz), theta (4.1–8 Hz), alpha (8.1–12 Hz), beta (12.1–32 Hz), and gamma (32.1–80 Hz)) was computed.

DAR

For each channel c , the DAR was calculated to get the ratio of the mean delta power to the mean alpha power as
D A R c = P c ( f ) f = 1 ,   ,   4   H z P c ( f ) f = 8 ,   ,     12   H z
The global DAR was calculated by averaging the ratios across all N EEG channels:
D A R = 1 N c = 1 N D A R c

BSI

The pairwise-derived BSI was calculated by taking the absolute difference of the mean spectral power between the homologous left c L and right c R EEG channels. The BSI was calculated over the 1–25 Hz range [17] by
B S I c p = | P c R ( f ) P c L ( f ) P c R ( f ) + P c L ( f ) | f = 1 , , 25   H z
The BSI across all channel pairs c p was averaged to get global BSI:
B S I = 2 N c p = 1 N / 2 B S I c p
The range of BSI is 0 to 1, where 0 means perfect symmetry, and 1 represents no symmetry between channel pairs.
Since BSI does not take the direction of asymmetry into account, the direction of asymmetry was calculated by not taking the absolute of Equation (3):
B S I d i r c p = P c R ( f ) P c L ( f ) P c R ( f ) + P c L ( f ) f = 1 , , 25   H z
The values were averaged over all channel pairs:
B S I d i r = 2 N c p = 1 N / 2 B S I d i r c p
B S I d i r ranges from −1 to +1, where 0 means perfect symmetry, positive values mean higher power in the right hemisphere compared to the left hemisphere and vice versa for the negative values. B S I d i r for all subjects was multiplied by −1. This way, the positive values characterized higher power in the affected hemisphere compared to the unaffected hemisphere, and vice versa for the negative values [17].
In addition to BSI for the 1–25 Hz band, separate BSI calculations were made for the delta, theta, alpha, beta, and gamma bands. The mid-line channels were not included in the BSI calculation.

Source Localization

The localization of electrical activity in the brain during rest was estimated in the LORETA-KEY software, version 20151222 [51] (freely available at www.uzh.ch/keyinst/loreta). For source localization, we used sLORETA, which is a linear inverse algorithm that estimates the distribution of cortical generators of the EEG in three-dimensions, with lowest localization error compared to several other linear inverse methods [51]. The sLORETA implementation uses a reference brain from the Montreal Neurological Institute (average MRI brain map (MNI-152)) [52] with cortical gray matter divided into 6239 voxels with a resolution of 5 mm.
For sLORETA analysis, the EEG was segmented into 8 s long epochs to obtain smooth power spectral density. The data from subjects with at least 2 minutes of clean EEG per session were used. For any session with clean data longer than 2 minutes, 15 epochs were randomly selected to equalize the length of data across sessions and subjects. The sLORETA was performed in the frequency domain where cross-spectral matrices for each subject were computed in the LORETA-KEY software for the same five frequency bands as in the power spectral analysis above. Subsequently, the cross-spectral matrices were averaged for each subject and used as input for the sLORETA.

2.5.2. SEPs

The EEG recorded during the median nerve stimulation of the non-paretic limb was analyzed. For artifact detection and correction, the PREPed EEG was high-pass filtered using a 2nd order Butterworth filter with a cutoff frequency of 1 Hz. The epochs were extracted from 100 ms before the stimulus to 150 ms after the stimulus and baseline corrected using the pre-stimulus period. After cleaning the data, the channels of the subjects with the right affected hemisphere were mirrored in the sagittal plane to model all subjects having a lesion(s) in the left hemisphere. For each subject, the number of epochs in each session was equalized to the session with a minimum number of epochs (for that subject), by randomly removing clean excess epochs. Finally, the epochs were averaged.

SEP Peaks

In this study, we analyzed the amplitudes and latencies of the N20 and N30 SEP peaks. The N20 amplitude was calculated from the central electrode. The most positive (P14) and the most negative (N20) peaks with respect to the stimulus were identified in the time window of 13–20 ms and 20–30 ms, respectively. The N30 amplitude was calculated from the frontal electrode, as done previously [39]. The most positive (P22) and the most negative (N30) peaks with respect to the stimulus were identified in the time window of 15–25 ms and 25–45 ms, respectively. Afterwards, manual inspection was performed, and the identified peaks were verified by an expert in SEP analysis. The N20 and N30 amplitudes were defined as the absolute difference in the amplitudes of these peaks.

Source Localization

Brain source modeling was performed in the 18 to 60 ms post-stimulus period using Brain Electrical Source Analysis (BESA) (BESA Research 5.3; MEGIS Software GmbH, Gräfelfing, Germany) software. First, the potential distributions over the scalp from preset voltage dipoles within the brain were calculated. Afterwards, the agreement between the recorded and the calculated field potentials was evaluated. For further analysis, the percentage of data that could not be explained by the model, termed as residual variance (RV), was required to be less than 10%. A 4-shell ellipsoidal model with a radius of 85 mm was used.
The models were first created on the pre-session grand-averages to get an indication of the location and number of sources. Afterwards, LORETA was used on individual pre-sessions to guide the estimation of the location and number of sources. LORETA is a current density model that produces blurred source images, requires no prior constraints, and has high accuracy [53]. Once the dipoles were placed, the model fit was obtained by fixing their locations but allowing their orientations to move freely. The model was then applied to the associated post-sessions. The source activation waveforms were exported to MATLAB 2015b (The MathWorks, Inc., Natick, MA, USA), and brain source strengths for the N30 peak were computed by calculating the area under the curve (AUC) in the 20 to 60 ms post-stimulus period.

2.5.3. Artifact Removal

The same parameters were used for marking epochs as bad for both the resting-state EEG and SEPs. The resting-state EEG was segmented into 0.5 s long epochs. An epoch was considered bad if, for any channel, (i) the amplitude exceeded 100 μV, (ii) peak-to-peak amplitude was more than 150 μV in any sliding window of 200 ms width with a step size of 100 ms, (iii) the amplitude was greater than 100 μV in a step-function with a sliding window 200 ms wide with a step size of 50 ms, (iv) sample-to-sample difference was more than 50 μV, or (v) the amplitude was less than 2 μV for 125 ms (SEPs) or 150 ms (resting-state EEG) (i.e., flat-lined data). For SEPs, the above conditions were not applied to −2 ms to 2 ms period corresponding to the stimulus artifact. Afterwards, all epochs were manually verified to classify into good and bad epochs. The epochs with step-like artifacts in the frontal channels were not removed as they were related to eye-blinks and eye-movements [54].
The high-pass filtered EEG was downsampled to 512 Hz and epoched in a similar way as described above. The bad epochs and bad channels (from PREP) were removed from the data, and adaptive mixture independent component analysis was employed to decompose EEG into maximally independent components i.e., which are spatially fixed and temporally discrete [55]. This algorithm was used since its performance is better than that of other ICA algorithms [56]. The ICA weights obtained were applied to band-pass (0.5–1000 Hz) filtered PREPed data with bad epochs removed. For resting-state EEG, a Kaiser windowed FIR filter (β = 5.653) with an order of 7420 corresponding to the transition bandwidth of 1 Hz was used. For SEPs, a 2nd order Butterworth filter was employed. Afterwards, all the independent components (ICs) were manually categorized into the brain or non-brain components resulting from activities of muscle or eye, or noise of channel or mains. The ICs were categorized based on their spatial distribution (scalp topography), time course, spectrograms, event-related potential (ERP) images, and current dipole models using recommendations from [57,58] and the website https://labeling.ucsd.edu/. The bad ICs were removed, followed by interpolation of the noisy channels to give cleaned datasets.

2.6. Statistics

The data are reported as mean ± SD unless otherwise indicated. The statistical significance threshold was set at p < 0.05.
Two-way, repeated-measures ANOVAs were used with factors of intervention (Control and Chiropractic) and session (pre and post) to identify the effects of the intervention on the SEPs’ amplitudes and latencies, brain source strengths, DAR and BSI. Tukey’s HSD was used to perform pairwise comparisons. The statistical procedure was performed in MATLAB 2015b (The MathWorks, Inc., Natick, MA, USA).
Non-parametric cluster-based permutation tests [59] were used with a cluster-threshold of 0.05 to evaluate the differences between the interventions based on the global power spectrum of the resting-state. The clusters were defined as two or more channel-frequency pairs, each with p < 0.05 from the paired t-test (two-tailed). The maximum of cluster-level statistics, obtained by adding the t-values within each cluster, 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, which was approximated by the Monte Carlo method with 5000 permutations.
The statistical procedure for source localization of resting-state EEG was performed in the LORETA-KEY software using non-parametric mapping [60], which utilized Fisher’s random permutation test with 5000 randomizations to control for multiple comparison problem. The paired two-tailed t-test was used to find differences in the current sources across the different frequency bands. The tests were used to identify differences between the pre-sessions; the post- and pre-chiropractic sessions; and the post- and pre-control sessions.

3. Results

Out of nineteen chronic stroke survivors recruited, two were excluded as they were unable to complete the experiment—one was uncomfortable with the median nerve stimulation, and the other could not follow the instructions of the experimenter. Hence, the analyses were performed on the remaining 17 subjects (53 ± 12 years old). From the questions to evaluate the success of subject blinding, out of the 17 subjects, only two felt that one of the sessions was not an active session, and one of these was correct in the identification of the order of the interventions (chiropractic or control) he received.

3.1. Resting-State EEG

The analyses were performed on 16 subjects who had a minimum of 2 minutes of clean resting-state EEG in every session.

3.1.1. Spatio-Spectral Power

For all frequency bands, no differences in the power were seen between the pre-intervention sessions (no clusters) (Figure 1A), post- and pre-control intervention (1 positive cluster p = 0.055) (Figure 1B), and post- and pre-chiropractic spinal interventions (no clusters) (Figure 1C). There was a trend towards higher grand-averaged absolute power in all frequency bands after the chiropractic intervention compared to post-control intervention, which was similar or lower than the pre-control intervention. The difference between the two is shown in Figure 1D.

3.1.2. DAR

The ANOVA revealed that there was no interaction between the interventions and sessions over DAR values (all p > 0.05). Descriptive statistics are given in Table 2.

3.1.3. BSI

No interaction between interventions and sessions was found in any of the frequency band’s BSI values and directions (all p > 0.05). Descriptive statistics are given in Table 2.

3.1.4. Source Localization

For all frequency bands, the sLORETA analysis showed that the pre-sessions were similar. A slight increase in activity was seen after the control intervention, however, neither of the interventions had any significant effect on brain activity (all p > 0.05).

3.2. SEPs

The number of epochs per subject used for the analysis were 799 ± 124 ([min, max] = [437, 956]).

3.2.1. SEP Peaks

The N20 amplitude was not affected by either intervention (all p > 0.05) (Figure 2A). The ANOVA (Table 3) showed a significant interaction of the intervention and session for the N30 SEP peak amplitude (p = 0.03). The post hoc analysis showed that the N30 amplitude was increased by 39% after the chiropractic spinal manipulation (p = 0.007; 95% CI = [0.11, 0.78]). The N30 amplitude was similar for the pre-intervention sessions (control vs chiropractic: p = 0.12; 95% CI = [−0.09, 0.99]) and there was no change after the control intervention (post vs pre: p = 0.95; 95% CI = [−0.56, 0.39]). Figure 2B shows the distribution of N30 amplitude across the four sessions. None of the SEPs latencies were significantly changed following either intervention (all p > 0.05) (Figure 3).

3.2.2. Source Localization

The analyses were performed on 16 subjects because for one subject one of their sessions had an RV that was greater than 10%. The LORETA solution revealed five distinct areas during the 20–60 ms post-stimulus period; contralateral primary somatosensory cortex (SI), prefrontal cortex, cingulate, and bilateral secondary somatosensory cortices (SII). Therefore, a 5-source solution was assumed, and the dipoles were placed in these areas.
The ANOVA showed no interactions between the interventions and sessions for source strengths (all p > 0.05). Table 4 shows the coordinates of the brain sources and AUCs of their activity waveforms. The AUCs obtained are shown in Figure 4.

4. Discussion

In this randomized cross-over study, we evaluated the effects of a single session of chiropractic spinal manipulation versus a control intervention on stroke patients by investigating the change in the resting-state EEG, early somatosensory evoked potentials and their underlying brain sources. We found that chiropractic spinal manipulation increased the N30 amplitude of the somatosensory potential evoked by median nerve stimulation, but there were no changes in the N20 amplitude and latencies. There were no effects on the strength of the N30 peak brain sources by either intervention. There were no changes in the power of frequency bands, resting-state EEG source localization, delta-alpha ratio, and brain-symmetry index after either intervention. To the best of our knowledge, the current randomized study is the first of its kind to investigate these parameters in a stroke population following chiropractic spinal manipulation.

4.1. Resting-State

The current study showed a trend towards higher grand-averaged absolute power in all frequency bands after the chiropractic intervention compared to post-control intervention, which was similar or lower than the pre-control intervention. However, none of these increases were significant. In our previous study, performed on sub-clinical pain patients [37], we found no differences in the spectral power and in the sLORETA based source localization after the chiropractic spinal manipulation and control interventions. Similarly, in this study, neither of the interventions had a significant effect on the resting-state EEG, although there was a trend towards higher grand-average absolute power in all frequency bands. The probable reason for this is likely the non-uniformity in the type and location of the stroke in the subjects. This variation among subjects makes finding a statistical difference more difficult and may have required a larger sample size. Hence, there is a need for more research to explore longer chiropractic intervention periods (>1 session) and larger sample sizes or more homogenous subject groups in future studies of this kind. It is possible that this study was underpowered to show any potential changes following the chiropractic intervention. It is also possible, of course, that chiropractic spinal manipulation does not affect the resting-state spectral power.
Since DAR has been shown to be a predictive parameter of recovery in (sub-) acute stroke [21,62], it was used in this study to see if it was affected by the chiropractic spinal manipulation. It was found that DAR was not changed significantly with either the chiropractic spinal manipulation or control intervention. The DAR has been found to be increased in (sub-) acute stroke survivors compared to healthy persons [23,24]. However, it has been shown to be similar between chronic stroke survivors and healthy individuals [17]. Therefore, in this study of chronic stroke survivors, it is plausible that this parameter is not predictive of recovery in stroke, and therefore, it was not affected by the interventions.
Previously, in a study conducted in chronic stroke patients [17], BSI was found to be increased in the delta and theta bands but was unchanged in the alpha and beta bands. The same study also showed, using directional BSI, that the affected hemisphere had more resting-state power compared to the unaffected hemisphere, particularly in the delta and theta bands. In the present study, however, the BSI was not affected by the interventions. The possible reasons for this were inter-individual differences, the cross-over design of the study with less time between interventions, or more sessions of the treatment is required to affect the BSI. Since BSI is based on spectral power, it is also possible that chiropractic spinal manipulation does not affect the spectral power, and hence, does not affect BSI.

4.2. SEPs

In the current study, a 39% increase in the N30 SEP peak amplitude was observed following the chiropractic intervention. Previously, SEP studies looking at the effects of chiropractic spinal manipulation, compared to a control intervention, have shown that the N30 SEP amplitude decreased by 17–30% [32,35,36,39]. Lelic et al. [39] showed that this decrease in N30 SEP peak amplitude occurred in the prefrontal cortex, where there was a decrease of approximately 20% in the activity. Most of these studies were conducted in subclinical pain patients.
The N30 SEP peak is believed to be generated at the motor, premotor, and prefrontal cortices [63,64,65], and reflects sensorimotor integration [66]. Effective sensorimotor integration is required for learning new motor skills as well as for recovering from injuries such as stroke [67,68]. The prefrontal cortex, in particular, is of interest since it is associated with cognitive thinking and decision-making using the somatosensory input and information from other internal and external sources [69]. Early SEP peak amplitudes are known to be severely diminished or even absent in stroke populations [7,10,11,13,14]. It is, therefore, interesting that the chiropractic intervention showed an increase in the N30 SEP peak amplitude, considering previous studies in subclinical subjects have consistently reported a decrease in N30 SEP peak amplitude after chiropractic spinal manipulation [32,36,39,40]. The increase observed in the present study may reflect an improvement in early sensorimotor function that is related to the observed increase in strength found in a recent study of chronic stroke survivors after a single session of chiropractic spinal manipulation [41]. In that study [41], Holt et al. investigated the effect of chiropractic spinal manipulation on motor recovery in a chronic stroke population and found that, on average, plantarflexion muscle strength increased by approximately 64% after a single session of chiropractic care.
The source strength of the N30 SEP peak was not affected by the interventions. The possible reason for no significant changes in the source strength of the N30 SEP peak can be that the BESA analysis was performed on a head model of the healthy brain which did not cater for the number, location, and size of the lesion(s) of the participants in this current study. Using individual MRI to obtain a realistic head model of each subject could have improved the accuracy of source modeling, as suggested by previous simulation studies [70,71,72]. Therefore, this limitation should be kept in mind when interpreting the results.

5. Study Considerations

There are several reasons that no changes in the brain activity were seen in other parameters apart from the N30 amplitude after the chiropractic spinal manipulation intervention. These include that chiropractic spinal manipulation does not alter these parameters, or it may be because there are differences in the brain morphology of the stroke patients and non-uniformity of the type of stroke and affected brain regions. The enrollment criteria of future studies could be modified to include more homogenous patients with respect to stroke type and location. This would, of course, make the recruitment of patients more difficult, and the results of that study may not be generalizable to the majority of the stroke population.
This study was an exploratory study with a sample size of 19. Previous studies based on stroke patients had sample sizes similar to this study (10 to 21 participants) [17,20,25,62,73]. However, we recognize that the sample size used may not be large enough, and this study may, therefore, be underpowered to detect changes following spinal manipulation that occur in this population. Future studies that further explore the potential changes following chiropractic care for this population should consider increasing the sample size.
It is also important to keep in mind that previous work has shown that some neuroplastic effects of chiropractic spinal manipulation can, in some instances, remain for at least a week [74]. This may have affected the current study because it could have reduced or nullified the changes in this study, where both interventions were performed within a week. Future studies should control for this factor by having more than a week between the interventions, or they should use a parallel-group design.

6. Conclusions

The findings of the study suggest that in chronic stroke survivors, a single session of chiropractic spinal manipulation increases the cortical activity related to the somatosensory evoked N30 potential, which may alter the early sensorimotor functionality. Future studies should investigate the long-term effects of chiropractic care in this population group to understand better the impact of chiropractic care on cortical activity in stroke survivors.

Author Contributions

Conceptualization, M.S.N., I.K.N., D.L., K.H., A.M.D. and H.H.; methodology, M.S.N.; software, M.S.N.; validation, M.S.N., I.K.N. and K.H.; formal analysis, M.S.N.; investigation, M.S.N, I.K.N., R.B.N., K.H. and I.A.; resources, I.K.N., I.A., A.M.D. and H.H.; data curation, M.S.N., I.K.N, K.H., I.A.; writing—original draft preparation, M.S.N.; writing—review and editing, M.S.N., I.K.N., D.L., R.B.N., K.H., I.A., A.M.D. and H.H.; visualization, M.S.N..; supervision, I.K.N., D.L., A.M.D. and H.H.; project administration, M.S.N., I.K.N., A.M.D. and H.H.; funding acquisition, A.M.D. and H.H. All authors have read and agreed to the published version of the manuscript.

Funding

The grant proposal for this study was reviewed by the Australian Spinal Research Foundation to support the facilitation of funding from the United Chiropractic Association. This study also received co-funding from Hamblin Chiropractic Research Fund Trust, The Rubicon Group, the Scottish Chiropractic Association, Halsa Chiropractic, and Physiotherapy and through donations to the Centre for Chiropractic Research Supporters Program at the New Zealand College of Chiropractic.

Acknowledgments

We would like to thank all participants who took part in the study.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. Benjamin, E.J.; Virani, S.S.; Callaway, C.W.; Chamberlain, A.M.; Chang, A.R.; Cheng, S.; Chiuve, S.E.; Cushman, M.; Delling, F.N.; Deo, R.; et al. Heart disease and stroke statistics - 2018 update: A report from the American Heart Association. Circulation 2018, 137, E67–E492. [Google Scholar] [CrossRef] [PubMed]
  2. Thrift, A.G.; Thayabaranathan, T.; Howard, G.; Howard, V.J.; Rothwell, P.M.; Feigin, V.L.; Norrving, B.; Donnan, G.A.; Cadilhac, D.A. Global stroke statistics. Int. J. Stroke 2017, 12, 13–32. [Google Scholar] [CrossRef] [PubMed]
  3. Gustavsson, A.; Svensson, M.; Jacobi, F.; Allgulander, C.; Alonso, J.; Beghi, E.; Dodel, R.; Ekman, M.; Faravelli, C.; Fratiglioni, L.; et al. Cost of disorders of the brain in Europe 2010. Eur. Neuropsychopharmacol. 2011, 21, 718–779. [Google Scholar] [CrossRef] [Green Version]
  4. Pollock, A.; George, B.S.; Fenton, M.; Firkins, L. Top 10 research priorities relating to life after stroke - consensus from stroke survivors, caregivers, and health professionals. Int. J. Stroke 2014, 9, 313–320. [Google Scholar] [CrossRef] [PubMed]
  5. Meyer, S.; De Bruyn, N.; Lafosse, C.; Van Dijk, M.; Michielsen, M.; Thijs, L.; Truyens, V.; Oostra, K.; Krumlinde-Sundholm, L.; Peeters, A.; et al. Somatosensory Impairments in the Upper Limb Poststroke: Distribution and Association With Motor Function and Visuospatial Neglect. Neurorehabilit. Neural Repair 2016, 30, 731–742. [Google Scholar] [CrossRef] [PubMed]
  6. Haupt, W.F.; Chopan, G.; Sobesky, J.; Liu, W.C.; Dohmen, C. Prognostic value of somatosensory evoked potentials, neuron-specific enolase, and S100 for short-term outcome in ischemic stroke. J. Neurophysiol. 2016, 115, 1273–1278. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Zeman, B.D.; Yiannikas, C. Functional prognosis in stroke: Use of somatosensory evoked potentials. J. Neurol. Neurosurg. Psychiatry 1989, 52, 242–247. [Google Scholar] [CrossRef] [Green Version]
  8. Keren, O.; Ring, H.; Solzi, P.; Pratt, H.; Groswasser, Z. Upper limb somatosensory evoked potentials as a predictor of rehabilitation progress in dominant hemisphere stroke patients. Stroke 1993, 24, 1789–1793. [Google Scholar] [CrossRef] [Green Version]
  9. Lee, S.Y.; Lim, J.Y.; Kang, E.K.; Han, M.K.; Bae, H.J.; Paik, N.J. Prediction of good functional recovery after stroke based on combined motor and somatosensory evoked potential findings. J. Rehabilitation Med. 2010, 42, 16–20. [Google Scholar] [CrossRef] [Green Version]
  10. Al-Rawi, M.A.W.; Hamdan, F.B.; Abdul-Muttalib, A.K. Somatosensory evoked potentials as a predictor for functional recovery of the upper limb in patients with stroke. J. Stroke Cerebrovasc. Dis. Off. J. Natl. Stroke Assoc. 2009, 18, 262–268. [Google Scholar] [CrossRef]
  11. Triccas, L.T.; Meyer, S.; Mantini, D.; Camilleri, K.; Falzon, O.; Camilleri, T.; Verheyden, G. A systematic review investigating the relationship of electroencephalography and magnetoencephalography measurements with sensorimotor upper limb impairments after stroke. J. Neurosci. Methods 2019, 311, 318–330. [Google Scholar] [CrossRef] [PubMed]
  12. Yoon, H.S.; Cha, Y.J.; Sohn, M.K.; You, J.H. Effect of rehabilitation on the somatosensory evoked potentials and gait performance of hemiparetic stroke patients. Technol. Health Care 2018, 26, S145–S150. [Google Scholar] [CrossRef] [Green Version]
  13. Gott, P.S.; Karnaze, D.S.; Fisher, M. Assessment of median nerve somatosensory evoked potentials in cerebral ischemia. Stroke 1990, 21, 1167–1171. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Feys, H.; Van Hees, J.; Bruyninckx, F.; Mercelis, R.; De Weerdt, W. Value of somatosensory and motor evoked potentials in predicting arm recovery after a stroke. J. Neurol. Neurosurg. Psychiatry 2000, 68, 323–331. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Guillery, R.W.; Sherman, S.M. The thalamus as a monitor of motor outputs. In Proceedings of the Philosophical Transactions of the Royal Society B: Biological Sciences; Royal Society: London, UK, 2002; Volume 357, pp. 1809–1821. [Google Scholar]
  16. Urbin, M.A.; Hong, X.; Lang, C.E.; Carter, A.R. Resting-state functional connectivity and its association with multiple domains of upper-extremity function in chronic stroke. Neurorehabilit. Neural Repair 2014, 28, 761–769. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Saes, M.; Meskers, C.G.M.; Daffertshofer, A.; de Munck, J.C.; Kwakkel, G.; van Wegen, E.E.H. How does upper extremity Fugl-Meyer motor score relate to resting-state EEG in chronic stroke? A power spectral density analysis. Clin. Neurophysiol. 2019, 130, 856–862. [Google Scholar] [CrossRef]
  18. Carter, A.R.; Shulman, G.L.; Corbetta, M. Why use a connectivity-based approach to study stroke and recovery of function? NeuroImage 2012, 62, 2271–2280. [Google Scholar] [CrossRef] [Green Version]
  19. Andraus, M.E.C.; Alves-Leon, S.V. Non-epileptiform EEG abnormalities: An overview. Arq. de Neuro-Psiquiatr. 2011, 69, 829–835. [Google Scholar] [CrossRef] [Green Version]
  20. Van Putten, M.J.A.M.; Tavy, D.L.J. Continuous quantitative EEG monitoring in hemispheric stroke patients using the brain symmetry index. Stroke 2004, 35, 2489–2492. [Google Scholar] [CrossRef] [Green Version]
  21. Finnigan, S.; van Putten, M.J.A.M. EEG in ischaemic stroke: Quantitative EEG can uniquely inform (sub-)acute prognoses and clinical management. Clin. Neurophysiol. Off. J. Int. Fed. Clin. Neurophysiol. 2013, 124, 10–19. [Google Scholar] [CrossRef]
  22. Britton, J.W.; Frey, L.C.; Hopp, J.L.; Korb, P.; Koubeissi, M.Z.; Lievens, W.E.; Pestana-Knight, E.M.; Louis, E.K.S. Electroencephalography (EEG): An Introductory Text and Atlas of Normal and Abnormal Findings in Adults, Children, and Infants; Allen Press, Inc.: Lawrence, KS, USA, 2016. [Google Scholar]
  23. Finnigan, S.; Wong, A.; Read, S. Defining abnormal slow EEG activity in acute ischaemic stroke: Delta/alpha ratio as an optimal QEEG index. Clin. Neurophysiol. Off. J. Int. Fed. Clin. Neurophysiol. 2016, 127, 1452–1459. [Google Scholar] [CrossRef] [Green Version]
  24. Sheorajpanday, R.V.A.; Nagels, G.; Weeren, A.J.T.M.; van Putten, M.J.A.M.; De Deyn, P.P. Reproducibility and clinical relevance of quantitative EEG parameters in cerebral ischemia: A basic approach. Clin. Neurophysiol. 2009, 120, 845–855. [Google Scholar] [CrossRef] [PubMed]
  25. Anastasi, A.A.; Falzon, O.; Camilleri, K.; Vella, M.; Muscat, R. Brain symmetry index in healthy and stroke patients for assessment and prognosis. Stroke Res. Treat. 2017, 2017, 8276136. [Google Scholar] [CrossRef] [PubMed]
  26. Li, S. Spasticity, Motor Recovery, and Neural Plasticity after Stroke. Front. Neurol. 2017, 8, 120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Haavik, H.; Murphy, B. The role of spinal manipulation in addressing disordered sensorimotor integration and altered motor control. J. Electromyogr. Kinesiol. 2012, 22, 768–776. [Google Scholar] [CrossRef] [PubMed]
  28. Niazi, I.K.; Türker, K.S.; Flavel, S.; Kinget, M.; Duehr, J.; Haavik, H. Changes in H-reflex and V-waves following spinal manipulation. Exp. Brain Res. 2015, 233, 1165–1173. [Google Scholar] [CrossRef]
  29. Pickar, J.G.; Bolton, P.S. Spinal manipulative therapy and somatosensory activation. J. Electromyogr. Kinesiol. 2012, 22, 785–794. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Haavik, H.; Murphy, B. Subclinical neck pain and the effects of cervical manipulation on elbow joint position sense. J. Manip. Physiol. Ther. 2011, 34, 88–97. [Google Scholar] [CrossRef]
  31. Taylor, H.H.; Murphy, B.A. Altered cortical integration of dual somatosensory input following the cessation of a 20 min period of repetitive muscle activity. Exp. Brain Res. 2007, 178, 488–498. [Google Scholar] [CrossRef]
  32. Haavik-Taylor, H.; Murphy, B. Cervical spine manipulation alters sensorimotor integration: A somatosensory evoked potential study. Clin. Neurophysiol. 2007, 118, 391–402. [Google Scholar] [CrossRef]
  33. Haavik, H.; Murphy, B.A. Transient modulation of intracortical inhibition following spinal manipulation. Chiropr. J. Aust. 2007, 37, 106–116. [Google Scholar]
  34. Taylor, H.H.; Murphy, B. Altered Sensorimotor Integration With Cervical Spine Manipulation. J. Manip. Physiol. Ther. 2008, 31, 115–126. [Google Scholar] [CrossRef] [PubMed]
  35. Taylor, H.H.; Murphy, B. Altered Central Integration of Dual Somatosensory Input After Cervical Spine Manipulation. J. Manip. Physiol. Ther. 2010, 33, 178–188. [Google Scholar] [CrossRef] [PubMed]
  36. Taylor, H.H.; Murphy, B. The Effects of Spinal Manipulation on Central Integration of Dual Somatosensory Input Observed After Motor Training: A Crossover Study. J. Manip. Physiol. Ther. 2010, 33, 261–272. [Google Scholar] [CrossRef]
  37. Navid, M.S.; Lelic, D.; Niazi, I.K.; Holt, K.; Mark, E.B.; Drewes, A.M.; Haavik, H. The effects of chiropractic spinal manipulation on central processing of tonic pain - a pilot study using standardized low-resolution brain electromagnetic tomography (sLORETA). Sci. Rep. 2019, 9, 6925. [Google Scholar] [CrossRef] [Green Version]
  38. Marshall, P.; Murphy, B. The effect of sacroiliac joint manipulation on feed-forward activation times of the deep abdominal musculature. J. Manip. Physiol. Ther. 2006, 29, 196–202. [Google Scholar] [CrossRef]
  39. Lelic, D.; Niazi, I.K.; Holt, K.; Jochumsen, M.; Dremstrup, K.; Yielder, P.; Murphy, B.; Drewes, A.M.; Haavik, H. Manipulation of Dysfunctional Spinal Joints Affects Sensorimotor Integration in the Prefrontal Cortex: A Brain Source Localization Study. Neural Plast. 2016, 2016, 3704964. [Google Scholar] [CrossRef] [Green Version]
  40. Haavik, H.; Niazi, I.K.; Jochumsen, M.; Sherwin, D.; Flavel, S.; Türker, K.S. Impact of spinal manipulation on cortical drive to upper and lower limb muscles. Brain Sci. 2017, 7, 2. [Google Scholar] [CrossRef] [Green Version]
  41. Holt, K.; Niazi, I.K.; Nedergaard, R.W.; Duehr, J.; Amjad, I.; Shafique, M.; Anwar, M.N.; Ndetan, H.; Turker, K.S.; Haavik, H. The effects of a single session of chiropractic care on strength, cortical drive, and spinal excitability in stroke patients. Sci. Rep. 2019, 9, 2673. [Google Scholar] [CrossRef]
  42. Guérit, J.M.; Witdoeckt, C.; de Tourtchaninoff, M.; Ghariani, S.; Matta, A.; Dion, R.; Verhelst, R. Somatosensory evoked potential monitoring in carotid surgery. I. Relationships between qualitative SEP alterations and intraoperative events. Electroencephalogr. Clin. Neurophysiol. 1997, 104, 459–469. [Google Scholar] [CrossRef]
  43. Fugl-Meyer, A.R. Post-stroke hemiplegia: Assessment of physical properties. Scand. J. Rehabilitation Med. 1980, 12, 85–93. [Google Scholar]
  44. Cooperstein, R.; Gleberzon, B.J. Technique Systems in Chiropractic; Churchill Livingstone: London, UK, 2004; ISBN 0443074135. [Google Scholar]
  45. Triano, J.J.; Budgell, B.; Bagnulo, A.; Roffey, B.; Bergmann, T.; Cooperstein, R.; Gleberzon, B.; Good, C.; Perron, J.; Tepe, R. Review of methods used by chiropractors to determine the site for applying manipulation. Chiropr. Man. Ther. 2013, 21, 36. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Delorme, A.; Makeig, S. EEGLAB: An open source toolbox for analysis of single-trial EEG dynamics including independent component analysis. J. Neurosci. Methods 2004, 134, 9–21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Lopez-Calderon, J.; Luck, S.J. ERPLAB: An open-source toolbox for the analysis of event-related potentials. Front. Hum. Neurosci. 2014, 8, 213. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Oostenveld, R.; Fries, P.; Maris, E.; Schoffelen, J.M. FieldTrip: Open source software for advanced analysis of MEG, EEG, and invasive electrophysiological data. Comput. Intell. Neurosci. 2011, 2011, 1–9. [Google Scholar] [CrossRef] [PubMed]
  49. Navid, M.S.; Niazi, I.K.; Lelic, D.; Drewes, A.M.; Haavik, H. The Effects of Filter’s Class, Cutoff Frequencies, and Independent Component Analysis on the Amplitude of Somatosensory Evoked Potentials Recorded from Healthy Volunteers. Sensors (Basel, Switzerland) 2019, 19, 2610. [Google Scholar] [CrossRef] [Green Version]
  50. Bigdely-Shamlo, N.; Mullen, T.; Kothe, C.; Su, K.-M.; Robbins, K.A. The PREP pipeline: Standardized preprocessing for large-scale EEG analysis. Front. Neuroinformatics 2015, 9, 16. [Google Scholar] [CrossRef]
  51. Pascual-Marqui, R.D. Standardized low-resolution brain electromagnetic tomography (sLORETA): Technical details. Methods Find. Exp. Clin. Pharmacol. 2002, 24 (Suppl. D), 5–12. [Google Scholar]
  52. Mazziotta, J.; Toga, A.; Evans, A.; Fox, P.; Lancaster, J.; Zilles, K.; Woods, R.; Paus, T.; Simpson, G.; Pike, B.; et al. A probabilistic atlas and reference system for the human brain: International Consortium for Brain Mapping (ICBM). Philos. Trans. R. Soc. Lond. Ser. B Biol. Sci. 2001, 356, 1293–1322. [Google Scholar] [CrossRef]
  53. Pascual-Marqui, R.D.; Esslen, M.; Kochi, K.; Lehmann, D. Functional imaging with low-resolution brain electromagnetic tomography (LORETA): A review. In Proceedings of the Methods and Findings in Experimental and Clinical Pharmacology; Prous Science: Barcelona, Spain, 2002; Volume 24, pp. 91–95. [Google Scholar]
  54. Luck, S.J. An Introduction to the Event-related Potential Technique; MIT Press: Cambridge, MA, USA, 2014; ISBN 9780262525855. [Google Scholar]
  55. Onton, J. High-frequency broadband modulation of electroencephalographic spectra. Front. Hum. Neurosci. 2009, 3, 61. [Google Scholar] [CrossRef] [Green Version]
  56. Delorme, A.; Palmer, J.; Onton, J.; Oostenveld, R.; Makeig, S. Independent EEG sources are dipolar. PLoS ONE 2012, 7, e30135. [Google Scholar] [CrossRef] [PubMed]
  57. Jung, T.P.; Makeig, S.; Humphries, C.; Lee, T.W.; Mckeown, M.J.; Iragui, V.; Sejnowski, T.J. Removing electroencephalographic artifacts by blind source separation. Psychophysiology 2000, 37, 163–178. [Google Scholar] [CrossRef] [PubMed]
  58. Chaumon, M.; Bishop, D.V.M.; Busch, N.A. A practical guide to the selection of independent components of the electroencephalogram for artifact correction. J. Neurosci. Methods 2015, 250, 47–63. [Google Scholar] [CrossRef]
  59. Maris, E.; Oostenveld, R. Nonparametric statistical testing of EEG- and MEG-data. J. Neurosci. Methods 2007, 164, 177–190. [Google Scholar] [CrossRef] [PubMed]
  60. Nichols, T.E.; Holmes, A.P. Nonparametric permutation tests for functional neuroimaging: A primer with examples. Hum. Brain Mapp. 2002, 15, 1–25. [Google Scholar] [CrossRef] [Green Version]
  61. Allen, M.; Poggiali, D.; Whitaker, K.; Marshall, T.R.; Kievit, R.A. Raincloud plots: A multi-platform tool for robust data visualization. Wellcome Open Res. 2019, 4, 63. [Google Scholar] [CrossRef] [Green Version]
  62. Finnigan, S.P.; Walsh, M.; Rose, S.E.; Chalk, J.B. Quantitative EEG indices of sub-acute ischaemic stroke correlate with clinical outcomes. Clin. Neurophysiol. 2007, 118, 2525–2532. [Google Scholar] [CrossRef]
  63. Cebolla, A.M.; Palmero-Soler, E.; Dan, B.; Cheron, G. Frontal phasic and oscillatory generators of the N30 somatosensory evoked potential. NeuroImage 2011, 54, 1297–1306. [Google Scholar] [CrossRef]
  64. Waberski, T.D.; Buchner, H.; Perkuhn, M.; Gobbelé, R.; Wagner, M.; Kücker, W.; Silny, J. N30 and the effect of explorative finger movements: A model of the contribution of the motor cortex to early somatosensory potentials. Clin. Neurophysiol. 1999, 110, 1589–1600. [Google Scholar] [CrossRef]
  65. Balzamo, E.; Marquis, P.; Chauvel, P.; Régis, J. Short-latency components of evoked potentials to median nerve stimulation recorded by intracerebral electrodes in the human pre- and postcentral areas. Clin. Neurophysiol. 2004, 115, 1616–1623. [Google Scholar] [CrossRef]
  66. Rossi, S.; Della Volpe, R.; Ginanneschi, F.; Ulivelli, M.; Bartalini, S.; Spidalieri, R.; Rossi, A. Early somatosensory processing during tonic muscle pain in humans: Relation to loss of proprioception and motor “defensive” strategies. Clin. Neurophysiol. 2003, 114, 1351–1358. [Google Scholar] [CrossRef]
  67. Andrew, D.; Haavik, H.; Dancey, E.; Yielder, P.; Murphy, B. Somatosensory evoked potentials show plastic changes following a novel motor training task with the thumb. Clin. Neurophysiol. 2015, 126, 575–580. [Google Scholar] [CrossRef] [PubMed]
  68. Treleaven, J. Sensorimotor disturbances in neck disorders affecting postural stability, head and eye movement control. Man. Ther. 2008, 13, 2–11. [Google Scholar] [CrossRef] [PubMed]
  69. Passmore, S.R.; Murphy, B.; Lee, T.D. The origin, and application of somatosensory evoked potentials as a neurophysiological technique to investigate neuroplasticity. J. Can. Chiropr. Assoc. 2014, 58, 170–183. [Google Scholar] [PubMed]
  70. Leahy, R.M.; Mosher, J.C.; Spencer, M.E.; Huang, M.X.; Lewine, J.D. A study of dipole localization accuracy for MEG and EEG using a human skull phantom. Electroencephalogr. Clin. Neurophysiol. 1998, 107, 159–173. [Google Scholar] [CrossRef]
  71. Cuffin, B.N. EEG localization accuracy improvements using realistically shaped head models. IEEE Trans. Bio-Med. Eng. 1996, 43, 299–303. [Google Scholar] [CrossRef]
  72. Waberski, T.D.; Buchner, H.; Lehnertz, K.; Hufnagel, A.; Fuchs, M.; Beckmann, R.; Rienäcker, A. Properties of advanced headmodelling and source reconstruction for the localization of epileptiform activity. Brain Topogr. 1998, 10, 283–290. [Google Scholar] [CrossRef]
  73. Wang, L.; Yu, C.; Chen, H.; Qin, W.; He, Y.; Fan, F.; Zhang, Y.; Wang, M.; Li, K.; Zang, Y.; et al. Dynamic functional reorganization of the motor execution network after stroke. Brain 2010, 133, 1224–1238. [Google Scholar] [CrossRef] [Green Version]
  74. Haavik, H.; Özyurt, M.G.; Niazi, I.K.; Holt, K.; Nedergaard, R.W.; Yilmaz, G.; Türker, K.S. Chiropractic Manipulation Increases Maximal Bite Force in Healthy Individuals. Brain Sci. 2018, 8, 76. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Resting-state frequency-domain analysis. Scalp topographies of the difference of grand averaged power from (A) the pre-intervention sessions (chiropractic-control), (B) post- and pre-control sessions, (C) post- and pre-chiropractic sessions, and (D) between interventions ((C) minus (B)). The cluster-based permutation tests showed non-significant differences in comparisons in (AC), however, the absolute power in all frequency bands across the scalp was higher after the chiropractic intervention.
Figure 1. Resting-state frequency-domain analysis. Scalp topographies of the difference of grand averaged power from (A) the pre-intervention sessions (chiropractic-control), (B) post- and pre-control sessions, (C) post- and pre-chiropractic sessions, and (D) between interventions ((C) minus (B)). The cluster-based permutation tests showed non-significant differences in comparisons in (AC), however, the absolute power in all frequency bands across the scalp was higher after the chiropractic intervention.
Brainsci 10 00253 g001
Figure 2. Somatosensory evoked potentials (SEPs) amplitude. Dots represent (A) N20 and (B) N30 amplitudes from all sessions of the analyzed subjects. Boxplots show the median, 25th and 75th percentiles. The error bars represent mean ± 95% CI. The distribution plots show the density distribution estimated by a Gaussian kernel with SD of 1.5. (A) The N20 amplitude was not affected by either intervention. (B) The N30 amplitude was similar for the pre-intervention sessions. After the chiropractic spinal manipulation, the N30 amplitude increased significantly, however, it was not changed after the control intervention. The figure was made using the code provided by Allen et al. [61].
Figure 2. Somatosensory evoked potentials (SEPs) amplitude. Dots represent (A) N20 and (B) N30 amplitudes from all sessions of the analyzed subjects. Boxplots show the median, 25th and 75th percentiles. The error bars represent mean ± 95% CI. The distribution plots show the density distribution estimated by a Gaussian kernel with SD of 1.5. (A) The N20 amplitude was not affected by either intervention. (B) The N30 amplitude was similar for the pre-intervention sessions. After the chiropractic spinal manipulation, the N30 amplitude increased significantly, however, it was not changed after the control intervention. The figure was made using the code provided by Allen et al. [61].
Brainsci 10 00253 g002
Figure 3. SEPs latency. Dots represent (A) P14, (B) N20, (C) P22, and (D) N30 latencies from all sessions of the analyzed subjects. Boxplots show the median, 25th and 75th percentiles. The error bars represent mean ± 95% CI. The distribution plots show the density distribution estimated by a Gaussian kernel with SD of 1.5. None of the SEPs latencies were significantly affected by either intervention.
Figure 3. SEPs latency. Dots represent (A) P14, (B) N20, (C) P22, and (D) N30 latencies from all sessions of the analyzed subjects. Boxplots show the median, 25th and 75th percentiles. The error bars represent mean ± 95% CI. The distribution plots show the density distribution estimated by a Gaussian kernel with SD of 1.5. None of the SEPs latencies were significantly affected by either intervention.
Brainsci 10 00253 g003
Figure 4. Area under the curve. The dots represent adjusted area under curve (AUC) with respect to the baselines set at 0 (black vertical line) (i.e. post-AUC minus pre-AUC). There were no effects on the strength of brain sources by either intervention. Abbreviations are same as in Table 4.
Figure 4. Area under the curve. The dots represent adjusted area under curve (AUC) with respect to the baselines set at 0 (black vertical line) (i.e. post-AUC minus pre-AUC). There were no effects on the strength of brain sources by either intervention. Abbreviations are same as in Table 4.
Brainsci 10 00253 g004
Table 1. Patients’ characteristics.
Table 1. Patients’ characteristics.
No.Age (Years)Type of StrokeArea InvolvedAffected HemisphereFM ScoreTime Since Event (Months)
154IschemiaMCALeft5524
251IschemiaACALeft5718
368HemorrhageMCALeft4160
475IschemiaACARight7612
536IschemiaMCALeft6418
661HemorrhageMCALeft835
733IschemiaMCALeft555
848IschemiaMCALeft6424
956IschemiaMCALeft635
1058IschemiaMCALeft7625
1141HemorrhageMCARight6420
1262IschemiaMCARight7216
1346IschemiaMCARight5413
1433HemorrhageMCARight7146
1551IschemiaMCARight4623
1666IschemiaMCARight7812
17 *63HemorrhageACARight2150
1858IschemiaMCALeft6316
19 *38HemorrhageMCARight583
Note. * = Patients excluded; MCA = Middle cerebral artery; ACA = Anterior cerebral artery; FM = Fugl-Meyer Score.
Table 2. DAR and BSI descriptive statistics.
Table 2. DAR and BSI descriptive statistics.
InterventionSessionDARBSI1-25 HzBSIdeltaBSIthetaBSIalphaBSIbetaBSIgamma
ValuedirValuedirValuedirValuedirValuedirValuedir
ControlPre2.37 ± 1.700.44 ± 0.170.08 ± 0.240.48 ± 0.170.11 ± 0.250.48 ± 0.150.11 ± 0.270.47 ± 0.210.03 ± 0.260.41 ± 0.160.01 ± 0.220.46 ± 0.100.07 ± 0.22
Post2.12 ± 1.550.42 ± 0.160.08 ± 0.230.47 ± 0.170.11 ± 0.240.48 ± 0.160.12 ± 0.270.44 ± 0.210.02 ± 0.250.39 ± 0.14−0.01 ± 0.210.47 ± 0.140.04 ± 0.28
ChiropracticPre2.48 ± 2.200.38 ± 0.110.11 ± 0.140.43 ± 0.150.15 ± 0.140.41 ± 0.140.15 ± 0.160.40 ± 0.170.04 ± 0.180.35 ± 0.08−0.01 ± 0.140.40 ± 0.090.02 ± 0.21
Post2.05 ± 1.140.44 ± 0.210.09 ± 0.190.45 ± 0.180.11 ± 0.190.49 ± 0.240.11 ± 0.190.51 ± 0.250.06 ± 0.220.38 ± 0.210.01 ± 0.170.41 ± 0.160.03 ± 0.16
Table 3. ANOVA–N30 Amplitude.
Table 3. ANOVA–N30 Amplitude.
Predictordf_numdf_denSS_numSS_denFp Value
(Intercept)116501.38112.1071.560.00
intervention1160.587.611.210.30
session1160.572.293.990.06
intervention × session1161.203.375.710.03
Note. df_num indicates degrees of freedom numerator. df_den indicates degrees of freedom denominator. SS_num indicates sum of squares numerator. SS_den indicates sum of squares denominator.
Table 4. Coordinates (Talairach) and areas under curve (AUCs) of brain sources.
Table 4. Coordinates (Talairach) and areas under curve (AUCs) of brain sources.
RegionControlChiropractic
XYZAUC PreAUC PostXYZAUC PreAUC Post
SI26 ± 11−31 ± 938 ± 7897 ± 492941 ± 48227 ± 6−29 ± 1136 ± 71001 ± 589949 ± 473
Pre Frontal3 ± 1344 ± 316 ± 5349 ± 208383 ± 226−4 ± 1844 ± 611 ± 5376 ± 276418 ± 260
Cingulate−22 ± 27−5 ± 10−6 ± 4537 ± 354474 ± 383−18 ± 31−4 ± 11−6 ± 7530 ± 412494 ± 281
Cont. SII27 ± 8−52 ± 9−14 ± 3449 ± 236515 ± 31427 ± 12−49 ± 13−14 ± 2511 ± 328518 ± 286
Ips. SII−28 ± 10−52 ± 13−9 ± 12362 ± 172348 ± 172−28 ± 8−48 ± 15−12 ± 11336 ± 191336 ± 172
Note. Abbreviations: AUC = Area under curve; SI = Primary somatosensory cortex; Cont. SII = Contralateral secondary somatosensory cortex; Ips. SII = Ipsilateral secondary somatosensory cortex.

Share and Cite

MDPI and ACS Style

Navid, M.S.; Niazi, I.K.; Lelic, D.; Nedergaard, R.B.; Holt, K.; Amjad, I.; Drewes, A.M.; Haavik, H. Investigating the Effects of Chiropractic Spinal Manipulation on EEG in Stroke Patients. Brain Sci. 2020, 10, 253. https://0-doi-org.brum.beds.ac.uk/10.3390/brainsci10050253

AMA Style

Navid MS, Niazi IK, Lelic D, Nedergaard RB, Holt K, Amjad I, Drewes AM, Haavik H. Investigating the Effects of Chiropractic Spinal Manipulation on EEG in Stroke Patients. Brain Sciences. 2020; 10(5):253. https://0-doi-org.brum.beds.ac.uk/10.3390/brainsci10050253

Chicago/Turabian Style

Navid, Muhammad Samran, Imran Khan Niazi, Dina Lelic, Rasmus Bach Nedergaard, Kelly Holt, Imran Amjad, Asbjørn Mohr Drewes, and Heidi Haavik. 2020. "Investigating the Effects of Chiropractic Spinal Manipulation on EEG in Stroke Patients" Brain Sciences 10, no. 5: 253. https://0-doi-org.brum.beds.ac.uk/10.3390/brainsci10050253

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop