Research | Open | Published:
Resting and active motor thresholds versus stimulus–response curves to determine transcranial magnetic stimulation intensity in quadriceps femoris
Journal of NeuroEngineering and Rehabilitationvolume 11, Article number: 40 (2014)
Transcranial magnetic stimulation (TMS) is a widely-used investigative technique in motor cortical evaluation. Recently, there has been a surge in TMS studies evaluating lower-limb fatigue. TMS intensity of 120-130% resting motor threshold (RMT) and 120% active motor threshold (AMT) and TMS intensity determined using stimulus–response curves during muscular contraction have been used in these studies. With the expansion of fatigue research in locomotion, the quadriceps femoris is increasingly of interest. It is important to select a stimulus intensity appropriate to evaluate the variables, including voluntary activation, being measured in this functionally important muscle group. This study assessed whether selected quadriceps TMS stimulus intensity determined by frequently employed methods is similar between methods and muscles.
Stimulus intensity in vastus lateralis, rectus femoris and vastus medialis muscles was determined by RMT, AMT (i.e. during brief voluntary contractions at 10% maximal voluntary force, MVC) and maximal motor-evoked potential (MEP) amplitude from stimulus–response curves during brief voluntary contractions at 10, 20 and 50% MVC at different stimulus intensities.
Stimulus intensity determined from a 10% MVC stimulus–response curve and at 120 and 130% RMT was higher than stimulus intensity at 120% AMT (lowest) and from a 50% MVC stimulus–response curve (p < 0.05). Stimulus intensity from a 20% MVC stimulus–response curve was similar to 120% RMT and 50% MVC stimulus–response curve. Mean stimulus intensity for stimulus–response curves at 10, 20 and 50% MVC corresponded to approximately 135, 115 and 100% RMT and 180, 155 and 130% AMT, respectively. Selected stimulus intensity was similar between muscles for all methods (p > 0.05).
Similar optimal stimulus intensity and maximal MEP amplitudes at 20 and 50% MVC and the minimal risk of residual fatigue at 20% MVC suggest that a 20% MVC stimulus–response curve is appropriate for determining TMS stimulus intensity in the quadriceps femoris. The higher selected stimulus intensities at 120-130% RMT have the potential to cause increased coactivation and discomfort and the lower stimulus intensity at 120% AMT may underestimate evoked responses. One muscle may also act as a surrogate in determining optimal quadriceps femoris stimulation intensity.
Transcranial magnetic stimulation (TMS) is a safe non-invasive technique employed to investigate motor cortical function. A rapidly changing magnetic field is produced by a coil placed over the target area of the brain and this causes electromagnetic induction to generate an electrical current in the brain. When sufficiently strong, this electrical current causes direct and trans-synaptic depolarization, and stimulation, of the pyramidal tract axons.
Selection of suitable TMS intensity is an important concern for researchers and clinicians. While being non-invasive, stimulation of the brain may be uncomfortable, particularly at high stimulus intensities. Thus, reducing the number of stimuli necessary to determine stimulus intensity and selecting the minimum intensity necessary to appropriately measure the desired parameters is beneficial to both investigators and subjects. The latter point has been largely absent in the literature despite several studies finding either similar or contradictory results when two different stimulus intensities were employed [1–3]. The majority of recent research has been conducted on clinical populations, and thus, recommendations are generally directed towards investigations in clinical populations or for clinical purposes [4, 5]. International Federation of Clinical Neurophysiology (IFCN) practical guidelines  discuss different methods of determining cortical motor threshold in relaxed muscle (RMT, resting motor threshold) and subsequent implications for stimulus intensity. These practical guidelines state that optimal intensity for TMS should correspond to the transition from the rising slope to the flat portion of the sigmoid stimulus–response (stimulator intensity-elicited motor-evoked potential (MEP) amplitude) curve and that this optimal intensity corresponds approximately to 140% RMT or 170% cortical motor threshold determined during voluntary muscular contraction (AMT, active motor threshold) . Stimulus–response curves are not routinely used for diagnostic purposes despite providing a direct means to determine stimulus intensity to elicit maximal MEP responses. This type of method has recently been employed by several research groups in the applied exercise sciences [2, 3, 6–8] while several other studies have determined stimulus intensity from RMT or AMT [9–12]. It remains to be determined if commonly employed selection of TMS intensity as determined by RMT, AMT and stimulus–response curves in this applied field result in selection of similar TMS intensities. Furthermore, practical guidelines for TMS intensity determination are normally based on investigations in upper-limb muscles. Data from lower-limb muscles are limited despite the functional importance of the lower limbs, specifically in regards to locomotion.
Studies utilizing TMS to investigate fatigue or acute exercise interventions in lower-limb muscles have used various methods to determine stimulus intensity. The most common of these has been RMT (the lowest intensity necessary to elicit MEPs, usually of at least 0.05 mV in amplitude, in at least one half of a given number of stimuli in the relaxed muscle) [9, 11, 13–15]. Another common method is AMT (the lowest intensity necessary to elicit detectable MEPs or MEPs of a pre-determined amplitude in at least one half of a given number of stimuli during weak voluntary contraction) [10, 12, 16–18]. More recently, numerous studies have selected a stimulus intensity to evoke MEP responses of a certain size in the target muscle during voluntary contraction [2, 6–8, 19, 20]. Some studies are unclear about the intensity chosen for TMS  or whether intensity determination was performed with the muscle in the relaxed or contracted state . Other studies based stimulus intensity on the intensity chosen to stimulate another muscle group  or simply selected maximal stimulator output .
Each of these methods produces a unique set of concerns. Cortical excitability is intrinsically linked to voluntary contraction intensity. While cortical excitability is low at rest, it increases rapidly as contraction intensity increases from rest [25, 26]. Whether determination of stimulus intensity in relaxed muscle (as with RMT) is appropriate for conducting measures in contracting muscle is unknown. Similarly, it remains to be determined whether selecting stimulus intensity at a different contraction level than that employed during evaluation is appropriate.
An additional complexity when evaluating leg muscles (e.g. knee extensors, knee flexors, plantar flexors) is that, unlike the elbow flexors, there is not a single dominant muscle. Whether it is appropriate to use a single muscle as a surrogate for all muscles within a muscle group (e.g. rectus femoris [RF] for the quadriceps femoris) when determining stimulus intensity remains to be investigated, especially since muscles and muscle groups may respond differently to TMS. This is a pertinent issue given both the functional importance of the quadriceps femoris and its increasing prevalence in studies utilizing TMS in the evaluation of fatigue [2, 7–9, 27].
Fatigue of the quadriceps is increasingly being evaluated in both healthy and clinical populations. An important measure in fatigue evaluation is voluntary activation (VA) [28, 29]. Evaluation of cortical VA utilizes superimposed twitches (SIT) evoked by TMS delivered during moderate- to high-intensity voluntary contractions (i.e. ≥50% maximal voluntary force [MVC]) [20, 30, 31]. Evoked MEP responses at ~50% MVC are theoretically maximal due to the firing of almost all motoneurons and maximal corticospinal excitability [20, 25, 31]. Since a key component of VA is the requirement that the muscle is driven maximally, maximal MEP amplitude is believed to be essential to ensure that SIT, and by extension VA, is not underestimated. Recently, quadriceps femoris studies have begun using TMS-induced antagonist coactivation as a criterion in the selection of TMS intensity [6–9] since this may cause SIT underestimation, and thus underestimate the development of central fatigue.
A comparison of selected stimulus intensity between published studies is impossible due to the use of different methods and equipment and different study aims. Thus, the primary objective of this study was to compare different methods of determining TMS intensity for the purposes of fatigue evaluation in the quadriceps femoris on selected stimulus intensity. Because of the use of voluntary contractions ≥50% MVC to determine VA and because maximal MEP responses have been observed to occur during contractions of approximately 50% MVC, a stimulus–response curve at 50% MVC was used as a baseline for comparison with other methods (i.e. this method most closely resembles fatigue evaluation). By using the same stimulator, coil and stimulation site, this protocol permits the isolation of differences between methods of stimulus intensity determination. The secondary objective was to determine whether selected stimulus intensity is similar for each of the three superficial quadriceps muscles.
Eight healthy active men participated in this study (means ± standard deviation: age, 30 ± 8 years; height, 181 ± 5 cm; body mass, 73 ± 4 kg). Subjects were informed of the experimental protocol and all associated risks prior to giving written informed consent as part of a medical inclusion. All procedures conformed to the Declaration of Helsinki and were approved by the Comité de Protection des Personnes Sud-Est 1, France.
Each subject completed one familiarization session and one experimental session. During the familiarization session, subjects were introduced to all procedures conducted in the experimental session and repeated trials until they performed all tests consistently and as directed. The largest MVC from the familiarization session was used to calculate contraction intensities and the reproducibility of MVC between sessions was verified.
Force and electromyographic recordings
Knee extensor force was measured during voluntary and evoked contractions by a calibrated force transducer (Meiri F2732 200 daN, Celians, Montauban, France) with amplifier that was attached by a non-compliant strap to the right leg immediately proximal to the malleoli of the ankle joint. Subjects were seated upright in a custom-built chair with both hips and right knee at 90° of flexion. The force transducer was fixed to the chair such that force was measured in direct line to the applied force. Electromyographic (EMG) activity of the right knee extensors (RF, vastus lateralis [VL] and vastus medialis [VM]) and flexors (biceps femoris, [BF]) was recorded.
EMG activity was recorded with a pair of self-adhesive surface (10-mm recording diameter) electrodes (Meditrace 100, Covidien, Mansfield, USA) in bipolar configuration with a 30-mm interelectrode distance and the reference on the patella. Low impedance (<5 kΩ) between electrodes was obtained by shaving, gently abrading the skin with sandpaper and then cleaning it with isopropyl alcohol. Signals were analogue-to-digitally converted at a sampling rate of 2000 Hz by PowerLab system (16/30—ML880/P, ADInstruments, Bella Vista, Australia) and octal bio-amplifier (ML138, ADInstruments) with bandpass filter (5–500 Hz) and analyzed offline using Labchart 7 software (ADInstruments).
Femoral nerve stimulation
Single electrical stimuli of 1-ms duration were delivered via constant-current stimulator (DS7A, Digitimer, Welwyn Garden City, Hertfordshire, UK) to the right femoral nerve via a 30-mm diameter surface cathode in the femoral triangle (Meditrace 100, Covidien, Mansfield, USA) and 50 x 90 mm rectangular anode (Durastick Plus, DJO Global, Vista, USA) on the gluteus maximus. Single stimuli were delivered incrementally until plateaus in maximal M-wave (Mmax) and twitch amplitude were reached. Three supramaximal stimuli at 130% of the intensity to produce maximal Mmax and twitch responses (52 ± 9 mA) were delivered at rest.
Transcranial magnetic stimulation
Single-pulses (0.1-ms rise time; 1-ms duration) were manually delivered by TMS to elicit MEPs and twitches in the right knee extensors. The contralateral motor cortex was stimulated by a magnetic stimulator (Magstim 2002, The Magstim Company Ltd, Whitland, UK) with 110-mm double-cone coil (maximum output of 1.4 T) to induce a postero-anterior current. The coil was manually controlled by an experienced investigator throughout the protocol. Subjects wore a cervical collar during all TMS measures to stabilize the head and neck.
Determination of coil position
Subjects wore a latex swim cap on which lines were drawn between the preauricular points and from nasion to inion to identify the vertex. Every centimeter from 1 cm anterior to 3 cm posterior to the vertex was demarcated along the nasal-inion line and also to 2 cm over the left motor cortex. At each point a stimulus was delivered at 70% maximal stimulator output during brief voluntary contraction of the knee extensors at 10% MVC. Target force was displayed on a screen and subjects were provided with real-time visual feedback during all voluntary contractions throughout the protocol. The coil was positioned at the site evoking the largest VL (39.5 ± 19.2% Mmax), RF (75.9 ± 26.7% Mmax) and VM (45.0 ± 21.3% Mmax) MEP amplitudes and SIT with minimal BF MEP amplitude. This coil position was drawn directly onto the swim cap and used throughout the protocol. Coil position was also verified before the delivery of each stimulus.
Determination of stimulus intensity
Four methods of determining stimulus intensity were investigated in the following order: 1) RMT: Beginning at 30% of maximal stimulator output and increasing by 5% increments to 80%, subjects received 10 stimuli at each stimulus intensity with the knee extensors completely relaxed. Stimuli were delivered at 10-s intervals. 2) AMT/stimulus–response curve at 10% MVC: Subjects performed brief voluntary contractions (~2-3 s) of the knee extensors with TMS delivered 10 times at 20, 25, 30, 35 and then 40% of maximal stimulator output. Subjects then performed brief contractions with TMS delivered 4 consecutive times at each of the following randomly-ordered stimulus intensities: 50, 60, 70 and 80% of maximal stimulator output. All stimuli were delivered at 15-s intervals. 3) Stimulus–response curve at 20% MVC: Subjects performed brief contractions (~2-3 s) of the knee extensors with TMS delivered 4 consecutive times at each of the following randomly-ordered stimulus intensities: 20, 30, 40, 50, 60, 70 and 80% maximal stimulator output. Stimuli were delivered at 15-s intervals. 4) Stimulus–response curve at 50% MVC: Similar to the stimulus–response curve at 20% MVC except that stimuli were delivered at 20-s intervals. During voluntary contractions, TMS was always delivered once the subject had contracted to the appropriate force level and the force had stabilized  and 10 min rest was provided between each of the four methods.
Peak-to-peak MEP and Mmax amplitudes were measured offline for each individual response. Individual MEP and Mmax amplitudes were then averaged and MEP amplitudes were normalized to Mmax amplitudes evoked in relaxed muscle. Data collected from a similar group of subjects in our laboratory indicated Mmax amplitudes were similar at rest and at the contraction intensities employed in this study (i.e. up to 50% MVC, unpublished observations, 2012). RMT was determined as the lowest stimulus intensity producing at least 5 MEPs of at least 0.05 mV from 10 stimuli. RMT was also determined from 6 and 8 stimuli (minimum of 3 and 4 MEPs, respectively). Stimulus intensities of 120 and 130% RMT were determined for comparison with methods used in other lower-limb studies [9, 11, 13, 14, 27, 30, 33]. AMT was determined by visual identification of MEPs above background EMG from contractions at 10% MVC  and corresponded to the lowest stimulus intensity producing MEPs in at least half the contractions. Classically, fixed thresholds are used to determine the presence of a MEP (i.e. 0.2 mV at 10% MVC ); however, the large variability in background EMG activity for the three measured quadriceps muscles rendered this method impractical. AMT was also determined from 6 and 8 stimuli (minimum of 3 and 4 MEPs, respectively). The stimulus intensity of 120% AMT was determined for comparison because of its use in other lower-limb studies [10, 12]. Stimulus–response curves at 10, 20 and 50% MVC were used to determine stimulus intensity by identifying the minimum stimulus intensity to evoke maximal MEP amplitude (i.e. the lowest intensity resulting in an increase of less than 5% MEP amplitude at higher stimulus intensities). Individual MEPs from a typical stimulus–response curve at 20% MVC for one subject are presented in Figure 1. Antagonist MEP amplitude was examined to verify that this stimulus intensity did not elicit increased TMS-induced coactivation. For the 10% MVC stimulus–response curve, only the first 4 stimuli at 20, 30 and 40% maximal stimulator output were considered. Where a plateau was not reached, MEP amplitude at 80% maximal stimulator output was compared to the estimated maximal MEP amplitude from Boltzmann modeling (see next paragraph). If mean MEP amplitude was greater or equal to the maximal modeled MEP amplitude, 80% was accepted as being part of the plateau and selected as the appropriate stimulus intensity. Otherwise, a plateau was determined to not have occurred and the data was excluded from analyses.
MEP amplitude from stimulus–response curves were modeled with a Boltzmann sigmoidal function  using the equation:
where MEPmax is the estimated maximal MEP amplitude, S is the stimulus intensity, S 50 is the stimulus intensity required to produce a response equal to half MEPmax and k is the slope parameter (inversely proportional to maximal function steepness). To eliminate the effects of background EMG in the modeling process, an amplitude of 0 mV was assigned to all responses in which there was no discernible MEP.
Statistical analyses were performed with Statistica (version 8, Tulsa, USA). The Shapiro-Wilk test was used to verify data normality. One-way repeated measures analyses of variance (ANOVA) were used to evaluate the method of stimulus determination (120 and 130% RMT, 120% AMT and stimulus–response curves), any difference between muscles and the effect of contraction intensity on Boltzmann parameters. One-way repeated measures ANOVA were also used to compare AMT and RMT determined from 6, 8 and 10 stimuli. When ANOVA revealed significant interactions, the Newman-Keuls post-hoc test was used to identify differences. The effect size as determined from the ANOVA were calculated as ω2 to reduce potential bias associated with the small sample size . Dependent t-tests were used to compare Boltzmann and linear relationships for the coefficient of determination of MEP amplitude. Statistical significance was set at p < 0.05. All data are expressed as means ± standard deviation except Figure 2 where values are expressed as means ± standard error of the mean.
Selected stimulus intensity
One subject did not reach a plateau in MEP amplitude in RF with the 10% MVC stimulus–response curve and was thus excluded from all relevant analyses.
Neither AMT nor RMT were different whether determination occurred with the first 6, 8 or 10 responses at each stimulus intensity for any muscle (p > 0.05). Therefore all subsequent analyses were conducted based upon AMT and RMT determined from 10 stimuli at each stimulus intensity. Selected TMS intensity determined by RMT, AMT and stimulus–response curves are presented in Figure 2. Stimulus intensities determined from RMT (120 and 130%) and stimulus–response curves at 10% MVC were higher than the intensity determined by stimulus–response curve at 50% MVC (VL: F(5,35) = 8.54, p < 0.001, ω2 = 0.48; RF: F(5,30) = 8.13, p < 0.001, ω2 = 0.50; VM: F(5,35) = 7.69, p < 0.001, ω2 = 0.45). Stimulus intensity at 120% AMT was lower than stimulus intensity determined from stimulus–response curves at both 10 and 20% MVC (p < 0.05). Table 1 presents the selected stimulus intensities from the stimulus–response curves as a percentage of the stimulus intensity to elicit both RMT and AMT to contextualize the differences between these methods. There was also no difference in selected intensity between muscles for any method (RMT: F(2,14) = 2.62, p = 0.11, ω2 = 0.16; AMT: F(2,14) = 1.21, p = 0.33, ω2 = 0.02; 10% MVC: F(2,12) = 1.00, p = 0.40, ω2 = 0.00; 20% MVC: F(2,14) = 1.15, p = 0.35, ω2 = 0.02; 50% MVC: F(2,14) = 0.778, p = 0.48, ω2 = 0) nor difference in normalized MEP amplitude at the selected stimulus intensity between 10, 20 and 50% MVC stimulus–response curves (VL: F(2,14) = 3.23, p = 0.07, ω2 = 0.21; RF: F(2,12) = 2.48, p = 0.13, ω2 = 0.16; VM: F(2,14) = 2.81, p = 0.09, ω2 = 0.18) (Table 2,). Raw BF MEP amplitudes at the selected stimulus intensities were 0.51 ± 0.54, 0.53 ± 0.41 and 0.53 ± 0.41 mV for VL, 0.42 ± 0.47, 0.43 ± 0.41 and 0.54 ± 0.31 mV for RF and 0.40 ± 0.45, 0.45 ± 0.40 and 0.59 ± 0.42 mV for VM for 10, 20 and 50% MVC stimulus response curves, respectively. At the stimulus intensity selected by AMT determined during contractions at 10% MVC, raw BF MEP amplitudes were 0.30 ± 0.41, 0.28 ± 0.42 and 0.29 ± 0.42 mV for VL, RF and VM, respectively. A single stimulus–response curve at 50% MVC is presented in Figure 3.
In VL, RF and VM, Mmax amplitudes were 16.2 ± 4.1 mV, 7.4 ± 1.8 mV and 17.0 ± 6.7 mV, respectively. Central drive as indicated by RMS · Mmax-1 for VL (0.0046 ± 0.0014), RF (0.0039 ± 0.0007) and VM (0.0053 ± 0.0025) at 10% MVC and VL (0.0088 ± 0.0024), RF (0.0086 ± 0.0019) and VM (0.0100 ± 0.0039) at 20% MVC were similar (F(2,14) = 1.32, p = 0.30, ω2 = 0.05, and F(2,14) = 0.660, p = 0.53, ω2 = 0, respectively). At 50% MVC, RMS · Mmax-1 for RF (0.0376 ± 0.0160) was greater than for both VL (0.0237 ± 0.0094) and VM (0.0264 ± 0.0115) (F(2,14) = 8.36, p = 0.004, ω2 = 0.00).
Boltzmann sigmoidal curves
Boltzmann curves from a typical subject are presented in Figure 4. Boltzmann curves provided a significantly better fit for the relationship between MEP amplitude and stimulator intensity than a linear relationship for stimulus–response curves at 10, 20 and 50% MVC for all muscles (p < 0.05). As contraction intensity increased, S 50 decreased in all muscles (VL: F(2,14) = 33.1, p < 0.001, ω2 = 0.79; RF: F(2,14) = 55.6, p < 0.001, ω2 = 0.87; VM: F(2,14) = 32.5, p < 0.001, ω2 = 0.79). Few differences were observed in MEPmax · Mmax-1 (only RF lower at 10% MVC; VL: F(2,14) = 1.88, p = 0.19, ω2 = 0.09; RF: F(2,14) = 3.88, p = 0.046, ω2 = 0.25; VM: F(2,14) = 2.40, p = 0.13, ω2 = 0.14) and k (only VL lower at 10% MVC; VL: F(2,14) = 7.50, p = 0.006, ω2 = 0.43; RF: F(2,14) = 1.62, p = 0.23, ω2 = 0.07; VM: F(2,14) = 0.911, p = 0.42, ω2 = 0). Results from modeling the stimulus–response curve data with the Boltzmann equation are presented in Table 3.
The main findings of this study are that (i) commonly-used stimulus intensities based upon RMT and a stimulus–response curve at 10% MVC are higher than those when determined using stimulus–response curves at 20 and 50% MVC and AMT and (ii) selected stimulus intensity, as determined by all methods, is similar between the three quadriceps muscles investigated. Because a stimulus–response curve performed at 20% MVC resulted in selection of a similar stimulus intensity to a stimulus–response curve at 50% MVC and because a stimulus–response curve at 20% MVC has a lower risk of inducing fatigue with repeated submaximal contractions, the present study indicates that this method is suitable for determining optimal stimulus intensity.
Comparison of methods
Resting motor threshold
In evaluation of the lower limbs to investigate fatigue or the effect of an exercise intervention, RMT has often been used to determine stimulus intensity. Most frequently this has been at 120 [11, 13, 14, 33] and 130% RMT [9, 27, 30]. The present study found that the use of these RMT intensities results in selection of higher stimulus intensities than a stimulus–response curve at 50% MVC and that stimulus intensity at 130% RMT is significantly greater than that from a stimulus–response curve at 20% MVC. No studies in the lower limbs were found to employ the suggested IFCN equivalent of 140% RMT , an intensity higher than the intensity at the transition from the rising slope to the plateau of the stimulus–response curves in the present study (Table 1).
There are several concerns about using RMT to determine optimal stimulus intensity in fatigue studies. The most important is whether it is appropriate to determine stimulus intensity in the relaxed muscle when evaluation of TMS-related parameters is conducted during muscular contraction. The rapid increase in cortical excitability from rest to even very weak contraction [25, 26] and the differential results in MEP evolution evaluated after fatiguing contractions when assessed in relaxed (i.e. decreased MEP amplitude/area [37–39]) and contracting (i.e. no change or increased MEP amplitude/area [10, 40, 41]) muscle present conceptual difficulties. More practically, increased stimulus intensity is associated with greater subject discomfort and this is important when recruiting healthy subjects and critical when evaluating clinical populations. If RMT is used to select stimulus intensity, no more than 6 stimuli should be delivered at each stimulus intensity since more stimuli do not better identify RMT contrary to the accepted standard of 10 stimuli at each intensity [5, 42]. It has also been reported that extremely high stimulus intensities are often required to determine RMT due to low cortical excitability at rest and that in some subjects RMT cannot be determined . This difficulty has also occurred in our laboratory. Given that high stimulus intensities may be required to evoke a MEP and the variable nature of MEP responses , particularly in the relaxed muscle , it may be difficult to identify an appropriate coil position.
Magnetic stimulation of the motor cortex with a double-cone coil permits more precise localization of specific brain areas than with a circular coil. It does not, however, permit localization with pin-point accuracy. Barker  detailed the induced electrical field and its rate of change with different coil types. Given that the motor cortex is not divided into discrete sections corresponding to individual muscles  and the imprecise area of stimulation with TMS, other muscle groups will inevitably be stimulated. Awiszus et al.  discussed the problem of high-intensity electrical muscle stimulation stimulating both agonist and antagonist muscles of the upper limb and these findings can likely be applied to transcranial motor cortical stimulation. To our knowledge, Todd et al.  were the first to specifically address this with a criterion in the determination of stimulus intensity (i.e. “a small MEP” in the antagonist). Figure 3 illustrates the 50% MVC stimulus–response curve of one subject. A plateau in quadriceps MEP amplitude corresponds to increased BF MEP amplitude and decreased SIT. In this subject, 120% RMT equated to 72, 66 and 78% maximal stimulator output in VL, RF and VM, respectively. This indicates that in some subjects, at 120% RMT, coactivation becomes apparent. Coactivation is problematic in the study of fatigue since quantification of cortical VA is essential. At higher stimulus intensities, such as those derived in the relaxed muscle from RMT, SIT during voluntary contraction may be underestimated because of increased contribution of antagonistic muscles  without corresponding augmentation of quadriceps femoris MEP amplitude.
Active motor threshold
Selected stimulus intensity at 120% AMT is significantly lower than stimulus–response curves at 10 and 20% MVC. All lower-limb studies employing AMT as a basis for TMS intensity determination utilized intensities much lower than the IFCN comparison equivalent of 170% AMT , recommendations much closer to a 10% MVC stimulus–response curve in this study (Table 1). As with RMT, the use of 6, 8 or 10 stimuli at each stimulus intensity when determining AMT did not affect the stimulus intensity selected.
Background EMG activity varies between individuals and also between muscles at a given contraction intensity; in some cases normal peak-to-peak amplitudes vary by >500% between subjects in the same muscle. Thus, the appropriateness of the common use of a fixed MEP amplitude to determine the presence of a MEP in evaluating AMT at different contraction intensities and in different subjects and/or muscles must be investigated. Boltzmann modeling indicates high inter-subject variability in evolution from no evoked MEP response to a maximal one (i.e. k; see Table 3). Some subjects demonstrated what could be characterized as a threshold from which no response immediately became a maximal one while in other subjects MEP amplitude gradually increased to maximum as stimulus intensity increased. Comparison with stimulus–response curves indicates that using AMT to determine stimulus intensity may result in submaximal MEP responses that are situated on the rising part of the Boltzmann curve. Unlike the use of maximal muscular responses to neural stimulation allowing serial or between-subject comparisons, comparison of submaximal evoked responses may introduce additional confounding factors. It remains to be established whether submaximal and maximal MEP responses and their evolution (e.g. with fatigue) are similar, particularly since preliminary indications from upper-  and lower-limb  studies suggest this may not always be the case. The evaluation of cortical VA may also be affected by the use of stimulus intensities derived from AMT (e.g. 120%). Stimulus intensity at 120% AMT was non-significantly lower than that determined from a 50% MVC stimulus–response curve and this might result in delivery of TMS at a submaximal intensity during contractions between 50 and 100% MVC and result in underestimated SIT. The effect on estimated resting twitch, calculated from the linear regression of three SITs from three different contraction intensities in this range and acceptable if r > 0.9 [48, 49], and subsequent estimation of cortical VA are unknown.
Generally, AMT is evaluated in voluntary contractions at 5 or 10% MVC in the upper limbs [50–53]. In lower limbs, Kalmar and Cafarelli  and Hilty et al.  used 3% MVC and found higher AMT than in Weier et al.  and the present study, the latter two having employed contractions at 10% MVC. This is consistent with Boltzmann modeling showing decreased stimulus intensity to evoke a MEP of half maximal amplitude (i.e. S 50; see Table 3) as contraction intensity increases.
All stimulus–response curves demonstrated a Boltzmann sigmoidal relationship, thus permitting the use of a stimulus–response curve to identify maximal MEP amplitudes and directly determine optimal diagnostic TMS stimulus intensity . Modeling of data indicated that estimated maximal MEP amplitude was lower at 10% MVC than at other contraction intensities although this was only significant in RF. Stimulus intensity to evoke a MEP of half maximal amplitude also decreased as contraction intensity increased. Determining stimulus intensity during contractions at 50% MVC would appear to be appropriate since evoked MEP responses at this contraction intensity are theoretically maximal [20, 25, 31] and both this and higher contraction intensities are used to determine cortical VA. A concern, however, is that an extended series of such contractions may produce measurable effects of fatigue, and consequently, that residual effects of fatigue may be present during a subsequent protocol as reported in a recent study . The lack of difference between stimulus intensity as determined by stimulus–response curves at 20 and 50% MVC and the similar maximal MEP amplitudes as determined by Boltzmann modeling suggest that in the quadriceps femoris, a stimulus–response curve at 20% MVC is appropriate to determine TMS intensity when the aim is to evaluate fatigue-related parameters such as VA.
Comparison of muscles
Studies determining stimulus intensity during contractions have often used normalized MEP amplitude or area of a given size as criteria [7, 8, 19, 20]. For example, Sidhu et al.  selected an intensity that produced the largest RF MEP with the stipulations that this must be at least 50% Mmax and that antagonist BF MEP amplitude be less than 10% raw RF MEP amplitude. In VL and VM in the present study, only 2 of 8 and 3 of 8 subjects, respectively, satisfied the requirement that MEP amplitude be ≥50% Mmax. In the case where several quadriceps muscles are examined, the latter criterion is ambiguous. BF may often be greater than 10% raw MEP amplitude in at least one knee extensor muscle and for one subject in the present study this was the case in all muscles at almost all stimulus intensities evaluated and also at almost all coil positions evaluated in the determination of optimal coil position.
There was no difference in selected stimulus intensity between muscles as determined by any method. This suggests that one muscle could be used as a surrogate for other quadriceps muscles. RF alone has frequently been used to determine quadriceps stimulus intensity [6, 8, 19, 20]. When RF is normalized, MEP amplitude is larger than for either VL or VM due to consistently smaller Mmax in the RF and little differences in raw MEP amplitude. The presentation of normalized RF MEP amplitudes instead of VL and VM may give the impression of eliciting greater corticospinal drive to the quadriceps muscles. In the present study, this was not due to a greater RF contribution since RMS · Mmax-1 was only greater than that of VL and VM at 50% MVC and normalized RF MEPs are larger than VL and VM at all contraction intensities. RF may not be an ideal surrogate because of the difficulty in recording clear M waves in this muscle. Furthermore, RF is the sole biarticular muscle of the quadriceps femoris, and thus, may not be representative of the muscle group.
An important limitation to the protocol is that it was not conducted on a second day to investigate the day-to-day variability of the methods employed in this study. Further investigations are required to establish whether the different methods employed to evaluate TMS parameters with fatigue are reproducible on different days. The present study also used a maximal response in contracting muscle as a reference point to evaluate multiple fatigue-related TMS parameters since this provides important insights into the manifestation and development of fatigue; however, recent studies suggest that TMS responses elicited by a submaximal stimulus intensity may also further understanding of fatigue mechanisms [1–3]. This reinforces the necessity of selecting an appropriate method to determine TMS intensity directly related to the parameters being investigated. In the context of the evaluation of cortical voluntary activation and corticospinal excitability with fatigue, maximal responses as investigated in this study are pertinent. In other research and diagnostic areas employing TMS, this may not be the case, and methods such as RMT and AMT may be the methods of choice for determining optimal stimulus intensity. Further studies must also determine the specific relevance of TMS-induced maximal and submaximal responses in both healthy and clinical populations in the context of fatigue, including the manner in which this affects measures of cortical voluntary activation and both excitatory and inhibitory mechanisms.
Percentages of AMT and RMT have often been employed to determine TMS intensity in studies evaluating fatigue; however, these methods do not accurately identify the minimum stimulus intensity to elicit MEPs of maximal amplitude in the quadriceps femoris. Thus, they may be inappropriate for cortical excitability and voluntary activation assessment. The potential for increased coactivation and discomfort at 120 and 130% RMT and possible underestimation of evoked responses at 120% AMT preclude their use. There are minor differences between selected stimulus intensity (lower at 50% MVC for VL only) from stimulus–response curves at 20 and 50% MVC. Both MEP amplitude at selected stimulus intensity and estimated maximal MEP amplitude determined from these stimulus–response curves are similar. This indicates that a stimulus–response curve performed at 20% MVC is a suitable method of determining TMS stimulus intensity while reducing the risk of inducing fatigue compared to methods at a higher percentage of MVC. From the present study, it is also concluded that determining stimulus intensity from a single muscle is acceptable in the quadriceps femoris.
Active motor threshold
Analysis of variance
International Federation of Clinical Neurophysiology
Slope parameter (Boltzmann modeling)
Maximal M-wave amplitude
Estimated maximal MEP amplitude (Boltzmann modeling)
Maximal voluntary force
Resting motor threshold
Stimulus intensity (Boltzmann modeling)
Stimulus intensity to produce a response half that of MEPmax (Boltzmann modeling)
Transcranial magnetic stimulation
McNeil CJ, Giesebrecht S, Gandevia SC, Taylor JL: Behaviour of the motoneurone pool in a fatiguing submaximal contraction. J Physiol 2011, 589: 3533-3544. 10.1113/jphysiol.2011.207191
Temesi J, Arnal PJ, Davranche K, Bonnefoy R, Levy P, Verges S, Millet GY: Does central fatigue explain reduced cycling after complete sleep deprivation? Med Sci Sports Exerc 2013, 45: 2243-2253. 10.1249/MSS.0b013e31829ce379
Temesi J, Rupp T, Martin V, Arnal PJ, Feasson L, Verges S, Millet GY: Central fatigue assessed by transcranial magnetic stimulation in ultratrail running. Med Sci Sports Exerc 2013. doi:10.1249/MSS.0000000000000207
Lefaucheur JP, Andre-Obadia N, Poulet E, Devanne H, Haffen E, Londero A, Cretin B, Leroi AM, Radtchenko A, Saba G, Thai-Van H, Litre CF, Vercueil L, Bouhassira D, Ayache SS, Farhat WH, Zouari HG, Mylius V, Nicolier M, Garcia-Larrea L: Recommandations franc¸aises sur l’utilisation de la stimulation magnétique transcrânienne répétitive (rTMS): règles de sécurité et indications thérapeutiques. Neurophysiol Clin 2011, 41: 221-295. 10.1016/j.neucli.2011.10.062
Groppa S, Oliviero A, Eisen A, Quartarone A, Cohen LG, Mall V, Kaelin-Lang A, Mima T, Rossi S, Thickbroom GW, Rossini PM, Ziemann U, Valls-Sole J, Siebner HR: A practical guide to diagnostic transcranial magnetic stimulation: report of an IFCN committee. Clin Neurophysiol 2012, 123: 858-882. 10.1016/j.clinph.2012.01.010
Rupp T, Jubeau M, Wuyam B, Perrey S, Levy P, Millet GY, Verges S: Time-dependent effect of acute hypoxia on corticospinal excitability in healthy humans. J Neurophysiol 2012, 108: 1270-1277. 10.1152/jn.01162.2011
Klass M, Roelands B, Levenez M, Fontenelle V, Pattyn N, Meeusen R, Duchateau J: Effects of noradrenaline and dopamine on supraspinal fatigue in well-trained men. Med Sci Sports Exerc 2012, 44: 2299-2308. 10.1249/MSS.0b013e318265f356
Sidhu SK, Bentley DJ, Carroll TJ: Locomotor exercise induces long-lasting impairments in the capacity of the human motor cortex to voluntarily activate knee extensor muscles. J Appl Physiol 2009, 106: 556-565.
Goodall S, Gonzalez-Alonso J, Ali L, Ross EZ, Romer LM: Supraspinal fatigue after normoxic and hypoxic exercise in humans. J Physiol 2012, 590: 2767-2782. 10.1113/jphysiol.2012.228890
Iguchi M, Shields RK: Cortical and segmental excitability during fatiguing contractions of the soleus muscle in humans. Clin Neurophysiol 2012, 123: 335-343. 10.1016/j.clinph.2011.06.031
Sammut R, Thickbroom GW, Wilson SA, Mastaglia FL: The origin of the soleus late response evoked by magnetic stimulation of human motor cortex. Electroencephalogr Clin Neurophysiol 1995, 97: 164-168. 10.1016/0924-980X(94)00324-Z
Weier AT, Pearce AJ, Kidgell DJ: Strength training reduces intracortical inhibition. Acta Physiol 2012, 206: 109-119. 10.1111/j.1748-1716.2012.02454.x
Lentz M, Nielsen JF: Post-exercise facilitation and depression of M wave and motor evoked potentials in healthy subjects. Clin Neurophysiol 2002, 113: 1092-1098. 10.1016/S1388-2457(02)00031-7
Tallent J, Goodall S, Hortobagyi T, St Clair Gibson A, French DN, Howatson G: Repeatability of corticospinal and spinal measures during lengthening and shortening contractions in the human tibialis anterior muscle. PLoS One 2012, 7: e35930. 10.1371/journal.pone.0035930
Mileva KN, Bowtell JL, Kossev AR: Effects of low-frequency whole-body vibration on motor-evoked potentials in healthy men. Exp Physiol 2009, 94: 103-116. 10.1113/expphysiol.2008.042689
Hilty L, Lutz K, Maurer K, Rodenkirch T, Spengler CM, Boutellier U, Jancke L, Amann M: Spinal opioid receptor-sensitive muscle afferents contribute to the fatigue-induced increase in intracortical inhibition in healthy humans. Exp Physiol 2011, 96: 505-517.
Kalmar JM, Cafarelli E: Central excitability does not limit postfatigue voluntary activation of quadriceps femoris. J Appl Physiol 2006, 100: 1757-1764. 10.1152/japplphysiol.01347.2005
Krishnan C, Dhaher Y: Corticospinal responses of quadriceps are abnormally coupled with hip adductors in chronic stroke survivors. Exp Neurol 2012, 233: 400-407. 10.1016/j.expneurol.2011.11.007
Fernandez-del-Olmo M, Rodriguez FA, Marquez G, Iglesias X, Marina M, Benitez A, Vallejo L, Acero RM: Isometric knee extensor fatigue following a Wingate test: peripheral and central mechanisms. Scand J Med Sci Sports 2013, 23: 57-65. 10.1111/j.1600-0838.2011.01355.x
Sidhu SK, Bentley DJ, Carroll TJ: Cortical voluntary activation of the human knee extensors can be reliably estimated using transcranial magnetic stimulation. Muscle Nerve 2009, 39: 186-196. 10.1002/mus.21064
McKay WB, Tuel SM, Sherwood AM, Stokic DS, Dimitrijevic MR: Focal depression of cortical excitability induced by fatiguing muscle contraction: a transcranial magnetic stimulation study. Exp Brain Res 1995, 105: 276-282. 10.1007/BF00240963
Hollge J, Kunkel M, Ziemann U, Tergau F, Geese R, Reimers CD: Central fatigue in sports and daily exercises. A magnetic stimulation study. Int J Sports Med 1997, 18: 614-617. 10.1055/s-2007-972691
Verin E, Ross E, Demoule A, Hopkinson N, Nickol A, Fauroux B, Moxham J, Similowski T, Polkey MI: Effects of exhaustive incremental treadmill exercise on diaphragm and quadriceps motor potentials evoked by transcranial magnetic stimulation. J Appl Physiol 2004, 96: 253-259.
Gibbons CE, Pietrosimone BG, Hart JM, Saliba SA, Ingersoll CD: Transcranial magnetic stimulation and volitional quadriceps activation. J Athl Train 2010, 45: 570-579. 10.4085/1062-6050-45.6.570
Taylor JL, Allen GM, Butler JE, Gandevia SC: Effect of contraction strength on responses in biceps brachii and adductor pollicis to transcranial magnetic stimulation. Exp Brain Res 1997, 117: 472-478. 10.1007/s002210050243
Ugawa Y, Terao Y, Hanajima R, Sakai K, Kanazawa I: Facilitatory effect of tonic voluntary contraction on responses to motor cortex stimulation. Electroencephalogr Clin Neurophysiol 1995, 97: 451-454. 10.1016/0924-980X(95)00214-6
Goodall S, Ross EZ, Romer LM: Effect of graded hypoxia on supraspinal contributions to fatigue with unilateral knee-extensor contractions. J Appl Physiol 2010, 109: 1842-1851. 10.1152/japplphysiol.00458.2010
Gandevia SC: Spinal and supraspinal factors in human muscle fatigue. Physiol Rev 2001, 81: 1725-1789.
Gruet M, Temesi J, Rupp T, Levy P, Millet GY, Verges S: Stimulation of the motor cortex and corticospinal tract to assess human muscle fatigue. Neuroscience 2013, 231: 384-399.
Goodall S, Romer LM, Ross EZ: Voluntary activation of human knee extensors measured using transcranial magnetic stimulation. Exp Physiol 2009, 94: 995-1004. 10.1113/expphysiol.2009.047902
Todd G, Taylor JL, Gandevia SC: Measurement of voluntary activation of fresh and fatigued human muscles using transcranial magnetic stimulation. J Physiol 2003, 551: 661-671. 10.1113/jphysiol.2003.044099
Gruet M, Temesi J, Rupp T, Millet GY, Verges S: Effect of different approaches to target force on transcranial magnetic stimulation responses. Muscle Nerve 2013, 48: 430-432. 10.1002/mus.23786
Ross EZ, Middleton N, Shave R, George K, Nowicky A: Corticomotor excitability contributes to neuromuscular fatigue following marathon running in man. Exp Physiol 2007, 92: 417-426. 10.1113/expphysiol.2006.035972
Sacco P, Hope PA, Thickbroom GW, Byrnes ML, Mastaglia FL: Corticomotor excitability and perception of effort during sustained exercise in the chronic fatigue syndrome. Clin Neurophysiol 1999, 110: 1883-1891. 10.1016/S1388-2457(99)00144-3
Devanne H, Lavoie BA, Capaday C: Input–output properties and gain changes in the human corticospinal pathway. Exp Brain Res 1997, 114: 329-338. 10.1007/PL00005641
Olejnik S, Algina J: Generalized eta and omega squared statistics: measures of effect size for some common research designs. Psychol Methods 2003, 8: 434-447.
Gandevia SC, Petersen N, Butler JE, Taylor JL: Impaired response of human motoneurones to corticospinal stimulation after voluntary exercise. J Physiol 1999, 521: 749-759. 10.1111/j.1469-7793.1999.00749.x
Khedr EM, Galal O, Said A, Abd-elsameea M, Rothwell JC: Lack of post-exercise depression of corticospinal excitability in patients with Parkinson’s disease. Eur J Neurol 2007, 14: 793-796. 10.1111/j.1468-1331.2007.01858.x
Milanovic S, Filipovic SR, Blesic S, Ilic TV, Dhanasekaran S, Ljubisavljevic M: Paired-associative stimulation can modulate muscle fatigue induced motor cortex excitability changes. Behav Brain Res 2011, 223: 30-35. 10.1016/j.bbr.2011.04.013
Klass M, Levenez M, Enoka RM, Duchateau J: Spinal mechanisms contribute to differences in the time to failure of submaximal fatiguing contractions performed with different loads. J Neurophysiol 2008, 99: 1096-1104. 10.1152/jn.01252.2007
Sogaard K, Gandevia SC, Todd G, Petersen NT, Taylor JL: The effect of sustained low-intensity contractions on supraspinal fatigue in human elbow flexor muscles. J Physiol 2006, 573: 511-523. 10.1113/jphysiol.2005.103598
Rossini PM, Barker AT, Berardelli A, Caramia MD, Caruso G, Cracco RQ, Dimitrijevic MR, Hallett M, Katayama Y, Lucking CH, de Noordhout ALM, Marsden CD, Murray NMF, Rothwell JC, Swash M, Tomberg C: Non-invasive electrical and magnetic stimulation of the brain, spinal cord and roots: basic principles and procedures for routine clinical application. Report of an IFCN committee. Electroencephalogr Clin Neurophysiol 1994, 91: 79-92. 10.1016/0013-4694(94)90029-9
Kiers L, Cros D, Chiappa KH, Fang J: Variability of motor potentials evoked by transcranial magnetic stimulation. Electroencephalogr Clin Neurophysiol 1993, 89: 415-423. 10.1016/0168-5597(93)90115-6
Darling WG, Wolf SL, Butler AJ: Variability of motor potentials evoked by transcranial magnetic stimulation depends on muscle activation. Exp Brain Res 2006, 174: 376-385. 10.1007/s00221-006-0468-9
Barker AT: The history and basic principles of magnetic nerve stimulation. Electroencephalogr Clin Neurophysiol 1999,51(Suppl):3-21.
Nudo RJ, Milliken GW, Jenkins WM, Merzenich MM: Use-dependent alterations of movement representations in primary motor cortex of adult squirrel monkeys. J Neurosci 1996, 16: 785-807.
Awiszus F, Wahl B, Meinecke I: Influence of stimulus cross talk on results of the twitch-interpolation technique at the biceps brachii muscle. Muscle Nerve 1997, 20: 1187-1190. 10.1002/(SICI)1097-4598(199709)20:9<1187::AID-MUS17>3.0.CO;2-Q
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: 1036-1044. 10.1152/japplphysiol.00103.2006
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: 1199-1209. 10.1152/japplphysiol.01246.2007
Cirillo J, Lavender AP, Ridding MC, Semmler JG: Motor cortex plasticity induced by paired associative stimulation is enhanced in physically active individuals. J Physiol 2009, 587: 5831-5842. 10.1113/jphysiol.2009.181834
Kidgell DJ, Pearce AJ: Corticospinal properties following short-term strength training of an intrinsic hand muscle. Hum Mov Sci 2010, 29: 631-641. 10.1016/j.humov.2010.01.004
Sale MV, Ridding MC, Nordstrom MA: Cortisol inhibits neuroplasticity induction in human motor cortex. J Neurosci 2008, 28: 8285-8293. 10.1523/JNEUROSCI.1963-08.2008
Todd G, Flavel SC, Ridding MC: Low-intensity repetitive transcranial magnetic stimulation decreases motor cortical excitability in humans. s 2006, 101: 500-505.
We sincerely acknowledge the assistance of Dr Léonard Féasson and Dr Pascal Edouard for conducting the medical inclusions and Régis Bonnefoy for his technical expertise. John Temesi was supported by a doctoral research grant from the Rhône-Alpes Region. Mathieu Gruet was supported by a grant from the Fonds de Dotation Recherche en Santé Respiratoire, appel d’offres 2011 for research on the brain during exercise.
The authors declare that they have no competing interests.
JT, MG, TR, SV and GYM all contributed to the project conception, data collection and analysis and manuscript preparation. All authors read and approved the final manuscript.