Skip to main content

The reliability and validity of rapid transcranial magnetic stimulation mapping for muscles under active contraction

Abstract

Rapid mapping is a transcranial magnetic stimulation (TMS) mapping method which can significantly reduce data collection time compared to traditional approaches. However, its validity and reliability has only been established for upper-limb muscles during resting-state activity. Here, we determined the validity and reliability of rapid mapping for non-upper limb muscles that require active contraction during TMS: the masseter and quadriceps muscles. Eleven healthy participants attended two sessions, spaced two hours apart, each involving rapid and ‘traditional’ mapping of the masseter muscle and three quadriceps muscles (rectus femoris, vastus medialis, vastus lateralis). Map parameters included map volume, map area and centre of gravity (CoG) in the medial-lateral and anterior-posterior directions. Low to moderate measurement errors (%SEMeas = 10–32) were observed across muscles. Relative reliability varied from good-to-excellent (ICC = 0.63–0.99) for map volume, poor-to-excellent (ICC = 0.11–0.86) for map area, and fair-to-excellent for CoG (ICC = 0.25–0.8) across muscles. There was Bayesian evidence of equivalence (BF’s > 3) in most map outcomes between rapid and traditional maps across all muscles, supporting the validity of the rapid mapping method. Overall, rapid TMS mapping produced similar estimates of map parameters to the traditional method, however the reliability results were mixed. As mapping of non-upper limb muscles is relatively challenging, rapid mapping is a promising substitute for traditional mapping, however further work is required to refine this method.

Highlights

We assessed the validity and reliability of rapid mapping for muscles that require active contraction during TMS: masseter and quadriceps.

Rapid TMS mapping showed equivalence in map parameters with the “gold standard”, however relative reliability was mixed.

Rapid TMS mapping is a promising substitute for traditional TMS mapping, but reliability of some parameters remains questionable.

Peer Review reports

Introduction

Transcranial magnetic stimulation (TMS) is a valuable tool for studying the structure and function of the primary motor cortex (M1) [1]. Using TMS, a “map” of the corticomotor representation of a target muscle can be obtained. These maps provide information about corticomotor excitability as well as the location of the corticomotor representation [2]. These measures are important for understanding neuroplasticity in healthy and clinical populations, including how corticomotor representations of muscles are affected by neurological disorders [3], motor task training [4; 5; 6] and pain [7; 8; 9; 10].

Traditionally, mapping is performed by delivering TMS pulses at multiple predefined sites, organised in  1 cm spaced rows and columns, with multiple stimuli delivered at each 1cm2 area [11; 12; 2]. However, the traditional mapping method can be lengthy because the experimenter must be spatially precise where TMS is delivered and cover the entirety of the surface representation. The recent development of neuronavigation software has allowed mapping to be as quick as two minutes using the “rapid mapping method” [13; 14]. In this method, TMS is delivered pseudorandomly within a specified area (e.g., a 6 × 6 cm2 area) rather than at predetermined grid sites. The motor evoked potential (MEP) amplitudes recorded at each stimulation location are then used to construct an estimated brain map via interpolation. This method has been shown to produce corticomotor maps that (a) have acceptable measurement error (absolute reliability), (b) have acceptable consistency between sessions (relative reliability), and (c) produce map parameter estimates similar to the traditional method (validity) in muscles of the upper limb and hand [15, 16, 13; 17; 18; 14].

Thus far, the rapid mapping technique has only been validated during resting muscle activity. However, certain muscles require active contraction to obtain reliable MEPs due to the depth of cortical representations, with two examples being the masseter and quadriceps muscles [11; 19]. Assessing the reliability and validity of rapid mapping for such muscles is vital for many areas of research. For example, mapping of the masseter and quadriceps has been useful for understanding the neurophysiology of chronic musculoskeletal pain conditions such as temporomandibular disorder [20] and patellofemoral pain [19], wherein alterations in motor representations of masseter and quadriceps muscles have been demonstrated respectively. Unfortunately, the time-consuming nature of traditional TMS mapping and the need to maintain muscle contraction can lead to muscle fatigue, potentially reducing data quality. Rapid TMS mapping can significantly reduce the acquisition time and overcome the shortcomings of the traditional method. However, whether corticomotor maps for the masseter and quadriceps obtained using the rapid method exhibit acceptable reliability and validity remains unknown.

Thus, this study aimed to determine the absolute reliability, relative reliability, and validity of corticomotor maps for the masseter and quadriceps muscles using rapid TMS mapping in healthy adults.

Methods

Participants

Eleven healthy participants (five females and six males, mean age:29 ± 8 years) were recruited using advertisements emailed throughout Neuroscience Research, Australia. Upon arrival, the participants completed a TMS safety screen [21]. Participants were excluded if they experienced any pain, were pregnant, had any major medical, neurological, or psychiatric conditions, reported a previous history of lower limb injuries or surgeries, or were taking psychoactive medication at the time of testing.

Experimental protocol

In line with a previous study, participants attended two sessions spaced two hours apart [13]. As we were interested in the reliability of the rapid mapping method itself, having a long intersession interval might lead to systematic changes in corticomotor excitability, making it difficult to disentangle whether differences between sessions are due to poor reliability or simply systematic changes in corticomotor excitability across time. For each session, maps were collected for the masseter and quadriceps muscles in a randomised order, using both traditional and rapid mapping methods (Fig. 1). One researcher (NC) conducted the masseter maps for all participants, while another (WJC) conducted the quadriceps maps for all participants. The researchers in this study are experienced in TMS mapping and have been involved in previous TMS masseter mapping and lower limb muscle mapping studies [20; 22].

Fig. 1
figure 1

Experimental design

Electromyography

Participants were seated in a comfortable chair and viewed a monitor that displayed electromyographic (EMG) feedback. Bipolar surface electrodes were used to record motor evoked potentials (MEPs) from the right masseter muscle using a belly tendon montage with active (muscle belly) and reference (tendon) electrodes placed along the mandibular angle, and the ground electrode was placed on the right acromion process [20]. For corticomotor representations of the quadriceps muscles, bipolar surface electrodes were placed over the belly of the right rectus femoris (RF), vastus lateralis (VL), and vastus medialis oblique (VMO), with the ground electrode placed at the right tibial shaft. The electrode for measuring RF muscle activity was placed 15 cm proximal to the superior border of the patella on the reference line (a line drawn between the centre of the patella and the anterior superior ischial spine) [23]. The electrode for measuring VL activity was placed 15 cm proximal to the superior border of the patella in a line 20° lateral to the reference line [24]. The electrode for measuring VMO muscle activity was placed 45° medial to the reference line, 4 cm proximal to the superior patellar border, and 4 cm lateral to the reference line [24]. EMG signals were amplified (x2000) using a CED 1902, filtered (20–1000 Hz), and digitally sampled at 2000 Hz using a Power 1401 Data Acquisition System, and Spike2 software (CED Limited, Cambridge, UK). Participants were seated with the hips and knees at 90 degrees of flexion during the entire mapping procedure. Prior to the TMS setup, participants performed three 3s maximal muscle contractions for the masseter and RF muscles to calculate the average root mean square (RMS) EMG during maximum voluntary contraction (MVC). The contraction target for the subsequent hotspot, threshold and mapping procedures were 20% MVC for the masseter [20] and 10% MVC for the RF muscle [25]. To maintain 10% MVC of the RF muscle, participants were instructed to perform isometric knee extension against a theraband. The participants received live visual feedback from the EMG signals (which was rectified and smoothed for display purposes) to maintain their contraction within a range ± 2.5% of the target contraction level.

TMS mapping procedures

While participants maintained the target contraction, the scalp sites evoking the largest MEP (the “hotspot”) for the masseter and RF muscles at a given TMS intensity were first determined. The active motor threshold (AMT) was determined using the TMS motor threshold assessment tool [26]. The AMT was defined as the minimum intensity required to evoke discernible MEPs. An MEP was identified if the EMG waveform was visibly larger in amplitude relative to the background EMG 5-15ms after the TMS pulse for the masseter [27; 28], and 26–46 ms for the RF muscle [29]. To minimise setup time, the hotspot locations and AMT values for the masseter and RF muscles in Session 1 were recorded and used in Session 2. Variations in these values within a two-hour intersession interval were expected to be negligible [5, 30]. The RF, VL, and VM were mapped concurrently in both sessions, as has been done previously [31].

Single monophasic stimuli were delivered using Magstim Rapid2 (Magstim Ltd., UK) and a 70 mm figure-of-eight coil at 120% AMT. An angle of 90° between the anterior-posterior line and medial-lateral line was used to induce a current in the lateral-to-medial direction [32; 33]. The Neural Navigator (Neurosoft, Russia) was used to track the positions of the TMS coil and the participant head in a 3-dimensional (3D) space. A model of the participant’s head was created by registering the position of anatomical landmarks (nasion, nose tip, and left and right ear tragus) into the software using a standard head-shape template.

TMS mapping was performed using a fixed 6 × 7 cm (7 rows and 8 columns) grid centred around the anticipated hotspot for the masseter (2 cm anterior and 6.5 cm lateral to the Cz) [34] and for the quadriceps muscles (0.5 cm anterior and 2 cm lateral to the Cz) [35]. For the traditional mapping method, three stimuli were delivered (interstimulus interval:4 s [13]) at each grid site (a total of 168 stimuli). For the rapid mapping method, 126 stimuli (interstimulus interval:2 s) were delivered using the pseudorandom walk method over the same grid as the traditional mapping. We based the number of stimuli from previous studies that used a stimulus density of 3 pulses per cm2 [13; 14]. This would also match the stimulus density of the traditional mapping method allowing for a fairer comparison. TMS stimuli were delivered in trains of three stimuli per grid site for the traditional mapping method and seven stimuli per train for the rapid mapping method. This was to provide breaks for participants in between contractions, thus minimizing muscle fatigue. As such, the interval between trains varied as the next train started when participants were ready for the next contraction. The neuronavigational display was monitored to ensure adequate coverage of the grid, and adjacent positions were not consecutively stimulated.

Data processing

The MEP amplitude at each grid site was calculated as the average RMS EMG amplitude of MEP traces, subtracting the background RMS EMG 55 to 5ms prior to the TMS pulse [36; 37], using a custom MATLAB script (MATLAB 7, MathWorks Inc., USA). A fixed MEP onset/offset window was used (26-46ms post-TMS for the quadriceps muscles, 7-15ms post-TMS for the masseter muscle) based on MEP latencies from previous studies using TMS mapping in the quadriceps and masseter muscles [20, 29] and confirmed using onsets/offsets that were manually determined for each participant. This fixed MEP window method has been shown to yield highly consistent TMS map parameters when compared to manually determined MEP onsets/offsets for each trial [20].

Maps were processed using a custom MATLAB Script. For the traditional mapping method, MEP amplitudes (mean of each location) were superimposed over the respective grid sites to construct a topographical map. For the rapid mapping method, the stimulation coordinates were first transformed in 3D space so that a 2-dimension (2D) plane fitted through the stimulation coordinates was parallel to the transverse plane. The 2D plane was then divided into 2744 (49 × 56) partitions, and triangular linear interpolation was used to approximate the MEP amplitude based on the nearest acquired MEP data [14]. Because there were 56 stimulation locations for the traditional mapping method, the interpolated plane was then divided into 56 equally sized sub-grids representing the same locations that were stimulated for the traditional mapping method. The mean MEP amplitude of each sub-grid was calculated, and these values were used to construct a topographical map. Figure 2 shows a schematic representation of the two mapping methods.

Fig. 2
figure 2

Schematic of map processing for the rapid and traditional mapping methods. For both methods, TMS was delivered over a 6 × 7 cm2 fixed area over the scalp. Traditional mapping consisted of three stimuli at 56 stimulation sites. The mean MEP amplitude at each site was calculated, and these values were used to create a map. The rapid mapping delivered 126 stimuli (3 stimuli/cm2). The 3D stimulation coordinates were transformed such that the 2D plane fitted through the coordinates was parallel to the transverse plane. This plane was then divided into 49 × 56 partitions, and MEP amplitude at each partition was interpolated using the nearest acquired MEP data. The mean MEP amplitude of each 7 × 8 sub-grid of partitions (representing the same locations as the traditional method) was then determined, and these values were used to create a map. Note the maps in this figure were generated from simulated data

The following map parameters were obtained for each map: (1) map volume, defined as the sum of all MEP amplitudes on the grid (interpolated grid in the case of rapid mapping) that were greater than 10% of the maximum MEP amplitude; (2) map area, defined as the number of sites (measured in cm2) on the grid which exhibited an MEP amplitude that was greater than 10% of the maximum MEP amplitude, and (3) centre of gravity (CoG), defined as the amplitude-weighted centre of the map in either the anterior-posterior direction or the medial-lateral direction, using the following formula: (CoG=\(\:\sum\:{V}_{i}x{\:X}_{i\:}/\sum\:{V}_{i}\), \(\:\sum\:{V}_{i}x{\:Y}_{i\:}/\sum\:{V}_{i}\); Vi=mean MEP amplitude at each site with the coordinates Xi, Yi) [38, 12].

Statistical analysis

Bland-Altman plots were used to assess homoscedasticity and the agreement between sessions and to visually inspect any systematic between-session differences and outliers for each map parameter of each target muscle [39; 40; 41; 42; 43]. When heteroscedasticity was identified, log transformation was performed and further analyses were conducted using log-transformed data.

Absolute reliability

Absolute Reliability was assessed using R, version 4.0.3 (R Development Core Team, Vienna, Austria) [44]. Absolute reliability, reflecting the measurement error and within-individual variability, was determined as the standard error of the measurement (SEMeas; a smaller SEMeas value indicates better measurement accuracy and lower within-individual variability) [39; 45; 42; 43]. The following formula was used: SEMeas=√mean squared error, where the mean squared error (or ‘residual error’) was obtained from one-way repeated measures analysis of variance applied to test and retest measurements [46]. SEMeas%, was also calculated (using the formula: SEMeas%= SEMeas/pooled mean*100%, where pooled mean was obtained from both testing sessions), to indicate the relative size of the measurement error [47; 42]. Additionally, the smallest detectable change (SDCindiv) at the individual level was derived from the SEMeas using the following formula: SDCindiv= SEMeas*1.96*√2 (for an individual), and minimal detectable change (MDC) for the group using the following formula: SDCgroup= SDCindiv/√n; n = sample size [48; 49; 42; 50]. Currently, there is no consensus on acceptable SDCs for TMS map parameters.

Relative reliability

The relative reliability was assessed using the MATLAB (Version 2021b) f_ICC function. To assess the relative reliability, the intraclass correlation coefficient (ICC) was evaluated using a single-rating, absolute-agreement, two-way mixed-effects model [ICC (2,1)] [39; 51; 43]. An ICC ≤ 0.2 indicates poor, 0.21 to 0.4 fair, 0.41 to 0.6 moderate, 0.61 to 0.8 good and > 0.8 excellent relative reliability [46].

Validity of rapid mapping method

To assess the validity of the rapid mapping method, we tested its equivalence against the traditional method [13] using equivalence Bayesian paired samples t-tests (paired-sample, equivalence interval: -0.25–0.25, non-overlapping hypotheses, Cauchy scale = 0.707) on JASP (Version 0.12.2.0, JASP Team, 2020) to compare the means of map volume, map area and CoG of the two sessions between the rapid and traditional mapping methods. A Bayes factor (BF) > 3 would provide at least moderate evidence that the rapid and traditional map parameters were inside the equivalence region, whereas a BF between 1 and 3 provided anecdotal/weak evidence that rapid and traditional map parameters were inside the equivalence region [52].

Results

Map characteristics

Owing to time constraints, one male participant did not complete the mapping protocols for the quadriceps muscles, and one female participant did not complete the mapping protocols for the masseter muscle. This led to ten participants (six males, four females) in the masseter muscle analysis, and ten participants (five females, five males) in the quadriceps muscle analysis. Table 1 shows the data for TMS procedural outcomes, including MVC, hotspot, AMT, and mapping times of both the rapid and traditional mapping methods. Rapid mapping was completed on average in  6–7 min, which was less than half the time required for traditional mapping ( 17 min). Table 2 lists the data for the map parameters of the rapid and traditional mapping methods for each muscle. Figure 3 shows examples of maps for each muscle from a representative participant obtained using the rapid and traditional mapping method.

Table 1 Measures relevant to rapid and traditional TMS mapping procedures for the masseter and rectus femoris muscles (mean, standard deviation and ranges)
Table 2 TMS map parameters (means and standard deviation) for the masseter and quadriceps muscles acquired by the rapid and traditional TMS mapping procedure
Fig. 3
figure 3

Map examples for each muscle obtained using the rapid and traditional mapping method. The axes represent the anterior-posterior and medial-lateral distance (in cm) relative to the vertex/Cz. Red and blue areas represent areas of higher and lower excitability respectively (measured in mV)

Reliability and SDCindv of rapid TMS mapping for the masseter and quadriceps muscles

The results for the measurement error and the minimal detectable change at the individual level (SDCindv) are presented in Table 3. Measurement errors for rapid mapping were low for the map area (%SEMeas=10%) but moderate for the map volume and the CoG location coordinates (%SEMeas=10–32%) for the masseter and quadriceps muscles, with the exception of the map volume for the RF muscle (%SEMeas=10%). The results for relative reliability are presented in Table 4. For map volume, the relative reliability of the rapid mapping procedure ranged from good to excellent for the quadriceps and excellent for the masseter muscles. For map area, relative reliability of rapid mapping was poor for the masseter, and ranged from moderate-excellent for the quadriceps muscles. For CoG, the relative reliability for rapid mapping ranged from fair to moderate for the quadriceps muscles and good to excellent for the masseter.

Table 3 Absolute reliability results for the rapid mapping procedure
Table 4 Intraclass correlation coefficient (ICC) data (and confidence intervals) for TMS map parameters of the masseter and quadriceps muscles acquired by the rapid mapping procedure

Validity of rapid TMS mapping for the masseter and quadriceps muscles

Table 5 presents the results of the Bayes factors for the equivalence tests between the rapid and traditional map parameters. There was moderate evidence of equivalence over non-equivalence between rapid and traditional mapping methods across all muscles (all BF’s > 3) with the exception of map area for the masseter, CoG Y for the masseter muscle and CoG X for the vastus medialis.

Table 5 Bayes factors (BF) from equivalence tests between the map parameters of the rapid and traditional mapping methods. The bold font represents substantial evidence for equivalence over non-equivalence (BF > 3)

Discussion

This study assessed the absolute reliability, relative reliability, and validity of the rapid TMS mapping method for the masseter and quadriceps muscles in healthy adults. The results showed low measurement error for map area, and mostly moderate measurement error for map volume and CoG. The relative reliability varied from good-to-excellent for map volume, poor-to-excellent for map area and fair-to-excellent for CoG coordinates. There was moderate Bayesian evidence of equivalence in the map parameters between the rapid and traditional maps in all muscles, supporting the validity of the rapid mapping method. Taken together, our findings provide reference values for minimal detectable changes in the masseter and quadriceps muscles in this population. While our results suggest the rapid mapping method produces similar estimates of map parameters to the traditional method, the reliability results of this method were mixed.

Absolute reliability of rapid TMS mapping

As there are no accepted norms for measurement error (SEMeas) or SDC, an arbitrary cut-off value of < 10% for SEMeas% has been proposed to reflect low measurement error [53; 42]. Our results showed that rapid mapping of the masseter and quadriceps muscles had moderate measurement errors for map volume and CoG location coordinates (%SEMeas=10–32%), except for map volume of the RF muscle (SEMeas%=10%). Consistent with previous reports on SDCs of TMS mapping measures [18; 54] and other TMS measures [42], TMS map measures for the masseter and quadriceps muscles had sizable SDCindv, precluding evaluative use for tracking changes within an individual. Nevertheless, our data will enable future studies with similar cohort characteristics to estimate how much an individual (SDCindv) or group (SDCindv divided by the square root of their sample size) would have to change in TMS map measures to be considered a real change exceeding the measurement errors [55]. As SDCs of TMS map measures for the masseter and quadriceps muscles become sufficiently low in modest sample sizes, they can be used for evaluative purposes to detect changes at the group level [42].

Research assessing the reliability of TMS mapping that requires muscle contractions during the procedure is scarce [56] and no study has examined the measurement error (absolute reliability). Higher TMS intensities are needed for mapping the facial and lower limb muscles at rest because of the depth of the M1 representation and the strength of corticospinal connections, which may not be tolerated by participants [56]. Thus, muscle contraction is necessary to allow mapping of these muscles with lower TMS intensities and better procedure tolerability [20; 29; 57]. The moderate measurement error of the TMS map for the masseter and quadriceps muscles may be attributed to between-session variation in the level of muscle contraction. Although 10% of MVC was set for the RF and 20% for the masseter muscles tailored to each participant and visual feedback of EMG was provided to maintain the target contraction during mapping, between-session variations in MVC within individuals may have influenced corticomotor representations of the mapped muscles.

Relative reliability of rapid TMS mapping

Our results suggest that rapid mapping yields good-to-excellent relative reliability for the estimation of map volume. This aligns with previous studies demonstrating the excellent relative reliability of rapid mapping of the upper limb muscles [15, 16, 13]. Relative reliability results for the other outcomes were mixed: the quadriceps muscles showed moderate to excellent relative reliability for the map area but fair relative reliability for CoG, while the masseter muscle showed good to excellent relative reliability for CoG but poor relative reliability for the map area. These findings are contrary to those from other studies demonstrating excellent relative reliability for the map area and CoG estimates for the rapid mapping of upper limb muscles [17; 14]. One explanation for this could be the low number of stimuli required for rapid mapping of the quadriceps and masseter muscles. In this study, we chose 126 stimuli to match the three pulses per cm2 in traditional mapping; however, achieving equally distributed pulses across the grid space using the pseudorandom walk method is not always possible. Non-upper limb muscles have corticomotor representations that differ from other muscles, for example, a different number of discrete peaks of corticomotor excitability [58; 19]. It is conceivable that a higher number of pulses would increase the researcher’s ability to more evenly distribute TMS pulses within the grid. This may produce a more reliable interpolation of the map and improve the identification of active sites, thus producing more accurate estimates of the map area and CoG location coordinates [56]. This may be further aided by the use of robotic neuronavigated TMS maps [59] to ensure that the grid space contains a sufficient number and distribution of TMS pulses. In any case, despite the unexpected results for the map area and the CoG coordinates, rapid TMS mapping appears to produce stable measures of map volume for both masseter and knee muscles.

Validity of rapid TMS mapping

Another aim of our study was to determine whether rapid mapping produces similar map estimates to the “gold-standard” traditional approach [12] where stimuli are delivered at 1 cm spaced points around the grid space. We used the same gold standard as Cavaleri et al. [13] where stimuli were delivered with a 4 s interstimulus interval. Our results showed that map volume, map area and most CoG assessments of the four muscles were equivalent between the rapid and traditional mapping methods. Although three comparisons yielded insufficient evidence, two of these (Map Area for Masseter, CoG X for VM) showed close to moderate evidence of equivalence (BF = 2.94–2.97). Overall, this suggests that while some map estimates obtained using the rapid mapping approach yield questionable relative reliability, they nonetheless produce similar map estimates to the “gold standard”. These results are mostly consistent with previous studies that showed equivalence between rapid and traditional maps [13; 14]. However, one previous study compared rapid and traditional mapping using null hypothesis testing [14], which does not appropriately assess equivalence, whereas the other used a frequentist approach for assessing equivalence [13], which requires corrections for multiple comparisons. A major strength of our study is the use of a Bayesian approach to assess equivalence. This approach allows researchers to quantify the strength of evidence for equivalence (rather than the accept/reject null hypothesis approach) [52] and does not require significance correction for multiple comparisons [60]. In summary, our findings support the validity of the rapid method and can be used instead of the traditional method.

Limitations and future research suggestions

This study has several limitations. Although our sample size of 10 is similar to a previous reliability study on the rapid mapping method [14], this could have been optimized using a sample size calculation. In situations where sample size is not calculated, it has been recommended to provide confidence intervals or posterior Bayesian distributions [61] – the former is shown in Table 4 and latter has been added to the supplementary file. Another limitation is the lack of trial exclusion based on coil orientation/angle, as these data were not available from our neuronavigation system. Trials with inaccurate coil positioning may have influenced the reliability of the rapid mapping data. However, we note that exclusion due to coil angle has previously been found to occur in a small percentage of trials (3.3%) [14]. Moreover, the easy-to-maintain coil orientation of 90° to the anterior-posterior line along with the large amount of TMS mapping experience of the experimenters in our study is likely to contribute to even fewer errors. Nonetheless, we advise future studies to assess rapid mapping of masseter and quadriceps muscles accounting for incorrect orientation and angles. Another consideration is coil orientation. Whereas the present study induced a lateral-medial current for all muscles, some studies oriented the coil that generated anterior-posterior direction currents for mapping the quadriceps considering that this coil orientation could stimulate the corticospinal pathway to the lower limb muscles more efficiently [29]. However, it should be noted that this rationale is based on studies in the lower leg and foot muscles [62] and other studies demonstrated that the coil should be oriented to generate medial-lateral currents for optimal stimulation for the lower leg muscle [63; 64]. As altering coil orientation can influence the motor threshold and stimulate different populations of neurons [63; 64], further investigation is needed to examine the optimal coil orientation for mapping M1 corticomotor representations of the quadriceps muscles. Lastly, the study did not assess reliability and validity using different ISIs, as has been done previously [13]. Future studies are encouraged to re-assess the reliability and validity of rapid mapping of the knee/masseter muscles with different ISIs, as the optimal ISI may differ for different muscles.

Conclusions

The process of acquiring TMS-evoked responses from the masseter and quadriceps muscles can be relatively challenging, as MEPs are typically collected while these muscles are under active contraction, which can lead to fatigue, especially during lengthy mapping sessions. Establishing methodologies that reduce map acquisition time for these muscles may improve the quality and participant tolerability of future experiments using masseter or quadriceps mapping. Here, we showed that maps produced using the rapid mapping method had modest absolute reliability, equivalence with the traditional mapping method for most outcomes/muscles, good to excellent relative reliability for map volume, but mixed relative reliability for other measures. Thus while rapid TMS mapping could be a promising substitute for traditional mapping of the non-upper limb muscles, further work is required to refine this methodology.

Data availability

The data that support the findings of this study are available from the corresponding author, NC, upon reasonable request.

References

  1. Schieber MH. Constraints on somatotopic organization in the primary motor cortex. J Neurophysiol. 2001;86(5):2125–43.

    Article CAS PubMed Google Scholar

  2. Wilson S, Thickbroom G, Mastaglia F. Transcranial magnetic stimulation mapping of the motor cortex in normal subjects: the representation of two intrinsic hand muscles. J Neurol Sci. 1993;118(2):134–44.

    Article CAS PubMed Google Scholar

  3. Lüdemann-Podubecká J, Nowak DA. Mapping cortical hand motor representation using TMS: a method to assess brain plasticity and a surrogate marker for recovery of function after stroke? Neurosci Biobehavioral Reviews. 2016;69:239–51.

    Article Google Scholar

  4. Kleim JA, Kleim ED, Cramer SC. Systematic assessment of training-induced changes in corticospinal output to hand using frameless stereotaxic transcranial magnetic stimulation. Nat Protoc. 2007;2(7):1675–84.

    Article CAS PubMed Google Scholar

  5. Lu S, Baad-Hansen L, Zhang Z, Svensson P. One hour jaw muscle training does not evoke plasticity in the corticomotor control of the masseter muscle. Arch Oral Biol. 2013;58(10):1483–90.

    Article PubMed Google Scholar

  6. Ngomo S, Leonard G, Moffet H, Mercier C. Comparison of transcranial magnetic stimulation measures obtained at rest and under active conditions and their reliability. J Neurosci Methods. 2012;205(1):65–71.

    Article PubMed Google Scholar

  7. Burns E, Chipchase LS, Schabrun SM. Altered function of intracortical networks in chronic lateral epicondylalgia. Eur J Pain. 2016;20(7):1166–75.

    Article CAS PubMed Google Scholar

  8. Chang WJ, Buscemi V, Liston MB, McAuley JH, Hodges PW, Schabrun SM. Sensorimotor cortical activity in Acute Low Back Pain: a cross-sectional study. J Pain. 2019;20(7):819–29.

    Article PubMed Google Scholar

  9. Chowdhury NS, Chang W-J, Millard SK, Skippen P, Bilska K, Seminowicz DA, Schabrun SM. The Effect of Acute and Sustained Pain on Corticomotor excitability: a systematic review and Meta-analysis of group-and individual-level data. J Pain. 2022;23(10):1680–96.

    Article PubMed Google Scholar

  10. Tsao H, Danneels LA, Hodges PW. ISSLS Prize winner: smudging the motor brain in young adults with recurrent low back pain. Spine (Phila Pa 1976). 2011;36(21):1721–7.

    Article PubMed Google Scholar

  11. McMillan AS, Watson C, Walshaw D. Transcranial magnetic-stimulation mapping of the cortical topography of the human masseter muscle. Arch Oral Biol. 1998;43(12):925–31.

    Article CAS PubMed Google Scholar

  12. Wassermann EM, McShane LM, Hallett M, Cohen LG. Noninvasive mapping of muscle representations in human motor cortex. Electroencephalogr Clin Neurophysiol. 1992;85(1):1–8.

    Article CAS PubMed Google Scholar

  13. Cavaleri R, Schabrun SM, Chipchase LS. The reliability and validity of rapid transcranial magnetic stimulation mapping. Brain Stimul. 2018;11(6):1291–5.

    Article PubMed Google Scholar

  14. Van De Ruit M, Perenboom MJ, Grey MJ. TMS brain mapping in less than two minutes. Brain Stimul. 2015;8(2):231–9.

    Article PubMed Google Scholar

  15. Cavaleri R, Chipchase LS, Massé-Alarie H, Schabrun SM, Shraim MA, Hodges PW. Corticomotor reorganization during short‐term visuomotor training in the lower back: a randomized controlled study. Brain Behav. 2020;10(8):e01702.

    Article PubMed PubMed Central Google Scholar

  16. Cavaleri R, Chipchase LS, Summers SJ, Schabrun SM. Repetitive transcranial magnetic stimulation of the primary motor cortex expedites recovery in the transition from acute to sustained experimental pain: a randomised, controlled study. Pain. 2019;160(11):2624–33.

    Article CAS PubMed Google Scholar

  17. Jin F, Bruijn SM, Daffertshofer A. Accounting for stimulations that do not elicit motor-evoked potentials when mapping cortical representations of multiple muscles. Front Hum Neurosci 2022;16.

  18. Jonker ZD, van der Vliet R, Hauwert CM, Gaiser C, Tulen JH, van der Geest JN, Donchin O, Ribbers GM, Frens MA, Selles RW. TMS motor mapping: comparing the absolute reliability of digital reconstruction methods to the golden standard. Brain Stimul. 2019;12(2):309–13.

    Article PubMed Google Scholar

  19. Te M, Baptista AF, Chipchase LS, Schabrun SM. Primary motor cortex organization is altered in persistent patellofemoral pain. Pain Med. 2017;18(11):2224–34.

    Article PubMed Google Scholar

  20. Chowdhury NS, Skippen P, Si E, Chiang AKI, Millard SK, Furman AJ, Chen S, Schabrun SM, Seminowicz DA. The reliability of two prospective cortical biomarkers for pain: EEG peak alpha frequency and TMS corticomotor excitability. J Neurosci Methods. 2023;385:109766.

    Article CAS PubMed Google Scholar

  21. Rossi S, Antal A, Bestmann S, Bikson M, Brewer C, Brockmöller J, Carpenter LL, Cincotta M, Chen R, Daskalakis JD, Di Lazzaro V, Fox MD, George MS, Gilbert D, Kimiskidis VK, Koch G, Ilmoniemi RJ, Lefaucheur JP, Leocani L, Lisanby SH, Miniussi C, Padberg F, Pascual-Leone A, Paulus W, Peterchev AV, Quartarone A, Rotenberg A, Rothwell J, Rossini PM, Santarnecchi E, Shafi MM, Siebner HR, Ugawa Y, Wassermann EM, Zangen A, Ziemann U, Hallett M. Safety and recommendations for TMS use in healthy subjects and patient populations, with updates on training, ethical and regulatory issues: Expert guidelines. Clin Neurophysiol. 2021;132(1):269–306.

    Article PubMed Google Scholar

  22. Moukhaiber N, Summers SJ, Opar D, Imam J, Thomson D, Chang W-J, Andary T, Cavaleri R. The effect of theta burst stimulation over the primary motor cortex on experimental hamstring pain: a randomized, controlled study. J Pain. 2023;24(4):593–604.

    Article PubMed Google Scholar

  23. Stensdotter AK, Hodges PW, Mellor R, Sundelin G, Hager-Ross C. Quadriceps activation in closed and in open kinetic chain exercise. Med Sci Sports Exerc. 2003;35(12):2043–7.

    Article PubMed Google Scholar

  24. Coqueiro KRR, Bevilaqua-Grossi D, Bérzin F, Soares AB, Candolo C, Monteiro-Pedro V. Analysis on the activation of the VMO and VLL muscles during semisquat exercises with and without hip adduction in individuals with patellofemoral pain syndrome. J Electromyogr Kinesiol. 2005;15(6):596–603.

    Article PubMed Google Scholar

  25. 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(1):1–13.

    Article Google Scholar

  26. Awiszus F, Borckardt J. TMS motor threshold assessment tool (MTAT 2.0). USA: Brain Stimulation Laboratory, Medical University of South Carolina; 2011.

    Google Scholar

  27. Macaluso GM, Pavesi G, Bonanini M, Mancia D, Gennari PU. Motor-evoked potentials in masseter muscle by electrical and magnetic stimulation in intact alert man. Arch Oral Biol. 1990;35(8):623–8.

    Article CAS PubMed Google Scholar

  28. Pavesi G, Macaluso GM, Tinchelli S, Medici D, Gemignani F, Mancia D. Magnetic motor evoked potentials (MEPs) in masseter muscles. Electromyogr Clin Neurophysiol. 1991;31(5):303–9.

    CAS PubMed Google Scholar

  29. Davies JL. Using transcranial magnetic stimulation to map the cortical representation of lower-limb muscles. Clin Neurophysiol Pract. 2020;5:87–99.

    Article PubMed PubMed Central Google Scholar

  30. Temesi J, Ly SN, Millet GY. Reliability of single-and paired-pulse transcranial magnetic stimulation for the assessment of knee extensor muscle function. J Neurol Sci. 2017;375:442–9.

    Article PubMed Google Scholar

  31. Cavaleri R. Exploring early corticomotor reorganisation. 2022.

  32. Guggisberg AG, Dubach P, Hess CW, Wüthrich C, Mathis J. Motor evoked potentials from masseter muscle induced by transcranial magnetic stimulation of the pyramidal tract: the importance of coil orientation. Clin Neurophysiol. 2001;112(12):2312–9.

    Article CAS PubMed Google Scholar

  33. Kesar TM, Stinear JW, Wolf SL. The use of transcranial magnetic stimulation to evaluate cortical excitability of lower limb musculature: challenges and opportunities. Restor Neurol Neurosci. 2018;36(3):333–48.

    PubMed PubMed Central Google Scholar

  34. Costa YM, Exposto FG, Kothari M, Castrillon EE, Conti PCR, Bonjardim LR, Svensson P. Masseter corticomotor excitability is decreased after intramuscular administration of nerve growth factor. Eur J Pain 2019;0(ja).

  35. Schwenkreis P, Pleger B, Cornelius B, Weyen U, Dertwinkel R, Zenz M, Malin J-P, Tegenthoff M. Reorganization in the ipsilateral motor cortex of patients with lower limb amputation. Neurosci Lett. 2003;349(3):187–90.

    Article CAS PubMed Google Scholar

  36. Chowdhury NS, Chiang AK, Millard SK, Skippen P, Chang W-J, Seminowicz DA, Schabrun SM. Combined transcranial magnetic stimulation and electroencephalography reveals alterations in cortical excitability during pain. Elife. 2023;12:RP88567.

    Article CAS PubMed PubMed Central Google Scholar

  37. Tsao H, Tucker KJ, Hodges PW. Changes in excitability of corticomotor inputs to the trunk muscles during experimentally-induced acute low back pain. Neuroscience. 2011;181:127–33.

    Article CAS PubMed Google Scholar

  38. Uy J, Ridding M, Miles T. Stability of maps of Human Motor Cortex Made with Transcranial magnetic stimulation. Brain Topogr. 2002;14(4):293–7.

    Article PubMed Google Scholar

  39. Atkinson G, Nevill AM. Statistical methods for assessing measurement error (reliability) in variables relevant to Sports Medicine. Sports Med. 1998;26(4):217–38.

    Article CAS PubMed Google Scholar

  40. Damron LA, Dearth DJ, Hoffman RL, Clark BC. Quantification of the corticospinal silent period evoked via transcranial magnetic stimulation. J Neurosci Methods. 2008;173(1):121–8.

    Article PubMed Google Scholar

  41. Martin Bland J, Altman D, Statistical methods for assessing, agreement between two methods of clinical measurement. Lancet. 1986;327(8476):307–10.

    Article Google Scholar

  42. Schambra H, Ogden RT, Martínez-Hernández I, Lin X, Chang YB, Rahman A, Edwards D, Krakauer J. The reliability of repeated TMS measures in older adults and in patients with subacute and chronic stroke. Front Cell Neurosci 2015;9(335).

  43. Weir JP. Quantifying test-retest reliability using the intraclass correlation coefficient and the SEM. J Strength Cond Res. 2005;19(1):231–40.

    PubMed Google Scholar

  44. Core Team R. R: a language and environment for statistical computing. Vienna: R Foundation for statistical computing; 2013.

    Google Scholar

  45. Beaulieu LD, Flamand VH, Masse-Alarie H, Schneider C. Reliability and minimal detectable change of transcranial magnetic stimulation outcomes in healthy adults: a systematic review. Brain Stimul. 2017;10(2):196–213.

    Article PubMed Google Scholar

  46. Matamala JM, Howells J, Dharmadasa T, Trinh T, Ma Y, Lera L, Vucic S, Burke D, Kiernan MC. Inter-session reliability of short-interval intracortical inhibition measured by threshold tracking TMS. Neurosci Lett. 2018;674:18–23.

    Article CAS PubMed Google Scholar

  47. Lexell JE, Downham DY. How to assess the reliability of measurements in rehabilitation. Am J Phys Med Rehabil. 2005;84(9):719–23.

    Article PubMed Google Scholar

  48. Beckerman H, Roebroeck ME, Lankhorst GJ, Becher JG, Bezemer PD, Verbeek ALM. Smallest real difference, a link between reproducibility and responsiveness. Qual Life Res. 2001;10(7):571–8.

    Article CAS PubMed Google Scholar

  49. de Vet HC, Terwee CB, Knol DL, Bouter LM. When to use agreement versus reliability measures. J Clin Epidemiol. 2006;59(10):1033–9.

    Article PubMed Google Scholar

  50. Schuck P, Zwingmann C. The ‘smallest real difference’ as a measure of sensitivity to change: a critical analysis. Int J Rehabil Res. 2003;26(2):85–91.

    PubMed Google Scholar

  51. Trevethan R. Intraclass correlation coefficients: clearing the air, extending some cautions, and making some requests. Health Serv Outcomes Res Method. 2017;17(2):127–43.

    Article Google Scholar

  52. van Doorn J, van den Bergh D, Böhm U, Dablander F, Derks K, Draws T, Etz A, Evans NJ, Gronau QF, Haaf JM. The JASP guidelines for conducting and reporting a bayesian analysis. Psychon Bull Rev. 2021;28:813–26.

    Article PubMed Google Scholar

  53. Flansbjer UB, Holmback AM, Downham D, Patten C, Lexell J. Reliability of gait performance tests in men and women with hemiparesis after stroke. J Rehabil Med. 2005;37(2):75–82.

    Article PubMed Google Scholar

  54. Nazarova M, Novikov P, Ivanina E, Kozlova K, Dobrynina L, Nikulin VV. Mapping of multiple muscles with transcranial magnetic stimulation: absolute and relative test-retest reliability. Hum Brain Mapp. 2021;42(8):2508–28.

    Article PubMed PubMed Central Google Scholar

  55. Terwee CB, Bot SD, Boer MR, Windt DA, Knol DL, Dekker J. Quality criteria were proposed for measurement properties of health status questionnaires. J Clin Epidemiol 2007;60.

  56. Sondergaard RE, Martino D, Kiss ZHT, Condliffe EG. TMS Motor Mapping Methodology and reliability: a structured review. Front NeuroSci 2021;15.

  57. Ward S, Bryant AL, Pietrosimone B, Bennell KL, Clark R, Pearce AJ. Cortical motor representation of the rectus femoris does not differ between the left and right hemisphere. J Electromyogr Kinesiol. 2016;28:46–52.

    Article PubMed Google Scholar

  58. Chang W-J, Bennell KL, Hodges PW, Hinman RS, Liston MB, Schabrun SM. Combined exercise and transcranial direct current stimulation intervention for knee osteoarthritis: protocol for a pilot randomised controlled trial. BMJ open. 2015;5(8):e008482.

    Article PubMed PubMed Central Google Scholar

  59. Kahl CK, Giuffre A, Wrightson JG, Kirton A, Condliffe EG, MacMaster FP, Zewdie E. Active versus resting neuro-navigated robotic transcranial magnetic stimulation motor mapping. Physiological Rep. 2022;10(12):e15346.

    Article Google Scholar

  60. Gelman A, Tuerlinckx F. Type S error rates for classical and bayesian single and multiple comparison procedures. Comput Stat. 2000;15(3):373–90.

    Article Google Scholar

  61. Dziak JJ, Dierker LC, Abar B. The interpretation of statistical power after the data have been gathered. Curr Psychol. 2020;39:870–7.

    Article PubMed Google Scholar

  62. 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(5):858–82.

    Article CAS PubMed PubMed Central Google Scholar

  63. Richter L, Neumann G, Oung S, Schweikard A, Trillenberg P. Optimal Coil Orientation for Transcranial magnetic stimulation. PLoS ONE. 2013;8(4):e60358.

    Article CAS PubMed PubMed Central Google Scholar

  64. Terao Y, Ugawa Y, Sakai K, Uesaka Y, Kohara N, Kanazawa I. Transcranial stimulation of the leg area of the motor cortex in humans. Acta Neurol Scand. 1994;89(5):378–83.

    Article CAS PubMed Google Scholar

Download references

Acknowledgements

No further acknowledgements.

Funding

This work was supported by 1R61NS113269-01 from The National Institutes of Health to SMS.

Author information

Authors and Affiliations

Authors

Contributions

N.S: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Project Administration, Writing – Original Draft, Writing – Reviewing and Editing. W.C.: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Project Administration, Writing – Original Draft, Writing – Reviewing and Editing. R.C: Conceptualization, Resources, Software and Writing – Reviewing and Editing.A.C: Conceptualization, Methodology and Writing – Reviewing and Editing. S.S: Conceptualization, Project Administration, Resources, Software, Supervision and Writing – Reviewing and Editing.

Corresponding authors

Correspondence to Nahian S. Chowdhury or Wei-Ju Chang.

Ethics declarations

Ethics approval and consent to participate

The study procedures adhered to the Declaration of Helsinki and were approved by the Human Research Ethics Committee of the University of New South Wales (HC200555). All the participants provided written informed consent.

Consent for publication

The authors give their consent for the publication of the material in this paper.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chowdhury, N.S., Chang, WJ., Cavaleri, R. et al. The reliability and validity of rapid transcranial magnetic stimulation mapping for muscles under active contraction. BMC Neurosci 25, 43 (2024). https://doi.org/10.1186/s12868-024-00885-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12868-024-00885-w

Keywords