Abstract
During unilateral contractions, muscular activation can be detected in both active and resting limbs. In healthy populations, the latter is referred to as physiological mirror activity (pMA). The study of pMA holds implications for clinical applications as well as the understanding of bilateral motor control. However, the underlying mechanisms of pMA remain to be fully resolved. A commonality of prevailing explanatory approaches is the concept of shared neural input. With this study, we, therefore, aimed to investigate neural input in the form of multiple analyses of surface electromyography (sEMG) recordings in the frequency domain. For this purpose, 14 healthy, right-handed males aged 18â35Â years were recruited. All participants performed a pinch-force task with the dominant hand in a blockwise manner. In total, 9 blocks of 5 contractions each were completed at 80% of maximum force output. Muscle activity was recorded via sEMG of the first dorsal interosseous muscle of the active and resting hand. We analyzed (1) spectral features as well as (2) intermuscular coherence (IMC). Our results demonstrate a blockwise increase in median frequency, mean frequency, and peak frequency in both hands. Frequency ratio analyses revealed a higher low-frequency component in the resting hand. Although we were able to demonstrate IMC on an individual level, results varied greatly and grand-averaged IMC failed to reach significance. Taken together, our findings imply an overlap of spectral properties between active and passive hands during repeated unilateral contractions. Combined with evidence from previous studies, this suggests a common neural origin between active and resting hands during unilateral contractions possibly resulting from a reduction in interhemispheric inhibition due to high force demands. Nevertheless, the exploratory nature of this study necessitates the classification of our results through follow-up studies.
Similar content being viewed by others
Introduction
Physiological mirror activity (pMA) describes unintentional muscle activity in resting limbs during unilateral contractions in neurologically healthy populations, that can be observed through surface electromyographic (sEMG) recordings. Inherent properties of pMA have been primarily studied by way of time-domain analyses of sEMG. To date, it is well known and often replicated, that pMA amplitude changes as a function of the amount of applied force during unilateral movements1,2,3,4,5,6,7,8. Factors influencing pMA amplitude apart from force requirements include increased movement frequency9 and rhythmicity10, increased cognitive load induced by simultaneously provided visual stimuli during task execution11, executive and attentional processes12,13, as well as central and peripheral fatigue resulting from repetitive exhaustive contractions14,15,16,17,18. Recently, it was further demonstrated that pMA is not time-locked to the onset of voluntarily contracting hand muscles but starts with varying and dynamically changing latencies, which are inversely correlated with respective amplitudes18.
Nevertheless, unresolved aspects of pMA remain, most prominently its neural origin15. Numerous explanatory approaches exist, for example, the theory of motor overflow19,20 or motor irradiation17, which assumes that pMA results from an ongoing modulation of interhemispheric communication during unilateral contractions with increasing force demands. These modulations lead to a shift from interhemispheric inhibition (IHI) to interhemispheric facilitation (IHF), which in turn induces bilateral activation of motor-relevant brain regions4,19. Alternative viewpoints suggest that ipsilateral projections of the corticospinal tract (CST) lead to bilateral activation in limbs both contralateral and ipsilateral to the cortical origin of the motor command15. A unifying feature of these theories is the question of the neural drive underlying active and resting limbs, i.e., the extent to which both effectors share neural input. In addition to time-domain analyses, frequency-domain analyses of sEMG allow insights into the neural drive to muscles by providing information on signal power distribution21, central and peripheral fatigue22, as well as fibre type characteristics23. Furthermore, quantifiable properties of sEMG frequency content, referred to as spectral features, have proven to be useful in the description and classification of physiological and pathophysiological muscle activities24. Among the commonly studied frequency features are median frequency (MDF), mean frequency (MNF), peak frequency (PKF), and frequency ratio (FR). MDF, MNF, and PKF are mainly used to assess fatigue-related effects22, whereas FR can be used to study the relative contribution of low and high frequency components to total signal power25. Each of these spectral features covers a particular aspect of the sEMG frequency content, thus aiding in the assessment of sEMG signals. Frequency analyses can be performed within a single signal, but can also be extended to multiple related signals. A well-established method in this context is intermuscular coherence (IMC). IMC enables the study of functional binding between muscles26. Such binding between muscles was first thought to reflect a common neural drive underlying functionally related muscle activities27. It was later demonstrated, that common oscillatory input to a set of motoneurons is present at the output of these motoneurons21 and more importantly emphasized, compared to other non-linear inputs, in the common neural drive to the muscle28,29. IMC can, therefore, be used to reveal the presence or absence of a common descending drive between functionally related muscles26.
In summary, a systematic analysis of (1) frequency content in active and resting hands during unilateral contractions, as well as (2) shared frequency content between active and resting hands has yet to be conducted. Such multimodal analyses in the frequency domain represent a promising extension to the study of pMA. Therefore, in this study, we aimed to investigate both the spectral features and IMC of active and resting limbs to contribute to the uncovering of their physiological origin. Here, we focussed on pMA of the upper extremities, particularly the first dorsal interosseous (FDI), as it is the best-studied muscle in relation to this phenomenon. Regarding spectral features, we hypothesized differences in terms of reduced MNF, MDF, and PKF as well as a higher proportion of low-frequency components in the resting compared to the active hand. This was mainly motivated by the fact that pMA has fundamentally lower amplitudes in the noncontracting muscle compared with the contracting muscle4,5. Furthermore, we hypothesized shared neural drive of both muscles and thus the presence of IMC based on contemporary theories on the origin of pMA. Perspectively, the results of this study apply to patients as well as healthy individuals, as the study of pMA has implications not only for clinical applications but also for understanding bimanual motor control in general.
Materials and methods
This study partly represents a reanalysis of previously published sEMG data18.
Ethical approval
The study was endorsed by the Ethics Committee at the Medical Faculty of Leipzig University (ref. no. 429/15-ff). All participants gave written informed permission to participate in the tests and were compensated for involvement. All methods were carried out in accordance with the Declaration of Helsinki.
Participants
Twenty-four healthy male participants were enrolled in this study. Since a detailed description of frequency-related parameters requires the presence of pMA, only participants with robust, visible pMA across all trials were analyzed in this study. Specifically, only participants with observable pMA, i.e. EMG activity exceeding mean baseline activityâ+â2 standard deviations for a time window of at least 10 ms, over each of the 45 contractions for at least 1 s of the 3-s contraction durations were selected to ensure actual analysis of pMA. Accordingly, fourteen healthy male participants were included in this study (age: (meanâ±âstandard deviation) 27.2â±â4.6 years). Only male participants were recruited to account for possible gender-related differences in brain structure and function30, which may have influenced the motor behavior under investigation. All participants were right-handed according to the Oldfield handedness inventory31 (laterality quotient: 86.4â±â18.9) and none of them had any history of playing musical instruments. None of the participants followed an organized training regime but performed general physical activities (3.9â±â2.1 h per week). In addition, participants were advised to prevent alcohol and caffeine consumption 24 h before the experimental sessions due to its well-known influence on force production and the functioning of the central nervous system32.
Behavioral experiment
Participants performed a unilateral force generation task. Therefore, the actively contracting (right) hand was used to operate a custom-built force sensor to perform an isometric pinch-task for which the thumb and index finger applied force to the force sensor while the resting (left) hand remained relaxed (see Fig. 1A). Force was applied on a single-point load cell (LAUMAS Elettronica Srl, Italy) with a net weight of 0.2 kg and platform dimensions of 25âÃâ30âÃâ35 mm. The maximum force capacity of the load cell is 50 kg with a combined error ofâ<â±â0.02%. Participants were instructed to focus solely on the active hand while there was no feedback from ongoing pMA to avoid deliberate inhibition of involuntary muscle activity33. The participants were not aware of the study interest at any time during the experiments.
Visual feedback during the pinch force task was provided on a PC using Presentation 16.5 (NeuroBehavioral Systems, Albany, USA). The screen showed a target field and a horizontal bar. The objective was to move the bar as rapidly and accurately as possible into the target field by exerting pressure on the force sensor. The target field, as well as the force required to reach it, were adjusted to individual maximum voluntary force values, representing 80% of maximum force.
A maximum force test was performed separately for both the right and left hands before the test. Participants exerted an individual maximum contraction 3 times (5 s duration for each repetition) with a resting period of 1 min between contractions to counteract fatigue and prevent injuries. The trial where participants exerted the highest force was taken as their individual maximum force. This trial was further used for the normalization of individual force levels during the behavioral task. The corresponding EMG amplitude to this trial was then taken as the individual maximum voluntary contraction (MVC) for further EMG MVC normalization. For this purpose, EMG activity over a window of 500 ms around the peak of the maximum contraction was averaged. The subsequent task used a block design consisting of a constant force level (80% of individual maximum force). Five isometric contractions were performed using the right hand within one block. One contraction lasted for 3 s, with 3 s of rest between contractions. As all participants performed 5 contractions per block, this resulted in a total duration of 30 s per block. Each block was followed by a resting period of 4:30 min to allow partial recovery. In total, each participant completed 9 blocks. Taken together, the time to finish the task was 45 min.
EMG recordings and analysis
sEMG was recorded on a wireless Desktop Direct Transmission System (NORAXON Inc., Scottsdale, USA). EMG signals were obtained from bilateral first dorsal interossei muscles (FDI) using bipolar surface electrodes (Ag/AgCl; diameter: 10 mm). Inter-electrode distances were consistent at 20 mm and electrodes were attached in parallel to muscle fiber orientation. This configuration allowed us to record sEMG activity over the voluntarily contracting FDI (FDIVol), as well as subliminal pMA over the homologous FDI (FDIpMA) of the relaxed hand. EMG data were recorded with a sampling frequency of 1500 Hz, band-pass filtered at 10â500 Hz, input impedanceâ>â100 MΩ, Common Mode Rejection Ratio (CMRR)â>â100 dB, and a gain of 500.
All EMG processing was performed using custom-written code implemented in MATLAB (v. R2021b, The MathWorks Inc., Natick, USA). Initially, all EMG data were decimated to 500Â Hz. For this purpose, we applied a type I Chebyshev filter (cutoff frequency of 200Â Hz) to the data before downsampling to avoid aliasing effects34, and subsequently, high-pass filtered EMG signals at 20Â Hz (4th order Butterworth filter) to remove low-frequency noise35. Additionally, both EMG signals (FDIVol and FDIpMA) were overlaid and time-locked to confirm the preservation of temporal relationships between both signals.
Time-domain analysis
Voluntary muscle onsets of FDIVol were manually defined for each burst through visual inspection of EMG traces performed by a single trained researcher. Participants were instructed to keep all muscles of the upper extremity as relaxed as possible in between every contraction. Onsets of FDIpMA were defined analogously to those of FDIVol. Mean EMG amplitudes of FDIVol, as well as FDIpMA, were determined by calculating the root mean square (RMS) values of each burst. To ensure comparable EMG amplitudes, all EMG amplitudes were normalized to individual MVC values (see âBehavioral experimentâ section), which were recorded separately for each hand at the beginning of each session. All bursts were subsequently pooled per block (BL 1âBL 9) for each participant and subjected to further statistical analyses.
Frequency-domain analysis: spectral features
Initially, power spectral densities (PSD) were estimated according to Welchâs method36. To characterize frequency-specific patterns of FDIVol and FDIpMA, four commonly reported frequency features were employed: mean frequency (MNF), median frequency (MDF), peak frequency (PKF), and frequency ratio (FR)25,37,38. Both MNF and MDF were estimated using the corresponding MATLAB functions, where MNF is the average frequency obtained by taking the sum of the product of the EMG power spectrum and the frequency and dividing both by the total sum of the power spectrum. MDF is defined as the frequency at which the EMG power spectrum is divided into two regions of equal amplitudes. PKF corresponds to the frequency of maximum energy content of the power spectrum and FR is calculated as a ratio between low- and high-frequency content of the EMG signals. Low- and high cutoff frequencies of FR can be defined either by experiment or MNF39. Here, we chose to set low and high frequency cutoffs after careful inspection of power spectral densities to account for the exploratory aspect of FDIpMA frequency content.
Frequency-domain analysis: intermuscular coherence
To analyze common input between active and resting hands, IMC was estimated between FDIVol and FDIpMA. Coherence is an extension of Pearsonâs correlation coefficient in the frequency domain and describes linear relationships between two input signals. Coherence is obtained through the normalization of the cross-spectra between two signals by their auto-spectra40. First, EMG signals were epoched per block and contractions yielding 45 trials (5 contractionsâÃâ9 blocks). All EMG bursts were then rectified by way of Hilbert transform. The amplitude of the Hilbert transform provides the envelope of the broadband EMG signal, gives similar results to full-wave rectification, and has been commonly employed to study intermuscular coherence41,42. In a final step, bursts were concatenated in a blockwise manner yielding 5 FDIVol and 5 corresponding FDIpMA bursts per block for each participant. This is done to reduce noise in the signals22. After further inspection of EMG signal quality, no bursts were excluded from further analysis, neither for FDIVol nor for FDIpMA. Intermuscular coherence was subsequently estimated between pairs of EMG data (FDIVol and FDIpMA) in this blockwise manner using Welchâs method with a Hanning window of 250 ms without overlap26. The significance of IMC results was determined by calculating confidence limits according to Rosenberg et al.43:
where α is the significance set to 5%, N is the number of epochs, and CL reflects the confidence limit above which observed coherence values are deemed significant. IMC estimates were then integrated over four relevant frequency bands: (1) α (8â12 Hz), (2) β (13â30 Hz), (3) γ (30â60 Hz), and (4) γⲠ(60â100 Hz) to yield IMC areas and used for further statistical analysis. These frequency bands were selected based on previous studies demonstrating that muscles receive oscillatory neuronal input spanning across these frequencies44,45. Analysis of the areas of coherence estimates is considered to be superior compared to analysis of peaks and frequencies of coherence estimates46,47,48.
Statistical analyses
All statistical analyses were conducted using the free and open-source statistical software program JASP (JASP Team (2021). JASP (Version 0.16.1.0) [Computer software]; Retrieved from https://jasp-stats.org). Since the majority of all variables were normally distributed, as assessed by ShapiroâWilk testing (αâ=â0.05), parametric analyses were implemented.
Amplitude
Blockwise amplitudes of FDIVol and FDIpMA were analyzed using a repeated-measures ANOVA with the within-subject factors HAND (FDIVol and FDIpMA), and BLOCK (BL 1âBL 9).
Frequency features
Concerning frequency features, blockwise MNF, MDF, PKF, and FR were analyzed using separate repeated-measures ANOVAs with the within-subject factors HAND (FDIVol and FDIpMA), and BLOCK (BL 1âBL 9).
Intermuscular coherence
Summed coherence values of each analyzed frequency band (α, β, γ, γâ²) were separately assessed using repeated-measures ANOVAs with the within-subject factor BLOCK (BL 1âBL 9).
For all analyses, the statistical threshold was set at pâ<â0.05. Effect sizes were expressed using partial eta squared (ηp2) for ANOVAs and Cohens d for Bonferroni corrected post-hoc comparisons. A GreenhouseâGeisser correction was implemented when the sphericity assumption of repeated-measures ANOVAs was violated.
Results
Amplitude
Repeated measures ANOVA indicated a significant effect for HAND (F(1, 13)â=â97.390, pâ<â0.001, ηp2â=â0.882) on sEMG amplitude (see Fig. 1). Post-hoc comparison showed that sEMG amplitudes of FDIVol were significantly higher compared to FDIpMA (mean difference (MD)â=â69.3%MVC, SEâ=â7.0%MVC, pâ<â0.001, dâ=â2.637). However, no significant effect was found for BLOCK (F(1.584, 20.598)â=â1.360, pâ=â0.274, ηp2â=â0.095) and no significant interaction effect HANDâÃâBLOCK (F(1.694, 22.016)â=â1.148, pâ=â0.327, ηp2â=â0.081) was observed.
Frequency features
Visual inspection of block-wise PSDs revealed a clear difference between FDIVol and FDIpMA in the low-frequency range up to 40 Hz (see Fig. 2). Therefore, the cutoff frequency between low and high-frequency components was set at 40 Hz, which was subsequently used to calculate FR. In addition, this cutoff is consistent with previous studies analyzing the relationship of low-frequency and high-frequency signal components in EMG37.
Regarding MNF, repeated-measures ANOVA revealed a significant effect for BLOCK (F(3.344, 43.467)â=â15.186, pâ<â0.001, ηp2â=â0.539; see Fig. 3A). Pairwise post-hoc comparisons indicated that from BL 4 onwards, MNF gradually increased compared to BL 1 (BL 1 vs. BL 4: MDâ=â4.73 Hz, SEâ=â1.31 Hz, pâ=â0.017, dâ=â0.964). No effect for HAND (F(1, 13)â=â3.971, pâ=â0.068, ηp2â=â0.234) or interaction effect HANDâÃâBLOCK (F(2.582, 33.564)â=â1.530, pâ=â0.228, ηp2â=â0.105) was found.
A similar relation was observed for MDF. Repeated-measures ANOVA showed a significant effect for BLOCK (F(3.581, 46.657)â=â13.301, pâ<â0.001, ηp2â=â0.506; see Fig. 3B). Pairwise post-hoc comparisons revealed that from BL 5 until BL 9, higher MDF compared to BL 1 were observed (BL 1 vs. BL 5: MDâ=â8.64 Hz, SEâ=â1.83 Hz, pâ<â0.001, dâ=â1.265). No effect for HAND (F(1, 13)â=â3.787, pâ=â0.074, ηp2â=â0.226) or interaction effect HANDâÃâBLOCK (F(2.763, 35.914)â=â1.361, pâ=â0.271, ηp2â=â0.095) was found.
When looking at PKF, repeated-measures ANOVA revealed a significant interaction effect HANDâÃâBLOCK (F(8, 104)â=â2.190, pâ=â0.034, ηp2â=â0.144) and a significant effect for BLOCK (F(8, 104)â=â7.043, pâ<â0.001, ηp2â=â0.351; see Fig. 3C). Simple main effect analysis showed that PKF of FDIVol was significantly higher than PKF of FDIpMA during BL 1 (89.8 Hz vs. 66.7 Hz, Fâ=â6.596, pâ=â0.023), BL 3 (112.4 Hz vs. 83.1 Hz, Fâ=â5.863, pâ=â0.031) and BL 4 (109.8 Hz vs. 82.2 Hz, Fâ=â7.026, pâ=â0.020). Pairwise post-hoc comparisons for BLOCK indicated that BL 6, BL 7, BL 8 and BL 9 had higher PKF compared to BL 1 (BL 1 vs. BL 6: MDâ=â37.02 Hz, SEâ=â6.228 Hz, pâ<â0.001, dâ=â1.589; BL 1 vs. BL 7: MDâ=â25.93 Hz, SEâ=â6.228 Hz, pâ=â0.002, dâ=â1.113; BL 1 vs. BL 8: MDâ=â25.77 Hz, SEâ=â6.228 Hz, pâ=â0.003, dâ=â1.106; BL 1 vs. BL 9: MDâ=â27.64 Hz, SEâ=â6.228 Hz, pâ<â0.001, dâ=â1.186). Again, no effect for HAND (F(1, 13)â=â2.230, pâ=â0.159, ηp2â=â0.156) was found.
In terms of FR, repeated-measures ANOVA revealed a significant effect for HAND (F(1, 13)â=â5.651, pâ=â0.033, ηp2â=â0.303) and BLOCK (F(2.383, 30.980)â=â5.469, pâ=â0.007, ηp2â=â0.296; see Fig. 3D). Post-hoc comparison of HAND indicated that FR of FDIVol was significantly lower compared to FDIpMA (MDâ=âââ0.053, SEâ=â0.022 Hz, pâ=â0.033, dâ=âââ0.653). Pairwise post-hoc comparisons of BLOCK revealed that from BL 7 onwards, FR gradually decreased compared to BL 1 (BL 1 vs. BL 7: MDâ=â0.039 Hz, SEâ=â0.011 Hz, pâ=â0.017, dâ=â0.967). However, no significant interaction effect HANDâÃâBLOCK (F(2.301, 29.914)â=â2.850, pâ=â0.067, ηp2â=â0.180) was found.
Intermuscular coherence
On an individual level, significant IMC was observable in all participants, although with a high inter-subject variability concerning the frequency band and block at which significant IMC was found (please see Fig. S1 for a detailed overview of individual IMC spectra). Critically, pooled IMC spectra between FDIVol and FDIpMA do not show significant grand-averaged coherence for any block in broadband sEMG (see Fig. 4).
Summed coherence values of the analyzed frequency bands did not show a blockwise modulation as indicated by repeated-measures ANOVAs with non-significant factor BLOCK (α-band: F(3.198, 41.576)â=â0.867, pâ=â0.472, ηp2â=â0.063; β-band: F(8, 104)â=â1.070, pâ=â0.390, ηp2â=â0.076; γ-band: F(8, 104)â=â0.727, pâ=â0.667, ηp2â=â0.053; γâ²-band: F(8, 104)â=â0.793, pâ=â0.610, ηp2â=â0.058).
Discussion
With this study, we aimed to investigate spectral features in active and resting hands during repeated unilateral contractions as well as IMC to potentially reveal common input between active and resting hands. Regarding frequency features of both FDIVol and FDIpMA, we found a significant increase in MNF, MDF, and PKF in both hands from blocks 4, 5, and 6, respectively, compared to block 1. We also observed a significant interaction between the factors HAND and BLOCK for PKF, where PKF was higher in FDIVol than FDIpMA during block 1 and block 4. No significant differences between both hands were found in either of these measures. For FR, results showed significantly higher values in FDIpMA compared to FDIVol as well as a gradual decrease in FR of both hands from block 7 to block 9. Intermuscular coherence was evident across participants on an individual level, however, grand-averaged IMC failed to show any significant pattern of common input between both hands. In addition, we integrated IMC across several frequency bands of interest to compare IMC content between blocks of unilateral contractions but again were unable to detect differences in IMC content between different blocks. All findings and their implications are discussed below.
As expected, amplitudes of FDIVol were significantly higher compared to FDIpMA. This is a robust finding as pMA amplitude, if evident, is usually only a fraction of the amplitude of the voluntarily contracting hand4,49. One aspect of our approach to improving the understanding of pMA was to study frequency features of pMA. Among the most commonly analyzed frequency features are MNF, MDF, and PKF. In this study, all three features were significantly increased from block 4 (MNF), block 5 (MDF), and block 6 (PKF) onwards compared to block 1 in both FDIVol and FDIpMA. Spectral features are subject to influence by a number of factors. One such factor is fatigue. Many studies report a decrease in MDF and MNF due to muscular fatigue, especially in isometric contractions25. This is attributed to the fact that peripheral fatigue leads to a downward shift in the frequency spectrum of sEMG signals due to slower muscle fiber conduction velocity, which results in a compression of the power spectrum and subsequently lower MDF and MNF values22. Nevertheless, it is worth noting that some studies show a broad range of MDF during muscular fatigue protocols, which somewhat contextualizes the strength of the association between MDF and muscular fatigue50,51. Importantly, the behavioral task employed here was designed to counteract potential fatigue effects by providing numerous rest periods between individual contractions and also between individual contraction blocks. Prolonged rest periods have been shown to be effective in preventing fatigue effects on sEMG signals during repeated contractions52. Thus, the increase in MDF and MNF observed in this study suggests that peripheral fatigue did not have a decisive effect on either parameter. Another factor that may modulate MDF and MNF is the force level which can be determined by MVC monitoring. However, the evidence is equivocal. Where some studies reported an increase of both parameters with increasing MVC50,51, others observed a decrease53 or no modulation between different levels of MVC54. Factors such as gender53, different muscle fiber compositions55, and electrode configurations25 have been cited as possible reasons for these varying results. An interesting finding in this context comes from Bilodeau et al.56 who were able to show that MDF and MNF increase in some muscles while they decrease in other muscles during contractions with increasing MVC. This implies a muscle-specific modulation of MDF and MNF or at least a considerable variability of these measures between different muscles. Similar findings demonstrate inconsistent patterns of PKF during force level modulations57. A further potential factor influencing MDF and MNF is the degree of motor recruitment. Generally, it is assumed that MVC-normalized sEMG signals allow conclusions to be drawn about motor recruitment, although this relationship has not been clearly established. Still, several studies show moderate58 to weak59 associations between sEMG amplitudes and the recruitment of motor units. The degree of motor unit recruitment may also influence frequency features. For example, a relationship between a linear increase in MDF and orderly recruitment of motor units has been demonstrated in animals60. Still, this relationship shows high inter-subject variability in humans59. Interestingly, although participants in this study were instructed to contract at 80% of their maximum force, our results show a gradual increase from roughly 60% MVC to 80% MVC from block 1 to block 9 in FDIVol and a comparable increase, albeit at lower force levels, for FDIpMA (Fig. 1). It is tempting to speculate that the observed increase in MDF and MNF is partly influenced by the increase in muscular activation, i.e., neural drive, as quantified by an increase in MVC in both FDIVol and FDIpMA. Thus, the increase in MDF and MNF observed in FDIVol and FDIpMA might be explained by progressive recruitment of new motor units with higher muscle fiber conduction velocities as a result of increasing muscle force60. Analysis of FR revealed significantly higher FR in FDIpMA compared to FDIVol as well as a significant decrease in FR from block 7 to block 9. In this study, a higher FR implies a higher proportion of low-frequency (20â40 Hz) signal components. This decrease in FR, similar to the observed changes in MDF, MNF, and PKF, can potentially be associated with increased recruitment of motor units. With increasing recruitment of motor units, the ratio between low- and high-frequency components of the EMG spectrum tilts in favor of the high-frequency components60,61. The fact that we observed an increased FR in FDIpMA predominantly in early blocks (Figs. 2, 3D) followed by an approximation of FR between both hands in later blocks may imply an increase in the spillover between FDIVol and FDIpMA related to the progressive recruitment of motor units in the active hand. Naturally, further evidence is needed to support this assumption. Surface EMG decomposition methods show promising results59 and should be considered in future studies investigating pMA mechanisms.
Regarding IMC, we observed two main findings. On an individual level, several participants showed significant IMC spanning from alpha to high gamma frequency bands. This finding is consistent with a previous study evaluating IMC between homologous hand muscles during rhythmic contractions. The authors found coherence between homologous hand muscles across participants although average patterns and their frequency bands were not analyzed62. Coherence was in part attributed to so-called neural cross-talk, i.e., patterns of interference between limb movements due to neural interactions between hemispheres, a process mainly mediated through the corpus callosum63. Indeed, studies on pMA in hand movements confirm an association between transcallosal fibers connecting bilateral primary motor cortices (M1) and pMA4. Specifically, decreased interhemispheric inhibition (IHI) from contralateral to ipsilateral primary motor cortices appears to be associated with increased pMA64, suggesting a governing role for the contralateral M1 in the suppression of pMA. However, inhibition appears to decrease with increasing force demands, as the ipsilateral M1 is progressively recruited in the process65. Isometric contractions above 30% of MVC have been shown to decrease IHI from contralateral to ipsilateral M1, which in turn allows an increase in corticospinal excitability66. In conjunction with our results, it can be speculated that the high contraction force (up to 80% MVC) caused a reduction in the inhibitory influence of the contralateral M1 on the ipsilateral M1, both inducing prominent pMA and serving as an explanation for the common input between the FDIVol and FDIpMA observed at the participant level. Critically, the patterns of IMC were highly variable between participants and, therefore, did not allow precise conclusions to be drawn about potential neuronal commonalities. Variability was evident both in terms of the blocks and the frequency bands within which IMC occurred. Blockwise grand-averaged IMC did not yield significant results, further emphasizing that precise patterns of common input between FDIVol and FDIpMA were not detected in this study. In general, these results do not equate to the absence of a common neural origin, but rather imply that the spectral correlates of this common input are highly individual. Future studies should pursue approaches in which pMA is maximized, for example, through fatigue protocols to increase the likelihood of robust patterns of shared input between active and resting muscles. Future studies should extend our results by corticomuscular coherence (CMC) analyses between electroencephalography and EMG of the active and resting hand67. CMC is particularly useful for inferring cortical origins in electrophysiological signal pairs and is therefore particularly useful for revealing further evidence regarding the neural origin of physiological mirror activity.
Limitations
A limitation of the present study is the selection of an all-male participant group. This was done to prevent variance issues due to possible sex differences in central-nervous structure and function30 as well as initial findings suggesting sex differences in pMA68. For this reason, our findings cannot be applied to female populations in general. Another limitation of our results relates to the age range of the participants. It is known that pMA follows an inverted U distribution as a function of age68. For this reason, our results need to be extended to other age ranges. A further influence on pMA results from the sporting activity of the participants. Previous research demonstrated that athletes who perform chronic, organized, high-volume training have higher pMA compared to non-athletes5. In this study, we selected participants without organized training backgrounds and low training volumes, corresponding to contemporary health guidelines. It will be interesting to investigate the influence of organized athletic activity at different performance levels on the magnitude as well as neural characteristics of pMA. Furthermore, we investigated only one force level (80% of maximum force) in this study. This was done in light of the positive correlation between pMA amplitude and force level. Nevertheless, future studies should compare different force levels to gain insight into modulatory distinctions of neural input between FDIVol and FDIpMA relating to force requirements. Ultimately, the present results apply to a single muscle only. Although the FDI is one of the most frequently studied muscles concerning pMA, further studies on neural mechanisms of pMA should also be carried out in other muscles, e.g., proximal muscles of the upper extremity as well as proximal and distal muscles of the lower extremity. This will allow more complete and robust conclusions to be drawn about the characteristics of pMA and its neural origins.
Conclusion
In summary, our results imply an overlap of spectral properties between active and passive hands during repeated unilateral contractions. This relationship may result from a reduction in interhemispheric inhibition resulting from high force levels as well as progressive recruitment of motor units. However, further studies are needed to support these assumptions. The investigation of pMA holds implications not only for clinical applications but also for understanding bimanual motor control. Our results provide an initial reference point for further investigations. Such studies may extend the inference of muscle activity from global (sEMG) to local (motor unit domain) parameters to understand the underlying processes more precisely. In this sense, the use of EMG decomposition methods is particularly promising.
Data availability
Data, in an anonymous format (according to data protection policy in the ethics agreement), is available at https://doi.org/10.6084/m9.figshare.19328747.v1.
References
Zijdewind, I. & Kernell, D. Bilateral interactions during contractions of intrinsic hand muscles. J. Neurophysiol. 85, 1907â1913. https://doi.org/10.1152/jn.2001.85.5.1907 (2001).
van Duinen, H., Renken, R., Maurits, N. M. & Zijdewind, I. Relation between muscle and brain activity during isometric contractions of the first dorsal interosseus muscle. Hum. Brain Mapp. 29, 281â299. https://doi.org/10.1002/hbm.20388 (2008).
Sehm, B., Perez, M. A., Xu, B., Hidler, J. & Cohen, L. G. Functional neuroanatomy of mirroring during a unimanual force generation task. Cereb. Cortex 20, 34â45. https://doi.org/10.1093/cercor/bhp075 (2010).
Sehm, B., Steele, C. J., Villringer, A. & Ragert, P. Mirror motor activity during right-hand contractions and its relation to white matter in the posterior midbody of the corpus callosum. Cereb. Cortex 26, 4347â4355. https://doi.org/10.1093/cercor/bhv217 (2016).
Maudrich, T. et al. Mirror electromyografic activity in the upper and lower extremity: A comparison between endurance athletes and non-athletes. Front. Hum. Neurosci. 11, 485. https://doi.org/10.3389/fnhum.2017.00485 (2017).
Todor, J. I. & Lazarus, J. A. Exertion level and the intensity of associated movements. Dev. Med. Child Neurol. 28, 205â212. https://doi.org/10.1111/j.1469-8749.1986.tb03856.x (1986).
Hopf, H. C., Schlegel, H. J. & Lowitzsch, K. Irradiation of voluntary activity to the contralateral side in movements of normal subjects and patients with central motor disturbances. Eur. Neurol. 12, 142â147 (1974).
Maudrich, T., Kenville, R., Lepsien, J., Villringer, A. & Ragert, P. Structural neural correlates of physiological mirror activity during isometric contractions of non-dominant hand muscles. Sci. Rep. 8, 9178. https://doi.org/10.1038/s41598-018-27471-5 (2018).
Uttner, I. et al. Mirror movements and the role of handedness: Isometric grip forces changes. Mot. Control 11, 16â28 (2007).
Tisseyre, J., Amarantini, D. & Tallet, J. Behavioural and cerebral asymmetries of mirror movements are specific to rhythmic task and related to higher attentional and executive control. Behav. Brain Res. 412, 113429. https://doi.org/10.1016/j.bbr.2021.113429 (2021).
Addamo, P. K., Farrow, M., Hoy, K. E., Bradshaw, J. L. & Georgiou-Karistianis, N. The influence of task characteristics on younger and older adult motor overflow. Q. J. Exp. Psychol. 62, 239â247. https://doi.org/10.1080/17470210802269217 (2009).
Tisseyre, J. et al. Mirror movements are linked to executive control in healthy and brain-injured adults. Neuroscience 379, 246â256. https://doi.org/10.1016/j.neuroscience.2018.03.027 (2018).
Tisseyre, J., Marquet-Doleac, J., Barral, J., Amarantini, D. & Tallet, J. Lateralized inhibition of symmetric contractions is associated with motor, attentional and executive processes. Behav. Brain Res. 361, 65â73. https://doi.org/10.1016/j.bbr.2018.12.034 (2019).
Liederman, J. & Foley, L. M. A modified finger lift test reveals an asymmetry of motor overflow in adults. J. Clin. Exp. Neuropsychol. 9, 498â510. https://doi.org/10.1080/01688638708410765 (1987).
Cincotta, M. & Ziemann, U. Neurophysiology of unimanual motor control and mirror movements. Clin. Neurophysiol. 119, 744â762. https://doi.org/10.1016/j.clinph.2007.11.047 (2008).
Post, M. et al. Reduced cortical activity during maximal bilateral contractions of the index finger. Neuroimage 35, 16â27. https://doi.org/10.1016/j.neuroimage.2006.11.050 (2007).
Post, M., Bayrak, S., Kernell, D. & Zijdewind, I. Contralateral muscle activity and fatigue in the human first dorsal interosseous muscle. J. Appl. Physiol. 105, 70â82. https://doi.org/10.1152/japplphysiol.01298.2007 (2008).
Maudrich, T. et al. Inverse relationship between amplitude and latency of physiological mirror activity during repetitive isometric contractions. Neuroscience 406, 300â313. https://doi.org/10.1016/j.neuroscience.2019.03.029 (2019).
Perez, M. A. & Cohen, L. G. Mechanisms underlying functional changes in the primary motor cortex ipsilateral to an active hand. J. Neurosci. 28, 5631â5640. https://doi.org/10.1523/jneurosci.0093-08.2008 (2008).
Yensen, R. A factor influencing motor overflow. Percept. Mot. Skills 20, 967â968. https://doi.org/10.2466/pms.1965.20.3.967 (1965).
Farina, D., Merletti, R. & Enoka, R. M. The extraction of neural strategies from the surface EMG: An update. J. Appl. Physiol. 117, 1215â1230 (2014).
McManus, L., De Vito, G. & Lowery, M. M. Analysis and biophysics of surface EMG for physiotherapists and kinesiologists: Toward a common language with rehabilitation engineers. Front. Neurol. https://doi.org/10.3389/fneur.2020.576729 (2020).
Amorim, C. F. & Marson, R. A. Application of Surface Electromyography in the Dynamics of Human Movement Vol. 10 (Acesso em, 2012).
Campanini, I., Cosma, M., Manca, M. & Merlo, A. Added value of dynamic EMG in the assessment of the equinus and the equinovarus foot deviation in stroke patients and barriers limiting its usage. Front. Neurol. 11, 583399 (2020).
Phinyomark, A., Thongpanja, S., Hu, H., Phukpattaranont, P. & Limsakul, C. The usefulness of mean and median frequencies in electromyography analysis. In Computational Intelligence in Electromyography AnalysisâA Perspective on Current Applications and Future Challenges 195â220 (2012).
Laine, C. M. & Valero-Cuevas, F. J. Intermuscular coherence reflects functional coordination. J. Neurophysiol. 118, 1775â1783. https://doi.org/10.1152/jn.00204.2017 (2017).
Farmer, S. F., Bremner, F. D., Halliday, D. M., Rosenberg, J. R. & Stephens, J. A. The frequency content of common synaptic inputs to motoneurones studied during voluntary isometric contraction in man. J. Physiol. 470, 127â155. https://doi.org/10.1113/jphysiol.1993.sp019851 (1993).
Negro, F. & Farina, D. Linear transmission of cortical oscillations to the neural drive to muscles is mediated by common projections to populations of motoneurons in humans. J. Physiol. 589, 629â637. https://doi.org/10.1113/jphysiol.2010.202473 (2011).
Negro, F. & Farina, D. Decorrelation of cortical inputs and motoneuron output. J. Neurophysiol. 106, 2688â2697. https://doi.org/10.1152/jn.00336.2011 (2011).
Grabowska, A. Sex on the brain: Are gender-dependent structural and functional differences associated with behavior?. J. Neurosci. Res. 95, 200â212. https://doi.org/10.1002/jnr.23953 (2017).
Oldfield, R. C. The assessment and analysis of handedness: The Edinburgh inventory. Neuropsychologia 9, 97â113 (1971).
Pesta, D. H., Angadi, S. S., Burtscher, M. & Roberts, C. K. The effects of caffeine, nicotine, ethanol, and tetrahydrocannabinol on exercise performance. Nutr. Metab. (Lond.) 10, 71. https://doi.org/10.1186/1743-7075-10-71 (2013).
Maudrich, T. et al. Voluntary inhibition of physiological mirror activity: An EEG-EMG study. Eneuro. https://doi.org/10.1523/eneuro.0326-20.2020 (2020).
De Luca, G. Fundamental concepts in EMG signal acquisition. (Delsys Inc, 2003).
Merletti, R. & Di Torino, P. Standards for reporting EMG data. J. Electromyogr. Kinesiol. 9, 3â4 (1999).
Welch, P. The use of fast Fourier transform for the estimation of power spectra: A method based on time averaging over short, modified periodograms. IEEE Trans. Audio Electroacoust. 15, 70â73 (1967).
Allison, G. & Fujiwara, T. The relationship between EMG median frequency and low frequency band amplitude changes at different levels of muscle capacity. Clin. Biomech. 17, 464â469 (2002).
Nazmi, N. et al. A review of classification techniques of EMG signals during isotonic and isometric contractions. Sensors 16, 1304 (2016).
Han, J.-S. et al. In Proc. of 6th international conference on soft computing (IIZUKA2000) 890â897.
Nojima, I., Watanabe, T., Saito, K., Tanabe, S. & Kanazawa, H. Modulation of EMGâEMG coherence in a choice stepping task. Front. Hum. Neurosci. 12, 50 (2018).
Boonstra, T. W. et al. Muscle networks: Connectivity analysis of EMG activity during postural control. Sci. Rep. 5, 17830. https://doi.org/10.1038/srep17830 (2015).
Kenville, R. et al. Intermuscular coherence between homologous muscles during dynamic and static movement periods of bipedal squatting. J. Neurophysiol. 124, 1045â1055. https://doi.org/10.1152/jn.00231.2020 (2020).
Rosenberg, J. R., Amjad, A. M., Breeze, P., Brillinger, D. R. & Halliday, D. M. The Fourier approach to the identification of functional coupling between neuronal spike trains. Prog. Biophys. Mol. Biol. 53, 1â31 (1989).
Conway, B. et al. Synchronization between motor cortex and spinal motoneuronal pool during the performance of a maintained motor task in man. J. Physiol. 489, 917â924 (1995).
Farmer, S. Rhythmicity, synchronization and binding in human and primate motor systems. J. Physiol. 509, 3â14 (1998).
Jaiser, S. R., Baker, M. R. & Baker, S. N. Intermuscular coherence in normal adults: Variability and changes with age. PLoSÂ One 11, e0149029. https://doi.org/10.1371/journal.pone.0149029 (2016).
Spedden, M. E. et al. The development of functional and directed corticomuscular connectivity during tonic ankle muscle contraction across childhood and adolescence. Neuroimage 191, 350â360. https://doi.org/10.1016/j.neuroimage.2019.02.054 (2019).
Omlor, W., Patino, L., Hepp-Reymond, M.-C. & Kristeva, R. Gamma-range corticomuscular coherence during dynamic force output. Neuroimage 34, 1191â1198 (2007).
Zijdewind, I., Butler, J. E., Gandevia, S. C. & Taylor, J. L. The origin of activity in the biceps brachii muscle during voluntary contractions of the contralateral elbow flexor muscles. Exp. Brain Res. 175, 526â535 (2006).
Doheny, E. P., Lowery, M. M., FitzPatrick, D. P. & OâMalley, M. J. Effect of elbow joint angle on forceâEMG relationships in human elbow flexor and extensor muscles. J. Electromyogr. Kinesiol. 18, 760â770. https://doi.org/10.1016/j.jelekin.2007.03.006 (2008).
Gerdle, B., Eriksson, N. E. & Brundin, L. The behaviour of the mean power frequency of the surface electromyogram in biceps brachii with increasing force and during fatigue. With special regard to the electrode distance. Electromyogr. Clin. Neurophysiol. 30, 483â489 (1990).
Ortega-Auriol, P. A., Besier, T. F., Byblow, W. D. & McMorland, A. J. C. Fatigue influences the recruitment, but not structure, of muscle synergies. Front. Hum. Neurosci. https://doi.org/10.3389/fnhum.2018.00217 (2018).
Kaplanis, P. A., Pattichis, C. S., Hadjileontiadis, L. J. & Roberts, V. C. Surface EMG analysis on normal subjects based on isometric voluntary contraction. J. Electromyogr. Kinesiol. 19, 157â171. https://doi.org/10.1016/j.jelekin.2007.03.010 (2009).
Cechetto, A. D., Parker, P. A. & Scott, R. N. The effects of four time-varying factors on the mean frequency of a myoelectric signal. J. Electromyogr. Kinesiol. 11, 347â354. https://doi.org/10.1016/s1050-6411(01)00010-4 (2001).
Farina, D., Fosci, M. & Merletti, R. Motor unit recruitment strategies investigated by surface EMG variables. J. Appl. Physiol. 1985(92), 235â247. https://doi.org/10.1152/jappl.2002.92.1.235 (2002).
Bilodeau, M., Arsenault, A. B., Gravel, D. & Bourbonnais, D. EMG power spectra of elbow extensors during ramp and step isometric contractions. Eur. J. Appl. Physiol. Occup. Physiol. 63, 24â28. https://doi.org/10.1007/bf00760796 (1991).
Gander, R. E. & Hudgins, B. S. Power spectral density of the surface myoelectric signal of the biceps brachii as a function of static load. Electromyogr. Clin. Neurophysiol. 25, 469â478 (1985).
Suzuki, H., Conwit, R. A., Stashuk, D., Santarsiero, L. & Metter, E. J. Relationships between surface-detected EMG signals and motor unit activation. Med. Sci. Sports Exerc. 34, 1509â1517. https://doi.org/10.1097/00005768-200209000-00018 (2002).
Del Vecchio, A., Negro, F., Felici, F. & Farina, D. Associations between motor unit action potential parameters and surface EMG features. J. Appl. Physiol. 1985(123), 835â843. https://doi.org/10.1152/japplphysiol.00482.2017 (2017).
Solomonow, M. et al. Electromyogram power spectra frequencies associated with motor unit recruitment strategies. J. Appl. Physiol. 1985(68), 1177â1185. https://doi.org/10.1152/jappl.1990.68.3.1177 (1990).
Broman, H., Bilotto, G. & De Luca, C. J. Myoelectric signal conduction velocity and spectral parameters: Influence of force and time. J. Appl. Physiol. 1985(58), 1428â1437. https://doi.org/10.1152/jappl.1985.58.5.1428 (1985).
Ridderikhoff, A., Daffertshofer, A., Peper, C. E. & Beek, P. J. Mirrored EMG activity during unimanual rhythmic movements. Neurosci. Lett. 381, 228â233. https://doi.org/10.1016/j.neulet.2005.02.041 (2005).
Swinnen, S. P. Intermanual coordination: From behavioural principles to neural-network interactions. Nat. Rev. Neurosci. 3, 348â359. https://doi.org/10.1038/nrn807 (2002).
Hübers, A., Orekhov, Y. & Ziemann, U. Interhemispheric motor inhibition: Its role in controlling electromyographic mirror activity. Eur. J. Neurosci. 28, 364â371. https://doi.org/10.1111/j.1460-9568.2008.06335.x (2008).
Uehara, K. & Funase, K. Contribution of ipsilateral primary motor cortex activity to the execution of voluntary movements in humans: A review of recent studies. J. Phys. Fit. Sports Med. 3, 297â306 (2014).
Cabibel, V., Muthalib, M., Teo, W. P. & Perrey, S. High-definition transcranial direct-current stimulation of the right M1 further facilitates left M1 excitability during crossed facilitation. J. Neurophysiol. 119, 1266â1272. https://doi.org/10.1152/jn.00861.2017 (2018).
Mima, T. & Hallett, M. Corticomuscular coherence: A review. J. Clin. Neurophysiol. 16, 501â511 (1999).
Koerte, I. et al. Mirror movements in healthy humans across the lifespan: Effects of development and ageing. Dev. Med. Child Neurol. 52, 1106â1112. https://doi.org/10.1111/j.1469-8749.2010.03766.x (2010).
Funding
Open Access funding enabled and organized by Projekt DEAL. This work was supported by the Max Planck Society.
Author information
Authors and Affiliations
Contributions
R.K. and T.M. designed the study, acquired the data, analyzed the data, and wrote the manuscript. All authors interpreted the data, contributed to the manuscript, reviewed it, approved the content of the final version, and agree to be accountable for all aspects of the work. All persons designated as authors qualify for authorship, and all those who qualify for authorship are listed.
Corresponding authors
Ethics declarations
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.
Supplementary Information
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/.
About this article
Cite this article
Kenville, R., Maudrich, T. Spectral properties of physiological mirror activity: an investigation of frequency features and common input between homologous muscles. Sci Rep 12, 15965 (2022). https://doi.org/10.1038/s41598-022-20413-2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-022-20413-2