Abstract
Background
Patients with knee osteoarthritis exhibit greater knee muscle co-contraction, and which in turn is a risk factor for disease progression.
Objective
The aims of this study were to assess the muscle co-contraction during walking and quadriceps exercises in patients with medial knee osteoarthritis, and to identify the factors most associated with the alterations in co-contraction.
Methods
The subjects were 12 women with medial knee osteoarthritis and 13 asymptomatic controls. Kinematic and kinetic parameters during walking, and co-contraction indices during walking, knee extension and straight leg raising exercises were measured. The factors associated with the alteration in co-contraction in the osteoarthritis group were determined by linear regression analyses.
Results
Compared with the control group, the osteoarthritis group showed significantly greater co-contraction of the vastus lateralis and biceps femoris (VLBF) in two quadriceps exercises as well as the VLBF, co-contraction of the vastus lateralis and lateral gastrocnemius (VLLG) in walking. In linear regression analyses, the VLBF and VLLG in walking were associated with knee extension moment and knee pain frequency in daily life. The VLBF in knee extension and straight leg raising exercise were associated with knee flexion angle in a quiet standing position and knee varus peak angle during the stance phase of walking.
Conclusions
Greater muscle co-contractions were observed in not only walking but quadriceps strengthening exercises in patients with knee OA. Besides, the factors of these patterns were different in different tasks. These findings may help the development of more effective rehabilitation strategies for patients with knee osteoarthritis.
Introduction
Muscle co-contraction refers to the simultaneous recruitment of synergistic muscles, agonist and antagonist muscles crossing a joint during dynamic tasks.1,2 Muscle co-contraction of agonist and antagonist muscles, such as the knee flexors and extensors crossing a joint of the lower limbs, is a common strategy adopted to reduce strain and shear forces at the joint, but it increases compressive forces and joint mechanical load. 3
Cumulative mechanical load is a risk factor for knee osteoarthritis (OA), a degenerative joint disease.4,5 Patients with knee OA exhibit greater muscle co-contraction around the knee joint during walking,6,7 which may be an attempt to stabilize the knee joint and avoid pain.6,8 However, this response inhibits the knee’s shock absorption mechanism and increases harmful loads on the joint, which can result in progression of the OA. 9 Therefore, the reduction of excessive co-contraction is important to reduce OA symptoms and slow the progression of OA.
Various factors have been reported to contribute to the greater muscle co-contraction during walking in patients with knee OA, including knee instability,10,11 medial compartment joint loading,12,13 internal knee extension moment, 14 the severity of OA,15–17 knee muscle strength, 18 gait speed, 19 and knee pain.14,18,20,21 Studies have examined the associations between greater co-contraction and knee instability,10,11 medial compartment joint loading, 13 internal knee extension moment, 14 the severity of OA, 15 and knee pain.14,20 However, to the best of our knowledge, no study has extensively evaluated multiple factors, so the most associated factor with the greater co-contraction in patients with knee OA remains unclear.
Quadriceps strengthening exercise programs are commonly prescribed for patients with knee OA to improve muscle strength,22,23 knee pain, physical function, and quality of life.22,24 However, there is no consensus on whether strengthening exercises reduce muscle co-contraction.18,25 The basic muscle strengthening exercises most often used for patients with knee OA are open kinetic chain (OKC) exercises, single-joint or single-muscle group exercise. 26 In healthy subjects, greater muscle co-contraction around the knee joint has rarely been observed during OKC exercises, whereas it has been reported during closed kinetic chain (CKC) exercises such as the squat and lunge.27,28 In patients with knee OA, greater muscle co-contraction has been observed while a CKC motion (descending a step), 6 but there has been no investigation of co-contraction during OKC quadriceps exercises. If muscle co-contraction increases during simple OKC quadriceps exercises, repeating these exercises may exacerbate the co-contraction pattern, with potential progression of the OA.
The aims of this study were (i) to investigate knee muscle co-contraction during walking and OKC quadriceps exercises in patients with knee OA and (ii) to identify the most associated factors with alterations in muscle co-contraction. Understanding the factors associated with alterations in muscle co-contraction in knee OA could lead to the development of more effective rehabilitation strategies.
Methods
Subjects
Characteristics of the knee OA and control groups.
JKOM, Japanese Knee Osteoarthritis Measure; K/L, Kellgren–Lawrence; *, p values of <0.05 were statistically significant. Values are expressed as mean (SD).
a>180° means varus knee.
The severity of knee OA in the OA group was K/L grade I in 5 patients, grade II in 2 patients, grade III in 3 patients, and grade IV in 2 patients. Their main symptom was knee pain. None reported feelings of knee instability.
Assessment of knee pain
The Japanese Knee Osteoarthritis Measure (JKOM) 34 was used to assess the subjects’ knee conditions in their daily lives. Each subject was asked to indicate the intensity and frequency of knee pain in these few days. Pain intensity was valuated using a visual analog scale (VAS) with a range of 0 to 100 mm, 0 indicating no pain and 100 indicating the worst possible pain. Pain frequency was assessed using a 5-point scale (5: always; 4: often; 3: sometimes; 2: rarely; 1: never). 35 Knee pain intensity during the walking trials was also assessed using the VAS.
Knee joint range of motion and knee angles
To measure range of motion (RoM), reflective markers were attached to the greater trochanters and lateral epicondyles of the femurs, the heads of the fibulae, and the lateral malleoli of the tibiae; then photographs were taken of the end angles of passive knee extension and flexion in the supine position, and of the knee flexion angle and knee varus angle in a quiet standing position. 36 The photographs were taken using a digital camera (IXY DIGITAL 510 IS; Canon, Tokyo, Japan). The knee angles were measured from the photographs using ImageJ 1.43u software (National Institutes of Health; MD, USA).
Measurement of maximum joint torque
The maximum isometric knee extension, knee flexion, and ankle plantar flexion torques were measured using a hand-held dynamometer (µTas MT-1; Anima, Tokyo, Japan) during maximal voluntary isometric contraction (MVIC). The knee extension torque was measured in a seated position with the knee flexed at approximately 45°. 19 The knee flexion torque measured with the subject prone, with the knee flexed at approximately 60°. 19 The ankle plantar flexion torque was measured with the subject supine, with the ankle in the neutral position. 21 During the measurement, the subject was asked to grip the edge of the bed with both hands. Isometric contraction for 3 s was performed twice, separated by a 1-min rest period. The torques, averaged over two trials, were normalized to the subject’s body weight (Nm/kg). Electromyogram (EMG) data were recorded during the MVIC trial for each muscle (as described later).
Tasks
The subjects performed three tasks: a walking task and then two OKC quadriceps exercise tasks. Prior to the measurements, the subjects practiced each task to become familiar with the motion and movements. For the walking task, the subjects were instructed to walk along a flat 10-m walkway at a self-selected gait speed. The task was performed three times, and the mean values were used in the analysis. The two OKC quadriceps exercise tasks were a knee extension task and a straight leg raising (SLR) task. For the knee extension task, the subjects sat on a height-adjustable chair without a backrest and armrests, with their lower legs vertical and their arms folded across their chest. After an oral cue, the subjects extended the knee joint from the starting position to its maximum angle over 3 s, held it there for 3 s, and then returned to the starting position over 3 s. For the SLR task, the subjects lay supine on a bed with both legs extended straight and their arms folded across their chest. After an oral cue, the subjects flexed the hip joint of one leg over 3 s until their toe touched a bar that was 10 cm above the toe at the starting position, held their leg there for 3 s, and then returned to the starting position for 3 s. A metronome set at 60 beat/min was used to time the motion speed. The quadriceps exercise tasks were performed twice, and the mean values were used in the analysis.
Motion data acquisition and analysis
Infrared-reflective markers were attached to 45 landmarks: the temples, lateral ends of the superior nuchal line, tragi, acromia, olecranon processes, styloid processes of the ulnae, superior edges of the iliac crests, anterior superior iliac spines, posterior superior iliac spines, greater trochanters, lateral and medial epicondyles of the left and right femurs, lateral and medial joint spaces of the left and right knee joints, lateral and medial condyles of the left and right tibiae, bilateral lateral and medial malleoli, bilateral first and fifth metatarsal heads, and bilateral calcaneal tuberosities. Kinematic data were collecting using a three-dimensional motion analysis system (Vicon MX; Vicon Motion Systems, Oxford, UK) with six infrared cameras at a sampling rate of 100 Hz. At the same time, three-dimensional ground reaction forces were collected using eight force plates (TF-400-A; Tec Gihan, Kyoto, Japan) at a sampling frequency of 1000 Hz. The stance phase during walking (from initial contact to toe-off of ipsilateral leg) was defined by a vertical ground reaction force >10 N. 37 For the quadriceps exercises, joint movement from initiation to termination was defined by the joint angle being greater than the mean plus three standard deviations of the joint angle for the baseline period prior to the oral cue.
The data processing and analysis were performed using BodyBuilder software (Vicon Motion Systems, Oxford, UK). The coordinates of the joint center, joint excursion, peak angle from initial contact to mid-stance of ipsilateral leg during walking, and the maximum, mean, and impulse of the knee joint moment in the stance phase were calculated as described previously38,39 (Supplemental Table). The joint moments were normalized to the subject’s body weight.
EMG acquisition and analysis
Surface EMG data were acquired using an EMG Master (Mediarea Support Business Union, Okayama, Japan), with a pass-band of 15–500 Hz and a sampling frequency of 1000 Hz. The skin surface was prepared by shaving and then cleaning with the skin preprocessing agent Skin Pure (Nihon Kohden, Tokyo, Japan). Bipolar surface circular Ag/AgCl electrodes (Blue-sensor M-00-S; Ambu, Ballerup, Denmark) were placed 30-mm apart, center to center, in line with the muscle fibers over the vastus medialis, vastus lateralis, semitendinosus, biceps femoris, lateral gastrocnemius, and medial gastrocnemius muscles. The electrodes were placed according to the SENIAM guidelines,
40
with a reference electrode placed over the radial styloid process. To analyze the co-contraction of muscle pairs, co-contraction indices (CCIs) were calculated from 100 ms prior to the initial contact to the toe-off of ipsilateral leg during walking, and from 100 ms prior to movement initiation to movement termination during the quadriceps exercises,41,42 for the following muscle pairs: vastus medialis and semitendinosus (VMST), vastus medialis and medial gastrocnemius (VMMG), vastus lateralis and biceps femoris (VLBF), and vastus lateralis and lateral gastrocnemius (VLLG) (Figure 1). The highest root mean square signals during each MVIC trial, averaged over a moving 50-ms window, were used for normalization of the values for each muscle.
43
The raw EMG data for each trial were subjected to full-wave rectification and then expressed as a percentage of the EMG value obtained during MVIC (%MVIC). Each CCI was determined as follows, using the method described by Lewek et al.
10
The definition of analytical range of the muscle co-contractions. The muscle co-contractions were determined using the equation described by Lewek et al.
7
The sample waveforms of vastus lateralis (VL) (dashed), biceps femoris (BF) (dotted), and co-contraction for VL and BF (VLBF) (solid) in a subject with OA. The upper panel (a) displays the waveforms during walking. Dashed vertical lines show the timing of initial contact, 100 ms prior to initial contact, and toe-off. The middle panel (b) displays the waveforms during knee extension exercise. Dashed dotted line shows the knee joint angle in the sagittal plane (full extension is 0° and the positive values indicate the knee joint flexion angle). Dashed vertical lines show the timing of knee motion initiation, 100 ms prior to knee motion initiation, and knee motion termination. The lower panel (c) displays the waveforms during straight leg raising exercise. Dashed dotted line shows the hip joint angle in the sagittal plane (the positive values indicate the hip joint flexion angle). Dashed vertical lines show the timing of hip motion initiation, 100 ms prior to hip motion initiation, and hip motion termination.

Here, the terms lower EMGi and higher EMGi refer to the muscle activity in the less active and more active muscle, respectively, in each muscle pair. 6
Statistical analyses
Statistical analyses were performed using IBM SPSS Statistics 26 (IBM Japan, Tokyo, Japan). The Shapiro–Wilk test was applied to assess variables for normality of their distributions. For normally distributed data, two-sample t-test (for homoscedastic data) or Welch’s test (for heteroscedastic data) were used. The Mann–Whitney U test was used if the data were not normally distributed.
In the OA group, Pearson’s correlation analysis (or Kendall’s correlation analysis if the data were not normally distributed) was used to evaluate the relationships between the CCIs and the following parameters: gait speed, knee pain intensity and frequency, maximum joint torques, RoM of passive knee extension, knee flexion angle in a quiet standing position, K/L grading system, knee varus peak angle during the stance phase of walking, and external knee varus and internal knee extension moments. Stepwise linear regression analyses were used to determine the most associated factors with alterations in CCIs in OA group. The independent variables were the parameters significantly correlated with CCIs in the correlation analyses. To avoid multicollinearity, if any two independent variables showed a correlation with a correlation coefficient was ≥0.70, the more statistically significant of the two variables was entered into the regression model. The significance level for all tests was <0.05.
Results
Subject characteristics
The subject’s characteristics, knee varus angle in a quiet standing position, JKOM total scores, and radiographic scores are shown in Table 1. The knee varus angle was significantly larger in the OA group than in the control group (p < 0.001). There were no significant differences between the two groups in age, height, body weight, or body mass index.
Knee pain, maximum joint torque, RoM, and knee angle
Knee pain, maximum joint torque, knee joint range of motion, and walking parameters.
VAS, visual analog scale; IC, initial contact; Mst, mid-stance; *, p-values of <0.05 were statistically significant. Values are expressed as mean (SD) or median [IQR].
aFull extension angle of the knee joint is 0° and the negative values indicate the limited knee extension.
bFull extension angle of the knee joint is 0° and the positive values indicate the knee joint flexion angle.
cSum of flexion and extension excursion.
The RoMs of passive knee extension and flexion were significantly smaller in the OA group than in the control group (p < 0.01; p < 0.05). The knee flexion angle in a quiet standing position showed significantly greater flexion in the OA group than in the control group (p = 0.001).
There were no significant differences in any of the maximum joint torques between the two groups.
Walking parameters
Table 2 summarizes the walking parameters in the stance phase. The gait speed was slower in the OA group than in the control group (p < 0.01). Knee extension excursion and total excursion (the sum of flexion and extension excursion, Supplemental Table) were significantly less in the OA group than in the control group (p < 0.05). Knee flexion excursion was also less in the OA group, but the difference was not statistically significant (p = 0.054). The knee varus peak angle during the stance phase was greater in the OA group than in the control group, but the difference was not statistically significant (p = 0.079).
The internal knee extension moment impulse was significantly greater in the OA group than in the control group (p = 0.001). The graph of the external knee varus moment in the stance phase showed a bimodal shape 44 ; we therefore compared the two peaks between the groups (Supplemental Table). The second peak, mean, and impulse of the external knee varus moment were significantly greater in the OA group (p < 0.05). No differences were found between the two groups for the maximum and mean internal knee extension moment and for the first peak of the external knee varus moment.
CCIs
CCIs for the stance phase of walking and quadriceps strengthening exercises.
CCIs, Co-contraction indices calculated from surface EMG data using the method described by Lewek et al. 7 ; VMST, vastus medialis and semitendinosus; VMMG, vastus medialis and medial gastrocnemius; VLBF, vastus lateralis and biceps femoris; VLLG, vastus lateralis and lateral gastrocnemius; *, p-values of <0.05 were statistically significant. Values are expressed as mean (SD) or median [IQR].
an = 12 for co-contraction indices for the stance phase of walking and n = 11 for co-contraction indices for quadriceps strengthening exercises (data of quadriceps exercises could not be collected from a subject due to the experimental problems).
bn = 13 for co-contraction indices for the stance phase of walking and n = 11 for co-contraction indices for quadriceps strengthening exercises (data of quadriceps exercises could not be collected from two subjects due to the experimental problems).
Correlation analyses for the OA group
The results of correlation analyses in the knee OA group.
VLBF, vastus lateralis and biceps femoris; VLBF during ext ex., vastus lateralis and biceps femoris during knee extension exercise; VLBF during SLR ex., vastus lateralis and biceps femoris during straight leg raising exercise; VLLG, vastus lateralis and lateral gastrocnemius; K/L, Kellgren–Lawrence; r, correlation coefficient; *, p-values of <0.05 were statistically significant.
Linear regression analyses for the OA group
The results of simple and multiple regression analyses in the knee OA group.
95% CI, 95% confidence interval; VLBF, vastus lateralis and biceps femoris; VLBF during knee extension ex., vastus lateralis and biceps femoris during knee extension exercise; VLBF during SLR ex., vastus lateralis and biceps femoris during straight leg raising exercise; VLLG, vastus lateralis and lateral gastrocnemius; R2, adjusted R2; B, unstandardized regression coefficients; β, standardized regression coefficients; *, p-values of <0.05 were statistically significant; ES, effect sizes with Cohen’s f2 (R2/1-R2). 39
Discussion
This study demonstrated relationships between greater VLBF and a higher maximum internal knee extension moment in walking, greater VLBF in knee extension exercises and a greater knee flexion angle in a quiet standing position, greater VLBF in SLR and a greater knee varus peak angle in stance phase of walking, and greater VLLG in walking and a higher knee pain frequency in daily life. To the best of our knowledge, this is the first study to establish the major factor among many factors associated with greater co-contraction in patients with knee OA.
In this study, most of subjects in the OA group (83%) had K/L grading systems I–III, indicating mild to moderate joint changes.12,21 The OA group showed gait characteristics typical of patients with knee OA, including reduced gait speed, decreased knee excursion in the sagittal plane, increased knee varus peak angle during the stance phase, and increased knee external varus and internal extension moment,46–51 as well as greater muscle co-contraction in the lateral knee muscles. This finding was consistent with previous reports that patients with relatively early stage of knee OA have greater co-contraction in the lateral knee muscles to protect medial joint.16,21 We observed associations between greater both VLLG and VLBF in walking and higher knee pain frequency in daily life. Chronic pain and fear of pain can induce the modulation of muscle activity, including increases in the activity of antagonist muscle.52–54 The association between knee pain and greater activity of the lateral gastrocnemius has also been reported. 20 The greater lateral muscle activities of knee may help to redistribute the medial load laterally. 55 These reports and the findings of this study suggest that greater VLLG and VLBF in walking may be muscle-guarding activity to reduce the medial knee pain that experienced by patients with knee OA in daily life. Patients with knee OA are typically asked about knee pain intensity; these findings suggest it may also be necessary to ask about pain frequency in the clinical interviews with patients with knee OA. VLBF in walking was significantly correlated with not only knee pain frequency in daily life (positive correlation) but also knee pain intensity during the walking trial (negative correlation) and the internal knee extension moments (positive correlation). These results may indicate that patients with OA increase muscle-guarding activity against the knee pain and stress experienced in daily life were attempting to reduce knee pain during the walking trial. However, alterations in antagonist muscle recruitment negatively affect tibiofemoral kinematics56,57 and increase the load on knee joint structures such as the surfaces of the tibiofemoral joint, muscles, and ligaments.12,58–60 The OA group walked with greater flexion in the knee joint in the stance phase than the control group. Greater VLBF and VLLG in walking may have interfered with the normal motion of knee extension. The tibiofemoral joint contact area decreases as the knee flexion angle increases. 61 Thus, increased knee flexion in walking increase the contact pressure on the tibiofemoral joint. The OA group also showed increased external varus and internal extension knee moments, indicating mechanical stress on the knee joint. In summary, our results suggest that the knee joints of the subjects in OA group were subjected to higher load arising from kinematic, kinetic, and muscular factors. Greater co-contraction may help stabilize the knee joint and relieve pain in the short term, but it may lead to abnormal knee motion and harmful loads on the knee joint, potentially resulting in OA progression in the long term.
The subjects with knee OA showed greater VLBF even in the OKC quadriceps exercises. This suggests that they may have become accustomed to controlling their knees with antagonist muscles, even though the quadriceps muscle is the prime mover muscle in these exercises. The linear regression analyses indicated that greater VLBF in knee extension exercise in the OA group was associated with a greater knee flexion angle in a quiet standing position. This malposture could result in alterations in muscle activity or flexibility.62,63 A standing posture with knee full extension needs little activity of the quadriceps and hamstring muscles because the stability of the knee joint depends on the passive stability system provided by, for example, the bones and ligaments.64,65 However, the activity of these muscles increases as the knee flexion angle increases.66,67 Appropriate joint stability and movement requires harmony between the passive (bones and ligaments), active (muscles), and neural systems. 68 A habitual standing posture with knee flexion, indicating a greater dependence on muscle activity in standing, might affect neuromuscular adaptation and result in alterations in muscle activity during knee extension exercises. Furthermore, the association between greater VLBF in SLR and a greater knee varus peak angle during the stance phase of walking may also be an indication of neuromuscular adaptation around the knee joint. Patients with knee OA who exhibit greater knee varus thrust may exert greater muscle co-contraction to keep the knee joint as straight as possible during walking. 13 Thus, habitual abnormal knee varus movement during walking might affect muscle activity in SLR exercise, the task of which keeps their knee straight. These results provide evidence in support of not recommending aggressive quadriceps strengthening exercises for patients with knee OA who have excessive malalignment or laxity of the knee. 69 More effective rehabilitation strategies may require different exercises to be prescribed according to the specific characteristics of patients with knee OA, such as their knee flexion angle in a quiet standing posture or the knee varus peak angle during the stance phase of walking. Our findings may also prompt further studies to develop more effective rehabilitation strategies for patients with knee OA.
This study had several limitations. First, after correlation analyses, additional regression analyses were performed to find a more direct relationship with CCIs. However, the number of subjects was small for regression analyses; further research is needed with a larger number of subjects. However, the results of the regression analyses showed large effect sizes (Cohen’s f2 > 0.35), 45 suggesting the results are meaningful. Nevertheless, further studies are needed with a larger number of subjects. Second, our sample population was limited to patients with medial knee OA who could walk without ambulatory aids. Finally, we did not confirm whether treatment for these factors reduced the co-contraction pattern. Future research is required to address these limitations.
Greater muscle co-contractions were observed in not only walking but also OKC quadriceps strengthening exercises in patients with knee OA. These patterns may be due to cumulative daily factors, including maximum internal knee extension moment and knee varus peak angle during the stance phase of walking, knee pain frequency, and knee flexion angle in a quiet standing position. Besides, the factor of these greater co-contractions may be different in different tasks. These findings could help the development of more effective rehabilitation strategies for patients with knee OA.
Supplemental Material
Supplemental Material - Factors associated with greater co-contraction during walking and quadriceps exercises in patients with knee osteoarthritis
Supplemental Material for Factors associated with greater co-contraction during walking and quadriceps exercises in patients with knee osteoarthritis by Yoshio Wakimoto, Masaya Anan, Kenji Tanimoto, Hiroka Hattori and Koichi Shinkoda in Physiotherapy Practice and Research
Footnotes
Author note
This work was performed at Hiroshima University, Hiroshima, Japan
Acknowledgements
The authors thank Dr. Masahide Yagi of Kyoto University for his valuable advices on data analysis.
Ethical approval
The study was undertaken with the approval of the Ethics Committee of Division of Physical Therapy and Occupational Therapy Sciences, Graduate School of Health Sciences, Hiroshima University (No. 1201).
Informed consent
Prior to the study, the purpose of study was explained to all subjects, who all signed informed consent.
Funding
The authors received no financial support for the research, authorship, and/or publication of this article.
Declaration of conflicting interests
The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Supplemental material
Supplemental material for this article is available online.
Appendix
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
