Research article | Open | Published:
The short-term recovery of corticomotor responses in elbow flexors
BMC Neurosciencevolume 20, Article number: 9 (2019)
The recovery of neurophysiological parameters at various time intervals following fatiguing exercise has been investigated previously. However, the repetition of neuromuscular assessments during the recovery period may have interfered with the true corticomotor excitability responses. In this experiment, fatiguing contractions were combined with a single post-fatigue assessment at varying time points. Ten participants undertook 5 bouts of 60-s maximal voluntary contractions (MVC) of the elbow flexors, separated by 20 min. Before and after each 60-s fatiguing exercise (FAT), participants performed a series of 6-s contractions at 100, 75 and 50% of their MVC during which transcranial magnetic, transmastoid electrical and brachial plexus electrical stimuli were used to elicit motor evoked potentials (MEP), cervicomedullary motor evoked potentials (CMEP) and compound muscle action potentials (Mmax) in the biceps brachii muscle, respectively. Post-FAT measurements were randomly performed 0, 15, 30, 60, or 120 s after each FAT.
MVC force declined to 65.1 ± 13.1% of baseline following FAT and then recovered to 82.7 ± 10.2% after 60 s. The MEP·Mmax−1 ratio recorded at MVC increased to 151.1 ± 45.8% and then returned to baseline within 60 s. The supraspinal excitability (MEP·CMEP−1) measured at MVC increased to 198.2 ± 47.2% and fully recovered after 30 s. The duration of post-MEP silent period recorded at MVC elongated by 23.4 ± 10.6% during FAT (all P < 0.05) but fully recovered after 15 s.
The current study represents the first accurate description of the time course and pattern of recovery for supraspinal and spinal excitability and inhibition following a short maximal fatiguing exercise in upper limb.
Neuromuscular fatigue (NMF) refers to any progressive exercise-induced change in the central (proximal to the neuromuscular junction) and/or peripheral (at or distal to the neuromuscular junction) systems that reduce the force production capacity of a voluntary or electrically-induced evoked contraction . In order to uncover the influence of exercise on the etiology of NMF development, transcranial magnetic stimulation (TMS) of the cerebral cortex in combination with electrical stimulation of the peripheral motor nerve (PNS) have been used extensively. With these non-invasive techniques, investigators can evaluate the contribution of central and peripheral mechanisms to the development of NMF in the exercised muscles . The specific patterns of NMF development and recovery are dependent upon the nature of the task, the subjects, and the muscle group tested [3,4,5]. Additionally, the various factors that contribute to fatigue recover at different rates. For instance, the gradual return of maximal voluntary force output (MVC) to baseline (which would take upwards of 10 min) has been attributed to peripheral factors while central factors such as central voluntary activation and corticospinal excitability have been shown to recover much more rapidly [6,7,8].
In the upper limb, previous studies have assessed the recovery of neurophysiological parameters following fatiguing single limb exercise by repeating neuromuscular evaluations at various time points from 2 s to 20 min post exercise [6, 7, 9,10,11,12,13]. Each neuromuscular evaluation involved multiple maximal or submaximal voluntary contractions; therefore, the repetitive contractions used throughout the recovery period may have altered the recovery pattern of neurophysiological parameters. Indeed, brief upper limb contractions (e.g. 2–6 s contractions at 50% MVC) have been shown to modulate the corticospinal pathway excitability for as long as 15 s  to 390 s  following the contractions. Some studies have avoided the use of additional contractions during the recovery period by analyzing recovery in the relaxed muscle; however, this method does not elucidate the recovery of central processes contributing to the performance of voluntary contractions [16,17,18].
The efficacy of the corticomotor pathway in the transfer of central commands (from cerebral cortex to the active muscles) can be evaluated through the electromyographic (EMG) responses to TMS and PNS, known as motor evoked potentials (MEP), and compound muscle action potentials (M wave), respectively [11, 16]. To further differentiate whether alterations in corticospinal excitability are occurring at the supraspinal or spinal motor neuron level, transmastoid electrical stimulation (TMES) of the descending corticospinal tract is used to produce cervico-medullary motor evoked potentials (CMEP) which reflect the excitability of the motoneuron pool [16, 17, 19]. Increases in MEP and CMEP amplitude or area relative to the Mwave, are indicative of heightened excitability of the corticospinal and spinal motoneurons, while the MEP·CMEP−1 ratio can be calculated to evaluate supraspinal excitability . An increase in silent period (SP, i.e. a period of EMG silence after MEP) is an indication of corticospinal inhibition [20,21,22,23].
MEP amplitude and area as well as the duration of silent period have been shown to increase during fatiguing contractions in upper limb muscles [7,8,9, 11]. When recovery measurements are taken on a relaxed muscle, the initial MEP facilitation is followed by a long-lasting MEP depression [16, 24, 25]. However, this depression is masked if MEP is measured during brief contractions. During contractions, the increased MEP displays an initial rapid recovery (depression). A complete return to baseline however has been shown to occur as quickly as 15 s following 3-min of intermittent elbow contractions , or as slowly as 10 min following repeated 22-s elbow flexor MVC’s . Similar to MEP, CMEP depression during the recovery period is only seen when the muscle is relaxed, whereas a return to baseline occurs within 15 s when CMEP is measured during contractions [16, 18, 26]. The increased SP during fatiguing upper limb exercise also typically recovers within 15–30 s of exercise termination [7, 8, 25, 27], although the exact time course of recovery for this measurement has yet to be established. Again, the repetitive post-fatigue evaluations may have prevented definitive conclusions on the recovery of EMG parameters to be drawn from these studies.
Thus, the current study aimed to specifically examine the recovery time courses of neurophysiological responses to a fatiguing task in the elbow flexors without the potential influence of repetitive post-fatigue assessment. We hypothesize that in the absence of repetitive post-fatigue assessments, the responses elicited by corticospinal stimuli including MEP, CMEP and SP would return to baseline level within a few seconds whereas the recovery of the maximal force output and voluntary activation would take much longer. An accurate estimate of the recovery of these neuromuscular responses (particularly within the first 60 s of recovery) could provide guidelines for the timing of future neuromuscular evaluations.
Ten recreationally active male participants (28.7 ± 6.1 year, 178.5 ± 5.2 cm, 75.1 ± 6.8 kg) volunteered to participate in this study. All participants completed the Physical Activity Readiness Questionnaire-Plus form, a checklist of contraindications for TMS , and signed the informed consent form before participating in the study. All the participants were right handed based on Edinburgh Handedness Inventory . None of the subjects had a history of musculoskeletal, neurological disease or were taking medications. The study protocol was conducted in accordance with declaration of Helsinki and was approved by the ethics board at the University of Calgary (REB16-0697).
The study was divided into two separate sessions: a familiarization and a testing session. During the familiarization session, the experimental procedure was explained, the written informed consent was obtained and participants were familiarized with the experimental protocol. During the actual testing session, participants were equipped with the EMG and stimulating electrodes. They were then seated in front of the elbow flexion ergometer with the hip, right shoulder and elbow positioned at 90°. The forearm was attached to a force transducer (Model LC101-2K, Omegadyne Inc., Sunbury, OH) with a velcro wrist strap. The force and EMG signals were monitored on a computer screen placed directly in front of the subject.
Before initiation of the neuromuscular evaluations, the optimal intensity for the PNS was determined (see below). The participants then undertook a warm-up protocol involving three repetitions of elbow flexion at 10%, 30%, 50% and one contraction at 70% of the MVC recorded during the familiarization session. Each warm up contraction was 5 s long with 5 s of rest in between contractions. The warm-up was followed by two 5-s isometric MVCs, with 2 min of rest in between. The determination of TMS and TMES intensities were performed at 20% MVC. Following this, neuromuscular function assessments including a series of TMS, TMES and PNS were elicited every 2-s while the participants performed continuous 6-s elbow flexor contractions at 100%, 75% and 50% of their MVC (Figs. 1, 2). For each assessment, the 75 and 50% of MVC values were immediately calculated after the MVC and displayed on a computer screen using the data acquisition software (LabChart 8), with no rest between contractions at each intensity. The fatigue protocol consisted of a 60-s elbow flexion MVC (FAT). The FAT was repeated 5 times every 20 min. Before and after each FAT, one neuromuscular function assessment was performed. The time-delay for the post-fatigue measurement was 0, 15, 30, 60, or 120 s, and was selected randomly for each participant.
EMG activity of the biceps and triceps brachii muscles was recorded with pairs of self-adhesive surface electrodes (10-mm recording diameter; Meditrace 100, Covidien, Mansfield, MA) in a bipolar configuration with 20-mm inter-electrode distance. A reference electrode was placed on the lateral epicondyle. Prior to placing the electrodes, the area of skin was shaved, abraded with sandpaper, and cleansed with an isopropyl alcohol swab to decrease skin resistance. An inter-electrode impedance of < 5 kΩ was obtained prior to recording to ensure an adequate signal-to-noise ratio. Signals were converted from analog to digital at a sampling rate of 2000 Hz using a PowerLab data acquisition system (16/35, ADInstruments, Bella Vista, Australia) and an octal bio-amplifier (ML138, ADInstruments; common mode rejection ratio = 85 dB, gain = 500) with a bandpass filter (5– 500 Hz) (ADInstruments). Data was then analyzed offline using LabChart 8 software (ADInstruments).
Transcranial magnetic stimulation (TMS)
The MEP responses of the biceps and triceps brachii muscles were elicited using a Magstim 200 stimulator (Magstim Company, UK) with a 110-mm double-cone coil (maximum output of 1.4 T) to preferentially stimulate the left motor cortex. Subjects wore a lycra swim cap on which intersecting lines were drawn to identify the vertex using the distance from nasion to inion and from left to right tragus. Every centimetre was demarcated along the nasal-inion line from the vertex to 2 cm posterior to the vertex, as well as laterally to 3 cm over the left motor cortex (6 positions). At each of these points, a stimulus was delivered at 50% maximal stimulator output during voluntary contractions at 20% MVC . The TMS hotspot was defined as the site where the largest MEP amplitude was evoked (on top of a contraction at 20% of MVC). The position of the coil was drawn directly onto the swim cap and kept constant throughout the protocol. Participants performed brief contractions at 20% MVC whereby a superimposed TMS was delivered at 30%, 40%, 50%, 60%, 70% and 80% of maximum stimulator output in random order. Four contractions were performed at each stimulus intensity with a 15-s rest interval between contractions. The intensity which showed the highest superimposed twitch (SIT), with largest MEP for biceps brachii and smallest MEP for triceps brachii, was chosen as TMS intensity for the rest of the testing session . The group mean stimulation intensity was 56 ± 9% of maximum stimulator output. The MEP amplitude recorded at this TMS intensity could be differentiated from the background EMG, during 100% MVC contractions. This stimulation intensity was then used for the remainder of the experiment.
Transmastoid electrical stimulation (TMES)
CMEP responses were evoked by passing a high-voltage electrical current between surface electrodes placed over the skin on the left (cathode) and right (anode) of the mastoid processes (stimulator Model DS7AH; Digitimer, Welwyn Garden City, Hertfordshire, UK). The stimuli intensity (pulse duration: 100 μs; 400 V maximum) was adjusted to produce CMEP amplitude that matched the MEP amplitude during brief 20% of MVC contraction. In other words, the amplitude of CMEP was matched with the MEP amplitude obtained from the TMS stimulus-responses curve explained above. The group mean stimulation intensity was 218 ± 92 mA.
Peripheral nerve stimulation (PNS)
Brachial plexus electrical stimulation was used to determine the size of the maximal compound muscle action potential (Mmax) of the biceps and triceps brachii. The stimulating electrodes (Ag–AgCl discs, 20-mm diameter) were placed on the supraclavicular fossa (cathode) and on the acromion process (anode). High-voltage percutaneous electrical stimuli (stimulator Model DS7AH; Digitimer) were delivered to the motor axons at the Erb’s point. The stimulation intensity (200-μs pulse duration; 400 V maximum) was increased in incremental steps (20 mA) until a plateau in Mmax and maximal twitch force was achieved at rest. The intensity was then increased by an additional 30% to ensure supramaximal stimulation. This stimuli intensity was used for the remainder of the experimental session. The group mean stimulation intensity was 232 ± 108 mA.
The force output of the elbow flexors and the background root mean square EMG (rmsEMG) of the biceps brachii muscle were assessed for the first, middle and last 5-s intervals of FAT in the 5 experimental conditions.
The force and the background EMG [rmsEMG of the biceps and triceps brachii muscles] were quantified over 500 ms prior to each MEP to ensure that changes in corticospinal responses measured at similar target forces were not due to differences in muscle activity prior to the stimulus. The maximum rmsEMG values were normalized to the amplitude of Mmax recorded during the same contraction to produce rmsEMG·Mmax−1 ratio (rmsEMG100, rmsEMG75, rmsEMG50).
The peak-to-peak amplitude and area under the signals were measured for each MEP, CMEP and Mmax signal recorded from biceps and triceps brachii muscles. Since the amplitude and area of the signals showed the same results, the area has been reported. In addition, the CMEPs recorded from triceps were not distinguishable from background EMG for majority of the participants, therefore these values were not reported. Onset of MEP, CMEP and Mmax were defined as the point at which the voltage trace became tangential to baseline in either the positive or negative direction. Because the Mmax can change as a result of the level of voluntary activation and fatigue, responses to biceps and triceps MEP and biceps CMEP were normalized to the subsequent Mmax recorded during the same contraction to produce MEP·Mmax−1 ratios (MEP100, MEP75, MEP50) and CMEP·Mmax−1 ratios (CMEP100, CMEP75, CMEP50), respectively. The MEP·CMEP−1 ratios were calculated for biceps at different contraction intensities (MEP·CMEP 100 −1 , MEP·CMEP 75 −1 and MEP·CMEP 50 −1 ) to identify the changes at the cortical level. The duration (ms) of silent period was assessed for biceps MEPs as the interval from the stimulus artefact to the return of the continuous EMG by visual inspection during 100, 75 and 50% of MVC contractions (SP100, SP75, SP50). The triceps silent period was not reported because of inconsistency of the data.
One of the aims of the present study was to measure voluntary activation using TMS. The amplitude of SITs evoked by TMS during contractions at 100, 75, and 50% MVC were calculated and the y-intercept of the linear regression between the SITs was used to quantify the estimated resting twitch . Although the TMS-evoked SITs demonstrated a linear regression at rest (r2 > 0.9), the post-FAT regression lines were not linear for some participants. This issue could be attributed to poor motor control during fatigue state. Specifically, several participants were unable to maintain the post-FAT MVC force at a plateau level. In these cases, the TMS was elicited at submaximal force (~ at 95% of peak MVC) which would result in overestimation of SIT amplitudes recorded at MVC. Accordingly, these VA values were removed from the manuscript.
Statistical analyses were computed using SPSS software (version 23.0; SPSS, Inc., Chicago, IL). Assumption of normality (Shapiro–Wilk test) and sphericity (Mauchley test) were tested for all of the dependent variables. If the assumption of sphericity was violated, the corrected value for non-sphericity with Greenhouse–Geisser correction was reported. Firstly, in order to compare the rate of decline in the MVC force and maximum rmsEMG across the five FAT trials, a two-way repeated measure analysis of variances (ANOVA) was performed, i.e. time (first, middle and last 5-s of each 60-s contraction) × five experimental conditions. Secondly, in order to ensure that there was no significant difference between the five randomized experimental conditions at baseline, one-way repeated measures ANOVA was run on the raw data for all variables. To measure the potential cumulative effect of the 5 × 60-s sustained MVCs performed in the experimental session, an additional one-way ANOVA was conducted on the five baseline measures (independent of the randomized post-fatigue trials). This cumulative effect was measured for all outcome variables. In addition, coefficient of variation (CV) was calculated for the MVC force, corresponding rmsEMG, MEP, CMEP and SP values using the five baseline measures (Table 1). Thirdly, in order to determine the effect of FAT on each variable, a paired t-test was used to compare the absolute baseline and 0 s values. Finally, the values obtained from measurements at each recovery time point were normalized to the corresponding baseline measure performed before each FAT trial. A one-way repeated measures ANOVA was conducted on five normalized data points to determine the influence of elapsed time on the recovery of MVC force, TMS-induced SIT, rmsEMG, corticospinal excitability and SP measures at 100, 75 and 50% of MVCs contractions. If a significant effect was obtained from the ANOVAs, Bonferroni corrected paired sample t-tests were performed to compare different time points. Significance was defined as P = 0.05.
No significant differences were observed between the five baseline measures for all variables when the randomization was applied (P > 0.2). The individual data for the MVC force and rmsEMG∙Mmax−1 are depicted in Fig. 3a. Similarly, no cumulative effect was evident for all outcome measures (P > 0.1) except for MVC force which demonstrated a trend towards lower force output at baseline #5 compared to #1 (P = 0.058). The CVs for the force output, rmsEMG·Mmax−1, MEP·Mmax−1, CMEP·Mmax−1 and SP measured at baseline are presented in Table 1. The CVs were also calculated for the MVC force and biceps rmsEMG recorded at the beginning of five FAT trails and were found to be 14.9% and 16.6%, respectively.
Force and rmsEMG during FAT trials
Force and the corresponding rmsEMG measured during the sustained 60-s MVC declined to 54.8 ± 12.2% (group mean: 335.7 ± 88.7 to 172.5 ± 52.8 N; time effect: F2,18 = 78.21, P < 0.001) and 83.9 ± 26.7% (1.06 ± 0.27 to 0.84 ± 0.29 mV; time effect: F2,18 = 5.18, P = 0.018), respectively (Fig. 3c, d). There was no significant condition or interaction effects between the five trials.
Force and rmsEMG
The MVC force significantly declined from pre- to post-FAT measurements (to 65.1 ± 13.1% of baseline, P < 0.001). MVC showed a significant recovery at 60 s (to 82.4 ± 10.1%, P = 0.005) compared to the 0 s time point (Fig. 4a) despite remaining below baseline even at the 120 s measurement (F4,36 = 7.95, P < 0.001). The rmsEMG100 recorded from biceps (500 ms before MEPs) was trending towards a significant decrease from baseline to 0 s, and then to recovery at 15 s (P = 0.052, Fig. 4b). The rmsEMG50 significantly declined from baseline to 0 s (P = 0.001) however ANOVA did not show any significant recovery during the post-FAT time points. No change was observed for rmsEMG75.
Biceps MEP100 demonstrated a significant facilitation from baseline to 0 s (increased to 151.1 ± 45.8% of baseline, P = 0.011) and then gradually returned towards baseline (F1.7,14.3 = 6.24, P = 0.014) until a significant recovery was displayed at 60 s (P = 0.038) (Fig. 5). The MEP75 did not show any statistical significance from baseline to 0 s, nor across the recovery time points. The MEP50 significantly declined from baseline to 0 s (to 84.6 ± 17.3%, P = 0.012). Although the recovery pattern for this parameter showed a trend towards significance (F4,36 = 2.56, P = 0.060), the Bonferroni corrected paired t-tests did not show any difference between post-FAT time points (Fig. 5). Since MEP values were normalized to Mmax, it is important to report that Mmax area recorded at 100, 75 and 50% of MVC did not show any significance difference between the time delays. The triceps MEP100 showed a trend toward a significant decrease (F1.6,10.1 = 4.26, P = 0.052), however Bonferroni corrected paired t-tests did not show any difference between the recovery time points (Table 2).
CMEP100 did not show any difference from baseline to 0 s (P = 0.091), however this measure declined significantly from 0 s to 15 s (to 69.3 ± 22.1%, P = 0.022). The CMEP75 and CMEP50 did now show any significant change from baseline to 0 s, nor across the recovery period (Fig. 5).
The MEP·CMEP 100 −1 indicated a significant increase in supraspinal excitability from baseline to 0 s (to 198.2 ± 47.2%, P < 0.001). During the recovery period, ANOVA showed a significant main effect for this parameter (F4,24 = 6.508, P < 0.001) where it demonstrated the highest values at 15 s, and then showed a significant return at 30 s (P = 0.045) (Fig. 5) The MEP·CMEP 75 −1 and MEP·CMEP 50 −1 ratios did not show any difference between the six time points.
While FAT resulted in a significant elongation of SP100 (to 123.4 ± 10.6%, P < 0.001) at 0 s, there was only a trend toward significance for SP75 (to 106.6 ± 9.1%, P < 0.055) and SP50 (+110.1% ± 10.6, P < 0.052). The SP100 returned to baseline at 15 s (P < 0.001) (Fig. 6).
The present study is the first to investigate the recovery time course of corticomotor responses in the elbow flexors after a fatiguing contraction without the interference of repeatitive recovery assessments. The main results supported our research hypotheses that (i) the impaired MVC force output partially recovered at 60 s, (2) the facilitated corticospinal and supraspinal excitability (indicated by increased MEP100 and MEP·CMEP 100 −1 , respectively) were recovered at 60 s and 30 s, respectively, and (3) the corticospinal inhibition, indicated by prolonged SP100, returned to baseline at 15 s. These findings confirm that measures of NMF and corticospinal excitability and inhibition are transient and recover at different paces. Thus, in order to achieve a precise estimation of the etiology of NMF, neurophysiological measurements should be performed immediately after cessation of the fatiguing task. Another important result of the present study is that the recovery pattern of neurophysiological responses is dependent upon the intensity of the voluntary contraction performed during the assessment.
MVC force and rmsEMG
Repeated neuromuscular assessments during the recovery period could modulate corticomotor excitability responses; therefore, the present experiment involved FAT trials combined with a single post-fatigue assessment at different time delays. The baseline MVC force output and the corresponding rmsEMG (Fig. 3a, b), as well as the force and rmsEMG at the beginning and termination of FAT contractions (i.e. the first and last 5 s of FATs, Fig. 3c, d) demonstrated similar values between the five randomized trails. These results confirm that the subjects were adequately recovered before the randomized FAT contractions.
Force declined to ~ 65% of baseline at 0 s which is similar to previous studies involving 60-s elbow flexor MVCs [25, 31]. At 60-s post-FAT, maximal force displayed a significant amount of recovery from its minimum value observed at 0 s, but never fully recovered to baseline over the course of measurements. Previous upper limb protocols involving sustained or intermitted isometric contractions (from 90 s to 3 min in duration) have demonstrated that the recovery of maximal voluntary force may take greater than 10 min to return to baseline [6,7,8]. The present study was not designed to directly discern the recovery pattern of central and/or peripheral mechanisms of fatigue. Yet the changes in rmsEMG100 (Fig. 4b) suggest that alteration in the descending central motor drive might have played a significant role in the reduction and recovery of force output at post-FAT level. On the other hand, those studies which investigated peripheral fatigue mechanisms in upper limb muscles have demonstrated that muscles contractile properties could take greater than 5 min to completely recover following fatiguing isometric contractions [6, 7, 10], very likely explaining why MVC was still 20% below initial values 120 s after exercise cessation. Further research with the single post-fatigue assessment paradigm is required to investigate the time course of full recovery of MVC force as well as central and peripheral mechanisms in upper limb muscles.
The MEP100 and MEP·CMEP 100 −1 ratio was increased at 0 s, whereas no change was observed for the CMEP100 at this time point (Fig. 5). These results confirm that the excitability of motor cortical cells is significantly increased during maximal contraction at the termination of FAT [9, 11, 20]. During the next 15 s, the CMEP100 declined significantly, however the MEP100 and MEP·CMEP 100 −1 ratio remained above the baseline (Fig. 5). Gandevia et al.,  also found an increase in MEP·Mmax−1 and MEP·CMEP−1 ratio immediately [2–25 s] after a 120-s MVC of the elbow flexors. Contrary to the present study however, these authors measured the corticospinal excitability on relaxed muscles and found a decrease in the CMEP·Mmax−1 ratio immediately after cessation of the task. Similar to our study, Taylor et al. [11, 32] and McNeil et al.  recorded MEP during maximal contraction and found MEP facilitation during a 120-s elbow flexor MVC. They attributed the enhanced MEP after sustained MVC to (1) a greater net excitability of the corticomotor cells as well as a decrease in refractoriness of these neurons which may increase descending excitatory volleys and/or (2) a decreased disynaptic inhibition at the spinal level. These mechanisms have been suggested to be compensatory pathways adopted by corticomotor circuitries to minimize deterioration of maximal muscle performance by increasing net excitatory motor drive to the spinal motoneurone pool. Nonetheless, contrary to the present work in which the recovery assessment began immediately following fatiguing exercise, the previously mentioned investigations initiated their post-fatigue cortical stimulations 15 s , 30 s  and 45 s  after the end of the sustained contraction, whereby MEP did not show any significant difference compared to baseline. Considering that an increase in MEP area during the sustained maximal elbow flexion is a well-established concept [8, 9, 11, 32], the failure of prior investigations in observing a residual MEP facilitation at 15, 30 and 45 s supports our argument that starting the post-fatigue assessments with a minimum delay provides a more accurate estimation of the corticospinal modulations.
Contrary to MEP100 facilitation, the corticospinal excitability recorded during submaximal contractions (i.e. MEP50) demonstrated a significant depression at 0 s. Very few studies have recorded the post-fatigue corticospinal responses during various maximal and submaximal contractions (e.g. contractions at 100, 75 and 50% of MVC); however, in line with this finding, Goodall et al.,  recently demonstrated a fatigue-induced depression in corticospinal excitability measured during submaximal contraction (10% MVC) in lower limb muscles. The reason for the disparity in MEP responses between maximal and submaximal contraction following FAT is unclear, however given that CMEP100 and CMEP50 did not demonstrate a significant alteration at 0 s (compared to baseline), it could be postulated that mediation at the motor cortical level was the primary determinant of the excitability of the entire corticospinal pathway (as indicated by the size of MEP). Indeed, supraspinal excitability assessed during MVC (MEP·CMEP 100 −1 ratio) showed significant facilitation at 0 and 15 s while no change was observed at 50% of MVC. The rmsEMG50 significantly declined from baseline to 0 s, thus it is plausible that a decrease in the magnitude of excitatory cortical volleys could have contributed in MEP50 depression. However, a trend towards a lower cortical discharge rate during MVC (indicated by rmsEMG100) could have artificially increased MEP100,as it reduces the likelihood of stimulating the corticomotoneuronal cells during the refractory period. Nonetheless, our results indicate that the mechanisms contributing to the increase in MEP100 and/or decrease in MEP50 after FAT are transient and recover quickly.
The CMEP100 initially demonstrated a non-significant facilitation at 0 s and then a significant depression at 15 s. As mentioned earlier, previous investigations that measured CMEP at rest demonstrated a long-lasting depression for this measure following maximal contractions (from 5 to 120 s). As such, it has been suggested that the intrinsic property of the spinal motoneurones becomes less responsive with fatigue [16, 18, 26, 34]. However, it has also been shown that the depression is typically masked when CMEP is recorded during voluntary contraction . Accordingly, we suggest that the deteriorated responsiveness of the spinal motoneurones at 0 s could have been masked by descending excitatory cortical input to the motoneurone pool. Considering that this descending excitatory cortical volley would also be present during the MVC at 15 s, it is unclear why CMEP100 decreased at this time point. It is unlikely that this phenomenon was due to a decline in the magnitude of excitatory input to the motoneurone pool. Indeed, it has been shown that withdrawal of excitatory input (disfacilitation) during maximal contraction would increase the responsiveness of the motoneurones and the amplitude of CMEP [35, 36]. Thus, the decrease in CMEP100 at 15 s could be associated with inhibitory signalling pathways diverging to the motoneurones, such as through the activation of spinal inhibitory interneurons . As mentioned earlier, Taylor et al.  suggest that the reduction of spinal inhibition could be a contributing factor in facilitation of MEP at the end of sustained contractions. We speculate that the restoration of inhibition from 0 to 15 s could not only partially recover (decrease) the facilitated MEP100 (Fig. 5), but also depress CMEP100 at this time point.
The SP100 was prolonged following FAT, however it fully recovered within 15 s of rest. The SP lengthens with fatigue leading to the conclusion that an augmented level of inhibition is built up in the corticomotor circuitries [8, 20]. The increase in duration of SP has been demonstrated following both sustained [9, 11] and intermittent isometric elbow flexor MVCs . The time course of complete recovery for this parameter has been reported as early as 30 s , 15 s , and 10 s  following the cessation of the isometric and intermittent contractions. Therefore, our data supports previous work indicating that SP recovers rapidly, and is not directly associated with motor cortical excitability, which increases after fatigue [8, 11]. SP is mediated through increasing activity of inhibitory interneurons releasing the neurotransmitter gamma-aminobutyric acid (GABA), and subsequent GABAB receptors activation [38, 39]. However, its mechanisms appear to be distinct from short interval intracortical inhibition which is mediated by GABAA inhibitory systems [33, 40, 41]. The data in the present study demonstrated only a trend toward significance for elongation of SP during submaximal contractions (i.e. SP75 and SP50). Although these results suggest that near maximal contractions are required to observe a prolongation of the SP [11, 20], Goodall et al.,  found a significant elongation of SP at 10% of MVC. Therefore, further research is necessary to evaluate mechanisms of SP recovery at different contraction intensities following fatigue.
There are several limitations to this study. As with all fatigue studies, the findings should be interpreted by taking the task specificity of fatigue mechanisms into consideration. To avoid contraction-induced fatigue, or modulations in corticomotor responses, only one neuromuscular fatigue assessment was performed before and after each fatiguing contraction. Previous investigations demonstrated that voluntary motor cortical outputs recovered within intervals varying from 15 s to 4 min [6,7,8], therefore we chose to record the neurophysiological parameters up to 120 s post-test. However, the recovery intervals chosen to represent this timeframe (i.e. 0, 15, 30, 60 and 120 s) were not ideal for the tracking of the changes in corticospinal excitability. Ideally neuromuscular assessments should have been performed after very short rest periods (e.g. 2, 5 and 10 s) to better analyze the patterns of corticospinal recovery. However, adding more fatiguing contractions was not feasible for this project. Finally, further experiments are required that include a control condition in which neurophysiological responses are quantified following repetitive neuromuscular evaluation (i.e. without fatiguing contractions).
The present study is the first to describe the time course of neuromuscular recovery following a short and maximal fatiguing exercise without the interference of repetitive post-fatigue assessments. We suggest that post-fatigue assessments should be initiated immediately following task cessation because corticospinal excitation and inhibition recover substantially within 30 s of recovery. Although the reasons are not clear, the present study also shows that the intensity of contraction influences corticospinal excitability and inhibition. Whether or not these findings apply to whole-body, dynamic exercises must be confirmed by using an innovative ergometer recently developed in our laboratory .
cervico-medullary motor evoked potentials
60-s elbow flexion MVC
motor evoked potentials
maximal voluntary force output
peripheral motor nerve
root mean square EMG·Mmax−1 ratio
transmastoid electrical stimulation
transcranial magnetic stimulation
MacIntosh BR, Rassier DE. What is fatigue? Can J Appl Physiol. 2002;27(1):42–55.
Gandevia SC. Spinal and supraspinal factors in human muscle fatigue. Physiol Rev. 2001;81(4):1725–89.
Babault N, Desbrosses K, Fabre MS, Michaut A, Pousson M. Neuromuscular fatigue development during maximal concentric and isometric knee extensions. J Appl Physiol. 2006;100(3):780–5.
Neyroud D, Rüttimann J, Mannion AF, Millet GY, Maffiuletti NA, Kayser B, Place N. Comparison of neuromuscular adjustments associated with sustained isometric contractions of four different muscle groups. J Appl Physiol. 2013;114(10):1426–34.
Morel B, Clémençon M, Rota S, Millet GY, Bishop DJ, Brosseau O, et al. Contraction velocity influence the magnitude and etiology of neuromuscular fatigue during repeated maximal contractions. Scand J Med Sci Sports. 2015;25(5):e432–41.
Hunter SK, Todd G, Butler JE, Gandevia SC, Taylor JL. Recovery from supraspinal fatigue is slowed in old adults after fatiguing maximal isometric contractions. J Appl Physiol. 2008;105(4):1199–209.
Szubski C, Burtscher M, Löscher WN. Neuromuscular fatigue during sustained contractions performed in short-term hypoxia. Med Sci Sports Exerc. 2007;39(6):948–54.
Taylor JL, Allen GM, Butler JE, Gandevia SC. Supraspinal fatigue during intermittent maximal voluntary contractions of the human elbow flexors. J Appl Physiol. 2000;89(1):305–13.
McNeil CJ, Martin PG, Gandevia SC, Taylor JL. The response to paired motor cortical stimuli is abolished at a spinal level during human muscle fatigue. J Physiol. 2009;587(Pt 23):5601–12.
Gandevia SC, Allen GM, Butler JE, Taylor JL. Supraspinal factors in human muscle fatigue: evidence for suboptimal output from the motor cortex. J Physiol. 1996;490(Pt 2):529–36.
Taylor JL, Butler JE, Allen GM, Gandevia SC. Changes in motor cortical excitability during human muscle fatigue. J Physiol. 1996;490(Pt 2):519–28.
Taylor JL, Butler JE, Gandevia SC. Changes in muscle afferents, motoneurons and motor drive during muscle fatigue. Eur J Appl Physiol. 2000;83(2–3):106–15.
Taylor JL, Petersen NT, Butler JE, Gandevia S. Corticospinal transmission after voluntary contractions. Adv Exp Med Biol. 2002;508:435–41.
Nørgaard P, Nielsen JF, Andersen H. Post-exercise facilitation of compound muscle action potentials evoked by transcranial magnetic stimulation in healthy subjects. Exp Brain Res. 2000;132(4):517–22.
Aboodarda SJ, Copithorne DB, Pearcey GEP, Button DC, Power KE. Changes in supraspinal and spinal excitability of the biceps brachii following brief, non-fatiguing submaximal contractions of the elbow flexors in resistance-trained males. Neurosci Lett. 2015;607:66–71.
Gandevia SC, Petersen N, Butler JE, Taylor JL. Impaired response of human motoneurones to corticospinal stimulation after voluntary exercise. J Physiol. 1999;521(Pt 3):749–59.
Taylor JL, Petersen NT, Butler JE, Gandevia SC. Interaction of transcranial magnetic stimulation and electrical transmastoid stimulation in human subjects. J Physiol. 2002;541(Pt 3):949–58.
Giesebrecht S, Martin PG, Gandevia SC, Taylor JL. Altered corticospinal transmission to the hand after maximum voluntary efforts. Muscle Nerve. 2011;43(5):679–87.
Todd G, Taylor JL, Gandevia SC. Reproducible measurement of voluntary activation of human elbow flexors with motor cortical stimulation. J Appl Physiol. 2004;97(1):236–42.
Sacco P, Thickbroom GW, Thompson ML, Mastaglia FL. Changes in corticomotor excitation and inhibition during prolonged submaximal muscle contractions. Muscle Nerve. 1997;20(9):1158–66.
Wilson SA, Lockwood RJ, Thickbroom GW, Mastaglia FL. The muscle silent period following transcranial magnetic cortical stimulation. J Neurol Sci. 1993;114(2):216–22.
Yacyshyn AF, Woo EJ, Price MC, McNeil CJ. Motoneuron responsiveness to corticospinal tract stimulation during the silent period induced by transcranial magnetic stimulation. Exp Brain Res. 2016;234(12):3457–63.
Inghilleri M, Berardelli A, Cruccu G, Manfredi M. Silent period evoked by transcranial stimulation of the human cortex and cervicomedullary junction. J Physiol. 1993;466:521–34.
Brasil-Neto JP, Pascual-Leone A, Valls-Solé J, Cammarota A, Cohen LG, Hallett M. Postexercise depression of motor evoked potentials: a measure of central nervous system fatigue. Exp Brain Res. 1993;93(1):181–4.
Sacco P, Thickbroom GW, Byrnes ML, Mastaglia FL. Changes in corticomotor excitability after fatiguing muscle contractions. Muscle Nerve. 2000;23(12):1840–6.
Butler JE, Taylor JL, Gandevia SC. Responses of human motoneurons to corticospinal stimulation during maximal voluntary contractions and ischemia. J Neurosci. 2003;23(32):10224–30.
Hunter SK, Butler JE, Todd G, Gandevia SC, Taylor JL. Supraspinal fatigue does not explain the sex difference in muscle fatigue of maximal contractions. J Appl Physiol. 2006;101(4):1036–44.
Rossi S, Hallett M, Rossini PM, Pascual-Leone A. Screening questionnaire before TMS: an update. Clin Neurophysiol. 2011;122(8):1686.
Oldfield RC. The assessment and analysis of handedness: the Edinburgh inventory. Neuropsychologia. 1971;9(1):97–113.
Temesi J, Gruet M, Rupp T, Verges S, Millet GY. Resting and active motor thresholds versus stimulus-response curves to determine transcranial magnetic stimulation intensity in quadriceps femoris. J Neuroeng Rehabil. 2014;11:40.
Todd G, Taylor JL, Gandevia SC. Measurement of voluntary activation of fresh and fatigued human muscles using transcranial magnetic stimulation. J Physiol. 2003;551(Pt 2):661–71.
Taylor JL, Butler JE, Gandevia SC. Altered responses of human elbow flexors to peripheral-nerve and cortical stimulation during a sustained maximal voluntary contraction. Exp Brain Res. 1999;127(1):108–15.
Goodall S, Howatson G, Thomas K. Modulation of specific inhibitory networks in fatigued locomotor muscles of healthy males. Exp Brain Res. 2018;236(2):463–73.
Petersen NT, Taylor JL, Butler JE, Gandevia SC. Depression of activity in the corticospinal pathway during human motor behavior after strong voluntary contractions. J Neurosci. 2003;23(22):7974–80.
Gandevia SC, Butler JE, Taylor JL. Viewpoint: Fatigue mechanisms determining exercise performance: integrative physiology is systems physiology. J Appl Physiol (1985). 2008;104(5):1546.
Martin PG, Gandevia SC, Taylor JL. Output of human motoneuron pools to corticospinal inputs during voluntary contractions. J Neurophysiol. 2006;95(6):3512–8.
Crone C, Nielsen J. Central control of disynaptic reciprocal inhibition in humans. Acta Physiol Scand. 1994;152(4):351–63.
Chen R, Lozano AM, Ashby P. Mechanism of the silent period following transcranial magnetic stimulation. Evidence from epidural recordings. Exp Brain Res. 1999;128(4):539–42.
Thompson SM, Gähwiler BH. Comparison of the actions of baclofen at pre- and postsynaptic receptors in the rat hippocampus in vitro. J Physiol. 1992;451:329–45.
Benwell NM, Mastaglia FL, Thickbroom GW. Differential changes in long-interval intracortical inhibition and silent period duration during fatiguing hand exercise. Exp Brain Res. 2007;179(2):255–62.
Werhahn KJ, Kunesch E, Noachtar S, Benecke R, Classen J. Differential effects on motorcortical inhibition induced by blockade of GABA uptake in humans. J Physiol. 1999;517(Pt 2):591–7.
Doyle-Baker D, Temesi J, Medysky ME, Holash RJ, Millet GY. An innovative ergometer to measure neuromuscular fatigue immediately after cycling. Med Sci Sports Exerc. 2018;50(2):375–87.
Conception and design of research: S.J.A., S.F., and G.Y.M.; Experimental acquisition: S.J.A., S.F., and G.Y.M.; Data analysis; S.J.A., S.F, and K.C.; Interpretation of the results: S.J.A., S.F., K.C. and G.Y.M.; Drafting, edition and revision of the manuscript; S.J.A., S.F., K.C. and G.Y.M. All authors read and approved the final manuscript.
Human Performance Laboratory, Faculty of Kinesiology, University of Calgary, Calgary, Canada.
This research was partially funded by the Eyes High Postdoctoral Scholars and Program for Undergraduate Research Experience (PURE) awards, University of Calgary, Alberta, Canada.
The authors declare that they have no competing interests.
Availability of data and materials
Raw data are available to all interested researchers by direct request of the corresponding author (Dr. Guillaume Y. Millet, Human Performance Laboratory, Faculty of Kinesiology, University of Calgary, 2500 University Dr NW, Calgary, Alberta, Canada T2N 1N4. E-mail: email@example.com)
Consent for publication
Ethics approval and consent to participate
The study protocol was conducted in accordance with declaration of Helsinki and was approved by the ethics board at the University of Calgary (REB16-0697).
This research was partially funded by the Eyes High Postdoctoral Scholars and Program for Undergraduate Research Experience (PURE) awards, University of Calgary, Alberta, Canada.
All participants completed the Physical Activity Readiness Questionnaire-Plus form, a checklist of contraindications for TMS , and signed the informed consent form before participating in the study.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.