Inter- and intra-session reliability of muscle activity patterns during cycling

The aim of this study was to determine the inter- and intra-session reliability of the temporal and magnitude components of activity in eight muscles considered important for the leg cycling action. On three separate occasions, 13 male non-cyclists and 11 male cyclists completed 6 min of cycling at 135, 150, and 165 W. Cyclists completed two additional 6-min bouts at 215 and 265 W. Surface electromyography was used to record the electrical activity of tibialis anterior, soleus, gastrocnemius medialis, gastrocnemius lateralis, vastus medialis, vastus lateralis, rectus femoris, and gluteus maximus. There were no differences (P > 0.05) in the muscle activity onset and offset or in the iEMG of any muscles between visits. There were also no differences (P > 0.05) between cyclists and non-cyclists in the variability of these parameters. Overall, standard error of measurement (SEM) and intra-class correlation analyses suggested similar reliability of both inter- and intra-session muscle activity onset and offset. The SEM of activity onset in tibialis anterior and activity offset in soleus, gastrocnemius lateralis and rectus femoris was markedly higher than in the other muscles. Intra-session iEMG was reliable (coefficient of variation (CV) = 5.3–13.5%, across all muscles), though a CV range of 15.8–43.1% identified low inter-session iEMG reliability. During submaximal cycling, the temporal components of muscle activity exhibit similar intra- and inter-session reliability. The magnitude component of muscle activity is reliable on an intra-session basis, but not on an inter-session basis.     

Quantifying thigh muscle co-activation during isometric knee extension contractions: Within- and between-session reliability

Muscle co-activation around the knee is important during ambulation and balance. The wide range of methodological approaches for the quantification of co-activation index (CI) makes comparisons across studies and populations difficult. The present study determined within- and between-session reliability of different methodological approaches for the quantification of the CI of the knee extensor and flexor muscles during maximum voluntary isometric contractions (MVICs). Eight healthy volunteers participated in two repeated testing sessions. A series of knee extension MVICs of the dominant leg with concomitant torque and electromyographic (EMG) recordings were captured. CI was calculated utilizing different analytical approaches. Intraclass correlation coefficient (ICC) showed that within-session measures displayed higher reliability (ICC > 0.861) and lower variability (Coefficient of variation; CV < 21.8%) than between-session measures (ICC < 0.645; CV > 24.2%). A selection of a 500 ms or larger window of RMS EMG activity around the PT delivered more reliable and less variable results than other approaches. Our findings suggest that the CI can provide a reliable measure for comparisons among conditions and is best utilized for within-session experimental designs.     

Algorithmic detection of the beginning and end of bipolar electrograms: Implications for novel methods to assess local activation time during atrial tachycardia

Activation mapping is required to effectively ablate atrial tachycardia (AT). Conventional tools to assess local activation time (LAT) are based upon the peak of the bipolar electrogram (B-EGM, LAT_Peak) and the maximal negative slope of the unipolar electrogram (U-EGM, LAT_Slope). Bipolar electrograms are influenced by wavefront direction, bipole orientation, and inter-electrode spacing causing ambiguity in peak detection, whereas unipolar electrograms are disturbed by the presence of far-field signals. We developed a new algorithm to detect the beginning and end of bipolar electrograms (t_begin and t_end). Then, we introduced new LAT methods related to the onset of B-EGMs (LAT_Onset), the center of mass of B-EGMs (LAT_CoM), and the slope of U-EGMs within a pre-defined window (LAT_Slope-hybrid).In total 3752 recordings from 31 AT patients were retrospectively analyzed. The signal-to-noise ratio (SNR) for B-EGMs was calculated to differentiate algorithmically high from low quality electrograms (HQ and LQ). In a subset of 328 B-EGMs, five experts validated the t_begin as determined by the algorithm by visual rating. The newly developed LAT methods were compared to the conventional LAT methods and to one another (Bland–Altman plots) in both HQ (n = 3003) and LQ EGMs (n = 749).The t_begin algorithm was accurate (deviation < ±10 ms) in 96 ± 4% of HQ and 91 ± 8% of LQ B-EGMs. BA plots revealed the following difference (bias) and variation in HQ and LQ EGMs respectively: (1) LAT_Onset vs. LAT_Peak: 27 ± 30 ms and 24 ± 62 ms; (2) LAT_CoM vs. LAT_Peak: 0 ± 16 ms and 2 ± 38 ms; (3) LAT_Slope-hybrid vs. LAT_Slope: 1 ± 32 ms and 15 ± 110 ms; (4) LAT_Onset vs. LAT_CoM: 22 ± 24 ms and 18 ± 22 ms; (5) LAT_Onset vs. LAT_Slope-hybrid: 16 ± 18 ms and 13 ± 22 ms; and (6) LAT_CoM vs. LAT_Slope-hybrid: 5 ± 20 ms and 4 ± 18 ms.In the present study, we introduced three new methods to assess local activation time in AT, based upon an algorithm detecting accurately the beginning and end of the B-EGM complex. BA analysis of the new methods showed similar variation in high and low quality EGMs, suggesting that they introduce less ambiguity than the conventional peak method. LAT_Onset consistently yielded an earlier activation moment. LAT_Slope-hybrid – by blanking far-field potentials – seems to be the optimal method for detection of the maximal negative slope in U-EGMs. Interestingly, LAT_CoM in B-EGMs coincided with the maximal negative slope in U-EGMs, suggesting its physiological sense and future use. The new LAT methods can be implemented in real-time mapping applications.     

Long-term stability of tensiomyography measured under different muscle conditions

Tensiomyography (TMG) is a technique utilised to measure mechanical and contractile properties of skeletal muscle. Aim of this study was to assess long-term stability of TMG across a variety of muscle conditions. Gastrocnemius Medialis of 21 healthy males was measured using TMG in rested conditions, after a warm-up, after a maximal voluntary contraction and after a fatigue protocol. Participants were re-tested on a second occasion 4 weeks apart. Among the parameters examined, Contraction Time, Sustain Time and Delay time exhibited a good level of absolute reliability (CV = 3.8–9.4%) and poor to excellent level of relative reliability (ICC = 0.56–0.92). On the other hand, relative reliability was good to excellent for muscle Displacement (ICC = 0.86–0.96), whereas its level of absolute reliability was questionable (CV = 8.0–14.8%). Minimum detectable change was less than 20% in most conditions for the aforementioned parameters. Half-relaxation Time yielded overall insufficient reliability. In general, the level of reliability tended to increase after the maximal voluntary contraction and the fatigue protocol were administered, probably because of more controlled conditions preceding the measurement. Information about the long-term stability of TMG across different muscle conditions is essential when intervention studies are undertaken with an exercising population, particularly athletes.     

Does gantry rotation time influence accuracy of volume computed tomography dose index (CTDIvol) in modern CT?

PurposeThe purpose of this study was to develop a gantry overrun corrected CTDI_vol (cCTDI_vol) dosimetry and evaluate the differences between the displayed CTDI_vol (dCTDI_vol), measured CTDI_vol (mCTDI_vol), and the cCTDI_vol.Methods and materialsThe each 8 rotation times between 275 and 1000 ms of two CT scanners were investigated. Rotation time (T_rot) and the beam-on time (T_beam) in axial scanning were measured accurately to determine the gantry overrun time (T_over) as T_beam − T_rot. Subsequently, mCTDI_vol was measured by using a 100 mm ionization chamber and CTDI phantoms. Furthermore, we introduced a gantry overrun correction factor (C_o = T_rot/T_beam) to obtain cCTDI_vol. Upon completion of the data acquisition, the dCTDI_vol and mCTDI_vol were compared with the cCTDI_vol.ResultsThe discrepancies of T_rot were 0.2 ± 0.2 ms as compared to the preset rotation times, and T_over was machine-specific and almost constant (22.4 ± 0.5 ms or 45.1 ± 0.3 ms) irrespective of the preset rotation time. Both dCTDI_vol and mCTDI_vol were increasingly overestimated compared to cCTDI_vol as the faster the preset rotation time was selected (1.7–23.5%).ConclusionThe rotation time influences the accuracy of CTDI_vol in modern CT, and should be taken into consideration when assessing the radiation output in modern CT.     

Intrarater reliability and agreement of linear encoder derived heel-rise endurance test outcome measures in healthy adults

A linear encoder measuring vertical displacement during the heel-rise endurance test (HRET) enables the assessment of work and maximum height in addition to the traditional repetitions measure. We aimed to compare the test-retest reliability and agreement of these three outcome measures. Thirty-eight healthy participants (20 females, 18 males) performed the HRET on two occasions separated by a minimum of seven days. Reliability was assessed by the intraclass correlation coefficient (ICC) and agreement by a range of measures including the standard error of measurement (SEM), coefficient of variation (CV), and 95% limits of agreement (LoA). Reliability for repetitions (ICC = 0.77 (0.66, 0.85)) was equivalent to work (ICC = 0.84 (95% CI 0.76, 0.89)) and maximum height (ICC = 0.85 (0.77, 0.90)). Agreement for repetitions (SEM = 6.7 (5.8, 7.9); CV = 13.9% (11.9, 16.8%); LoA = −1.9 ± 37.2%) was equivalent to work (SEM = 419 J (361, 499 J); CV = 13.1% (11.2, 15.8%); LoA = 0.1 ± 34.8%) with maximum height superior (SEM = 0.8 cm (0.6, 1.0 cm); CV = 6.6% (5.7, 7.9%); LoA = 1.3 ± 17.1%). Work and maximum height demonstrated acceptable reliability and agreement that was at least equivalent to the traditional repetitions measure.     

High-field MR spectroscopy in the multiparametric MRI evaluation of breast lesions

IntroductionThe aim of this work was to assess the role of 3T-MR spectroscopy (MRS) in the multi-parametric MRI evaluation of breast lesions, using a pattern-recognition based classification method.Methods291 patients (301 lesions, median 2.3 cm³) were enrolled in the study (age 18–85 y, mean 54.2 y). T1-TSE (TR/TE = 400/10 ms) and T2-STIR imaging (TR/TE = 5000/60 ms), dynamic-contrast-enhanced MRI (DCE-MRI), apparent diffusion coefficient (ADC) (b = 0–800 s/mm²), and single-voxel MRS (10 × 10 × 10 mm³, PRESS, TR/TE = 3000 ms/135 ms) were performed by means of a 3T scanner. MRS results were accepted if the FWHM of the water peak was ⩽45 Hz. Total choline (tCho) was considered detected if the signal-to-noise ratio (SNR) of the 3.2 ppm peak was ⩾2. A classifier-based analysis (support-vector-machines, SVM) was performed with 4-dimensional vectors including type of margin, DCE-MRI kinetic curve type, ADC mean value, and tCho SNR. A comparison with 3-dimensional vectors (without tCho SNR) was used to assess MRS impact on sensitivity, specificity, and positive-negative predictive values (PPV-NPV) for malignancy.Results228 lesions (180 malignant/48 benign) showed acceptable spectral quality. Comparison of classification results with histopathological examination of surgical specimens showed sensitivity = 93.7%, specificity = 84.9%, PPV = 95.2%, NPV = 81.5% without the inclusion of MRS in the SVM analysis. When MRS was included, the figures increased to 95.1%, 90.7%, 97.2%, and 85.0%, respectively.ConclusionsInclusion of 3T-MRS in the multi-parametric MRI evaluation of breast lesions improved the performance of the SVM-based classifier, showing a possible role of high-field MR spectroscopy in the differential diagnosis between benign and malignant breast lesions. Further research is however needed to confirm this initial evidence.     

The architecture and contraction time of intrinsic foot muscles

Although critical for effective human locomotion and posture, little data exists regarding the segmentation, architecture and contraction time of the human intrinsic foot muscles. To address this issue, the Abductor Hallucis (AH), Abductor Digiti Minimi (ADM), Flexor Digitorum Brevis (FDB) and Extensor Digitorum Brevis (EDB) were investigated utilizing a cadaveric dissection and a non-invasive whole muscle mechanomyographic (wMMG) technique. The segmental structure and architecture of formaldehyde-fixed foot specimens were determined in nine cadavers aged 60–80 years. The wMMG technique was used to determine the contraction time (Tc) of individual muscle segments, within each intrinsic foot muscle, in 12 volunteers of both genders aged between 19 and 24 years.While the pattern of segmentation and segmental –architecture (e.g. fibre length) and –Tc of individual muscle segments within the same muscle were similar, they varied between muscles. Also, the average whole muscle Tc of FDB was significantly (p < 0.05) shorter (faster) (Tc = 58 ms) than in all other foot muscles investigated (ADM Tc = 72 ms, EDB Tc = 72 ms and ABH Tc = 69 ms). The results suggest that the architecture and contraction time of the FDB reflect its unique direct contribution, through toe flexion, to postural stability and the rapid development of ground reaction forces during forceful activities such as running and jumping.     

Signal-to-noise ratio of diffusion weighted magnetic resonance imaging: Estimation methods and in vivo application to spinal cord

Diffusion tensor imaging (DTI) and tractographic reconstruction may be applied for in vivo clinical spinal cord studies. However, this structure represents a challenge to current acquisition and reconstruction strategies, due to its small size, motion artifacts, partial volume effects and low signal-to-noise-ratio (SNR). Aims of this work were to select the best approach for the estimate of SNR and to use it for spinal cord diffusion weighted (DW) sequence optimization.Seven methods for the estimate of SNR were compared on uniform phantom DW images, and the best performing approach (single ROI for signal and noise, difference of images—SNR_diff) was applied for the following in vivo sequence evaluations.Fifteen sequences with different parameters (voxel size, repetition (TR) and echo (TE) times) were compared according to SNR, resolution, fractional anisotropy (FA) and tractography performances on three healthy volunteers. In vivo optimization of DW sequences resulted in: axial sequence, with voxel size = 1.5 mm × 1.5 mm × 3.5 mm, TR = 3200 ms and TE = 89 ms, sagittal sequence with voxel size = 2.2 mm × 2.2 mm × 2 mm, TR = 3000 ms and TE = 84 ms.An objective method tested on phantom and a practical index for in vivo spinal cord DTI SNR estimation allowed to obtain axial and sagittal optimized sequences, providing excellent tractographic results, with acceptable acquisition times for in vivo clinical applications.     

Kinematic and electromyographic analysis of the Nordic Hamstring Exercise

The Nordic Hamstring Exercise (NHE) has been introduced as a training tool to improve the efficiency of eccentric hamstring muscle contraction. The aim of this study was to perform a biomechanical analysis of the NHE. Eighteen participants (20.4 ± 1.9 years) performed two sets of five repetitions each of the NHE and maximal eccentric voluntary contraction (MEVC) of the knee flexors on an isokinetic dynamometer whilst knee angular displacement and electrical activity (EMG) of biceps femoris were measured. EMG was on average higher during the NHE (134.3% of the MEVC). During the forward fall of the NHE, the angle at which a sharp increase in downward velocity occurred varied between 47.9 and 80.5 deg, while the peak knee angular velocity (pVelocity) varied between 47.7 and 132.8 deg s^?1. A significant negative correlation was found between pVelocity and peak EMG (r = ?0.62, p < 0.01) and EMG at 45 deg (r = ?0.75, p < 0.01) expressed as a percentage of peak MEVC EMG. Some of the variables analyzed exhibited good to excellent levels of intra- and inter-session reliability. This type of analysis could be used to indirectly monitor the level of eccentric strength of the hamstring muscles while performing the NHE and potentially any training- or injury-related changes.     

Electromyography responses of pediatric and young adult volunteers in low-speed frontal impacts

No electromyography (EMG) responses data exist of children exposed to dynamic impacts similar to automotive crashes, thereby, limiting active musculature representation in computational occupant biomechanics models. This study measured the surface EMG responses of three neck, one torso and one lower extremity muscles during low-speed frontal impact sled tests (average maximum acceleration: 3.8 g; rise time: 58.2 ms) performed on seated, restrained pediatric (n = 11, 8–14 years) and young adult (n = 9, 18–30 years) male subjects. The timing and magnitude of the EMG responses were compared between the two age groups. Two normalization techniques were separately implemented and evaluated: maximum voluntary EMG (MVE) and neck cross-sectional area (CSA). The MVE-normalized EMG data indicated a positive correlation with age in the rectus femoris for EMG latency; there was no correlation with age for peak EMG amplitudes for the evaluated muscles. The cervical paraspinous exhibited shorter latencies compared with the other muscles (2–143 ms). Overall, the erector spinae and rectus femoris peak amplitudes were relatively small. Neck CSA-normalized peak EMG amplitudes negatively correlated with age for the cervical paraspinous and sternocleidomastoid. These data can be useful to incorporate active musculature in computational models, though it may not need to be age-specific in low-speed loading environments.     

Tensiomyography of selected lower-limb muscles in professional soccer players

Tensiomyography is a non-invasive method of neuromuscular assessment used to measure muscle action characteristics, muscle tone, and muscle fiber type, and provides information on acute and chronic responses of muscle to different training loads. The aims of the present study were: to analyse differences in muscle response and mechanical characteristics of two major muscles of the lower extremity in a large group of Spanish soccer players according to playing position, and to provide group norms against which clinical findings may be compared. Data were collected from 78 professional soccer players (age 26.6 ± 4.4 years; height: 179.2 ± 5.3 cm; body mass: 75.8 ± 5.3 kg). Tensiomyography was recorded from the rectus femoris (RF) and biceps femoris (BF) muscles after 2 days without take part in any strenuous exercise or training. Five tensiomyographic parameters were analyzed: maximal displacement (D_m), contraction time (T_c), sustain time (T_s), delay time (T_d), and half-relaxation time (T_r). A good to excellent intra-session reliability was found for all contractile parameters (ICC ranged from 0.78 to 0.95). No significant differences between players of any position were observed in absolute values of BF. However, significant differences were observed for T_c, T_r and T_s between the different playing positions on RF (P < 0.05, effect size ranged from 1.3 to 1.6). Professional soccer players showed muscles with ability to rapidly generate force during contractions. The neuromuscular profile provided could help in identifying the normative data that are important for the different positions in order to optimize the training and recovery process of each individual player.     

Biceps femoris and semitendinosus tendon/aponeurosis strain during passive and active (isometric) conditions

The purpose of this study was to quantify strain and elongation of the long head of the biceps femoris (BFlh) and the semitendinosus (ST) tendon/aponeurosis. Forty participants performed passive knee extension trials from 90° of knee flexion to full extension (0°) followed by ramp isometric contractions of the knee flexors at 0°, 45° and 90° of knee flexion. Two ultrasound probes were used to visualize the displacement of BFlh and ST tendon/aponeurosis. Three-way analysis of variance designs indicated that: (a) Tendon/aponeurosis (passive) elongation and strain were higher for the BFlh than the ST as the knee was passively extended (p < 0.05), (b) contraction at each angular position was accompanied by a smaller BFlh tendon/aponeurosis (active) strain and elongation than the ST at higher levels of effort (p < 0.05) and (c) combined (passive and active) strain was significantly higher for the BFlh than ST during ramp contraction at 0° but the opposite was observed for the 45° and 90° flexion angle tests (p < 0.05). Passive elongation of tendon/aponeurosis has an important effect on the tendon/aponeurosis behavior of the hamstrings and may contribute to a different loading of muscle fibers and tendinous tissue between BFlh and ST.     

Acute changes in knee cartilage transverse relaxation time after running and bicycling

PurposeTo compare the acute effect of running and bicycling of an equivalent cumulative load on knee cartilage composition and morphometry in healthy young men. A secondary analysis investigated the relationship between activity history and the change in cartilage composition after activity.MethodsIn fifteen men (25.8±4.2 years), the vertical ground reaction force was measured to determine the cumulative load exposure of a 15-min run. The vertical pedal reaction force was recorded during bicycling to define the bicycling duration of an equivalent cumulative load. On separate visits that were spaced on average 17 days apart, participants completed these running and bicycling bouts. Mean cartilage transverse relaxation times (T₂) were determined for cartilage on the tibia and weight-bearing femur before and after each exercise. T₂ was measured using a multi-echo spin-echo sequence and 3T MRI. Cartilage of the weight bearing femur and tibia was segmented using a highly-automated segmentation algorithm. Activity history was captured using the International Physical Activity Questionnaire.ResultsThe response of T₂ to bicycling and running was different (p=0.019; mean T₂: pre-running=34.27 ms, pre-bicycling=32.93 ms, post-running=31.82 ms, post-bicycling=32.36 ms). While bicycling produced no change (−1.7%, p=0.300), running shortened T₂ (−7.1%, p<0.001). Greater activity history predicted smaller changes in tibial, but not femoral, T₂.ConclusionsChanges in knee cartilage vary based on activity type, independent of total load exposure, in healthy young men. Smaller changes in T₂ were observed after bicycling relative to running. Activity history was inversely related to tibial T₂, suggesting cartilage conditioning.     

Muscle activation patterns of knee flexors and extensors during passive and active movement of the spastic lower limb in chronic stroke patients

The aim of this study was to describe the characteristics of spasticity, quantified as muscle activity during stretch, during passive and active movement. For this cross sectional study 19 stroke patients with spasticity in the lower limb were recruited. Reflex activity was studied with surface electromyography of knee flexor and extensor muscles during passive and active movement of the lower leg.On both the affected and unaffected side, root mean square values of the knee extensor muscles, while stretched, were higher during active than during passive movement (p < 0.05). For the vastus lateralis (VL) the correlation was moderate (ρ = 0.54, p = 0.022), for the rectus femoris (RF) high (ρ = 0.83, p < 0.001). For the semitendinosus (ST) the correlation was low (ρ = 0.27) and not significant.During active movement the correlation between VL activity and activity of the antagonist ST, as an indicator for co-contraction of the affected muscles, was marked (ρ = 0.73, p = 0.001). A moderate negative correlation was found between reflex activity of RF during passive stretch and the active range of motion (ρ = ?0.51, p = 0.027).The results show that a passive stretch test alone is insufficient either as assessment method for spasticity during active motor tasks or as a measure for motor control.     

Surface electromyography reveals males have a slower patellar reflex than females

In humans the cross sectional area of spinal motor neurons at L3 is larger in males than in females. Since these contribute to the control of the quadriceps femoris muscle group and are involved in the patellar reflex (PR), gender differences in the PR are expected. We have investigated this possibility using a group of 28 young subjects (14 male and 14 female) aged 20–22 years. The PR was quantified by the muscle compound action potential (MCAP) from the surface electromyogram (sEMG) of the vastus lateralis muscle. We found that the PR latency in females (17 ± 0.19 ms), was significantly (p < 0.001) faster than in males (21 ± 0.37 ms). This 4 ms difference in latency could not be ascribed to differences in stature or thigh length. In conclusion, for the age range tested females posses a significantly faster patellar reflex than males. We suggest that the slower PR latency of male subjects may arise in part from their larger α-motorneurons: such that longer integration times are required for the summation of postsynaptic excitation to be sufficient to excite α-motorneurons.     

Single-electron photoreduction of the PM intermediate of cytochrome c oxidase

The P_M → F transition of the catalytic cycle of cytochrome c oxidase from bovine heart was investigated using single-electron photoreduction and monitoring the subsequent events using spectroscopic and electometric techniques. The P_M state of the oxidase was generated by exposing the oxidized enzyme to CO plus O₂. Photoreduction results in rapid electron transfer from heme a to oxoferryl heme a₃ with a time constant of about 0.3 ms, as indicated by transients at 605 nm and 580 nm. This rate is ∼ 5-fold more rapid than the rate of electron transfer from heme a to heme a₃ in the F → O transition, but is significantly slower than formation of the F state from the P_R intermediate in the reaction of the fully reduced enzyme with O₂ to form state F (70–90 μs). The ∼ 0.3 ms P_M → F transition is coincident with a rapid photonic phase of transmembrane voltage generation, but a significant part of the voltage associated with the P_M → F transition is generated much later, with a time constant of 1.3 ms. In addition, the P_M → F transition of the R. sphaeroides oxidase was also measured and also was shown to have two phases of electrogenic proton transfer, with τ values of 0.18 and 0.85 ms.     

The highest antagonistic coactivation of the vastus intermedius muscle among quadriceps femoris muscles during isometric knee flexion

Although the possibility that the vastus intermedius (VI) muscle contributes to flexion of the knee joint has been suggested previously, the detail of its functional role in knee flexion is not well understood. The purpose of this study was to examine the antagonist coactivation of VI during isometric knee flexion. Thirteen men performed 25–100% of maximal voluntary contraction (MVC) at 90°, 120°, and 150° knee joint angles. Surface electromyography (EMG) of the four individual muscles in the quadriceps femoris (QF) was recorded and normalized by the EMG signals during isometric knee extension at MVC. Cross-talk on VI EMG signal was assessed based on the median frequency response to selective cooling of hamstring muscles. Normalized EMG of the VI was significantly higher than that of the other synergistic QF muscles at each knee joint angle (all P < 0.05) with minimum cross-talk from the hamstrings to VI. There were significant correlations between the EMG signal of the hamstrings and VI (r = 0.55–0.85, P < 0.001). These results suggest that VI acts as a primary antagonistic muscle of QF during knee flexion, and that VI is presumably a main contributor to knee joint stabilization.