Skip to main content

Longitudinal estimation of intramuscular Tibialis Anterior coherence during subacute spinal cord injury: relationship with neurophysiological, functional and clinical outcome measures



Estimation of surface intramuscular coherence has been used to indirectly assess pyramidal tract activity following spinal cord injury (SCI), especially within the 15-30 Hz bandwidth. However, change in higher frequency (>40 Hz) muscle coherence during SCI has not been characterised. Thus, the objective of this study was to identify change of high and low frequency intramuscular Tibialis Anterior (TA) coherence during incomplete subacute SCI.


Fifteen healthy subjects and 22 subjects with motor incomplete SCI (American Spinal Injury Association Impairment Scale, AIS, C or D grade) were recruited and tested during 4 sessions performed at 2-week intervals up to 8 months after SCI. Intramuscular TA coherence estimation was calculated within the 10–60 Hz bandwidth during controlled maximal isometric and isokinetic foot dorsiflexion. Maximal voluntary dorsiflexion torque, gait function measured with the WISCI II scale, and TA motor evoked potentials (MEP) were recorded.


During subacute SCI, significant improvement in total lower limb manual muscle score, TA muscle strength and gait function were observed. No change in TA MEP amplitude was identified. Significant increase in TA coherence was detected in the 40–60 Hz, but not the 15–30 Hz bandwidth. The spasticity syndrome was associated with lower 15-30 Hz TA coherence during maximal isometric dorsiflexion and higher 10–60 Hz coherence during fast isokinetic movement (p < 0.05).


Longitudinal estimation of neurophysiological and clinical measures during subacute SCI suggest that estimation of TA muscle coherence during controlled movement provides indirect information regarding adaptive and maladaptive motor control mechanisms during neurorehabilitation.


Limited recovery up to 6–9 months following spinal cord injury (SCI) depends on both the severity and neurological level of the injury [1], especially after motor incomplete injury [1,2,3]. Currently, clinicians assume that optimal neurorehabilitation can be achieved during subacute incomplete SCI [3], including improvement in gait speed and walking distance [2, 3]. Neurophysiological measures of motor system function provide an objective basis to detect early recovery following SCI [4], although further characterisation of these measures are required, especially in relation to recovery of motor function [5]. Thus, there is a need to perform longitudinal studies of surrogate neurophysiological markers that can be measured quickly and which may be predictive of long-term follow up clinical outcomes. New surrogate measures with standard clinical outcome measures may provide indirect information regarding the development of adaptive neuroplasticity leading to the recovery of motor function [6].

Non-invasive neurophysiological measurement of motor evoked potentials (MEP) allows clinicians to obtain an objective evaluation of the integrity of motor pathways, such as the corticospinal tract, and has been used to predict the severity of SCI [7,8,9]. Longitudinal studies of corticospinal tract recovery during the subacute phase of SCI using this technique has previously identified a direct relationship between MEP amplitude and gait function [10, 11]. However, the reliability and validity of MEP parameters as biomarkers of recovery after SCI is not clear [12].

Electromyographic (EMG) muscle coherence is a mathematical index that estimates the degree of synchronization of activity from two independent muscle sites or signal sources, calculated in the frequency domain [13]. This non-invasive neurophysiological technique provides indirect information regarding the degree of common neural drive to motor neurones, and has been extensively used to assess pyramidal tract activity in healthy individuals [13,14,15,16,17,18,19]. Coherence estimation can be performed on the Tibialis Anterior (TA) [15], with corticospinal motor tract activity inferred from the beta band (10-40 Hz) after SCI, measured either at the low 10-20 Hz [20] or higher 24-40 Hz frequencies [21]. Furthermore, TA coherence estimation from two locations of the same muscle (also called intramuscular coherence estimation) after incomplete SCI subjects has been shown to reflect muscle strength, the severity of injury and the grade of spasticity symptoms [22].

More recently, coherence has been detected in the higher gamma (30-46 Hz) frequency band following SCI recorded from hand muscles in a non-human primate animal model, but without a change in the 15-30 Hz [23]. This fact suggests that higher frequency gamma band coherence could reflect activity mediated from either subcortical or propriospinal systems. However, beta and gamma band EMG coherence have also been related to static and dynamic motor tasks respectively [24]. Further characterisation of muscle coherence above 40 Hz has not been systematically characterised during subacute SCI, especially in relation to beta and standard gamma band activity.

Estimation of TA muscular coherence could help to indirectly identify adaptive neuroplasticity mechanisms related to top-down or spinal motor control, that have been demonstrated in animal models of SCI [25], and which may be present in humans [26]. The hypothesis of this study was that longitudinal estimation of a wider range of gamma-band coherence during the first months of incomplete SCI would provide information regarding the most appropriate test conditions required to highlight motor system neuroplasticity in relation to lower limb motor function and clinical outcome measures. Therefore, the objective of the present study was to monitor muscular TA coherence during controlled foot dorsiflexion in subjects during subacute incomplete SCI, in order to identify significant changes in the overall frequency content in the alpha (5-10 Hz), beta (10-40 Hz) and gamma (40-60 Hz) range, and to relate these changes to improvement in neurophysiological and clinical lower limb motor function. Preliminary findings of this study have been presented as an oral presentation at the International Congress on Neurorehabilitation [27].


Participants and ethical committee approval

This study was conducted at the “Hospital Nacional de Parapléjicos” in Toledo. Healthy subjects were recruited from the Castilla La Mancha region around Toledo. Ethical approval for the study was obtained from the Toledo Clinical Research Ethical Committee (#47, 07/05/2013) and all subjects were required to sign the informed consent before inclusion. The inclusion criteria for the SCI group included: age between 18 to 65 years, more than 1 month from the time of the SCI, injury severity graded between C-D according to the AIS score [28], neurological level between C2 y Th12 with a medically stable evolution and, a TA manual muscle strength score of >2 [29]. The exclusion criteria included: injury to the musculoskeletal or peripheral nervous system, neurological injury to supraspinal centres, epilepsy and pregnancy.

Experimental design

The experimental design was a descriptive prospective longitudinal design that included 4 repeated testing sessions performed at 2 week intervals during subacute SCI (1st-8th month after injury), which included tests for TA motor evoked potentials (MEP), muscle strength, gait function and lower limb spasticity syndrome. A group of non-injured subjects were recruited in one testing session to obtain reference data (main inclusion criteria, no neurological or musculoskeletal disorders and age between 18 and 65 years old.)

Experimental procedures

In addition to estimation of TA muscular coherence, maximal voluntary torque (MVT) of ankle dorsiflexion at 10° of plantarflexion, combined with Quadriceps (Q), Hamstrings (H), Triceps Surae (TS) and TA muscle scores performed to assess strength [29], and gait function using the walking index spinal cord injury scale (WISCI II) were measured [30]. Outcome measures related to the SCI spasticity syndrome [31] included muscle hypertonia with the modified Ashworth scale (MAS) [32], and the Penn scale [33]. SCI severity and neurological level was performed by a clinical consultant who had participated in the AIS training course with the European Multicenter Study about Spinal Cord Injury (EMSCI) network. Estimation of TA intramuscular coherence was performed by an engineer who was blinded to the experimental group allocations (MA) [22]. A trained physiotherapist (EBE) assessed neurophysiological outcomes (eg. MEP and TA EMG during controlled foot movement), and clinical measures (eg. muscle strength scores and lower limb spasticity tests).

Controlled foot dorsiflexion and Tibialis Anterior coherence estimation

The methods used to calculate TA muscular coherence have been published previously [22]. Subjects were seated in an isokinetic dynamometer (KinCom, Chattanooga Group Inc.). The trunk and pelvis were stabilized using straps and the hip, knee and ankle joints were flexed at 90°. A minimal TA muscle score of greater than 2 in the manual muscle score was set as the criteria for recruitment into the study, associated with a moderate muscle contraction level. Testing was performed on the most severely compromised lower limb following SCI, as long as the minimum muscle strength score was met or exceeded. Electromyographic activity was recorded using double differential surface electrodes, with a preamplifier gain of 10 V/V and an open filter bandwidth (Delsys Inc. Signal Conditioning Electrodes 3.1). Double differential electrodes were used to reduce electrical cross talk to a minimum between the recording electrodes. To record EMG intramuscular coherence, two surface electrodes were placed on two specific locations on a midline distal and proximal to the mid-point of the TA muscle belly, and separated by a distance of 10 cm to avoid electrical cross-talk [14, 17, 22, 34].

TA muscular coherence was estimated from EMG signals recorded during: i) two maximal isometric dorsiflexion contractions of the TA muscle which were maintained for 5 s, and separated by a 10 s rest period and, ii) ten isokinetic TA muscle contractions with the aim of recording a minimal signal of 3.5 s when combined together (see below), performed at 60°/s and 120°/s, with the ankle joint displaced from 30° plantarflexion to 0° dorsiflexion.

Tibialis Anterior motor evoked potentials

TA MEP were evoked with a transcranial magnetic stimulation device applied over the primary motor cortex (Magstim Rapid 2, Magstim Company Ltd) [35] with the subject seated on the dynamometer chair as described previously. Participants were provided visual biofeedback to help the subject reach 20% of the maximal dorsiflexion torque. A double-coned coil was oriented onto the cranium to induce a current stimulation in the posterior-to-anterior direction, with the optimal site located by varying the position of the coil from the vertex and increasing stimulus intensity until a MEP in the contralateral TA muscle was recorded at the lowest intensity. Ten MEP were recorded, with the threshold defined as the stimulus intensity at which five MEP of a minimal peak-to-peak amplitude of 100 mV [35].

Functional and clinical outcome measures

MVT during dorsiflexion was recorded using the isokinetic dynamometer. Individuals were explored in the supine position, including the: i) AIS scale [28]; ii) manual muscle score (0–5) of Q H, TS, TA and total muscle score (0–20 sum up of Q,H,G and TA scores); iii) muscle hypertonia detected during knee and ankle joint passive flexion and extension with the MAS [32]; iv) spasm frequency quantified with Penn scale [33].

Data analysis

TA electromyographic signals were sampled at a 10 KHz sampling frequency (MicroPlus 1401, Cambridge Electronic Design) and were subsequently down sampled to 2KHz using a low pass filter of 700 Hz to avoid aliasing (Matlab 7.11) [22]. TA intramuscular coherence estimation was performed with the Signal Processing Toolbox of Matlab 7.11 by estimating power spectral densities with Welch’s method [36]. Due to the methodological requirement of estimating coherence from a minimum length of EMG signal, several controlled movement tasks were performed so that when concatenated together a total muscle signal of at least 3.5 s could be analysed for each movement tasks. The signal was divided into 8 data segments using 50% overlapping segments performed with a Hamming window [37].

Coherence spectra were calculated against frequency for each subject for quality control and to identify line noise if present (eg. Fig. 1) and for all subjects (see Fig. 4). The coherence estimation for each frequency band (10–16, 15–30, 24–40 and 40-60 Hz) was calculated differently, with all coherence data points within each specific band for each movement averaged to obtain a final grand average coherence value (from 0 to 1). Inspection of individual TA EMG signals were not found to contain a marked level of line noise at 50 Hz, reducing the possibility that this frequency content could have contributed to false positive results associated with intramuscular coherence analysis (see Fig. 1). The velocity dependence of intramuscular TA coherence was calculated by calculating the ratio of its value during isokinetic TA muscle movement at 120°/s and at 60°/s.

Fig. 1
figure 1

Example Tibialis Anterior intramuscular coherence estimation made during maximal isometric foot dorsiflexion from a subject with SCI and presented for the 0-120 Hz frequency range. The two EMG measurements recorded simultaneously from the two recording sites from the Tibialis Anterior muscle are shown in Panels a and b. The Tibialis Anterior coherence estimation calculated between them is shown in Panel c. A peak in intramuscular coherence can be observed at 30 Hz. Notice that no marked line noise is recorded at 50 Hz. Surface EMG signals were recorded from the midline of the TA muscle, with point 1 located 5 cm proximally (1) and 5 cm distally (2) to the midpoint of the muscle belly

Amplified TA MEP were captured using an analogue-digital converter (MicroPlus 1401, Cambridge Electronic Design, UK) for subsequent analysis using a signal processing software package (Signal version 2.14. Cambridge Electronic Design, UK). MEP latency and amplitude were calculated when 7 out of 10 MEP’s were found to be consistent and showed low variability to the mean, with the response averaged from 5 signals. MEP amplitude was calculated from the highest to the lowest multiphase peak. Patients were diagnosed with the spasticity syndrome if a modified Ashworth score > 1 and/or a Penn score of ≥1 was diagnosed in 3 of the 4 of repeated testing sessions [22, 31, 38].

Statistical analysis

Two commercial software packages (SPSS, version 17.0 and SigmaStat, version 3.1, Systat software, Inc) were used to perform statistical analysis. Kolmogorov Smirnov test confirmed the non-normally distribution of the data and therefore non-parametric statistics were applied. With the exception of spectra data, used for quality data checking, all other data presented here are expressed as median values with 25th and 75th percentiles. In addition, this was performed to facilitate comparison with previous studies made by the group for coherence estimation [22, 27]. The Friedman test was performed for longitudinal analysis within patient groups, followed by application of the Bonferroni post-test used to identify differences among the sessions. The Mann-Whitney test was used to compare TA intramuscular coherence values between the healthy control and SCI groups, and between subjects with and without the SCI spasticity syndrome. Spearman correlations were applied to determine relationships between TA coherence, MEP, neuromuscular parameters (AIS scale, MVT, muscle score and spasticity measures) and gait function (WISCI II). Statistical significance was set at α ≤ 0.05.


Subject characteristics

Fifteen healthy subjects (8 males) with a median age of 26 years (24–29, 25th - 75th percentiles value range) and 22 subjects with motor incomplete SCI (18 males) with a median age of 39 years (36–55) were recruited for the follow-up study. No difference was identified for gender.

Initially fourteen subjects with SCI were diagnosed with the spasticity syndrome (modified Ashworth score > 1 and/or a Penn score of ≥1) and 8 were diagnosed without spasticity. However, one of the subjects (from the spastic group) was discharged prematurely after the first session and other individual from the non spastic group failed to complete the fourth session due to secondary complications. Finally, 13 SCI subjects diagnosed with spasticity syndrome and 7 SCI subjects without spasticity completed the study. No difference was identified in general for age, gender, AIS score, and time after SCI for the groups with and without the SCI spasticity syndrome. Furthermore, no significant differences were identified for general clinical characteristics, manual muscle, MVT during dorsiflexion of the foot or WISCI II scores. As expected, the Penn spasm and MAS scores were higher in individuals with SCI spasticity when compared SCI subjects without the spasticity syndrome (p < 0.001, Table 1).

Table 1 Individual SCI characteristics, and lower limb muscle and spasticity scores measured from the lower limb with the lowest total muscle strength score. Subjects without (1–7) and with the SCI spasticity syndrome (8–20) were recruited into the study

Change in clinical and functional scores during subacute SCI

The median testing period during subacute SCI ranged from 13 weeks (7–15, p25-p75)) for the initial session up to 19 weeks (14–22) for the fourth session, reflecting a repeated testing interval of approximately 2 weeks. Manual muscle strength scores generally increased after first session, for subjects with and without SCI spasticity from session 1 to 4 (p < 0.001; Fig. 2a). In contrast, dorsiflexion MVT only showed a specific increase for all subjects with SCI from session 3 (114 Nm 89–136) to 4 (139 Nm, 102–173; p < 0.05; Fig. 2b). Gait function also increased during subacute SCI for subjects with or without the spasticity syndrome (Fig. 2c).

Fig. 2
figure 2

Clinical and functional measures of motor function in all subjects with SCI during the four testing sessions of the subacute phase, and subcategorised into individuals with and without SCI spasticity syndrome. a. Total muscle (Quadriceps, Hamstring, Tibialis Anterior and Triceps Surae) strength score (0–20) [29]. b. Maximum voluntary dorsiflexion torque. Dotted line represents the median level of torque recorded from the noninjured control group. c. Gait function scale (WISCI II) (0–20) [30]. *: p < 0.05 with respect to session 1. Data are presented as median values with 25th and 75th percentiles

Tibialis Anterior motor evoked potentials

TA MEP records were evoked in subjects with SCI during session 1 (Fig. 3a) and 4 (Fig. 3b). However, no change in TA MEP amplitude (Fig. 3c) or latency (Fig. 3d) was observed during subacute SCI, for individuals with or without spasticity syndrome.

Fig. 3
figure 3

Tibialis Anterior motor evoked potentials recorded during subacute SCI. a. Example averaged TA MEP from ten individual records from a subject with SCI performed during session 1, and (b.) session 4. c. Group averaged TA MEP amplitude and (d.) latency recorded in subjects with SCI (n = 20), and those diagnosed with (n = 13) and without the spasticity syndrome (n = 7). Dotted line corresponds to the median of the voluntary control group data. Data are presented as median values with 25th and 75th percentiles

Correlation analysis of TA MEP with clinical and functional motor scores failed to reveal significant relationships with muscle strength scores (Q, H, TA, TS or total score), dorsiflexion MVT or the WISCI II scale, for subjects with or without the spasticity syndrome. However, the SCI spasticity group were characterised by correlations between TA MEP amplitude and the Penn spasm scale (rho = 0.41; p = 0.03), and TA MEP latency with TS muscle (rho = 0.33; p = 0.01) and total MAS score (rho = 0.28; p = 0.04).

TA muscle coherence during controlled activation

Analysis of the entire 10-60 Hz TA coherence bandwidth, recorded during maximum isometric dorsiflexion in the non-injured and SCI groups, showed no marked difference during the first testing session for individual frequency (Fig. 4 – Session 1). In contrast, during the last session a significant increase in the 50-60 Hz bandwidth was evident for subjects with SCI (Fig. 4 – Session 4). At the group level this increase in TA muscle coherence during subacute SCI was evident as a significant increase in the 40-60 Hz bandwidth during MVT dorsiflexion (p = 0.03, Fig. 5d), but not within the other studied bandwidths either during isometric (Fig. 5a-c) or for isokinetic movement (10-60 Hz). Furthermore, a weak correlation was observed between the 40-60 Hz TA coherence bandwidth measured during MVT dorsiflexion and the TA muscle score (Rho = 0.22, p = 0.05).

Fig. 4
figure 4

Intramuscular TA muscle coherence spectra (10–70 Hz). Session 1. Mean coherence spectra (as defined in the methods section) was calculated during maximal isometric dorsiflexion in the SCI group (n = 20) compared to the non-injured group (n = 15) during the first, and Session 2 last testing session 4 (see data analysis section in the methods). Data is presented as mean and standard error

Fig. 5
figure 5

Intramuscular TA coherence estimated following incomplete SCI during the 4 repeated testing sessions. a. Median TA coherence analysed within the 10–16 Hz, b. 15–30 Hz, c. 24–40 Hz and d. 40–60 Hz frequency range during maximal isometric dorsiflexion from 20 subjects with SCI. Dotted line corresponds to the median non-injured group coherence value. *: p ≤ 0.05 with respect to session 1 and #: p ≤ 0.05 with respect to the non-injured control data. Further methodological information can be found in the analysis section of the methods. Data presented as median values with 25th and 75th percentiles

With respect to SCI severity, subjects diagnosed with an AIS D injury were characterised with higher TA coherence within the 24-40 Hz, but not 10–16 Hz, 15–30 Hz, and 40-60 Hz bandwidths, when compared to the AIS C group, tested during isometric MVT for the entire subacute period (0.16, 0.11–0.19 vs 0.11, 0.09–0.15 respectively, p = 0.03).

TA muscle coherence and SCI spasticity syndrome

During the last testing session of subacute SCI (Session 4) intramuscular TA coherence estimation at high isokinetic movement speeds (120°/s) was higher for subjects with the SCI spasticity syndrome for all the studied bandwidths (Fig. 6a, 40-60 Hz: p < 0.01; 24-40 Hz: p < 0.01; 15-30 Hz: p = 0.02 and 10-16 Hz: p = 0.04). Likewise, the ratio of TA intramuscular coherence measured during isokinetic movement at 120/60 °/s also revealed higher values for the SCI spasticity group (Fig. 6b). Correlation analysis also revealed a negative relationship between velocity-dependent TA 40-60 Hz muscle coherence and specific SCI spasticity symptoms, specifically with the total lower limb (rho = −0.63; p = 0.02) and knee flexion MAS score (rho = −0.70; p = 0.01).

Fig. 6
figure 6

Velocity-dependent intramuscular TA coherence estimated during the last testing session of subacute SCI, in subjects diagnosed with and without the spasticity syndrome. a. Intramuscular TA coherence was estimated during isokinetic dorsiflexion of the foot at 120°/s, and (b.) expressed as a ratio of 120/60 °/s movement in subjects with and without SCI spasticity syndrome. Ankle joint movement was set at the same angular displacement for all subjects (see methods). *: p ≤ 0.05


For the first time, this study shows that limited recovery of lower limb muscle strength and general walking function during the first few months after incomplete SCI is accompanied by a significant increase in high-frequency (>40 Hz) band intramuscular TA coherence, without an increase in either TA motor evoked potential amplitude or beta band 15-30 Hz intramuscular coherence. Increase in high-frequency (>40 Hz) TA muscle coherence without an increase in TA MEP amplitude suggests that subcortical motor control mechanisms could play a role in the early recovery of lower limb muscle function during subacute SCI. Furthermore, the development of SCI spasticity was associated with a lower level of 15-30 Hz intramuscular TA coherence and TA motor evoked potentials, which highlights the relationship between this motor syndrome and corticospinal tract dysfunction.

Beta and gamma intramuscular TA coherence after spinal cord injury

Limited recovery of lower limb function, such as gait, observed in animal models of SCI has been shown to reflect neuroplasticity and new connections between corticospinal, reticulospinal and propriospinal motor control mechanisms [25]. Indeed, recovery of motor function in humans following incomplete SCI most probably includes similar neuroplasticity mechanisms active throughout the motor neuroaxis [1, 39, 40]. In the current study subjects with SCI showed clinical motor improvement, in addition to an increase in intramuscular gamma band TA 40–60 Hz coherence, but not in the beta band, which has traditionally been associated with corticospinal activity as assessed with transcraneal magnetic stimulation [41]. Intramuscular coherence within the gamma band has been proposed to reflect reticulospinal and propriospinal residual activity after spinal injury and other pathologies [23, 41, 42], although the contribution of corticospinal activity to standard gamma band coherence cannot be excluded [43, 44]. Interestingly in a previous study we have shown that TA gamma band coherence in motor incomplete SCI patients correlates with severe spasticity symptoms, including spasms and hypertonia during isometric dorsiflexion [22]. Taken together, these studies suggest that increased activity within reticulospinal and/or propriospinal control mechanisms, without an increase in TA MEP (see below), could be mediated both adaptive and maladaptive neuroplasticity of subcortical motor control pathways after SCI [45].

In our previous study, estimation of TA beta band coherence (15-30 Hz) after SCI was shown to correlate with residual muscle strength and gait function [22]. However, in the present study, 15–30 Hz coherence was generally stable in individuals with incomplete SCI and showed no change during subacute SCI. One explanation could be the relatively high number of subjects with the SCI spasticity syndrome (14/22), which may have blunted adaptive neuroplasticity of the corticospinal motor control system [22]. Another explanation could be that estimation of coherence activity in the 24–40 Hz may better reflect corticospinal activity after incomplete SCI [21], especially as we have shown previously that this band correlates with dorsiflexion isometric torque and is higher in patients with less severe SCI (AIS D) [22]. Nevertheless, no change in 24-40 Hz was detected in the present study for our cohort during subacute SCI.

Motor evoked potentials elicited by transcranial magnetic stimulation permits an objective assessment of SCI severity [8, 9, 46, 47], although the reliability of this technique in smaller patient cohorts has been questioned, especially in relation to reliability [12]. In other larger studies MEP amplitude and latency have been shown to correlate with functional and clinical motor outcome measures in specific cohorts [7, 46]. In the present study, TA MEP amplitude was decreased respect to the non-injured group, measured during subacute SCI, suggesting that longer follow-up times may be required to identify correlations with intramuscular TA coherence [11].

SCI spasticity syndrome and intramuscular coherence after spinal cord injury

The complex pathophysiology of the SCI spasticity syndrome may involve concomitant changes in cortical, subcortical and spinal motor control pathways [48, 49]. In this study, spasticity in subjects with SCI was characterised by a general decrease in intramuscular beta band TA coherence (15–30 Hz) during the isometric dorsiflexion. This reduction in potentially cortical motor control activity, together with correlations between TA MEP and severity of spasticity symptoms may suggest that maladaptive neuroplasticity of the corticospinal pathway plays a role in the development of this motor disorder.

In general, no difference was found with either the clinical scores such as the AIS, manual muscle, maximal voluntary dorsiflexion torque and gait scores, or the neurophysiological (TA muscle coherence or MEP’s) measures when they were compared between patients with and without spasticity during subacute SCI. However when TA coherence was measured for all the frequency bands studied (10–60 Hz), a significant velocity-dependent increase was noted for the SCI spasticity group, suggesting a generalised excitability of activity within the neuroaxis [22] including subcortical motor control systems [45]. Furthermore gamma intramuscular TA coherence was also higher at faster isokinetic speeds, but also correlated negatively with muscle hypertonia, similar to our previous study [22]. Finally, because low frequency coherence has been associated with activity in spinal networks (10-16 Hz) [21], increased TA coherence detected during isometric activation in subjects with spasticity suggests that hyperexcitability of spinal motor control mechanisms also play a role during faster movements.


Due to the small cohort size (n = 20) and limited follow up to 6 weeks during subacute SCI, caution should be made when we conclude that the increase in TA muscular coherence within the high-frequency (>40 Hz) band is not due to a corresponding increase in TA MEP amplitude. Longer follow up would also have clarified the impact of the development of SCI spasticity on TA muscle coherence. Furthermore, the relatively small and short duration TA muscle activation produced by patients with SCI may have affected the size of MEP, suggesting that standardisation of the TA MEP should be made to avoid methodological bias.


This study demonstrates that periodic estimation of high frequency intramuscular gamma band TA coherence (>40 Hz) in subjects with subacute incomplete SCI, may provide additional diagnostic information regarding the development of either adaptive or maladaptive motor control neuroplasticity during lower limb motor recovery. More information is required regarding the involvement of subcortical mechanisms that could mediate the limited lower limb motor recovery at this stage of rehabilitation, particularly as no evidence for a significant increase in either TA motor evoked potentials or 15-30 Hz beta band coherence identified in this small patient cohort. An important set of experimental studies using animal models have identified several secondary mechanisms of neuroplasticity organised at the subcortical and intraspinal level that could mediate adaptive neuronal motor control mechanisms following incomplete SCI, which suggest that the present clinical findings could be mediated by similar mechanisms such as the reticulospinal and/or propiospinal tracts. At the clinical level, evidence-based diagnosis of motor control adaptive neuroplasticity could be achieved by periodically testing TA intramuscular coherence during a simple motor task, but further studies should be performed to assess the clinical impact of this technique at the individual level.



American Spinal Injury Association Impairment Scale


electromyographic activity




modified Ashworth scale


motor evoked potentials


Maximal voluntary torque




spinal cord injury


Tibialis Anterior


Triceps Surae


  1. Fawcett JW, Curt A, Steeves JD, et al. Guidelines for the conduct of clinical trials for spinal cord injury as developed by the ICCP panel: spontaneous recovery after spinal cord injury and statistical power needed for therapeutic clinical trials. Spinal Cord. 2007;45(3):190–205.

    Article  CAS  PubMed  Google Scholar 

  2. Dobkin B, Apple D, Barbeau H, et al. Weight-supported treadmill vs over-ground training for walking after acute incomplete SCI. Neurology. 2006;66(4):484–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Yang JF, Musselman KE. Training to achieve over ground walking after spinal cord injury: a review of who, what, when, and how. J Spinal Cord Med. 2012;35(5):293–304.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Xie J, Boakye M. Electrophysiological outcomes after spinal cord injury. Neurosurg Focus. 2008;25(5):E11.

    Article  PubMed  Google Scholar 

  5. Steeves JD, Lammertse D, Curt A, et al. Guidelines for the conduct of clinical trials for spinal cord injury (SCI) as developed by the ICCP panel: clinical trial outcome measures. Spinal Cord. 2007;45(3):206–21.

    Article  CAS  PubMed  Google Scholar 

  6. Curt A, Schwab ME, Dietz V. Providing the clinical basis for new interventional therapies: refined diagnosis and assessment of recovery after spinal cord injury. Spinal Cord. 2004;42(1):1–6.

    Article  CAS  PubMed  Google Scholar 

  7. Curt A, Van Hedel HJ, Klaus D, Dietz V, Group E-SS. Recovery from a spinal cord injury: significance of compensation, neural plasticity, and repair. J Neurotrauma. 2008;25(6):677–85.

    Article  PubMed  Google Scholar 

  8. Calancie B, Alexeeva N, Broton JG, Suys S, Hall A, Klose KJ. Distribution and latency of muscle responses to transcranial magnetic stimulation of motor cortex after spinal cord injury in humans. J Neurotrauma. 1999;16(1):49–67.

    Article  CAS  PubMed  Google Scholar 

  9. Kobayashi M, Pascual-Leone A. Transcranial magnetic stimulation in neurology. Lancet Neurol. 2003;2(3):145–56.

    Article  PubMed  Google Scholar 

  10. Thomas SL, Gorassini MA. Increases in corticospinal tract function by treadmill training after incomplete spinal cord injury. J Neurophysiol. 2005;94(4):2844–55.

    Article  PubMed  Google Scholar 

  11. Petersen JA, Spiess M, Curt A, Dietz V, Schubert M, Group E-SS. Spinal cord injury: one-year evolution of motor-evoked potentials and recovery of leg motor function in 255 patients. Neurorehabil Neural Repair. 2012;26(8):939–48.

    Article  PubMed  Google Scholar 

  12. Boakye M, Harkema S, Ellaway PH, Skelly AC. Quantitative testing in spinal cord injury: overview of reliability and predictive validity. J Neurosurg Spine. 2012;17(1 Suppl):141–50.

    Article  PubMed  Google Scholar 

  13. Farmer SF, Bremner FD, Halliday DM, Rosenberg JR, Stephens JA. The frequency content of common synaptic inputs to motoneurones studied during voluntary isometric contraction in man. J Physiol. 1993;470:127–55.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Halliday DM, Rosenberg JR, Amjad AM, Breeze P, Conway BA, Farmer SF. A framework for the analysis of mixed time series/point process data--theory and application to the study of physiological tremor, single motor unit discharges and electromyograms. Prog Biophys Mol Biol. 1995;64(2–3):237–78.

    Article  CAS  PubMed  Google Scholar 

  15. Davey NJ, Ellaway PH, Baker JR, Friedland CL. Rhythmicity associated with a high degree of short-term synchrony of motor unit discharge in man. Exp Physiol. 1993;78(5):649–61.

    Article  CAS  PubMed  Google Scholar 

  16. Mills KR, Schubert M. Short term synchronization of human motor units and their responses to transcranial magnetic stimulation. J Physiol. 1995;483(Pt 2):511–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Conway BA, Halliday DM, Farmer SF, et al. Synchronization between motor cortex and spinal motoneuronal pool during the performance of a maintained motor task in man. J Physiol. 1995;489(Pt 3):917–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Nielsen JB, Conway BA, Halliday DM, Perreault MC, Hultborn H. Organization of common synaptic drive to motoneurones during fictive locomotion in the spinal cat. J Physiol. 2005;569(Pt 1):291–304.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Rosenberg JR, Amjad AM, Breeze P, Brillinger DR, Halliday DM. The Fourier approach to the identification of functional coupling between neuronal spike trains. Prog Biophys Mol Biol. 1989;53(1):1–31.

    Article  CAS  PubMed  Google Scholar 

  20. Barthelemy D, Willerslev-Olsen M, Lundell H, et al. Impaired transmission in the corticospinal tract and gait disability in spinal cord injured persons. J Neurophysiol. 2010;104(2):1167–76.

    Article  PubMed  Google Scholar 

  21. Norton JA, Gorassini MA. Changes in cortically related intermuscular coherence accompanying improvements in locomotor skills in incomplete spinal cord injury. J Neurophysiol. 2006;95(4):2580–9.

    Article  PubMed  Google Scholar 

  22. Bravo-Esteban E, Taylor J, Aleixandre M, et al. Tibialis Anterior muscle coherence during controlled voluntary activation in patients with spinal cord injury: diagnostic potential for muscle strength, gait and spasticity. J Neuroeng Rehabil. 2014;11:23.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Nishimura Y, Morichika Y, Isa T. A subcortical oscillatory network contributes to recovery of hand dexterity after spinal cord injury. Brain. 2009;132(Pt 3):709–21.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Gwin JT, Ferris DP. Beta- and gamma-range human lower limb corticomuscular coherence. Front Hum Neurosci. 2012;6:258.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Filli L, Engmann AK, Zorner B, et al. Bridging the gap: a reticulo-propriospinal detour bypassing an incomplete spinal cord injury. J Neurosci. 2014;34(40):13399–410.

    Article  CAS  PubMed  Google Scholar 

  26. Villiger M, Grabher P, Hepp-Reymond MC, et al. Relationship between structural brainstem and brain plasticity and lower-limb training in spinal cord injury: a longitudinal pilot study. Front Hum Neurosci. 2015;9:254.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Bravo-Esteban E, Taylor J, Aleixandre M, et al. Longitudinal estimation of intramuscular Tibialis Anterior coherence during subacute spinal cord injury: relationship with neurophysiological, clinical and functional measures. Paper presented at: Proceedings of the 2nd International Conference on NeuroRehabilitation (ICNR2014); 24–26 June, 2014; Aalborg.

  28. Maynard FM Jr, Bracken MB, Creasey G, et al. International Standards for Neurological and Functional Classification of Spinal Cord Injury. American Spinal Injury Association. Spinal Cord. 1997;35(5):266–74.

    Article  PubMed  Google Scholar 

  29. Medical Research Council of the UK. Aids to the investigation of peripheral nerve injuries. London: Pendragon House; 1976.

    Google Scholar 

  30. Dittuno PL, Ditunno JF Jr. Walking index for spinal cord injury (WISCI II): scale revision. Spinal Cord. 2001;39(12):654–6.

    Article  CAS  PubMed  Google Scholar 

  31. Gomez-Soriano J, Bravo-Esteban E. Perez-Rizo E, et al. Spinal Cord: Abnormal cutaneous flexor reflex activity during controlled isometric plantarflexion in human spinal cord injury spasticity syndrome; 2016.

    Google Scholar 

  32. Bohannon RW, Smith MB. Interrater reliability of a modified Ashworth scale of muscle spasticity. Phys Ther. 1987;67(2):206–7.

    Article  CAS  PubMed  Google Scholar 

  33. Penn RD, Savoy SM, Corcos D, et al. Intrathecal baclofen for severe spinal spasticity. N Engl J Med. 1989;320(23):1517–21.

    Article  CAS  PubMed  Google Scholar 

  34. Halliday DM, Conway BA, Christensen LO, Hansen NL, Petersen NP, Nielsen JB. Functional coupling of motor units is modulated during walking in human subjects. J Neurophysiol. 2003;89(2):960–8.

    Article  CAS  PubMed  Google Scholar 

  35. Barker AT, Jalinous R, Freeston IL. Non-invasive magnetic stimulation of human motor cortex. Lancet. 1985;1(8437):1106–7.

    Article  CAS  PubMed  Google Scholar 

  36. Welch P. The use of fast Fourier transform for the estimation of power spectra: a method based on time averaging over short, modified periodograms. Audio and electroacoustics, IEEE Transactions. 1967;15:4.

    Article  Google Scholar 

  37. Harris FJ. On the use of windows for harmonic analysis with the discrete Fourier transform. Proceedings IEEE. 1978;1

  38. Gómez-Soriano J, Castellote JM, Pérez-Rizo E, Esclarin A, Taylor JS. Voluntary ankle flexor activity and adaptive coactivation gain is decreased by spasticity during subacute spinal cord injury. Exp Neurol. 2010;224(2):507–16.

    Article  PubMed  Google Scholar 

  39. McKay WB. Neurophysiological characterization of the 'new anatomy' and motor control that results from neurological injury or disease. Clin Neurol Neurosurg. 2012;114(5):447–54.

    Article  PubMed  Google Scholar 

  40. Oudega M, Perez MA. Corticospinal reorganization after spinal cord injury. J Physiol. 2012;590(16):3647–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Fisher KM, Zaaimi B, Williams TL, Baker SN, Baker MR. Beta-band intermuscular coherence: a novel biomarker of upper motor neuron dysfunction in motor neuron disease. Brain. 2012;135(Pt 9):2849–64.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Nishimura Y, Isa T. Cortical and subcortical compensatory mechanisms after spinal cord injury in monkeys. Exp Neurol. 2012;235(1):152–61.

    Article  PubMed  Google Scholar 

  43. Barthelemy D, Willerslev-Olsen M, Lundell H, Biering-Sorensen F, Nielsen JB. Assessment of transmission in specific descending pathways in relation to gait and balance following spinal cord injury. Prog Brain Res. 2015;218:79–101.

    Article  PubMed  Google Scholar 

  44. Willerslev-Olsen M, Petersen TH, Farmer SF, Nielsen JB. Gait training facilitates central drive to ankle dorsiflexors in children with cerebral palsy. Brain. 2015;138(Pt 3):589–603.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Kumru H, Vidal J, Kofler M, Benito J, Garcia A, Valls-Sole J. Exaggerated auditory startle responses in patients with spinal cord injury. J Neurol. 2008;255(5):703–9.

    Article  PubMed  Google Scholar 

  46. Curt A, Keck ME, Dietz V. Functional outcome following spinal cord injury: significance of motor-evoked potentials and ASIA scores. Arch Phys Med Rehabil. 1998;79(1):81–6.

    Article  CAS  PubMed  Google Scholar 

  47. Diehl P, Kliesch U, Dietz V, Curt A. Impaired facilitation of motor evoked potentials in incomplete spinal cord injury. J Neurol. 2006;253(1):51–7.

    Article  PubMed  Google Scholar 

  48. Adams MM, Hicks AL. Spasticity after spinal cord injury. Spinal Cord. 2005;43(10):577–86.

    Article  CAS  PubMed  Google Scholar 

  49. Nielsen JB, Crone C, Hultborn H. The spinal pathophysiology of spasticity--from a basic science point of view. Acta Physiol (Oxf). 2007;189(2):171–80.

    Article  CAS  Google Scholar 

Download references


We are grateful for the continued support of the medical staff of the “Servicio de Rehabilitación”, and for the participation of the physiotherapists at the “Hospital Nacional de Parapléjicos” in recruiting subjects with subacute spinal cord injury. We also thank all the participants who participated in this study.


This project was funded by the “Consejería de Sanidad de Castilla–La Mancha” (03019–00, EQ02015), “Fundación Mutua Madrileña 2013” and the “Consolider-Ingenio 2010 Hyper” from the “Ministerio de Ciencia e Innovación” (CSDOOC-09-61,313 to E.B.E.).

Availability of data and materials

The datasets for the current study are available from the corresponding author on reasonable request.

Authors’ contributions

EBE participated in collecting the subject data, conducting data and statistical analysis, and helped to draft the manuscript. MA developed the software required to analyse EMG activity and muscle coherence analyses. CSM, helped collect the data and with GAM and IGA helped organize the data and to draft the manuscript. DT and JLP helped to draft the manuscript and supervise the research. JT and JGS conceived the study, interpreted the data, drafted the manuscript and critically revised it. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

This manuscript does not contain any individual person’s data in any form.

Ethics approval and consent to participate

Ethical approval was obtained from the Toledo Clinical Research Ethical Committee (# 47, 07/05/2013) and all subjects signed an informed consent sheet before recruitment.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Julian Taylor.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Bravo-Esteban, E., Taylor, J., Aleixandre, M. et al. Longitudinal estimation of intramuscular Tibialis Anterior coherence during subacute spinal cord injury: relationship with neurophysiological, functional and clinical outcome measures. J NeuroEngineering Rehabil 14, 58 (2017).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: