ORIGINAL RESEARCH article

Front. Physiol., 27 May 2025

Sec. Striated Muscle Physiology

Volume 16 - 2025 | https://doi.org/10.3389/fphys.2025.1593851

Longitudinal and transverse muscle stiffness change differently with knee osteoarthritis and do not align with stiffness sensation

Angela V. DieterichAngela V. Dieterich1Katrin SkerlKatrin Skerl2Filip PaskaliFilip Paskali3Leonardo Gizzi,Leonardo Gizzi4,5Mehrin AzanMehrin Azan2Gabriela F. CarvalhoGabriela F. Carvalho1Matthias KohlMatthias Kohl3Andreas Haueise,
Andreas Haueise1,6*
  • 1Ultrasound Lab, Department of Physiotherapy, Faculty of Health, Medical, and Life Sciences, Furtwangen University, Freiburg, Germany
  • 2Institute of Technical Medicine, Faculty of Health, Medical, and Life Sciences, Furtwangen University, Furtwangen, Germany
  • 3Institute of Precision Medicine, Faculty of Health, Medical, and Life Sciences, Furtwangen University, Villingen-Schwenningen, Germany
  • 4Fraunhofer Institute for Production Engineering and Automation, Department of Biomechatronics, Stuttgart, Germany
  • 5Institute for Modelling and Simulation of Biomechanical Systems IMSB, University of Stuttgart, Stuttgart, Germany
  • 6Department of Sport and Sport Science, University of Freiburg, Freiburg, Germany

Knee osteoarthritis (OA) is a prevalent joint condition associated with disability, pain, and stiffness in the muscles surrounding the knee. Myofascial and massage techniques are employed to alleviate these symptoms. Unclear is whether muscle stiffness is measurably increased in the painful muscles, and how measured muscle stiffness relates to perceived stiffness, pain, and functional impairment. Given muscle anisotropy, stiffness changes may occur in the longitudinal plane parallel to muscle fibers or perpendicularly in the transverse plane. Shear wave velocity (SWV) was measured in both scanning planes in 21 individuals with diagnosed knee OA, 21 sex- and age-matched controls, and 20 young controls, focusing on the gastrocnemius medialis and vastus lateralis muscles under four conditions: (a) resting state, (b) holding the shank against gravity, (c) double-leg stance, and (d) single-leg stance. Median stiffness measures, the ratio of longitudinal-to-transverse stiffness, and the ratio of single-leg stance-to-baseline stiffness were compared between groups using the Kruskal- Wallis and Pairwise Asymptotic Wilcoxon rank sum tests. Correlations with the Knee Injury and Osteoarthritis Outcome Score and the Tampa Scale of Kinesiophobia were examined. Longitudinal stiffness of the gastrocnemius medialis muscle was significantly lower in the OA group at double-leg (P = 0.033) and single-leg stance (P = 0.019), with tendencies toward lower median stiffness in both muscles across most tasks. Transverse stiffness of the vastus lateralis muscle was significantly higher in the OA group at baseline (P = 0.027), with tendencies toward higher median stiffness in both muscles across most tasks. Significant moderate to excellent correlations support the clinical relevance of both longitudinal and transverse stiffness measurements. Measured and perceived stiffness were not correlated. Study results suggest that knee OA may differentially affect muscle stiffness in the longitudinal and transverse planes and that transverse stiffness measures may have clinical relevance.

1 Introduction

Osteoarthritis (OA) of the knee is the most prevalent form of osteoarthritis (Long et al., 2022) characterized by progressive limitations on ambulatory function and joint mobility, significantly impacting the individual’s quality of life (Driban et al., 2020). The perception of joint stiffness is one of the major clinical symptoms of osteoarthritis (Bastick et al., 2016). Patients with knee OA commonly report pain and stiffness in the muscles surrounding the knee joint (Teo et al., 2019; NICE, 2023). Physiotherapy is the first-line treatment recommended for alleviating symptoms and functional impairment (Gibbs et al., 2023). In response to patients’ complaints and based on the assumption of objectively increased muscle stiffness, interventions such as massage (Bahns and Kopkow, 2023) or myofascial techniques (Dor and Kalichman, 2017) are often employed aiming to reduce muscle stiffness. However, it remains unclear whether muscle stiffness is objectively increased in the affected muscles, and how measured muscle stiffness relates to perceived stiffness, pain, and functional impairment.

Ultrasound shear wave elastography is the most frequently used modality to non-invasively estimate tissue stiffness, with muscle shear wave velocity (SWV) or the related shear modulus serving as the primary measures of stiffness (Davis et al., 2019; Wang et al., 2019). To date, two studies explored muscle stiffness using shear wave elastography around the knee joint in patients with knee OA compared to healthy individuals (Gökşen et al., 2021; Li et al., 2021). Gökşen et al. found no significant difference in SVV of the rectus femoris and the biceps femoris muscles between individuals with and without knee OA (Gökşen et al., 2021) while Li et al. reported an increased shear modulus in all three hamstring muscles (Li et al., 2021). In both studies blinding was not reported and manually selected small regions of muscle tissue were measured, introducing a potential risk of bias and raising concerns about representative measurements (Haueise et al., 2023). Additionally, it remains unclear whether the measured muscles were perceived as stiff by the patients. Thus, current literature does not answer the question whether subjective and objective muscle stiffness correspond. Moreover, although patients with musculoskeletal conditions frequently report sensations of tense or stiff muscles, evidence for objectively increased muscle stiffness is conflicting (Haueise et al., 2023). All studies included in the aforementioned review measured muscle stiffness only along muscle fiber direction, which is the standard procedure (Săftoiu et al., 2019). Previous research suggests that large samples (Kolding et al., 2018; Hvedstrup et al., 2020; Li et al., 2023) and eventually the summation of stiffness measures across different muscle parts (Li et al., 2023) may be necessary to detect measurably increased muscle stiffness in pain conditions. This contrasts with the prevailing clinical perception of palpably increased muscle stiffness. The question arises of whether the well-established longitudinal shear wave elastography measurements capture the phenomenon underlying perceived muscle stiffness sufficiently.

Muscle tissue is anisotropic, meaning its stiffness differs between the longitudinal (parallel to the direction of muscle fibers) and the transverse (cross-sectional) planes (Green et al., 2024; Gennisson et al., 2010; Ngo et al., 2024). Longitudinal muscle stiffness reflects both active and passive tension, increasing with muscle activation and stretch (Hug et al., 2015). Transverse muscle stiffness is thought to represent the mechanical properties in cross-sectional direction. Previous studies have indicated that muscle stiffness is lower in the transverse plane and is less influenced by muscle activation compared to longitudinal measurements (Gennisson et al., 2010; Eby et al., 2013; Ngo et al., 2024). Muscle palpation and massage techniques are typically applied transversely to the direction of the muscle fibers (Kolster, 2016). This traditional approach may reflect implicit knowledge and experience of altered transverse muscle properties. To our knowledge, transverse muscle stiffness has not been measured in individuals with knee OA.

The mechanical properties of muscles change with age (Kragstrup et al., 2011; Lieber and Ward, 2013). To distinguish changes of muscle stiffness due to ageing from those associated with pain in OA, two control groups were included, an age- and sex-matched group to isolate OA-specific effects, and a younger sex-matched group to capture age-related changes. The primary aim of the current study was to compare SWV of the most painful muscles around the knee joint in both longitudinal and transverse planes across three groups 1) individuals with diagnosed knee OA, 2) age- and sex-matched control participants, and 3) sex-matched young participants under both resting and active conditions. A previous study questions a relationship between perceived and measured muscle stiffness (Akagi and Kusama, 2015). To examine this and the clinical relevance of the stiffness measures, the secondary aim was to investigate the relationships between measured muscle stiffness and clinical symptoms including pain, stiffness, knee joint function, and kinesiophobia.

2 Methods

This cross-sectional observational study was planned and conducted according to the STROBE-checklist for cross-sectional studies (von Elm et al., 2007). Recruitment and enrolment of participants were in accordance with the Declaration of Helsinki. Ethical approval was obtained from the Ethics Committee of the University of Stuttgart (No 21–012) and each participant provided informed consent prior to participation.

2.1 Participants

Participants aged 45–68 years, either diagnosed with knee osteoarthritis or free of lower limb conditions and discomfort, were recruited through university posters, local newspaper advertisements, and a digital neighborhood platform. Eligibility was assessed via telephone screening based on the following inclusion and exclusion criteria. The group of participants with knee OA (labelled OA) included individuals with a medical diagnosis of knee OA irrespective the OA stage, and painful muscles around the knee joint for at least 6 months. Both healthy control groups comprised individuals without lower limb diagnoses or pain lasting longer than a week within the past 2 years, and were labelled old control and young control to emphasize the age contrast, not to judge age. The age-matched control group included healthy individuals matched for age and sex, whereas the young control group included sex-matched healthy individuals aged 18–30 years. Exclusion criteria for all three groups included systemic diseases affecting muscle performance (e.g. rheumatic diseases, COPD), a history of knee or hip trauma (e.g. ACL rupture), intake of medication potentially affecting muscle tone (e.g. muscle relaxants), and skin lesions at the lower limb. An a priori power analysis was conducted using G*Power (Version 3.1.9.7), based on a cross-sectional study including four groups (Alfuraih et al., 2020). The effect size (f = 0.48) was calculated in G*Power based on the vastus lateralis SWV between the healthy and active rheumatoid arthritis groups, as these yielded the largest differences. As G*Power does not support the Kruskal–Wallis test, a one-way ANOVA was used as its parametric equivalent. With α = 0.05 and a power (1 – β) = 0.80, the required sample size was n = 45 for group comparisons and n = 38 for pairwise Wilcoxon tests. The participants of the OA group were recruited first to allow for matching of the control groups. Demographic information and self-reported measures were collected from all participants using the German versions of the “Knee Injury and Osteoarthritis Outcome Score” (KOOS) (Kessler et al., 2003) and the “Tampa Scale of Kinesiophobia” (TSK) questionnaires (Rusu et al., 2014). The KOOS questionnaire includes subscores for pain, functional limitations, and symptoms, including two items (S6, S7) specifically assessing the intensity of perceived stiffness.

2.2 Procedure

Muscle stiffness was estimated using an Aixplorer Ultimate system (V.12, Supersonic Imagine, Aix-en-Provence, France) equipped with a linear transducer (SL15-4) with a 5 cm footprint. The system was configured for shear wave elastography, MSK preset with 90% gain, depth 2.0–2.5 cm, focus zone centered, sampling frequency 1.2–2.0 Hz, dynamic range 72 dB, a maximal measurement range of <10 m/s, low persistence, a standard setting between resolution and penetration, and 50% opacity.

After completing demographic information and questionnaires, participants were prepared for data collection. Individuals with knee OA were asked to point out and rank their two most painful muscles. The most painful muscle was then scanned using B-mode ultrasound to determine and mark the direction of the muscle fibers, a longitudinal and a transverse scanning location. The locations were chosen on top of the muscle belly, approximately in accordance with the SENIAM recommendations (Hermens et al., 2000). Adhesive foam tape was used to mark and frame the transducer locations to ensure consistency for repeated measurements (Dieterich et al., 2020). To standardize the potential influence of previous activity (Proske and Gandevia, 2012), the prepared participants walked for 5 min at comfortable pace on a treadmill. Data collection started with participants positioned supine or prone on a plinth, depending on their most painful muscle. The ultrasound probe was hosted in a custom-cut foam block, covered with sufficient ultrasound gel, and secured to the leg with hook-and-loop straps (Figure 1). The elastogram was configured to cover the maximum visualizable area of the target muscle, positioned at a depth of 1.2–1.5 cm just below the superficial fascia, to optimize area-based visualization of muscle tissue. SWV of the most painful muscle was recorded at four different conditions, 1) baseline, i.e., resting state, 2) holding the shank against gravity (lift shank), 3) double-leg stance and 4) single-leg stance (Table 1). Measurements were taken in longitudinal and in transverse transducer orientation. Trials were recorded as 30-s video clips except single-leg stance (20 s) to prevent exacerbation of pain. Dependent on recording duration and depth, each clip included 24–54 elastograms (images). All conditions except baseline were repeated three times. After recording SWV in the most painful muscle, measurements were taken from the second most painful muscle only during the more challenging tasks double-leg and single-leg stance. The same muscles considered as the first and second most painful muscles by a participant of the OA group were measured in a sex- and age matched participant from the old control group and a sex-matched participant from the young control group. Data were stored without any indication of group allocation.

Figure 1
www.frontiersin.org

Figure 1. The ultrasound transducer was hosted in a custom-cut foam block, here for recording the M. vastus lateralis.

Table 1
www.frontiersin.org

Table 1. Description of the recording conditions. All conditions were recorded in longitudinal and in transverse transducer orientation. Depending on the most painful muscle, participants were positioned supine or prone on the plinth.

2.3 Data processing

Elastography clips were exported as DICOM files and computer-analyzed using a custom script (MATLAB R2021b, The MathWorks, Natick, MA, United States) modified from Skerl et al. (Skerl et al., 2017) to extract stiffness measures from the complete elastograms within the clips. The SWE clips comprised greyscale B-mode ultrasound data at 11 Hz and stiffness data (elastography) at 1.2–2.0 Hz. Each clip contained 200–750 greyscale images with 24–54 superimposed elastograms. The elastogram’s position within the greyscale images was automatically detected and the elastogram’s upper 1 mm was discarded. The elastogram was segmented into 14 small (4 mm × 4 mm) squares (Figure 2) from which stiffness was extracted separately. For each square in each frame of a clip, the following data were extracted: maximum, mean and median SWV, standard deviation and interquartile range, percentage of colored pixels (representing stiffness information), and the number and percentage of pixels with a SWV >9.9 m/s (suggesting saturated measurements). As a measure of imaging quality control, the number of pixels containing colored elastography values was calculated for each square in each frame. Frames containing less than 50% colored pixels were excluded from further analyses. If more than 50% of the squares were excluded, the elastogram was discarded. Extracted data were exported for each clip separately for the statistical analyses. Group allocation was not disclosed prior to the statistical analysis.

Figure 2
www.frontiersin.org

Figure 2. GUI (graphical user interface) of the custom-programmed script for SWE data processing, example from the young control group, M. gastrocnemius medialis, single-leg stance, longitudinal transducer orientation.

2.4 Statistical analysis

Demographic data, the sum of the KOOS S6 and S7 items asking for stiffness sensations, and the results of the KOOS and TSK questionnaires were described per group by mean and standard deviation. Statistical analysis of the stiffness data was performed using R (Team, 2023) with the following packages: tidyr (Wickham et al., 2023b) and dplyr (Wickham et al., 2023a) for data handling, MKinfer (Kohl, 2023) and SPSS (V.29, IBM corporation) for statistical analysis, and ggplot2 (Wickham, 2016) and patchwork (Pedersen, 2023) for data visualization.

The extracted data included median stiffness values from squares (1–7) in two rows (A; B; Figure 2) over 20- or 30-s recording time. Measurements of SWV from the first and second most painful muscles were combined for analysis. Longitudinal and transverse measurements were evaluated separately for each muscle. Given that muscle function involves three-dimensional changes in muscle shape (Azizi et al., 2017), mechanical alteration in one plane likely relates to the mechanical properties in other planes. Consequently, the ratio of longitudinal to transverse SWV was calculated as a two-dimensional indicator of longitudinal muscle tension relative to transverse elasticity. Additionally, the ratio of SWV at baseline to that during single-leg stance was computed as an indicator of stiffness change with a high activation level. Stiffness measures and stiffness ratios were analyzed by group and muscle.

Reliability of the stiffness measures between repeated recordings was estimated using the ICC(3,3) in SPSS (IBM, Version 29, reliability analysis two-way mixed absolute agreement) and the standard error of measurement (SEM). The Shapiro-Wilk test demonstrated skewed data distributions for most testing conditions. Therefore, muscle stiffness and stiffness ratios were described by median and interquartile range (IQR). Group differences were examined using the Kruskal–Wallis Rank Sum test. Significance was set to P < 0.05. For significant group differences, Pairwise Asymptotic Wilcoxon Rank Sum Tests were used to identify the different groups. Within the knee OA group, the Spearman correlation coefficient ρ (rho) was used to assess the relationships between stiffness measures, stiffness sensations, the KOOS domains, the TSK scores, age, and BMI. Correlation coefficients ρ > 0.75 were judged good to excellent, ρ between 0.50 and 0.75 moderate to good, and ρ between 0.25 and 0.50 fair (Portney and Watkins, 2000).

3 Results

Sixty-two participants took part in the study (Table 2), resulting in 1240 trials of the most painful muscle and 740 trials of the second most painful muscle, which was only recorded during standing activities. In the OA group, the most painful muscle was identified as the vastus lateralis in 11 participants, the gastrocnemius medialis in 7 participants, the vastus medialis in 2 participants, and biceps femoris in 1 participant. For the second most painful muscle, the vastus lateralis was selected by 4 participants, the gastrocnemius medialis by 11 participants, the vastus medialis by 3 participants, and the biceps femoris by 3 participants. Due to the limited number of measurements of the vastus medialis and biceps femoris muscles, the presented results focus on the vastus lateralis and gastrocnemius medialis muscles. Data inclusion following the quality control described in 2.3 was on average 85% for longitudinal scanning and 74% for transverse scanning, varying by muscle. The highest data inclusion was observed for the gastrocnemius medialis muscle in the longitudinal plane (98%), while the lowest was for the vastus lateralis muscle in the transverse plane (64%). The reliability of stiffness measurements was higher in longitudinal than in transverse measurements (Table 3).

Table 2
www.frontiersin.org

Table 2. Demographic description and patient reported outcome measures (KOOS and Tampa scale of kinesiophobia) of the three groups; mean ± SD. Perceived stiffness = KOOS items S6 + S7.

Table 3
www.frontiersin.org

Table 3. Reliability (ICC(3,3) of median stiffness measures in the recording conditions that were repeated three-times, mean (95% confidence intervals) over squares.

3.1 Longitudinal muscle stiffness

Longitudinal measurements in both muscles indicated tendencies toward lower median SWV in the knee OA group (Figure 3; Table 4), with significance in the Kruskal Wallis Rank Sum Test observed for the gastrocnemius medialis muscle at double-leg stance, P = 0.033, and at single-leg stance, P = 0.019, groupwise comparisons see Table 4. At baseline and lift shank, the young control group demonstrated the lowest and the old control group the highest longitudinal stiffness. In the knee OA group, median longitudinal stiffness at baseline and lift shank was similar to the young control group.

Figure 3
www.frontiersin.org

Figure 3. Longitudinal shear wave velocity of the gastrocnemius medialis and vastus lateralis muscles per group at resting state and tasks of increasing activity; logarithmic scaling of the x-axis to enhance the visibility in the lower range. *P < 0.05. The high gastrocnemius medialis stiffness at stance is likely due to the stretch over two joints. Abbreviations: M., muscle; n, number of participants per group; OA, osteoarthritis.

Table 4
www.frontiersin.org

Table 4. Longitudinal shear wave velocity of the gastrocnemius medialis and vastus lateralis muscles per group at resting state and tasks of increasing activity, P values for group differences.

3.2 Transverse muscle stiffness

Transverse measurements demonstrated tendencies toward higher median SWV in the knee OA group acoss most tasks (Figure 4; Table 5), with significance observed in the Kruskal Wallis Rank Sum Test for the vastus lateralis muscle, P = 0.026. In the knee OA group, transverse SWV showed minimal or no increase with rising activity levels. At baseline and lift shank, the OA group demonstrated the highest and the young control group the lowest transverse stiffness.

Figure 4
www.frontiersin.org

Figure 4. Transverse shear wave velocity of the gastrocnemius medialis and vastus lateralis muscles per group at resting state and tasks of increasing activity. *P < 0.05. Abbreviations: M., muscle; n, number of participants per group; OA, osteoarthritis.

Table 5
www.frontiersin.org

Table 5. Transverse shear wave velocity of the gastrocnemius medialis and vastus lateralis muscles per task and group, P values for group differences.

3.3 Ratios of longitudinal-to-transverse stiffness and single-leg stance-to-baseline stiffness

In the knee OA group, the median ratio of longitudinal-to-transverse SWV was significantly lower for the gastrocnemius medialis muscle at single leg stance, P = 0.049, with tendencies toward a lower ratio in the OA group compared to the control groups in all but one measured condition (Table 6). Although not statistically significant (P = 0.051 for transverse SWV of the vastus lateralis muscle), the low single-leg stance-to-baseline ratio of transverse SWV suggests the least change in transverse muscle stiffness with activity in the knee OA group (Figure 5).

Table 6
www.frontiersin.org

Table 6. Ratio of longitudinal relative to transverse SWV of the gastrocnemius medialis and vastus lateralis muscles per group at resting state and tasks of increasing activity.

Figure 5
www.frontiersin.org

Figure 5. Ratio of shear wave velocity at single-leg stance relative to baseline of the longitudinal (left) and the transverse (right) shear wave velocity of the gastrocnemius medialis and vastus lateralis muscles per group; (left) in logarithmic scaling. Note the generally higher ratio for the gastrocnemius muscle due to higher longitudinal stiffness at stance.

3.4 Correlation of stiffness measures with patient reported outcomes

Longitudinal SWV demonstrated significant moderate to excellent correlations with the KOOS symptoms score (double-leg stance, ρ = 0.487 P = 0.047), the KOOS QoL score (double-leg stance, ρ = 0.724 P = 0.005), and with the Tampa Scale of Kinesiophobia (baseline, ρ = 0.668 P = 0.025). Transverse SWV was well to excellently correlated with the KOOS ADL score (baseline ρ = 0.812 P = 0.049; double-leg stance ρ = 0.679 P = 0.022), with the KOOS QoL score (double-leg stance ρ = 0.668 P = 0.025) and with the Tampa Scale of Kinesiophobia (baseline ρ = 0.856 P < 0.001; lift shank ρ = 0.748 P = 0.020) The consistently different signs of the correlations of longitudinal SWV with the KOOS pain score suggest separate consideration of the low activation conditions baseline and lift shank versus the standing tasks (Table 7). Inconsistencies across tasks and muscles request cautious interpretation of the correlation results.

Table 7
www.frontiersin.org

Table 7. Spearman correlation coefficients of longitudinal and transverse shear wave velocity of the gastrocnemius medialis and vastus lateralis muscles in the knee OA group at resting state and tasks of increasing activity.

Considering all participants of the OA group irrespective the most painful muscles, perceived stiffness correlated significantly with age, ρ = 0.491 P < 0.001, BMI ρ = 0.310 P = 0.046, and all KOOS scores: Pain, ρ = −0.526 P < 0.001, ADL ρ = −0.469 P = 0.002, Sport ρ = −0.454 P = 0.003, and QoL ρ = −0.490 P < 0.001; but not with measured muscle stiffness.

4 Discussion

The current study results suggest lower longitudinal stiffness in the painful medial gastrocnemius and vastus lateralis muscles of individuals with knee OA during standing positions compared to age- and sex-matched controls. Conversely, results suggest higher muscle stiffness in the transverse plane, with significant differences in the vastus lateralis muscle at baseline. Across all groups, transverse muscle stiffness increased only marginally with activity, with the smallest increase observed in the knee OA group. Significant moderate to excellent correlations suggest meaningful relationships between both longitudinal and transverse shear wave velocity measurements, the KOOS subscores, and kinesiophobia, but not with the sensation of stiffness.

We identified no comparable studies specifically examining muscle stiffness in individuals with knee OA. Gökşen et al. (2021) used shear wave elastography to investigate stiffness in the rectus femoris and biceps femoris muscles in individuals with early knee OA compared to healthy individuals. They observed highly variable stiffness in both groups and muscles, with a trend toward higher median stiffness in the rectus femoris and lower median stiffness in the biceps femoris muscles in participants with knee OA. Thirty-one participants with early-stage knee OA were measured in a prone position with 30°–45° knee flexion. Several small measurements zones were manually positioned within the elastogram by an examiner; blinding was not reported (Gökşen et al., 2021). In contrast, Li et al. (2021) found significantly elevated stiffness in the hamstring muscles of 50 study participants with knee OA and symptoms of pain, stiffness, and heightened sensitivity in the posterior knee region. Apparently, stiffness measures were recorded in a prone position with an extended knee (2021). The measurement zone was larger than that used in the study by Gökşen et al., but similarly, blinding of the examiner was not reported. In both studies, the manual and maybe unblinded positioning of small measurement zones within muscles of heterogeneous stiffness may have introduced bias (Haueise et al., 2023).

A further study examined 79 participants aged 60–80 years without knee pain or known arthritis and found that increased passive longitudinal stiffness of the quadriceps femoris muscles (summed measures of the vastus lateralis, vastus medialis and rectus femoris muscles) was associated with a higher risk of developing knee OA within 12 months (Li et al., 2023). Measurements were recorded in a supine position with 60° of knee flexion using a computed extraction of stiffness from a large muscle region. The authors suggest that increased quadriceps muscle stiffness may contribute to a higher risk of developing knee OA (Li et al., 2023). We acknowledge the methodological strength of this study, which suggests increased longitudinal muscle stiffness precedes the onset of knee OA. The examined sample and the study question differ from the here presented work.

The aforementioned knee OA studies report inconsistent results and are not comparable to the present study in terms of design or measurement methods. Participants were not asked to identify their most painful muscles nor were functionally relevant weight-bearing conditions included. In the present study, we recorded elastography clips comprising 24–54 elastograms, over which stiffness measures over the full width of the transducer’s footprint were averaged, resulting in highly reliable and representative measurements. The extraction of muscle stiffness has been performed fully computed and thus blinded for group allocation. The current results align with previous work that found equal or lower longitudinal stiffness of muscles with chronic pain (Dieterich et al., 2020; Friede et al., 2020; Klauser et al., 2022; Wang et al., 2022). Longitudinal muscle stiffness relates to the contractile tension built up during muscle activation (Nordez and Hug, 2010; Hug et al., 2015). Thus, lower active longitudinal stiffness in painful muscles suggests reduced muscle activation around the painful arthritic knee joint.

In the past years, transverse muscle stiffness was increasingly measured. Studies focussed on methodological and reliability assessments in healthy participants (Gennisson et al., 2010; Eby et al., 2013; Miyamoto et al., 2015; Chino et al., 2017; Ruby et al., 2019; Vuorenmaa et al., 2022; Ngo et al., 2024). Additionally, a few diagnostic studies have employed transverse stiffness measurements (Nicholls et al., 2020; Kolb et al., 2021a; Kolb et al., 2021b). Most of these studies measured muscle stiffness in a relaxed state. Some studies included assessments under increasing active load in human volunteers (Gennisson et al., 2010) or passive load in specimen (Eby et al., 2013; Ngo et al., 2024). The existing literature confirms lower (Gennisson et al., 2010; Eby et al., 2013; Vuorenmaa et al., 2022; Ngo et al., 2024) or equal (Kolb et al., 2021b) transverse compared to longitudinal stiffness measures in healthy individuals and only marginal increases of transverse stiffness under loading. In systemic sclerosis (Kolb et al., 2021b) and in myositis (Kolb et al., 2021a), increased passive transverse muscle stiffness has been identified as biomarker for disease. Correspondingly, the current study suggests higher transverse stiffness in the knee OA group, with a significant difference observed in the vastus lateralis muscle at baseline. A clinical relevance of transverse stiffness measurements is underlined by the correlation results.

Longitudinal and transverse shear wave elastography measure the velocity of shear waves induced by focused acoustic beams. Consequently, shear wave elastography relies on shear wave propagation along echogenic structures that are organized approximately parallel to the transducer’s footprint. The transverse mechanical properties of muscles are primarily influenced by the connective tissue between muscle fibers - the extracellular matrix (Purslow, 2020; Wheatley, 2020). The extracellular matrix forms a complex network of proteins and fibers that envelop muscle fibers, anchoring them to facilitate efficient force transmission during contraction and movement (Ramaswamy et al., 2011; Marcucci et al., 2019; Sinha et al., 2020; Minato et al., 2021; Wohlgemuth et al., 2023; Spadoni et al., 2024). Additionally, it adapts to the mechanical conditions of the muscle (Jozsa et al., 1988; Järvinen et al., 2002; Kjaer et al., 2006; Hyldahl et al., 2015; Gihring et al., 2020). Thus, transverse muscle stiffness may be interpreted as the elasticity of the connective tissue between muscle fibers.

Compared to longitudinal measurements, transverse measures of muscle stiffness are technically challenging due to an increased difficulty in tracking shear waves (Eby et al., 2013). This is indicated by areas of void pixels where stiffness information is absent. Since shear wave tracking depends on echogenic structures along which shear waves propagate, the orientation of the intramuscular connective tissue is of high interest. Purslow described the endomysial network as a “planar feltwork of quasi-randomly orientated, wavy fibers” lying in a plane parallel to the muscle fiber surface (2020). Similarly, the perimysium, which ensheathes several muscle fiber bundles, has been described as a lattice of wavy collagen fiber bundles arranged at +55° and −55° angles relative to the muscle fiber direction (Purslow, 2010). A famous cross-sectional electron micrograph of the endomysial and perimysial networks (Purslow and Trotter, 1994; Purslow, 2020; Figure 1) suggests structures which probably enable shear wave propagation approximately perpendicular to the muscle fibers. Variations in the orientation of the intramuscular connective tissue may increase the difficulty of shear wave transmission and tracking, potentially resulting in lower quality elastograms. However, Miyamoto et al. (2015), Eby et al. (2013) and Gennisson et al. (2010) reported negligible stiffness differences in measurements angled up to 45° and 90° relative to muscle fiber direction suggesting that transverse stiffness measures are robust even when the measurement direction is not perfectly aligned with the echogenic structures.

In the present study, the low single-leg stance-to-baseline ratio of transverse stiffness in the knee OA group suggests the least stiffness change with high-level contraction compared to the control groups. Also, a trend toward higher transverse stiffness in the knee OA group was observed, both pointing to reduced elasticity of the between-fiber connections. During contraction, muscle deforms in three dimensions while maintaining constant volume (Ryan et al., 2020). Modelling studies suggest that limiting the radial expansion of a muscle, i.e., restricting transverse force, constrains longitudinal contraction (Gindre et al., 2013; Azizi et al., 2017; Ryan et al., 2020). Higher transverse muscle stiffness may therefore contribute to reduced muscle contraction levels. With the longitudinal-to-transverse stiffness ratio, we combined longitudinal and transverse stiffness measurements as a two-dimensional descriptor. In the OA group, the longitudinal-to-transverse ratio was lowest for both muscles across most tasks, suggesting reduced longitudinal tension relative to increased between-fiber stiffness. Interestingly, palpation and massage techniques are usually applied transversely to muscle fiber direction, implying an intuitive understanding of the importance of between-fiber elasticity.

The current work has important limitations. Most notably, the selection of the most painful muscles for assessment resulted in four different muscles being evaluated, dividing the sample into subgroups and reducing the statistical power of the study. Measuring the most painful muscles was mandatory for answering the research question and represented the clinically most relevant approach. More data are needed to confirm the observed tendencies. The current work enables sample size calculations for larger studies. Further, the sensation of stiffness in OA is typically reported when starting movement. Probably, we diminished group differences by the standardized warm-up procedure prior to data acquisition. We standardized prior activity because it influences muscle length (Proske and Gandevia, 2012), and therefore likely also muscle tension and stiffness. Moreover, the technical challenges associated with transverse stiffness measurements resulted in fewer data being included in the statistical analysis. This affected predominantly the transverse measurements of the vastus lateralis muscle. Only the use of elastography clips comprising 24–54 elastograms, from which measurements were averaged, enabled the high data inclusion and good to excellent reliability. In light of these limitations, the present findings should be interpreted with caution. Nonetheless, the observed trends, statistically significant differences and correlations suggest that the inclusion of transverse muscle stiffness measurements may contribute to a more comprehensive understanding of muscle mechanical properties and hold potential clinical relevance. Future studies should expand the sample size to further validate the current observations. This study did not include an assessment of biomechanical factors that may also influence muscle stiffness, such as lower limb alignment. Future research should consider incorporating biomechanical measurements to evaluate their potential impact on muscle stiffness.

5 Conclusion

In individuals with knee osteoarthritis, the stiffness of painful muscles was not elevated along the muscle fibers; rather, it was reduced during weight-bearing conditions. The inclusion of transverse muscle stiffness measurements appears to offer additional clinically relevant insights and may contribute to a more comprehensive assessment of muscle mechanics in this population. The sensation of stiffness did not align with measured stiffness.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: The raw data analyzed for this study can be found in ZENODO, doi 10.5281/zenodo.15025467 and 10.5281/zenodo.15025478.

Ethics statement

The studies involving humans were approved by Ethics Committee of the University of Stuttgart, No 21-012. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

AD: Project administration, Supervision, Data curation, Formal Analysis, Methodology, Funding acquisition, Writing – review and editing, Conceptualization, Investigation, Writing – original draft, Resources. KS: Writing – original draft, Software, Supervision, Data curation, Writing – review and editing. FP: Data curation, Writing – review and editing, Formal Analysis, Writing – original draft, Visualization. LG: Methodology, Investigation, Resources, Writing – original draft, Funding acquisition, Writing – review and editing. MA: Data curation, Software, Writing – original draft, Writing – review and editing. GC: Conceptualization, Writing – original draft, Writing – review and editing, Formal Analysis. MK: Writing – original draft, Visualization, Writing – review and editing, Formal Analysis, Validation, Supervision. AH: Methodology, Formal Analysis, Project administration, Data curation, Writing – original draft, Investigation, Writing – review and editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This study was funded by the Ministry of Science, Research, and Arts (MWK) Baden-Württemberg, Germany, Innovative Projekte 2020.

Acknowledgments

The authors acknowledge support by the Open Access Publication Fund of the University of Freiburg. The authors want to thank Franziska Bubeck for her support in the initial steps of this project.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that Gen AI was used in the creation of this manuscript exclusively to improve the language style and correct writing of the manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Akagi R., Kusama S. (2015). Comparison between neck and shoulder stiffness determined by shear wave ultrasound elastography and a muscle hardness meter. Ultrasound Med. and Biol. 41 (8), 2266–2271. doi:10.1016/j.ultrasmedbio.2015.04.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Alfuraih A. M., Tan A. L., O'Connor P., Emery P., Wakefield R. J. (2020). Muscle stiffness in rheumatoid arthritis is not altered or associated with muscle weakness: a shear wave elastography study. (1439–7609. doi:10.1080/14397595.2019.1645374

CrossRef Full Text | Google Scholar

Azizi E., Deslauriers A. R., Holt N. C., Eaton C. E. (2017). Resistance to radial expansion limits muscle strain and work. Biomech. Model. Mechanobiol. 16 (5), 1633–1643. doi:10.1007/s10237-017-0909-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Bahns C., Kopkow C. (2023). Physiotherapy for patients with hip and knee osteoarthritis in Germany: a survey of current practice. BMC Musculoskelet. Disord. 24 (1), 424. doi:10.1186/s12891-023-06464-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Bastick A. N., Wesseling J., Damen J., Verkleij S. P., Emans P. J., Bindels P. J., et al. (2016). Defining knee pain trajectories in early symptomatic knee osteoarthritis in primary care: 5-year results from a nationwide prospective cohort study (CHECK). Br. J. Gen. Pract. 66 (642), e32–e39. doi:10.3399/bjgp15X688129

PubMed Abstract | CrossRef Full Text | Google Scholar

Chino K., Kawakami Y., Takahashi H. (2017). Tissue elasticity of in vivo skeletal muscles measured in the transverse and longitudinal planes using shear wave elastography. Clin. Physiology Funct. imaging 37, 394–399. doi:10.1111/cpf.12315

PubMed Abstract | CrossRef Full Text | Google Scholar

Davis L. C., Baumer T. G., Bey M. J., van Holsbeeck M. (2019). Clinical utilization of shear wave elastography in the musculoskeletal system. Ultrasonography 38 (1), 2–12. doi:10.14366/usg.18039

PubMed Abstract | CrossRef Full Text | Google Scholar

Dieterich A. V., Yavuz U. S., Petzke F., Nordez A., Falla D. (2020). Neck muscle stiffness measured with shear wave elastography in women with chronic nonspecific neck pain. J. Orthop. Sports Phys. Ther. 50 (4), 179–188. doi:10.2519/jospt.2020.8821

PubMed Abstract | CrossRef Full Text | Google Scholar

Dor A., Kalichman L. (2017). A myofascial component of pain in knee osteoarthritis. J. Bodyw. Mov. Ther. 21 (3), 642–647. doi:10.1016/j.jbmt.2017.03.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Driban J. B., Harkey M. S., Barbe M. F., Ward R. J., MacKay J. W., Davis J. E., et al. (2020). Risk factors and the natural history of accelerated knee osteoarthritis: a narrative review. BMC Musculoskelet. Disord. 21 (1), 332. doi:10.1186/s12891-020-03367-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Eby S. F., Song P., Chen S., Chen Q., Greenleaf J. F., An K.-N. (2013). Validation of shear wave elastography in skeletal muscle. J. Biomechanics 46 (14), 2381–2387. doi:10.1016/j.jbiomech.2013.07.033

PubMed Abstract | CrossRef Full Text | Google Scholar

Friede M. C., Klauser A., Fink C., Csapo R. (2020). Stiffness of the iliotibial band and associated muscles in runner’s knee: assessing the effects of physiotherapy through ultrasound shear wave elastography. Phys. Ther. Sport 45, 126–134. doi:10.1016/j.ptsp.2020.06.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Gennisson J.-L., Deffieux T., Macé E., Montaldo G., Fink M., Tanter M. (2010). Viscoelastic and anisotropic mechanical properties of in vivo muscle tissue assessed by supersonic shear imaging. Ultrasound Med. Biol. 36 (5), 789–801. doi:10.1016/j.ultrasmedbio.2010.02.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Gibbs A. J., Gray B., Wallis J. A., Taylor N. F., Kemp J. L., Hunter D. J., et al. (2023). Recommendations for the management of hip and knee osteoarthritis: a systematic review of clinical practice guidelines. Osteoarthr. Cartil. 31 (10), 1280–1292. doi:10.1016/j.joca.2023.05.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Gihring A., Gärtner F., Liu C., Hoenicka M., Wabitsch M., Knippschild U., et al. (2020). Influence of obesity on the organization of the extracellular matrix and satellite cell functions after combined muscle and thorax trauma in C57bl/6J mice. Front. Physiol. 11, 849. doi:10.3389/fphys.2020.00849

PubMed Abstract | CrossRef Full Text | Google Scholar

Gindre J., Takaza M., Moerman K. M., Simms C. K. (2013). A structural model of passive skeletal muscle shows two reinforcement processes in resisting deformation. J. Mech. Behav. Biomed. Mater 22, 84–94. doi:10.1016/j.jmbbm.2013.02.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Gökşen A., Can F., Yilmaz S. (2021). Does muscle stiffness predict early-onset knee osteoarthritis? J. Surg. Med. 5 (11), 1126–1130. doi:10.28982/josam.1014433

CrossRef Full Text | Google Scholar

Green M. A., Geng G., Qin E., Sinkus R., Gandevia S. C., Bilston L. E. (2024). Measuring anisotropic muscle stiffness properties using elastography. NMR Biomed. 26, 1387–1394. doi:10.1002/nbm.2964

PubMed Abstract | CrossRef Full Text | Google Scholar

Haueise A., Le Sant G., Eisele-Metzger A., Dieterich A. V. (2023). Is musculoskeletal pain associated with increased muscle stiffness? - evidence map and critical appraisal of muscle measurements using Shear Wave Elastography. Clin. Physiology Funct. Imaging 44 (3), 187–204. doi:10.1111/cpf.12870

PubMed Abstract | CrossRef Full Text | Google Scholar

Hermens H. J., Freriks B., Disselhorst-Klug C., Rau G. (2000). Development of recommendations for SEMG sensors and sensor placement procedures. J. Electromyogr. Kinesiol 10, 361–374. doi:10.1016/s1050-6411(00)00027-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Hug F., Tucker K., Gennisson J. L., Tanter M., Nordez A. (2015). Elastography for muscle biomechanics: toward the estimation of individual muscle force. Exerc. Sports Sci. Rev. 43 (3), 125–133. doi:10.1249/JES.0000000000000049

PubMed Abstract | CrossRef Full Text | Google Scholar

Hvedstrup J., Tørring Kolding L., Ashina M., Winther Schytz H. (2020). Increased neck muscle stiffness in migraine patients with ictal neck pain: a shear wave elastography study. Cephalgia 40 (6), 565–574. doi:10.1177/0333102420919998

PubMed Abstract | CrossRef Full Text | Google Scholar

Hyldahl R. D., Nelson B., Xin L., Welling T., Groscost L., Hubal M. J., et al. (2015). Extracellular matrix remodeling and its contribution to protective adaptation following lengthening contractions in human muscle. Faseb J. 29 (7), 2894–2904. doi:10.1096/fj.14-266668

PubMed Abstract | CrossRef Full Text | Google Scholar

Järvinen T. A., Józsa L., Kannus P., Järvinen T. L., Järvinen M. (2002). Organization and distribution of intramuscular connective tissue in normal and immobilized skeletal muscles. An immunohistochemical, polarization and scanning electron microscopic study. J. Muscle Res. Cell. Motil. 23 (3), 245–254. doi:10.1023/a:1020904518336

PubMed Abstract | CrossRef Full Text | Google Scholar

Jozsa L., Thöring J., Järvinen M., Kannus P., Lehto M., Kvist M. (1988). Quantitative alterations in intramuscular connective tissue following immobilization: an experimental study in the rat calf muscles. Exp. Mol. Pathol. 49 (2), 267–278. doi:10.1016/0014-4800(88)90039-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Kessler S., Lang S., Puhl W., Stöve J. (2003). The Knee Injury and Osteoarthritis Outcome Score--a multifunctional questionnaire to measure outcome in knee arthroplasty. Z Orthop. Ihre Grenzgeb 141 (3), 277–282. doi:10.1055/s-2003-40083

PubMed Abstract | CrossRef Full Text | Google Scholar

Kjaer M., Magnusson P., Krogsgaard M., Boysen Møller J., Olesen J., Heinemeier K., et al. (2006). Extracellular matrix adaptation of tendon and skeletal muscle to exercise. J. Anat. 208 (4), 445–450. doi:10.1111/j.1469-7580.2006.00549.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Klauser A. S., de-Koekkoek F., Schwabl C., Fink C., Friede M., Csapo R. (2022). Intraobserver assessment of shear wave elastography in tensor fasciae latae and gluteus maximus muscle: the importance of the hip abductor muscles in runners knee compared to healthy controls. J. Clin. Med. 11 (13), 3605. doi:10.3390/jcm11133605

PubMed Abstract | CrossRef Full Text | Google Scholar

Kohl M. (2023). MKinfer: inferential statistics. Available online at: https://www.stamats.de.

Google Scholar

Kolb M., Ekert K., Schneider L., Fritz J., Ioanoviciu S. D., Henes J., et al. (2021a). The utility of shear-wave elastography in the evaluation of myositis. Ultrasound Med. Biol. 47 (8), 2176–2185. doi:10.1016/j.ultrasmedbio.2021.04.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Kolb M., Peisen F., Ekert K., Xenitidis T., Fritz J., Ioanoviciu S. D., et al. (2021b). Shear wave elastography for assessment of muscular abnormalities related to systemic sclerosis. Acad. Radiol. 28 (8), 1118–1124. doi:10.1016/j.acra.2020.04.043

PubMed Abstract | CrossRef Full Text | Google Scholar

Kolding L. T., Do T. P., Ewertsen C., Winter Schytz H. (2018). Muscle stiffness in tension-type headache patients with pericranial tenderness: a shear wave elastography study. Cephalgia Rep. 1, 1–6. doi:10.1177/2515816318760293

CrossRef Full Text | Google Scholar

Kolster B. (2016). Massage: klassische massage, querfriktionen, funktionsmassage, faszienbehandlung. Berlin: Springer. doi:10.1007/978-3-662-47273-6

CrossRef Full Text | Google Scholar

Kragstrup T. W., Kjaer M., Mackey A. L. (2011). Structural, biochemical, cellular, and functional changes in skeletal muscle extracellular matrix with aging. Scand. J. Med. Sci. Sports 21 (6), 749–757. doi:10.1111/j.1600-0838.2011.01377.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Li F., Wang Z. Y., Zhang Z. J., Shen S. H., Guo J. Y., Guo Y. X., et al. (2021). In hamstring muscles of patients with knee osteoarthritis an increased ultrasound shear modulus indicates a permanently elevated muscle tonus. Front. Physiol. 12, 752455. doi:10.3389/fphys.2021.752455

PubMed Abstract | CrossRef Full Text | Google Scholar

Li Z., Leung K. L., Huang C., Huang X., Chung R., Fu S. N. (2023). Passive stiffness of the quadriceps predicts the incidence of clinical knee osteoarthritis in twelve months. Eur. J. Phys. Rehabil. Med. 59 (1), 65–74. doi:10.23736/s1973-9087.22.07634-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Lieber R. L., Ward S. R. (2013). Cellular Mechanisms of Tissue Fibrosis. 4. Structural and functional consequences of skeletal muscle fibrosis. Am. J. Physiology Cell. Physiology 305 (3), C241–C252. doi:10.1152/ajpcell.00173.2013

PubMed Abstract | CrossRef Full Text | Google Scholar

Long H., Liu Q., Yin H., Wang K., Diao N., Zhang Y., et al. (2022). Prevalence trends of site-specific osteoarthritis from 1990 to 2019: findings from the global burden of disease study 2019. Arthritis and Rheumatology 74 (7), 1172–1183. doi:10.1002/art.42089

PubMed Abstract | CrossRef Full Text | Google Scholar

Marcucci L., Bondì M., Randazzo G., Reggiani C., Natali A. N., Pavan P. G. (2019). Fibre and extracellular matrix contributions to passive forces in human skeletal muscles: an experimental based constitutive law for numerical modelling of the passive element in the classical Hill-type three element model. PLoS One 14 (11), e0224232. doi:10.1371/journal.pone.0224232

PubMed Abstract | CrossRef Full Text | Google Scholar

Minato K., Yoshimoto Y., Kurosawa T., Watanabe K., Kawashima H., Ikemoto-Uezumi M., et al. (2021). Measurement of lateral transmission of force in the extensor digitorum longus muscle of young and old mice. Int. J. Mol. Sci. 22 (1422-0067 (Electronic)), 12356. doi:10.3390/ijms222212356

PubMed Abstract | CrossRef Full Text | Google Scholar

Miyamoto N., Hirata K., Kanehisa H., Yoshitake Y. (2015). Validity of measurement of shear modulus by ultrasound shear wave elastography in human pennate muscle. PLOS ONE 10 (4), e0124311. doi:10.1371/journal.pone.0124311

PubMed Abstract | CrossRef Full Text | Google Scholar

Ngo H. H. P., Andrade R. J., Lancelot J., Loumeaud A., Cornu C., Nordez A., et al. (2024). Unravelling anisotropic nonlinear shear elasticity in muscles: towards a non-invasive assessment of stress in living organisms. J. Mech. Behav. Biomed. Mater. 150, 106325. doi:10.1016/j.jmbbm.2023.106325

PubMed Abstract | CrossRef Full Text | Google Scholar

NICE (2023). Osteoarthritis in over 16s: diagnosis and management. NICE Guidel. NG226, 1–42. Available online at: www.nice.org.uk/guidance/ng226.

Google Scholar

Nicholls J., Alfuraih A. M., Hensor E. M. A., Robinson P. (2020). Inter- and intra-reader reproducibility of shear wave elastography measurements for musculoskeletal soft tissue masses. Skelet. Radiol. 49 (5), 779–786. doi:10.1007/s00256-019-03300-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Nordez A., Hug F. (2010). Muscle shear elastic modulus measured using supersonic shear imaging is highly related to muscle activity level. J. Appl. Physiology 108, 1389–1394. doi:10.1152/japplphysiol.01323.2009

PubMed Abstract | CrossRef Full Text | Google Scholar

Pedersen T. L. (2023). Patchwork: the composer of plots. Available online at: https://CRAN.R-project.org/package=patchwork.

Google Scholar

Portney L. G., Watkins M. P. (2000). Foundations of clinical research. Upper Saddle River, New Jersey 07458: Prentice-Hall Health.

Google Scholar

Proske U., Gandevia S. C. (2012). The proprioceptive senses: their roles in signaling body shape, body position and movement, and muscle force. Physiol. Rev. 92, 1651–1697. doi:10.1152/physrev.00048.2011

PubMed Abstract | CrossRef Full Text | Google Scholar

Purslow P. P. (2010). Muscle fascia and force transmission. J. Bodyw. Mov. Ther. 14 (4), 411–417. doi:10.1016/j.jbmt.2010.01.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Purslow P. P. (2020). The structure and role of intramuscular connective tissue in muscle function. Front. Physiol. 11, 495. doi:10.3389/fphys.2020.00495

PubMed Abstract | CrossRef Full Text | Google Scholar

Purslow P. P., Trotter J. A. (1994). The morphology and mechanical properties of endomysium in series-fibred muscles: variations with muscle length. J. Muscle Res. Cell. Motil. 15 (3), 299–308. doi:10.1007/bf00123482

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramaswamy K. S., Palmer M. L., Van der Meulen J. H., Renoux A., Kostrominova T. Y., Michele D. E., et al. (2011). Lateral transmission of force is impaired in skeletal muscles of dystrophic mice and very old rats. J. Physiology 589, 1195–1208. doi:10.1113/jphysiol.2010.201921

PubMed Abstract | CrossRef Full Text | Google Scholar

Ruby L., Mutschler T., Martini K., Klingmüller V., Frauenfelder T., Rominger M. B., et al. (2019). Which confounders have the largest impact in shear wave elastography of muscle and how can they be minimized? An elasticity phantom, ex vivo porcine muscle and volunteer study using a commercially available system. Ultrasound Med. and Biol. 45 (10), 2591–2611. doi:10.1016/j.ultrasmedbio.2019.06.417

PubMed Abstract | CrossRef Full Text | Google Scholar

Rusu A. C., Kreddig N., Hallner D., Hülsebusch J., Hasenbring M. I. (2014). Fear of movement/(Re)injury in low back pain: confirmatory validation of a German version of the Tampa Scale for Kinesiophobia. BMC Musculoskelet. Disord. 15, 280. doi:10.1186/1471-2474-15-280

PubMed Abstract | CrossRef Full Text | Google Scholar

Ryan D. S., Domínguez S., Ross S. A., Nigam N., Wakeling J. M. (2020). The energy of muscle contraction. II. Transverse compression and work. Front. Physiol. 11, 538522. doi:10.3389/fphys.2020.538522

PubMed Abstract | CrossRef Full Text | Google Scholar

Săftoiu A., Gilja O. H., Sidhu P. S., Dietrich C. F., Cantisani V., Amy D., et al. (2019). The EFSUMB guidelines and recommendations for the clinical practice of elastography in non-hepatic applications: update 2018. Ultraschall der Medizin 40, 425–453. doi:10.1055/a-0838-9937

PubMed Abstract | CrossRef Full Text | Google Scholar

Sinha U., Malis V., Chen J. S., Csapo R., Kinugasa R., Narici M. V., et al. (2020). Role of the extracellular matrix in loss of muscle force with age and unloading using magnetic resonance imaging, biochemical analysis, and computational models. Front. Physiol. 11, 626. doi:10.3389/fphys.2020.00626

PubMed Abstract | CrossRef Full Text | Google Scholar

Skerl K., Cochran S., Evans A. (2017). First step to facilitate long-term and multi-centre studies of shear wave elastography in solid breast lesions using a computer-assisted algorithm. Int. J. Comput. Radiology Surg. 9, 1533–1542. doi:10.1007/s11548-017-1596-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Spadoni S., Todros S., Reggiani C., Marcucci L., Pavan P. G. (2024). The role of the extracellular matrix in the reduction of lateral force transmission in muscle bundles: a finite element analysis. Comput. Biol. Med. 1745, 108488–110534. doi:10.1016/j.compbiomed.2024.108488

PubMed Abstract | CrossRef Full Text | Google Scholar

Team R. C. (2023). in R: a language and environment for statistical computing. Editor R. C. Team (R Foundation for Statistical Computing). (Vienna). Available online at: https://www.R-project.org/

Google Scholar

Teo P. L., Hinman R. S., Egerton T., Dziedzic K. S., Bennell K. L. (2019). Identifying and prioritizing clinical guideline recommendations most relevant to physical therapy practice for hip and/or knee osteoarthritis. Jornal Orthop. Sports Phys. Ther. 49 (7), 501–512. doi:10.2519/jospt.2019.8676

PubMed Abstract | CrossRef Full Text | Google Scholar

von Elm E., Altman D. G., Egger M., Pocock S. J., Gøtzsche P. C., Vandenbroucke J. P., et al. (2007). The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement: guidelines for reporting observational studies. PLoS Med. 4 (10), e296. doi:10.1371/journal.pmed.0040296

PubMed Abstract | CrossRef Full Text | Google Scholar

Vuorenmaa A. S., Siitama E.A.-O., Mäkelä K. S. (2022). Inter-operator and inter-device reproducibility of shear wave elastography in healthy muscle tissues. J. Appl. Clin. Med. Phys. 23, e13717. doi:10.1002/acm2.13717

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang A. B., Perreault E. J., Royston T. J., Lee S. S. M. (2019). Changes in shear wave propagation within skeletal muscle during active and passive force generation. J. Biomech. 94, 115–122. doi:10.1016/j.jbiomech.2019.07.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang C., Cui L., Xu Y., Xue H., Sun Y., Fu Q., et al. (2022). Ultrasound and shear-wave elastography characteristics of the iliocapsularis muscle in pincer-type femoroacetabular impingement. Pm R. 14 (11), 1360–1367. doi:10.1002/pmrj.12707

PubMed Abstract | CrossRef Full Text | Google Scholar

Wheatley B. B. (2020). Investigating passive muscle mechanics with biaxial stretch. Front. Physiol. 11, 1021. doi:10.3389/fphys.2020.01021

PubMed Abstract | CrossRef Full Text | Google Scholar

Wickham H. (2016). ggplot2: elegant graphics for data analysis. New York: Springer.

Google Scholar

Wickham H., Francois R., Henry L., Müller K., Vaughan D. (2023a). Dplyr: a grammar of data manipulation. Available online at: https://CRAN.R-project.org/package=dplyr.

Google Scholar

Wickham H., Vaughan D., Girlich M. (2023b). Tidyr: tidy messy data. Available online at: https://CRAN.R-project.org/package=tidyr.

Google Scholar

Wohlgemuth R. P., Feitzinger R. M., Henricson K. E., Dinh D. T., Brashear S. E., Smith L. R. (2023). The extracellular matrix of dystrophic mouse diaphragm accounts for the majority of its passive stiffness and is resistant to collagenase digestion. Matrix Biol. Plus 18, 100131–100285. (Electronic)), 100131. doi:10.1016/j.mbplus.2023.100131

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: shear wave elastography, muscle, stiffness, knee osteoarthritis, anisotropy, age, patient reported outcome measures

Citation: Dieterich AV, Skerl K, Paskali F, Gizzi L, Azan M, Carvalho GF, Kohl M and Haueise A (2025) Longitudinal and transverse muscle stiffness change differently with knee osteoarthritis and do not align with stiffness sensation. Front. Physiol. 16:1593851. doi: 10.3389/fphys.2025.1593851

Received: 14 March 2025; Accepted: 06 May 2025;
Published: 27 May 2025.

Edited by:

Elisa Belluzzi, University of Padua, Italy

Reviewed by:

Naveen Kumar Potti, Naveen institute of manual therapy, India
Zongpan Li, University of Maryland, United States

Copyright © 2025 Dieterich, Skerl, Paskali, Gizzi, Azan, Carvalho, Kohl and Haueise. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Andreas Haueise, YW5kcmVhcy5oYXVlaXNlQGVtYWlsLnVuaS1mcmVpYnVyZy5kZQ==

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.