Next Article in Journal
Chinese Digital Arm (CDA): A High-Precision Digital Arm for Electrical Stimulation Simulation
Previous Article in Journal
Connectivity Analysis in EEG Data: A Tutorial Review of the State of the Art and Emerging Trends
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ankle Joint Angle Influences Relative Torque Fluctuation during Isometric Plantar Flexion

1
Department of Biomedical Engineering, McCormick School of Engineering, Northwestern University, Evanston, IL 60208, USA
2
Shirley Ryan AbilityLab (Formerly the Rehabilitation Institute of Chicago), Chicago, IL 60611, USA
3
Department of Physical Medicine and Rehabilitation, Feinberg School of Medicine, Northwestern University, Chicago, IL 60611, USA
4
Department of Biomedical Engineering, Newark College of Engineering, New Jersey Institute of Technology, Newark, NJ 07102, USA
*
Author to whom correspondence should be addressed.
Bioengineering 2023, 10(3), 373; https://doi.org/10.3390/bioengineering10030373
Submission received: 29 January 2023 / Revised: 9 March 2023 / Accepted: 13 March 2023 / Published: 18 March 2023

Abstract

:
The purpose of this study was to investigate the influence of changes in muscle length on the torque fluctuations and on related oscillations in muscle activity during voluntary isometric contractions of ankle plantar flexor muscles. Eleven healthy individuals were asked to perform voluntary isometric contractions of ankle muscles at five different contraction intensities from 10% to 70% of maximum voluntary isometric contraction (MVIC) and at three different muscle lengths, implemented by changing the ankle joint angle (plantar flexion of 26°-shorter muscle length; plantar flexion of 10°-neutral muscle length; dorsiflexion of 3°-longer muscle length). Surface electromyogram (EMG) signals were recorded from the skin surface over the triceps surae muscles, and rectified-and-smoothed EMG (rsEMG) were estimated to assess the oscillations in muscle activity. The absolute torque fluctuations (quantified by the standard deviation) were significantly higher during moderate-to-high contractions at the longer muscle length. Absolute torque fluctuations were found to be a linear function of torque output regardless of muscle length. In contrast, the relative torque fluctuations (quantified by the coefficient of variation) were higher at the shorter muscle length. However, both absolute and relative oscillations in muscle activities remained relatively consistent at different ankle joint angles for all plantar flexors. These findings suggest that the torque steadiness may be affected by not only muscle activities, but also by muscle length-dependent mechanical properties. This study provides more insights that muscle mechanics should be considered when explaining the steadiness in force output.

1. Introduction

Force (or torque) fluctuations in the low-frequency (<1 Hz) band are well-known to arise during voluntary sustained isometric contractions in human skeletal muscles [1]. Such fluctuations, often quantified by standard deviation (SD) or coefficient of variation (CoV), appear to be a useful biomarker to help us understand motor performance. For example, greater force fluctuations have been observed in individuals in aging [2,3,4], musculoskeletal injuries [5,6] and neurophysiology disorders [3], demonstrating significant correlation with motor function [1,7]. Accordingly, a systematic investigation into underlying mechanisms of force fluctuations is essential for understanding motor activities and for improving motor control ability in broad clinical populations.
Force fluctuations were found to be closely related to variability in the recruitment and rate-coding properties of the motor units [8]. Both experimental and simulation studies demonstrated a linear increase of SD for force (FSD) with increasing voluntary isometric contraction level [9]. Such increase in FSD can be explained by the presence of signal-dependent noise that appears to be associated with the statistical distribution of the motoneuron spike trains [10]. Previous studies have also suggested that CoV for force (FCoV) decreases non-linearly as the contraction level increases by up to approximately 20% of maximum voluntary isometric contraction (MVIC) [9]. Such reduction in FCoV can be explained, in part, by the characteristics of motor unit recruitment and motor unit discharge variability [7,11]. For instance, significant correlations between FCoV and CoV for motor unit discharge rate were found in younger individuals [7]. Collectively, low-frequency oscillations in force may be consequences of variability in a single motor unit’s firing rate or low-frequency oscillations in the neural drive to the muscle, commonly estimated by the cumulative spike trains or the rectified-and-smoothed EMG (rsEMG) [12,13,14].
Until now, many studies have suggested the role and significance of neurophysiological factors in the generation of force fluctuation but lacked the materials to address how muscle mechanics can modulate the force fluctuation. Considering muscle length-mediated changes in contractile properties (e.g., time-to-peak, half-relaxation time, or force-length characteristics) [15] as well as in afferent feedback [16], it is likely that muscle mechanics can also impact the steadiness of the force output. A few studies using the joint angle as a proxy for muscle length showed that shorter muscle length led to an increase in the plantar flexion torque variability, when compared with longer muscle length at the same force level [17,18]. This may be due to the finding that the twitch time of the recruited motor units are shortened at shorter muscle lengths [15], and more motor units are recruited [19]. Both factors may induce a decreased firing rate or larger variability in firing rate and thus lead to higher force fluctuation. Moreover, any changes in muscle length can result in changes in spinal cord excitability. At a more flexed ankle joint angle, there is less stretch on muscle spindles (Ia afferent) and thus reduced excitatory drive to the motoneuron pool, potentially affecting the steadiness of motor output [20].
The purpose of this study was to evaluate the influence of changes in muscle length on the plantar flexion torque variability and oscillations in muscle activities during submaximal voluntary isometric contractions in neurologically intact individuals. To assess the oscillations in muscle activities, we collected the surface EMG signals and obtained the SD for the rsEMG of triceps surae muscles. We hypothesized that torque variability will be affected by muscle length, presumably related to underlying changes in low-frequency oscillations in neural drive and/or muscle mechanics.

2. Materials and Methods

2.1. Participants

Eleven young healthy individuals (Age: 27 ± 1.5 yrs.; F/M: 5/6) participated in this study after providing consent. All procedures were approved by the Institutional Review Board in Northwestern University and complied with the Helsinki declaration.

2.2. Experimental Setup

Participants were asked to sit upright in a Biodex chair with the knee fully extended (Figure 1). The right foot was fixed onto a footplate attached to a 6-axis loadcell (Omega160, ATI Industrial Automation, Apex, NC, USA), secured via an adhesive foot strap. Single differential surface EMG electrodes (Bagnoli, Delsys Inc., Boston, MA, USA) were placed on the triceps surae muscles (i.e., middle part of the lateral gastrocnemius (LG), soleus (SOL), and tibialis anterior, and on the distal part of the MG muscle belly) with a ground electrode attached on the bony landmark (i.e., lateral malleolus). The locations of electrodes were determined after muscle palpations on the muscles of interest. Prior to electrode placements, the sites were cleaned with alcohol pads. EMG and plantar flexion torque signals were performed simultaneously at a sampling frequency of 2 kHz (NI USB-6259 BNC, National Instrument, Austin, TX, USA).

2.3. Experimental Protocol

Plantar flexion contractions with the right leg were performed at three ankle joint angles: ankle joint angles at 87°, at 100°, and at 116° (the angle between shank and sole of foot), corresponding to longer, neutral, and shorter muscle lengths, respectively. Each configuration was set by adjusting the foot-plate position. At each ankle joint angle, three MVICs were performed, and the averaged value of these was used to estimate the target submaximal contraction torque. Each subject was asked to perform isometric contractions at 10%, 20%, 30%, 50%, and 70% MVIC. Visual feedback was provided in real time involving reaching the target level from rest state (5–10 s) and then maintaining at the target level ± 3% MVIC for at least 6 s. Each contraction level was repeated three times. To minimize the effects from systematic muscle fatigue, the orders of the ankle joint angles and target contraction levels were randomized, and at least a 30 s break was provided between trials.

2.4. Data Analysis

The torque signals were low-pass filtered using a fourth-order Butterworth filter with a cutoff frequency at 6 Hz. Raw EMG signals were band-pass filtered using a fourth-order Butterworth filter with a passband of 20–450 Hz. The filtered EMG signals were then rectified and smoothed by applying a fourth-order Butterworth filter with a cutoff frequency of 2 Hz (hereafter called rsEMG).
To determine a representative stable measure of the variability of muscle contractions for each trial, a 4 s segment in the middle of each sustained contraction was chosen by calculating the minimum standard deviation in the filtered torque signals. Using this time segment, both absolute and relative variability for the filtered torque signals were quantified as SD and as the CoV, respectively. Hereafter, absolute and relative variability for torque is called TSD and TCoV. CoV was defined as the ratio of SD to the mean torque value of the segment. Both absolute and relative variability values for the rsEMG of each plantar flexor (MGSD, LGSD, and SOLSD; MGCoV, LGCoV, and SOLCoV) were calculated to examine the oscillations in muscle activities. Figure 2 shows representative detrended torque trials and rsEMG signals from a representative subject.

2.5. Statistical Analysis

A Kolmogorov–Smirnov test was performed to assess the normality of the data for each test. Since the normality test on the major outcomes rejected the null hypothesis, a log transformation was applied. Given the evaluation of their skewness and kurtosis after log transformation, the data distributions were acceptable to conduct recommended parametric tests [21]. Thus, a two-way, repeated-measures ANOVA was performed to examine the differences in the log-transformed torque measurement outcomes (TSD and TCoV) and oscillations of muscle activities at the five contraction levels and three ankle joint angles via the SPSS (IBM Corp., Armonk, NY, USA) with a significance level (α) of 0.05. When necessary, post hoc pairwise multiple comparisons with Bonferroni correction were used. In any case the sphericity assumption was violated, Greenhouse–Geisser would be used instead.
Stepwise regression analysis was applied to test what is the more dominant contributor to the changes in TCoV with contraction levels and ankle joint angles. To determine the dominant muscle that contributes to the changes in torque variability with contraction levels, the data were collapsed across the contraction levels, and the stepwise regression analysis was conducted at each ankle joint angle. When determining the dominant muscle that contributes to the changes in torque variability with the ankle joint angle at each contraction level, the data were collapsed across the ankle joint angles, and the analysis was conducted at each contraction level.

3. Results

3.1. Absolute Variability for Torque

The ANOVA analysis revealed significant main effects from both the contraction level (F (4, 40) = 142.458, p < 0.001) and ankle angle (F (2, 20) = 6.547, p = 0.006) on TSD. No significant interaction was also found between the contraction level and ankle angle (F (8, 80) = 1.325, p = 0.243). As shown in Figure 3a, further analysis suggested that TSD at the longer muscle length was observed to be significantly larger than that at the shorter muscle length when the contraction level is at 30% (dz = 0.274, p = 0.026), 50% (dz = 0.420, p = 0.029) and 70% MVIC (dz = 0.414, p = 0.006). As shown in Figure 3b, TSD increases linearly with the actual plantar flexion torque regardless of muscle length.

3.2. Absolute Variability for rsEMG

The average MGSD, LGSD, and SOLSD at different contraction levels and ankle angles are presented in Table 1. Significant effects of contraction level are indicated for all three muscles (MG: F (4, 40) = 240.271, p < 0.001; LG: F (4, 40) = 145.711, p < 0.001; and SOL: F (4, 40) = 72.440, p < 0.001). For the MG muscle, the ANOVA results indicated that there is no significant effect of the ankle angle (F (2, 20) = 1.240, p = 0.311) or no interaction (F (8, 80) = 0.589, p = 0.784). Similarly, LGSD was not affected by the ankle angle (F (2, 20) = 0.612, p = 0.552) or the interaction (F (8, 80) = 0.598, p = 0.777). For the SOL muscle, there are no significant effects from both the ankle angle (F (2, 20) = 1.168, p = 0.331) and the interaction (F (8, 80) = 1.101, p = 0.372).

3.3. Relative Variability for Torque

The statistical analysis suggested a significant main effect from the contraction level (F (4, 40) = 8.037, p < 0.001) and ankle joint angle (F (2, 20) = 11.177, p < 0.001) on TCoV, but there was no significant effect from the interaction between the contraction level and ankle joint angle (F (8, 80) = 1.536, p = 0.158). As shown in Figure 4a, TCoV at the ankle joint angle of 116° was significantly higher than that at the other two ankle joint angles at the contraction of 10% MVIC (Shorter vs. Longer: dz = 0.724, p = 0.031; Shorter vs. Neutral: dz = 0.743, p = 0.039) and 20% MVIC (Shorter vs. Longer: dz = 0.734, p < 0.001; Shorter vs. Neutral: dz = 0.818, p < 0.001). As shown in Figure 4b, TCoV at ankle joint angle of 116° tends to be the highest at the same plantar flexion torque compared with the other two ankle joint angles.

3.4. Relative Variability for rsEMG

The average MGCoV, LGCoV, and SOLCoV at each ankle angle are presented in Table 2. For the MG muscle, the ANOVA results indicated that there is no significant effect from the ankle joint angle (F (2, 20) = 1.529, p = 0.241) and no significant interaction between the ankle joint angle and the contraction level (F (8, 80) = 0.463, p = 0.879). Although it shows a significant dependence on the contraction level (F (4, 40) = 5.311, p = 0.002), it only shows a significant reduction from 10% to 50% MVIC (p = 0.036) at the longest muscle length and then remained relatively constant among other comparisons. For the LG muscle, no significant dependence was observed from the contraction level (F (4, 40) = 0.822, p = 0.519), ankle joint angle (F (2, 20) = 0.572, p = 0.573), or their interaction (F (3.788, 37.872) = 1.065, p = 0.396). For the SOL muscle, there is no significant effect from the contraction level (F (4, 40) = 1.434, p = 0.240) and the interaction between the ankle joint angle and the contraction level (F (3.413, 34.125) = 0.419 p = 0.765). It shows a significant effect from the ankle joint angle (F (2, 20) = 5.915, p = 0.010) with a significant reduction in SOLCoV at 20% and 70% MVIC when the ankle angle is changed from 116° to 87°, and a reduction at 30% MVIC when it is changed from 100° to 87°.

3.5. Contributions from Relative Variability of rsEMG to the Resultant Relative Torque Variability

It appears that individual muscle contributions to TCoV may be different at different ankle joint angles (Table 3). The stepwise regression model revealed that the MGCoV contribution to TCoV was significant at all tested ankle angles. The LGCoV significantly contributed to TCoV at both the neutral and shorter muscle lengths, whereas the SOLCoV contribution was significant only at the longer muscle length.
When considering the individual muscle contributions to TCoV at different contraction intensities (Table 4), MGCoV seems to be a more dominant contributor compared to the other two muscles, as supported by its significant contribution to TCoV from 20 to 70% MVIC. The contribution of both LGCoV and SOLCoV to TCoV was significant only at 30 and 50% MVIC. Interestingly, TCoV at 10% MVIC was not explained significantly by any of triceps surae.

4. Discussion

The main findings of this study are as follows: (1) absolute torque variability (TSD) increased linearly with increasing the contraction level, and a significantly higher TSD was observed during moderate-to-high submaximal contractions (i.e., 30–70% MVIC) at the longer muscle length compared to the shorter muscle length, while the oscillations in muscle activities (absolute variability for rsEMG of each muscle—MGSD, LGSD, and SOLSD) remained relatively constant at different ankle joint angles; (2) TSD was found to be linearly increased with actual torque output regardless of ankle joint angles; (3) relative torque variability (TCoV) was found to be significantly higher during low submaximal contractions (i.e., 10–20% MVIC) at the shorter muscle length compared to the other two muscle lengths, but it appears that there is no significant effect of ankle joint angle on TCoV during moderate-to-high submaximal contractions; and (4) the stepwise regression models suggested that among triceps surae, the MG may play an important role to control the steadiness of plantar flexion torque.
TSD increased monotonically with contraction level at each ankle joint angle, as supported by previous findings [22]. This linear increment with contraction level of TSD may be explained by the presence of the signal-dependent noise (SDN) model, which assumes that there is noise from the motor command and that the amount of noise scales with the magnitude of motor command [10]. Our study further emphasized that the presence of SDN may also hold in a joint with multiple synergists. Studies have demonstrated that TSD is highly correlated with the low-frequency oscillations in the common neural drive, which can be quantified via cumulative spike trains or rsEMG of the target muscles [12,13,14]. This assertion can also be supported by our findings that MGSD, LGSD, and SOLSD increase monotonically with the contraction level (Table 1).
The significantly larger TSD was found at the longest muscle length during the moderate-to-high contractions (i.e., 30–70% MVIC). Since there are no significant effects from the ankle joint angle on the oscillations in muscle activities (Table 1), it suggests that the significant changes in TSD at the longer muscle length is not likely affected by oscillations in common neural drive; instead, TSD appears to be a linear function of actual torque output (Figure 3b), which further implies that the increment in TSD at the longer muscle during moderate-to-high contractions may result from the higher torque output linked to the length-tension function [23,24,25], while the lack of significant differences in TSD at lower contractions may be due to the fact that the actual torque output in that range across the muscle lengths is similar. This may suggest that absolute torque variability is a function of actual torque output regardless of muscle mechanics.
Our results revealed that at all ankle joint angles, TCoV decreased from a higher value and remained constant afterwards with significance detected at both the neutral (100°) and shorter muscle length (116°), in agreement with previous findings [9]. There are no clear underlying mechanisms for this yet, but it is plausible that such nonlinear relationships can be explained, in part, by motor unit mechanics (i.e., firing patterns and contractile properties), as supported by the previous finding, which demonstrated a significant correlation between CoV for motor unit discharge rate and TCoV at low contraction intensities of the first dorsal interosseous muscles in older individuals [7]. Given that the influence of variability in discharge rate of a single motor unit is almost attenuated by convolution of motoneuron spike trains with motor unit twitches and summation of twitch forces, motor unit discharge variability may affect the steadiness in force output only at low contraction levels. At higher force levels thereafter, the more fused contractions due to more active motor units and their increased firing rates would result in higher absolute but lower relative force variability. However, rsEMG has limitations to provide information regarding an individual MU firing characteristics. To better understand the association between the firing characteristics and force variability, further studies are needed.
TCoV at low contraction levels (i.e., 10–20% MVIC) was significantly higher at the shorter muscle length (Figure 4a). Indeed, to generate the comparable forces at the shorter muscle length, it is most likely that more active motor units (MUs) are recruited and/or higher firing rates of the active MUs are required [19], leading to a more fused status and thus a smaller TCoV. The potential underlying mechanisms may involve the changes in twitch properties of the activated MUs (i.e., decreased half-relaxation time) but with a relatively constant firing rate at the shorter muscle length [15], which may lead to less fused status and thus a higher TCoV. Another possible mechanism may be related to the inhibited transmission of the low-frequency oscillations in force signals at the shorter muscle length while the muscle is under high slack (i.e., 10–20% MVIC). The lack of significance in the changes of TCoV at moderate-to-high contractions (i.e., 30–70% MVIC in Figure 4a) may be due to the fact that muscle slack uptake occurs while the muscle is with high active tension and thus can limit the effects from changes in muscle length at moderate-to-high contractions. As demonstrated in Table 3, the lack of contribution from any of triceps surae on TCoV at 10% MVIC may further suggest that the relative torque variability at low contractions can be highly affected by muscle mechanics instead of muscle activation. Future studies are required to better understand the potential effects of length-dependent mechanisms on force variability, which have received less attention so far.
Some animal studies demonstrated direct evidence on selective recruitment of different types of motor units, particularly in synergistic muscles [26]. It is then possible that SOL may be recruited first and contribute predominantly to ankle torque generation, thus affecting the changes in torque variability at low contraction levels. However, different from animal models, MG is a mixed muscle containing at least 40% slow-twitch fibers in humans [27]. Our results (Table 3 and Table 4) suggested that MG muscle may play an important role to control the steadiness of plantar flexion torque among triceps surae. Moreover, TCoV at low contraction levels (i.e., 10% MVIC) among different ankle angles was not explained significantly by any of triceps surae, which may provide indirect evidence that torque variability can be highly affected by muscle mechanics instead of by muscle activation strategy at a low contraction level, and thus the potential effects from the recruitment order among triceps surae on the different torque variabilities at different ankle angles may be limited.
Several studies indicated that most of the variability in the force signal during the steady contractions can be explained by fluctuations in the common modulation of the motor unit discharge rate, with low-frequency oscillations over time, indicated by rsEMG [13,28]. It has been proposed to use force variability as a potential measurement approach to evaluate motor function in clinical populations (i.e., stroke survivors) [13,28]. Our study suggested that not only the potential changes in motor function, but also the changes in muscle mechanics, can impact the changes in force variability. Considering that muscle mechanics disorders (i.e., muscle stiffness) can also be highly impacted in individuals with neurological disorders, care should be taken when explaining the changes in force steadiness in clinical populations. Further studies would be needed to understand the relative contribution from changes in muscle mechanics and motor function to the steadiness of the resultant force output.
There are several limitations in our study. Previous research has demonstrated that tibialis anterior (TA) muscle contracts simultaneously [29,30] to maintain the joint stability during isometric contractions. However, our further evaluation on the activation status of the TA muscle shows that the RMS EMG of TA was significantly smaller when compared with MG and LG; moreover, its amplitude was less than 4% of the RMS EMG of TA collected during maximum dorsiflexion, suggesting that the effects from the antagonist on the plantar flexion torque variability may be negligible. We also have limitations on the manipulations of changes in the ankle joint angle due to the constraints of our setup. Our results only demonstrated a significant difference in TCoV between the shorter and the longer/neutral muscle length but lack an observable difference when comparing between the neutral and longer muscle length. This may be due to the small interval of changes in the ankle joint angles. Another limitation is that the targeted population in this study only involves young individuals. Considering that aging and neurological impairment (i.e., stroke survivors) also affect the steadiness of force output [2,3,31] as well as mechanical properties of the muscle–tendon unit [32], it would be interesting to investigate the relative contribution of changes in the mechanical properties to the force variability in different populations. Lastly, it is possible that the higher TCoV observed at the shorter muscle length during lower contractions (10–20% MVIC) can be affected by the smaller torque output at a shorter muscle length. However, considering that the average torque value during 20% MVIC at a longer muscle length (12.9 N m) is comparable to that during 30% MVIC at a shorter muscle length (12.1 N m), our further evaluation suggested that TCoV at the shorter muscle length was ~20% greater compared to the longer muscle length. TCoV–torque relation also illustrated that TCoV at the shorter muscle length is likely higher at a given, comparable torque output (Figure 4b). Future studies are necessary to examine torque fluctuations at the matched actual torque.

5. Conclusions

This study demonstrated that the absolute torque variability increased during modest submaximal contractions at a longer muscle length, whereas the relative torque variability increased during low submaximal contractions at a shorter muscle length. These findings suggest that the torque steadiness may be affected by both neural drive and muscle mechanics. Future studies are needed to better describe how these neuromuscular properties can influence the variability of force output.

Author Contributions

Conceptualization, F.S., W.Z.R. and J.S.; methodology, F.S., W.Z.R. and J.S.; software, F.S. and J.S.; validation, F.S., W.Z.R. and J.S.; formal analysis, F.S., W.Z.R. and J.S.; investigation, F.S., W.Z.R. and J.S.; resources, W.Z.R. and J.S.; data curation, F.S., W.Z.R. and J.S.; writing—original draft preparation, F.S. and J.S.; writing—review and editing, W.Z.R. and J.S.; visualization, F.S., W.Z.R. and J.S.; supervision, W.Z.R. and J.S.; project administration, F.S., W.Z.R. and J.S.; funding acquisition, W.Z.R. and J.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by grants from the National Institute on Disability, Independent Living, and Rehabilitation Research (90SFGE0005), the Davee Foundation Stroke Research Seed Grant initiative and the Northwestern University Department of Neurology, Division of Stroke and Neurocritical Care, and the National Institutes of Health (R01HD089952).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Institutional Review Board (or Ethics Committee) of Northwestern University (Approval code: STU00206371; Approval date: 02/07/2018).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

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

Acknowledgments

We thank all participants in this study. We also thank Andres Cardona for technical preparation. Andres is informed about this acknowledgement, and he agreed to be mentioned.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Enoka, R.M.; Farina, D. Force steadiness: From motor units to voluntary actions. Physiol. J. 2021, 36, 114–130. [Google Scholar] [CrossRef] [PubMed]
  2. Vaillancourt, D.E.; Newell, K.M. Aging and the time and frequency structure of force output variability. J. Appl. Physiol. 2003, 94, 903–912. [Google Scholar] [CrossRef] [Green Version]
  3. Lodha, N.; Misra, G.; Coombes, S.A.; Christou, E.A.; Cauraugh, J.H. Increased force variability in chronic stroke: Contributions of force modulation below 1 Hz. PLoS ONE 2013, 8, e83468. [Google Scholar] [CrossRef] [PubMed]
  4. Kent-Braun, J.A.; Callahan, D.M.; Fay, J.L.; Foulis, S.A.; Buonaccorsi, J.P. Muscle weakness, fatigue, and torque variability: Effects of age and mobility status. Muscle Nerve 2014, 49, 209–217. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Smith, J.W.; Christensen, J.C.; Marcus, R.L.; LaStayo, P.C. Muscle force and movement variability before and after total knee arthroplasty: A review. World J. Orthop. 2014, 5, 69. [Google Scholar] [CrossRef]
  6. Wang, W.-E.; Roy, A.; Misra, G.; Archer, D.B.; Ribeiro-Dasilva, M.C.; Fillingim, R.B.; Coombes, S.A. Motor-evoked pain increases force variability in chronic jaw pain. J. Pain Res. 2018, 19, 636–648. [Google Scholar] [CrossRef]
  7. Tracy, B.L.; Maluf, K.S.; Stephenson, J.L.; Hunter, S.K.; Enoka, R.M. Variability of motor unit discharge and force fluctuations across a range of muscle forces in older adults. Muscle Nerve 2005, 32, 533–540. [Google Scholar] [CrossRef]
  8. Taylor, A.M.; Christou, E.A.; Enoka, R.M. Multiple features of motor-unit activity influence force fluctuations during isometric contractions. J. Neurophysiol. 2003, 90, 1350–1361. [Google Scholar] [CrossRef] [Green Version]
  9. Galganski, M.E.; Fuglevand, A.J.; Enoka, R.M. Reduced control of motor output in a human hand muscle of elderly subjects during submaximal contractions. J. Neurophysiol. 1993, 69, 2108–2115. [Google Scholar] [CrossRef]
  10. Jones, K.E.; Hamilton, A.F.d.C.; Wolpert, D.M. Sources of signal-dependent noise during isometric force production. J. Neurophysiol. 2002, 88, 1533–1544. [Google Scholar] [CrossRef]
  11. Dideriksen, J.L.; Negro, F.; Enoka, R.M.; Farina, D. Motor unit recruitment strategies and muscle properties determine the influence of synaptic noise on force steadiness. J. Neurophysiol. 2012, 107, 3357–3369. [Google Scholar] [CrossRef] [Green Version]
  12. Yoshitake, Y.; Shinohara, M. Oscillations in motor unit discharge are reflected in the low-frequency component of rectified surface EMG and the rate of change in force. Exp. Brain Res. 2013, 231, 267–276. [Google Scholar] [CrossRef]
  13. Yoshitake, Y.; Shinohara, M. Low-frequency component of rectified EMG is temporally correlated with force and instantaneous rate of force fluctuations during steady contractions. Muscle Nerve 2013, 47, 577–584. [Google Scholar] [CrossRef] [PubMed]
  14. Park, S.H.; Casamento-Moran, A.; Yacoubi, B.; Christou, E.A. Voluntary reduction of force variability via modulation of low-frequency oscillations. Exp. Brain Res. 2017, 235, 2717–2727. [Google Scholar] [CrossRef] [PubMed]
  15. Hali, K.; Zero, A.M.; Rice, C.L. Effect of ankle joint position on triceps surae contractile properties and motor unit discharge rates. Physiol. Rep. 2021, 8, e14680. [Google Scholar] [CrossRef] [PubMed]
  16. Becker, R.; Awiszus, F. Physiological alterations of maximal voluntary quadriceps activation by changes of knee joint angle. Muscle Nerve 2001, 24, 667–672. [Google Scholar] [CrossRef]
  17. Shinohara, M.; Yoshitake, Y.; Kouzaki, M.; Fukunaga, T. The medial gastrocnemius muscle attenuates force fluctuations during plantar flexion. Exp. Brain Res. 2006, 169, 15–23. [Google Scholar] [CrossRef]
  18. Sosnoff, J.J.; Voudrie, S.J.; Ebersole, K.T. The effect of knee joint angle on torque control. J. Mot. Behav. 2009, 42, 5–10. [Google Scholar] [CrossRef]
  19. Pasquet, B.; Carpentier, A.; Duchateau, J. Change in muscle fascicle length influences the recruitment and discharge rate of motor units during isometric contractions. J. Neurophysiol. 2005, 94, 3126–3133. [Google Scholar] [CrossRef] [Green Version]
  20. Johansson, H.; Sjölander, P.; Sojka, P. A sensory role for the cruciate ligaments. Clin. Orthop. Relat. Res. 1991, 268, 161–178. [Google Scholar]
  21. George, D.; Mallery, P. IBM SPSS Statistics 26 Step by Step: A Simple Guide and Reference; Routledge: New York, NY, USA, 2019. [Google Scholar]
  22. Baweja, H.S.; Patel, B.K.; Martinkewiz, J.D.; Vu, J.; Christou, E.A. Removal of visual feedback alters muscle activity and reduces force variability during constant isometric contractions. Exp. Brain Res. 2009, 197, 35–47. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Ishikawa, M.; Pakaslahti, J.; Komi, P. Medial gastrocnemius muscle behavior during human running and walking. Gait Posture 2007, 25, 380–384. [Google Scholar] [CrossRef] [PubMed]
  24. Arnold, E.M.; Delp, S.L. Fibre operating lengths of human lower limb muscles during walking. Philos. Trans. R. Soc 2011, 366, 1530–1539. [Google Scholar] [CrossRef] [Green Version]
  25. Rubenson, J.; Pires, N.J.; Loi, H.O.; Pinniger, G.J.; Shannon, D.G. On the ascent: The soleus operating length is conserved to the ascending limb of the force–length curve across gait mechanics in humans. J. Exp. Biol. 2012, 215, 3539–3551. [Google Scholar] [CrossRef] [Green Version]
  26. Kanda, K.; Burke, R.; Walmsley, B. Differential control of fast and slow twitch motor units in the decerebrate cat. Exp. Brain Res. 1977, 29, 57–74. [Google Scholar] [CrossRef]
  27. Saltin, B. Skeletal muscle adaptability: Significant for metabolism and performance in Handbook of Physiology. Skeletal Muscle Bethesda. Am. Physiol. Soc. 1983, 10, 555–663. [Google Scholar]
  28. Chen, Y.-T.; Zhang, C.; Liu, Y.; Magat, E.; Verduzco-Gutierrez, M.; Francisco, G.E.; Zhou, P.; Zhang, Y.; Li, S. The effects of botulinum toxin injections on spasticity and motor performance in chronic stroke with spastic hemiplegia. Toxics 2020, 12, 492. [Google Scholar] [CrossRef] [PubMed]
  29. Riemann, B.L.; Lephart, S.M. The sensorimotor system, part II: The role of proprioception in motor control and functional joint stability. J. Athl. Train. 2002, 37, 80. [Google Scholar] [PubMed]
  30. Le, P.; Best, T.M.; Khan, S.N.; Mendel, E.; Marras, W.S. A review of methods to assess coactivation in the spine. J. Electromyogr. Kinesiol. 2017, 32, 51–60. [Google Scholar] [CrossRef]
  31. Christou, E.A. Aging and variability of voluntary contractions. Exerc. Sport Sci. Rev. 2011, 39, 77. [Google Scholar] [CrossRef] [Green Version]
  32. McCrum, C.; Leow, P.; Epro, G.; König, M.; Meijer, K.; Karamanidis, K. Alterations in leg extensor muscle-tendon unit biomechanical properties with ageing and mechanical loading. Front. Physiol. 2018, 9, 150. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Schematic diagram of experimental setup.
Figure 1. Schematic diagram of experimental setup.
Bioengineering 10 00373 g001
Figure 2. (a) Typical force trials of 30%MVIC for 4 s at longer muscle length from a representative subject. (b) From the same segment in (a), rsEMG from MG muscle; (c) rsEMG from LG muscle; (d) rsEMG from SOL muscle. All are detrended signals.
Figure 2. (a) Typical force trials of 30%MVIC for 4 s at longer muscle length from a representative subject. (b) From the same segment in (a), rsEMG from MG muscle; (c) rsEMG from LG muscle; (d) rsEMG from SOL muscle. All are detrended signals.
Bioengineering 10 00373 g002
Figure 3. (a) Absolute torque variability at different contraction intensities and ankle joint positions. (b) Relationships between absolute torque variability and actual plantar flexion torque at different ankle joint positions. Asterisk indicates a significant difference (* p < 0.050) between the longer muscle length and the shorter one. TSD indicates the absolute torque variability. Error bar indicates SEM.
Figure 3. (a) Absolute torque variability at different contraction intensities and ankle joint positions. (b) Relationships between absolute torque variability and actual plantar flexion torque at different ankle joint positions. Asterisk indicates a significant difference (* p < 0.050) between the longer muscle length and the shorter one. TSD indicates the absolute torque variability. Error bar indicates SEM.
Bioengineering 10 00373 g003
Figure 4. (a) Relative torque variability at different contraction levels and joint angles. (b) Relative torque variability at different plantar flexion torque and joint angles. Asterisk indicates a significant difference (* p < 0.050, ** p < 0.001) detected between different muscle lengths. TCoV indicates the relative torque variability. Error bars indicate SEM.
Figure 4. (a) Relative torque variability at different contraction levels and joint angles. (b) Relative torque variability at different plantar flexion torque and joint angles. Asterisk indicates a significant difference (* p < 0.050, ** p < 0.001) detected between different muscle lengths. TCoV indicates the relative torque variability. Error bars indicate SEM.
Bioengineering 10 00373 g004
Table 1. Absolute variability of MG, LG, and SOL (MGSD, LGSD, and SOLSD) at different contraction intensities and ankle joint angles. L: Longer muscle length (87°); N: Neutral muscle length (100°); S: Shorter muscle length (116°). Data shows as Mean ± SEM. Unit: µV.
Table 1. Absolute variability of MG, LG, and SOL (MGSD, LGSD, and SOLSD) at different contraction intensities and ankle joint angles. L: Longer muscle length (87°); N: Neutral muscle length (100°); S: Shorter muscle length (116°). Data shows as Mean ± SEM. Unit: µV.
MusclePosition10%
MVIC
20%
MVIC
30%
MVIC
50%
MVIC
70%
MVIC
MGL0.62 ± 0.130.85 ± 0.171.18 ± 0.211.99 ± 0.402.56 ± 0.51
N0.52 ± 0.090.87 ± 0.161.16 ± 0.242.04 ± 0.492.78 ± 0.64
S0.69 ± 0.120.98 ± 0.191.31 ± 0.272.10 ± 0.403.18 ± 0.70
LGL0.42 ± 0.050.65 ± 0.070.94 ± 0.101.67 ± 0.272.37 ± 0.40
N0.48 ± 0.090.59 ± 0.091.06 ± 0.241.52 ± 0.272.23 ± 0.39
S0.49 ± 0.080.69 ± 0.111.06 ± 0.201.75 ± 0.312.84 ± 0.46
SOLL0.30 ± 0.030.30 ± 0.050.50 ± 0.060.78 ± 0.111.01 ± 0.14
N0.35 ± 0.060.35 ± 0.080.57 ± 0.130.70 ± 0.140.84 ± 0.13
S0.31 ± 0.060.31 ± 0.050.48 ± 0.100.69 ± 0.141.02 ± 0.19
Table 2. Relative variability of MG, LG, and SOL (MGCoV, LGCoV,and SOLCoV) at different contraction intensities and ankle joint angles. L: Longer muscle length (87°); N: Neutral muscle length (100°); S: Shorter muscle length (116°). Data shows as Mean ± SEM. Unit: %.
Table 2. Relative variability of MG, LG, and SOL (MGCoV, LGCoV,and SOLCoV) at different contraction intensities and ankle joint angles. L: Longer muscle length (87°); N: Neutral muscle length (100°); S: Shorter muscle length (116°). Data shows as Mean ± SEM. Unit: %.
MusclePosition10%
MVIC
20%
MVIC
30%
MVIC
50%
MVIC
70%
MVIC
MGL9.98 ± 0.688.73 ± 0.508.41 ± 0.628.74 ± 0.628.25 ± 0.55
N8.81 ± 0.707.77 ± 0.527.52 ± 0.438.07 ± 0.618.13 ± 0.62
S8.66 ± 9.688.03 ± 0.507.60 ± 0.678.28 ± 0.568.42 ± 0.55
LGL8.30 ± 0.378.01 ± 0.728.04 ± 0.418.41 ± 0.768.77 ± 0.96
N8.17 ± 0.147.77 ± 0.618.46 ± 0.807.47 ± 0.628.05 ± 0.68
S8.19 ± 0.487.46 ± 0.517.64 ± 0.518.16 ± 0.388.41 ± 0.53
SOLL7.53 ± 0.406.76 ± 0.356.98 ± 0.337.31 ± 0.467.38 ± 0.31
N7.78 ± 0.367.46 ± 0.477.08 ± 0.367.36 ± 0.678.83 ± 1.17
S9.02 ± 0.938.30 ± 0.378.63 ± 0.588.80 ± 0.829.29 ± 0.68
Table 3. Stepwise regression model to identify possible predictors of the relative variability for torque (TCoV) out of the relative variabilities for EMG burst (MGCoV, LGCoV, and SOLCoV) across the contraction levels at each ankle joint angle. Data are shown as the estimated coefficient ±95% confidence interval. The significant predictors are in boldface type and the most important variable is highlighted. * p < 0.05. ** p < 0.001.
Table 3. Stepwise regression model to identify possible predictors of the relative variability for torque (TCoV) out of the relative variabilities for EMG burst (MGCoV, LGCoV, and SOLCoV) across the contraction levels at each ankle joint angle. Data are shown as the estimated coefficient ±95% confidence interval. The significant predictors are in boldface type and the most important variable is highlighted. * p < 0.05. ** p < 0.001.
MGLGSOL
Longer* 0.14 ± 0.080.04 ± 0.07* 0.18 ± 0.13
Neutral* 0.08 ± 0.07** 0.13 ± 0.070.04 ± 0.06
Shorter* 0.19 ± 0.11* 0.18 ± 0.14−0.01 ± 0.09
Table 4. Stepwise regression model to identify possible predictors of the relative variability for torque (TCoV) out of the relative variabilities for EMG burst (MGCoV, LGCoV, and SOLCoV) across the ankle joint angles at each contraction level. Data are shown as the estimated coefficient ±95% confidence interval. The significant predictors are in boldface type and the most important variable is highlighted. * p < 0.05. ** p < 0.001.
Table 4. Stepwise regression model to identify possible predictors of the relative variability for torque (TCoV) out of the relative variabilities for EMG burst (MGCoV, LGCoV, and SOLCoV) across the ankle joint angles at each contraction level. Data are shown as the estimated coefficient ±95% confidence interval. The significant predictors are in boldface type and the most important variable is highlighted. * p < 0.05. ** p < 0.001.
% MVICMGLGSOL
100.03 ± 0.150.17 ± 0.280.04 ± 0.16
20* 0.22 ± 0.140.04 ± 0.130.08 ± 0.16
30** 0.19 ± 0.08* 0.15 ± 0.09* 0.12 ± 0.10
50* 0.10 ± 0.09* 0.11 ± 0.09* 0.11 ± 0.08
70** 0.19 ± 0.100.04 ± 0.090.03 ± 0.07
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Shi, F.; Rymer, W.Z.; Son, J. Ankle Joint Angle Influences Relative Torque Fluctuation during Isometric Plantar Flexion. Bioengineering 2023, 10, 373. https://doi.org/10.3390/bioengineering10030373

AMA Style

Shi F, Rymer WZ, Son J. Ankle Joint Angle Influences Relative Torque Fluctuation during Isometric Plantar Flexion. Bioengineering. 2023; 10(3):373. https://doi.org/10.3390/bioengineering10030373

Chicago/Turabian Style

Shi, Fandi, William Zev Rymer, and Jongsang Son. 2023. "Ankle Joint Angle Influences Relative Torque Fluctuation during Isometric Plantar Flexion" Bioengineering 10, no. 3: 373. https://doi.org/10.3390/bioengineering10030373

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop