Acousto-elasticity of transversely isotropic incompressible soft tissues: characterization of skeletal striated muscle

Unravelling anisotropic nonlinear shear elasticity in muscles: Towards a non-invasive assessment of stress in living organisms

Acoustoelasticity theory describes propagation of shear waves in uniaxially stressed medium and allows the retrieval of nonlinear elastic coefficients of tissues. In transverse isotropic medium such as muscles the theory leads to 9 different configurations of propagating shear waves (stress axis vs. fibers axis vs. shear wave polarization axis vs. shear wave propagation axis). In this work we propose to use 4 configurations to quantify these nonlinear parameters ex vivo and in vivo. Ex vivo experiments combining ultrasound shear wave elastography and mechanical testing were conducted on iliopsoas pig muscles to quantify three third-order nonlinear coefficients A, H and K that are possibly linked to the architectural structure of muscles. In vivo experiments were performed with human volunteers on biceps brachii during a stretching exercise on an ergometer. A combination of the third order nonlinear elastic parameters was assessed. The knowledge of this nonlinear elastic parameters paves the way to quantify in vivo the local forces produced by muscle during exercise, contraction or movements.

Rotational 3D shear wave elasticity imaging: Effect of knee flexion on 3D shear wave propagation in in vivo vastus lateralis

Skeletal muscle is a complex tissue, exhibiting not only direction-dependent material properties (commonly modeled as a transversely isotropic material), but also changes in observed material properties due to factors such as contraction and passive stretch. In this work, we evaluated the effect of muscle passive stretch on shear wave propagation along and across the muscle fibers using a rotational 3D shear wave elasticity imaging system and automatic analysis methods. We imaged the vastus lateralis of 10 healthy volunteers, modulating passive stretch by imaging at 8 different knee flexion angles (controlled by a BioDex system). In addition to demonstrating the ability of this acquisition and automatic processing system to estimate muscle shear moduli over a range of values, we evaluated potential higher order biomarkers for muscle health that capture the change in muscle stiffness along and across the fibers with changing knee flexion. The median within-subject variability of these biomarkers is found to be <16%, suggesting promise as a repeatable clinical metric. Additionally, we report an unexpected observation: that shear wave signal amplitude along the fibers increases with increasing flexion and muscle stiffness, which is not predicted by transversely isotropic (TI) material simulations. This observation may point to an additional potential biomarker for muscle health or inform other material modeling choices for muscle.

Load sharing between synergistic muscles characterized by a ligand-binding approach and elastography

The skeletal muscle contraction is determined by cross-bridge formation between the myosin heads and the actin active sites. When the muscle contracts, it shortens, increasing its longitudinal shear elastic modulus ([Formula: see text]). Structurally, skeletal muscle can be considered analogous to the molecular receptors that form receptor-ligand complexes and exhibit specific ligand-binding dynamics. In this context, this work aims to apply elastography and the ligand-binding framework to approach the possible intrinsic mechanisms behind muscle synergism. Based on the short-range stiffness principle and the acoustic-elasticity theory, we define the coefficient [Formula: see text], which is directly related to the fraction saturation of molecular receptors and links the relative longitudinal deformation of the muscle to its [Formula: see text]. We show that such a coefficient can be obtained directly from [Formula: see text] estimates, thus calculating it for the biceps brachii, brachioradialis, and brachialis muscles during isometric elbow flexion torque (τ) ramps. The resulting [Formula: see text] curves were analyzed by conventional characterization methods of receptor-ligand systems to study the dynamical behavior of each muscle. The results showed that, depending on muscle, [Formula: see text] exhibits typical ligand-binding dynamics during joint torque production. Therefore, the above indicates that these different behaviors describe the longitudinal shortening pattern of each muscle during load sharing. As a plausible interpretation, we suggested that this could be related to the binding kinetics of the cross-bridges during their synergistic action as torque increases. Likewise, it shows that elastography could be useful to assess contractile processes at different scales related to the change in the mechanical properties of skeletal muscle.

Decoupling Uniaxial Tensile Prestress and Waveguide Effects From Estimates of the Complex Shear Modulus in a Cylindrical Structure Using Transverse-Polarized Dynamic Elastography

Dynamic elastography, whether based on magnetic resonance, ultrasound, or optical modalities, attempts to reconstruct quantitative maps of the viscoelastic properties of biological tissue, properties altered by disease and injury, by noninvasively measuring mechanical wave motion in the tissue. Most reconstruction strategies that have been developed neglect boundary conditions, including quasi-static tensile or compressive loading resulting in a nonzero prestress. Significant prestress is inherent to the functional role of some biological tissues currently being studied using elastography, such as skeletal and cardiac muscle, arterial walls, and the cornea. In the present article a configuration, inspired by muscle elastography but generalizable to other applications, is analytically and experimentally studied. A hyperelastic polymer phantom cylinder is statically elongated in the axial direction while its response to transverse-polarized vibratory excitation is measured. We examine the interplay between uniaxial prestress and waveguide effects in this muscle-like tissue phantom using computational finite element simulations and magnetic resonance elastography measurements. Finite deformations caused by prestress coupled with waveguide effects lead to results that are predicted by a coordinate transformation approach that has been previously used to simplify reconstruction of anisotropic properties using elastography. Here, the approach estimates material viscoelastic properties that are independent of the nonhomogeneous prestress conditions without requiring advanced knowledge of those stress conditions.

Nonlinear propagation of quasiplanar shear wave beams in soft elastic media with transverse isotropy

Model equations are developed for shear wave propagation in a soft elastic material that include effects of nonlinearity, diffraction, and transverse isotropy. A theory for plane wave propagation by Cormack [J. Acoust. Soc. Am. 150, 2566 (2021)] is extended to include leading order effects of wavefront curvature by assuming that the motion is quasiplanar, which is consistent with other paraxial model equations in nonlinear acoustics. The material is modeled using a general expansion of the strain energy density to fourth order in strain that comprises thirteen terms defining the elastic moduli. Equations of motion for the transverse displacement components are obtained using Hamilton's principle. The coupled equations of motion describe diffraction, anisotropy of the wave speeds, quadratic and cubic plane wave nonlinearity, and quadratic nonlinearity associated with wavefront curvature. Two illustrative special cases are investigated. Spatially varying shear vertical wave motion in the fiber direction excites a quadratic nonlinear interaction unique to transversely isotropic soft solids that results in axial second harmonic motion with longitudinal polarization. Shear horizontal wave motion in the fiber plane reveals effects of anisotropy on third harmonic generation, such as beam steering and dependence of harmonic generation efficiency on the propagation and fiber directions.

Axial stress determines the velocity of shear wave propagation in passive but not active muscles in vivo

Ultrasound shear wave elastography can be used to characterize mechanical properties of unstressed tissue by measuring shear wave velocity (SWV), which increases with increasing tissue stiffness. Measurements of SWV have often been assumed to be directly related to the stiffness of muscle. Some have also used measures of SWV to estimate stress, since muscle stiffness and stress covary during active contractions, but few have considered the direct influence of muscle stress on SWV. Rather, it is often assumed that stress alters the material properties of muscle, and in turn, shear wave propagation. The objective of this study was to determine how well the theoretical dependency of SWV on stress can account for measured changes of SWV in passive and active muscles. Data were collected from six isoflurane-anesthetized cats; three soleus muscles and three medial gastrocnemius muscles. Muscle stress and stiffness were measured directly along with SWV. Measurements were made across a range of passively and actively generated stresses, obtained by varying muscle length and activation, which was controlled by stimulating the sciatic nerve. Our results show that SWV depends primarily on the stress in a passively stretched muscle. In contrast, the SWV in active muscle is higher than would be predicted by considering only stress, presumably due to activation-dependent changes in muscle stiffness. Our results demonstrate that while SWV is sensitive to changes in muscle stress and activation, there is not a unique relationship between SWV and either of these quantities when considered in isolation.NEW & NOTEWORTHY Ultrasound shear wave elastography may be an inexpensive way to measure muscle stress in passive muscle. Here, using a cat model we directly measured shear wave velocity (SWV), muscle stress, and muscle stiffness. Our results show that SWV depends primarily on the stress in a passively stretched muscle. In contrast, the SWV in active muscle is higher than would be predicted by considering only stress, presumably due to activation-dependent changes in muscle stiffness.

Noninvasive measurement of local stress inside soft materials with programmed shear waves

Mechanical stresses across different length scales play a fundamental role in understanding biological systems' functions and engineering soft machines and devices. However, it is challenging to noninvasively probe local mechanical stresses in situ, particularly when the mechanical properties are unknown. We propose an acoustoelastic imaging-based method to infer the local stresses in soft materials by measuring the speeds of shear waves induced by custom-programmed acoustic radiation force. Using an ultrasound transducer to excite and track the shear waves remotely, we demonstrate the application of the method by imaging uniaxial and bending stresses in an isotropic hydrogel and the passive uniaxial stress in a skeletal muscle. These measurements were all done without the knowledge of the constitutive parameters of the materials. The experiments indicate that our method will find broad applications, ranging from health monitoring of soft structures and machines to diagnosing diseases that alter stresses in soft tissues.

Biaxial Tensile Prestress and Waveguide Effects on Estimates of the Complex Shear Modulus Using Optical-Based Dynamic Elastography in Plate-Like Soft Tissue Phantoms

Dynamic elastography attempts to reconstruct quantitative maps of the viscoelastic properties of biological tissue, properties altered by disease and injury, by noninvasively measuring mechanical wave motion in the tissue. Most reconstruction strategies that have been developed neglect boundary conditions, including quasi-static tensile or compressive loading resulting in a nonzero prestress. Significant prestress is inherent to the functional role of some biological tissues, such as skeletal and cardiac muscle, arterial walls, and the cornea. In the present article a novel configuration, inspired by corneal elastography but generalizable to other applications, is studied. A polymer phantom layer is statically elongated via an in-plane biaxial normal stress while the phantom's response to transverse vibratory excitation is measured. We examine the interplay between biaxial prestress and waveguide effects in this plate-like tissue phantom. Finite static deformations caused by prestressing coupled with waveguide effects lead to results that are predicted by a novel coordinate transformation approach previously used to simplify reconstruction of anisotropic properties. Here, the approach estimates material viscoelastic properties independent of the nonzero prestress conditions without requiring advanced knowledge of those stress conditions.

Non-destructive mapping of stress and strain in soft thin films through sound waves

Measuring the in-plane mechanical stress in a taut membrane is challenging, especially if its material parameters are unknown or altered by the stress. Yet being able to measure the stress is of fundamental interest to basic research and practical applications that use soft membranes, from engineering to tissues. Here we present a robust non-destructive technique to measure directly in-situ stress and strain in soft thin films without the need to calibrate material parameters. Our method relies on measuring the speed of elastic waves propagating in the film. Using optical coherence tomography, we verify our method experimentally for a stretched rubber membrane, a piece of cling film (about 10 μm thick), and the leather skin of a traditional Irish frame drum. We find that our stress predictions are highly accurate and anticipate that our technique could be useful in applications ranging from soft matter devices to biomaterial engineering and medical diagnosis.

The combined importance of finite dimensions, anisotropy, and pre-stress in acoustoelastography

Dynamic elastography, whether based on magnetic resonance, ultrasound, or optical modalities, attempts to reconstruct quantitative maps of the viscoelastic properties of biological tissue, properties that are altered by disease and injury, by noninvasively measuring mechanical wave motion in the tissue. Most reconstruction strategies that have been developed neglect boundary conditions, including quasistatic tensile or compressive loading resulting in a nonzero prestress. Significant prestress is inherent to the functional role of some biological tissues currently being studied using elastography, such as skeletal and cardiac muscle, arterial walls, and the cornea. In the present article, we review how prestress alters both bulk mechanical wave motion and wave motion in one- and two-dimensional waveguides. Key findings are linked to studies on skeletal muscle and the human cornea, as one- and two-dimensional waveguide examples. This study highlights the underappreciated combined acoustoelastic and waveguide challenge to elastography. Can elastography truly determine viscoelastic properties of a material when what it is measuring is affected by both these material properties and unknown prestress and other boundary conditions?

Optical coherence elastography for assessing the influence of intraocular pressure on elastic wave dispersion in the cornea

The cornea is a highly specialized organ that relies on its mechanical stiffness to maintain its aspheric geometry and refractive power, and corneal diseases such as keratoconus have been linked to abnormal tissue stiffness and biomechanics. Dynamic optical coherence elastography (OCE) is a clinically promising non-contact and non-destructive imaging technique that can provide measurements of corneal tissue stiffness directly in vivo. The method relies on the concepts of elastography where shear waves are generated and imaged within a tissue to obtain mechanical properties such as tissue stiffness. The accuracy of OCE-based measurements is ultimately dependent on the mathematical theories used to model wave behavior in the tissue of interest. In the cornea, elastic waves propagate as guided wave modes which are highly dispersive and can be mathematically complex to model. While recent groups have developed detailed theories for estimating corneal tissue properties from guided wave behavior, the effects of intraocular pressure (IOP)-induced prestress have not yet been considered. It is known that prestress alone can strongly influence wave behavior, in addition to the associated non-linear changes in tissue properties. This present study shows that failure to account for the effects of prestress may result in overestimations of the corneal shear moduli, particularly at high IOPs. We first examined the potential effects of IOP and IOP-induced prestress using a combination of approximate mathematical theories describing wave behavior in thin plates with observations made from data published in the OCE literature. Through wave dispersion analysis, we deduce that IOP introduces a tensile hoop stress and may also influence an elastic foundational effect that were observable in the low-frequency components of the dispersion curves. These effects were incorporated into recently developed models of wave behavior in nearly incompressible, transversely isotropic (NITI) materials. Fitting of the modified NITI model with ex vivo porcine corneal data demonstrated that incorporation of the effects of IOP resulted in reduced estimates of corneal shear moduli. We believe this demonstrates that overestimation of corneal stiffness occurs if IOP is not taken into consideration. Our work may be helpful in separating inherent corneal stiffness properties that are independent of IOP; changes in these properties and in IOP are distinct, clinically relevant issues that affect the cornea health.

Acoustoelasticity in transversely isotropic soft tissues: Quantification of muscle nonlinear elasticity

Recent developments in the field of elastography aim at developing the quantification of new mechanical properties of tissues, that are complementary to the shear modulus, which is characteristic of the linear elastic properties of a quasi-incompressible medium. In this context, measurement of the elastic nonlinearity of tissues was recently proposed based on acoustoelasticity. Up to now, most of the experimental applications of acoustoelasticity theory using Landau formalism in human tissues have assumed isotropy. However, this strong hypothesis does not hold in all human tissues, such as muscles that are generally considered as transversely isotropic (TI). In this work, after reviewing the constraints imposed by TI symmetry on the linear and nonlinear elastic properties of TI media, the acoustoelasticity theory in TI incompressible media is developed and implemented experimentally on a TI polyvinyl alcohol phantom and on ex vivo muscular tissues. Based on this theory and on the evolutions of the shear wave speed, with respect to uniaxial static stress, the nonlinear elastic parameter A is experimentally quantified. The estimations of A in ex vivo bovine and porcine muscles are on the order of hundreds of kPa. This work paves the way for more thorough muscle mechanical properties characterization as well as for the development of a potential new biomarker.

Non-invasive measurement of the internal pressure of a pressurized biological compartment using Lamb waves

In this study, we propose a mechanical analysis for estimating the internal pressure of a finitely deformed spherical compartment from Lamb wave measurements. The proposed analysis produces a dispersion relation associating Lamb wave speed with pressure using limited material parameters (only a strain stiffening term). The analysis was validated on ultrasound bladder vibrometry (UBV) experiments collected from 9 ex vivo porcine bladders before and after formalin cross-linking. Estimated pressures were compared with pressures measured directly by a pressure transducer. The proposed analysis proved broadly effective at estimating pressure from UBV based Lamb wave without calibration as demonstrated by the observed concordance between estimated and measured pressures (Lins CCC = 0.82 (0.66-0.91). Theoretical limitations and potential refinements to improve the accuracy and generality of the approach are discussed.

Plane nonlinear shear wave propagation in transversely isotropic soft solids

Nonlinear wave equations are obtained for the two plane shear wave modes in a transversely isotropic soft solid. The material is modeled using a general expansion of the strain energy density up to fourth order in strain. Whereas, in an isotropic soft solid, leading order nonlinearity for plane wave propagation appears at cubic order in strain, elastic anisotropy in a transversely isotropic material introduces nonlinear effects at quadratic order, including interaction between the modes of a wave with two displacement components. Expressions for second harmonic generation in an elliptically polarized wave field illustrate the low efficiency of nonlinear interactions between the two displacement components, which results from the disparity between propagation speeds of the two shear wave modes. Coupled wave equations with up to cubic nonlinearity are presented and then approximated to describe linearly polarized waves by neglecting interaction between modes. Evolution equations are obtained for linearly polarized progressive waves, and explicit expressions are given in terms of elastic moduli and propagation direction for the coefficients of leading order nonlinearity. Expressions are presented for up to third harmonic generation from a time-harmonic source.