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

Quantitative Viscoelastic Response (QVisR): Direct Estimation of Viscoelasticity With Neural Networks

We present a machine learning method to directly estimate viscoelastic moduli from displacement time-series profiles generated by viscoelastic response (VisR) ultrasound excitations. VisR uses two colocalized acoustic radiation force (ARF) pushes to approximate tissue viscoelastic creep response and tracks displacements on-axis to measure the material relaxation. A fully connected neural network is trained to learn a nonlinear mapping from VisR displacements, the push focal depth, and the measurement axial depth to the material elastic and viscous moduli. In this work, we assess the validity of quantitative VisR (QVisR) in simulated materials, propose a method of domain adaption to phantom VisR displacements, and show in vivo estimates from a clinically acquired dataset.

Quantitative Evaluation of the Lumbar Multifidus Muscle by Shear Wave Dispersion in Healthy Adults: A Preliminary Study

OBJECTIVES

To investigate the application value of shear wave dispersion (SWD) in healthy adults with the lumbar multifidus muscle (LMM), to determine the range of normal reference values, and to analyze the influences of factors on the parameter.

METHODS

Ninety-five healthy volunteers participated in the study, from whom 2-dimensional, shear wave elastography (SWE), and SWD images of the bilateral LMM were acquired in three positions (prone, standing, and anterior flexion). Subcutaneous fat thickness (SFH), SWE velocity, and SWD slope were measured accordingly for analyses.

RESULTS

The mean SWD slope of the bilateral LMM in the prone position was as follows: left: 14.8 ± 3.1 (m/second)/kHz (female) and 13.0 ± 2.5 (m/second)/kHz (male); right: 14.8 ± 3.7 (m/second)/kHz (female) and 14.2 ± 3.4 (m/second)/kHz (male). In the prone position, there was a weak negative correlation between the bilateral LMM SWD slope of activity level 2 and level 1 (β = -1.5 (2 versus 1, left), -1.9 (2 versus 1, right), all P < .05), and between the left SWD slope of activity level 3 and level 1 (β = -2.3 [3 versus 1, left], P < .05). The correlation between SWE velocity and SWD slope value changed with the position: there was a weak positive correlation in the prone position (r = 0.3 [left], 0.37 [right], both P < .05), and a moderate positive correlation in the standing and anterior flexed positions (r = 0.49-0.74, both P < .001). SFH was moderately negatively correlated with bilateral SWD slope values in the anterior flexion (left: r = -0.4, P = .01; right: r = -0.7, P < .01).

CONCLUSIONS

SWD imaging can be used as an adjunct tool to aid in the assessment of viscosity in LMM. Further, activity level, and position are influencing factors that should be considered in clinical practice.

Matrix Nonlinear Viscoelasticity Regulates Skeletal Myogenesis through MRTF Nuclear Localization and Nuclear Mechanotransduction

Mechanically sensitive tissues (e.g., skeletal muscles) greatly need mechanical stimuli during the development and maturation. The extracellular matrix (ECM) mediates these signals through nonlinear viscoelasticity of collagen networks that are predominant components of the ECM. However, the interactions between cells and ECM form a feedback loop, and it has not yet been possible to determine the degree to which, if any, of the features of matrix nonlinear viscoelasticity affect skeletal muscle development and regeneration. In this study, a nonlinear viscoelastic feature (i.e., strain-enhanced stress relaxation (SESR)) in normal skeletal muscles is observed, which however is almost absent in diseased muscles from Duchenne muscular dystrophy mice. It is recapitulated such SESR feature in vitro and separated the effects of mechanical strain and ECM viscoelasticity on myoblast response by developing a collagen-based hydrogel platform. Both strain and stress relaxation induce myogenic differentiation and myotube formation by C2C12 myoblasts, and myogenesis is more promoted by applying SESR. This promotion can be explained by the effects of SESR on actin polymerization-mediated myocardin related transcription factor (MRTF) nuclear localization and nuclear mechanotransduction. This study represents the first attempt to investigate the SESR phenomenon in skeletal muscles and reveal underlying mechanobiology, which will provide new opportunities for the tissue injury treatments.

Feasibility of Phase Velocity Imaging Using Multi Frequency Oscillation-Shear Wave Elastography

OBJECTIVE

To assess viscoelasticity, a pathologically relevant biomarker, shear wave elastography (SWE) generally uses phase velocity (PV) dispersion relationship generated via pulsed acoustic radiation force (ARF) excitation pulse. In this study, a multi-frequency oscillation (MFO)- excitation pulse with higher weight to higher frequencies is proposed to generate PV images via the generation of motion with energy concentrated at the target frequencies in contrast to the broadband frequency motion generated in pulsed SWE (PSWE).

METHODS

The feasibility of MFO-SWE to generate PV images at 100 to 1000 Hz in steps of 100 Hz was investigated by imaging 6 and 70 kPa inclusions with 6.5 and 10.4 mm diameter and ex vivo bovine liver with and without the presence of an aberration layer and chicken muscle ex vivo, and 4T1 mouse breast tumor, in vivo with comparisons to PSWE.

RESULTS

MFO-SWE-derived CNR was statistically higher than PSWE for 6 kPa (both with and without aberration) and 70 kPa (with aberration) inclusions and derived SNR of the liver was statistically higher than PSWE at higher frequency (600-1000 Hz). Quantitatively, at 600-1000 Hz, MFO-SWE improved CNR of inclusions (without and with) aberration on an average by (8.2 and 156)% and of the tumor by 122%, respectively, and improved SNR of the liver (without and with) aberration by (20.2 and 51.5)% and of chicken muscle by 72%, respectively compared to the PSWE.

CONCLUSIONS AND SIGNIFICANCE

These results indicate the advantages of MFO-SWE to improve PV estimation at higher frequencies which could improve viscoelasticity quantification and feature delineation.

Draft Guidance for Industry Duchenne Muscular Dystrophy, Becker Muscular Dystrophy, and Related Dystrophinopathies - Developing Potential Treatments for the Entire Spectrum of Disease

Background

Duchenne muscular dystrophy (DMD) and related dystrophinopathies are neuromuscular conditions with great unmet medical needs that require the development of effective medical treatments.

Objective

To aid sponsors in clinical development of drugs and therapeutic biological products for treating DMD across the disease spectrum by integrating advancements, patient registries, natural history studies, and more into a comprehensive guidance.

Methods

This guidance emerged from collaboration between the FDA, the Duchenne community, and industry stakeholders. It entailed a structured approach, involving multiple committees and boards. From its inception in 2014, the guidance underwent revisions incorporating insights from gene therapy studies, cardiac function research, and innovative clinical trial designs.

Results

The guidance provides a deeper understanding of DMD and its variants, focusing on patient engagement, diagnostic criteria, natural history, biomarkers, and clinical trials. It underscores patient-focused drug development, the significance of dystrophin as a biomarker, and the pivotal role of magnetic resonance imaging in assessing disease progression. Additionally, the guidance addresses cardiomyopathy's prominence in DMD and the burgeoning field of gene therapy.

Conclusions

The updated guidance offers a comprehensive understanding of DMD, emphasizing patient-centric approaches, innovative trial designs, and the importance of biomarkers. The focus on cardiomyopathy and gene therapy signifies the evolving realm of DMD research. It acts as a crucial roadmap for sponsors, potentially leading to improved treatments for DMD.

Muscle ultrasound in hereditary muscle disease

In this review we summarise the key techniques of muscle ultrasound as they apply to hereditary muscle disease. We review the diagnostic utility of muscle ultrasound including its role in guiding electromyography and muscle biopsy sampling. We summarize the different patterns of sonographic muscle involvement in the major categories of genetic muscle disorders and discuss the limitations of the technique. We hope to encourage others to adopt ultrasound in their care for patients with hereditary muscle diseases.

Parametric Analysis of SV Mode Shear Waves in Transversely Isotropic Materials Using Ultrasonic Rotational 3-D SWEI

Ultrasonic rotational 3-D shear wave elasticity imaging (SWEI) has been used to induce and evaluate multiple shear wave modes, including both the shear horizontal (SH) and shear vertical (SV) modes in in vivo muscle. Observations of both the SH and SV modes allow the muscle to be characterized as an elastic, incompressible, transversely isotropic (ITI) material with three parameters: the longitudinal shear modulus μL , the transverse shear modulus μT , and the tensile anisotropy χE . Measurement of the SV wave is necessary to characterize χE , but the factors that influence SV mode generation and characterization with ultrasonic SWEI are complicated. This work uses Green's function (GF) simulations to perform a parametric analysis to determine the optimal interrogation parameters to facilitate visualization and quantification of SV mode shear waves in muscle. We evaluate the impact of five factors: μL , μT , χE , fiber tilt angle [Formula: see text], and F-number of the push geometry on SV mode speed, amplitude, and rotational distribution. These analyses demonstrate that the following hold: 1) as μL increases, SV waves decrease in amplitude so are more difficult to measure in SWEI imaging; 2) as μT increases, the SV wave speeds increase; 3) as χE increases, the SV waves increase in speed and separate from the SH waves; 4) as fiber tilt angle [Formula: see text] increases, the measurable SV waves remain approximately the same speed, but change in strength and in rotational distribution; and 5) as the push beam geometry changes with F-number, the measurable SV waves remain approximately the same speed, but change in strength and rotational distribution. While specific SV mode speeds depend on the combinations of all parameters considered, measurable SV waves can be generated and characterized across the range of parameters considered. To maximize measurable SV waves separate from the SH waves, it is recommended to use an F/1 push geometry and [Formula: see text].

Shear-Wave Elastography and Viscosity PLUS for the Assessment of Peripheric Muscles in Healthy Subjects: A Pre- and Post-Contraction Study

Viscosity is a novel parameter, recently introduced in the use of elastographic techniques, correlating to shear-wave dispersion. The purpose of this study was to provide normal reference viscosity values for the peripheral muscles in healthy volunteers. This prospective study included 38 subjects who underwent US examinations between November 2021 and January 2022. Measurements were taken on the calf and the deltoid muscles in both pre- and post-contraction states. The age range was 21–29 years, with a median of 26 years. The SWE and ViPLUS values in the deltoid muscles were significantly higher than in the soleus muscles in both pre- and post-contraction sets (p = 0.002). There were statistically significant differences between the pre- and post-contraction values for both the SWE and ViPLUS values in the subgroup analysis. The ICC estimates and the 95% confidence intervals were based on a mean rating (k = 2), an absolute agreement, and a two-way random-effects model, demonstrating excellent agreement between the measurements taken by the two examiners.

Quantitative Estimation of Mechanical Anisotropy Using Acoustic Radiation Force (ARF)-Induced Peak Displacements (PD): In Silico and Experimental Demonstration

Elastic degree of anisotropy (DoA) is a diagnostically relevant biomarker in muscle, kidney, breast, and other organs. Previously, elastic DoA was qualitatively assessed as the ratio of peak displacements (PD) achieved with the long-axis of a spatially asymmetric Acoustic Radiation Force Impulse (ARFI) excitation point spread function (PSF) aligned along versus across the axis of symmetry (AoS) in transversely isotropic materials. However, to better enable longitudinal and cross-sectional analyses, a quantitative measure of elastic DoA is desirable. In this study, qualitative ARFI PD ratios are converted to quantitative DoA, measured as the ratio of longitudinal over transverse shear elastic moduli, using a model empirically derived from Field II and finite element method (FEM) simulations. In silico, the median absolute percent error (MAPE) in ARFI-derived shear moduli ratio (SMR) was 1.75%, and predicted SMRs were robust to variations in transverse shear modulus, Young's moduli ratio, speed of sound, attenuation, density, and ARFI excitation PSF dimension. Further, ARFI-derived SMRs distinguished two materials when the true SMRs of the compared materials differed by as little as 10%. Experimentally, ARFI-derived SMRs linearly correlated with the corresponding ratios measured by Shear Wave Elasticity Imaging (SWEI) in excised pig skeletal muscle ( [Formula: see text], MAPE = 13%) and in pig kidney, in vivo ( [Formula: see text], MAPE = 5.3%). These results demonstrate the feasibility of using the ARFI PD to quantify elastic DoA in biological tissues.

Performance of Passive Muscle Stiffness in Diagnosis and Assessment of Disease Progression in Duchenne Muscular Dystrophy

The aim of this study was to evaluate the performance of passive muscle stiffness in diagnosing and assessing disease progression in Duchenne muscular dystrophy (DMD). Boys with DMD and age-matched controls were recruited. Shear wave elastography (SWE) videos were collected by performing dynamic stretching of the gastrocnemius medius (GM). At ankle angles from plantar flexion (PF) 30° to dorsiflexion (DF) 20°, the shear modulus of the GM was measured for each 10° of ankle movement. Shear modulus at each ankle angle was compared between the DMD and control group. Correlation between passive muscle stiffness and motor function grading was also analyzed. A total of 26 patients with DMD and 20 healthy boys were enrolled. At multiple stretch levels, passive muscle stiffness of the GM was significantly higher in patients with DMD than in those in the control group (all p values <0.05). The shear modulus of GM at an ankle angle of DF 10° had the largest area under the receiver operating characteristic curve in differentiating DMD patients from normal subjects (AUC = 0.902, 95% confidence interval: 0.814-0.990). Motor function grading was a significant determinant of passive muscle stiffness at an ankle angle of DF 10° (B = 21.409, t = 3.372, p = 0.003). Passive muscle stiffness may potentially serve as a useful non-invasive tool to monitor disease progression in DMD patients.

Full Characterization of in vivo Muscle as an Elastic, Incompressible, Transversely Isotropic Material Using Ultrasonic Rotational 3D Shear Wave Elasticity Imaging

Using a 3D rotational shear wave elasticity imaging (SWEI) setup, 3D shear wave data were acquired in the vastus lateralis of a healthy volunteer. The innate tilt between the transducer face and the muscle fibers results in the excitation of multiple shear wave modes, allowing for more complete characterization of muscle as an elastic, incompressible, transversely isotropic (ITI) material. The ability to measure both the shear vertical (SV) and shear horizontal (SH) wave speed allows for measurement of three independent parameters needed for full ITI material characterization: the longitudinal shear modulus μL , the transverse shear modulus μT , and the tensile anisotropy χE . Herein we develop and validate methodology to estimate these parameters and measure them in vivo, with μL = 5.77±1.00 kPa, μT = 1.93±0.41 kPa (giving shear anisotropy χμ = 2.11±0.92 ), and χE = 4.67±1.40 in a relaxed vastus lateralis muscle. We also demonstrate that 3D SWEI can be used to more accurately characterize muscle mechanical properties as compared to 2D SWEI.

Electronic Point Spread Function Rotation Using a Three-Row Transducer for ARFI-Based Elastic Anisotropy Assessment: In Silico and Experimental Demonstration

Degree of anisotropy (DoA) of mechanical properties has been assessed as the ratio of acoustic radiation force impulse (ARFI)-induced peak displacements (PDs) achieved using spatially asymmetric point spread functions (PSFs) that are rotated 90° to each other. Such PSF rotation has been achieved by manually rotating a linear array transducer, but manual rotation is cumbersome and prone to misalignment errors and higher variability in measurements. The purpose of this work is to evaluate the feasibility of electronic PSF rotation using a three-row transducer, which will reduce variability in DoA assessment. A Siemens 9L4, with 3×192 elements, was simulated in Field II to generate spatially asymmetric ARFI PSFs that were electronically rotated 63° from each other. Then, using the finite element method (FEM), PD due to the ARFI excitation PSFs in 42 elastic, incompressible, transversely isotropic (TI) materials with shear moduli ratios of 1.0-6.0 were modeled. Finally, the ratio of PDs achieved using the two rotated PSFs was evaluated to assess elastic DoA. DoA increased with increasing shear moduli ratios and distinguished materials with 17% or greater difference in shear moduli ratios (Wilcoxon, ). Experimentally, the ratio of PDs achieved using ARFI PSF rotated 63° from each other distinguished the biceps femoris muscle from two pigs, which had median shear moduli ratios of 4.25 and 3.15 as assessed by shear wave elasticity imaging (SWEI). These results suggest that ARFI-based DoA assessment can be achieved without manual transducer rotation using a three-row transducer capable of electronically rotating PSFs by 63°. It is expected that electronic PSF rotation will facilitate data acquisitions and improve the reproducibility of elastic anisotropy assessments.

Usability and Pitfalls of Shear-Wave Elastography for Evaluation of Muscle Quality and Its Potential in Assessing Sarcopenia: A Review

Sarcopenia is age-related progressive and generalized loss of skeletal muscle mass and strength. Its prevalence is rising, which poses a burden for society because it increases disability and dependency and therefore raises health care costs. Muscle mass quality, however-an essential part of sarcopenia-is not easily diagnosable yet. Recent interest has risen for ultrasonographic evaluation of muscle. This review introduces muscle elastography as a possible, easy and cheap tool to evaluate qualitative muscle parameters. Basic principles of muscle elastography are described, as well as different elastography techniques and some technical considerations. Furthermore, a proposal for practical guidelines is offered and factors influencing muscle stiffness are highlighted.

In vivo Ultrafast Quantitative Ultrasound and Shear Wave Elastography Imaging on Farm-Raised Duck Livers during Force Feeding

Shear wave elastography (speed and dispersion), local attenuation coefficient slope and homodyned-K parametric imaging were used for liver steatosis grading. These ultrasound biomarkers rely on physical interactions between shear and compression waves with tissues at both macroscopic and microscopic scales. These techniques were applied in a context not yet studied with ultrasound imaging, that is, monitoring steatosis of force-fed duck livers from pre-force-fed to foie gras stages. Each estimated feature presented a statistically significant trend along the feeding process (p values <10^-3). However, whereas a monotonic increase in the shear wave speed was observed along the process, most quantitative ultrasound features exhibited an absolute maximum value halfway through the process. As the liver fat fraction in foie gras is much higher than that seen clinically, we hypothesized that a change in the ultrasound scattering regime is encountered for high-fat fractions, and consequently, care has to be taken when applying ultrasound biomarkers to grading of severe states of steatosis.

Viscoelastic Response Ultrasound Derived Relative Elasticity and Relative Viscosity Reflect True Elasticity and Viscosity: In Silico and Experimental Demonstration

Viscoelastic response (VisR) ultrasound characterizes the viscoelastic properties of tissue by fitting acoustic radiation force (ARF)-induced displacements in the region of ARF excitation to a 1-D mass-spring-damper (MSD) model. Elasticity and viscosity are calculated separately but relative to the applied ARF amplitude. We refer to these parameters as "relative elasticity (RE)" and "relative viscosity (RV)." We herein test the hypothesis that RE and RV linearly correlate to true elasticity and viscosity in tissue. VisR imaging was simulated in 144 homogeneous viscoelastic materials with varying elasticities and viscosities. Derived RE linearly correlated with material elasticity and varied by an average of 2.52% when the material viscosity changed from 0.1 to 1.3 Pa · s. Derived RV linearly correlated with material viscosity but varied by an average of 102.5% when material elasticity changed from 3.33 to 20 kPa. The effect of elasticity on RV measurement was compensated using the slope of the linear relationship between RV and natural frequency ( ωtextn ). After compensation, RV [Formula: see text] (elasticity compensated RV) linearly correlated with material viscosity and varied by less than 1.00% on average when the modeled shear elastic modulus changed from 3.3 to 20 kPa. In addition to elasticity compensation, variation in ARF amplitude over depth was compensated, yielding REDC and [Formula: see text]. REDC and [Formula: see text] successfully contrasted elastic and viscous inclusions, respectively, in three simulated phantoms. Experimentally, in the homogeneous oil-in-gelatin phantoms and excised livers, REDC linearly correlated with shear wave dispersion ultrasound vibrometry (SDUV) derived shear elastic modulus, and [Formula: see text] linearly correlated with SDUV-derived shear viscosity. In excised livers containing viscoelastic oil-in-gelatin inclusions, the inclusions were successfully contrasted from the liver background by both REDC and [Formula: see text]. These results suggest that RE and RV are relevant for qualitatively assessing the elastic and viscous properties of tissue.

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.