# The Effects of Spinal Manipulation on Motor Unit Behavior

Lucien Robinault, Aleš Holobar, Sylvain Cremoux, Usman Rashid, Imran Khan Niazi, Kelly Holt, Jimmy Lauber, Heidi Haavik

► To cite this version:

Lucien Robinault, Aleš Holobar, Sylvain Cremoux, Usman Rashid, Imran Khan Niazi, et al. The Effects of Spinal Manipulation on Motor Unit Behavior. Brain Sciences, 2021, 11 (1), pp.105. 10.3390/brainsci11010105. [https://doi.org/10.3390/brainsci11010105](https://doi.org/10.3390/brainsci11010105). [hal-03424918](https://HAL Id: hal-03424918)

Submitted on 21 Apr 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
Article

# The Effects of Spinal Manipulation on Motor Unit Behavior

Lucien Robinault<fnref n="1" />, Aleš Holobar<fnref n="2" /> ID, Sylvain Crémoux<fnref n="3" />,<fnref n="4" /> ID, Usman Rashid<fnref n="5" /> ID, Imran Khan Niazi <fnref n="5" />,<fnref n="6" />,<fnref n="7" />,* ID,
Kelly Holt<fnref n="6" /> ID, Jimmy Lauber<fnref n="1" /> ID and Heidi Haavik<fnref n="6" /> ID

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Robinault, L.; Holobar, A.; Crémoux, S.; Rashid, U.; Niazi, I.K.; Holt, K.; Lauber, J.; Haavik, H. The Effects of Spinal Manipulation on Motor Unit Behavior. *Brain Sci.* **2021**, *11*, 105. https://doi.org/10.3390/brainsci11010105

Received: 16 December 2020

Accepted: 12 January 2021

Published: 14 January 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** Over recent years, a growing body of research has highlighted the neural plastic effects of spinal manipulation on the central nervous system. Recently, it has been shown that spinal manipulation improved outcomes, such as maximum voluntary force and limb joint position sense, reflecting improved sensorimotor integration and processing. This study aimed to further evaluate how spinal manipulation can alter neuromuscular activity. High density electromyography (HD sEMG) signals from the tibialis anterior were recorded and decomposed in order to study motor unit changes in 14 subjects following spinal manipulation or a passive movement control session in a crossover study design. Participants were asked to produce ankle dorsiflexion at two force levels, 5% and 10% of maximum voluntary contraction (MVC), following two different patterns of force production ("ramp" and "ramp and maintain"). A significant decrease in the conduction velocity (p = 0.01) was observed during the "ramp and maintain" condition at 5% MVC after spinal manipulation. A decrease in conduction velocity suggests that spinal manipulation alters motor unit recruitment patterns with an increased recruitment of lower threshold, lower twitch torque motor units.

**Keywords:** high-density surface electromyography; chiropractic; electromyography decomposition; motor unit

## 1. Introduction

Over recent years, a growing body of research has investigated the neurophysiological effects of spinal manipulation [1,2].

These studies revealed that spinal manipulation improves outcomes, such as maximum voluntary contraction and joint position sense [3–7], and that these changes are associated with central plastic adaptations in sensorimotor integration at the cortical level [3,7–11].

Two previous studies have reported that participants receiving a single session of spinal manipulation [5], or 12 weeks of chiropractic care [6], showed an improvement in elbow or ankle joint position sense, respectively, in comparison to control participants. This suggests that spinal manipulation has an impact on the integration and processing of somatosensory information from the limbs. Other studies have used electroencephalography (EEG) to measure sensory integration adaptation at the cortical level after spinal manipulation and have reported that early cortical somatosensory evoked potentials (SEPs), such

<footnote n="1">¹ Laboratoire d'Automatique, de Mécanique et d'Informatique Industrielles et Humaines (LAMIH UMR CNRS 8201), Université Polytechnique Hauts-de-France, F-59313 Valenciennes, France; lucien.robinault@uphf.fr (L.R.); jimmy.lauber@uphf.fr (J.L.)</footnote>

<footnote n="2">² Faculty of Electrical Engineering and Computer Science, University of Maribor, SI-2000 Maribor, Slovenia; ales.holobar@um.si</footnote>

<footnote n="3">³ Centre de Recherche Cerveau et Cognition, Université de Toulouse, UPS, 31052 Toulouse, France; sylvain.cremoux@cnrs.fr</footnote>

<footnote n="4">⁴ The Brain and Cognition Research Center (CerCo), CNRS UMR 5549, 31052 Toulouse, France</footnote>

<footnote n="5">⁵ Health and Rehabilitation Research Institute, AUT University, Auckland 0627, New Zealand; usman.rashid@aut.ac.nz</footnote>

<footnote n="6">⁶ Centre for Chiropractic Research, New Zealand College of Chiropractic, Auckland 1060, New Zealand; kelly.holt@nzchiro.co.nz (K.H.); heidi.haavik@nzchiro.co.nz (H.H.)</footnote>

<footnote n="7">⁷ Department of Health Science and Technology, Aalborg University, 9220 Aalborg, Denmark</footnote>

<footnote>* Correspondence: imran.niazi@nzchiro.co.nz; Tel.: +64-9-526-6789; Fax: +64-9-526-6788</footnote>
as N20 and N30, significantly decreased after cervical spinal manipulation in a relatively healthy population [12–14]. These early cortical SEP peaks are thought to represent the arrival of the afferent signals at the primary somatosensory cortex and early sensorimotor integration, respectively [15,16]. Interestingly, in a chronic stroke population, who are known to have a reduced amplitude of early SEP peaks [17], a recent study showed an increase in N30 SEP peak amplitude following chiropractic care [18]. These observed changes in SEP peak amplitude after chiropractic care indicate that spinal manipulation can alter cortical somatosensory processing and results in neural plastic changes in the central nervous system (CNS) [11].

Motor control changes following spinal manipulation have also previously been documented [4]. One study showed that spinal manipulation altered cortical motor control of two upper limb muscles in a muscle-specific manner [14]. This study used a paired-pulse transcranial magnetic stimulation (TMS) protocol to explore specific, central corticomotor facilitatory, and inhibitory neural pathways to the two upper limb target muscles. After spinal manipulation, there was a decrease in short interval intracortical inhibition, an increase in intracortical facilitation, and a shortening of the cortical silent period in the flexor muscle, whereas the extensor muscle showed opposite changes, with a decrease in intracortical facilitation and a lengthening of the cortical silent period [14]. Another TMS study measuring stimulus response curves found that spinal manipulation led to changes in corticospinal excitability, as measured by a significantly larger motor evoked potential for both upper and lower limb muscles, and larger amplitudes of movement-related cortical potential (MRCP) components with no change in spinal F-wave amplitude or persistence [10], suggesting these changes occurred at the supra-spinal level [16].

Overall, these studies suggest that spinal manipulation alters central corticomotor processing. Interestingly, these central motor adaptations also influence motor behavior, as several studies have shown significant increases in maximum voluntary plantarflexion force production after a single session of spinal manipulation [3,7,19]. The first of these studies showed that the force increase was associated with a large increase in the cortical-based V-wave and a small, but significant, decrease in the spinal H-reflex threshold, respectively, reflecting supraspinal and spinal cord adaptive changes in α-motor neuron excitability [19]. Two follow-up studies in elite athletes [3] and chronic stroke patients [3] also showed increased ankle plantarflexion force production with a large V-wave increase following spinal manipulation, but no changes in the H-reflex. Overall, these results suggest that spinal manipulation alters motor preparation and commands from the supraspinal regions to more efficiently control force production. Combined, these findings suggest that spinal manipulation has an impact on central cortical processing that improves the accuracy with which the brain is aware of limb position and alters the way the brain controls the upper and lower limbs, leading to more efficient motor control [3–5,7]. Yet, how these central neural adaptations translate at the neuromuscular level remains to be studied.

Neuromuscular changes can be studied via the assessment of motor unit (MU) activity. The gold standard to study MU properties is the decomposition of intramuscular electromyographic (EMG) signals [20,21], but this technique is invasive and limits the area of the muscle that can be studied [22]. To overcome these drawbacks, high-density surface EMG (HD sEMG) has been developed, together with advanced signal processing techniques [23–29]. Numerous variables can be estimated after decomposition of the HD sEMG signals into MU action potentials (MUAP), such as MU firing rate, instantaneous MU firing rate gradient, and changes in the propagation velocity of the MUAP across the muscle fibers [30–32]. The MU firing rate refers to the instantaneous rate of coding of the MU [33]. The maximum firing rate value is known as the peak firing rate, while the instantaneous MU firing rate gradient refers to the derivative of the firing rate. The conduction velocity refers to the speed at which the neural impulse propagates from the motor unit through the muscle fibers [16].

The previous studies that have shown changes in ankle plantarflexion force, V-waves, and the H-reflex following spinal manipulation [3,7,19] have been conducted using a small
number of surface EMG electrodes. This enables few conclusions to be made about changes in MU behavior associated with spinal manipulation [3,7,19]. To better understand how spinal manipulation influences MU behavior and neuromuscular control, HD sEMG can be used in place of standard surface electrodes [23,34]. The changes in motor control following spinal manipulation that have been observed in previous studies could be due to a shift towards activation of lower threshold MU's [4]. If this is the case, and because conduction velocity of a MU is highly correlated to its twitch torque production, one would expect to see lower conduction velocity and higher peak firing rate values and an increase in instantaneous MU firing rate gradient after spinal manipulation [35]. Thus, to further explore the changes in motor unit behavior that occur after spinal manipulation, this study aimed to investigate neuromuscular adaptation of the tibialis anterior (TA) muscle before and after a single session of spinal manipulation in a subclinical pain population.

## 2. Materials and Methods

### 2.1. *Experimental Design and Setting*

The experiment was conducted as a randomized, crossover study, where each subject participated in both control and intervention groups (see Figure 1).

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Illustration of the study design.</figcaption></figure>

The experiment took place at Aalborg University Hospital, Aalborg, Denmark. Before the experiment, the subjects were allocated randomly to either receive the chiropractic intervention or sham intervention first. Participants were then assessed before and after receiving the appropriate intervention. After a wash out period of at least 24 h, the participants were retested before and after receiving the alternate intervention. The study was conducted according to the Declaration of Helsinki. The North Denmark Region Committee on Health Research Ethics (N-20140027) approved the study. All participants gave written informed consent before participating in the study.

### 2.2. Participants

Fourteen participants (10 males and 4 females; age: 28.6 ± 2.8 years) were recruited via convenience sampling for this study. Participants were included if they presented with a history of mild back pain or neck pain for which they had not previously sought treatment and had evidence of mild spinal dysfunction based on a chiropractic assessment. Participants fulfilling these criteria can be referred to as having subclinical pain [8]. Subclinical pain (SCP) refers to recurring spinal dysfunction, such as mild pain, ache, and/or stiffness, with or without a history of known spinal trauma. Individuals with SCP do not have constant symptoms and have not yet sought treatment of their spinal problems. There is an increasing interest in SCP in the literature [5,8,36–43] because individuals that fall into this category provide an opportunity to explore neurophysiologic dysfunction without the confounding effect of current pain, nor the confounding effects of previous drug use or surgery. Changes in pain is known to alter sensory processing and motor control [44], as would various drugs and surgery. Studying this group may, therefore, help to gain a better understanding of the features of those developing more serious chronic pain states
and may also help characterizing which subgroups of those developing pain problems that
may respond to different sorts of treatment options.

On the day of the experiment, the participants were excluded from the study if they presented with a pain level greater than 2 on a 10-points pain scale [45]. Sample size calculations and predicted effect sizes were based on changes observed in a previous study that investigated similar neurophysiological changes before and after spinal manipulation (SM) in previous study Niazi et al. [19]. For an effect size of 0.5, 11 subjects were needed to be able to reject the null hypothesis that the population means of the experimental and control groups are equal with probability (power) 0.9. To allow for drop out during the trial and relative uncertainty relating to power outcomes, we aimed to enroll 15 subjects in the trial but were able to recruit 14.

2.3. *Spinal Manipulation and Sham Intervention*

The chiropractor assessed the entire spine and both sacroiliac joints for segmental dysfunction (also called vertebral subluxations by chiropractors). The clinical indicators that were used to assess joint function included assessing for tenderness to palpation of the relevant joints, manually palpating for restricted range of motion, assessing for palpable asymmetric muscle tension, and any abnormal or blocked joint play or end-feel of the joints. All of these biomechanical characteristics are used by chiropractors as clinical indicators of joint dysfunction [46].

For the spinal manipulation intervention, a New Zealand registered chiropractor checked participants for spinal dysfunction known as chiropractic or vertebral subluxations and adjusted them, where necessary. Vertebral subluxations are recognized as biomechanical lesions of the vertebral column by the World Health Organization (ICD-10-CM code M99.1) [47]. Due to the past two decades of basic science research, they are now viewed 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 (CNS's) ability to self-regulate, self-organize, adapt, repair, and heal [48]. It is thought that such "maladaptation" of body posture may be beneficial and is used to avoid further pain from the region (pain adaptation concept) [49]; however, when maintained for a long period of time, they become maladaptive or harmful. The clinical indicators that were used to assess joint function, and to identify vertebral subluxations, included assessing for tenderness to palpation of the relevant joints, manually palpating for restricted range of motion, assessing for palpable asymmetric muscle tension, and any abnormal or blocked joint play or end-feel of the joints. These clinical indicators are routinely used by chiropractors when analyzing the spine and have previously been shown to be reliable for the identification of chiropractic subluxations when used as a multidimensional battery of tests [46,49].

During the chiropractic intervention, the chiropractor manipulated the areas of joint dysfunction using high velocity, low amplitude thrusts. All of the spinal manipulations carried out in this study were manual high-velocity, low-amplitude thrusts to the spine or pelvic joints [50]. This is a standard manipulation technique used by chiropractors, and is also referred to as spinal adjustments. The mechanical properties of this intervention have been investigated; and, although the actual force applied to the subject's spine depends on the chiropractor, the patient, and the spinal location of the manipulation, the general shape of the force-time history of spinal manipulations is very consistent [51], and the duration of the thrust is always less than 200 milliseconds [52]. This high-velocity type of manipulation was specifically chosen because previous research has shown that reflex EMG activation is observed after high-velocity, low-amplitude manipulations but not lower-velocity mobilizations [53]. This manipulation technique has also been previously used in studies that have investigated the neurophysiological effects of spinal manipulation [4].

The control intervention was intended to act as a control treatment session, as well as to
act as a physiological control for possible changes occurring due to the cutaneous, muscular,
or vestibular input that would occur with the type of passive and active movements
involved in preparing a subject/patient for a manipulation. During the control intervention, the subject's head and/or spine and peripheral joints were moved passively and actively in ways similar to what was done by the chiropractor that provided the actual chiropractic intervention. Thus, the control intervention involved the subjects being moved into the manipulation setup positions similar to how the chiropractor would normally setup a subject prior to applying a manipulative thrust. Loading a joint, as is done prior to spinal manipulation, has been shown to alter paraspinal proprioceptive firing in anesthetized cats [54] and, therefore, was carefully avoided by ending the movement prior to end range of motion when passively moving the subjects. No spinal manipulation was performed during any control intervention.

For each subject, the chiropractor performing the real and control interventions was the same and had at least 10 years of clinical experience.

## 2.4. Procedure

The foot of the participant was placed in an ankle ergometer, at a neutral angle, with a strap to fix the foot in place at the metatarsal level. At the beginning and the end of the experimental protocol, the participants were asked to perform 3 dorsiflexion maximal voluntary contractions (MVCs) of 3 s, with 2 min rest in between each MVC.

The participants were then asked to produce dorsiflexion following a triangular "ramp" or trapezoidal "ramp and maintain" force production pattern, as described in Figure 2. The "ramp" condition consisted of a 3 s increase of force production to reach the target force intensity, then 3 s ramping down to rest. The "ramp and maintain" condition consisted of a ramping up phase of 1 s to the target force intensity, then maintaining this intensity for 6 s, followed by a ramping down phase of 1 s to rest. The force intensity was either 5% or 10% of the MVC recorded at the beginning of the experiment. The reasons to use low % MVC was to specifically investigate the low threshold single motor units because only the threshold of the H-reflex (i.e., reflecting low threshold motor units) was shown to change following a spinal manipulation intervention by Niazi et al. [19].

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Force pattern for the subject to follow. (a) "Ramp" condition; (b) "ramp and maintain" condition.</figcaption></figure>

The purpose of the "ramp" condition was to find the firing threshold of the detected MUs and the "ramp & maintain" condition was used to "isolate" MU's with different firing thresholds [55-57]. The force intensity of 5% and 10% was chosen as the HD sEMG decomposition techniques used are reliable at these force levels [58] and because previous research suggested, via indirect measures, that spinal manipulation is more likely to have an impact on lower threshold MUs [19].

Before and after the interventions, the subjects executed 3 contractions following the "ramp" condition pattern at the first intensity of contractions, and then at the second intensity of contractions. The first force intensity was randomly determined for each participant using computer-generated random allocation. The subjects then executed
3 contractions following the "ramp and maintain" condition pattern at each intensity of contractions. Contractions were separated by a 10 s rest period. Conditions were separated by a 5-min break.

## 2.5. *Equipment*

### 2.5.1. Ankle Ergometer

Force was recorded from a force transducer [19] mounted on a custom-made pedal. The force data were also displayed on a screen in front of the subject and used in real time to provide the target pattern and force level that the subject had to develop. The force was sampled at 2000 Hz and acquired with 'Mr. Kick' software (SMI, Aalborg University, Aalborg, Denmark).

### 2.5.2. HD EMG

EMG signals from the tibialis anterior (TA) muscle were recorded with semi-disposable adhesive 64 sensor grid matrix electrodes (ELSCH064R3, OT Bioelettronica, Torino, Italy) coupled to an amplifier (EMG-USB2+, OT Bioelettronica). The grid was applied on the dominant leg of the subject after the skin had been abraded with a paste (Meditec-Every, Parma, Italy). The position of the *missing* electrode was at the most distal and medial position on the leg, with the columns of the grid in the direction of the muscle fibers, as shown in Figure 3. The grid was fixed on the skin by adhesive tape, and the reference electrode was placed on the ankle of the subject. The EMG signal was acquired at a frequency of 2048 Hz, gain of 2000, band-pass filter 10–500 Hz, and analog to digital (A/D) converted with a resolution of 12 bits. All the signals were acquired via OTBiolab+ software (OT Bioelettronica) [59]. An example of recorded EMG signals can be seen in Figure 4.

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Placement of the ELSCH064R3 electrode grid; note the missing electrode on the top right corner.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Representative electromyographic (EMG) signals from the first row of electrodes of the ELSCH064R3 (OT-Bioelectronica, Italy) grid, acquired on subject 2, before control intervention during a “ramp” condition task at 10% maximum voluntary contraction (MVC).</figcaption></figure>

## 2.6. *Data Processing*

Recorded sEMG signals were processed using DEMUSE software (Maribor, Slovenia), [60] where they were off-line bandpass filtered (10–500 Hz), then decomposed with the convolution kernel compensation (CKC) technique [34,61–64]. For each contraction, 60 of the 63 acquired monopolar signals were used for the decomposition. The signals were treated separately for each contraction of each condition, each recording having a small buffer time before and after the actual contraction (the total duration of each recording was between 12 s to 15 s, and had no influence on the decomposition process). This decomposition technique has been validated extensively with both simulated and experimental signals, and the detailed analysis procedure has been described in previous studies [65,66]. Once decomposed, the firing pattern of different MUs were identified, manually inspected and edited by an experienced operator using the DEMUSE tool interface [60]. Movement artifacts and abnormally high (>30 Hz) or low (<5 Hz) instantaneous firing rate (IFR) values were discarded. Duplicates of detected MUs were also discarded. The signals that presented a pulse to noise ratio (PNR) of 30 dB or above were kept [67]. Finally, MUs with an extremely low number of firing (<15) were discarded. MUs were not tracked across the contractions, therefore, labels of MUs that were detected during the “ramp” or “ramp and maintain” contractions did not correspond to those obtained during another contraction for individual subjects. For example, MU 1 detected for one contraction was not the same as MU 1 detected for another contraction.

From those firing patterns, the following variables were computed:

*   **Conduction velocity (CV):** The speed which an action potential travels along the membrane of a skeletal muscle fiber [68]. The estimation of the conduction velocity was calculated using the techniques described by Farina et al. [69,70], which use a maximum likelihood multiple-channel method, using the Newton method for efficient optimization. Coefficient of correlation (CC) was used to measure the goodness of match of MUAPs on neighboring HDsEMG channels and, thus, the accuracy of CV estimation. Only the CV values with corresponding CC ≥0.7 were used for the statistical analysis.
* Peak firing rate (PFR): The PFR of a MU was defined as the maximum instantaneous firing rate value for the MU. The instantaneous firing rate (IFR) of each MU was calculated:

$$IFR_n = \frac{1}{(t_{n+1} - t_n)}, \qquad (1)$$

with $IFR_n$ in pps, and $t_{n+1}$ and $t_n$ in seconds.

An example of the MU firing rate is shown in Figure 5.

<figure><img src="image_8.png" /><figcaption>Figure 5. Motor unit (MU) firing time of the detected MUs for subject 2. (a) During the “ramp” condition at 10% MVC after control intervention. (b) During the “ramp and maintain” condition at 5% MVC after control intervention. A force pattern is depicted in orange.</figcaption></figure>

* Instantaneous MU firing rate gradient (IFRG): The IFRG is the slope of linear approximation of MU firing rate increase from recruitment to its peak IFR value (Figure 6).
The following technique was used to compute this value: a polynomial of the second order was fitted on the instantaneous firing rate curve of the MU, using the `polyfit()` function from MATLAB 2018a (MathWorks, Natick, MA, United States). Then, a linear function was fitted on the first value and the peak value of this polynomial. The slope of the linear expression defined the IFRG. In each participant, IFRG values were averaged across contraction repetitions and across condition and intensity before being statistically analyzed.

<figure><img src="image_9.png" /><figcaption>**Figure 6.** Representation of the instantaneous motor unit (MU) firing rate gradient (IFRG). The blue line depicts the 2nd order polynomial fitted to the IFR, pictured in red dots. From this polynomial, the peak IFR was retrieved (dotted green line) and then used to compute the linear expression from the first detected IFR value to the peak IFR (yellow line).</figcaption></figure>

## 2.7. Statistical Analysis

The statistical analysis was performed in R (R Foundation for Statistical Computing, Vienna, Austria) version 3.5.1 using *lme4* package version 1.1-21 [71,72] and *emmeans* package version 1.3.4. [73] Tukey's HSD method was used to perform pairwise contrasts. Linear mixed regression or generalized linear mixed regression models were setup to investigate the influence of the chiropractic intervention on the following variables compared to the control: (i) average peak firing rate of the MUs (PFR), (ii) average instantaneous firing rate gradient of the MUs (IFRG), and (iii) average conduction velocity (CV) of the muscle fibers associated with the MUs. For each dependent variable, a separate mixed model was setup to estimate its mean and standard error across the 2 intervention sessions (Intervention: chiropractic, control) under 2 intensity levels (Intensity: 5%, 10%) in the 2 conditions (Condition: "ramp", "ramp and maintain"). Each model was adjusted for the baseline (pre-intervention) values of the outcome variables and included a saturated fixed effects structure consisting of the independent variables and all their possible interactions. To cater for correlated repeated measures data from the same participants, the models also included random intercept effects for participants. PFR, IFRG, and CV were coded as continuous variables. Intensity, Intervention, and Condition were entered as categorical variables. A *p*-value of less than 0.05 was considered significant for all statistical tests.

## 3. Results

Fourteen participants (10 males and 4 females; age: 28.6 ± 2.8 years) participated in this study. Their baselines characteristics can be found in Table 1. Two participants (one male and one female) completed only a control session as they opted out of the study following their first assessment due to scheduling issues. Their results have, however, been included in the data processing and analysis.
<table><caption>**Table 1.** Baseline measurements before Intervention (mean ± SE).</caption><thead><tr><th></th><th>Control Group<br>(N = 14)</th><th>Intervention Group<br>(N = 12)</th></tr></thead><tbody><tr><td>Conduction velocity (m.s<sup>-1</sup>)</td><td>5.11 ± 0.62</td><td>5.12 ± 0.61</td></tr><tr><td>Peak firing rate (pps)</td><td>15.79 ± 4.04</td><td>15.45 ± 3.9</td></tr><tr><td>Instantaneous MU firing rate gradient (pps.s<sup>-1</sup>)</td><td>1.33 ± 1.48</td><td>1.37 ± 1.51</td></tr></tbody></table>

Across all the subjects, a total of 2240 (Pre) and 2329 (Post) motor units were decomposed from the Control intervention session and total of 1954 (Pre-) and 1962 (Post-) motor units (MUs) were decomposed during from spinal manipulation session. In control sessions, there were 14 participants; thus, the average was 2240/14 = 160, whereas, in the experiment session, the average was 1954/12 = 162.83. Thus, the average number of MU was higher in the intervention session. The number of MU discarded because they were showcasing less than 15 firing represent 4.7% of the total MU detected (399 MU discarded for a total of 8485).

Tables 2–4 show pair-wise contrasts across sessions for the two conditions: control and intervention for instantaneous MU firing rate gradient (IFRG) of MUs, peak pulse/min of Mus, and conduction velocity (CV).

<table><caption>**Table 2.** Estimated Group differences in the conduction velocity (CV) across groups.</caption><thead><tr><th>Condition, Intensity</th><th>Control—Intervention<br>Difference (Mean ± SE)</th><th>df</th><th>T Ratio</th><th>p Value</th></tr></thead><tbody><tr><td>Ramp and maintain, 5%</td><td>0.3 ± 0.11 m.s<sup>-1</sup></td><td>53.7</td><td>2.64</td><td>0.01</td></tr><tr><td>Ramp and maintain, 10%</td><td>0.2 ± 0.12 m.s<sup>-1</sup></td><td>53.8</td><td>1.65</td><td>0.11</td></tr><tr><td>Ramp, 5%</td><td>0.07 ± 0.11 m.s<sup>-1</sup></td><td>53.8</td><td>0.66</td><td>0.51</td></tr><tr><td>Ramp, 10%</td><td>0.03 ± 0.12 m.s<sup>-1</sup></td><td>54.4</td><td>0.28</td><td>0.78</td></tr></tbody></table>

<table><caption>**Table 3.** Estimated group differences in the average instantaneous MU firing rate gradient (IFRG) of MUs across groups.</caption><thead><tr><th>Condition, Intensity</th><th>Ratio Control/Intervention<br>Difference (Mean ± SE)</th><th>z-Value</th><th>p-Value</th></tr></thead><tbody><tr><td>Ramp and maintain, 5%</td><td>0.99 ± 0.15</td><td>-0.05</td><td>0.96</td></tr><tr><td>Ramp and maintain, 10%</td><td>0.95 ± 0.15</td><td>-0.32</td><td>0.75</td></tr><tr><td>Ramp, 5%</td><td>1.27 ± 0.2</td><td>1.52</td><td>0.13</td></tr><tr><td>Ramp, 10%</td><td>1.174 ± 0.19</td><td>0.98</td><td>0.33</td></tr></tbody></table>

<table><caption>**Table 4.** Estimated group differences in the average peak firing rate (PFR) of MUs across groups.</caption><thead><tr><th>Condition, Intensity</th><th>Control—Intervention<br>Difference (Mean ± SE)</th><th>z-Value</th><th>p-Value</th></tr></thead><tbody><tr><td>Ramp and maintain, 5%</td><td>0.286 ± 0.4 pps</td><td>0.715</td><td>0.47</td></tr><tr><td>Ramp and maintain, 10%</td><td>-0.24 ± 0.47 pps</td><td>-0.508</td><td>0.61</td></tr><tr><td>Ramp, 5%</td><td>0.14 ± 0.42 pps</td><td>0.33</td><td>0.74</td></tr><tr><td>Ramp, 10%</td><td>0.1 ± 0.52 pps</td><td>0.2</td><td>0.84</td></tr></tbody></table>

### 3.1. Conduction Velocity

The pair-wise contrasts across sessions for the control and chiropractic intervention indicated a statistically significant difference in the conduction velocity between control and spinal manipulation for “ramp and maintained” at 5% MVC only (0.3 ± 0.11 m.s⁻¹; *p* = 0.01). The changes were not significant for the other conditions. Figure 7 shows the estimated values for the “ramp and maintain” condition between the control and spinal manipulation sessions.
<figure><img src="image_10.png" /><figcaption>**Figure 7.** Distribution plot of the CV after intervention (control vs. spinal manipulation) for “ramp & maintain” condition adjusted to baseline.</figcaption></figure>

## 3.2. *Instantaneous MU Firing Rate Gradient & Peak Firing Rate*

No significant changes were found for the instantaneous MU firing rate gradient, and peak firing rate variables. Tables 3 and 4 summarize the results, respectively, for instantaneous MU firing rate gradient and peak firing rate. The following figures display the distribution plots of the results for instantaneous MU firing rate gradient for the condition “ramp & maintain” in Figure 8 and for peak firing rate for the condition “ramp & maintain” in Figure 9.
<figure><img src="image_11.png" /><figcaption>**Figure 8.** Magnified group distribution plot of the IFRG after intervention (control vs. spinal manipulation) for “ramp & maintain” condition adjusted to baseline. (Unmagnified version is in Appendix A)</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>**Figure 9.** Group distribution plot of the PFR after intervention (control vs. spinal manipulation) for “ramp & maintain” condition adjusted to baseline.</figcaption></figure>

## 4. Discussion

The results revealed a main effect of spinal manipulation on the conduction velocity of the MUs. This supports previous studies that have suggested that spinal manipulation has an impact on sensory and motor neural processing that impacts motor control [3-7,74]. This study aimed to evaluate the influence of spinal manipulation on MU activity and therefore, neuromuscular control. To achieve this aim, we used HD sEMG with advanced signal processing techniques to obtain information regarding MU behavior. Based on previous research involving spinal manipulation, we expected a decrease in conduction velocity and/or an increase in instantaneous MU firing rate gradient and peak firing rate following the chiropractic intervention. Each of these changes could explain neuromuscular adaptions previously documented after spinal manipulation [1-3,7-11].

The results showed significant changes only in the “ramp & maintain” condition at 5% of the MVC, where the conduction velocity decreased in the spinal manipulation session as compared to the control session. It is known that conduction velocity of an MU is correlated to its twitch torque production capacity [35]; therefore, the decrease of conduction velocity after spinal manipulation would reflect the recruitment of MU
producing, on average, lower twitch torque. This decrease in conduction velocity after spinal manipulation suggests that spinal manipulation increases the activation of lower threshold MUs which are better suited for the production of low force and precise action during low contraction level of the TA muscle.

A recent study [74] explored the effects of a single spinal manipulation session on maximum voluntary contractions (MVC) of the ankle dorsiflexors, HD EMG, intramuscular EMG, and near-infrared spectroscopy (NIRS) from the TA muscle in 25 participants with low level recurring spinal dysfunction. They found a significant increase in conduction velocity in both the isometric steady-state contractions at 10% of MVC and during the isometric ramp contractions at 10% of MVC compared to the control intervention [74]. The increased MVC in TA strength and increase in conduction velocity, without changes in motor unit discharge rate, suggests that the MVC strength changes observed following spinal manipulation in this study were mostly due to increased recruitment of larger, higher threshold motor units. It is interesting to note that this recent study found an increase in conduction velocity at 10% MVC following the spinal manipulation session, whereas, in the current study, using the HD sEMG decomposition methodology, we found a decrease in the conduction velocity at 5% MVC [74]. Both studies recorded EMG during similar “ramp” and “ramp & maintain” conditions at both 5% and 10% MVC [74]. This apparent contradiction in finding may be due to the difference in the way the conduction velocity across the TA was calculated. In a recent study [74], the conduction velocity was calculated using algorithms of cross-correlation and a maximum likelihood multi-channel approach adapted from Farina et al. [75]. This method involved iteratively cross-correlating signals in a multi-channel matrix to find the electrodes whose signals had the least mean square error while calculating the similarity between them. In the current study, we used the HDsEMG decomposition methodology and calculated the conduction velocity using the techniques described by Farina et al. in 2000 and 2001, respectively [69,70], which use a maximum likelihood multiple-channel method, using the Newton method for efficient optimization. Coefficient of correlation (CC) was used to measure the goodness of match of MUAPs on neighboring HDsEMG channels and, thus, the accuracy of CV estimation. Only the CV values with corresponding CC ≥0.7 were used for the statistical analysis. These methodological differences may explain the differences in conduction velocity changes following spinal manipulation. Alternatively, it could be individual variation in responses to the spinal manipulations themselves that could explain these differences. However, the most likely explanation is that the two different methods for calculating CV are using different strategies. One method is measuring CV from decomposed single Mus, and the other is measuring the CV of all MUs under the HDsEMG electrode grid. Regardless, it appears that spinal manipulation alters the central corticomotor control of the TA muscle, possibly by increased recruitment of larger, higher threshold motor units to enable greater levels of MVC, and/or by increases the activation of lower threshold MUs which are better suited for the production of low force and precise action as required by the 5% “ramp” and “ramp & maintain” conditions required of the subjects in this study.

Further studies should explore these finding in higher percentage of MVC contractions, and carefully consider any significant changes of the IFRG induced by spinal manipulation. Indeed, as it is developed in the theories on MU recruitment, be it the dominant theory of Kernell [76,77] or its main alternative proposed by De Luca and Contessa [78], the IFRG is correlated to the MUs excitation threshold; thus, the IFRG could yield valuable information. No significant results concerning the IFRG could be found in our study, but regarding the data acquired, it seems that the choice of the type of force production pattern could impact the IFRG as we can see opposite changes between the condition “ramp” and “ramp & maintain” were the IFRG seems to, respectively, increase or decrease. Therefore, the choice of the force production patterns should be carefully considered in future studies. It is interesting to note this trend for a task-specific change in IFRG between “ramp” and “ramp & maintain” conditions following spinal manipulation. A previous study has shown that spinal manipulation can have complete opposite neuromuscular changes to
two different functional muscles [14]. This study found that, after spinal manipulation, there was a decrease in short interval intracortical inhibition, an increase in intracortical facilitation, and a shortening of the cortical silent period in the flexor muscle, whereas the extensor muscle showed opposite changes, with a decrease in intracortical facilitation and a lengthening of the cortical silent period [14]. It may be that spinal manipulation can also alter corticospinal excitability to the same muscle in a different direction depending on the task required (i.e., “ramp” and “ramp & maintain” in the current study). This should be explored further in future studies.

Multiple previous studies have documented strength increases following a spinal ma-
nipulation intervention. One study demonstrated a 16% increase in plantar flexor muscle
strength after a single session of spinal manipulation in subclinical pain participants [19].
Another study found an 8% increase in plantar flexor muscle strength following a single
session of spinal manipulation in elite taekwondo athletes [3]. A third study was conducted
in a group that had lost their ability to cortically activate their muscles, i.e., chronic stroke
patients with ongoing plantar flexor muscle weakness which found that spinal manipu-
lation significantly increase their plantarflexion muscle strength by on average 64.2% [7].
This is rather remarkable and suggests spinal manipulation may be of great value for both
athletes and those who have lost their cortical ability to activate their muscles. This should
be further explored in future studies. The current study finding that showed changes
in the conduction velocity MU level suggests that spinal manipulation may induce
task dependent changes and also help to elucidate the exact mechanisms of how spinal
manipulation can impact human motor control.

5. Limitation

There are multiple limitations to consider for this study. The force production patterns were chosen in order to "isolate" lower threshold MUs that were believed to be impacted by spinal manipulation [19,55-57]. However, it is important to keep in mind that high-threshold motor units are also very likely to change following spinal manipulation [10] which the current study design in unable to elucidate. The current study design only allowed us to extract the lower threshold MU activity patterns that were executed at low force levels. In addition, note that the non-significant change in instantaneous MU firing rate gradient between "ramp" and "ramp & maintain" condition are only trends as no significant differences in the alterations from spinal manipulation were seen between the two patterns. This may be because both conditions were the same type of contraction, isometric, and executed at low force. The lack of any significant difference in the results between the different conditions may, therefore, be that they were too similar. Each individual was also adjusted according to his or her spinal dysfunction pattern, and this may have influenced the homogeneity of the results and, thus, the outcome. Each subject was adjusted at different spinal levels so the degree to which their particular spinal dysfunction pattern, and thus the spinal manipulations they received for it, was impacting the motor control of the tibialis anterior may have varied. In addition, we have to take into account that the spinal manipulation provided in this study was only a single session. Some neurophysiological changes observed after spinal manipulation have been shown to require multiples sessions over time [6].

Gender may have influenced the outcome of this study, although gender has been
shown to influence motor control more in an elderly population [79] than the participants
of the current study. This could be followed up in future studies, as the current sample size
does not allow for such factorial analysis.

6. Conclusions

This study shows that spinal manipulation provides at least short-term benefits in low
intensity movement, via the improvement of the recruitment of low threshold motor unit
and low twitch torque that are better suited for low force precise tasks. Further research
is required to investigate the longer term and potential functional effects of spinal manip-
ulation and chiropractic care in a variety of patients who would benefit from improved
muscle function. Future studies should also use higher force levels to investigate the effects
of spinal manipulation on motor unit behavior as EMG decomposition techniques have
improved and now also allow for this type of analysis.

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

**Funding:** Slovenian Research Agency (projects J2-1731 and L7-9421 and Programme funding P2-0041).

**Institutional Review Board Statement:** The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the North Denmark Region Committee on Health Research Ethics (N-20140027).

**Informed Consent Statement:** Informed consent was obtained from all subjects involved in the study.

**Data Availability Statement:** The data are not publicly available because of ethics committee restriction.

**Acknowledgments:** We would like to acknowledge NCMIC Foundation (Iowa, USA) and French government grant to University of Valenciennes for supporting the lead author Lucien with his PhD scholarship.

**Conflicts of Interest:** The authors declare no conflict of interest.

## Appendix A

<figure><img src="image_13.png" /><figcaption>**Figure A1.** Unmagnified distribution plot of the IFRG after intervention (control vs. spinal manipulation) for “ramp & Maintain” condition adjusted to baseline (Unmagnified Figure 8).</figcaption></figure>

## References

1. Jenkins, H.J.; Downie, A.S.; Moore, C.S.; French, S.D. Current evidence for spinal X-ray use in the chiropractic profession: A narrative review. *Chiropr. Man. Therap.* **2018**, *26*, 48. [CrossRef] [PubMed]

2. McGill, S.M. Stability: From biomechanical concept to chiropractic practice. *J. Can. Chiropr. Assoc.* **1999**, *43*, 75–88.
3. Christiansen, T.L.; Niazi, I.K.; Holt, K.; Nedergaard, R.W.; Duehr, J.; Allen, K.; Marshall, P.; Türker, K.S.; Hartvigsen, J.; Haavik, H. The effects of a single session of spinal manipulation on strength and cortical drive in athletes. *Eur. J. Appl. Physiol. 2018*, *118*, 737–749. [CrossRef] [PubMed]

4. 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. [CrossRef]

5. 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. [CrossRef]

6. Holt, K.R.; Haavik, H.; Lee, A.C.L.; Murphy, B.; Elley, C.R. Effectiveness of Chiropractic Care to Improve Sensorimotor Function Associated With Falls Risk in Older People: A Randomized Controlled Trial. *J. Manip. Physiol. Ther. 2016*, *39*, 267–278. [CrossRef] [PubMed]

7. Holt, K.; Niazi, I.K.; Nedergaard, R.W.; Duehr, J.; Amjad, I.; Shafique, M.; Anwar, M.N.; Ndetan, H.; Türker, 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. [CrossRef]

8. Daligadu, J.; Haavik, H.; Yielder, P.C.; Baarbe, J.; Murphy, B. Alterations in Cortical and Cerebellar Motor Processing in Subclinical Neck Pain Patients Following Spinal Manipulation. *J. Manip. Physiol. Ther. 2013*, *36*, 527–537. [CrossRef]

9. Haavik, H.; Murphy, B. Transient modulation of intracortical inhibition following spinal manipulation. *Chiropr. J. Aust.* **2007**, *37*, 106–116.

10. Haavik, H.; Niazi, I.; Jochumsen, M.; Sherwin, D.; Flavel, S.; Türker, K. Impact of Spinal Manipulation on Cortical Drive to Upper and Lower Limb Muscles. *Brain Sci.* **2016**, *7*, 2. [CrossRef]

11. 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. [CrossRef] [PubMed]

12. Mauguiere, F. Somatosensory evoked potentials: Normal responses, abnormal waveforms, and clinical applications in neurological disease. In *Electroencephalography. Basic Principles, Clinical Applications, and Related Fields*; Lippincott Williams & Wilkins: Philadelphia, PA, USA, 2005; pp. 1067–1119.

13. Taylor, H.H.; Murphy, B. Altered Central Integration of Dual Somatosensory Input After Cervical Spine Manipulation. *J. Manip. Physiol. Ther. 2010*, *33*, 178–188. [CrossRef] [PubMed]

14. Taylor, H.H.; Murphy, B. Altered Sensorimotor Integration with Cervical Spine Manipulation. *J. Manip. Physiol. Ther. 2008*, *31*, 115–126. [CrossRef]

15. Mendonça, G.V.; Teodósio, C.; Mouro, M.; Freitas, S.; Mil-Homens, P.; Pezarat-Correia, P.; Vila-Chã, C. Improving the Reliability of V-Wave Responses in the Soleus Muscle. *J. Clin. Neurophysiol. 2019*, *36*, 97–103. [CrossRef] [PubMed]

16. Mills, K. *Oxford Textbook of Clinical Neurophysiology*; Oxford University Press: Oxford, UK, 2017.

17. 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. [CrossRef] [PubMed]

18. 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. [CrossRef] [PubMed]

19. 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. [CrossRef]

20. Holobar, A.; Glaser, V.; Gallego, J.A.; Dideriksen, J.L.; Farina, D. Non-invasive characterization of motor unit behaviour in pathological tremor. *J. Neural Eng.* **2012**, *9*, 056011. [CrossRef]

21. Radeke, J.; van Dijk, J.P.; Holobar, A.; Lapatki, B.G. Electrophysiological method to examine muscle fiber architecture in the upper lip in cleft-lip patients. *J. Orofac. Orthop.* **2014**, *75*, 51–61. [CrossRef]

22. Perry, J.; Easterday, C.S.; Antonelli, D.J. Surface Versus Intramuscular Electrodes for Electromyography of Superficial and Deep Muscles. *Phys. Ther. 1981*, *61*, 7–15. [CrossRef]

23. Farina, D.; Holobar, A.; Merletti, R.; Enoka, R.M. Decoding the neural drive to muscles from the surface electromyogram. *Clin. Neurophysiol.* **2010**, *121*, 1616–1623. [CrossRef] [PubMed]

24. Gazzoni, M.; Farina, D.; Merletti, R. A new method for the extraction and classification of single motor unit action potentials from surface EMG signals. *J. Neurosci. Methods 2004*, *136*, 165–177. [CrossRef]

25. Hogrel, J.-Y. Use of surface EMG for studying motor unit recruitment during isometric linear force ramp. *J. Electromyogr. Kinesiol.* **2003**, *13*, 417–423. [CrossRef]

26. Holobar, A.; Zazula, D. Multichannel Blind Source Separation Using Convolution Kernel Compensation. *IEEE Trans. Signal Process.* **2007**, *55*, 4487–4496. [CrossRef]

27. Merletti, R. *Surface Electromyography*, 1st ed.; Wiley-IEEE Press: Piscataway, NJ, USA, 2016.

28. Kleine, B.U.; van Dijk, J.P.; Lapatki, B.G.; Zwarts, M.J.; Stegeman, D.F. Using two-dimensional spatial information in decomposition of surface EMG signals. *J. Electromyogr. Kinesiol.* **2007**, *17*, 535–548. [CrossRef]

29. Nawab, S.H.; Chang, S.-S.; De Luca, C.J. High-yield decomposition of surface EMG signals. *Clin. Neurophysiol.* **2010**, *121*, 1602–1615. [CrossRef] [PubMed]

30. Vila-Chã, C.; Falla, D.; Correia, M.V.; Farina, D. Adjustments in Motor Unit Properties during Fatiguing Contractions after Training. *Med. Sci. Sports Exerc.* **2012**, *44*, 616–624. [CrossRef]
31. Duchateau, J.; Semmler, J.; Enoka, R. Training adaptations in the behavior of human motor units. *J. Appl. Physiol.* **2006**, *101*, 1766-1775. [CrossRef]

32. Vila-Chã, C.; Falla, D.; Farina, D. Motor unit behavior during submaximal contractions following six weeks of either endurance or strength training. *J. Appl. Physiol.* **2010**, *109*, 1455-1466. [CrossRef]

33. Deschenes, M. Short review: Rate coding and motor unit recruitment patterns. *J. Appl. Sport Sci. Res.* **1989**, *3*, 33-39. [CrossRef]

34. Merletti, R.; Holobar, A.; Farina, D. Analysis of motor units with high-density surface electromyography. *J. Electromyogr. Kinesiol.* **2008**, *18*, 879-890. [CrossRef] [PubMed]

35. Andreassen, S.; Arendt-Nielsen, L. Muscle fibre conduction velocity in motor units of the human anterior tibial muscle: A new size principle parameter. *J. Physiol.* **1987**, *391*, 561-571. [CrossRef] [PubMed]

36. Andrew, D.; Yielder, P.; Haavik, H.; Murphy, B. The effects of subclinical neck pain on sensorimotor integration following a complex motor pursuit task. *Exp. Brain Res.* **2018**, *236*, 1-11. [CrossRef] [PubMed]

37. Baarbé, J.; Yielder, P.; Haavik, H.; Debison-Larabie, C.; Holmes, M.; Murphy, B. Enhanced cerebellar disinhibition when cervical manipulation precedes motor learning in patients with subclinical neck pain. *J. Chiropr. Med.* **2015**, *14*, 115-116.

38. Baarbé, J.K.; Holmes, M.W.R.; Murphy, H.E.; Haavik, H.; Murphy, B.A. Influence of Subclinical Neck Pain on the Ability to Perform a Mental Rotation Task: A 4-Week Longitudinal Study with a Healthy Control Group Comparison. *J. Manip. Physiol. Ther.* **2016**, *39*, 23-30. [CrossRef]

39. Baarbé, J.K.; Yielder, P.; Haavik, H.; Holmes, M.W.R.; Murphy, B.A. Subclinical recurrent neck pain and its treatment impacts motor training-induced plasticity of the cerebellum and motor cortex. *PLoS ONE* **2018**, *13*, e0193413. [CrossRef]

40. Farid, B.; Yielder, P.; Holmes, M.; Haavik, H.; Murphy, B.A. Association of Subclinical Neck Pain with Altered Multisensory Integration at Baseline and 4-Week Follow-up Relative to Asymptomatic Controls. *J. Manip. Physiol. Ther.* **2018**, *41*, 81-91. [CrossRef]

41. Lee, H.; Nicholson, L.L.; Adams, R.D.; Bae, S.-S. Proprioception and Rotation Range Sensitization Associated with Subclinical Neck Pain. *Spine* **2005**, *30*, E60-E67. [CrossRef]

42. Lee, H.-Y.; Wang, J.-D.; Yao, G.; Wang, S.-F. Association between cervicocephalic kinesthetic sensibility and frequency of subclinical neck pain. *Man. Ther.* **2008**, *13*, 419-425. [CrossRef]

43. Paulus, I.; Brumagne, S. Altered interpretation of neck proprioceptive signals in persons with subclinical recurrent neck pain. *J. Rehabil. Med.* **2008**, *40*, 426-432. [CrossRef]

44. 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. [CrossRef]

45. Hawker, G.A.; Mian, S.; Kendzerska, T.; French, M. Measures of adult pain: Visual Analog Scale for Pain (VAS Pain), Numeric Rating Scale for Pain (NRS Pain), McGill Pain Questionnaire (MPQ), Short-Form McGill Pain Questionnaire (SF-MPQ), Chronic Pain Grade Scale (CPGS), Short Form-36 Bodily Pain Scale (SF). *Arthritis Care Res.* **2011**, *63*, S240-S252. [CrossRef]

46. 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. Therap.* **2013**, *21*, 36. [CrossRef] [PubMed]

47. WHO. *WHO Guidelines on Basic Training and Safety in Chiropractic*; WHO: Geneva, Switzerland, 2005.

48. The Rubicon Group Definition and Position Statement on the Chiropractic Subluxation. 2017. Available online: http://www.therubicongroup.org/#/policies/ (accessed on 10 December 2020).

49. Holt, K.; Russell, D.; Cooperstein, R.; Young, M.; Sherson, M.; Haavik, H. Interexaminer reliability of a multidimensional battery of tests used to assess for vertebral subluxations. *Chiropr. J. Aust.* **2018**, *46*, 100-117.

50. Cooperstein, R.; Gleberzon, B. *Technique Systems in Chiropractic*; Churchill-Livingstone: New York, NY, USA, 2004.

51. Hessell, B.; Herzog, W.; Conway, P.; McEwen, M. Experimental measurement of the force exerted during spinal manipulation using the Thompson technique. *J. Manip. Physiol. Ther.* **1990**, *13*, 448-453.

52. Herzog, W.; Conway, P.; Zhang, Y.; Gál, J.; Guimaraes, A. Reflex responses associated with manipulative treatments on the thoracic spine: A pilot study. *J. Manip. Physiol. Ther.* **1995**, *18*, 233-236.

53. Herzog, W. Mechanical, physiologic, and neuromuscular considerations of chiropractic treatments. *Adv. Chiropr.* **1996**, *3*, 69-285.

54. Pickar, J.; Wheeler, J. Response of muscle proprioceptors to spinal manipulative-like loads in the anesthetized cat. *J. Manip. Physiol. Ther.* **2001**, *24*, 2-11. [CrossRef]

55. Masakado, Y.; Nielsen, J.B. Task-and Phase-related Changes in Cortico-muscular Coherence. *Keio J. Med.* **2008**, *57*, 50-56. [CrossRef]

56. Duchateau, J.; Enoka, R.M. Human motor unit recordings: Origins and insight into the integrated motor system. *Brain Res.* **2011**, *1409*, 42-61. [CrossRef]

57. Conway, B.A.; Reid, C.; Halliday, D.M. Low Frequency Cortico-Muscular Coherence During Voluntary Rapid Movements of the Wrist Joint. *Brain Topogr.* **2003**, *16*, 221-224. [CrossRef] [PubMed]

58. Farina, D.; Holobar, A. Characterization of Human Motor Units from Surface EMG Decomposition. *Proc. IEEE.* **2016**, *104*, 353-373. [CrossRef]

59. OTBioelettronica. *OT BioLab+*; OTBioelettronica: Turin, Italy, 2017.

60. Holobar, A. DEMUSE® Tool. 2020. Available online: https://demuse.feri.um.si/ (accessed on 10 March 2019).
61. Holobar, A.; Zazula, D. Correlation-based decomposition of surface electromyograms at low contraction forces. *Med. Biol. Eng. Comput.* **2004**, *42*, 487–495. [CrossRef] [PubMed]

62. Holobar, A.; Zazula, D. On the selection of the cost function for gradient-based decomposition of surface electromyograms. In *Proceedings of the 2008 30th Annual International Conference of the IEEE Engineering in Medicine and Biology Society*, Vancouver, BC, Canada, 20–25 August 2008; pp. 4668–4671. [CrossRef]

63. Holobar, A.; Farina, D.; Gazzoni, M.; Merletti, R.; Zazula, D. Estimating motor unit discharge patterns from high-density surface electromyogram. *Clin. Neurophysiol.* **2009**, *120*, 551–562. [CrossRef] [PubMed]

64. Minetto, M.A.; Holobar, A.; Botter, A.; Farina, D. Discharge Properties of Motor Units of the Abductor Hallucis Muscle during Cramp Contractions. *J. Neurophysiol.* **2009**, *102*, 1890–1901. [CrossRef]

65. Grassi, B.; Quaresima, V. Near-infrared spectroscopy and skeletal muscle oxidative function in vivo in health and disease: A review from an exercise physiology perspective. *J. Biomed. Opt.* **2016**, *21*, 091313. [CrossRef] [PubMed]

66. Gross, A.; Langevin, P.; Burnie, S.J.; Bédard-Brochu, M.S.; Empey, B.; Dugas, E.; Faber-Dobrescu, M.; Andres, C.; Graham, N.; Goldsmith, C.H.; et al. Manipulation and mobilisation for neck pain contrasted against an inactive control or another active treatment. In *Cochrane Database of Systematic Reviews*; Cochrane Back and Neck Group, Ed.; John Wiley & Sons: Hoboken, NJ, USA, 2015. [CrossRef]

67. Holobar, A.; Minetto, M.A.; Farina, D. Accurate identification of motor unit discharge patterns from high-density surface EMG and validation with a novel signal-based performance metric. *J. Neural Eng.* **2014**, *11*, 016008. [CrossRef]

68. Beck, R. Muscle fiber conduction velocity. In *Wiley Encyclopedia of Biomedical Engineering*; John Wiley & Sons: Hoboken, NJ, USA, 2006.

69. Farina, D.; Fortunato, E.; Merletti, R. Noninvasive estimation of motor unit conduction velocity distribution using linear electrode arrays. *IEEE Trans. Biomed. Eng.* **2000**, *47*, 380–388. [CrossRef]

70. Farina, D.; Muhammad, W.; Fortunato, E.; Meste, O.; Merletti, R.; Rix, H. Estimation of single motor unit conduction velocity from surface electromyogram signals detected with linear electrode arrays. *Med. Biol. Eng. Comput.* **2001**, *39*, 225–236. [CrossRef]

71. Bates, D.; Mächler, M.; Bolker, B.; Walker, S. Fitting Linear Mixed-Effects Models using lme4. *arXiv* **2014**, arXiv:1406.5823. Available online: http://arxiv.org/abs/1406.5823 (accessed on 5 March 2020).

72. Jochumsen, M.; Cremoux, S.; Robinault, L.; Lauber, J.; Arceo, J.C.; Navid, M.S.; Nedergaard, R.W.; Rashid, U.; Haavik, H.; Niazi, I.K. Investigation of Optimal Afferent Feedback Modality for Inducing Neural Plasticity with A Self-Paced Brain-Computer Interface. *Sensors* **2018**, *18*, 3761. [CrossRef] [PubMed]

73. Lenth, R.; Singmann, H.; Love, J. Emmeans: Estimated marginal means, aka least-squares means. *R Package Version* **2018**, *1*, 3.

74. Niazi, I.K.; Kamavuako, E.N.; Holt, K.; Janjua, T.A.M.; Kumari, N.; Amjad, I.; Haavik, H. The Effect of Spinal Manipulation on the Electrophysiological and Metabolic Properties of the Tibialis Anterior Muscle. *Healthcare* **2020**, *8*, 548. [CrossRef] [PubMed]

75. Farina, D.; Pozzo, M.; Merlo, E.; Bottin, A.; Merletti, R. Assessment of average muscle fiber conduction velocity from surface EMG signals during fatiguing dynamic contractions. *IEEE Trans. Biomed. Eng.* **2004**, *51*, 1383–1393. [CrossRef] [PubMed]

76. Kernell, D. The Limits of Firing Frequency in Cat Lumbosacral Motoneurones Possessing Different Time Course of Afterhyperpolarization. *Acta Physiol. Scand.* **1965**, *65*, 87–100. [CrossRef]

77. Kernell, D. Principles of Force Gradation in Skeletal Muscles. *Neural Plast.* **2003**, *10*, 69–76. [CrossRef]

78. De Luca, C.J.; Contessa, P. Biomechanical benefits of the onion-skin motor unit control scheme. *J. Biomech.* **2015**, *48*, 195–203. [CrossRef]

79. Gawel, M.; Kostera-Pruszczyk, A. Effect of Age and Gender on the Number of Motor Units in Healthy Subjects Estimated by the Multipoint Incremental MUNE Method. *J. Clin. Neurophysiol.* **2014**, *31*, 272–278. [CrossRef]