Ultrasound elastography in the assessment of post-stroke muscle stiffness: a systematic review

Anisotropic mechanical properties Quantification in skeletal muscle using magnetic resonance elastography and diffusion tensor imaging

Skeletal muscle contains a highly hierarchical structure, leading to anisotropic mechanical properties, with varying morphological responses to mechanical loadings from different directions. However, this feature is rarely studied in clinical studies, mainly due to the challenges in quantifying muscle anisotropic mechanical properties in vivo. The aim of the current study was to quantify the anisotropic mechanical properties of skeletal muscle using an integrated approach combining multi-frequency magnetic resonance elastography (MRE) and diffusion tensor imaging (DTI). Muscle fascicle orientation was determined through DTI tractography. Direct inversion of the curl-based wave equation was used to quantify three complex-valued moduli (μ⊥∗, μ‖∗, and E‖∗) assuming muscle as an incompressible transversely isotropic material. This approach was evaluated on one ex vivo muscle sample by comparing MRE-derived moduli to rheometry measurements, and further assessed in vivo in the ankle plantarflexors of nine able-bodied subjects. Consistency in the anisotropic ratio was observed between rheometry and MRE measurements in the ex vivo muscle sample, though discrepancies were noted in absolute shear moduli values. In vivo, the anisotropy of skeletal muscle was observed by the relationship of μ⊥∗≠1/3E‖∗ and μ‖∗≠1/3E‖∗ at different MRE driving frequencies with higher parallel shear modulus (μ‖∗) than the perpendicular shear modulus (μ⊥∗). This study demonstrated a promising approach for quantifying the muscle anisotropic mechanical properties in vivo, which can be useful in various clinical applications.

Simultaneous measurement of tensile and shear elasticity and anisotropy in human skeletal muscle tissue using steered ultrasound shear waves

Load-bearing skeletal muscle tissues are reinforced by intricate networks of protein fibers aligned in preferential orientations, imparting direction-dependent mechanical properties (anisotropy). Characterizing this anisotropy in vivo is essential for understanding both normal and pathological muscle function, as well as structural integrity. However, current noninvasive techniques are limited in their ability to accurately measure the mechanical properties of anisotropic tissues such as skeletal muscle. Here, we used an innovative angle-resolved ultrasound elastography method, recently developed by our team, to simultaneously quantify tensile and shear elasticity and anisotropy, enabling comprehensive assessment of muscle biomechanics. We fully characterized the mechanical properties of the biceps brachii in fourteen healthy young adults under passive and active axial loadings, revealing distinct shear and tensile mechanical behaviors both along and across muscle fibers. Notably, tensile and shear moduli along the main fiber orientation were found to be uncoupled, while cross-muscle fiber measurements exhibited a consistent modulus ratio of 3.4 ± 0.2, regardless of axial loading conditions or intensities. These findings highlight the anisotropic nature of skeletal muscle and provide valuable insights into its in vivo mechanical behavior. Both shear and tensile anisotropy increased with muscle axial physiological loading, indicating that our method is sensitive to changes in anisotropic tissue mechanics. Lastly, we demonstrated that angle-resolved ultrasound shear wave elastography provides direct estimates of shear and tensile properties, offering significant promise for clinical applications, including neuromuscular disease diagnostics and monitoring, biomechanical modeling for predicting tissue responses to loading and therapies, and tissue engineering. STATEMENT OF SIGNIFICANCE: : Conventional ultrasound shear wave elastography techniques overlook the anisotropy of skeletal muscles, leading to incomplete tissue mechanical characterization. In this study, we leveraged an innovative angle-resolved elastography method to assess tensile and shear elasticity, along with their anisotropic factors, of human muscle in vivo. For the first time, we reveal the intricate relationships between tensile and shear elasticities during active and passive physiological loading. This technique enhances our understanding of muscle mechanics and has promising clinical implications for muscle health and neuromuscular disease management, where tissue structural and mechanical properties are often altered. Additionally, it could help in developing constitutive models for muscle tissue and contribute to the design of tissue-engineered materials.

Feasibility of shear wave elastography to assess upper limb spasticity in patients after brain injury

PURPOSE

To assess the feasibility of Shear Wave Elastography (SWE) for quantifying upper limb spasticity in people after brain injury.

METHODS

This observational study included 52 patients with upper limb spasticity after brain injury. Participants underwent a two-week rehabilitation treatment. The modified Ashworth scale (MAS), upper extremity section of the Fugl-Meyer motor scale (FMMS-UE), Barthel index (BI) and SWE were assessed before and after a two-week rehabilitation. The shear modulus was used to evaluate the stiffness of the biceps brachii at the transverse and longitudinal plane of 0°and 90°elbow flexion.

RESULTS

In our study, paretic side shear modulus values were significantly larger than non-paretic sides (all p-values < 0.05), with significant reductions only on the paretic side (all p-values < 0.001) after 2-week rehabilitation treatment. At baseline, MAS and shear modulus had a significant positive correlation. After rehabilitation, only the longitudinal plane shear modulus showed a correlation with MAS. No significant correlations were found between shear modulus and FMMS-UE or BI scores, except for a negative correlation between shear modulus and BI at 0° transverse plane flexion.

CONCLUSIONS

Our findings indicate that SWE effectively quantifies upper limb spasticity and can assess the impact of rehabilitation treatment.

Shear wave elastography in the assessment of gastrocnemius spastic muscle elasticity: influences of ankle position and muscle contraction

BACKGROUND

Following upper motor neuron syndromes (UMNS), intrinsic viscoelastic muscle properties such as elastic stiffness may be altered, which leads to muscle hyper-resistance to passive mobilization. So far, no gold standard assessment of hyper-resistance, whether clinical or instrumental, is available. Shear wave elastography (SWE) has been increasingly used for non-invasive evaluation of elastic stiffness of skeletal muscles in people with hyper-resistance.

AIM

Our study aimed to evaluate the validity of SWE at ankle neutral resting position (ANRP). Additional objectives included assessing the influence of ankle position, muscle contraction, and laterality on elastic stiffness, and evaluating the reproducibility and responsiveness of SWE.

DESIGN

Observational cohort study.

SETTING

Outpatients of the physical and rehabilitation medicine department of a university hospital in Brussels.

POPULATION

Thirty hemiparetic patients following UMNS with hyper-resistance in gastrocnemii muscles.

METHODS

Elastic stiffness was quantified by shear wave velocity (SWV) measurements of gastrocnemii muscles. A higher SWV corresponds to a higher elastic stiffness. Measurements were performed on the affected and less-affected limbs in ANRP, in passive dorsiflexion and during isometric contraction. Assessments were performed 3 times. Criterion validity, reproducibility, and responsiveness were evaluated. A linear mixed model was used to study position and laterality effect.

RESULTS

In ANRP, reproducibility was excellent and SWV was significantly higher in the affected limb than in the less-affected limb. This laterality effect disappeared in passive dorsiflexion and was even reversed during isometric contraction. SWV was significantly higher on both sides in passive dorsiflexion and during contraction than in ANRP.

CONCLUSIONS

Our results suggest that SWE measurements in ANRP are reliable and may provide a more valid measurement of gastrocnemii elastic stiffness following UMNS.

CLINICAL REHABILITATION IMPACT

SWE may be a useful clinical tool as an extension of the physical exam for longitudinal monitoring of passive muscle elastic stiffness, to assist with treatment decisions and to better quantify the therapeutic effect of procedures to reduce muscle overactivity. However, a standardized protocol should be used. ANRP seems to be the most valid position for assessing gastrocnemius elastic stiffness in neurological populations. This should be kept in mind for the choice of positioning in further studies.

Effect of Ankle-Foot Orthosis on Paretic Gastrocnemius and Tibialis Anterior Muscle Contraction of Stroke Survivors During Walking: A Pilot Study

Ankle-foot orthoses (AFOs) have been commonly prescribed for stroke survivors with foot drop, but their impact on the contractions of paretic tibialis anterior (TA) and medial gastrocnemius (MG) has remained inconclusive. This study thus investigated the effect of AFOs on these muscle contractions in stroke survivors. The contractions of paretic TA and MG muscles were assessed in twenty stroke patients and compared between walking with and without AFOs, using a novel wearable dynamic ultrasound imaging and sensing system. The study found an increase in TA muscle thickness throughout a gait cycle (p > 0.05) and a significant increase in TA muscle surface mechanomyography (sMMG) signals during the pre- and initial swing phases (p < 0.05) when using an AFO. MG muscle thickness generally decreased with the AFO (p > 0.05), aligning more closely with trends seen in healthy adults. The MG surface electromyography (sEMG) signal significantly decreased during the initial and mid-swing phases when wearing an AFO (p < 0.05). The TA-MG co-contraction index significantly decreased during initial and mid-swing phases with the AFO (p < 0.05). These results suggest that AFOs positively influenced the contraction patterns of paretic ankle muscles during walking in stroke patients, but further research is needed to understand their long-term effects.

Aging-associated differences between perioral and trunk muscle characteristics

The physiological and functional changes from aging affect the systemic and swallowing muscles. While ultrasound is used to examine muscle characteristics, no previous studies have examined the relationship between age and perioral and trunk muscle characteristics, including stiffness. This study aimed to investigate the relationship between age and the characteristics of perioral and trunk muscles, including quality, quantity, and stiffness, using ultrasound. A total of 215 participants were recruited for this cross-sectional study: 98 younger adults (20-64 years) and 117 older adults (≥ 65 years). Muscle characteristics were assessed along with data on age, sex, body mass index, and number of teeth. The correlation between age and muscle characteristics was examined. Multiple linear regression was performed with muscle characteristics as dependent variables and age, body mass index, and number of teeth as explanatory variables. We found that only trunk muscle stiffness was associated with age, with a stronger correlation in quality and quantity than perioral muscles. Masseter muscle quality and quantity showed a higher correlation with age than those of the geniohyoid muscle. In conclusion, characteristics of the trunk and perioral muscles exhibit differences in age-related changes, highlighting the importance of maintaining both trunk and perioral muscle activity in older adults to prevent functional decline.

Quantifying Plantar Flexor Muscles Stiffness During Passive and Active Force Generation Using Shear Wave Elastography in Individuals With Chronic Stroke

OBJECTIVES

This study aims to investigate the mechanical properties of paretic and healthy plantar flexor muscles and assesses the spatial distribution of stiffness between the gastrocnemius medialis (GM) and lateralis (GL) during active force generation.

METHODS

Shear wave elastography measurements were conducted on a control group (CNT, n=14; age=59.9±10.6 years; BMI=24.5±2.5 kg/m2) and a stroke survivor group (SSG, n=14; age=63.2±9.6 years; BMI=23.2±2.8 kg/m2). Shear modulus within the GM and GL was obtained during passive ankle mobilization at various angles of dorsiflexion (P0 =0°; P1=10°; P2=20°; P3=-20° and P4=-30°) and during different levels (30%, 50%, 70%, 100%) of maximal voluntary contraction (MVC). Muscle activations of GM, GL, soleus and tibialis anterior were also evaluated.

RESULTS

The results revealed a significant increase in passive stiffness within the paretic plantar flexor muscles under high tension during passive mobilization (p<0.05). Yet, during submaximal and maximal MVC, the paretic plantar flexors exhibited decreased active stiffness levels (p<0.05). A notable discrepancy was found between the stiffness of the GM and GL, with the GM demonstrating greater stiffness from 0° of dorsiflexion in the SSG (p<0.05), and from 10° of dorsiflexion in the CNT (p<0.05). No significant difference in stiffness was observed between the GM and GL muscles during active condition.

CONCLUSION

Stroke affects the mechanical properties differently depending on the state of muscle activation. Notably, the distribution of stiffness among synergistic plantar flexor muscles varied in passive condition, while remaining consistent in active condition.

Acute effects of static and proprioceptive neuromuscular facilitation stretching on hamstrings muscle stiffness and range of motion: a randomized cross-over study

This study aimed to compare the acute effects of static stretching (SS) and proprioceptive neuromuscular facilitation (PNF) stretching on hamstrings flexibility and shear modulus. Sixteen recreationally active young volunteers participated in a randomized cross-over study. Participants underwent an aerobic warm-up (WU), followed by either SS or PNF stretching. Range of motion (RoM) during passive straight leg raise and active knee extension, as well as shear modulus of the biceps femoris (BF) and semitendinosus (ST) muscles, were measured at baseline, post-WU, and post-stretching. Both stretching techniques significantly increased RoM, with no differences observed between SS and PNF (p < 0.001; η2 = 0.59-0.68). However, only PNF stretching resulted in a significant decrease in BF shear modulus (time×stretching type interaction: p = 0.045; η2 = 0.19), indicating reduced muscle stiffness. No changes in ST shear modulus were observed after either stretching technique. There was no significant correlation between changes in RoM and shear modulus, suggesting that the increase in RoM was predominantly due to changes in stretch tolerance rather than mechanical properties of the muscles. These findings suggest that both SS and PNF stretching can effectively improve hamstring flexibility, but PNF stretching may additionally reduce BF muscle stiffness. The study highlights the importance of considering individual muscle-specific responses to stretching techniques and provides insights into the mechanisms underpinning acute increases in RoM.

Spatial distribution of stiffness between and within muscles in paretic and healthy individuals during prone and standing positions

This study investigated the inter- and intramuscular variability of plantar flexors stiffness during prone and standing positions at different muscle lengths in healthy and paretic individuals. To access tissue stiffness, shear wave elastography (SWE) measurements were carried out on two groups: control group (CG; n=14; age 43.9±9.6 years; body mass index [BMI]=24.5±2.5 kg/m2) and stroke survivor group (SSG; n=14; age 43.9±9.6 years; BMI=24.5±2.5 kg/m2). Shear Modulus (μ, kPa) within three plantar flexors (the gastrocnemius medialis [GM], gastrocnemius lateralis [GL], and soleus [SOL]) was obtained during two conditions: prone and standing position, at different angles of dorsiflexion (0°, 10°, and 20°). Measurements were also performed in different proximo-distal regions of each muscle. Muscle activation of the GM, GL, SOL, and tibialis anterior were evaluated during the two conditions. Results showed a high spatial stiffness variability between and within plantar flexors during dorsiflexion. The highest stiffness was observed in the GM, especially in the distal region at 20° in healthy and paretic muscles. In the prone position, the paretic muscle exhibits greater stiffness compared to the healthy muscle (p < 0.05). In contrast, in the standing position, an increase of stiffness in the healthy muscle compared to the paretic muscle was observed (p < 0.05). Thus, mechanical properties are differently affected by stroke depending on active and passive states of ankle muscles during dorsiflexion. In addition, the modification of ankle muscle state change stiffness distribution between and within plantar flexors.

Anisotropy and reproducibility of ultrasound shear wave elastography in patella tendons with and without tendinopathy

PURPOSE

Ultrasound shear wave elastography (SWE) is a tool that can be utilized to assess biomechanical properties of tendons. Anisotropy, an ultrasound imaging artifact has been commonly cited as a potential source of error in the accuracy and reproducibility of SWE. The aim of the study was to assess reproducibility in performing SWE of patella tendons and differences in SWE and anisotropy between normal patella tendons and patellar tendinopathy.

METHODS

After obtaining the Institutional Review Board approval and written informed consent, we prospectively measured the shear wave velocity (SWV) of patella tendons with and without tendinopathy in 25 volunteers. SWVs were measured in three anatomic planes: longitudinal, perpendicular transverse, and tilted transverse with the probe tilted 15-30° from the perpendicular transverse plane by three operators with varied levels of experience. Anisotropy coefficient (A) was calculated by formula of A = (SWVLongitudinal - SWVTransverse) / SWVTransverse.

RESULTS

Differences in SWV and anisotropy coefficient between normal tendons and tendons with tendinopathy were significant (p < 0.05). The intra- and inter-observer reproducibility in performing SWE were moderate to good (intraclass correlation coefficient: 0.81-0.95). The mean difference of 95% Bland-Altman limits of agreement for measuring tendon SWV ranged -0.08 to 0.41 (upper 0.08 to 1.14, lower -1.22 to -0.22) between senior and junior operators.

CONCLUSION

The results of this study suggest that SWE and anisotropy coefficient are feasible tools to differentiate patellar tendinopathy from normal patella tendons. The reproducibility of performing SWE of patella tendons is moderate to good.

Whole-Body Photobiomodulation Therapy for Fibromyalgia: A Feasibility Trial

Effective treatment for fibromyalgia (FM) is lacking and further treatment options are needed. Photobiomodulation therapy (PBMT) represents one potential treatment option. Whilst favourable findings have been reported using localised PBMT, no investigations have established the value of whole-body PBMT for the complete set of symptom domains in FM. A single-arm feasibility study was conducted in accordance with CONSORT (Consolidated Standards of Reporting Trials) guidelines. A non-probability sampling method was used to access individuals with FM. The primary outcome measure was identified as the Revised Fibromyalgia Impact Questionnaire (FIQR). Forty-nine participants were screened and twenty-one trial participants entered the trial. Nineteen participants completed the intervention (18 whole-body PBMT sessions over approximately six weeks). Descriptive statistics and qualitative analysis was undertaken to represent feasibility outcomes. Acceptability of the trial device and processes were established. Outcome measures towards efficacy data were guided by core and peripheral OMERACT (outcomes measures in rheumatological clinical trials) domains, utilising a combination of participant-reported and performance-based outcome measures. Data for the embedded qualitative component of the trial were captured by participant-reported experience measures and audio-recorded semi-structured interviews. Positive changes were observed for FM-specific quality of life, pain, tenderness, stiffness, fatigue, sleep disturbance, anxiety, depression and cognitive impairment. Patient global assessment revealed improvements at 6 weeks, with continued effect at 24 weeks. FM-specific quality of life at 24 weeks remained improved compared with baseline scores. The findings provided evidence to support a full-scale trial and showed promise regarding potential efficacy of this novel non-invasive treatment in an FM population.

Utilization of shear wave elastography to quantify and predict response to upper extremity botulinum toxin injections in patients with cerebral palsy: A pilot study

OBJECTIVE

Shear wave elastography (SWE) was used to quantify change in upper extremity muscle stiffness in patients with unilateral spastic cerebral palsy (USCP) following botulinum toxin A (BTX-A) therapy. We hypothesized that SWE measures would decrease following ultrasound-guided BTX-A injection, and correlate with functional improvement.

METHODS

SWE measures of BTX-A treated muscles were recorded immediately pre-injection, and at 1-, 3- and 6-months post-injection. At the same timepoints, functional assessment was performed using the Modified Ashworth Scale (MAS), and passive and active range of motion (PROM and AROM) measures. Correlation of SWE with MAS, PROM and AROM, as well as the relationship between change in SWE and change in MAS, PROM and AROM was determined using Spearman's rank correlation coefficient and generalized estimating equation modeling.

RESULTS

16 muscles were injected and longitudinally assessed. SWE and MAS scores decreased following BTX-A injection (p = 0.030 and 0.004, respectively), reflecting decreased quantitative and qualitative muscle stiffness. Decreased SWE reached statistical significance at 1- and 3-months, and 1-, 3- and 6-months for MAS. When comparing relative change in SWE to relative change in AROM, larger change in SWE strongly correlated with positive change in AROM (p-value range:<0.001-0.057). BTX-A responders also demonstrated lower baseline SWE (1.4 m/s) vs. non-responders (1.9 m/s), p = 0.035.

CONCLUSION

Ultrasound-guided BTX-A injections in patients with USCP resulted in decreased quantitative and qualitative muscle stiffness. Strong correlation between change in SWE and AROM, as well as the significant difference in baseline SWE for BTX-A responders and non-responders, suggests SWE may provide a useful tool to predict and monitor BTX-A response.

Correlation between Body Composition and Inter-Examiner Errors for Assessing Lumbar Multifidus Muscle Size, Shape and Quality Metrics with Ultrasound Imaging

Ultrasound imaging (US) is widely used in several healthcare disciplines (including physiotherapy) for assessing multiple muscle metrics such as muscle morphology and quality. Since measuring instruments are required to demonstrate their reliability, accuracy, sensitivity, and specificity prior to their use in clinical and research settings, identifying factors affecting their diagnostic accuracy is essential. Since previous studies analyzed the impact of sociodemographic but not body composition characteristics in US errors, this study aimed to assess whether body composition metrics are correlated with ultrasound measurement errors. B-mode images of the lumbar multifidus muscle at the fifth lumbar vertebral level (L5) were acquired and analyzed in 49 healthy volunteers by two examiners (one experienced and one novel). Cross-sectional area, muscle perimeter and mean echo intensity were calculated bilaterally. A multivariate correlation matrix was calculated for assessing the inter-examiner differences with body composition metrics. Results demonstrated excellent reliability (intraclass correlation coefficient, ICC > 0.9) for assessing the muscle cross-sectional area and perimeter, and good reliability for assessing the muscle shape and mean echo intensity (ICC > 0.7). Inter-examiner errors for estimating muscle size were correlated with participants' age (p value, p < 0.01), weight (p < 0.05), total and trunk lean mass (both, p < 0.01) and water volume (p < 0.05). Greater shape descriptors and mean brightness disagreements were correlated with older ages (p < 0.05) and total lean mass (p < 0.05). No correlations between age and body composition metrics were found (p > 0.05). This study found US to be a reliable tool for assessing muscle size, shape and mean brightness. Although aging showed no correlations with body composition changes in this sample, it was the main factor correlated with US measurement errors.

Morphological and Mechanical Properties of Lower-Limb Muscles in Type 2 Diabetes: New Potential Imaging Indicators for Monitoring the Progress of DPN

The aim of this study was to explore changes in morphological and mechanical properties of lower-limb skeletal muscles in patients with diabetes with and without diabetic peripheral neuropathy (DPN) and seek to find a potential image indicator for monitoring the progress of DPN in patients with type 2 diabetes mellitus (T2DM). A total of 203 patients with T2DM, with and without DPN, were included in this study. Ultrasonography and ultrasound shear wave imaging (USWI) of the abductor hallux (AbH), tibialis anterior (TA), and peroneal longus (PER) muscles were performed for each subject, and the shear wave velocity (SWV) and cross-sectional area (CSA) of each AbH, TA, and PER were measured. The clinical factors influencing AbH_CSA and AbH_SWV were analyzed, and the risk factors for DPN complications were investigated. AbH_CSA and AbH_SWV in the T2DM group with DPN decreased significantly (P < 0.05), but no significant differences were found in the SWV and CSA of the TA and PER between the two groups. Toronto Clinical Scoring System (CSS) score and glycosylated hemoglobin (HbA1c) were independent predictors of AbH_CSA and AbH_SWV. As AbH_SWV and AbH_CSA decreased, Toronto CSS score and HbA1c increased and incidence of DPN increased significantly. In conclusion, the AbH muscle of T2DM patients with DPN became smaller and softer, while its morphological and mechanical properties were associated with the clinical indicators related to the progression of DPN. Thus, they could be potential imaging indicators for monitoring the progress of DPN in T2DM patients.