Application of the novel estimation method by shear wave elastography using vibrator to human skeletal muscle

In recent years, non-invasive measurement of tissue stiffness (hardness) using ultrasound elastography has attracted considerable attention. It has been used to evaluate muscle stiffness in the fields of rehabilitation, sports, and orthopedics. However, ultrasonic diagnostic devices with elastography systems are expensive and clinical use of such devices has been limited. In this study, we proposed a novel estimation method for vibration-based shear wave elastography measurement of human skeletal muscle, then determined its reproducibility and reliability. The coefficient of variation and correlation coefficient were used to determine reproducibility and reliability of the method by measuring the shear wave velocities in konjac phantom gels and agar phantom gels, as well as skeletal muscle. The intra-day, day-to-day, and inter-operator reliabilities were good when measuring the shear wave velocities in phantom gels. The intra-day and day-to-day reliabilities were good when measuring the shear wave velocities in skeletal muscle. The findings confirmed adequate reproducibility and reliability of the novel estimation method for vibration-based shear wave elastography. Therefore, the proposed measurement method may be a useful tool for evaluation of muscle stiffness.

Scientific Reports | (2020) 10:22248 | https://doi.org/10.1038/s41598-020-79215-z www.nature.com/scientificreports/ which comprises a local mechanical property of the constituent material. Hence, quantification of elastic modulus distribution could improve the accuracy of muscle stiffness measurement. A method known as ultrasound strain elastography has been used to assess skeletal muscle stiffness. In this method, an ultrasonic diagnostic probe is held to the body surface and minimal repetitive pressure is applied to tissues. This pressure generates real-time strain images using a transducer. Accordingly, when a consistent amount of pressure is applied, softer tissues demonstrate a greater amount of deformation and therefore experience more strain, compared with stiffer tissues. However, tissue stiffness is measured as a ratio of the amount of actual strain to target strain, indicating that quantitative measurement using Young's modulus remains impossible, although semiquantitative measurement is possible. At least three cycles of compression and decompression are recommended for optimal measurement 4 ; however, extensive tissue preload through repeated compressions can alter tissue elasticity. Notably, this measurement procedure requires specific technical skills, such that the measurement accuracy depends on the examiner's experience [5][6][7][8] .
Recently, an approach that utilizes shear wave elastography has attracted attention as a new method for quantitative estimation of tissue stiffness. Shear wave propagation speed is highly dependent on tissue stiffness; Doppler ultrasonography can be used to measure tissue stiffness because it can measure small propagationrelated displacement fluctuations with relatively high accuracy 9,10 . In ultrasound shear wave elastography, shear waves can be launched in a controlled manner by external mechanical vibrations [11][12][13][14][15][16][17][18][19][20][21][22] and acoustic radiation force impulses [23][24][25][26][27][28] . The use of an external vibrating source is much less common among commercial ultrasound shear wave elastography systems because it is relatively difficult to image the target tissue at appropriate depth; moreover, the required devices can be cumbersome 11,22 . Most commercial ultrasound shear wave elastography systems (e.g., Philips, GE Healthcare, Siemens, Canon, Mindray, Samsung, and Supersonic) use acoustic radiation force impulse-based techniques to generate ultrasonic shear waves in tissues 27,28 . This method requires strong ultrasonic waves, but undesirable temperature elevation can occur if bone and surrounding tissues are located near the focal point 29 . One major limitation of this method is that the detection of shear waves generally requires ultrasound imaging equipment with a high frame rate; these systems are expensive. Therefore, it is difficult to introduce such systems in general clinical practice, rehabilitation facilities, and sports medicine clinics.
Accordingly, external mechanical vibration has been proposed as a cost-effective alternative method for shear-wave induction, with the goal of enabling wider adoption of shear wave elastography 14,17,20,22 . Recently, there has been considerable progress in developing the technological foundation required for external mechanical vibration-based shear wave elastography imaging 16,19,30 . Vibro-elastography has found broad clinical applications in prostate cancer detection and liver fibrosis staging 20,31,32 . To the best of our knowledge, there is no widely available method that can visualize the shear wave wavefront in real time by means of conventional external vibration-based measurement. A shear wave that is reflected from the tissue boundary produces a standing wave in the vicinity of the boundary, which reduces the accuracy of shear wave velocity estimation. The complicated propagation of the shear wave inside tissue also diminishes measurement accuracy. Therefore, we presume that confirmation of the propagation wavefront (i.e., propagation direction, reflection, refraction, and complicated wavefront) of the shear wave and measurement of the shear wave velocity will facilitate quantitative measurements with good reproducibility.
Color Doppler shear wave imaging is a novel wave front imaging method in which ultrasound color flow images are used to detect the wave fronts of shear waves; this approach allows direct, real-time observation of the wave fronts of shear waves [33][34][35] . Shear wave elastography uses continuous shear waves generated by a mechanical vibrator attached to the tissue surface. The addition of a very simple device system to a traditional ultrasound device allows quantitative visualization of tissue stiffness; this offers an easy-to-use method for evaluation of the viscoelastic properties of superficial tissues. However, the adaptation of this system to skeletal muscle has not yet been investigated in detail. Therefore, this study was performed to examine the reproducibility and reliability of a novel estimation method for vibration-based shear wave elastography measurement of human skeletal muscle.
Evaluation of reproducibility when applying the novel estimation method to vibration-based shear wave elastography measurement of skeletal muscle. Figure 3 shows the CV and ICC values of the intra-day (CV, < 9.8; ICC, > 0.92) and day-to-day (CV, < 9.6; ICC, > 0.90) reliabilities of shear wave velocity estimation of skeletal muscle stiffness. It also shows example color flow images and shear wave velocity maps for the biceps brachii, flexor carpi radialis, semitendinosus, biceps femoris, medial gastrocnemius, and tibialis anterior.

Discussion
Our study provides evidence of excellent reliability when using our novel estimation method for vibration-based shear wave elastography measurements of phantom elasticity and skeletal muscle stiffness. These observations are clinically important because skeletal muscle stiffness is observed in various conditions, such as myofascial pain syndrome 36 , neck and shoulder pain 27 , and low back pain 37 . Thus, our results suggest that multiple patient groups who experience muscle stiffness may benefit from our proposed method, which is a useful tool for evaluating diseases characterized by muscle stiffness and determining treatment effects; this method is expected to help reduce costs and promote introduction of vibration-based shear wave elastography in general clinical practice. The proposed estimation method for vibration-based shear wave elastography demonstrated excellent intraday, day-to-day, and inter-operator reliabilities when measuring phantom elasticity. Konjac gel is used as a phantom for ultrasound imaging 38,39 and elastography 34 . The mean shear wave velocity in konjac phantom gel was 2.2 m/s, which is reasonably consistent with the findings in a previous study 34 . The CV calculated from the mean and standard deviation of five measurements during one trial was 0.38%-1.38%, indicating excellent   Table 1. Intra-and inter-operator reliabilities of agar phantom gel elasticity. The intra-class correlation coefficient (ICC), 95% confidence interval (CI), and coefficient of variation (CV) were calculated for each trial. Phantom gels with different concentrations of agar (1%, 2%, 3%, 4%, and 8%) were measured twice per day for 2 days with a between-measurement interval of 1 day (three measurements per trial). www.nature.com/scientificreports/ repeatability for measurements of phantom elasticity. Therefore, triplicate measurements during each trial were performed in subsequent experiments. Intra-day and day-to-day reliabilities were measured to evaluate temporal variation. These results revealed that our proposed method was able to reliably measure elasticity at two different time points (intra-day CV, 0.2%-0.6%; day-to-day CV, 0.6%-0.7%). Furthermore, inter-operator reliability was measured to evaluate operator variation; the findings showed that two independent operators produced reliable elasticity measurements using our novel method (inter-operator CV, 0.5%-0.7%).
To verify the effectiveness of our proposed method, agar gels were used because agar is economical and because elasticity can be easily manipulated by modifying the agar concentration 34,40 . Several attempts were made to produce different concentrations of agar phantom gels, which were prepared based on the recipes described in Table 2. These recipes ensured that the respective components did not separate, powdered agar was dissolved in the solution, and the agar was uniformly solidified. Our results revealed that, as the concentration of agar increased, shear wave velocity increased; moreover, the agar concentration was positively correlated with shear wave velocity. In addition, the CV and ICC measurements were consistent among different agar concentrations, www.nature.com/scientificreports/ revealing good intra-day, day-to-day, and inter-operator reliabilities. Therefore, our novel estimation method for shear wave elastography may be reliable for measurement of elasticity.
The most striking finding of this study was that our estimation method for shear wave elastography had excellent intra-day and day-to-day reliabilities when measuring the shear wave velocities in the right biceps brachii, flexor carpi radialis, semitendinosus, biceps femoris, medial gastrocnemius, and tibialis anterior. This finding suggests that our method can reliably assess muscle elasticity.
Other authors have proposed alternative techniques to measure muscle stiffness in healthy skeletal muscle by means of ultrasound elastography 27,[41][42][43][44][45] . It is difficult to compare shear wave velocities measured in our study with values measured in previous studies because of inconsistencies in probe positioning, variations in ultrasound elastography techniques, and differences in the methods used to report stiffness. Available shear wave elastography measurement parameters include shear wave speed, elastic modulus, shear elastic modulus, Young's modulus, and elasticity [44][45][46][47] . Notably, Ballyns et al. 41 measured biceps brachii shear wave velocity by means of external vibrations. Although their probe and vibrator positions, frequencies, and analytical methods differed from our method, the overall approaches were similar in their study and our study; thus, we presumed that their results could be readily compared with our findings. In the study by Ballyns et al., measurements were performed at 60, 110, 160, and 200 Hz; measurements at 276 Hz (the frequency used in our study) were not directly reported. However, the data provided in their publication implied that biceps brachii shear wave velocity at 276 Hz would be approximately 4 m/s, similar to our measured values of 3.02-5.12 m/s (see Supplementary  Fig. S1 online). Therefore, the shear wave velocity measured by the proposed method is similar to the value reported in a previous study that involved external vibration measurements.
Our findings are comparable with those of previous studies that have shown good/excellent reliability in measuring muscle elasticity by means of shear wave elastography in resting muscles and during active contraction [48][49][50] . Lacourpaille et al. 48 investigated the reliability of elastic modulus in resting muscles by means of shear wave elastography using supersonic shear wave imaging. They reported good intra-session (n = 20) and inter-day (n = 21) reliabilities (biceps brachii ICC, 0.868 and 0.832, respectively; medial gastrocnemius ICC, 0.950 and 0.922, respectively; and tibialis anterior ICC, 0.811 and 0.809, respectively). Depending on the muscle, the mean temporal variability in the CV over 10 measurements ranged from 2.9% to 5.0%. Dubois et al. 49 examined the intra-operator repeatabilities of a measurement protocol to assess shear modulus in the resting lower limb muscles of 10 healthy participants who underwent shear wave elastography using supersonic shear wave imaging. They reported poor repeatability in the biceps femoris (CV, 18%), semimembranosus (CV, 15%), semitendinosus (CV, 20%), and medial gastrocnemius (CV, 9%) muscles. In that study, the CV of shear modulus was calculated and measured six times. Koo et al. 50 performed shear wave elastography using supersonic shear wave imaging to estimate muscle elasticity in the resting tibialis anterior of 20 healthy participants, then calculated the CV of shear elastic modulus in triplicate. They observed excellent test-retest reliability (ICC , 0.942).
In this study, we confirmed adequate reproducibility and reliability using our novel estimation method for vibration-based shear wave elastography measurement of skeletal muscle. To verify inter-day and day-to-day reliabilities, we measured participants' major surface muscles (i.e., biceps brachii, flexor carpi radialis, semitendinosus, biceps femoris, medial gastrocnemius, and tibialis anterior) twice per day for 2 days, with a betweenmeasurement interval of 1 day. For all muscles, the intra-day (CV, < 9.8; ICC, > 0.92) and day-to-day (CV, < 9.6; ICC, > 0.90) reliabilities of the novel estimation method for skeletal muscle were good. Compared with the results of prior studies 48-50 , our findings indicate that the proposed method is a highly reliable approach for estimating major muscle elasticity.
Good reliability in our study may have been achieved through standardization of the scanning site with reference to body landmarks, imaging parameters, and both location and size of the region of interest. In addition, the positions of the probe and excitation point were adjusted to ensure that shear waves propagated in the direction of muscle fibers. Muscle stiffness at rest is perhaps the simplest measurement to capture using ultrasound elastography, but this method must be performed in a standardized manner. In particular, the orientation of the transducer probe affects measurement reliability 46 . Because our method can directly show shear wave wavefronts in real time, the observer can easily identify the optimal probe orientation and/or excitation point, then reproduce shear wave wavefronts that propagate uniformly in the region of interest; these approaches facilitate good reliability when measuring muscle stiffness.
Limitations. There were several notable limitations in the present study. First, we did not investigate correlations with other quantitative methods in this study, because we previously determined the effectiveness of the new estimation method by comparison with values measured using other instruments. Moreover, we previously Table 2. Recipes for agar phantom gel preparation. Phantom gels with different concentrations of agar (1%, 2%, 3%, 4%, and 8%) were mixed with sterilized milk.

Concentration (%)
Powdered agar (g) Sterilized milk (ml) www.nature.com/scientificreports/ showed that shear wave velocity measurements in phantom gels, supraclavicular muscles, and trapezius muscles acquired using the new method were consistent with those acquired by means of other instruments 15,33,34,51,52 .
Comparison of shear wave velocity among skeletal muscle sites and quantitative methods is an important consideration for clinical application of novel skeletal muscle stiffness measurement methods and devices. Second, all participants were asymptomatic young male volunteers and elasticity was measured in static muscles. Thus, the findings may not be generalizable to a variety of patients or muscles. Nevertheless, the findings provide important support for further validation of this technique's reliability under different conditions. In particular, this technique should be performed in asymptomatic patients and the findings should be compared with measurements in patients who exhibit various symptoms or pathologies.
Third, only the elasticity of the middle portion of the muscle was examined, on the basis of the recommended placement of surface electrodes for standardized electromyography measurement. Although we focused on sites that are often measured in clinical practice, our results were highly reproducible and reliable in all parts of skeletal muscle measured in this study. Moreover, the results of this study provide important guidance for clinical application of our method. In the future, we plan to investigate areas other than the central section of the muscle, as well as muscle-tendon transition zones, to determine differences in shear elastic modulus among these sites.
Finally, only the elasticity of major surface muscles was examined, largely because these muscles are frequently measured in clinical practice. As noted above, our results were highly reproducible and reliable in the measured muscles. In the future, we plan to investigate other muscles (i.e., deep tissues) to support clinical implementation of our method. As shown in the photograph of the hamstring (Fig. 3), tissues can be readily measured close to the body surface; however, the measurement of deep tissues may be difficult. A previous study 51 showed that stable measurement of deeper tissues may be achieved by modifying the shear wave frequency of the actuator with respect to shear wave tissue absorption. However, the feasibility of this approach must be confirmed through additional investigations.

Conclusions
In this experimental study, we assessed phantom gels and volunteers' skeletal muscles by means of vibration-based shear wave elastography. Our findings confirmed adequate reproducibility and reliability of a novel estimation method for shear wave elastography measurement of skeletal muscle at sites commonly assessed in clinical practice. Depending on the results of additional validation studies in study populations with different ethnicities and disease states, the proposed measurement method may constitute a useful tool for evaluation of muscle stiffness in both clinical and research settings.

Methods
Experimental system. Figure 4 illustrates the experimental system. Our novel estimation method for shear wave elastography was implemented using ultrasound scanners (Xario SSA-660A, Toshiba, Japan). The shear wave imaging method was previously described; quantitative evaluation confirmed a good correlation between the small displacement estimation method and shear wave velocity measured using our new method 15,33,34 . Notably, shear wave velocity is derived from stiffness by assuming that shear viscosity is negligible. Thus, comparative experiments between these two methods for human supraspinatus muscle, trapezius muscle, and an agar gel phantom have been performed; shear wave velocities measured using our novel estimation method were consistent with measurements performed by means of acoustic radiation force impulses 34,51,52 . In our novel estimation method, the shear wave frequency is selected from among several frequencies that satisfy the shear wave frequency conditions. The spatial resolution of a shear wave map decreases upon selection of a low frequency. However, determination of shear wave displacement amplitude is difficult because of the limited A short (45 mm) and lightweight (30 g) linear vibration motor, designed to manipulate an electric toothbrush, was adopted for the actuators by generating vibrations with a displacement amplitude of > 400 μm and frequency range of 230-300 Hz. Displacement of the linear vibration motor was measured using a microelectromechanical system accelerometer (ADXL001, Analog Devices, USA) at an accelerometer sampling frequency of 200 kHz; the displacement was derived by two-step integration. Figure 5 shows the vibration displacement when positioned over skeletal muscle. The harmonic distortion was measured by Fourier analysis. Power spectrum analysis showed that the second and third harmonic components of displacement were − 44.5 dB and − 79.9 dB at a vibration frequency of 276 Hz and a displacement amplitude of 300 μm. This result was considered negligible and the harmonic distortion was therefore ignored. Hence, this actuator was presumed to generate a sinusoidal waveform that was suitable for the novel estimation method.
To reconstruct shear wave maps, the color flow image on the ultrasound scanner was recorded as a moving picture for 2 s. Then, the color flow image was fed into a computer by means of a video capture device. Fourier analysis and a directional filter were applied; shear wave phase, velocity, and propagation maps were reconstructed. The central processing unit time required for image reconstruction was < 1 s. In our experimental setup, shear wave maps were reconstructed at intervals of 4 s until the stop button on the computer screen was clicked.
Participants. Fourteen healthy male college students voluntarily participated in the study. Written informed consent was obtained prior to the first measurement and all study procedures related to human participants were approved by the Human Research Ethics Committee of Nihon Fukushi University. All methods were performed in accordance with the relevant guidelines and regulations. Participants were excluded if they: (1) had lower-extremity contracture; (2) had undergone surgery on their back or lower extremities; (3) had neurological disorders; (4) took hormones or muscle-affecting drugs; and/or (5) regularly engaged in competitive sports or resistance, aerobic, or flexibility training. The participants' mean (± standard deviation) age, height, body mass, and body mass index were 21.7 ± 0.9 years, 170.7 ± 6.6 cm, 63.5 ± 11.2 kg, and 21.8 ± 3.3 kg/m 2 , respectively. Participants were asked to maintain their normal dietary habits and refrain from vigorous physical activity the day before and immediately before the experiment. Experimental design. Shear wave velocities were measured in phantom gels (konjac and agar) and skeletal muscles (right biceps brachii, flexor carpi radialis, semitendinosus, biceps femoris, medial gastrocnemius, and tibialis anterior) twice per day for 2 days, with a between-measurement interval of 1 day. During measurement of phantom gels, two investigators measured shear wave velocity. Phantom gels with different concentrations of agar (1%, 2%, 3%, 4%, and 8%) were mixed with sterilized milk ( Table 2). The probe of an ultrasonic diagnostic device was placed directly over the measurement point and the excitation point of the vibrator was positioned next to the probe (Fig. 6). Before measurement, the probe position and excitation point were carefully adjusted, while the wave fronts of shear waves were monitored. In addition, the investigator arranged his posture to ensure maintenance of consistent application force to the vibrator and probe during each measurement. In most instances, the investigator's elbows were placed on the bed to maintain a stable posture. This positioning was maintained and 3-5 measurements were performed. During measurements of skeletal muscle, the probe position and excitation point were parallel to the long axis of the muscle fibers; this positioning was adjusted to ensure that shear waves clearly propagated in the muscle fiber direction. The elasticity of the middle portion of www.nature.com/scientificreports/ the muscle was examined and adjusted with reference to the recommended placement of surface electrodes for standardized electromyography measurement 53,54 . The measurement position was determined based on a body surface landmark. If the positions of the probe and vibrator varied between measurements, the measurement location and shear wave propagation would have differed, resulting in inconsistent measurements. Therefore, the positions of the probe and vibrator were marked on the body surface; the probe and vibrator were positioned at the same site for each measurement. Furthermore, a B-mode image acquired during the first measurement was used as a reference image. Measurements were performed on B-mode images; characteristic points of the fascia, bone, nerves, and blood vessels were kept consistent as much as possible.
Statistical analysis. Statistical analyses were performed using the Statistical Package for the Social Sciences version 15.0 (SPSS Inc., USA). Distribution normality was assessed using the Shapiro-Wilk test. CV values were calculated to indicate measurement reproducibility and reliability 55 . To investigate the relationship between agar gel concentration and shear wave velocity, Pearson's correlation analysis was performed. Data were compared using a paired t-test. Differences between conditions were considered statistically significant when P values were < 0.05. ICC (1,1) and ICC (2,1) values with 95% confidence intervals were calculated to determine intra-day, day-to-day, and inter-operator reliabilities. Guidelines from Landis and Koch 56 were used to interpret reliability values as follows: 0.0-0.2 indicated slight agreement; 0.21-0.40 indicated fair agreement; 0.41-0.60 indicated moderate agreement; 0.61-0.80 indicated substantial agreement; and 0.81-1.0 indicated near-perfect or perfect agreement.

Data availability
All data needed to evaluate the conclusions in the paper are present in the paper, supplementary material, and/ or references cited within.