Skip to main content

Are there two forms of isometric muscle action? Results of the experimental study support a distinction between a holding and a pushing isometric muscle function

Abstract

Background

In isometric muscle function, there are subjectively two different modes of performance: one can either hold isometrically – thus resist an impacting force – or push isometrically – therefore work against a stable resistance. The purpose of this study is to investigate whether or not two different isometric muscle actions – the holding vs. pushing one (HIMA vs PIMA) – can be distinguished by objective parameters.

Methods

Ten subjects performed two different measuring modes at 80% of MVC realized by a special pneumatic system. During HIMA the subject had to resist the defined impacting force of the pneumatic system in an isometric position, whereby the force of the cylinder works in direction of elbow flexion against the subject. During PIMA the subject worked isometrically in direction of elbow extension against a stable position of the system. The signals of pressure, force, acceleration and mechanomyography/-tendography (MMG/MTG) of the elbow extensor (MMGtri/MTGtri) and the abdominal muscle (MMGobl) were recorded and evaluated concerning the duration of maintaining the force level (force endurance) and the characteristics of MMG-/MTG-signals. Statistical group differences comparing HIMA vs. PIMA were estimated using SPSS.

Results

Significant differences between HIMA and PIMA were especially apparent regarding the force endurance: During HIMA the subjects showed a decisively shorter time of stable isometric position (19 ± 8 s) in comparison with PIMA (41 ± 24 s; p = .005). In addition, during PIMA the longest isometric plateau amounted to 59.4% of the overall duration time of isometric measuring, during HIMA it lasted 31.6% (p = .000). The frequency of MMG/MTG did not show significant differences. The power in the frequency ranges of 8–15 Hz and 10–29 Hz was significantly higher in the MTGtri performing HIMA compared to PIMA (but not for the MMGs). The amplitude of MMG/MTG did not show any significant difference considering the whole measurement. However, looking only at the last 10% of duration time (exhaustion), the MMGtri showed significantly higher amplitudes during PIMA.

Conclusion

The results suggest that under holding isometric conditions muscles exhaust earlier. That means that there are probably two forms of isometric muscle action. We hypothesize two potential reasons for faster yielding during HIMA: (1) earlier metabolic fatigue of the muscle fibers and (2) the complexity of neural control strategies.

Peer Review reports

Background

The human organism basically has two operation modes of muscle action: the static and the dynamic mode [1]. The dynamic working method can be categorized into concentric and eccentric muscle action. The static mode corresponds to the isometric muscle action [1,2,3]. Within concentric muscle action, we are able to overcome a resistance. In contrast, the eccentric muscle action enables us to decelerate an object. While isometric muscle action includes no gross joint motion, nevertheless, energy is consumed executing it. During performance of isometric muscle action in pilot studies, we have observed that two types of isometric muscle action can be executed. One can resist an object or one can push against it – although in both cases no motion is carried out.

The literature regarding different forms of isometric muscle action is very poor. Garner et al. [4] suggest that two forms of isometric action could exist based on different activity of electromyography (EMG) during concentric and eccentric muscle action. They hypothesize that the eccentrically loaded isometric muscle action (i.e., holding isometric muscle action (HIMA)) shows lower EMG-amplitude than the concentrically loaded one (i.e., pushing isometric muscle action (PIMA)). However, their readings of EMG-activity during muscle contractions up to 50% of the MVC of the soleus and gastrocnemius muscles showed no differences between HIMA and PIMA. Therefore, they reject the hypothesis and postulate that the modes of isometric muscle action are identical. Semmler et al. [5] also considered the possibility of two forms of isometric muscle action, but also found no statistical differences in EMG-activity.

Enoka and Rudroff et al. [6,7,8] used the term “force task” to describe a subject pulling or pushing actively against a stable structure (PIMA). A “position task” was defined as maintaining a constant position of the freely movable limb, while holding an inertial load (equivalent to the load during the “force task”), which represents HIMA. In summary, the two forms of isometric muscle action represents a position task (HIMA) that requires position control and a force task (PIMA) that requires the control of force.

Previous research suggests some objective differences between both hypothesized isometric actions. The most commonly found distinction is the “time to task failure” (i.e., the duration of maintaining either the position or the force task). Hunter et al. [9] performed experiments measuring fatigue of the elbow flexor muscles at 15% of the MVC in both of the isometric control tasks. The measurements showed a significantly longer duration time during the force task (PIMA) compared to the position task (HIMA). These findings were partially confirmed by Rudroff et al. [6, 7]. Similar results were only observed at lower intensities (20 and 30% MVC) [6, 7], when the arm was positioned horizontally. With higher force levels (45 and 60% of MVC) no differences between HIMA and PIMA were noted [7]. When the forearm was positioned vertically, no differences were found at any intensities [6].

Furthermore, the average EMG (AEMG) at exhaustion was greater during the force task (PIMA) compared to the position task (HIMA), despite a comparable rate of AEMG increase in both tasks. In contrast, the mean arterial pressure was higher at exhaustion during HIMA, while heart rate was not different [9].

Rudroff et al. [7] found significant differences in the spectral analysis of the power in the frequency band of 10 to 29 Hz. This bandwidth could be of particular interest, as the physiological frequencies of MMG-oscillations are found around 10 Hz. Thus, this could be another item to distinguish the hypothesized two forms of isometric muscle action.

Furthermore, Rudroff et al. [8] found a higher glucose uptake during HIMA compared to the PIMA in young men but not in older men.

Given that inconsistent findings regarding the two forms of isometric muscle action (HIMA vs. PIMA) we sought to examine this issue with a different methodological approach, which probably mimics real life situations better. Because there is freedom to vary position and force simultaneously in real life setting, our methodology aim to replicate that freedom.

It is hypothesized that two different forms of isometric muscle action (HIMA and PIMA) can be distinguished from each other based on the force endurance, the amplitudes of the mechanical myofascial oscillations (recorded by mechanomyography/-tendography; MMG/MTG) and the power spectral density (PSD). We hypothesize, due to previous research of Enoka et al. and pilot investigations at our department, that during PIMA the duration time is significantly longer than during HIMA. Furthermore, the mean amplitude of the whole duration time and the last 10% of duration time as well as the PSD of the MMG/MTG are hypothesized to differ significantly between HIMA and PIMA. Since it has been shown that EMG and MMG can synchronize and are therefore interdependent (e.g., [10, 11]), and furthermore Garner et al. [4] as well as Semmler et al. [5] did not found differences concerning the frequency of EMG comparing HIMA and PIMA, it is not expected to find differences in the mean frequency of the MMG/MTG either.

Methods

Participants

Ten healthy subjects (m = 5, f = 5) volunteered to participate in the study. The triceps brachii muscle on their dominant side was measured. All subjects were students of the University of Potsdam (studying sports therapy or sport for teaching profession), except for one participant, who was a high school student of age 18. Female subjects were aged averaged 24.4 years (±1.95), weighed 59.94 kg (±3.87) and were 169 cm tall (±1.58). They reached an averaged maximal voluntary isometric force (MVC) of 18.25 Nm (±1.96). The five male subjects were averagely 24 years old (±4.95), weighed 71.8 kg (±5.93) and were 182 cm tall (±3.81). They reached an averaged MVC of 31.166 Nm (±9.78). All subjects were right-handed. The study was conducted according to the declaration of Helsinki and local ethical permission was given. All subjects were informed in detail and gave their informed written consent to participate.

Explanation of the terms holding and pushing isometric muscle action

The Holding Isometric Muscle Action (HIMA) designates the mode of isometry, during which the subject shall resist an applied external load in a static position. Accordingly, the subject has to apply as much force as is exerted from the exterior. Consequently, the position of holding isometric muscle action can only be maintained if the external load does not exceed the individual maximal force. If supramaximal forces are applied to the subject, the muscle action merges from isometric into eccentric muscle action. The Pushing resp. Pulling Isometric Muscle Action (PIMA) describes the mode of isometric muscle action, during which the subject pushes against or pulls at a resistance. However, the produced force does not overcome the given resistance. Since a constant muscle length is maintained during both modes, a lengthening (during HIMA) or shortening (during PIMA) of the muscle is prevented.

Realization of the holding and pushing isometric muscle action

The holding and the pushing isometric muscle actions were realized by means of the pneumatic system. Due to the compressibility of air the subject can push isometrically against the push rod to perform the PIMA. Thereby, the subject has to control the force and the position, since the push rod yields if the adjusted force is overcome by the subject, or extends, if the subject gives in. The pneumatic system is adjusted in a way that the pressure amounts to 80% of the pressure at the MVC of the subject when the push rod is pushed halfway into the cylinder. If the participant’s force is not increased or decreased, the pressure and the push rod stay stable. Thus, the subject’s task is to first push against the push rod until 80% of the MVC is reached and then to stay in this position.

However, the push rod can also push against the subject with a defined pressure. Again due to the compressibility of air, the subject is able to hold the push rod in a stable position, if the subject maintains just as much force as the pneumatic system applies externally. Thus, the maximal force of the subject must not be overcome. If the subject yields, the push rod moves out of the cylinder and if the subject pushes against the push rod, it is pressed into the cylinder and consequently in both situations the force, as well as the position, are no longer constant. The pneumatic system is adjusted individually in a way, that the subject reaches 80% of its MVC when the push rod is extended halfway. If the subject yields slightly, the pressure tends to decrease first, but is then readjusted by the system and increases automatically until 80% of the MVC are reached again. The instructed task for the subject is, to hold the out-coming push rod in the given position of 90° elbow flexion.

Setting

The setup was constructed especially for the measurements of the two isometric muscle actions (Fig. 1). The subject sits on a chair, while the angles between the subject’s leg and trunk, arm and trunk as well as the elbow angle measure 90°. The subject has contact to the pneumatic system via an interface, which is placed proximal to the ulnar styloid process. It was made sure, that the positioning was kept stable and reproducible during the whole measuring period. The interface consists of a shell made of a thermal deformable polymer material, which is commonly used in rehabilitation technology. Between the interface and push rod of the pneumatic cylinder a strain gauge is fixed (model: MLMZ 2000 N 36), in order to record the reaction force between the subject and the pneumatic system. Furthermore, an acceleration sensor (comp. Biovision) is attached to the strain gauge to detect the acceleration along the longitudinal acting force vector. In order to record the mechanic muscle oscillations, piezoelectric MMG-sensors (model: Shadow SH 4001) are applied painlessly with tape to the lateral head of the triceps brachii muscle, its tendon and the abdominal external oblique muscle. Due to the transversal muscle oscillation, the lateral positioning of the sensors on the skin above the muscle belly has an influence on the signal [12,13,14,15]. This fact, as well as the difficulty to standardize the contact pressure, has to be considered. That is one reason for the limited comparability of amplitudes between the subjects. Nevertheless, the amplitudes are comparable intrapersonally, because the sensors were not detached between the measurements of HIMA and PIMA. The MMG-signals are conducted across an amplifier (Nobels preamp booster pre-1) to an A/D-converter (comp. Biovision/National Instruments, 14-bit) and subsequently are recorded by the software NI DIAdem 10.2 (National Instruments on Sony Vaio: PCG- 61111 M, Windows 7). The sampling rate was 1000 Hz, therefore, the Shannon theorem was fulfilled. The subjects were able to control the produced pressure level via biofeedback at the beginning of the measurement on the laptop. Afterwards, the rater gave oral information, if the pressure level deviated.

Fig. 1
figure 1

Setting for measuring HIMA and PIMA with the pneumatic system. The MMG/MTG of triceps brachii muscle and its tendon as well as the MMG of the abdominal external oblique muscle are recorded with piezo-sensors. The pressure is regulated with the control unit. All signals pass the A/D-converter and are recorded with DIAdem 10.2 on the measuring laptop

Measuring procedure

In total 12 measurements were done: At first, two measurements with maximal intensity were performed against a stable resistance to determine the MVC of the subject. After this, each subject performed ten more measurements. The first six measurements were made for 15 s at 80% of the MVC. Thereby, three consecutive measurements for HIMA and PIMA were completed. The initial muscle action was randomized. Then the three measurement of each task were done in blocked order. The resting time between the trials was 60s. The endurance measurements followed. During this task, the subject had to maintain the position at 80% of the MVC for as long as possible performing either the HIMA or the PIMA. The initial muscle action was again randomized, followed by alternating HIMA and PIMA tasks. In this way, the influential factor of fatigue was minimized. The resting period was 120 s between endurance trials.

The abort criterion for all test trials was a deviation of the push rod of more than 3 cm from the middle position. Based on subjects’ varying arm lengths (range = 24.5…30.0 cm), 3 cm represents a change of 1.27°. This was considered to be irrelevant here. During the HIMA any change of the push rod in direction into the cylinder (pushing against it) was a failure criterion, because consequently the muscle action has turned from HIMA into PIMA. Failed trials were repeated.

Data processing

The MMG-/MTG-, force-, pressure- and ACC-signals were recorded with the software DIAdem 10.2 (National Instruments) for all trials and were used for different considerations. The values of the strain gauge and of the pressure sensor during the given intensities were compared between the HIMA and PIMA to ensure the same force/pressure level. Torque was not calculated from the strain gauges’ force values as torque behaves proportionally to the force performing intraindividual comparisons.

First of all, the isometric intervals of all raw signals were cut for further data analysis. They were excised based on the pressure signal (Fig. 2).

Fig. 2
figure 2

Exemplary curve of pressure signal during HIMA including the isometric parts. The different parts of the pressure signal were used inter alia for the cut of the other signals. The following parts can be differentiated: the pressure increase (gray), the total isometric phase (black), single isometric phases (orange) and the longest single isometric phase (red). The serrated deviations of the line show moments of jerky yielding

Isometric phases of the pressure signal for analyzing the force endurance

From the moment of reaching the given pressure level, the first deviation of the pressure signal (marked by a short visible pressure loss) was set at the starting point of the (first) isometric phase. The end of the first single isometric phase was appointed to the moment immediately before the second jerky yielding of pressure. These were so apparent in the pressure signal that a visual observation suffices. The next single isometric phase began immediately after the second deviation and this was continued (Fig. 2) until the subject falls below the given pressure level. Phases shorter than 1 s were not taken into account, because stable isometric condition were not built up. The sum of the duration time of all single isometric phases results in the total isometric phase. Furthermore, the longest single isometric phase of all signals was also used for specific considerations. Further preparations of data differ depending on the subsequent analysis.

Frequency of MMG-/MTG-signals

To evaluate the frequency of the MMG- and MTG-signals different methods were used: 1.) Power Spectral Density (Software DIAdem), 2.) Continuous Wavelet Transform (using Software Python), 3.) a Python script, which calculates the average of the time intervals between the maxima [16]. For 1.) and 2.) raw data were used. For 3.) the signal was filtered with a Butterworth filter (filtering degree 5, cut-off frequency: 20 Hz).

Furthermore, the power in the frequency range of 8 to 15 Hz, which is the commonly known physiological frequency of MMG, and the frequency band of 10 to 29 Hz, was investigated. The latter was considered by Hunter et al. [9] concerning the EMG power.

Amplitude of MMG-/MTG-signals

The data of the signals were filtered with a Butterworth filter (filtering degree 5, cut-off frequency: 20 Hz). Afterwards the arithmetic mean values of the amplitude of the total isometric phase and of the longest isometric phase were calculated using Excel (Microsoft Office). For another research question, just the last 10% of the duration time of the total isometric phase of the endurance measurements were included. This consideration of amplitude at exhaustion was done by Hunter et al. [9] concerning the AEMG. Regarding this analysis, two test persons had to be excluded, because their measurements at PIMA were stopped by the rater prior to complete exhaustion due to prolonged trial duration. For a further consideration, the coefficient of variation (CV) of the mean amplitude of the MMG-/MTG-signals between the three 15 s-measurements during the HIMA and the PIMA was calculated.

Group statistics

All group comparisons between HIMA and PIMA were done with SPSS Statistics (IBM). The significance level was set at α = 0.05 . Firstly, the data were tested by the Kolmogorov-Smirnov-Test regarding their normal distribution. If the normal distribution was fulfilled, the paired t-test was utilized. Otherwise, the comparison was done with the Wilcoxon-Test for paired samples.

Since we included several items in the testing procedure, the problem of multiple testing could be of relevance. This topic is controversially discussed in the literature (e.g., [17,18,19]). Rothman [20] claims that there are no adjustments needed for multiple comparisons. We basically support this opinion and agree with the conclusion of O’Brien [19]. Nevertheless, the problem of multiple comparisons seems to be considered particularly in clinical trials [21]. Since the present study is not a clinical trial, even proponents of adjustments for multiple comparisons might agree with our decision to set no adjustments here.

Results

Figure 3 displays exemplar raw signals of the MMG of the triceps brachii muscle during measurements at 80% of the MVC in the HIMA vs. PIMA.

Fig. 3
figure 3

Exemplary signals of MMG of the triceps muscle. Displayed are the isometric phases of offset raw MMG signals of the triceps muscle during the holding isometric muscle action (HIMA; above) and during the pushing isometric muscle action (PIMA; below) at 80% of the MVC. The zoomed areas show a one-second sequence

Force levels

The force level of the 15 s-measurements during HIMA and PIMA shows no significant difference. The mean force value of the holding isometric muscle action amounts to 99.31% of the PIMA. The mean force value of the strain gauge during HIMA was 55.19 N (±18.78) (men (m) = 67.18 N (±20.07), female (f) = 43.21 N (±5.65)). During the PIMA the arithmetic mean value of the strain gauge was 56.25 N (±20.94) (m = 69.24 N (±22.53), f = 43.25 N (±7.51)). The signals were distributed normally. The paired t-test was not significant (t(9) = 0.860; p = .412).

During the endurance measurements the force levels of HIMA resp. PIMA differ, although the pressure level does not. The mean force level of the n = 10 subjects during the HIMA endurance trials (55.14 N (±18.83)) amounts to 96.35% of the force level during the PIMA (57.23 N (±19.12); t(9) = 3,087; p = .013).

Amplitude

The mean amplitudes of the MMG and MTG signals did not differ significantly between PIMA and HIMA (range of p = 0.069…0.765) regarding the total isometric measuring time, the longest isometric phase during the 15 s-measurements as well as during the endurance measurements.

Considering only the last 10% of the time duration of the endurance trials (exhaustion), the amplitude of the MMG of the triceps brachii muscle showed significantly higher amplitudes at exhaustion during PIMA comparing to HIMA (n = 8; W = −2.521, p = .012). The MMGobl and MTGtri did not show any substantial differences.

During PIMA the CV of the amplitude of the MMGobl between the three 15 s-measurements (n = 10) was significantly higher comparing to HIMA (t(29) = 2523, p = .017) (Fig. 4). In eight of ten subjects the CV of mean amplitude of MMGobl (not in MMGtri or MTGtri) is higher during the PIMA.

Fig. 4
figure 4

Coefficient of variation (CV) of the mean amplitude of MMGobl. Illustrated is the CV of the MMG amplitude of the abdominal external oblique muscle between the three 15 s-measurements of all n = 10 subjects during HIMA (orange) and PIMA (blue) sorted by the difference of the CV between PIMA and HIMA (dashed line)

Frequency

Figure 5 shows the exemplary wavelet spectra of MMG and MTG-signals of one measurement. Since the HIMA and PIMA did not show any differences concerning these spectra, only the three spectra of one isometric mode are illustrated. The mean frequency shows no significant differences in the paired t-Test between the PIMA and the HIMA (p > .05).

Fig. 5
figure 5

Schematic diagram of the Wavelet Spectra of the MMG/MTG-signals. The Wavelet spectra of the MMG resp. MTG of the triceps brachii muscle and its tendon (top left/right) as well as of the abdominal external oblique muscle (bottom left) during isometric muscle action at 80% of the MVC show the stochastically distributed oscillations in the low frequency area around approximately 15 Hz. The right scale displays the intensity of power. The left side of each figure displays the mean of the power over the time. This can be seen analogues to the Power Spectral Density (PSD)

Power

Looking at the specific frequency ranges of 8 to 15 Hz and 10 to 29 Hz, respectively, the MTG-signal of the tendon of the triceps brachii muscle shows higher power values in HIMA than in PIMA. In the frequency range of 10 to 29 Hz the t-test for paired samples displays a significant difference (t(9) = −2.284, p = .048). This frequency range was chosen following Hunter et al. [9]. Looking at the common frequency range of MMG of 8 to 15 Hz, the significance becomes even clearer with p = .037 resulting from the Wilcoxon-Test for paired samples (power values were non-parametrically). The signals of the triceps brachii muscle as well as the abdominal external oblique muscle do not show any differences in either frequency ranges.

Endurance time

Using the pneumatic system, the single isometric phases of the duration measurements can be exactly identified, because any deviation of the push rod can be seen in the pressure signal (Fig. 2). Table 1 Therefore, the total isometric phase consists of several shorter single isometric phases. Table 1 shows the parameters’ raw values of the endurance trials. As can be seen in the last line and in Fig. 6, all items show significant differences between the PIMA and HIMA, respectively. During the pushing isometric muscle action, the subjects were able to stabilize a pure isometric position of the push rod (longest single isometric phase – without any deviation of the push rod) for approximately 41.4 s, which amounts to 59.4% of the total time duration of the isometric plateau (Fig. 7). During the holding isometric muscle action, the longest single isometric phase is 19.1 s, which relates to 31.6% of the total duration time. Thus, the pure isometric position in the PIMA could be maintained twice as long as in the HIMA.

Table 1 Overview of endurance parameters of the fatiguing measurements comparing PIMA vs. HIMA
Fig. 6
figure 6

Items of the endurance measurements comparing the PIMA (blue) vs. HIMA (gray). The bars display the duration time [s] of different items of maintaining 80% of MVC comparing the pushing isometric muscle action (PIMA; blue) vs. the holding isometric muscle action (HIMA, gray). *p < .05, ** p < .01

Fig. 7
figure 7

Relation between longest isometric phase and total isometric phase. Displayed are the arithmetic mean and the standard deviation comparing PIMA (blue) and HIMA (gray). ***p < .000

Discussion

Content-related discussion

The present results provide findings, which partially support the hypothesis of two different forms of isometric muscle action. Overall, the authors assume the findings to be sufficient to distinct between the holding and the pushing isometric muscle action (HIMA vs. PIMA), respectively. Nevertheless, one also has to consider, that there are some items, which do not support the differentiation between such two forms.

The parameters, which show statistical differences between HIMA and PIMA, are the force endurance, the power in the frequency range of 8 to 15 Hz and 10 to 29 Hz of the MTG of the triceps tendon, the mean amplitude of the MMG of the triceps muscle at exhaustion and potentially the higher variability of the MMG of the abdominal external oblique. Therefore, the findings are generally consistent with them of Enoka and Rudroff et al. [6, 7, 9] as in the present study a longer force endurance appeared during PIMA. However, in contrast to their findings, the present study shows a significant difference between the HIMA and the PIMA regarding the endurance time even under higher intensities (80% of MVC) and with a vertically positioned forearm. One has to consider primarily, that their measurements were made in a different setting (no use of pneumatics, different forearm position), which could explain the different outcomes.

The oscillations of the MMG/MTG show inconclusive results regarding the hypothesis. The results do not display a statistically relevant difference between both isometric muscle actions concerning their frequency. This was expected based on the findings of other investigators (e.g., [4, 5]), who did not find differences concerning the frequency of the EMG.

The main question for us to discuss is how the lower force endurances during the HIMA, compared to the PIMA, can be explained. This also includes the findings concerning the oscillation characteristics (frequency and amplitude), since they may be interdependent. In general, different physiological processes could be responsible for faster muscular fatiguing. Weineck [22] lists several reasons, for example a central fatigue or a fatigue caused by reduced hormone production. These two possibilities can be excluded due to the short measuring time. The measurement duration also speaks against a faster synaptic fatigability in the HIMA. Furthermore, studies have shown that the EMG-activity is lower during eccentric muscle action (e.g., [23,24,25]). Assuming that the HIMA is closer to the eccentric muscle action (see below), this would also speak against a synaptic fatiguing. Thus, two potential sources of the faster yielding at submaximal static force endurance remain:

  1. (1)

    Metabolic fatigue of the muscle fibers (e.g., oxygen supply, removal of lactate)

  2. (2)

    Complexity of the neural control strategies.

The subsequent discussion is based on the hypothetical idea, that during the holding isometric muscle action (HIMA) the neuronal control strategies and muscle physiological aspects could be similar to the lengthening mode – analog to the eccentric activity – but without performing a change in length or joint angle. Vice versa, the mechanisms of the pushing isometric muscle action (PIMA) could be similar to the shortening mode – comparable to the concentric activity – without making a change in length or joint angle. On the one hand, this hypothesis is based on the subjective feeling during the performance of either the HIMA or the PIMA. On the other hand, this assumption could also be supported by the findings of Hunter et al. [9], who found greater AEMG at exhaustion during PIMA compared to HIMA. It is well known that during concentric muscle action, the amplitude of EMG is greater than during the eccentric one. So perhaps this could be a link between the proposed isometric muscle actions and the common eccentric and concentric ones. Garner et al. [4] also proposed this hypothesis. Furthermore, the study of Grabiner et al. [26] showed a higher muscle activation (EMG) if the subject expects a concentric contraction compared to if it anticipates an eccentric one. These different a priori activations of the CNS might have also occurred in the current investigation. The subjects were instructed to either hold the external resistance without allowing the push rod to overcome them (i.e., they were expecting an external resistance, which could force the subject into eccentric action, if the impacting force was not resisted sufficiently). Or the subjects were instructed to push against the push rod, which included the task to contract the muscle as in concentric muscle action, but not more than to a specific level of force. There are different approaches, which could support the hypothesis that HIMA is closer to eccentric muscle action, while PIMA corresponds to the concentric one instead. The following discussion grounds on this hypothesis.

Metabolic fatigue of the muscle fibers as reason for lower force endurance during HIMA

For the first mentioned aspect as a potential reason for the lower force endurance during HIMA – the metabolic fatigue (1) – a local lack of metabolic exchange under the conditions of intramuscular compression during isometric action could be presumed. This consideration is speculative as the current investigation did not measure substrate transport.

According to several authors [28,29,30] an ischemia of muscles can occur at intensities as low as 5–30% of MVC because the muscle is compressing its own blood vessels. Jensen et al. [31] reported a reduction of oxygenation at 30–40 mmHg, which corresponds to a contraction level of 20% of MVC. With intensities of 80% of MVC it is most likely that the blood vessels are also compressed, resulting in a probable ischemia. In this case it is plausible that not only the oxygenation, but also the exchange of substrates might be limited or not possible at all. The inadequate oxygenation due to the muscles’ compression would rather implicate an anaerobic metabolism during the present muscle tasks. Thereby, inter alia, lactate would be produced and would increase during the isometric muscle action. According to Katz et al. [32] this can amount up to 15 times of the initial value of lactate. Under these tensed up conditions during static muscle action, one has to ask, how an exchange and transport of substrates can be enabled at all over an interval of 50–60s? And why does the duration of force maintenance differ between HIMA and PIMA, even though the generated force, and thus the compression of muscle fibers, should remain the same? One possible explanation on how the exchange and transport of substrates is maintained could be based on the neuromuscular oscillations.

In the present study, the amplitudes of the MMG of the triceps muscle show higher values at exhaustion during PIMA compared to HIMA. EMG and MMG are different, but causally linked, methods with similar behavior at exhaustion (i.e., amplitude increase at higher intensities [11, 33,34,35,36,37,38]). Thus, one could expect related results. The investigations of Rudroff et al., Hunter et al. and the present study all show higher values at exhaustion concerning the mean amplitude during PIMA compared to HIMA in the AEMG [6,7,8,9] and in the MMG of triceps muscle, respectively (present study). The amplitude of the MMG reflects the fine mechanic oscillations of the muscle fibers. We presume that they are possibly necessary to maintain any of the exchange and transport of substrates during isometric muscle action. Because under particular circumstances the transport of fluids can be supported by vibration [39, 40]. The muscle tissue may possibly use a similar mechanism during reduced flow of substrates due to compression. Therefore, we hypothesize that higher oscillations – equivalent to vibrations – in the MMG of the triceps brachii muscle during the PIMA would probably support the substrate flow between the cell and interstitium as well as the flow within these areas. Conversely, lactate would probably be removed slower, and therefore would accumulate more during HIMA. Even though the lactate production is lower during eccentric muscle action [41, 42], the ratio between production and removal could potentially be more detrimental during the HIMA. Thus, the muscular fatigue could begin sooner during the HIMA and consequently the subject could yield earlier. Further investigation including measurements of lactate or O2 remain.

In fact, the hypothesis of exchange and transport of substrate could be a fundamental physiological reason for the obligatory requirement of muscular oscillations during isometric action in general. This could be an area of future research.

In the present measurements, it was impressively apparent, that during HIMA the position could not be maintained stable by the subjects for as long as during PIMA. The holding position was interrupted earlier by sudden, quick, but short motions of the push rod (and thus, short-time lengthening of the muscle). But immediately after these abrupt releasing actions, the position could again be maintained at the same level of force, but with a slightly changed elbow angle. This was merely a subjective observation, but maybe the exchange and transport of substances could be expedited by the jerky yielding of muscle tension (compared to a very short intensive oscillation). In this way, the metabolic environment might be improved again for a short time.

In view of these considerations it has to be taken into account, that only the MMG of the triceps brachii muscle shows a higher mean amplitude at exhaustion. The other sensor locations do not. Since the abdominal external oblique muscle is probably only utilized to stabilize the trunk, the different tasks of HIMA and PIMA may not be reflected there.

The question remains, why the triceps tendon does not show the same output as its muscle? Even though both methods, the MMG and MTG, measure mechanic oscillations, the mechanisms behind these oscillations are different. The working muscle generates oscillations actively, whereas the tendon just resonates in a passive way during force transmission. From this point of view, it is conceivable that at higher intensities, the muscle (motor) vibrates more intensively. In contrast, a more tightened tendon (like a passive rope) would swing with lower amplitudes (and possibly higher frequencies).

Almost no basic knowledge exists about the MTG. Based on our investigation we know that the oscillations of the MTG and MMG show similar frequency ranges and are able to generate coherent behavior during isometric muscle action [43, 44]. The authors are not aware of any other studies about the MTG. Thus, the MTG as a method is not yet established in science. Since the same sensors are used for MMG and the tendon is a passive structure which directly is connected to the muscle, the oscillations of the muscle theoretically has to be transmitted to the tendon. This is what we saw in previous investigations [43, 44]. Anyhow, we suppose that the tendon of triceps muscle occupies a kind of superposition, since all three heads, and a large number of motor units of the triceps muscle come together and form the oscillations. This, in turn, could be an explanation for the different outcomes regarding the power of the MTG resp. MMG signals. The MTG of the triceps brachii muscle show significantly higher power in the frequency range of 8 to 15 Hz (as well as 10 to 29 Hz) during HIMA compared to PIMA. We remember that the mean amplitude of MMG of triceps brachii was higher during PIMA at exhaustion. Normally, higher amplitudes indicate higher power. In this case one may ask, how these findings can be brought together?

Complexity of the neural control strategies as a possible reason for lower force endurance during HIMA

Among others, Semmler et al. [5] have shown that the synchronization of motor units (EMG) is higher during eccentric muscle action than during the concentric one. Possibly – and still under the assumption that the HIMA is closer to the processes of eccentric muscle action – a higher synchronization would enable the muscle fibers to oscillate in a narrower frequency range than during the PIMA. This could be the reason for the higher power of the MTG of the triceps in the low frequency range of 8 to 15 Hz (and 10 to 29 Hz) during the HIMA. While the PIMA is potentially wider spread, including other frequencies, and therefore results in lower power in the specific frequency range of 8 to 15 Hz. Rudroff et al. [7] found contrary results concerning the EMG. Their findings show that the power during the force task (similar to PIMA) was higher in 10 to 29 Hz than during the position task and not vice versa as found here.

Fang et al. [25] investigated the movement-related cortical potential (MRCP) via EEG comparing eccentric and concentric muscle action. They found that – although the EMG in their study showed lower amplitude during eccentric muscle action – the MRCP still resulted in greater values. The investigators postulated that eccentric muscle action is more difficult to perform compared to the concentric one, and therefore results in the higher brain activity. These findings are supported by the review of Enoka [27], who summarized, that the eccentric muscle action involves more complex control strategies of the nervous system than the concentric one. This was underpinned by Duchateau & Baudry [45], who added, that the unique mechanisms of eccentric muscle action are still unknown.

Bringing all of this together, the first paradoxically appearing result of significantly higher power in small frequency ranges during HIMA in the MTG, and conversely the greater amplitude at exhaustion during PIMA in the MMG of the triceps muscle could rely on the assumption, that eccentric muscle action is based on more complex neural control strategies. Thus, the present investigation indicates that HIMA could be more difficult to perform compared to PIMA due to a more complex adjustment of the neuromuscular system. Indicators for this are the shorter endurance time and the higher power of the MTG oscillations in a small frequency range of 8 to 15 Hz as a hint for a greater synchronization of muscle activation during HIMA. In turn, this possibly influences the fatigue mechanism, but further research is needed to test this hypothesis.

Methodological limitations

In summary, we have to note two main methodological limitations within this study. Firstly, the force level shows a significant difference between HIMA and PIMA during the endurance measurements. We suppose that the reason for this lies in the stick and slip effect of the push rod. Since the push rod is not frictionless, jerky movements can occur during the slip out of the cylinder. Throughout the measurements, we controlled the force level of the subjects using the pressure signal, not the strain gauge’s signal. The reason for this was, that we had to use the pressure values to adjust the force of the pneumatic system individually for each subject to realize both measuring modes. While the pressure values do not show a statistical distinction between the HIMA and the PIMA, the values of the strain gauge do solely during the endurance measurements. Thereby, the stick and slip effect might appear to play a greater role when the piston is driving out of the push rod, which occurred more frequently during the endurance measurements compared to the 15 s trials. Therefore, prior to each yielding of the subject, the force level had to decrease. While the strain gauge captures this decline immediately, the push rod does not drive out instantly because of the friction. This may explain the different PIMA and HIMA pressure and force value behaviors. Anyway, the force needed to be maintained was lower during the HIMA, which should have made the task easier for subjects. As a result, this limitation may further support the results. A second limitation probably lies in identifying the maximal voluntary isometric contraction (MVC). As the setting shows, the subjects were fixed neither at their backs nor at their feet (Fig. 1). A higher MVC could possibly be generated with a different kind of setting, for instance by using a counterfort.

Conclusion

The results seem to be sufficient for the authors to partially support the hypothesis of a differentiation of the isometric muscle action into two forms: The holding vs. pushing isometric muscle action (HIMA vs. PIMA). The results suggest that under holding isometric conditions muscles exhaust earlier. Further investigations have to be done on these findings. If the results will be supported, it has to be discussed how these postulated – but in motor science still not established – two different forms of muscle action can be included in the current theoretical construct of muscle activity. Currently, measurements concerning HIMA and PIMA with maximal intensities are in progress at our department. Furthermore, to obtain more detailed information concerning the assumption that HIMA is closer to eccentric muscle action, we have conducted first pilot measurements using EEG and MMG during HIMA and PIMA. Perhaps the hypothesis, that the amplitude in the EEG is higher during HIMA – analogues to eccentric – can be supported by this investigation. The evaluation remains ongoing. All this will probably deliver further hints for the control mechanisms of these specific muscle actions and the integration into the construct of operating modes of musculature.

The supposed differentiation of the two forms could probably bring benefits for diagnostics, prevention, and rehabilitation in the future. Since the eccentric muscle action implies specific control strategies of the nervous system, the supposed HIMA might require them as well. It is known, that during eccentric muscle action the vulnerability of the musculoskeletal system is at its highest [46,47,48]. Since the eccentric muscle action is based on a more complex neural control strategy, this also emphasizes the importance of a functionally intact neuromuscular system, while performing eccentric muscle action. It has to be investigated, whether or not HIMA is more sensitive for disorders of the neuromuscular system compared to PIMA, because of the higher vulnerability of eccentric muscle action. If the differences between the holding and the pushing isometric mode could be underpinned by further investigations, this differentiation would possibly lead to a more specific clinical muscle testing, e.g., used for neurological diagnosis.

In this regard, the authors also see potential in measuring two interacting neuromuscular systems to compare the two isometric muscle actions with each other [44]. This requires an even higher sensomotoric control of the neuromuscular system. Thus, it could deliver additional information about the control strategies during the two isometric muscle actions. The interaction, inter alia, will depend on the ability to synchronize the mechanical oscillations of both neuromuscular systems [44]. Results of these measurements of interacting partners concerning both isometric muscle actions will be presented soon.

Abbreviations

ACC:

Acceleration sensor

CV:

Coefficient of variation

f:

Female

HIMA:

Holding isometric muscle action

m:

Men

MMG:

Mechanomyography

MMGobl:

Mechanomyography of the abdominal external oblique muscle

MMGtri:

Mechanomyography of triceps brachii muscle

MTG:

Mechanotendography

MTGtri:

Mechanotendography of the tendon of the triceps brachii muscle

MVC:

Maximal voluntary contraction

PIMA:

Pushing isometric muscle action

SD:

Standard deviation

References

  1. Wick D. Biomechanik sportlicher Bewegungen - Lehrbuch biomechanischer Grundlagen in Theorie und Praxis des Sports. Balingen: Spitta-Verlag; 2005.

    Google Scholar 

  2. Schnabel G, Harre HD, Krug J. Trainingslehre – Trainingswissenschaft. Leistung – Training – Wettkampf. 2nd ed. Aachen: Meyer & Meyer; 2011.

    Google Scholar 

  3. Komi PV. Kraft und Schnellkraft im Sport. Enzyklopädie der Sportmedizin. Köln: Deutscher Ärzte Verlag; 1994.

    Google Scholar 

  4. Garner JC, Blackburn T, Weimar W, Campbell B. Comparison of electromyographic activity during eccentrically versus concentrically loaded isometric contractions. J Electromyogr Kinesiol. 2008;18(3):466–71.

    Article  PubMed  Google Scholar 

  5. Semmler JG, Kornatz KW, Dinenno DV, Zhou S, Enoka RM. Motor unit synchronization is enhanced during slow lengthening contractions of a hand muscle. J Physiol. 2002;545(2):681–95.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Rudroff T, Barry BK, Stone AL, Barry CJ, Enoka RM. Accessory muscle activity contributes to the variation in time to task failure for different arm postures and loads. J Appl Physiol. 2007;102:1000–6.

    Article  PubMed  Google Scholar 

  7. Rudroff T, Justice JN, Holmes MR, Matthews SD, Enoka RM. Muscle activity and time to task failure differ with load compliance and target force for elbow flexor muscles. J Appl Physiol. 2011;110:125–36.

    Article  PubMed  Google Scholar 

  8. Rudroff T, Kalliokoski KK, Block DE, Gould JR, Klingensmith III WC, Enoka RM. PET/CT imaging of age- and task-associated differences in muscle activity during fatiguing contractions. J Appl Physiol. 2013;114:1211–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Hunter SK, Ryan DL, Ortega JD, Enoka RM. Task differences with the same load torque alter the endurance. J Neurophysiol. 2002;88:3087–96.

    Article  PubMed  Google Scholar 

  10. Coates MJ, Baker SN, Fitzgerald WJ. A complex wavelet-transform approach to power and coherence measurement from non-stationary data. J Physiol. 1997;501:36.

    Google Scholar 

  11. Barry DT, Geiringer S, Ball RD. Acoustic myography: a noninvasive monitor of motor unit fatigue. Muscle Nerve. 1985;8:189–94.

    Article  CAS  PubMed  Google Scholar 

  12. Beck T. Applications of Mechanomyography for Examining Muscle Function. Editor Travis. W. Beck; 2010.

  13. Cescon C, Madeleine P, Farina D. Longitudinal and transverse propagation of surface mechanomyographic waves generated by single motor unit activity. Med Biol Eng Comput. 2008;46:871–7.

    Article  PubMed  Google Scholar 

  14. Kaczmarek P, Celichowski J, Kasiński A. Experimentally verified model of mechanomyograms recorded during single motor unit contractions. J Electromyogr Kinesiol. 2005;15:617–30.

    Article  CAS  PubMed  Google Scholar 

  15. Farina D, Li X, Madeleine P. Motor unit acceleration maps and interference mechanomyographic distribution. J Biomech. 2008;41:2843–9.

    Article  PubMed  Google Scholar 

  16. Pikovsky A, Rosenblum M, Kurths J. Synchronization. A universal concept in nonlinear sciences. Cambridge: University Press; 2003.

    Google Scholar 

  17. Bender R, Langer S, Ziegler A. Multiples Testen. Dtsch Med Wochenschr. 2002;127:T4–7.

    Google Scholar 

  18. Proschan MA, Waclawiw MA. Practical guidelines for multiplicity adjustment in clinical trials. Control Clin Trials. 2000;21:527–39.

    Article  CAS  PubMed  Google Scholar 

  19. O’Brien PC. The appropriateness of analysis of variance and multiple comparison procedures. Biometrics. 1983;39:787–94.

    Article  PubMed  Google Scholar 

  20. Rothman KJ. No adjustments are needed for multiple comparisons. Epidemiology. 1990;1:43–6.

    Article  CAS  PubMed  Google Scholar 

  21. Bender R, Lange S. Adjusting for multiple testing – when and how? J Clin Epidemiol. 2001;54:343–9.

    Article  CAS  PubMed  Google Scholar 

  22. Weineck J. Sportbiologie. 9th ed. Balingen: Spitta Verlag; 2004.

    Google Scholar 

  23. Baudry S, Klass M, Pasquet B, Duchateau J. Age-related fatigability of the ankle dorsiflexor muscles during concentric and eccentric contractions. Eur J Appl Physiol. 2007;100:515–25.

    Article  PubMed  Google Scholar 

  24. Christou EA, Shinohara M, Enoka RM. Fluctuations in acceleration during voluntary contractions lead to greater impairment of movement accuracy in older adults. J Appl Physiol. 2003;95:373–84.

    Article  PubMed  Google Scholar 

  25. Fang Y, Siemionow V, Sahgal V, Xiong F, Yue GH. Greater movement-related cortical potential during human eccentric versus concentric muscle contractions. J Neurophysiol. 2001;86:1764–72.

    CAS  PubMed  Google Scholar 

  26. Grabiner MD, Owings TM, George MR, Enoka RM. Eccentric contractions are specified a priori by the CNS. Proc. VXth Congr. Int. Soc. Biomech. Jyva¨skyla¨ Finland July 2–6, 1995:338–9.

  27. Enoka RM. Eccentric contractions require unique activation strategies by the nervous system. J Appl Physiol. 1996;81:2339–46.

    CAS  PubMed  Google Scholar 

  28. Mense S. Neurobiologische Grundlagen von Muskelschmerz. Schmerz. 1999;13:3–17.

    Article  CAS  PubMed  Google Scholar 

  29. Järvholm U, Styf J, Suurkula M, Herberts P. Intramuscular pressure and muscle blood flow in supraspinatus. Eur J Appl Physiol. 1988;58:219–24.

    Article  Google Scholar 

  30. Sjøgaard G, Savard G, Juel C. Muscle blood flow during isometric activity and its relation to muscle fatigue. Eur J Appl Physiol. 1988;57:327–35.

    Article  Google Scholar 

  31. Jensen B, Jørgensen K, Hargens A, Nielsen P, Nicolaisen T. Physiological response to submaximal isometric contractions of the paravertebral muscles. Spine. 1999;24(22):2332.

    Article  CAS  PubMed  Google Scholar 

  32. Katz A, Sahlin K, Henriksson J. Muscle ammonia metabolism during isometric contraction in humans. Am J Physiol Cell Physiol. 1986;250(6):C834–40.

    CAS  Google Scholar 

  33. Akataki K, Mita K, Watakabe M, Itoh K. Mechanomyogram and force relationship during voluntary isometric ramp contractions of the biceps brachii muscle. Eur J Appl Physiol. 2001;84:19–25.

    Article  CAS  PubMed  Google Scholar 

  34. Orizio C, Perini R, Veicsteinas A. Muscular sound and force relationship during isometric contraction in man. Eur J Appl Physiol. 1989;58:528–33.

    Article  CAS  Google Scholar 

  35. Orizio C, Gobbo M, Diemont B, Exposito F, Veicsteinas A. The surface mechanomyography as a tool to describe the influence of fatigue on biceps brachii motor unit activation strategy. Historical basis and novel evidence. Eur J Appl Physiol. 2003;90:326–36.

    Article  PubMed  Google Scholar 

  36. Yoshitake Y, Moritani T. The muscle sound properties of different muscle fiber types during voluntary and electrically induced contractions. J Electromyogr Kinesiol. 1999;9:209–17.

    Article  CAS  PubMed  Google Scholar 

  37. Stokes MJ, Dalton PA. Acoustic myography for investigating human skeletal muscle fatigue. J Appl Physiol. 1991;71:1422–6.

    CAS  PubMed  Google Scholar 

  38. Evetovich T, Housh T, Stout J, Johnson G, Smith D, Ebersole K. Mechanomyographic responses to concentric isokinetic muscle concentric. Eur J Appl Physiol. 1997;75:166–9.

    Article  CAS  Google Scholar 

  39. Thiria B, Zhang J. Ratcheting fluid with geometric anisotropy. Appl Phys Lett. 2015;106:054106. doi:10.1063/1.4906927.

    Article  Google Scholar 

  40. John K, Thiele U. Kletternde Tropfen. Phys J. 2013;12(2):33–8.

    CAS  Google Scholar 

  41. Linnamo V, Bottas R, Komi PV. Force and EMG power spectrum durgin and after eccentric and concentric fatigue. J Electromyogr Kinesiol. 2000;10:293–300.

    Article  CAS  PubMed  Google Scholar 

  42. Durand RJ, Castracane VD, Hollander DB, Tryniecki JL, Bamman MM, O’Neal S, et al. Hormonal Responses from Concentric and Eccentric Muscle Contraction. Med Sci Sports Exerc. 2003;35(6):937–43.

  43. Schaefer L. Synchronisationsphänomene myotendinöser Oszillationen interagierender neuromuskulärer Systeme – mit Betrachtung einer Hypothese bezüglich unterschiedlicher Qualitäten isometrischer Muskelaktion. Doctoral dissertation, The University of Potsdam, Germany. 2014. Available: http://opus.kobv.de/ubp/volltexte/2014/7244/.

  44. Schaefer L, Torick A, Matuschek H, Holschneider M, Bittmann F. Synchronization of muscular oscillations between two subjects during isometric interaction. Eur J Trans Myol - Basic Appl Myol. 2014;24(3):195–202.

    Google Scholar 

  45. Duchateau J, Baudry S. Insights to the neural control of eccentric contractions. Review. J Appl Physiol. 2013;116:1418–25. doi:10.1152/japplphysiol.00002.2013.

  46. Kieb M, Lorbach O, Engelhardt M. Muskelverletzungen: diagnostik und behandlungen. Orthopade. 2010;39:1098–107.

    Article  CAS  PubMed  Google Scholar 

  47. Fousekis K, Tsepis E, Poulmedis P, et al. Intrinsic risk factors of non-contact quadriceps and hamstring strains in soccer: a prospective study of 100 professional players. Br J Sports Med. 2011;45:709–14.

    Article  PubMed  Google Scholar 

  48. Witchalls J, Blanch P, Waddington G, Adams R. Intrinsic functional deficits associated with increased risk of ankle injuries: a systematic review with meta-analysis. Br J Sports Med. 2011. doi:10.1136/bjsports-2011-090137.

    PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

No funding was used for this study or the production of this article.

Availability of data and materials

The datasets used and/or analyzed during the current study available from the corresponding author on reasonable request.

Authors’ contributions

Both authors (LS and FB) contributed to the conception and design of the study. LS perfomed the acquisition and analysis of data. In the process of interpretation and discussion of data both authors were participated. LS drafts the manuscript and it was critically revised by FB. Both authors gave their final approval of the version to be published.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

The person seen in Fig. 1 gave his written consent for publication. Further consent for publication is not applicable.

Ethics approval and consent to participate

The study was conducted according to the declaration of Helsinki. The study was approved by the ethics committee of the University of Potsdam (Germany). All participants gave their informed consent to participate.

Publisher’s Note

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

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Laura V. Schaefer.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Schaefer, L.V., Bittmann, F.N. Are there two forms of isometric muscle action? Results of the experimental study support a distinction between a holding and a pushing isometric muscle function. BMC Sports Sci Med Rehabil 9, 11 (2017). https://doi.org/10.1186/s13102-017-0075-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13102-017-0075-z

Keywords