Abstract
BACKGROUND:
A method of measurement of voluntary activation (VA, percent of full muscle recruitment) during isometric and isokinetic concentric contractions of the quadriceps femoris (QF) at 60
OBJECTIVE:
This study aimed to quantify the test-retest minimal real difference (MRD) of VA during isometric (ISOM) and isokinetic concentric contractions of QF (100
METHODS:
VA was measured through the interpolated twitch technique. Pairs of electrical stimuli were delivered to the QF at 40
RESULTS:
The VA (median, range) was 84.20% (38.2–99.9%) in ISOM and 94.22% (33.8–100%) in ISOK. MRD was 0.78 and 1.12 logit for ISOM and ISOK, respectively. As an example, in terms of percent VA these values correspond to a change from 76% to 95% and from 79% to 98% in ISOM and in ISOK, respectively.
CONCLUSIONS:
The provided MRD values allow to detect significant individual changes in VA, as expected after training and rehabilitation programs.
Introduction
Voluntary activation (VA) is defined as the degree of recruitment of muscles engaged in a voluntary contraction effort [1]. VA is usually measured as a percentage of the maximum force that could be provided by a muscle if all its fibres contract at tetanic frequency. A 100% VA is rarely achieved for the quadriceps femoris (QF), and for most muscles, even among trained healthy participants [2]. The most commonly used method for VA measurement is the interpolated twitch technique (ITT) which consists of stimulating a representative portion of the muscle belly through an electric shock, causing some initial discomfort but no pain. The stimulus is delivered both at rest and during maximal contraction effort. If the shock does not generate extra force during contraction, all muscle fibres belonging to the sample reached by the electric shock can be claimed to be tetanized [1, 3]. Otherwise, the contraction/rest ratio of the peak forces during the “twitches” evoked by the shock provides an estimate of the percentage of muscle fibres not (or not fully) engaged by the volitional drive: VA is the complement to this ratio. Therefore, the ITT provides an estimation of how voluntary effort translates into muscle force production [3].
A less common procedure, known as the central activation ratio (CAR), implies that the peak moment is evoked with supramaximal stimulation of the femoral nerve, making the resting twitch unnecessary [4, 5, 6, 7, 8]. The two methods provide quite superimposable results. In previous research on various muscle groups VA has been measured through ITT during isometric contractions (ISOM) [2, 9, 10, 11] and during both muscle shortening (“concentric contraction”) [9, 12] and lengthening (“eccentric contraction”) [13, 14, 15, 16]. Given that muscle force depends both on muscle length and velocity of shortening (or lengthening), most studies adopted isokinetic dynamometers Here, the activation of the QF during lower limb extension is of concern. The extant literature agrees (despite initial controversies [16, 17, 18]) that VA of QF is incomplete during ISOM (around 80–95%) [10, 19, 20, 21]. In contrast, the debate is ongoing for isokinetic (ISOK) concentric and eccentric contractions. Some authors concluded that maximal ISOK concentric VA did not differ significantly from ISOM VA either at slow (20
In a previous paper [1] a method of measurement of VA during ISOM and ISOK knee extension (velocities of 60 and 120
A key property of measurement is its test-retest reliability. This parameter for VA has been evaluated in isometric conditions at 90
Another neglected issue was the nonlinearity of percentage measures, flawed by floor and/or ceiling effects. In the present study, measures of VA and their changes were analysed both as raw percentages (VA) and as values transformed (VA
The primary endpoint of the present study was to determine the MRD of VA
Methods
The trial was conducted from September 2019 to April 2020 in a research laboratory on human movement, within a hospital department of neurorehabilitation.
Participants
Twenty-five healthy participants (18–45 years old, 12 women) took part in the experiments. The inclusion criteria were: i) ability to sign the informed consent form; ii) ability to understand the instructions and to complete the motor task; and iii) to be recreationally active.
Exclusion criteria were: i) pregnancy; ii) a history of epilepsy (to avoid the risk of seizures triggered by the stimuli); iii) implanted electro-sensitive devices; iv) any neurologic or orthopaedic condition limiting the articular mobility or muscular strength of the lower limbs; and v) current treatment with oral anticoagulant or antiplatelet therapy (to avoid the risk of muscle haemorrhage). None of the participants were familiar with the testing method.
Ethics
All participants gave written informed consent before participating. The study was approved by the ethics committee of the institution (project code 24C721_2017 date of approval 14
Procedures
Test sequencing
Participants took part in the experiment during two sessions, test (T0) and retest (T1), with a 14-day interval. The two sessions took place at about the same time of the day and participants were asked to maintain the same lifestyle habits between the sessions. The 14-day interval was chosen as it is close to the time frame usually adopted for re-testing after rehabilitation or sports training. Both lower limbs were tested. The order in which the lower limbs were tested was randomised (see below).
Trial sequencing
In a separate way for both men and women, for the first enrolled participant the limb to be tested first (dominant vs. non-dominant) was randomly selected. For subsequent participants, the sequence in limb testing was alternated. Then, those participants who had their dominant limb tested as first at T0 started retest with the non-dominant limb, and vice-versa.
Anthropometric measurements
All anthropometric parameters were measured at T0 and at T1 at the beginning of the testing session. Participants were tested for foot dominance using the Waterloo footedness questionnaire-revised [27]. Height and weight were measured.
Instruments
All tests were performed using a Cybex Humac Norm
Range of knee motion and rotation speeds
Isometric tests were conducted with the knee held at 40
Electrical stimulation of quadriceps muscle
VA levels were determined according to the ITT method. Electric shocks were delivered through large percutaneous rubber electrodes (120
Schematic representation of the experimental protocol for one limb. The experimental session starts with isometric (ISOM) habituation in which 4 submaximal repetitions were requested. Then, one ISOM repetition at 100%MVC was asked for, in order to calculate the parameters of the electrical doublets to be delivered in the test (represented by the vertical gray bar). After 3-min rest, the participants performed two ISOM repetitions separated by a 3-min rest. In each ISOM repetition, three electrical doublets (represented by arrows) were delivered: the first one was delivered at rest, the second one was delivered 3–4 s after the beginning of the effort initiation, and the third one was delivered at rest 2–3 s after the contraction ended. After the ISOM test, a 5-min rest separated the habituation trials, made by two isokinetic knee extension-flexion contractions at 100
Joint moments were displayed online on a PC-screen. The participant was instructed to neglect the screen and to focus on perceiving her/his maximal effort. Verbal encouragement was provided to the participants to elicit maximal effort.
To achieve habituation at the beginning of each session participants performed 2–4 submaximal isometric contractions at increasing effort at the pre-determined joint angle (40
After a 3-minute pause, ISOM contractions were requested. During ISOM (knee immobilised at 40
A 5-minute pause separated the end of ISOM testing from the beginning of isokinetic habituation: two extension-flexion repetitions at maximal effort at the chosen angular velocity (100
Surface electromyography (sEMG) recording
The sEMG was recorded at rest to ascertain the fully relaxed state of the vastus lateralis and the homolateral biceps femoris muscle. Two pairs of silver chloride surface electrodes were applied to each muscle according to the SENIAM guidelines [35]. Low impedance (
2.3.9.1. Data acquisition and analysis
A signal acquisition system (CED 1401, Cambridge Electronic Designed Limited-Ced, Cambridge) was used to simultaneously record knee moment, angle, and angular velocity through the isokinetic dynamometer. These variables were sampled at 500 Hz. A customised software (Spike 2, version 8, Cambridge Electronic Design Limited-CED, Cambridge) connected to the acquisition system controlled the electric stimulator in eliciting the stimulus at the desired knee angle. The average moment applied to the transducer when the lower limb was fully extended was used to offset the gravitational moment at all angles by proper computation. Moment signals were off-line filtered (low-pass FIR filter, 1.3 Hz).
The VA was determined by the following formula:
Tracings from a representative participant (woman, 35 years, 53 kg, 163 m). The gravity-corrected moment (Nm), on the ordinate, is given as a function of time (s), on the abscissa. From top to bottom, in the left column the panels refer to a maximal isometric contraction at 40
where IT and RT are the peak moments caused by the electric doublets during the voluntary contraction and at rest, respectively. During ISOM contractions, pre-contraction RT was adopted in the computation of VA. Indeed, like the twitch during contraction the post-contraction RT may exhibit some potentiation, but it may also be diminished by fatigue [36, 37]. The pre-contraction twitch is preferred for consistency with ISOK computation of VA [1]. During ISOK contractions, the time course of the moments was visually analysed off-line to detect any gross abnormalities. In particular, IT had to occur during the descending part of the bell-shaped moment curve (see ref. 1 for details) (Fig. 2).
2.3.9.2. Algebraic considerations
The goal of this study was the computation of the minimal real difference (MRD) of VA, at two time points. The MRD is the minimal change that cannot be attributed to random error, only (usually, at
where
2.3.10. Statistical analyses
For each participant, the best repetition out of two ISOM contractions and out of three ISOK contractions were selected for statistical analysis [39, 40, 41]. The normality of data distribution was assessed using the Shapiro-Wilk’s test. When applicable, means (SD) were computed and Student’s
Means, medians, first lower quartiles (Q1) and third upper quartiles (Q3) for age, weight, height, body mass index (BMI) of the dominant and the non-dominant limbs of the twenty-five participants
The ordinate gives the voluntary activation (VA) values of the participants (
The left ordinate gives the voluntary activation levels after logit transformation (VA
The absolute agreement between VA
2.3.11. Software
Computations, statistics, and graphic representations were performed using MATLAB
Means, medians, 5
From left to right, the columns report the values of intraclass correlation coefficients (ICC
Twenty-five subjects (12 women), median age 27.3 (24.0
In agreement with the classic muscle force-velocity relationship [45], higher moment peaks were recorded during ISOM with respect to ISOK, both at T0 and T1. Means, medians, and 5
Tests on raw VA data are given in the Appendix, Table A1.
In the present healthy sample, complete activation was never achieved during ISOM, either on the dominant or on the nondominant limb by any participant. During ISOK at T0 one woman and three men reached VA 100% with the dominant limb, while four women and four men reached VA 100% with the nondominant limb. During ISOK at T1 two women and three men reached VA 100% with the dominant limb, while nine women and four men reached VA 100% with the nondominant limb. Among these participants, one woman and one man, only, performed 100% VA during all isokinetic tests.
The distribution of VA
Bland-Altman plot of the differences between the VA
It can be seen that the logit transformation provided a wider spread of measures (hence, a higher discrimination), compared to the original VA observations. The distribution of the data was tested with Shapiro Wilk’s test (Appendix: Table A2). The repeated ANOVA model on VA
In contrast, no differences were recorded depending on side dominance (
The values of ICC
Figure 5 shows a Bland-Altman plot on VA
The results indicate that the proposed method of isokinetic measurement of VA on a sample of healthy participants through the ITT method may be reliable and clinically acceptable. To the best of the authors’ knowledge, this is the first study proposing aggregate data and test-retest reliability evaluation of quadriceps muscle VA under isokinetic conditions.
Amount and reliability of activation
These results confirm that VA of the quadriceps muscle is incomplete during maximal ISOM contractions [2, 4, 19, 20, 23]. In contrast, some participants in this study displayed complete activation (i.e. VA 100%) during maximal ISOK at 40
Comparison with previous studies; reasons for discrepancies
In this sample VA levels did not differ significantly between ISOK and ISOM, although mean and median values were slightly higher in isokinetic conditions (see Table 2). Previous findings by Babault reported slightly higher VA levels in isokinetic concentric contractions at 120
A possible explanation for these discrepancies lies in the higher VA levels reported for QF at joint angles closer to 90
In this sample, as far as between-limb differences are concerned, the dominant and non-dominant limbs displayed superimposable VA
The present study revealed statistically significant differences between women and men, with women demonstrating higher VA
Of relevance here, in the above cited studies and in most currently used isokinetic VA measurement techniques, the electrical stimulus was triggered when the joint reaches a certain angle during its movement. This technique of stimulation introduces an element of variability in the time of elicitation of the interpolated twitch along the curve because its latency depends on many factors: type of contraction (isometric vs. isokinetic), joint angular velocity, muscle conduction velocity, electrode positioning, time for excitation-contraction coupling, and tendon stiffness [1]. Therefore, comparisons are difficult when different angular velocities and knee joint angles are involved. In contrast, the stimulation technique performed in the present study is characterised by the computation of the instant of stimulation in order to have the peak of the interpolated twitch to occur consistently at exactly the desired knee angle [1].
Advantages of logit linearisation
In the present study, only the analysis performed on transformed data was able to detect a sex-related difference in VA
a difference in muscle fibre structure. Although the overall distribution of fast and slow fibre types in the vastus lateralis muscle is similar between women and men, a significant difference exists when the total area occupied by each fibre type is considered. In a study by Staron et al. [52], slow fibres were observed to occupy a greater area in the vastus lateralis muscle of women, whereas fast IIA fibres occupied a greater area in the vastus lateralis muscle of men. According to Henneman’s size principle of motor unit recruitment, the fast motor units also have the highest recruitment thresholds [53]. As already hypothesised by Behm et al. [2], fast-twitch predominant muscles may be more difficult to fully activate during voluntary contractions. Therefore, when compared to women, the quadriceps femoris muscle of men might display lower mean and median values in voluntary activation levels. The results may reflect a (hardly avoidable) mechanical artefact. Since VA is calculated as the complement to the ratio of IT and RT, the modification of one of these two parameters would alter the resulting VA value. In particular, an anterior tilt of the pelvis during knee extension would imply a shortening of the rectus femoris muscle and, as a consequence, a loss of knee extensor moment. The subsequent reduction of IT would lead to an increase in VA. Anatomical resting pelvic tilt does not display significant betweensex differences [54, 55]. Nevertheless, women are known to have significantly lower abdominal muscle thickness than men [56], hence a potentially lower capacity for preventing pelvic tilt during maximal knee-extension effort. In addition, women might demonstrate later onsets in anticipatory postural adjustments of trunk muscles compared to men in single lower limb lift movements [57]. Again, these differences may eventually lead to a pelvic tilt causing lower IT peaks and higher VA estimates in women.
The suggestion of using logit units, unfamiliar to most clinicians, may appear a major limitation of the method. However, the simple computation of the MRD (second Equation in the Methods) can be easily performed on any spreadsheet or hand calculator. This value represents a threshold that is easy to interpret for decisions on individuals. Its critical advantage on percentage scores is that the measure of change becomes independent of the baseline values. In other words, a 1-logit change means the same amount of change of “activation” at whatever level of the 0–100% scale. It is not so for raw percent changes, which are “compressed” (and therefore, underestimated) when the VA values are approaching the extremes of the acceptable range [46]. The linearization problem, often overlooked in the biomedical literature, arises with VA but not with moments: these are intrinsically linear (i.e., a 1-moment change means the same change whatever was the baseline level). Moments, however, were not the topic of the present article.
The provided MRD values found here can be used to exemplify the logit-to-percentages transformation of change values. The MRD was 0.78 and 1.12 logits for ISOM and ISOK, respectively. A 0.78 logit change might reflect a change in VA from 28% to 70% (see the second Equation in the Methods) and from 76% to 95%. Likewise, a 1.12 logit change might reflect a change in VA from 79% to 98% and from 30% to 85%. In clinical practice, mirror reasoning may apply. The typical question may be checking whether a given observed percentage change does reveal or not a “real” change. In the latter example above, baseline and retest ISOK VA values of 30% and 85% correspond to
On the size of the observed MRD
In clinical terms, the MRD values of activation, for both ISOM and ISOK contexts, seem rather large. Stated otherwise, a large change must be observed before it can be considered as “real”. This may reflect a) limitations of the present study as well as b) difficulties intrinsic to the measurement of this variable.
In this particular study, some test-retest instability (thus increasing the MRD) may be due to the small sample size and to the presence of outliers. On the other hand, there were no previously available data allowing for a-priori estimation of sample size. In addition, even smaller sample sizes are common in the literature on isokinetic quadriceps dynamometry [10, 12, 46]. Some instability might also stem from the absence of substantial habituation-learning processes, prevented by the experimental design. On the other hand, the design aimed to reproduce clinical contexts, in which no substantial practice is usually allowed. With respect to activation itself, this seems to be intrinsically unstable. In fact, VA is a psychological (not a biological) variable which may be influenced by numerous extraneous factors (learning, fatigue, fear, motivation, attention, etc.). In addition, in human physiology the maximum voluntary muscle activation is close to its limit (i.e., 100%). The “crowding” of values hinders the precision of measurement, hence reliability studies.
Last, a “real” change, i.e., a statistically significant one, is not necessarily a clinically important one. In general, clinical “significance” requires statistical significance, but the cut-off is higher and based on external criteria (e.g., mortality risk, the incidence of a given disease, success probability in an examination, or whatever) [58]. In the case of VA
The study limits cannot be overemphasised. The results refer to isometric contraction and to one isokinetic velocity (100
Author contributions
CONCEPTION: Luigi Tesio.
PERFORMANCE OF WORK: Luigi Catino, Chiara Malloggi, Valeria Cerina and Stefano Scarano.
INTERPRETATION OR ANALYSIS OF DATA: Luigi Catino, Chiara Malloggi, Stefano Scarano and Luigi Tesio.
PREPARATION OF MANUSCRIPT: Luigi Catino and Stefano Scarano.
REVISION FOR IMPORTANT INTELLECTUAL CONTENT: Luigi Catino Chiara Malloggi, Valeria Cerina, Stefano Scarano, Viviana Rota and Luigi Tesio.
SUPERVISION: Luigi Tesio.
Ethical considerations
All participants gave written informed consent before participating. The study was approved by the ethics committee of the institution (project code 24C721_ 2017, date of approval 14
Funding
The study was supported by Istituto Auxologico Italiano, IRCCS [AMULET, “ricerca corrente” project code 24C721_2017].
Footnotes
Acknowledgments
The authors have no acknowledgments.
Conflict of interest
The authors have no conflicts of interest to report.
Appendix
Results from A) Shapiro-Wilk’s, B) and C) Wilcoxon signed rank tests on percent voluntary activation (VA) data for the isometric contractions (ISOM) and the isokinetic concentric contractions (ISOK) for both limb sides (dominant and non-dominant) during the test (T0) and the retest (T1).
VA
Shapiro-Wilk’s test
Bonferroni corrected
ISOM_T0
0.000
ISOM_T1
0.005
ISOK_T0
0.000
ISOK_T1
0.000
ISOM_dominant
0.000
ISOM_non-dominant
0.000
ISOK_dominant
0.000
ISOK_non-dominant
0.000
Results from the Shapiro-Wilk test on logit-transformed voluntary activation (percent voluntary activation, VA, converted to the logit of the raw values, VA
VA
Shapiro-Wilk test
Bonferroni corrected
ISOM_T0
0.30
ISOM_T1
0.00
ISOK_T0
0.24
ISOK_T1
0.05
ISOM_dominant
0.54
ISOM_non-dominant
0.00
ISOK_dominant
0.11
ISOK_non-dominant
0.10
Results from the repeated ANOVA model for VA
VA
Repeated ANOVA model
R
0.06
P
Model
0.08
Sex
0.01
Test (ISOM vs. ISOK)
0.18
Lower limb side dominance
0.29
Time point
0.27
Dominance#Time point
0.49
Sex#Time point
0.89
