Measurement of passive skeletal muscle mechanical properties in vivo: recent progress, clinical applications, and remaining challenges

Test-retest reliability and follow-up of muscle magnetic resonance elastography in adults with and without muscle diseases

BACKGROUND

We investigated the potential of magnetic resonance elastography (MRE) stiffness measurements in skeletal muscles as an outcome measure, by determining its test-retest reliability, as well as its sensitivity to change in a longitudinal follow-up study.

METHODS

We assessed test-retest reliability of muscle MRE in 20 subjects with (n = 5) and without (n = 15) muscle diseases and compared this to Dixon proton density fat fraction (PDFF) and volume measurements. Next, we measured MRE muscle stiffness in 21 adults with Becker muscular dystrophy (BMD) and 21 age-matched healthy controls at baseline, and after 9 and 18 months. We compared two different methods of analysing MRE data in this study: 'Method A' used the stiffness maps generated by the Philips MRE software, and 'Method B' applied a custom-made procedure based on wavelength measurements on the MRE images.

RESULTS

Intraclass correlation coefficients (ICC) of muscle stiffness ranged from good (0.83 for left vastus medialis, P < 0.001) to poor (0.19 for right rectus femoris, P = 0.212) for the examined thigh muscles with Method A, but we did not find a significant test-retest reliability with Method B (P > 0.050 for all). The ICC of muscle PDFF and volume measurements was excellent (>0.90; P < 0.001) for all muscles. At baseline, the average stiffness of all thigh muscles was significantly lower in adults with BMD than in controls for both Method A (-0.2 kPa, P = 0.025) and Method B (-0.6 kPa, P < 0.001). Regardless of which method was used, there was no significant difference in the evolution of muscle stiffness in patients and controls over 18 months.

CONCLUSIONS

Test-retest reliability of muscle MRE using a simple 2D technique was suboptimal, and did not reliably measure muscle stiffness changes in adults with BMD as compared with controls over 18 months. While the results provide motivation for testing more advanced 3D MRE methods, we conclude that the simple 2D MRE implementation used in this study is not suitable as an outcome measure for characterizing thigh muscle in clinical trials.

The Biomechanics of Musculoskeletal Tissues during Activities of Daily Living: Dynamic Assessment Using Quantitative Transmission-Mode Ultrasound Techniques

The measurement of musculoskeletal tissue properties and loading patterns during physical activity is important for understanding the adaptation mechanisms of tissues such as bone, tendon, and muscle tissues, particularly with injury and repair. Although the properties and loading of these connective tissues have been quantified using direct measurement techniques, these methods are highly invasive and often prevent or interfere with normal activity patterns. Indirect biomechanical methods, such as estimates based on electromyography, ultrasound, and inverse dynamics, are used more widely but are known to yield different parameter values than direct measurements. Through a series of literature searches of electronic databases, including Pubmed, Embase, Web of Science, and IEEE Explore, this paper reviews current methods used for the in vivo measurement of human musculoskeletal tissue and describes the operating principals, application, and emerging research findings gained from the use of quantitative transmission-mode ultrasound measurement techniques to non-invasively characterize human bone, tendon, and muscle properties at rest and during activities of daily living. In contrast to standard ultrasound imaging approaches, these techniques assess the interaction between ultrasound compression waves and connective tissues to provide quantifiable parameters associated with the structure, instantaneous elastic modulus, and density of tissues. By taking advantage of the physical relationship between the axial velocity of ultrasound compression waves and the instantaneous modulus of the propagation material, these techniques can also be used to estimate the in vivo loading environment of relatively superficial soft connective tissues during sports and activities of daily living. This paper highlights key findings from clinical studies in which quantitative transmission-mode ultrasound has been used to measure the properties and loading of bone, tendon, and muscle tissue during common physical activities in healthy and pathological populations.

Upper extremity soft tissue stiffness measured with shear wave elastography: A between-day reliability study

Background

Imbalances in soft tissue stiffness (STS) of the pectoralis and common wrist extensors (CWEs) can adversely alter upper extremity function. Shear wave elastography (SWE) has the capacity to provide precise, repeatable, objective data on STS. The purpose of this study was to determine the between-day intrarater reliability for the pectoralis and CWE tendon stiffness as measured by SWE.

Methods

STS measured by the shear wave modulus (kilopascals (kPa)) was captured bilaterally on the pectoralis major, pectoralis minor, and CWE tendon using 2D SWE ultrasound imaging (GE Logiq S8, 9 L transducer) on two separate days. Within examiner intraclass correlation coefficients (ICC (3, 3)) with 95% confidence intervals, standard error of the measure, and minimal detectable changes were calculated.

Results

The investigators recruited 10 healthy participants (mean age = 23.50  ±  1.43). Between day intrarater reliability values were obtained for the dominant pectoralis major ICC  =  0.70 (-0.15, -0.93), nondominant pectoralis major ICC  =  0.63 (-0.30, -0.91), dominant pectoralis minor ICC  = 0 .72(-0.04, -0.93), nondominant pectoralis minor ICC  = 0.75(0.08, -0.94), dominant CWE tendon ICC  =  0.75(-0.08, -0.94), and nondominant CWE tendon ICC  =  0.75(-0.12, -0.94).

Discussion

Our data demonstrate acceptable reliability but future studies should include methods that control for variables known to effect STS.

Chinese Ultrasound Doctors Association Guideline on Operational Standards for 2-D Shear Wave Elastography Examination of Musculoskeletal Tissues

The Ultrasound Physician Branch of the Chinese Medical Doctor Association sought to develop evidence-based recommendations on the operational standards for 2-D shear wave elastography examination of musculoskeletal tissues. A consensus panel of 22 Chinese musculoskeletal ultrasound experts reviewed current scientific evidence and proposed a set of 12 recommendations for 13 key issues, including instruments, operating methods, influencing factors and image interpretation. A final consensus was reached through discussion and voting. On the basis of research evidence and expert opinions, the strength of recommendation for each proposition was assessed using a visual analog scale, while further emphasizing the best available evidence during the question-and-answer session. These expert consensus guidelines encourage facilitation of the standardization of clinical practices for collecting and reporting shear wave elastography data.

Using ultrasound elastography to assess non-invasive, non-pharmacological interventions for musculoskeletal stiffness: a systematic review and meta-analysis

PURPOSE

To evaluate the current evidence regarding the use of ultrasound elastography for assessing non-invasive, non-pharmacological interventions for eliciting changes in musculoskeletal stiffness.

METHODS

A systematic search of MEDLINE, CINAHL, EMBASE, and Web of Science databases was performed in accordance with Preferred Reporting Items for Systematic Review and Meta-Analysis (PRISMA) guidelines. Information on measurement and intervention procedures was extracted. Bias was assessed using Cochrane Risk of Bias or Risk of Bias In Non-randomised Studies of Interventions (ROBINS-I) tools for studies with true or quasi-experimental designs, respectively. Analyses were conducted for adequately powered subgroups based on intervention type, measurement site, and population assessed.

RESULTS

Twenty-one studies were included in the review. Overall risk of bias was low for true experimental studies and moderate for quasi-experimental studies. Subgroup analyses indicated a large overall effect for interventions involving manual physiotherapy and taping/splinting for reducing masseter muscle stiffness in patients with masticatory muscle disorders (g = 1.488, 95% CI = 0.320-2.655, p = 0.013). Analyses for other intervention types and patient groups were underpowered.

CONCLUSION

Ultrasound elastography demonstrates clinical applicability for assessing non-invasive, non-pharmacological interventions for musculoskeletal stiffness. However, the comparative efficacy of these interventions for modulating tissue stiffness remains inconclusive.

Effect of durations and pressures of cupping therapy on muscle stiffness of triceps

Cupping therapy has been used for the alleviation of muscle soreness in athletes. However, clinical studies of cupping therapy show conflicting results. Lack of standardized guidelines of the dose-response relationship of cupping therapy, such as appropriate cupping duration and negative pressure, limits the adoption of cupping therapy in clinical practice. The objectives of this study were to investigate the effect of various pressures and durations of cupping therapy on reducing muscle stiffness. The 2 × 2 factorial design with the repeated measures and counterbalanced design was used to test four cupping protocols, including two negative pressures at -225 and -300 mmHg and two durations at 5 and 10 min, in 12 healthy young people. B-mode and elastographic ultrasound was used to assess muscle stiffness of the triceps before and after cupping therapy. The region of interest of elastographic image was divided into the superficial and deep layers for assessing the effect of cupping therapy on stiffness of various depths of the triceps. Normalized stiffness was calculated as a ratio of pre-cupping stiffness divided by post-cupping stiffness of each participant. The two-way analysis of variance (ANOVA) was used to examine the main effects of the pressure and duration factors and the interaction effect between the pressure and duration factors. The results showed that there were no interactions between the pressure and duration factors (overall layer p = 0.149, superficial layer p = 0.632, and deep layer p = 0.491). The main effects of duration of the overall, superficial and deep layers were p = 0.538, p = 0.097 and p = 0.018, respectively. The results showed that 10-min cupping at -300 mmHg is more effective on reducing stiffness of the deep layer of the triceps compared to 5-min cupping (p = 0.031). This study provides the first evidence that the dose of cupping therapy could significantly affect changes of triceps stiffness and the deep layer of the muscle is more sensitive to cupping therapy compared to the superficial and overall layers.

Noninvasive Assessment of In Vivo Passive Skeletal Muscle Mechanics as a Composite Material Using Biomedical Ultrasound

OBJECTIVE

This study develops a biomedical ultrasound imaging method to infer microstructural information (i.e., tissue level) from imaging mechanical behavior of skeletal muscle (i.e., organ level).

METHODS

We first reviewed the constitutive model of skeletal muscle by regarding it as a transversely isotropic (TI) hyperelastic composite material, for which a theoretical formula was established among shear wave speed, deformation, and material parameters (MPs) using the acoustoelasticity theory. The formula was evaluated by finite element (FE) simulations and experimentally examined using ultrasound shear wave imaging (SWI) and strain imaging (SI) on in vivo passive biceps brachii muscles of two healthy volunteers. The imaging sequence included 1) generation of SW in multiple propagation directions while resting the muscle at an elbow angle of 90; 2) generation of SW propagating along the myofiber direction during continuous uniaxial muscle extension by passively changing the elbow angle from 90 to 120. Ultrasound-quantified SW speeds and muscle deformations were fitted by the theoretical formula to estimate MPs of in vivo passive muscle.

RESULTS

Estimated myofiber stiffness, stiffness ratio of myofiber to extracellular matrix (ECM), ECM volume ratio all agreed with literature findings.

CONCLUSION

The proposed mathematical formula together with our in-house ultrasound imaging method enabled assessing microstructural material properties of in vivo passive skeletal muscle from organ-level mechanical behavior in an entirely noninvasive way.

SIGNIFICANCE

Noninvasive assessment of both micro and macro properties of in vivo skeletal muscle will advance our understanding of complex muscle dynamics and facilitate treatment and rehabilitation planning.

Impact of chronic obstructive pulmonary disease on passive viscoelastic components of the musculoarticular system

Chronic obstructive pulmonary disease (COPD) produces skeletal muscle atrophy and weakness, leading to impairments of exercise performance. The mechanical work needed for movement execution is also provided by the passive tension developed by musculoarticular connective tissue. To verify whether COPD affects this component, the passive viscoelastic properties of the knee joint were evaluated in 11 patients with COPD and in 11 healthy individuals. The levels of stiffness and viscosity were assessed by means of the pendulum test, consisting in a series of passive leg oscillations. In addition, to explore the contribution of passive tension in the mechanical output of a simple motor task, voluntary leg flexion-extension movements were performed. Patients with COPD showed a statistically significant reduction in stiffness and viscosity compared to controls. Voluntary execution of flexion-extension movements revealed that the electromyographic activity of the Rectus Femoris and Biceps Femoris was lower in patients than in controls, and the low viscoelastic tension in the patients conditioned the performance of active movements. These results provide novel insights on the mechanism responsible for the movement impairments associated with COPD.

MRI and muscle imaging for idiopathic inflammatory myopathies

Although idiopathic inflammatory myopathies (IIM) are a heterogeneous group of diseases nearly all patients display muscle inflammation. Originally, muscle biopsy was considered as the gold standard for IIM diagnosis. The development of muscle imaging led to revisiting not only the IIM diagnosis strategy but also the patients' follow-up. Different techniques have been tested or are in development for IIM including positron emission tomography, ultrasound imaging, ultrasound shear wave elastography, though magnetic resonance imaging (MRI) remains the most widely used technique in routine. Whereas guidelines on muscle imaging in myositis are lacking here we reviewed the relevance of muscle imaging for both diagnosis and myositis patients' follow-up. We propose recommendations about when and how to perform MRI on myositis patients, and we describe new techniques that are under development.

Bonding and Thermal-Mechanical Property of Gradient NiCoCrAlY/YSZ Thermal Barrier Coatings with Millimeter Level Thickness

The thermal insulation properties of thermal barrier coatings (TBCs) can be significantly improved with increasing the coating thickness. However, due to the weak bonding of high-thickness TBCs, the low reliability and short lifetime greatly limits their application under some severe operating conditions. In this study, a novel and high-efficiency synchronous dual powder feeding method is used to deposit a series of gradient NiCoCrAlY/YSZ coatings with millimeter level thickness. The tensile bonding strengths and residual stress state of coatings are evaluated in order to explore the effect of thickness on the bonding strength of coatings. The results suggested that, due to some micro-convex structure at the “GC/TC” interface and inside “GC” layer, the bonding strength of 1000-μm-thickness gradient NiCoCrAlY/YSZ TBCs with the 4:6 and 2:8 hybrid ratios is over 44 MPa compared to the common TBCs. The fracture position gradually shifts from NiCoCrAlY bond coat to NiCoCrAlY/YSZ transition zone and finally to the YSZ top coat owing to the different position of residual stress concentrations. After thermal cycling tests, the 1000-μm-thickness gradient coating exhibits a higher thermal cycling life. Some coarse cracks initiate and propagate at the bottom region of TBCs, which is mainly due to thermal expansion mismatch stress that finally results in the failure of the gradient coating between the “BC” layer and the substrate.

Quantitative Muscle Ultrasonography Using 2D Textural Analysis: A Novel Approach to Assess Skeletal Muscle Structure and Quality in Chronic Kidney Disease

Chronic kidney disease (CKD) is characterized by progressive reductions in skeletal muscle function and size. The concept of muscle quality is increasingly being used to assess muscle health, although the best means of assessment remains unidentified. The use of muscle echogenicity is limited by an inability to be compared across devices. Gray level of co-occurrence matrix (GLCM), a form of image texture analysis, may provide a measure of muscle quality, robust to scanner settings. This study aimed to identify GLCM values from skeletal muscle images in CKD and investigate their association with physical performance and strength (a surrogate of muscle function). Transverse images of the rectus femoris muscle were obtained using B-mode 2D ultrasound imaging. Texture analysis (GLCM) was performed using ImageJ. Five different GLCM features were quantified: energy or angular second moment (ASM), entropy, homogeneity, or inverse difference moment (IDM), correlation, and contrast. Physical function and strength were assessed using tests of handgrip strength, sit to stand-60, gait speed, incremental shuttle walk test, and timed up-and-go. Correlation coefficients between GLCM indices were compared to each objective functional measure. A total of 90 CKD patients (age 64.6 (10.9) years, 44% male, eGFR 33.8 (15.7) mL/minutes/1.73 m2) were included. Better muscle function was largely associated with those values suggestive of greater image texture homogeneity (i.e., greater ASM, correlation, and IDM, lower entropy and contrast). Entropy showed the greatest association across all the functional assessments (r = -.177). All GLCM parameters, a form of higher-order texture analysis, were associated with muscle function, although the largest association as seen with image entropy. Image homogeneity likely indicates lower muscle infiltration of fat and fibrosis. Texture analysis may provide a novel indicator of muscle quality that is robust to changes in scanner settings. Further research is needed to substantiate our findings.

Quantitative Ultrasound to Assess Skeletal Muscles in Post Stroke Spasticity

Quantitative ultrasound (QUS) techniques such as pixel intensity, ultrasound strain, and shear wave elastography have made it possible to identify the echogenicity (brightness) and mechanical properties (stiffness) of normal and pathological tissues. These techniques can be utilized as an alternative diagnosis tool to assess post stroke spasticity. Current clinical assessment methods include the Modified Ashworth Scale (MAS) and the Modified Tardieu Scale (MTS), which can result in inconsistencies due to their subjective nature. QUS provides robust approaches to assessing muscle stiffness associated with post stroke spasticity. Computer-aided pixel count quantifies tissue echogenicity in grayscale image. A strain ratio in ultrasound strain imaging compares the stiffness and movement (lengthening or shortening) of a spastic muscle with nonspecific muscle. In addition, shear wave elastography provides the shear wave velocity of an affected muscle that directly associated with the muscle stiffness before and after treatment for spasticity. This article reviews the theory behind these aforementioned concepts and discuss the relations between QUS and skeletal muscles in post stroke spasticity.

Skeletal muscle stiffness as measured by magnetic resonance elastography after chronic spinal cord injury: a cross-sectional pilot study

Skeletal muscle stiffness is altered after spinal cord injury (SCI). Assessing muscle stiffness is essential for rehabilitation and pharmaceutical interventions design after SCI. The study used magnetic resonance elastography to assess the changes in stiffness after chronic SCI compared to matched able-bodied controls and determine its association with muscle size, spasticity, and peak torque in persons with SCI. Previous studies examined the association between muscle stiffness and spasticity, however, we are unaware of other studies that examined the effects of muscle composition on stiffness after SCI. Ten participants (one female) with chronic SCI and eight (one female) matched able-bodied controls participated in this cross-sectional study. Magnetic resonance elastography was utilized to monitor stiffness derived from shear waves propagation. Modified Ashworth scale was used to evaluate spasticity scores in a blinded fashion. Peak isometric and isokinetic torques were measured using a biodex dynamometer. Stiffness values were non-significantly lower (12.5%; P = 0.3) in the SCI group compared to able-bodied controls. Moreover, stiffness was positively related to vastus lateralis whole muscle cross-sectional area (CSA) (r² = 0.64, P < 0.005) and vastus lateralis absolute muscle CSA after accounting for intramuscular fat (r² = 0.78, P < 0.0007). Stiffness was also positively correlated to both isometric (r²= 0.55-0.57, P < 0.05) and isokinetic peak (r²= 0.46-0.48, P < 0.05) torques. Our results suggest that larger clinical trial is warranted to confirm the preliminary findings that muscle stiffness is altered after SCI compared to healthy controls. Stiffness appeared to be influenced by infiltration of intramuscular fat and modestly by the spasticity of the paralyzed muscles. The preliminary data indicated that the relationship between muscle stiffness and peak torque is not altered with changing the frequency of pulses or angular velocities. All study procedures were approved by the Institutional Review Board at the Hunter Holmes McGuire VA Medical Center, USA (IRB #: 02314) on May 3, 2017.

Fascial or Muscle Stretching? A Narrative Review

Stretching exercises are integral part of the rehabilitation and sport. Despite this, the mechanism behind its proposed effect remains ambiguous. It is assumed that flexibility increases, e.g., action on muscle and tendon, respectively, but this is not always present in the stretching protocol of the exercises used. Recently, the fasciae have increased popularity and seems that they can have a role to define the flexibility and the perception of the limitation of the maximal range of motion (ROM). Deep fascia is also considered a key element to transmit load in parallel bypassing the joints, transmitting around 30% of the force generated during a muscular contraction. So, it seems impossible dividing the action of the muscles from the fasciae, but they have to be considered as a “myofascial unit”. The purpose of this manuscript is to evaluate the mechanical behavior of muscles, tendons, and fasciae to better understand how they can interact during passive stretching. Stress-strain values of muscle, tendon and fascia demonstrate that during passive stretching, the fascia is the first tissue that limit the elongation, suggesting that fascial tissue is probably the major target of static stretching. A better understanding of myofascial force transmission, and the study of the biomechanical behavior of fasciae, with also the thixotropic effect, can help to design a correct plan of stretching.

Biceps brachii muscle hardness assessed by a push-in meter in comparison to ultrasound strain elastography

This study investigated the relationship between push-in meter (PM) and ultrasound strain elastography (USE) for biceps brachii (BB) muscle hardness. BB hardness of 21 young men was assessed by PM and USE during rest and isometric contractions of six different intensities (15, 30, 45, 60, 75, 90% of maximal voluntary contraction: MVC) at 30°, 60° and 90° elbow flexion. Muscle hardness (E) was calculated from the force-displacement relationship in PM, and strain ratio (SR) between an acoustic coupler (elastic modulus: 22.6 kPa) and different regions of interest (ROIs) in BB was calculated and converted to Young's modulus (YM) in USE. In resting muscle, E was 26.1 ± 6.4 kPa, and SR and YM for the whole BB was 0.88 ± 0.4 and 30.8 ± 12.8 kPa, respectively. A significant (p < 0.01) correlation was evident between E and logarithmical transformed SR (LTSR) for the ROI of whole BB (r = - 0.626), and E and converted YM (r = 0.615). E increased approximately ninefold from resting to 90% MVC, and E and LTSR (r = - 0.732 to - 0.880), and E and converted YM for the SR above 0.1 were correlated (r = 0.599-0.768, p < 0.01). These results suggest that muscle hardness values obtained by PM and USE are comparable.

Finite-Element analysis of a lateral femoro-tibial impact on the total knee arthroplasty

BACKGROUND AND OBJECTIVE

Total knee arthroplasty (TKA) is a routine surgery performed to treat patients with severe knee osteoarthritis. The success of a TKA depends strongly on the initial stability of the prosthetic components and its long-term osseointegration due to the optimal distribution of mechanical stresses in the surrounding bones under the effect of the different biomechanical loads applied to the Femur-TKA-Tibia system. The purpose of this study is to analyze the level and the distribution of the induced stresses in a Femur-TKA-Tibia system subjected to combined triaxial forces, which mimic a femoral mechanical shock.

METHODS

In this study, complex TKA system implanted in both femoral and tibial bones has been analyzed numerically with a three-dimensional finite-element method. A virtual model is designed to examine in silico the effect of the combined triaxial forces acting on this prosthesis in femoral region. Anatomical three-dimensional finite-element models of both femoral and tibial bones were constructed to calculate the interfacial stresses around the TKA components. The 3D finite-element processing program ABAQUS was used to perform the analysis.

RESULTS

The stresses propagated in the bone regions adjacent to the TKA osseointegrated components, and the decreased in their magnitude to the outer region. These stresses reached the highest level in the cortical bone areas that are right next to the proximal upper attachment portions of the TKA osseointegrated components. The magnitude of the stresses in the tibial component is higher than that in the femoral component. Finally, it is very important to emphasize the role of the polyethylene articulating spacer in the shock absorption of bone support sections. Thus, this component should be preserved mechanically from the impact of high shocks in order to maintain healthy TKA systems.

CONCLUSIONS

Optimizing TKA model by controlling the biomechanical stresses distributed within its both components and supporting bones is a valid approach to achieving favorable long-term outcomes. The 3D finite-element analysis provides an effective pre-operative method for planning patient-specific TKA prostheses, and for designing future models that preserves the biomechanical function of the Femur-TKA-Tibia system.

Exploration of New Contrasts, Targets, and MR Imaging and Spectroscopy Techniques for Neuromuscular Disease - A Workshop Report of Working Group 3 of the Biomedicine and Molecular Biosciences COST Action BM1304 MYO-MRI

Neuromuscular diseases are characterized by progressive muscle degeneration and muscle weakness resulting in functional disabilities. While each of these diseases is individually rare, they are common as a group, and a large majority lacks effective treatment with fully market approved drugs. Magnetic resonance imaging and spectroscopy techniques (MRI and MRS) are showing increasing promise as an outcome measure in clinical trials for these diseases. In 2013, the European Union funded the COST (co-operation in science and technology) action BM1304 called MYO-MRI (www.myo-mri.eu), with the overall aim to advance novel MRI and MRS techniques for both diagnosis and quantitative monitoring of neuromuscular diseases through sharing of expertise and data, joint development of protocols, opportunities for young researchers and creation of an online atlas of muscle MRI and MRS. In this report, the topics that were discussed in the framework of working group 3, which had the objective to: Explore new contrasts, new targets and new imaging techniques for NMD are described. The report is written by the scientists who attended the meetings and presented their data. An overview is given on the different contrasts that MRI can generate and their application, clinical needs and desired readouts, and emerging methods.

Magnetic resonance elastography (MRE) shows significant reduction of thigh muscle stiffness in healthy older adults

Determining the effect of ageing on thigh muscle stiffness using magnetic resonance elastography (MRE) and investigate whether fat fraction and muscle cross-sectional area (CSA) are related to stiffness. Six healthy older adults in their eighth and ninth decade and eight healthy young men were recruited and underwent a 3 T MRI protocol including MRE and Dixon fat fraction imaging. Muscle stiffness, fat fraction and muscle CSA were calculated in ROIs corresponding to the four quadriceps muscles (i.e. vastus lateralis (VL), vastus medialis (VM), vastus intermedius (VI), rectus femoris (RF)), combined quadriceps, combined hamstrings and adductors and whole thigh. Muscle stiffness was significantly reduced (p < 0.05) in the older group in all measured ROIs except the VI (p = 0.573) and RF (p = 0.081). Similarly, mean fat fraction was significantly increased (p < 0.05) in the older group over all ROIs with the exception of the VI (p = 0.059) and VL muscle groups (p = 0.142). Muscle CSA was significantly reduced in older participants in the VM (p = 0.003) and the combined quadriceps (p = 0.001), hamstrings and adductors (p = 0.008) and whole thigh (p = 0.003). Over the whole thigh, stiffness was significantly negatively correlated with fat fraction (r = − 0.560, p = 0.037) and positively correlated with CSA (r = 0.749, p = 0.002). Stepwise regression analysis revealed that age was the most significant predictor of muscle stiffness (p = 0.001). These results suggest that muscle stiffness is significantly decreased in healthy older adults. Muscle fat fraction and muscle CSA are also significantly changed in older adults; however, age is the most significant predictor of muscle stiffness.

Theoretical elastic tensile behavior of muscle fiber bundles in traumatic loading events

BACKGROUND

The mechanical characterization of skeletal muscle under high-rate loading regimes is important for predicting traumatic injuries due to traffic accidents and contact sports. However, it is difficult to perform dynamic mechanical tests at rates relevant to such rapid loading events.

METHODS

In the present study, a series of stress relaxation tests were conducted on rabbit hind-limb muscle fiber bundles using a custom tensile tester. Using relatively moderate loading conditions compared to those typically associated with traumatic injuries, the passive stress-decaying mechanical properties of muscle fiber bundles were characterized. In addition, stress relaxation responses to various ramp-hold stretches were theoretically predicted by a custom-built code.

FINDINGS

The results showed that the muscle fiber bundles exhibit greater stress relaxation at higher loading rates and greater stretch magnitudes. Based on these results, the data points representing the "elastic" stress-strain tensile behavior typical of traumatic injury were extrapolated using curve fitting. The theoretical model revealed rate-dependent characteristics of the muscle fiber bundles under traumatic loading conditions, which would result in tensile strengths of 300-500 kPa at the maximum engineering strain of 54%. This strength is on the order of magnitude as the maximum isometric stress of an active muscle contraction.

INTERPRETATION

The proposed numerical model is expected to serve as a powerful research tool to investigate injury mechanisms of the skeletal muscle. Moreover, the elastic response that was theoretically predicted here will be useful in the development of effective countermeasures to prevent traumatic injuries due to rapid loading events.

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

The goal of the study was to investigate the quantitative impact of region of interest (ROI), software choice, muscle fiber orientation and preload tension on shear wave velocity (SWV). First, SWV was assessed in an isotropic elasticity phantom and ex vivo porcine muscle using a commercially available clinical ultrasound system. Secondly, SWV was acquired in relaxed and stretched calf muscles of healthy volunteers (dorsal extension of the talocrural joint), for both parallel and transverse probe direction to the fibers, as well as for different ROIs and software versions. The effect of intermediate probe-fiber alignments was also analyzed. Finally, the impact of confounding factors on SWV reproducibility was minimized with a second force-controlled volunteer study, in which the calf was isometrically loaded, and fiber orientation and ROI were well-defined. 2046 in vivoSWE images were acquired to analyze SWV reproducibility with different confounder settings. In healthy volunteers, the main variance-contributing factors were in order of importance muscle tension, fiber orientation, horizontal ROI size and insertion depth. Regression analysis showed significantly reduced SWV with increasing insertion depth for each study material. Parallel probe-fiber orientation, muscle stretch and increasing horizontal ROI size led to significantly higher SWV. Based on the results of the study, we provide recommendations to minimize the impact of confounders in musculoskeletal elastography and discuss the main confounding mechanisms and trade-offs between confounding variables. Coefficients of variation can be significantly reduced with a controlled protocol, if the confounders are clinically taken into account.

Multiscale modeling of the neuromuscular system: Coupling neurophysiology and skeletal muscle mechanics

Mathematical models and computer simulations have the great potential to substantially increase our understanding of the biophysical behavior of the neuromuscular system. This, however, requires detailed multiscale, and multiphysics models. Once validated, such models allow systematic in silico investigations that are not necessarily feasible within experiments and, therefore, have the ability to provide valuable insights into the complex interrelations within the healthy system and for pathological conditions. Most of the existing models focus on individual parts of the neuromuscular system and do not consider the neuromuscular system as an integrated physiological system. Hence, the aim of this advanced review is to facilitate the prospective development of detailed biophysical models of the entire neuromuscular system. For this purpose, this review is subdivided into three parts. The first part introduces the key anatomical and physiological aspects of the healthy neuromuscular system necessary for modeling the neuromuscular system. The second part provides an overview on state-of-the-art modeling approaches representing all major components of the neuromuscular system on different time and length scales. Within the last part, a specific multiscale neuromuscular system model is introduced. The integrated system model combines existing models of the motor neuron pool, of the sensory system and of a multiscale model describing the mechanical behavior of skeletal muscles. Since many sub-models are based on strictly biophysical modeling approaches, it closely represents the underlying physiological system and thus could be employed as starting point for further improvements and future developments. This article is categorized under: Physiology > Mammalian Physiology in Health and Disease Analytical and Computational Methods > Computational Methods Models of Systems Properties and Processes > Organ, Tissue, and Physiological Models.

Magnetic resonance elastography of skeletal muscle deep tissue injury

The current state-of-the-art diagnosis method for deep tissue injury in muscle, a subcategory of pressure ulcers, is palpation. It is recognized that deep tissue injury is frequently preceded by altered biomechanical properties. A quantitative understanding of the changes in biomechanical properties preceding and during deep tissue injury development is therefore highly desired. In this paper we quantified the spatial-temporal changes in mechanical properties upon damage development and recovery in a rat model of deep tissue injury. Deep tissue injury was induced in nine rats by two hours of sustained deformation of the tibialis anterior muscle. Magnetic resonance elastography (MRE), T₂ -weighted, and T₂ -mapping measurements were performed before, directly after indentation, and at several timepoints during a 14-day follow-up. The results revealed a local hotspot of elevated shear modulus (from 3.30 ± 0.14 kPa before to 4.22 ± 0.90 kPa after) near the center of deformation at Day 0, whereas the T₂ was elevated in a larger area. During recovery there was a clear difference in the time course of the shear modulus and T₂ . Whereas T₂ showed a gradual normalization towards baseline, the shear modulus dropped below baseline from Day 3 up to Day 10 (from 3.29 ± 0.07 kPa before to 2.68 ± 0.23 kPa at Day 10, P < 0.001), followed by a normalization at Day 14. In conclusion, we found an initial increase in shear modulus directly after two hours of damage-inducing deformation, which was followed by decreased shear modulus from Day 3 up to Day 10, and subsequent normalization. The lower shear modulus originates from the moderate to severe degeneration of the muscle. MRE stiffness values were affected in a smaller area as compared with T₂ . Since T₂ elevation is related to edema, distributing along the muscle fibers proximally and distally from the injury, we suggest that MRE is more specific than T₂ for localization of the actual damaged area.

Contemporary image-based methods for measuring passive mechanical properties of skeletal muscles in vivo

Skeletal muscles' primary function in the body is mechanical: to move and stabilize the skeleton. As such, their mechanical behavior is a key aspect of their physiology. Recent developments in medical imaging technology have enabled quantitative studies of passive muscle mechanics, ranging from measurements of intrinsic muscle mechanical properties, such as elasticity and viscosity, to three-dimensional muscle architecture and dynamic muscle deformation and kinematics. In this review we summarize the principles and applications of contemporary imaging methods that have been used to study the passive mechanical behavior of skeletal muscles. Elastography measurements can provide in vivo maps of passive muscle mechanical parameters, and both MRI and ultrasound methods are available (magnetic resonance elastography and ultrasound shear wave elastography, respectively). Both have been shown to differentiate between healthy muscle and muscles affected by a broad range of clinical conditions. Detailed muscle architecture can now be depicted using diffusion tensor imaging, which not only is particularly useful for computational modeling of muscle but also has potential in assessing architectural changes in muscle disorders. More dynamic information about muscle mechanics can be obtained using a range of dynamic MRI methods, which characterize the detailed internal muscle deformations during motion. There are several MRI techniques available (e.g., phase-contrast MRI, displacement-encoded MRI, and "tagged" MRI), each of which can be collected in synchrony with muscle motion and postprocessed to quantify muscle deformation. Together, these modern imaging techniques can characterize muscle motion, deformation, mechanical properties, and architecture, providing complementary insights into skeletal muscle function.

Passive changes in muscle length

This review, the first in a series of minireviews on the passive mechanical properties of skeletal muscles, seeks to summarize what is known about the muscle deformations that allow relaxed muscles to lengthen and shorten. Most obviously, when a muscle lengthens, muscle fascicles elongate, but this is not the only mechanism by which muscles change their length. In pennate muscles, elongation of muscle fascicles is accompanied by changes in pennation and changes in fascicle curvature, both of which may contribute to changes in muscle length. The contributions of these mechanisms to change in muscle length are usually small under passive conditions. In very pennate muscles with long aponeuroses, fascicle shear could contribute substantially to changes in muscle length. Tendons experience moderate axial strains even under passive loads, and, because tendons are often much longer than muscle fibers, even moderate tendon strains may contribute substantially to changes in muscle length. Data obtained with new imaging techniques suggest that muscle fascicle and aponeurosis strains are highly nonuniform, but this is yet to be confirmed. The development, validation, and interpretation of continuum muscle models informed by rigorous measurements of muscle architecture and material properties should provide further insights into the mechanisms that allow relaxed muscles to lengthen and shorten.

Ultrasound Elastography to Assess Botulinum Toxin A Treatment for Post-stroke Spasticity: A Feasibility Study

The aim of the study was to evaluate the feasibility of using ultrasound elastography to assess the effect of botulinum toxin A (BoNT-A) treatment on post-stroke spasticity of the upper limb. In this prospective study, we performed ultrasound elastography (strain imaging and acoustic radiation force impulse-based shear wave elastography) of the spastic biceps brachii muscle in seven patients (five men and two women, mean age: 45y) who underwent BoNT-A injection treatment for post-stroke spasticity of the upper limb. We measured ultrasound elasticity parameters including axial strain ratio (SR) (SR = muscle strain/reference strain), longitudinal SR and shear wave velocity of the biceps muscle immediately before and 17-30 d (mean: 22 d) after BoNT-A injection. Statistical analyses included a two-tailed paired t-test to examine the difference in ultrasound elasticity parameters of bilateral biceps muscles before and after BoNT-A treatment and a Spearman rank correlation coefficient (r_s) to analyze the correlation of ultrasound elasticity parameters to clinical assessment with the Modified Ashworth Scale (MAS) and Tardieu Scale (TS). The difference in ultrasound elasticity parameters before and after BoNT-A treatment was significant (p < 0.05) in the treated spastic muscle. The correlation of ultrasound elasticity parameters with MAS and the angle of catch and range of motion in TS were also significant (r_s = 0.55-0.95, p < 0.05). Our results suggest that ultrasound elastography is feasible in assessment of the effectiveness of BoNT-A treatment for post-stroke spasticity of the upper limb.

A Systematic Review of Continuum Modeling of Skeletal Muscles: Current Trends, Limitations, and Recommendations

Finite elasticity theory has been commonly used to model skeletal muscle. A very large range of heterogeneous constitutive laws has been proposed. In this review, the most widely used continuum models of skeletal muscles were synthetized and discussed. Trends and limitations of these laws were highlighted to propose new recommendations for future researches. A systematic review process was performed using two reliable search engines as PubMed and ScienceDirect. 40 representative studies (13 passive muscle materials and 27 active muscle materials) were included into this review. Note that exclusion criteria include tendon models, analytical models, 1D geometrical models, supplement papers, and indexed conference papers. Trends of current skeletal muscle modeling relate to 3D accurate muscle representation, parameter identification in passive muscle modeling, and the integration of coupled biophysical phenomena. Parameter identification for active materials, assumed fiber distribution, data assumption, and model validation are current drawbacks. New recommendations deal with the incorporation of multimodal data derived from medical imaging, the integration of more biophysical phenomena, and model reproducibility. Accounting for data uncertainty in skeletal muscle modeling will be also a challenging issue. This review provides, for the first time, a holistic view of current continuum models of skeletal muscles to identify potential gaps of current models according to the physiology of skeletal muscle. This opens new avenues for improving skeletal muscle modeling in the framework of in silico medicine.

In Vivo Viscoelastic Response (VisR) Ultrasound for Characterizing Mechanical Anisotropy in Lower-Limb Skeletal Muscles of Boys with and without Duchenne Muscular Dystrophy

Our group has previously found that in silico, mechanical anisotropy may be interrogated by exciting transversely isotropic materials with geometrically asymmetric acoustic radiation force excitations and then monitoring the associated induced displacements in the region of excitation. We now translate acoustic radiation force-based anisotropy assessment to human muscle in vivo and investigate its clinical relevance to monitoring muscle degeneration in Duchenne muscular dystrophy (DMD). Clinical anisotropy assessments were performed using Viscoelastic Response ultrasound, with a degree of anisotropy reflected by the ratios of Viscoelastic Response relative elasticity (RE) or relative viscosity (RV) measured with the asymmetric radiation force oriented parallel versus perpendicular to muscle fiber alignment. In vivo results from rectus femoris and gastrocnemius muscles of boys aged ∼7.9-10.4 y indicate that RE and RV anisotropy ratios in rectus femoris muscles of boys with DMD were significantly higher than those of healthy control boys (RE: DMD = 1.51 ± 0.87, control = 0.99 ± 0.69, p = 0.04, Wilcoxon rank sum test; RV: DMD = 1.04 ± 0.71, control = 0.74 ± 0.22, p = 0.02). In the gastrocnemius muscle, only the RV anisotropy ratio was significantly higher in dystrophic than control patients (DMD = 1.23 ± 0.35, control = 0.88 ± 0.31, p = 0.04). In the dystrophic rectus femoris muscle, the RE anisotropy ratio was inversely correlated (slope = -0.03/lbf, r = -0.43, p = 0.07, Pearson correlation) with quantitative muscle testing functional output measures but was not correlated with quantitative muscle testing in the dystrophic gastrocnemius. These results suggest that Viscoelastic Response RE and RV measures reflect differences in mechanical anisotropy associated with functional impairment with dystrophic degeneration that are relevant to monitoring DMD clinically.

Soft tissue rheology and its implications for elastography: Challenges and opportunities

Magnetic resonance elastography and related shear wave ultrasound elastography techniques can be used to estimate the mechanical properties of soft tissues in vivo by using the relationships between wave propagation and the elastic properties of materials. These techniques have found numerous clinical and research applications, tracking changes in tissue properties as a result of disease or other interventions. Most dynamic elastography approaches estimate tissue elastic (or viscoelastic) properties from a simplified version of the equations for the propagation of acoustic waves through a homogeneous linear (visco)elastic medium. However, soft tissue rheology is complex and departs significantly from this idealized picture. In particular, soft tissues are nonlinearly viscoelastic, inhomogeneous and often anisotropic, and their apparent stiffness can vary with the current loading state. All of these features have implications for the reliability and reproducibility of elastography measurements, from data acquisition to analysis and interpretation. New developments in inversion algorithms for elastography are beginning to offer solutions to account for the complex rheology of tissues, including inhomogeneity and anisotropy. There remains considerable potential to further refine elastography to capture the full spectrum of tissue rheology, and thus to better understand the underlying tissue microstructural changes in a broad range of clinical disorders.

Measurement of large strain properties in calf muscles in vivo using magnetic resonance elastography and spatial modulation of magnetization

It is important to measure the large deformation properties of skeletal muscle in vivo in order to understand and model movement and the force-producing capabilities of muscle. As muscle properties are non-linear, an understanding of how the deformation state affects the measured shear moduli is also useful for clinical applications of magnetic resonance elastography (MRE) to muscle disorders. MRE has so far only been used to measure the linear viscoelastic (small strain) properties of muscles. This study aims to measure the shear moduli of human calf muscles under varying degrees of strain using MRE. Nine healthy adults (four males; age range, 25-38 years) were recruited, and the storage modulus G' was measured at three ankle angle positions: P0 (neutral), P15 (15° plantarflexed) and P30 (30° plantarflexed). Spatial modulation of magnetization (SPAMM) was used to measure the strain in the calf associated with the ankle rotations between P0 to P15 and P0 to P30. SPAMM results showed that, with plantarflexion, there was a shortening of the medial gastrocnemius and soleus muscles, which resulted in an expansion of both muscles in the transverse direction. Strains for each ankle rotation were in the range 3-9% (in compression). MRE results showed that this shortening during plantarflexion resulted in a mean decrease in G' in the medial gastrocnemius (p = 0.013, linear mixed model), but not in the soleus (p = 0.47). This study showed that MRE is a viable technique for the measurement of large strain deformation properties in vivo in soft tissues by inducing physiological strain within the muscle during imaging.

Quantitative Ultrasound Imaging to Assess the Biceps Brachii Muscle in Chronic Post-Stroke Spasticity: Preliminary Observation

We prospectively investigated the feasibility of using quantitative ultrasound imaging (QUI) to assess the biceps brachii muscle (BBM) in individuals with chronic post-stroke spasticity. To quantify muscle echogenicity and stiffness, we measured QUI parameters (gray-scale pixel value and shear wave velocity [SWV, m/s]) of the BBM in three groups: 16 healthy BBMs; 12 post-stroke, non-spastic BBMs; and 12 post-stroke, spastic BBMs. The QUI results were compared with the Modified Ashworth Scale and Tardieu Scale. A total of 20 SWVs were measured in each BBM, once at elbow in 90° flexion and again at maximally achievable extension using acoustic radiation force impulse imaging. BBM pixel value was measured in gray-scale images captured at 90° elbow flexion using ImageJ software. Statistical analyses included analysis of variance for examining the difference in SWV and pixel values among the three groups; Bonferroni correction for testing the difference in SWV and pixel values in a paired group; t-test for examining the difference in SWV values measured at two elbow angles; and Pearson correlation coefficient for analyzing the correlation of QUI to Modified Ashworth Scale and Tardieu Scale. SWV significantly differed between spastic BBMs and non-spastic or healthy BBMs. For pixel values, each of the three groups significantly differed from the others at elbow 90° flexion. The difference in SWV measured between the two elbow angles was also significant (p <0.01). A strong negative correlation was found between SWV and passive range of motion (R² = -0.88, p <0.0001) in spastic upper limbs. These results suggest that the use of QUI is feasible in quantitative assessment of spastic BBM.

Ultrasound Strain Imaging to Assess the Biceps Brachii Muscle in Chronic Poststroke Spasticity

OBJECTIVES

The aim of the study was to assess the feasibility of ultrasound strain imaging in characterizing the biceps brachii muscle in chronic poststroke spasticity.

METHODS

We prospectively analyzed strain imaging data from bilateral biceps brachii muscles in 8 healthy volunteers and 7 patients with poststroke chronic spasticity. Axial deformations of the biceps brachii muscle and overlying subcutaneous tissue were produced by external compression using a sandbag (1.0 kg) attached to a transducer. The lengthening and shortening of the biceps brachii muscle and subcutaneous tissue were produced by manual passive elbow extension (from 90° to 0°) and flexion (from 0° to 90°), respectively. We used offline 2-dimensional speckle tracking to estimate axial and longitudinal strain ratios (biceps brachii strain/subcutaneous tissue strain), and the longitudinal tissue velocity of the biceps brachii muscle. Statistical analyses included analysis of variance for testing differences in strain imaging parameters among healthy, nonspastic, and spastic biceps brachii muscles, the Bonferroni correction for further testing differences in US strain imaging among paired groups (healthy versus spastic, nonspastic versus spastic, and healthy versus nonspastic), and the Pearson correlation coefficient for assessing the intraobserver reliability of performing strain imaging in stroke survivors.

RESULTS

The differences in strain imaging parameters between healthy and spastic and between nonspastic and spastic biceps brachii muscles were significant at both 90° elbow flexion and maximal elbow extension (P < .01). There was no significant difference in axial strain ratios at 90° of elbow flexion or longitudinal tissue velocities between healthy and nonspastic muscles (P > .05). The intraobserver reliability of performing strain imaging in stroke survivors was good (r = 0.85; P < .01).

CONCLUSIONS

Ultrasound strain imaging seems to be feasible for characterizing the biceps brachii muscle in chronic poststroke spasticity.

Muscle Shear Wave Elastography in Inclusion Body Myositis: Feasibility, Reliability and Relationships with Muscle Impairments

Degenerative muscle changes may be associated with changes in muscle mechanical properties. Shear wave elastography (SWE) allows direct quantification of muscle shear modulus (MSM). The aim of this study was to evaluate the feasibility and reliability of SWE in the severely disordered muscle as observed in inclusion body myositis. To explore the clinical relevance of SWE, potential relationships between MSM values and level muscle impairments (weakness and ultrasound-derived muscle thickness and echo intensity) were investigated. SWE was performed in the biceps brachii at 100°, 90°, 70° and 10° elbow flexion in 34 patients with inclusion body myositis. MSM was assessed before and after five passive stretch-shortening cycles at 4°/s from 70° to 10° elbow angle and after three maximal voluntary contractions to evaluate potential effects of muscle pre-conditioning. Intra-class correlation coefficients and standard errors of measurements were >0.83 and <1.74 kPa and >0.64 and <1.89 kPa for within- and between-day values, respectively. No significant effect of passive loading-unloading and maximal voluntary contractions was found (all p values  >0.18). MSM correlated to predicted muscle strength (all Spearman correlation coefficients (ρ) > 0.36; all p values < 0.05). A significant correlation was found between muscle echo intensity and muscle shear modulus at 70° only (ρ = 0.38, p <0.05). No correlation was found between muscle thickness and MSM (all ρ values > 0.23 and all p values > 0.25, respectively). Within- and between-day reliability of muscle SWE was satisfactory and moderate, respectively. SWE shows promise for assessing changes in mechanical properties of the severely disordered muscle. Further investigations are required to clarify these findings and to refine their clinical value.

Ultrasound Strain Imaging in Assessment of Biceps Muscle Stiffness and Dynamic Motion in Healthy Adults

We prospectively evaluated the feasibility of using ultrasound strain imaging (USI) to assess biceps brachii muscle (BBM) stiffness and dynamic motion in 10 healthy adults. The BBM axial deformation was produced by external compression with a sandbag (1.0 kg) tied onto the transducer. The BBM lateral movement was produced by manual passive elbow flexion and extension. By use of 2-D speckle tracking, captured 5-s real-time ultrasound data of BBM were processed to estimate axial strain, representing muscle stiffness, and lateral strain and tissue velocity, representing muscle dynamic motion. Axial (lateral) strain ratio was defined as BBM strain divided by subcutaneous soft tissue strain. There was no significant difference in lateral strain or tissue velocity between the left and right BBM (lateral strain ratio: 4.69 ± 0.07 vs. 4.51 ± 0.08 for extension, 4.82 ± 0.09 vs. 4.69 ± 0.11 for flexion; tissue velocity: 1.58 ± 0.32 cm/s vs. 1.78 ± 0.85 cm/s for extension, -2.03 ± 0.63 vs. -2.03 ± 0.59 for flexion; all p values > 0.05) or between men and women (lateral strain ratio: 4.52 ± 0.06 vs. 4.67 ± 0.1 for extension, 4.71 ± 0.11 vs. 4.83 ± 0.09 for flexion; tissue velocity, cm/s: 1.76 ± 0.76 vs. 1.66 ± 0.65 for extension, -2.21 ± 0.65 vs. -1.88 ± 0.52 for flexion, all p values > 0.05). The difference in axial stain between men and women was significant (axial strain ratio: 3.09 ± 0.43 vs. 3.52 ± 0.26, p = 0.02). Inter- and intra-observer reliability in performing USI of the BBM was good (all intra-class correlation coefficients [ICCs] >0.75). Our results suggest that USI seems to be feasible for and reproducible in estimating BBM mechanical properties and motion dynamics in healthy adults.

Multiscale modeling methods in biomechanics

More and more frequently, computational biomechanics deals with problems where the portion of physical reality to be modeled spans over such a large range of spatial and temporal dimensions, that it is impossible to represent it as a single space-time continuum. We are forced to consider multiple space-time continua, each representing the phenomenon of interest at a characteristic space-time scale. Multiscale models describe a complex process across multiple scales, and account for how quantities transform as we move from one scale to another. This review offers a set of definitions for this emerging field, and provides a brief summary of the most recent developments on multiscale modeling in biomechanics. Of all possible perspectives, we chose that of the modeling intent, which vastly affect the nature and the structure of each research activity. To the purpose we organized all papers reviewed in three categories: 'causal confirmation,' where multiscale models are used as materializations of the causation theories; 'predictive accuracy,' where multiscale modeling is aimed to improve the predictive accuracy; and 'determination of effect,' where multiscale modeling is used to model how a change at one scale manifests in an effect at another radically different space-time scale. Consistent with how the volume of computational biomechanics research is distributed across application targets, we extensively reviewed papers targeting the musculoskeletal and the cardiovascular systems, and covered only a few exemplary papers targeting other organ systems. The review shows a research subdomain still in its infancy, where causal confirmation papers remain the most common. WIREs Syst Biol Med 2017, 9:e1375. doi: 10.1002/wsbm.1375 For further resources related to this article, please visit the WIREs website.

Biomimetic collagen/elastin meshes for ventral hernia repair in a rat model

Ventral hernia repair remains a major clinical need. Herein, we formulated a type I collagen/elastin crosslinked blend (CollE) for the fabrication of biomimetic meshes for ventral hernia repair. To evaluate the effect of architecture on the performance of the implants, CollE was formulated both as flat sheets (CollE Sheets) and porous scaffolds (CollE Scaffolds). The morphology, hydrophylicity and in vitro degradation were assessed by SEM, water contact angle and differential scanning calorimetry, respectively. The stiffness of the meshes was determined using a constant stretch rate uniaxial tensile test, and compared to that of native tissue. CollE Sheets and Scaffolds were tested in vitro with human bone marrow-derived mesenchymal stem cells (h-BM-MSC), and finally implanted in a rat ventral hernia model. Neovascularization and tissue regeneration within the implants was evaluated at 6weeks, by histology, immunofluorescence, and q-PCR. It was found that CollE Sheets and Scaffolds were not only biomechanically sturdy enough to provide immediate repair of the hernia defect, but also promoted tissue restoration in only 6weeks. In fact, the presence of elastin enhanced the neovascularization in both sheets and scaffolds. Overall, CollE Scaffolds displayed mechanical properties more closely resembling those of native tissue, and induced higher gene expression of the entire marker genes tested, associated with de novo matrix deposition, angiogenesis, adipogenesis and skeletal muscles, compared to CollE Sheets. Altogether, this data suggests that the improved mechanical properties and bioactivity of CollE Sheets and Scaffolds make them valuable candidates for applications of ventral hernia repair.

STATEMENT OF SIGNIFICANCE

Due to the elevated annual number of ventral hernia repair in the US, the lack of successful grafts, the design of innovative biomimetic meshes has become a prime focus in tissue engineering, to promote the repair of the abdominal wall, avoid recurrence. Our meshes (CollE Sheets and Scaffolds) not only showed promising mechanical performance, but also allowed for an efficient neovascularization, resulting in new adipose and muscle tissue formation within the implant, in only 6weeks. In addition, our meshes allowed for the use of the same surgical procedure utilized in clinical practice, with the commercially available grafts. This study represents a significant step in the design of bioactive acellular off-the-shelf biomimetic meshes for ventral hernia repair.

A MRI-Compatible Combined Mechanical Loading and MR Elastography Setup to Study Deformation-Induced Skeletal Muscle Damage in Rats

Deformation of skeletal muscle in the proximity of bony structures may lead to deep tissue injury category of pressure ulcers. Changes in mechanical properties have been proposed as a risk factor in the development of deep tissue injury and may be useful as a diagnostic tool for early detection. MRE allows for the estimation of mechanical properties of soft tissue through analysis of shear wave data. The shear waves originate from vibrations induced by an external actuator placed on the tissue surface. In this study a combined Magnetic Resonance (MR) compatible indentation and MR Elastography (MRE) setup is presented to study mechanical properties associated with deep tissue injury in rats. The proposed setup allows for MRE investigations combined with damage-inducing large strain indentation of the Tibialis Anterior muscle in the rat hind leg inside a small animal MR scanner. An alginate cast allowed proper fixation of the animal leg with anatomical perfect fit, provided boundary condition information for FEA and provided good susceptibility matching. MR Elastography data could be recorded for the Tibialis Anterior muscle prior to, during, and after indentation. A decaying shear wave with an average amplitude of approximately 2 μm propagated in the whole muscle. MRE elastograms representing local tissue shear storage modulus G_d showed significant increased mean values due to damage-inducing indentation (from 4.2 ± 0.1 kPa before to 5.1 ± 0.6 kPa after, p<0.05). The proposed setup enables controlled deformation under MRI-guidance, monitoring of the wound development by MRI, and quantification of tissue mechanical properties by MRE. We expect that improved knowledge of changes in soft tissue mechanical properties due to deep tissue injury, will provide new insights in the etiology of deep tissue injuries, skeletal muscle damage and other related muscle pathologies.

Shear waves elastography for assessment of human Achilles tendon's biomechanical properties: an experimental study

INTRODUCTION

Achilles tendon is the most frequently ruptured tendon, but its optimal treatment is increasingly controversial. The mechanical properties of the healing tendon should be studied further. Shear waves elastography (SWE) measures the shear modulus, which is proven to be correlated to elastic modulus in animal tendons. The aim of our study was to study whether the shear moduli of human cadaveric Achilles tendon, given by SWE, were correlated with the apparent elastic moduli of those tendons given by tensile tests.

MATERIALS AND METHODS

Fourteen cadaveric lower-limbs were studied. An elastographic study of the Achilles tendon (AT) was first done in clinical-like conditions. SWE was performed at three successive levels (0, 3 and 6cm from tendon insertion) with elastographic probe oriented parallel to tendon fibers, blindly, for three standardized ankle positions (25° plantar flexion, neutral position, and maximal dorsal flexion). The mean shear moduli were collected through blind offline data-analysis. Then, AT with triceps were harvested. They were subjected to tensile tests. A continuous SWE of the Achilles tendon was performed simultaneously. The apparent elastic modulus was obtained from the experimental stress-strain curve, and correlation with shear modulus (given by SWE) was studied.

RESULTS

Average shear moduli of harvested AT, given by SWE made an instant before the tensile tests, were significantly correlated with shear moduli of the same AT made at the same level, previously in clinical-like condition (p<0.05), only in neutral position. There was a statistical correlation (p<0.005) and a correlation coefficient R² equal to 0.95±0.05, between shear moduli (SWE) and apparent elastic moduli (tensile tests), for 11 tendons (3 tendons were inoperable due to technical error), before a constant disruption in the correlation curves.

DISCUSSION

We demonstrated a significant correlation between SWE of Achilles tendon performed in clinical-like conditions (in neutral position) and SWE performed in harvested tendon. We also found a correlation between SWE performed on harvested tendon and apparent elastic moduli obtained with tensile tests (for 11 specimens). As a consequence, we can suppose that SWE of AT in clinical-like conditions is related to tensile tests. To our knowledge, the ability of SWE to reliably assess biomechanical properties of a tendon or muscle was, so far, only demonstrated in animal models.

CONCLUSION

SWE can provide biomechanical information of the human AT non-invasively.

Ultrasound shear wave velocity in skeletal muscle: A reproducibility study

PURPOSE

The purpose of the study was threefold: to assess the reliability of shear wave velocities (SWV) measurements in normal skeletal muscles; to evaluate intra- and inter-operator reproducibility of measurements for a specific site of the muscle and for the mean value in the whole muscle.

MATERIALS AND METHODS

Two sets of measurements were performed at three weeks intervals of each other on 16 volunteers by two radiologists on medial gastrocnemius and tibialis anterior muscles. Each muscle was evaluated in 5 different sites, with three measurements for each site in the transverse and longitudinal planes. Reliability of SWV measurements was assessed by means of intraclass correlation coefficient (ICC).

RESULTS

Reliability of the three independent SWV measurements was excellent, slightly better in the longitudinal plane. Inter/intra-operator reproducibility per site was fair to good in the longitudinal plane and poor to fair in the transverse plane. For global values of the whole muscle, ICC showed good agreement in the longitudinal plane and fair agreement in the transverse plane.

CONCLUSION

Quantitative SWV measurements are reliable when performed in rigorous conditions. In conditions that mirror clinical practice, inter/intra-operator reproducibility is moderate, better for longitudinal compared to transverse plane.

Muscular hypertonicity: a suspected contributor to rheumatological manifestations observed in ambulatory practice

OBJECTIVE

The objective of this retrospective study of non-inflammatory rheumatic disease patients was to investigate if the individuals clinically identified with muscular hypertonicity (MHT) had increased clinical manifestations compared with those of age- and gender-matched patients with the same disorders.

MATERIAL AND METHODS

The MHT status was clinically identified in the rheumatologist's myofascial protocol examination as relatively increased passive resistance of relaxed muscle on a slow gentle stretch. Clinical and laboratory data were abstracted on a pre-coded form, including symptom and physical examination features, serum assays, and medications.

RESULTS

The 19 MHT cases complained of greater subjective stiffness (p=0.010) and tiredness (p=0.018) at initial encounters and increased aching pain (p=0.049) and were prescribed more (p=0.003) mild narcotic analgesics than the 19 comparison patients. The cases had higher (p=0.027) serum creatine kinase levels, and patients with diffuse MHT had greater frequency of heavy (30+pack-years) cigarette smoking (p=0.002) than comparison subjects. Narcotic usage was also greater in cases with diffuse involvement.

CONCLUSION

Non-inflammatory rheumatic disease patients with MHT had an overall similar profile as that of comparison patients but had greater musculoskeletal complaints, and those with diffuse involvement had greater narcotic usage. Further research, including quantitative measurements of muscle stiffness, are required to determine whether MHT is a documented entity associated with increased rheumatological manifestations.

Three-parameter shear wave inversion in MR elastography of incompressible transverse isotropic media: Application to in vivo lower leg muscles

PURPOSE

To develop and demonstrate MR elastography (MRE) for the measurement of three independent viscoelastic constants of skeletal muscle according to the theory of linear elasticity of incompressible materials with transverse isotropy (TI).

METHODS

Three-dimensional multifrequency MRE was applied to soleus, gastrocnemius, and tibialis anterior muscles in 10 healthy volunteers. The rotational wave fields were solved for complex-valued viscoelastic parameters μ12, μ13, and E3 corresponding to two shear moduli (within the planes of isotropy and symmetry of TI materials) and Young's modulus (along the principal fiber axis).

RESULTS

Anisotropy was represented by the inequality μ12  < μ13  < 1/3E3 considering storage and loss properties of the soleus and gastrocnemius muscles, whereas storage shear moduli of tibialis were indistinguishable. Storage moduli were: 1.06 ± 0.12, 1.33 ± 0.10, 6.92 ± 0.95 kPa (soleus); 0.90 ± 0.11, 1.30 ± 0.15, 8.22 ± 1.37 kPa (gastrocnemius); 1.26 ± 0.16, 1.27 ± 0.11, 9.29 ± 1.42 kPa (tibialis), for μ12, μ13, and E3, respectively. The muscles were different in their μ12 and E3 values, whereas μ13 was less sensitive to the muscle type. Leg differences were observed in the soleus and gastrocnemius muscles.

CONCLUSION

Recovery of the full elasticity tensor in incompressible TI materials is feasible by three-dimensional inversion of the time-harmonic shear wave equation. The method is potentially useful for the clinical evaluation of skeletal muscle anisotropy.