Next Article in Journal
Research on Mechanical Equipment Fault Diagnosis Method Based on Deep Learning and Information Fusion
Previous Article in Journal
Multi-Camera-Based Human Activity Recognition for Human–Robot Collaboration in Construction
Previous Article in Special Issue
A Multifunctional Network with Uncertainty Estimation and Attention-Based Knowledge Distillation to Address Practical Challenges in Respiration Rate Estimation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Investigation of Surface EMG Shorts-Derived Training Load during Treadmill Running

1
Faculty of Computing, Engineering and the Built Environment, Ulster University, Derry BT48 7JL, UK
2
School of Cardiovascular and Metabolic Health, College of Medical, Veterinary and Life Sciences, University of Glasgow, Glasgow G12 8QQ, UK
3
School of Health and Medical Sciences, University of Southern Queensland, Ipswich, QLD 4305, Australia
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(15), 6998; https://doi.org/10.3390/s23156998
Submission received: 8 June 2023 / Revised: 1 August 2023 / Accepted: 5 August 2023 / Published: 7 August 2023
(This article belongs to the Special Issue Wearable Sensors for Health and Physiological Monitoring)

Abstract

:
The purpose of this study was two-fold: (1) to determine the sensitivity of the sEMG shorts-derived training load (sEMG-TL) during different running speeds; and (2) to investigate the relationship between the oxygen consumption, heart rate (HR), rating of perceived exertion (RPE), accelerometry-based PlayerLoadTM (PL), and sEMG-TL during a running maximum oxygen uptake ( V ˙ O2max) test. The study investigated ten healthy participants. On day one, participants performed a three-speed treadmill test at 8, 10, and 12 km·h−1 for 2 min at each speed. On day two, participants performed a V ˙ O2max test. Analysis of variance found significant differences in sEMG-TL at all three speeds (p < 0.05). A significantly weak positive relationship between sEMG-TL and % V ˙ O2max (r = 0.31, p < 0.05) was established, while significantly strong relationships for 8 out of 10 participants at the individual level (r = 0.72–0.97, p < 0.05) were found. Meanwhile, the accelerometry PL was not significantly related to % V ˙ O2max (p > 0.05) and only demonstrated significant correlations in 3 out of 10 participants at the individual level. Therefore, the sEMG shorts-derived training load was sensitive in detecting a work rate difference of at least 2 km·h−1. sEMG-TL may be an acceptable metric for the measurement of internal loads and could potentially be used as a surrogate for oxygen consumption.

1. Introduction

Monitoring of the training load (TL) is an integral part of the training process in numerous sports and exercise regimes. The TL informs coaches of the physical and psychological stressors imposed on athletes during training and is categorised into two theoretical constructs: internal and external [1]. The internal TL captures the individual response to a given training stimulus and consists of objective physiological measures such as the heart rate (HR) or subjective psychophysiological methods such as the session rating of perceived exertion (sRPE) [2]. In contrast, the external load is the physical work completed by the individual. It includes measures such as global positioning systems (GPS), inertial sensors, time motion analysis (TMA), and local positioning systems (LPS) [3,4,5]. Combining the external TL with the internal response, or physiological changes in response to the load, is known as the dose–response relationship and is considered an important component in assessing the effectiveness of the training process [6]. It is suggested that adequate TL prescription could mitigate the chances of soft tissue injury, reduce illness, and increase game readiness [7,8].
It is recommended that a variety of internal and external load metrics are captured simultaneously as they are usually not interchangeable [9]. For example, HR-based measures of internal load demonstrate moderate–strong positive relationships with GPS metrics, such as jogging and total distance (r = 0.66), but no other GPS metrics, such as high-intensity running, sprinting, and maximum speed [4]. Contextual factors such as the upcoming opponent and previous game result (win, loss, and draw) can impact subjective TL methods, such as the RPE and sRPE, due to the psychological element that they inherit [4]. Thus, other objective internal TL monitoring approaches are recommended to capture different parameters of the individual response to sports and exercise.
Surface electromyography (sEMG) is a useful method for the evaluation of muscle activity, particularly in dynamic exercise. It has gained extensive usage in studying the electrophysiological aspects of muscle contractions and force generation within the realms of sports and exercise medicine [10]. Normally, the connection of wires to diagnostic equipment restricts free movement, which is a limiting factor in exercise-based sEMG protocols. The integration of sEMG sensors into clothing fabric, also referred to as “textile sEMG electrodes” [11,12,13,14], is potentially a novel technique for the objective measurement of the sEMG-derived training load (sEMG-TL). This would offer a convenient solution to the use of sEMG during dynamic exercise or sports activities in the field [14]. Strive™, Athos™, and Myontec™ provide wearable sEMG solutions claiming to capture sEMG-TL during exercise. Athos™ defines sEMG-TL as the total muscle activation from all sensors, divided by a scaling factor. While promising results exist for textile sEMG sensors in capturing muscular responses during isokinetic and functional exercise [11,14,15,16], research on sEMG-TL during dynamic exercise is limited.
Previous research has shown the strong validity and reliability of textile sEMG sensors embedded into compression shorts compared to a gold-standard BIOPAC system during knee flexion and extension, performed at different intensities [15]. Current research has used textile sEMG sensors to study muscle activity during controlled, closed-chain exercise [15,17], cycling [16], low-intensity daily locomotive activity [18], and functional exercise [11,14]. With claims that wearable sEMG shorts can be used in sports and exercise to provide a measure of TL, more research is warranted for assessing the sEMG derived TL metric given by the shorts during exercise.
The association between the integrated EMG (iEMG) of the quadriceps, hamstrings, and calf muscles and oxygen consumption suggests that muscle activation may be a representation of individuals’ physical exertion [19,20,21]. Thus, sEMG-TL may serve as a possible surrogate for oxygen consumption during specific treadmill running protocols, such as a V ˙ O2max test. To date, one study has demonstrated how sEMG shorts may be used as a viable predictor of lactate threshold work rates during an incremental cycling protocol [16]. This demonstrates the versatility of the use of sEMG shorts during exercise and they may be a viable tool in capturing the TL during dynamic whole-body exercise, such as running.
CatapultTM utilises a triaxial accelerometer to calculate the PlayerLoadTM (PL), an external load variable expressed as the square root of the sum of squared instantaneous acceleration changes in X, Y, and Z vectors, which has been validated in various sports and exercise research protocols [22,23].
External loads, such as PL, are often combined with internal loads such as HR measures to understand how prescribed training loads relate to individual athlete responses [24,25]. However, as previous research shows that accelerometry’s energy expenditure predictions decrease as the slope of uphill running increases due to inertial sensors’ limitations in detecting ground slope changes, a weak relationship between PL and oxygen consumption would be expected [26]. In contrast, sEMG-TL may better account for the changes in running slope during a treadmill running test as it accounts for the internal response.
Therefore, this study aims to (1) investigate the sensitivity of sEMG-TL across different running speeds; (2) explore and compare the relationships of sEMG-TL, HR, RPE, and PL with oxygen consumption during a maximal running test; and (3) examine the associations of sEMG-TL with HR, RPE, and PL. Understanding these relationships will contribute to advancing future research in the field of sports and exercise TL monitoring. Moreover, the findings will have practical implications, guiding the potential application of sEMG shorts as a tool for the monitoring of training loads in sports and exercise settings.

2. Materials and Methods

2.1. Participants

Six healthy, recreationally active males and four healthy, recreationally active females volunteered to participate in this study (Table 1). Written informed consent was obtained from all subjects involved in the study following a full explanation of the study’s research objectives, procedures, risks, and benefits. Ethical approval was obtained from the University of Glasgow Institutional Ethics Committee, and data collection procedures were in accordance with the Declaration of Helsinki.

2.2. Experimental Design

Participants visited the laboratory on two consecutive days. On the first day, they undertook preliminary anthropometric tests, including body fat percentage recording using an air displacement plethysmograph (BodPod®, COSMED, Rome, Italy). Recording of the body fat percentage was deemed necessary as excess subcutaneous fat has been shown to weaken sEMG signals and dampen the peak amplitude [27]. The participants’ thigh length and hip circumference were measured to appropriately fit them with sEMG shorts (AthosTM, Redwood City, CA, USA). A neoprene vest that involved holding a triaxial accelerometer unit (OptimEye™ S5, Catapult Innovations, Melbourne, Australia) between the scapulae was also fitted to each participant. Participants were familiarised with the wearable technology and a stationary sEMG calibration protocol, which was undertaken to configure the sEMG thresholds for each muscle group. The sEMG sensors were positioned to measure four muscle groups on each leg (inner and outer quadriceps, hamstrings, and glutes). Following the calibration protocol, participants completed a 3-speed treadmill running test. This test was conducted to assess whether the sEMG-TL could detect small changes in workload, as, to date, no research has explored the sEMG-TL.
On the second day, participants completed a treadmill running maximum oxygen uptake ( V ˙ O2max) test, during which various internal loads (sEMG-TL, RPE, HR, V ˙ O2) and external loads (accelerometery) were continuously monitored.

2.3. Methodology

On the first day, participants wore the sEMG shorts and were allocated an Athos unit (Figure 1) and completed the sEMG calibration protocol, which included four specific movements: seated leg extension, supine leg raise, prone knee flexion, and prone hip extension. Each movement was repeated four times while the researcher manually applied varying passive forces (low, medium, and high) and provided resistance to induce an isometric maximal voluntary contraction (MVC) (Figure 2). Following the sEMG calibration, each participant was fitted with a Polar H10 heart rate chest monitor (Polar Electro Oy, Kempele, Finland) and a triaxial accelerometer. Participants were then required to undertake a 5-min jogging warm-up on a treadmill (Woodway, ELG 70 Weiss, Weil am Rhein, Germany). The warm-up was important to increase the stability of the silicone overlaid electrode-to-skin interface, and to promote perspiration to increase skin conductivity. The 3-speed treadmill test then commenced, which required the participants to run at 8 (low), 10 (mod), and 12 (high) km⋅h−1 for 2 min at each speed.
On the second day, participants completed an incremental treadmill exercise test to determine V ˙ O2max. Males ran at 10 km⋅h−1, while females ran at 8 km⋅h−1, both on a 1% gradient increase every 1 min until reaching volitional exhaustion. Breath-by-breath pulmonary gas exchange was measured using a metabolic cart (MedGraphicsTM, Gloucester, UK) to measure V ˙ O2 (mL·kg−1·min−1), which was averaged over consecutive 60 s intervals. At the end of each 60 s interval, the participants’ heart rate (bpm) and RPE (Borg’s Category Ratio 6–20 RPE scale) were also measured.

2.4. Data Collection and Processing

The Athos™ unit electronics subsystem is responsible for signal reception, transmission, and conditioning. Within the electronics subsystem is a signal conditioning module that filters the sEMG data. Sampled sEMG signals were captured at 1 kHz. The anti-aliasing filter applied prior to sampling prevented high-frequency noise greater than 500 Hz from aliasing into the sEMG spectrum. Filtering included a linear band-pass with a cut-off frequency from 10 Hz to 500 Hz (23 dB frequencies with centre frequency at 120 Hz); a notch (removal of 60 Hz noise) filter, rectification, and linear envelope were applied. The linear envelope was down-sampled by a factor of 25 and smoothed using a 16-sample root-mean-square (RMS) transformation along with signal conversion from analogue to digital. sEMG signals were averaged at a 10 s window. sEMG recordings were considered viable when impedance was below 5 kΩ. Poor contact quality was deemed if loss of the contact signal occurred over 10% of the time. These filtering and signal processes were predetermined by the Athos™.
The sEMG peak amplitude, which was captured through the calibration protocol, provided individualised reference points for each muscle group, which the raw integrated sEMG output (area under the curve of the rectified EMG signal) could be normalised against. The integrated sEMG for each muscle group was measured as a percentage of the MVC. The accumulation of the normalised integrated sEMG values across all muscle groups relative to the MVC peak sEMG amplitude represented the sEMG-TL reported in arbitrary units (a.u.). A single arbitrary unit was equivalent to one muscle activating at 100% of the MVC for 1-s.
The triaxial accelerometer sampled at 100 Hz. The data were accumulated to generate a measure of the external load (PL). The equation for PL can be seen below:
PL = a x i a x i 1 2 + a y i a y i 1 2 + a z i a z i 1 2
where a x is the mediolateral acceleration, a y the anteroposterior acceleration, and a z the vertical acceleration.

2.5. Statistical Analysis

The data were analysed using the statistical software package SPSS (version 26.0). Since we had a small sample size, determination of the distribution of the data for each variable (sEMG-TL, HR, RPE, PL, and % V ˙ O2max) using a Shapiro–Wilk test was performed. Results showed that the distributions for all variables (sEMG-TL, HR, RPE, PL, and %VO2max) departed significantly from normality (W = 0.92–0.96, p < 0.05), as well as the 3-speed test data (W = 0.82, p < 0.05). Descriptive statistics were determined for the participants’ sEMG-TL data from the 3-speed test and for all variables during each minute in the V ˙ O2max test. The sensitivity of the sEMG-TL in detecting 2 km⋅h−1 changes in running intensity during the 3-speed treadmill test was assessed by a Kruskal–Wallis one-way analysis of variance by ranks test, and pairwise comparisons were performed using a Dunn–Bonferroni post hoc test. Statistical significance was set at p < 0.05. Raw accelerometry data and V ˙ O2 (mL·kg−1·min−1) were averaged at every 60 s interval during the V ˙ O2max test. V ˙ O2 was normalised to V ˙ O2max (% V ˙ O2max). Thus, a non-parametric Spearman’s correlation test was applied to determine the correlations between all variables in the composite data set, while separate correlations were applied on the individual data sets for sEMG-TL, HR, RPE, PL, and % V ˙ O2max. The significance level was set at p < 0.05. Correlations were classified as previously suggested: 0 = no correlation, (0 < rs < 0.2) = very weak correlation, (0.2 ≤ rs < 0.4) = weak correlation, (0.4 ≤ rs < 0.6) = moderate correlation, (0.6 ≤ rs < 0.8) = strong correlation, (0.8 ≤ rs < 1.0) = very strong correlation, and 1 = perfect correlation [28].

3. Results

3.1. Three-Speed Treadmill Test

sEMG-TL is presented as a median (interquartile range (IQR)) for each running speed. Figure 3 illustrates the spread of data between each running speed. A low running speed resulted in the lowest median sEMG-TL, 254.17 a.u. (IQR = 148.50 to 204.27 a.u.). A high running speed resulted in the highest sEMG-TL, 236.67 a.u. (IQR = 192.10 to 295.48 a.u.), and the sEMG-TL for a moderate running speed was between those of the low and high speeds, 208.07 (IQR = 183.51 to 244.29 a.u.). The Kruskal–Wallis one-way analysis of variance by ranks test revealed a significant difference in sEMG-TL among the three running speeds, χ2(2) = 11.50, p < 0.05. Post hoc pairwise comparisons using the Dunn–Bonferroni test demonstrated that the sEMG-TL at a low running speed was significantly different from both the moderate and high running speeds (adjusted p < 0.05), but no statistical difference was established between moderate and high running speeds (adjusted p > 0.05).

3.2. Treadmill V ˙ O2max Test

There were 74 individual recordings included in the analysis for each dependent variable. Median (IQR) loads are presented in Table 2 for each measurement every 1 min during the V ˙ O2max test. Generally, the median load for each variable increased at every 1 min interval during the V ˙ O2max test.
Figure 4 presents illustrates the Spearman’s rank correlation coefficients between each load variable for the composite data set. While HR and RPE showed significant strong positive correlations (rs = 0.70, p = 0.002) and rs = 0.91, p = 0.001, respectively) with % V ˙ O2max, sEMG-TL revealed a very weak correlation (rs = 0.31, p = 0.02). In contrast, no correlation was established between PL and % V ˙ O2max (rs = 0.04, p = 0.74). sEMG-TL showed weak significant correlations between RPE (rs = 0.39, p = 0.01), HR (rs = 0.24, p = 0.03) and PL (rs = 0.42, p = 0.02).
Correlation coefficients at the individual level are presented in Table 3. Most notably, sEMG-TL demonstrated very strong significant positive correlations with % V ˙ O2max in nine out of ten participants, while only three participants showed very strong non-significant positive correlations between PL and % V ˙ O2max. In addition, Table 3 presents very strong significant positive correlations between RPE and sEMG-TL in seven out of ten participants, while three participants showed non-significant strong correlations. Similarly, HR and sEMG-TL were very strongly significantly correlated in eight out of ten participants, while two participants showed strong non-significant positive correlations.

4. Discussion

The purpose of this study was to (1) investigate the sensitivity of sEMG-TL across different running speeds; (2) explore and compare the relationships of sEMG-TL, HR, RPE, and PL with oxygen consumption during a maximal running test; and (3) examine the associations of sEMG-TL with HR, RPE, and PL. An important feature of this study is the novelty of the textile sEMG sensors embedded in compression shorts to provide an sEMG-based internal TL metric as a method for the objective monitoring of the muscular load during dynamic exercise. The results show that sEMG-TL exhibited significant differences among three different running intensities, distinguished by increments of 2 km⋅h−1 in running velocity. As expected, higher running velocities were associated with an increased sEMG-TL. Figure 3 illustrates the distribution of participants’ sEMG-TL for each speed. A significant difference was shown between low–moderate and low–high, but, surprisingly, not moderate–high speeds. While the sEMG-TL was higher, the larger spread of the data was likely attributed to a non-significant difference. In fact, the spread of the data was greater as the speed increased. This may imply that the sEMG sensors become less accurate at higher intensities. However, it is important to also consider the individual biomechanical and physiological factors that could result in interindividual differences that affect the sEMG signal. Overall, these findings reflect the increased excitability of the muscles in response to the escalating demands imposed by running at higher velocities. Previous findings support these findings, showing increased sEMG signal amplitudes in running compared to low-speed jogging [29]. These findings indicate that sEMG-TL is sensitive in detecting small changes in intensity during exercise.
Considering the second aim of the study, a significant yet somewhat modest association between sEMG-TL and % V ˙ O2max was established (Figure 4). However, at the individual level, very strong significant correlations between sEMG-TL and % V ˙ O2max were established in ninety percent of the participants (Table 3). The interplay of participants’ unique biomechanical and physiological characteristics likely resulted in a weaker relationship between sEMG-TL and % V ˙ O2max in the pooled data set. Previous research identified that when analysing grouped data, the median frequency values of integrated EMG were insensitive to changes in muscle recruitment, metabolite accumulation, and fatigue associated with the increases in work intensity during a cycle ergometer V ˙ O2max test [30]. However, the integrated EMG results were highly individualistic between some subjects and muscles, which could help to explain the findings in the present study. Previous studies have also reported similar results regarding the intra-individual variability based on the sEMG signal amplitude [31,32]. During the V ˙ O2max test, the velocity of the treadmill was maintained at a constant level while the gradient of the treadmill increased by 1% every minute. As more force output is required every minute during the test, there is a higher demand for oxygen by the oxidative pathway, which metabolises to produce adenosine triphosphate (ATP) for energy. This was reflected through the V ˙ O2 measured from the participant. The higher force output requirements and associated fatigue as the test progresses will stimulate the muscles to increase the electrical discharge to compensate for the decreasing muscle strength, resulting in an increased sEMG amplitude [33]. During fatiguing exercise, the more frequent recruitment of larger motor unit action potentials is necessary to sustain the required force output during the test, which increases the sEMG amplitudes and, in part, explains the association between sEMG-TL and % V ˙ O2max [33].
HR and RPE demonstrated very strong significant correlations with % V ˙ O2max both in the composite data set and at the individual level for all participants. These findings align with the previous literature, suggesting that both HR and RPE are reliable indicators of intensity and load in the context of load monitoring [2,34,35]. It is important to emphasise that the load monitoring approaches employed in this study captured distinct parameters of the load from the participants, as exemplified through the varying degrees of correlation with % V ˙ O2max. HR, being a measure of the heart’s response to exercise, appears to be particularly effective in representing the oxygen output during physical activity, as anticipated, given its direct reflection of the cardiovascular system’s response to exertion. These results support the use of HR and RPE for load monitoring, as they highlight the utility of HR and RPE as practical measures to assess % V ˙ O2max during exercise. By capturing complementary aspects of the participants’ physiological responses, HR and RPE provide valuable insights into the intensity of the exercise and are more strongly related to % V ˙ O2max than sEMG-TL and PL.
PL demonstrated no correlation with % V ˙ O2max for the grouped data. This highlights the different parameters that internal and external player monitoring systems capture during exercise [2,25,36]. The type of treadmill running test used in the present study may have negatively influenced the relationship between PL and % V ˙ O2max. It has been previously shown that sEMG from the quadriceps and hamstrings may better reflect energy expenditure during uphill and downhill locomotion compared to accelerometry [37]. The metabolic rate during running is determined by two factors: (a) the rate of muscle force development, such as during the stance phase in running, where the tensile properties of soft tissue produce force; and (b) the volume of active leg muscle [38,39]. Thus, as the gradient of the treadmill increases, the sEMG-TL may be a better representation of the load experienced by the individual as PL only captures acceleration.
In regard to the third aim of the study, the current findings established a significant weak association between sEMG-TL and RPE in the composite data set (Figure 4). Previous research has illustrated similar results [40]. Fontes et al. found a moderate positive correlation between RPE and sEMG during exhaustive constant-load cycling bouts. Previous research has also indicated that the recruitment of additional muscle fibres is associated with increased RPE [41]. In addition, type II muscle fibre recruitment has been shown to increase the sEMG amplitude, which occurs at greater intensities nearing the end of the V ˙ O2max test, which in part may explain this relationship. At the individual level, sEMG-TL demonstrated very strong significant correlations with RPE in seventy percent of participants. Establishing a relationship between RPE and sEMG-TL is an important finding as this could help with exercise prescription and in designing routines for specific muscle loads.
In addition, a significant, very weak correlation was found between HR and sEMG-TL in the composite data, yet very strong significant correlations were established in eighty percent of participants when the HR and sEMG-TL correlations were analysed at the individual level. While there is no direct causal relationship between sEMG and HR, they are often positively associated during exercise [42,43]. As muscle activity increases, indicated by sEMG-TL, HR increases as well to meet the metabolic demands.

Strengths and Limitations

Most of the previous research conducted on textile sEMG sensors embedded in shorts for the monitoring of muscle activity involves isokinetic exercises, general locomotion, or low-level functional exercise [14,15,16,17,18]. In general, most studies investigate specific muscle groups’ sEMG responses, whereas this study used sEMG shorts to derive a TL metric, which opens new research possibilities in monitoring sEMG-TL in different sports and exercise environments. Additionally, sEMG-TL appears to closely correlate with % V ˙ O2max at the individual level; thus, the shorts could serve as a potential surrogate for the measurement of V ˙ O2 during incremental running.
Individual biomechanical and physiological characteristics likely impact the relationship between sEMG-TL and the other load variables [30,31,32]. A limiting factor in the current research was the lack of control of contextual factors, such as the participants’ training experience, fitness level, running technique, and muscle fibre type composition. For example, previous research has shown that athletes who are fitter and inherit better technique recruit muscle at different rates compared to non-athletes, especially in higher threshold motor units [44]. Thus, ensuring sample homogeneity, such as participants with similar fitness characteristics, could reduce confounding variables and enhance the internal validity of the study. In addition, a weakness of this study was that it entailed a small sample size. A larger sample size could potentially enhance the statistical power and generalisability of the results.
Another limitation of the current study is that, while HR and RPE are global measurements of internal load, sEMG represents a portion of superficial muscle (upper legs, only), yet the participants likely experienced physiological stress elsewhere in the body. Therefore, caution should be taken when using the data to inform exercise loading parameters where exercise includes more muscles other than the upper leg. In the future, research might consider using a full-body sEMG compression garment that encompasses many more sEMG sensors to better reflect the whole-body muscle load.

5. Conclusions

This study aimed to investigate the sensitivity of sEMG-TL across different running speeds, explore its relationship with oxygen consumption, and examine its associations with HR, RPE, and PL. The novelty of using textile sEMG sensors embedded in compression shorts provided an objective method for the monitoring of the muscular load during dynamic exercise.
The results revealed significant differences in sEMG-TL among different running intensities, with higher velocities associated with increased sEMG-TL. The study also established a significant yet weak association between sEMG-TL and % V ˙ O2max at the group level, but strong individual correlations were found. HR and RPE showed very strong significant correlations with % V ˙ O2max, supporting their reliability for the assessment of exercise loads. However, PL demonstrated no correlation with   % V ˙ O2max, indicating the limitations of accelerometry in capturing loads during treadmill running with gradient increases in the uphill slope. Additionally, a significant weak association was found between sEMG-TL and RPE at the composite level, while showing more substantial individual-level correlations. The study’s strengths include the use of sEMG shorts to derive a TL metric and its potential application for individual exercise prescription. However, factors such as the small sample size and the need to control contextual factors were identified as study limitations.
In conclusion, this research emphasises the sensitivity of sEMG-TL in detecting small changes in exercise intensity and highlights the potential for future research to be conducted with training load monitoring using sEMG. Future research should consider larger sample sizes and use full-body sEMG compression garments to better reflect the whole-body muscle load during exercise.

Author Contributions

Conceptualisation, K.A., V.P. and A.W.; methodology, V.P., K.A., T.R. and J.C.; formal analysis, K.A.; writing—original draft preparation, K.A. and V.P.; writing—review and editing, V.P., T.R. and S.P.B.; supervision, V.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Institutional Review Board (or Ethics Committee) of IRCCS Fondazione Don Gnocchi ONLUS (number 13663_oss).

Informed Consent Statement

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

Data Availability Statement

The data generated in this study are not available in any public network database. For anyone who would like to have access to the data, please contact the author of the study, K.A., email: [email protected].

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. McLaren, S.J.; Macpherson, T.W.; Coutts, A.J.; Hurst, C.; Spears, I.R.; Weston, M. The Relationships Between Internal and External Measures of Training Load and Intensity in Team Sports: A Meta-Analysis. Sports Med. 2018, 48, 641–658. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Impellizzeri, F.M.; Marcora, S.M.; Coutts, A.J. Internal and External Training Load: 15 Years On. Int. J. Sports Physiol. Perform. 2019, 14, 270–273. [Google Scholar] [CrossRef]
  3. Manchado, C.; Pueo, B.; Chirosa-Rios, L.J.; Tortosa-Martínez, J. Time–Motion Analysis by Playing Positions of Male Handball Players during the European Championship 2020. Int. J. Environ. Res. Public Health 2021, 18, 2787. [Google Scholar] [CrossRef]
  4. Sobolewski, E.J. The Relationships between Internal and External Load Measures for Division I College Football Practice. Sports 2020, 8, 165. [Google Scholar] [CrossRef] [PubMed]
  5. Alt, P.S.; Baumgart, C.; Ueberschär, O.; Freiwald, J.; Hoppe, M.W. Validity of a Local Positioning System during Outdoor and Indoor Conditions for Team Sports. Sensors 2020, 20, 5733. [Google Scholar] [CrossRef] [PubMed]
  6. Impellizzeri, F.M.; Shrier, I.; McLaren, S.J.; Coutts, A.J.; McCall, A.; Slattery, K.; Jeffries, A.C.; Kalkhoven, J.T. Understanding Training Load as Exposure and Dose. Sports Med. 2023, 53, 1–13. [Google Scholar] [CrossRef] [PubMed]
  7. Gabbett, T.J. The Training—Injury Prevention Paradox: Should Athletes Be Training Smarter and Harder? Br. J. Sports Med. 2016, 50, 273–280. [Google Scholar] [CrossRef] [Green Version]
  8. Drew, M.K.; Finch, C.F. The Relationship Between Training Load and Injury, Illness and Soreness: A Systematic and Literature Review. Sports Med. 2016, 46, 861–883. [Google Scholar] [CrossRef]
  9. Gonçalves, L.G.C.; Kalva-Filho, C.A.; Nakamura, F.Y.; Rago, V.; Afonso, J.; de Bedo, B.L.S.; Aquino, R. Effects of Match-Related Contextual Factors on Weekly Load Responses in Professional Brazilian Soccer Players. Int. J. Env. Res. Public. Health 2020, 17, 5163. [Google Scholar] [CrossRef]
  10. De Luca, C.J. The Use of Surface Electromyography in Biomechanics. J. Appl. Biomech. 1997, 13, 135–163. [Google Scholar] [CrossRef] [Green Version]
  11. Finni, T.; Hu, M.; Kettunen, P.; Vilavuo, T.; Cheng, S. Measurement of EMG Activity with Textile Electrodes Embedded into Clothing. Physiol. Meas. 2007, 28, 1405–1419. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Guo, L.; Sandsjö, L.; Ortiz-Catalan, M.; Skrifvars, M. Systematic Review of Textile-Based Electrodes for Long-Term and Continuous Surface Electromyography Recording. Text. Res. J. 2020, 90, 227–244. [Google Scholar] [CrossRef] [Green Version]
  13. Kim, S.; Lee, S.; Jeong, W. EMG Measurement with Textile-Based Electrodes in Different Electrode Sizes and Clothing Pressures for Smart Clothing Design Optimization. Polymers 2020, 12, 2406. [Google Scholar] [CrossRef] [PubMed]
  14. Colyer, S.L.; McGuigan, P.M. Textile Electrodes Embedded in Clothing: A Practical Alternative to Traditional Surface Electromyography When Assessing Muscle Excitation during Functional Movements. J. Sports Sci. Med. 2018, 17, 101. [Google Scholar]
  15. Lynn, S.K.; Watkins, C.M.; Wong, M.; Balfany, K.; Feeney, D.F. Validity and Reliability of Surface Electromyography Measurements from a Wearable Athlete Performance System. J. Sports Sci. Med. 2018, 17, 205. [Google Scholar]
  16. Snarr, R.L.; Tolusso, D.V.; Hallmark, A.V.; Esco, M.R. Validity of Wearable Electromyographical Compression Shorts to Predict Lactate Threshold during Incremental Exercise in Healthy Subjects. J. Strength. Cond. Res. 2021, 35, 702–708. [Google Scholar] [CrossRef] [PubMed]
  17. Svensson, F.; Aasa, U.; Strong, A. Textile Electromyography Electrodes Reveal Differences in Lower Limb Muscle Activation during Loaded Squats When Comparing Fixed and Free Barbell Movement Paths. Front. Sports Act. Living 2022, 4, 1021323. [Google Scholar] [CrossRef]
  18. Tikkanen, O.; Haakana, P.; Pesola, A.J.; Häkkinen, K.; Rantalainen, T.; Havu, M.; Pullinen, T.; Finni, T. Muscle Activity and Inactivity Periods during Normal Daily Life. PLoS ONE 2013, 8, e52228. [Google Scholar] [CrossRef] [Green Version]
  19. Kyröläinen, H.; Belli, A.; Komi, P.V. Biomechanical Factors Affecting Running Economy. Med. Sci. Sports Exerc. 2001, 33, 1330–1337. [Google Scholar] [CrossRef]
  20. Bigland-Ritchie, B.; Woods, J.J. Integrated Electromyogram and Oxygen Uptake during Positive and Negative Work. J. Physiol. 1976, 260, 267. [Google Scholar] [CrossRef]
  21. Watanabe, K.; Narouei, S. Association between Oxygen Consumption and Surface Electromyographic Amplitude and Its Variation within Individual Calf Muscles during Walking at Various Speeds. Sensors 2021, 21, 1748. [Google Scholar] [CrossRef] [PubMed]
  22. Howe, S.T.; Aughey, R.J.; Hopkins, W.G.; Cavanagh, B.P.; Stewart, A.M. Sensitivity, Reliability and Construct Validity of GPS and Accelerometers for Quantifying Peak Periods of Rugby Competition. PLoS ONE 2020, 15, e0236024. [Google Scholar] [CrossRef] [PubMed]
  23. Nicolella, D.P.; Torres-Ronda, L.; Saylor, K.J.; Schelling, X. Validity and Reliability of an Accelerometer-Based Player Tracking Device. PLoS ONE 2018, 13, e0191823. [Google Scholar] [CrossRef] [PubMed]
  24. Marynowicz, J.; Kikut, K.; Lango, M.; Horna, D.; Andrzejewski, M. Relationship Between the Session-RPE and External Measures of Training Load in Youth Soccer Training. J. Strength. Cond. Res. 2020, 34, 2800–2804. [Google Scholar] [CrossRef] [PubMed]
  25. Fox, J.L.; O’Grady, C.J.; Scanlan, A.T. The Relationships Between External and Internal Workloads During Basketball Training and Games. Int. J. Sports Physiol. Perform. 2020, 15, 1081–1086. [Google Scholar] [CrossRef]
  26. Chang, C.H.; Lin, K.C.; Ho, C.S.; Huang, C.C. Accuracy of the Energy Expenditure during Uphill Exercise Measured by the Waist-Worn ActiGraph. J. Exerc. Sci. Fit. 2019, 17, 62–66. [Google Scholar] [CrossRef]
  27. Kuiken, T.A.; Lowery, M.M.; Stoykov, N.S. The Effect of Subcutaneous Fat on Myoelectric Signal Amplitude and Cross-Talk. Prosthet. Orthot. Int. 2003, 27, 48–54. [Google Scholar] [CrossRef] [Green Version]
  28. Akoglu, H. User’s Guide to Correlation Coefficients. Turk. J. Emerg. Med. 2018, 18, 91–93. [Google Scholar] [CrossRef]
  29. Gazendam, M.G.J.; Hof, A.L. Averaged EMG Profiles in Jogging and Running at Different Speeds. Gait Posture 2007, 25, 604–614. [Google Scholar] [CrossRef]
  30. Miura, H.; Araki, H.; Matoba, H.; Kitagawa, K. Relationship among Oxygenation, Myoelectric Activity, and Lactic Acid Accumulation in Vastus Lateralis Muscle during Exercise with Constant Work Rate. Int. J. Sports Med. 2000, 21, 180–184. [Google Scholar] [CrossRef]
  31. Martens, J.; Daly, D.; Deschamps, K.; Fernandes, R.J.P.; Staes, F. Intra-Individual Variability of Surface Electromyography in Front Crawl Swimming. PLoS ONE 2015, 10, e0144998. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Guidetti, L.; Rivellini, G.; Figura, F. EMG Patterns during Running: Intra- and Inter-Individual Variability. J. Electromyogr. Kinesiol. 1996, 6, 37–48. [Google Scholar] [CrossRef] [PubMed]
  33. Linssen, W.H.J.P.; Stegeman, D.F.; Joosten, E.M.G.; Binkhorst, R.A.; Merks, M.J.H.; Laak, H.J.T.; Notermans, S.L.H. Fatigue in Type I Fiber Predominance: A Muscle Force and Surface EMG Study on the Relative Role of Type I and Type II Muscle Fibers. Muscle Nerve 1991, 14, 829–837. [Google Scholar] [CrossRef]
  34. Bourdon, P.C.; Cardinale, M.; Murray, A.; Gastin, P.; Kellmann, M.; Varley, M.C.; Gabbett, T.J.; Coutts, A.J.; Burgess, D.J.; Gregson, W.; et al. Monitoring Athlete Training Loads: Consensus Statement. Int. J. Sports Physiol. Perform. 2017, 12, 161–170. [Google Scholar] [CrossRef]
  35. Halson, S.L. Monitoring Training Load to Understand Fatigue in Athletes. Sports Med. 2014, 44 (Suppl. S2), 139. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Curtis, R.M.; Huggins, R.A.; Benjamin, C.L.; Sekiguchi, Y.; Adams, W.M.; Arent, S.M.; Jain, R.; Miller, S.J.; Walker, A.J.; Casa, D.J. Contextual Factors Influencing External and Internal Training Loads in Collegiate Men’s Soccer. J. Strength. Cond. Res. 2020, 34, 374–381. [Google Scholar] [CrossRef]
  37. Tikkanen, O.; Kärkkäinen, S.; Haakana, P.; Kallinen, M.; Pullinen, T.; Finni, T. EMG, Heart Rate, and Accelerometer as Estimators of Energy Expenditure in Locomotion. Med. Sci. Sports Exerc. 2014, 46, 1831–1839. [Google Scholar] [CrossRef]
  38. Kipp, S.; Grabowski, A.M.; Kram, R. What Determines the Metabolic Cost of Human Running across a Wide Range of Velocities? J. Exp. Biol. 2018, 221, jeb184218. [Google Scholar] [CrossRef] [Green Version]
  39. Bastien, G.J.; Willems, P.A.; Schepens, B.; Heglund, N.C. Effect of Load and Speed on the Energetic Cost of Human Walking. Eur. J. Appl. Physiol. 2005, 94, 76–83. [Google Scholar] [CrossRef]
  40. Fontes, E.B.; Smirmaul, B.P.C.; Nakamura, F.Y.; Pereira, G.; Okano, A.H.; Altimari, L.R.; Dantas, J.L.; De Moraes, A.C. The Relationship between Rating of Perceived Exertion and Muscle Activity during Exhaustive Constant-Load Cycling. Int. J. Sports Med. 2010, 31, 683–688. [Google Scholar] [CrossRef] [Green Version]
  41. Ni, W.; Zhang, Y.; Li, X.; Wang, X.; Wu, Y.; Liu, G. A Study on the Relationship between RPE and SEMG in Dynamic Contraction Based on the GPR Method. Electronics 2022, 11, 691. [Google Scholar] [CrossRef]
  42. Seals, D.R.; Enoka, R.M. Sympathetic Activation Is Associated with Increases in EMG during Fatiguing Exercise. J. Appl. Physiol. 1989, 66, 88–95. [Google Scholar] [CrossRef] [PubMed]
  43. Fallentin, N.; Sidenius, B.; Jorgensen, K. Blood Pressure, Heart Rate and EMG in Low Level Static Contractions. Acta Physiol. Scand. 1985, 125, 265–275. [Google Scholar] [CrossRef] [PubMed]
  44. Chapman, A.R.; Vicenzino, B.; Blanch, P.; Hodges, P.W. Patterns of Leg Muscle Recruitment Vary between Novice and Highly Trained Cyclists. J. Electromyogr. Kinesiol. 2008, 18, 359–371. [Google Scholar] [CrossRef]
Figure 1. AthosTM unit anterior view (a) and posterior view (set of contacts) (b). Exterior right leg (c) and interior left leg (d) view of sEMG shorts. Note: sEMG dry electrodes and electrode leads are composed of an inkjet-printed conductive polymer comprising an ether-based conductive thermoplastic polyurethane material. The electrodes are overlaid with a soft conductive silicone, which increases the stability of the electrode–skin interface.
Figure 1. AthosTM unit anterior view (a) and posterior view (set of contacts) (b). Exterior right leg (c) and interior left leg (d) view of sEMG shorts. Note: sEMG dry electrodes and electrode leads are composed of an inkjet-printed conductive polymer comprising an ether-based conductive thermoplastic polyurethane material. The electrodes are overlaid with a soft conductive silicone, which increases the stability of the electrode–skin interface.
Sensors 23 06998 g001
Figure 2. Movements for the sEMG calibration protocol to establish sEMG amplitude thresholds. Movements include prone knee flexion (a), prone hip extension (b), seated knee extension (c), and supine leg raise (d).
Figure 2. Movements for the sEMG calibration protocol to establish sEMG amplitude thresholds. Movements include prone knee flexion (a), prone hip extension (b), seated knee extension (c), and supine leg raise (d).
Sensors 23 06998 g002
Figure 3. Boxplot showing the distribution of sEMG-TL across different running speeds. sEMG-TL = surface electromyography training Load; a.u. = arbitrary units; low = 8 km⋅h−1; mod = 10 km⋅h−1; high = 12 km⋅h−1; black line = median, and black dots = 1 participant with a very high sEMG-TL.
Figure 3. Boxplot showing the distribution of sEMG-TL across different running speeds. sEMG-TL = surface electromyography training Load; a.u. = arbitrary units; low = 8 km⋅h−1; mod = 10 km⋅h−1; high = 12 km⋅h−1; black line = median, and black dots = 1 participant with a very high sEMG-TL.
Sensors 23 06998 g003
Figure 4. Scatterplot matrix of associations between variables: % V ˙ O2max = percentage of maximum oxygen uptake; sEMG-TL = surface electromyography training load; RPE = rating of perceived exertion; HR = heart rate; PL = PlayerLoad. Each blue dot corresponds to individual measurements at each one-minute stage during the treadmill running test. Solid lines are the least-squares derived best-fitting lines. * p < 0.05; ** p < 0.01.
Figure 4. Scatterplot matrix of associations between variables: % V ˙ O2max = percentage of maximum oxygen uptake; sEMG-TL = surface electromyography training load; RPE = rating of perceived exertion; HR = heart rate; PL = PlayerLoad. Each blue dot corresponds to individual measurements at each one-minute stage during the treadmill running test. Solid lines are the least-squares derived best-fitting lines. * p < 0.05; ** p < 0.01.
Sensors 23 06998 g004
Table 1. Participant characteristics.
Table 1. Participant characteristics.
CharacteristicFemaleMale
Age (years)23.5 ± 1.526.1 ± 3.7
Body mass (kg)59.3 ± 3.678.3 ± 9.4
Height (m)165.8 ± 8.0180.4 ± 8.3
Body fat (%)30.9 ± 6.317.4 ± 5.5
V ˙ O2max (mL⋅kg −1⋅min−1)40.8 ± 2.649.6 ± 6.0
Note: sEMG-TL = surface electromyography training load; RPE = rating of perceived exertion; HR = heart rate; PL = PlayerLoad. Values are presented as mean ± SD.
Table 2. Medians (IQR) for variables during each minute of the treadmill V ˙ O2max test.
Table 2. Medians (IQR) for variables during each minute of the treadmill V ˙ O2max test.
No. of
Participants
Minute   of   V ˙ O2maxsEMG-TL (a.u)HR (bpm)RPE (a.u.)PL (a.u.)
10185.85 (63.77, 100.37)138 (128, 156)8 (7, 12)15.00 (12.00, 16.75)
10281.42 (69.40, 105.46)166 (151, 170)11 (9, 14)15.00 (12.75, 16,25)
10382.79 (71.40, 108.62)163 (159, 178)13 (11, 14)15.50 (13.55, 16,25)
10487.07 (71.96, 109.15)167 (159, 176)15 (11, 17)15.50 (13.50, 16.25)
9578.44 (72.76, 107.05)170 (164, 181)14 (12, 17)16.00 (12.00, 17.00)
7694.26 (72.30, 109.51)177 (170, 194)16 (15, 19)16.00 (12.00, 17.00)
77102.89 (74.31, 112.42)181 (152, 192)16 (9, 19)16.00 (13.00, 17.00)
6894.52 (68.57, 110.88)182 (164, 186)18 (13, 19)16.50 (12.75, 18.75)
39112.69 (94.35, 116.31)186 (180, 189)18 (17, 18)17.00 (12.00, 18.50)
110113.152042018.00
111101.002042019.00
Note: sEMG-TL = surface electromyography training load; RPE = rating of perceived exertion; HR = heart rate; PL = PlayerLoad.
Table 3. Individual correlations between variables for each participant.
Table 3. Individual correlations between variables for each participant.
Participant 12345Participant12345
1PL 2
sEMG-TL0.33 −0.09
% V ˙ O2max0.420.83 ** −0.570.80
RPE0.460.81 **0.98 ** −0.380.92 *0.96 **
HR0.480.78 *0.98 **0.97 ** −0.760.610.94 *0.82
3PL 4
sEMG-TL0.49 0.67
% V ˙ O2max0.420.81 * 0.490.92 **
RPE0.400.670.95 ** 0.620.99 **0.96 **
HR0.430.84 **0.99 **0.95 ** 0.380.81 *0.96 **0.87 **
5PL 6
sEMG-TL−0.07 0.66 *
% V ˙ O2max−0.160.90 ** 0.85 **0.84 **
RPE−0.040.99 **0.89 ** 0.79 **0.88 **0.99 **
HR−0.150.94 **0.99 **0.93 ** 0.82 **0.74 **0.98 **0.96 **
7PL 8
sEMG-TL0.30 0.91 *
% V ˙ O2max0.310.97 * 0.97 **0.96 **
RPE0.000.560.73 0.96 **0.830.90 *
HR0.340.950.99 **0.77 0.95 *0.99 **0.99 **0.88
9PL 10
sEMG-TL−0.60 0.88 **
% V ˙ O2max−0.83 *0.85 ** 0.740.88 **
RPE−0.74 *0.96 **0.95 ** 0.750.88 **0.99 **
HR−0.83 *0.86 **0.99 **0.96 ** 0.79 *0.95 **0.97 **0.97 **
Note: % V ˙ O2max = percentage of maximum oxygen uptake; sEMG-TL = surface electromyography training load; RPE = rating of perceived exertion; HR = heart rate; * p < 0.05; ** p < 0.01.
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

Ashcroft, K.; Robinson, T.; Condell, J.; Penpraze, V.; White, A.; Bird, S.P. An Investigation of Surface EMG Shorts-Derived Training Load during Treadmill Running. Sensors 2023, 23, 6998. https://doi.org/10.3390/s23156998

AMA Style

Ashcroft K, Robinson T, Condell J, Penpraze V, White A, Bird SP. An Investigation of Surface EMG Shorts-Derived Training Load during Treadmill Running. Sensors. 2023; 23(15):6998. https://doi.org/10.3390/s23156998

Chicago/Turabian Style

Ashcroft, Kurtis, Tony Robinson, Joan Condell, Victoria Penpraze, Andrew White, and Stephen P. Bird. 2023. "An Investigation of Surface EMG Shorts-Derived Training Load during Treadmill Running" Sensors 23, no. 15: 6998. https://doi.org/10.3390/s23156998

APA Style

Ashcroft, K., Robinson, T., Condell, J., Penpraze, V., White, A., & Bird, S. P. (2023). An Investigation of Surface EMG Shorts-Derived Training Load during Treadmill Running. Sensors, 23(15), 6998. https://doi.org/10.3390/s23156998

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