Nonlinear elastic behavior of phantom materials for elastography

Biomechanics of the Human Osteochondral Unit: A Systematic Review

The damping system ensured by the osteochondral (OC) unit is essential to deploy the forces generated within load-bearing joints during locomotion, allowing furthermore low-friction sliding motion between bone segments. The OC unit is a multi-layer structure including articular cartilage, as well as subchondral and trabecular bone. The interplay between the OC tissues is essential in maintaining the joint functionality; altered loading patterns can trigger biological processes that could lead to degenerative joint diseases like osteoarthritis. Currently, no effective treatments are available to avoid degeneration beyond tissues' recovery capabilities. A thorough comprehension on the mechanical behaviour of the OC unit is essential to (i) soundly elucidate its overall response to intra-articular loads for developing diagnostic tools capable of detecting non-physiological strain levels, (ii) properly evaluate the efficacy of innovative treatments in restoring physiological strain levels, and (iii) optimize regenerative medicine approaches as potential and less-invasive alternatives to arthroplasty when irreversible damage has occurred. Therefore, the leading aim of this review was to provide an overview of the state-of-the-art-up to 2022-about the mechanical behaviour of the OC unit. A systematic search is performed, according to PRISMA standards, by focusing on studies that experimentally assess the human lower-limb joints' OC tissues. A multi-criteria decision-making method is proposed to quantitatively evaluate eligible studies, in order to highlight only the insights retrieved through sound and robust approaches. This review revealed that studies on human lower limbs are focusing on the knee and articular cartilage, while hip and trabecular bone studies are declining, and the ankle and subchondral bone are poorly investigated. Compression and indentation are the most common experimental techniques studying the mechanical behaviour of the OC tissues, with indentation also being able to provide information at the micro- and nanoscales. While a certain comparability among studies was highlighted, none of the identified testing protocols are currently recognised as standard for any of the OC tissues. The fibril-network-reinforced poro-viscoelastic constitutive model has become common for describing the response of the articular cartilage, while the models describing the mechanical behaviour of mineralised tissues are usually simpler (i.e., linear elastic, elasto-plastic). Most advanced studies have tested and modelled multiple tissues of the same OC unit but have done so individually rather than through integrated approaches. Therefore, efforts should be made in simultaneously evaluating the comprehensive response of the OC unit to intra-articular loads and the interplay between the OC tissues. In this regard, a multidisciplinary approach combining complementary techniques, e.g., full-field imaging, mechanical testing, and computational approaches, should be implemented and validated. Furthermore, the next challenge entails transferring this assessment to a non-invasive approach, allowing its application in vivo, in order to increase its diagnostic and prognostic potential.

Characterization of the concentration of agar-based soft tissue mimicking phantoms by impact analysis

In various medical fields, a change of soft tissue stiffness is associated with its physio-pathological evolution. While elastography is extensively employed to assess soft tissue stiffness in vivo, its application requires a complex and expensive technology. The aim of this study is to determine whether an easy-to-use method based on impact analysis can be employed to determine the concentration of agar-based soft tissue mimicking phantoms. Impact analysis was performed on soft tissue mimicking phantoms made of agar gel with a mass concentration ranging from 1% to 5%. An indicator Δt is derived from the temporal variation of the impact force signal between the hammer and a small beam in contact with the sample. The results show a non-linear decrease of Δt as a function of the agar concentration (and thus of the sample stiffness). The value of Δt provides an estimation of the agar concentration with an error of 0.11%. This sensitivity of the impact analysis based method to the agar concentration is of the same order of magnitude than results obtained with elastography techniques. This study opens new paths towards the development of impact analysis for a fast, easy and relatively inexpensive clinical evaluation of soft tissue elastic properties.

A Force-Matched Approach to Large-Strain Nonlinearity in Elasticity Imaging for Breast Lesion Characterization

OBJECTIVE

Ultrasound elasticity imaging is a class of ultrasound techniques with applications that include the detection of malignancy in breast lesions. Although elasticity imaging traditionally assumes linear elasticity, the large strain elastic response of soft tissue is known to be nonlinear. This study evaluates the nonlinear response of breast lesions for the characterization of malignancy using force measurement and force-controlled compression during ultrasound imaging.

METHODS

54 patients were recruited for this study. A custom force-instrumented compression device was used to apply a controlled force during ultrasound imaging. Motion tracking derived strain was averaged over lesion or background ROIs and matched with compression force. The resulting force-matched strain was used for subsequent analysis and curve fitting.

RESULTS

Greater median differences between malignant and benign lesions were observed at higher compressional forces (p-value < 0.05 for compressional forces of 2-6N). Of three candidate functions, a power law function produced the best fit to the force-matched strain. A statistically significant difference in the scaling parameter of the power function between malignant and benign lesions was observed (p-value = 0.025).

CONCLUSIONS

We observed a greater separation in average lesion strain between malignant and benign lesions at large compression forces and demonstrated the characterization of this nonlinear effect using a power law model. Using this model, we were able to differentiate between malignant and benign breast lesions.

SIGNIFICANCE

With further development, the proposed method to utilize the nonlinear elastic response of breast tissue has the potential for improving non-invasive lesion characterization for potential malignancy.

Bubble dynamics and speed of jets for needle-free injections produced by thermocavitation

Significance

The number of injections administered has increased dramatically worldwide due to vaccination campaigns following the COVID-19 pandemic, creating a problem of disposing of syringes and needles. Accidental needle sticks occur among medical and cleaning staff, exposing them to highly contagious diseases, such as hepatitis and human immunodeficiency virus. In addition, needle phobia may prevent adequate treatment. To overcome these problems, we propose a needle-free injector based on thermocavitation.

Aim

Experimentally study the dynamics of vapor bubbles produced by thermocavitation inside a fully buried 3D fused silica chamber and the resulting high-speed jets emerging through a small nozzle made at the top of it. The injected volume can range from ∼ 0.1 to 2    μ L per shot. We also demonstrate that these jets have the ability to penetrate agar skin phantoms and ex-vivo porcine skin.

Approach

Through the use of a high-speed camera, the dynamics of liquid jets ejected from a microfluidic device were studied. Thermocavitation bubbles are generated by a continuous wave laser (1064 nm). The 3D chamber was fabricated by ultra-short pulse laser-assisted chemical etching. Penetration tests are conducted using agar gels (1%, 1.25%, 1.5%, 1.75%, and 2% concentrations) and porcine tissue as a model for human skin.

Result

High-speed camera video analysis showed that the average maximum bubble wall speed is about 10 to 25 m/s for almost any combination of pump laser parameters; however, a clever design of the chamber and nozzle enables one to obtain jets with an average speed of ∼ 70    m / s . The expelled volume per shot (0.1 to 2    μ l ) can be controlled by the pump laser intensity. Our injector can deliver up to 20 shots before chamber refill. Penetration of jets into agar of different concentrations and ex-vivo porcine skin is demonstrated.

Conclusions

The needle-free injectors based on thermocavitation may hold promise for commercial development, due to their cost and compactness.

Assessing engineered tissues and biomaterials using ultrasound imaging: In vitro and in vivo applications

Quantitative assessment of the structural, functional, and mechanical properties of engineered tissues and biomaterials is fundamental to their development for regenerative medicine applications. Ultrasound (US) imaging is a non-invasive, non-destructive, and cost-effective technique capable of longitudinal and quantitative monitoring of tissue structure and function across centimeter to sub-micron length scales. Here we present the fundamentals of US to contextualize its application for the assessment of biomaterials and engineered tissues, both in vivo and in vitro. We review key studies that demonstrate the versatility and broad capabilities of US for clinical and pre-clinical biomaterials research. Finally, we highlight emerging techniques that further extend the applications of US, including for ultrafast imaging of biomaterials and engineered tissues in vivo and functional monitoring of stem cells, organoids, and organ-on-a-chip systems in vitro.

Hybrid Nanoparticles of Citrate-Coated Manganese Ferrite and Gold Nanorods in Magneto-Optical Imaging and Thermal Therapy

The development of nanomaterials has drawn considerable attention in nanomedicine to advance cancer diagnosis and treatment over the last decades. Gold nanorods (GNRs) and magnetic nanoparticles (MNPs) have been known as commonly used nanostructures in biomedical applications due to their attractive optical properties and superparamagnetic (SP) behaviors, respectively. In this study, we proposed a simple combination of plasmonic and SP properties into hybrid NPs of citrate-coated manganese ferrite (Ci-MnFe₂O₄) and cetyltrimethylammonium bromide-coated GNRs (CTAB-GNRs). In this regard, two different samples were prepared: the first was composed of Ci-MnFe₂O₄ (0.4 wt%), and the second contained hybrid NPs of Ci-MnFe₂O₄ (0.4 wt%) and CTAB-GNRs (0.04 wt%). Characterization measurements such as UV-Visible spectroscopy and transmission electron microscopy (TEM) revealed electrostatic interactions caused by the opposing surface charges of hybrid NPs, which resulted in the formation of small nanoclusters. The performance of the two samples was investigated using magneto-motive ultrasound imaging (MMUS). The sample containing Ci-MnFe₂O₄_CTAB-GNRs demonstrated a displacement nearly two-fold greater than just using Ci-MnFe₂O₄; therefore, enhancing MMUS image contrast. Furthermore, the preliminary potential of these hybrid NPs was also examined in magnetic hyperthermia (MH) and photoacoustic imaging (PAI) modalities. Lastly, these hybrid NPs demonstrated high stability and an absence of aggregation in water and phosphate buffer solution (PBS) medium. Thus, Ci-MnFe₂O₄_CTAB-GNRs hybrid NPs can be considered as a potential contrast agent in MMUS and PAI and a heat generator in MH.

Microbubble dynamics and jetting near tissue-phantom biointerfaces

Precise excitation of cavitation is a promising mechanism for microsurgery procedures and targeted drug delivery enhancement. The underlying phenomenon of interest, jetting behaviour of oscillating cavitation bubbles, occurs due to near-surface interactions between the boundary, liquid, and bubble. Within this study we measured boundary effects on the cavitation bubble dynamics and morphology, with an emphasis on observation and measurement of jetting behaviour near tissue-phantom biointerfaces. An important mechanism of boundary poration has been observed using time-resolved optical microscopy and explained for different tissue-phantom surface densities and Young's modulus. Below a critical distance to the boundary, around γ = 1.0, the resulting jets penetrated the tissue-phantom, resulting in highly localized few micrometer diameter jets.

Nonlinear Characterization of Tissue Viscoelasticity With Acoustoelastic Attenuation of Shear Waves

Shear-wave elastography (SWE) measures shear-wave speed (SWS), which is related to the underlying shear modulus of soft tissue. SWS in soft tissue changes depending on the amount of external strain that soft tissue is subjected to due to the acoustoelastic (AE) phenomenon. In the literature, variations of SWS as a function of applied uniaxial strain were used for nonlinear characterization, assuming soft tissues to be elastic, although soft tissues are indeed viscoelastic in nature. Hence, nonlinear characterization using SWS alone is insufficient. In this work, we use SWS together with shear-wave attenuation (SWA) during incremental quasi-static compressions in order to derive biomechanical characterization based on the AE theory in terms of well-defined storage and loss moduli. As part of this study, we also quantify the effect of applied strain on measurements of SWS and SWA since such confounding effects need to be taken into account when using SWS and/or SWA, e.g., for staging a disease state, while such effects can also serve as an additional imaging biomarker. Our results from tissue-mimicking phantoms with varying oil percentages and ex vivo porcine liver experiments demonstrate the feasibility of our proposed methods. In both experiments, SWA was observed to decrease with applied strain. For 10% compression in ex vivo livers, shear-wave attenuation decreased, on average, by 28% (93 Np/m), while SWS increased, on average, by 20% (0.26 m/s).

Axially- and torsionally-polarized radially converging shear wave MRE in an anisotropic phantom made via Embedded Direct Ink Writing

Magnetic Resonance Elastography (MRE) is a non-invasive imaging method to quantitatively map the shear viscoelastic properties of soft tissues. In this study, Embedded Direct Ink Writing is used to fabricate a muscle mimicking anisotropic phantom that may serve as a standard for imaging studies of anisotropic materials. The technique allowed us to obtain a long shelf life silicone-based phantom expressing transverse isotropic mechanical properties. Another goal of the present investigation is to introduce a torsionally-polarized, radially-converging shear wave actuation method for MRE. The implemented design for this novel setup was first validated via its application to isotropic and homogeneous gelatin phantoms. Then, a comparison of the resulting complex wave images from axially- and torsionally-polarized MRE on the developed anisotropic phantom and on a skeletal muscle murine sample is presented, highlighting the value of using multiple actuation and motion encoding polarization directions when studying anisotropic materials.

Measuring viscoelastic parameters in Magnetic Resonance Elastography: a comparison at high and low magnetic field intensity

Magnetic Resonance Elastography (MRE) is a non-invasive imaging technique which involves motion-encoding MRI for the estimation of the shear viscoelastic properties of soft tissues through the study of shear wave propagation. The technique has been found informative for disease diagnosis, as well as for monitoring of the effects of therapies. The development of MRE and its validation have been supported by the use of tissue-mimicking phantoms. In this paper we present our new MRE protocol using a low magnetic field tabletop MRI device at 0.5 T and sinusoidal uniaxial excitation in a geometrical focusing condition. Results obtained for gelatin are compared to those previously obtained using high magnetic field MRE at 11.7 T. A multi-frequency investigation is also provided via a comparison of commonly used rheological models: Maxwell, Springpot, Voigt, Zener, Jeffrey, fractional Voigt and fractional Zener. Complex shear modulus values were comparable when processed from images acquired with the tabletop low field scanner and the high field scanner. This study serves as a validation of the presented tabletop MRE protocol and paves the way for MRE experiments on ex-vivo tissue samples in both normal and pathological conditions.

A Novel Theranostic Platform: Integration of Magnetomotive and Thermal Ultrasound Imaging With Magnetic Hyperthermia

OBJECTIVE

Nanotheranostic systems integrate therapeutic and diagnostic procedures using nanotechnology. This type of approach has enabled the development of methods for early detection and treatment of different pathologies. Magnetic hyperthermia (MH) has been proposed as an alternative or complementary method of cancer therapy. However, challenges such as delivering and localizing the magnetic nanoparticles (MNPs) within tissues and monitoring the temperature during the treatment hinder this technique to be effectively translated into a clinical routine. Therefore, in this study a theranostic platform has been proposed and examined to address two main issues, localizing MNPs and real-time temperature monitoring, for preclinical MH.

METHODS

The system integrates magnetomotive (MMUS) and thermal ultrasound imaging with MH. An ultrasound device was used to acquire MMUS images to detect MNPs, and ultrasound thermometry to monitor the temperature. This platform was designed such that a single coil generated the magnetic field for MMUS and MH. The feasibility of the system was examined using a tissue mimicking phantom containing an inclusion filled with zinc substituted magnetite NPs.

RESULTS

These MNPs were effectively used as contrast agent for MMUS and to generate heat during MH. In addition to localizing MNPs, real-time two-dimensional temperature maps were obtained with substantial concordance (ρ_c > 0.97) with invasive measurements using fiber optic thermometer. The heating rate was proportional to the displacements in MMUS (r = 0.92).

CONCLUSION

Ultrasound thermometry was successfully used to monitor the temperature during MH. In addition, it was shown that acquiring MMUS images prior to MH can qualitatively predict the temperature distribution of the MNP-laden regions.

Multiangle Long-Axis Lateral Illumination Photoacoustic Imaging Using Linear Array Transducer

Photoacoustic imaging (PAI) combines optical contrast with ultrasound spatial resolution and can be obtained up to a depth of a few centimeters. Hand-held PAI systems using linear array usually operate in reflection mode using a dark-field illumination scheme, where the optical fiber output is attached to both sides of the elevation plane (short-axis) of the transducer. More recently, bright-field strategies where the optical illumination is coaxial with acoustic detection have been proposed to overcome some limitations of the standard dark-field approach. In this paper, a novel multiangle long-axis lateral illumination is proposed. Monte Carlo simulations were conducted to evaluate light delivery for three different illumination schemes: bright-field, standard dark-field, and long-axis lateral illumination. Long-axis lateral illumination showed remarkable improvement in light delivery for targets with a width smaller than the transducer lateral dimension. A prototype was developed to experimentally demonstrate the feasibility of the proposed approach. In this device, the fiber bundle terminal ends are attached to both sides of the transducer’s long-axis and the illumination angle of each fiber bundle can be independently controlled. The final PA image is obtained by the coherent sum of subframes acquired using different angles. The prototype was experimentally evaluated by taking images from a phantom, a mouse abdomen, forearm, and index finger of a volunteer. The system provided light delivery enhancement taking advantage of the geometry of the target, achieving sufficient signal-to-noise ratio at clinically relevant depths.

C-Elastography: In Vitro Feasibility Phantom Study

C-Elastography (CE) is a new ultrasound technique that locally maps the non-linear elasticity of soft tissue using low-frequency (150-250 Hz) shear waves generated by the acoustic radiation force (ARF). CE is based on a recent finding that the magnitude of the ARF in an isotropic tissue-like solid is related linearly to a third-order modulus of elasticity, C, which is responsible for the coupling between deviatoric and volumetric constitutive behaviors. The main objective of the work described here was to examine the feasibility of using and performance of C-elastography in differentiating and characterizing soft tissue via a pilot study on ex vivo tissue and tissue-mimicking inclusions cast in a gelatin block. In this vein, the CE technique deploys a combination of ultrasound motion sensing and 3-D visco-elastodynamic simulation to estimate the non-linear modulus C. As ultrasound focusing inherently confines the ARF to a small region, CE provides the means for measuring C within O(mm3) volumes. Equipped with such data analysis, we performed in vitro CE experiments on agar-based, xenograft and normal breast tissue samples embedded in a gelatin matrix. The compound C-elastograms indicate marked (and sharp) C-contrast, with average values of 1.9 and 5.6 at push points inside the featured soft and hard inclusions, respectively.

Walled vessel-mimicking phantom for ultrasound imaging using 3D printing with a water-soluble filament: design principle, fluid-structure interaction (FSI) simulation, and experimental validation

The geometry and stiffness of a vessel are pertinent to blood dynamics and vessel wall mechanical behavior and may alter in diseased conditions. Ultrasound-based ultrafast Doppler (uDoppler) imaging and shear wave imaging (SWI) techniques have been extensively exploited for the assessment of vascular hemodynamics and mechanics. Their performance is conventionally validated on vessel-mimicking phantoms (VMPs) prior to their clinical use. Compared with commercial ones, customized VMPs are favored for research use because of their wider range of material properties, more complex lumen geometries, or wall structures. Fused deposition modeling (FDM) 3D printing technique with plastic filaments is a promising method for making VMPs with a complex vessel lumen. However, it may require a toxic solvent or a long dissolution time currently. In this paper, we present a safe, efficient and geometrically flexible method where FDM 3D printing with a water-soluble polyvinyl alcohol (PVA) filament is exploited to fabricate a walled three-branch VMP (VMP-I). As a key step in fabrication, to avoid dissolution of the PVA-printed vessel core by the solution of the tissue-mimicking material, paraffin wax was used for isolation. Paraffin wax is easy to coat (i.e. without any special equipment), of satisfactory thickness (∼0.1 mm), chemically stable, and easy to remove after fabrication, thus making the proposed method practicable for ultrasound imaging studies. VMP-I was examined by B-mode imaging and power Doppler imaging (PDI) to verify complete dissolution of PVA-printed vessel core in its lumen, confirming good fabrication quality. The flow velocities in VMP-I were estimated by uDoppler imaging with a -0.8% difference, and the shear wave propagation speeds for the same phantom were estimated by SWI with a -6.03% difference when compared with fluid-structure interaction (FSI) simulation results. A wall-less VMP of a scaled and simplified coronary arterial network (VMP-II) was additionally fabricated and examined to test the capability of the proposed method for a complex lumen geometry. The proposed fabrication method for customized VMPs is foreseen to facilitate the development of ultrasound imaging techniques for blood vessels.

Predictive Model for Designing Soft-Tissue Mimicking Ultrasound Phantoms With Adjustable Elasticity

The use of mechanically representative phantoms is important for experimental validation in ultrasound (US) imaging, elastography, and image registration. This article proposes a model to predict the elastic modulus of a soft tissue-mimicking phantom based on two very easily controllable parameters: gelatin concentration and refrigeration duration. The model has been validated on small- and large-scale phantoms; it provides a good prediction of the elastic modulus in both cases (error < 16.2%). The tissue-mimicking phantom is made following a low-cost and simple fabrication procedure using commercial household gelatin with psyllium hydrophilic mucilloid fiber to obtain echogenicity. A large range of elastic properties was obtained (15-100kPa) by adjusting the gelatin concentration between 5% and 20% (g/mL) and the refrigeration time of the sample between 2 and 168 h, allowing to mimic normal and pathological human soft tissues. The phantom's acoustic properties (velocity, attenuation, and acoustic impedance) are also assessed using the American Institute of Ultrasound in Medicine (AIUM) standard.

Mechanical characterisation of an organic phantom candidate for breast tissue

The aim of this study is to analyse the mechanical properties of agarose, which can be used as a mechanical phantom for breast tissues. In general, tissue mimicking materials may be used to achieve a better understanding of the structure and properties of tissues and organs; this work shares these objectives. The phantom can be used as a tissue surrogate with realistic mechanical behaviour for biomechanical applications. To validate agarose as a suitable mechanical phantom for breast tissues, indentation tests were performed in homogeneous, rectangular agarose blocks. Blocks with different agarose concentrations were moulded and tested. An estimation of the material stiffness was then compared with experimental data on different breast tissues’ types found in literature. The phantom stiffness increased consistently with agarose concentration and stress. The results show that agarose-based mechanical phantoms of stiffer tissues require higher agarose concentrations (0.5% and 0.6%). In contrast, normal tissues can be mimicked with 0.3% and 0.4% of agarose. In addition, it was observed that preconditioning affects the mechanical properties of the gel, as it had already been shown in literature for breast tissues.

Evaluation of the Effect of Tissue Compression on the Results of Shear Wave Elastography Measurements

Shear wave imaging is considered to be more precise and less operator dependent when compared with strain imaging. It enables quantitative and reproducible data (Young's modulus of the imaged tissue). However, results of shear wave imaging can be affected by a variety of different factors. The aim of this study is to evaluate the effect of the pressure applied by the ultrasound probe during examination on the measured values of Young's modulus. The effect of the tissue compression on the results of the real-time shear wave elastography was evaluated via the gelatine phantom measurements, via the ex vivo experiments with pig liver, and via the in vivo measurements of the thyroid gland stiffness on healthy volunteers. The results of our measurements confirmed that the measured value of Young's modulus increases with the increasing pressure applied on the imaged object. The highest increase was observed during the ex vivo experiments (400%), and the lowest increase was detected in the case of the phantom measurements (8%). A two- to threefold increase in Young's modulus was observed between the minimum and maximum pressure in the case of the in vivo elastography measurements of thyroid gland. The Veronda-Westman theoretical model was used for the description of the tissue nonlinearity. We conclude that tissue compression by the force exerted on the probe can significantly affect the results of the real-time shear wave elastography measurements. Minimum pressure should be used when measuring the absolute value of Young's modulus of superficial organs.

Arterial Phantoms with Regional Variations in Wall Stiffness and Thickness

Regional wall stiffening and thickening are two common pathological features of arteries. To account for these two features, we developed a new arterial phantom design framework to facilitate the development of vessel models that contain a lesion segment whose wall stiffness and thickness differ from those of other segments. This new framework is based on multi-part injection molding principles that sequentially casted the lesion segment and the flank segments of the vessel model using molding parts devised with computer-aided design tools. The vessel-mimicking material is created from polyvinyl alcohol cryogel, and its acoustic properties are similar to those of arteries. As a case demonstration, we fabricated a stenosed three-segment phantom composed of a central lesion segment (5.1-mm diameter, 1.95-mm wall thickness, 212.6-kPa elastic modulus) and two flank segments (6.0-mm diameter, 1.5-mm wall thickness, 133.7-kPa elastic modulus). B-mode imaging confirmed the difference in thickness between the lesion segment and flank segments of the phantom. Also, Doppler-based vessel wall displacement analysis revealed that when pulsatile flow was fed through the phantom (carotid pulse; 27 mL/s peak flow rate), the lesion segment distended less compared with the flank segments. Specifically, the three-beat averaged peak wall displacement in the lesion segment was measured as 0.28 mm, and it was significantly smaller than that of the flank segments (0.60 mm). It is anticipated that this new multi-segment arterial phantom can serve as a performance testbed for the evaluation of local arterial stiffness estimation algorithms.

Application of Acoustoelasticity to Evaluate Nonlinear Modulus in Ex Vivo Kidneys

Currently, dynamic elastography techniques estimate the linear elastic shear modulus of different body tissues. New methods that investigate other properties of soft tissues such as anisotropy, viscosity, and shear nonlinearity would provide more information about the structure and function of the tissue and might provide a better contrast than tissue stiffness and hence provide more effective diagnostic tools for some diseases. It has previously been shown that shear wave velocity in a medium changes due to an applied stress, a phenomenon called acoustoelasticity (AE). Applying a stress to compress a medium while measuring the shear wave velocity versus strain provides data with which the third-order nonlinear shear modulus can be estimated. To evaluate the feasibility of estimating , we evaluated ten ex vivo porcine kidneys embedded in 10% porcine gelatin to mimic the case of a transplanted kidney. Under assumptions of an elastic incompressible medium for AE measurements, the shear modulus was quantified at each compression level and the applied strain was assessed by measuring the change in the thickness of the kidney cortex. Finally, was calculated by applying the AE theory. Our results demonstrated that it is possible to estimate a nonlinear shear modulus by monitoring the changes in strain and due to kidney deformation. The magnitudes of are higher when the compression is performed progressively and when using a plate attached to the transducer. Nevertheless, the values obtained for are similar to those previously reported in the literature for breast tissue.

NADPH oxidase mediates microtubule alterations and diaphragm dysfunction in dystrophic mice

Skeletal muscle from mdx mice is characterized by increased Nox2 ROS, altered microtubule network, increased muscle stiffness, and decreased muscle/respiratory function. While microtubule de-tyrosination has been suggested to increase stiffness and Nox2 ROS production in isolated single myofibers, its role in altering tissue stiffness and muscle function has not been established. Because Nox2 ROS production is upregulated prior to microtubule network alterations and ROS affect microtubule formation, we investigated the role of Nox2 ROS in diaphragm tissue microtubule organization, stiffness and muscle/respiratory function. Eliminating Nox2 ROS prevents microtubule disorganization and reduces fibrosis and muscle stiffness in mdx diaphragm. Fibrosis accounts for the majority of variance in diaphragm stiffness and decreased function, implicating altered extracellular matrix and not microtubule de-tyrosination as a modulator of diaphragm tissue function. Ultimately, inhibiting Nox2 ROS production increased force and respiratory function in dystrophic diaphragm, establishing Nox2 as a potential therapeutic target in Duchenne muscular dystrophy.

Quantitative characterization of viscoelastic behavior in tissue-mimicking phantoms and ex vivo animal tissues

Viscoelasticity of soft tissue is often related to pathology, and therefore, has become an important diagnostic indicator in the clinical assessment of suspect tissue. Surgeons, particularly within head and neck subsites, typically use palpation techniques for intra-operative tumor detection. This detection method, however, is highly subjective and often fails to detect small or deep abnormalities. Vibroacoustography (VA) and similar methods have previously been used to distinguish tissue with high-contrast, but a firm understanding of the main contrast mechanism has yet to be verified. The contributions of tissue mechanical properties in VA images have been difficult to verify given the limited literature on viscoelastic properties of various normal and diseased tissue. This paper aims to investigate viscoelasticity theory and present a detailed description of viscoelastic experimental results obtained in tissue-mimicking phantoms (TMPs) and ex vivo tissues to verify the main contrast mechanism in VA and similar imaging modalities. A spherical-tip micro-indentation technique was employed with the Hertzian model to acquire absolute, quantitative, point measurements of the elastic modulus (E), long term shear modulus (η), and time constant (τ) in homogeneous TMPs and ex vivo tissue in rat liver and porcine liver and gallbladder. Viscoelastic differences observed between porcine liver and gallbladder tissue suggest that imaging modalities which utilize the mechanical properties of tissue as a primary contrast mechanism can potentially be used to quantitatively differentiate between proximate organs in a clinical setting. These results may facilitate more accurate tissue modeling and add information not currently available to the field of systems characterization and biomedical research.

Anatomically realistic ultrasound phantoms using gel wax with 3D printed moulds

Here we describe methods for creating tissue-mimicking ultrasound phantoms based on patient anatomy using a soft material called gel wax. To recreate acoustically realistic tissue properties, two additives to gel wax were considered: paraffin wax to increase acoustic attenuation, and solid glass spheres to increase backscattering. The frequency dependence of ultrasound attenuation was well described with a power law over the measured range of 3–10 MHz. With the addition of paraffin wax in concentrations of 0 to 8 w/w%, attenuation varied from 0.72 to 2.91 dB cm⁻¹ at 3 MHz and from 6.84 to 26.63 dB cm⁻¹ at 10 MHz. With solid glass sphere concentrations in the range of 0.025–0.9 w/w%, acoustic backscattering consistent with a wide range of ultrasonic appearances was achieved. Native gel wax maintained its integrity during compressive deformations up to 60%; its Young’s modulus was 17.4  ±  1.4 kPa. The gel wax with additives was shaped by melting and pouring it into 3D printed moulds. Three different phantoms were constructed: a nerve and vessel phantom for peripheral nerve blocks, a heart atrium phantom, and a placental phantom for minimally-invasive fetal interventions. In the first, nerves and vessels were represented as hyperechoic and hypoechoic tubular structures, respectively, in a homogeneous background. The second phantom comprised atria derived from an MRI scan of a patient with an intervening septum and adjoining vena cavae. The third comprised the chorionic surface of a placenta with superficial fetal vessels derived from an image of a post-partum human placenta. Gel wax is a material with widely tuneable ultrasound properties and mechanical characteristics that are well suited for creating patient-specific ultrasound phantoms in several clinical disciplines.

Three-dimensional Ultrasound Elasticity Imaging on an Automated Breast Volume Scanning System

Ultrasound elasticity imaging has demonstrated utility in breast imaging, but it is typically performed with handheld transducers and two-dimensional imaging. Two-dimensional (2D) elastography images tissue stiffness of only a plane and hence suffers from errors due to out-of-plane motion, whereas three-dimensional (3D) data acquisition and motion tracking can be used to track out-of-plane motion that is lost in 2D elastography systems. A commercially available automated breast volume scanning system that acquires 3D ultrasound data with precisely controlled elevational movement of the 1D array ultrasound transducer was employed in this study. A hybrid guided 3D motion-tracking algorithm was developed that first estimated the displacements in one plane using a modified quality-guided search method, and then performed an elevational guided-search for displacement estimation in adjacent planes. To assess the performance of the method, 3D radiofrequency echo data were acquired with this system from a phantom and from an in vivo human breast. For both experiments, the axial displacement fields were smooth and high cross-correlation coefficients were obtained in most of the tracking region. The motion-tracking performance of the new method was compared with a correlation-based exhaustive-search method. For all motion-tracking volume pairs, the average motion-compensated cross-correlation values obtained by the guided-search motion-tracking method were equivalent to those by the exhaustive-search method, and the computation time was about a factor of 10 lesser. Therefore, the proposed 3D ultrasound elasticity imaging method was a more efficient approach to produce a high quality of 3D ultrasound strain image.

Direct Error in Constitutive Equation Formulation for Plane stress Inverse Elasticity Problem

We present a new computational formulation for inverse problems in elasticity with full field data. The formulation is a variant of an error in the constitutive equation formulation, but allows direct solution for the modulus field and accommodates discontinuous strain fields. The development of the formulation is motivated by the relatively poor performance of current direct formulations, reported so far in literature, in dealing with discontinuities in the strain and material property distribution. The formulation relies on minimizing the error in the constitutive equation, and a momentum equation constraint. Numerical results on model problems show that the formulation is capable handling discontinuous, and noisy strain fields, and also converging with mesh refinement for continuous and discontinuous material property distributions. The application to reconstruct the elastic modulus distribution in solid breast tumors is shown.

Recovering vector displacement estimates in quasistatic elastography using sparse relaxation of the momentum equation

We consider the problem of estimating the 2D vector displacement field in a heterogeneous elastic solid deforming under plane stress conditions. The problem is motivated by applications in quasistatic elastography. From precise and accurate measurements of one component of the 2D vector displacement field and very limited information of the second component, the method reconstructs the second component quite accurately. No a priori knowledge of the heterogeneous distribution of material properties is required. This method relies on using a special form of the momentum equations to filter ultrasound displacement measurements to produce more precise estimates. We verify the method with applications to simulated displacement data. We validate the method with applications to displacement data measured from a tissue mimicking phantom, and in-vivo data; significant improvements are noticed in the filtered displacements recovered from all the tests. In verification studies, error in lateral displacement estimates decreased from about 50% to about 2%, and strain error decreased from more than 250% to below 2%.

Improvement of Lesion Detection by Complete Angular Compound Ultrasonic Elastography

Quasi-static ultrasound elastography is an emerging diagnostic imaging modality for determining the stiffness of pathologically changed soft tissues, which do not show significant differences in acoustic impedance for B-mode imaging. Although some methods were applied to improve the signal-to-noise ratio (SNRe) and contrast-to-noise ratio (CNRe) of the constructed elastogram, nonuniform strain distribution at the internal boundary of a hard inclusion, even with the uniform displacement on the surface, is an inherent mechanical effect and results in distortion at the detected lesion boundary. To overcome such stress concentrations, a new elastographic modality was proposed, where the elastograms from different angles throughout 360° were compounded. The strain field and subsequent ultrasound images were calculated using the finite element method (FEM) and Field II, respectively, from which the elastograms were constructed. The performance of complete angular compound elastography with varied interval angles, lesion sizes, and ratios of Young's moduli of the lesion to the background was simulated and compared with that of conventional axial strain elastography. It is found that viewing the lesion from only about 10 angles (interval of 36°) would significantly improve the image quality of elastogram (increasing SNRe by at least 13% and CNRe by at least 5.8 dB), reduce the lesion distortion in the lateral direction, and enhance the sensitivity, resolution, and accuracy of lesion detection. A preliminary phantom study showed similar improvements. Altogether, complete angular compound elastography improves the elastogram quality and reduces the mechanical effects in lesion detection.

Inferring spatial variations of microstructural properties from macroscopic mechanical response

Disease alters tissue microstructure, which in turn affects the macroscopic mechanical properties of tissue. In elasticity imaging, the macroscopic response is measured and is used to infer the spatial distribution of the elastic constitutive parameters. When an empirical constitutive model is used, these parameters cannot be linked to the microstructure. However, when the constitutive model is derived from a microstructural representation of the material, it allows for the possibility of inferring the local averages of the spatial distribution of the microstructural parameters. This idea forms the basis of this study. In particular, we first derive a constitutive model by homogenizing the mechanical response of a network of elastic, tortuous fibers. Thereafter, we use this model in an inverse problem to determine the spatial distribution of the microstructural parameters. We solve the inverse problem as a constrained minimization problem and develop efficient methods for solving it. We apply these methods to displacement fields obtained by deforming gelatin-agar co-gels and determine the spatial distribution of agar concentration and fiber tortuosity, thereby demonstrating that it is possible to image local averages of microstructural parameters from macroscopic measurements of deformation.

Quantitative photoacoustic elastography in humans

We report quantitative photoacoustic elastography (QPAE) capable of measuring Young’s modulus of biological tissue in vivo in humans. By combining conventional PAE with a stress sensor having known stress–strain behavior, QPAE can simultaneously measure strain and stress, from which Young’s modulus is calculated. We first demonstrate the feasibility of QPAE in agar phantoms with different concentrations. The measured Young’s modulus values fit well with both the empirical expectation based on the agar concentrations and those measured in an independent standard compression test. Next, QPAE was applied to quantify the Young’s modulus of skeletal muscle in vivo in humans, showing a linear relationship between muscle stiffness and loading. The results demonstrated the capability of QPAE to assess the absolute elasticity of biological tissue noninvasively in vivo in humans, indicating its potential for tissue biomechanics studies and clinical applications.

Comparison between shear wave dispersion magneto motive ultrasound and transient elastography for measuring tissue-mimicking phantom viscoelasticity

Several methods have been developed over the last several years to analyze the mechanical properties of soft tissue. Elastography, for example, was proposed to evaluate soft tissue stiffness in an attempt to reduce the need for invasive procedures, such as breast biopsies; however, its qualitative nature and the fact that it is operator-dependent have proven to be limitations of the technique. Quantitative shearwave- based techniques have been proposed to obtain information about tissue stiffness independent of the operator. This paper describes shear wave dispersion magnetomotive ultrasound (SDMMUS), a new shear-wave-based method in which a viscoelastic medium labeled with iron oxide nanoparticles is displaced by an external tone burst magnetic field. As in magnetomotive ultrasound (MMUS), SDMMUS uses ultrasound to detect internal mechanical vibrations induced by the interaction between a magnetic field and magnetic nanoparticles. These vibrations generated shear waves that were evaluated to estimate the viscoelastic properties of tissue-mimicking phantoms. These phantoms were manufactured with different concentrations of gelatin and labeled with iron oxide nanoparticles. The elasticity and viscosity obtained with SDMMUS agreed well with the results obtained by traditional ultrasound-based transient elastography.

Quantitative micro-elastography: imaging of tissue elasticity using compression optical coherence elastography

Probing the mechanical properties of tissue on the microscale could aid in the identification of diseased tissues that are inadequately detected using palpation or current clinical imaging modalities, with potential to guide medical procedures such as the excision of breast tumours. Compression optical coherence elastography (OCE) maps tissue strain with microscale spatial resolution and can delineate microstructural features within breast tissues. However, without a measure of the locally applied stress, strain provides only a qualitative indication of mechanical properties. To overcome this limitation, we present quantitative micro-elastography, which combines compression OCE with a compliant stress sensor to image tissue elasticity. The sensor consists of a layer of translucent silicone with well-characterized stress-strain behaviour. The measured strain in the sensor is used to estimate the two-dimensional stress distribution applied to the sample surface. Elasticity is determined by dividing the stress by the strain in the sample. We show that quantification of elasticity can improve the ability of compression OCE to distinguish between tissues, thereby extending the potential for inter-sample comparison and longitudinal studies of tissue elasticity. We validate the technique using tissue-mimicking phantoms and demonstrate the ability to map elasticity of freshly excised malignant and benign human breast tissues.

A simulation technique for 3D MR-guided acoustic radiation force imaging

PURPOSE

In magnetic resonance-guided focused ultrasound (MRgFUS) therapies, the in situ characterization of the focal spot location and quality is critical. MR acoustic radiation force imaging (MR-ARFI) is a technique that measures the tissue displacement caused by the radiation force exerted by the ultrasound beam. This work presents a new technique to model the displacements caused by the radiation force of an ultrasound beam in a homogeneous tissue model.

METHODS

When a steady-state point-source force acts internally in an infinite homogeneous medium, the displacement of the material in all directions is given by the Somigliana elastostatic tensor. The radiation force field, which is caused by absorption and reflection of the incident ultrasound intensity pattern, will be spatially distributed, and the tensor formulation takes the form of a convolution of a 3D Green's function with the force field. The dynamic accumulation of MR phase during the ultrasound pulse can be theoretically accounted for through a time-of-arrival weighting of the Green's function. This theoretical model was evaluated experimentally in gelatin phantoms of varied stiffness (125-, 175-, and 250-bloom). The acoustic and mechanical properties of the phantoms used as parameters of the model were measured using independent techniques. Displacements at focal depths of 30- and 45-mm in the phantoms were measured by a 3D spin echo MR-ARFI segmented-EPI sequence.

RESULTS

The simulated displacements agreed with the MR-ARFI measured displacements for all bloom values and focal depths with a normalized RMS difference of 0.055 (range 0.028-0.12). The displacement magnitude decreased and the displacement pattern broadened with increased bloom value for both focal depths, as predicted by the theory.

CONCLUSIONS

A new technique that models the displacements caused by the radiation force of an ultrasound beam in a homogeneous tissue model theory has been rigorously validated through comparison with experimentally obtained 3D displacement data in homogeneous gelatin phantoms using a 3D MR-ARFI sequence. The agreement of the experimentally measured and simulated results demonstrates the potential to use MR-ARFI displacement data in MRgFUS therapies.

Characterization of biomechanical properties of agar based tissue mimicking phantoms for ultrasound stiffness imaging techniques

Pathological changes of the body have been observed to change the mechanical properties of soft tissue types which can be imaged by ultrasound elastography. Though initial clinical results using ultrasound elastography in detection of tumors are promising, quantification of signal to noise ratio, resolution and strain image patterns are the best achieved under a controlled study using phantoms with similar biomechanical properties of normal and abnormal tissues. The purpose of this work is to characterize the biomechanical properties of agar based tissue mimicking phantoms by varying the agar concentration from 1.7 to 6.6% by weight and identify the optimum property to be used in classification of cancerous tissues. We performed quasi-static uniaxial compression test under a strain rate of 0.5mm/min up to 15% strain and measured Young's modulus of phantom samples which are from 50kPa to 450kPa. Phantoms show nonlinear stress-strain characteristics at finite strain which were characterized using hyperelastic parameters by fitting Neo-Hookean, Mooney Rivlin, Ogden and Veronda Westmann models. We also investigated viscoelastic parameters of the samples by conducting oscillatory shear rheometry at various precompression levels (2-5%). Loss modulus values are always less than storage modulus which represents the behavior of soft tissues. The increase in agar concentration increases the shear modulus of the samples as well as decreases the linear viscoelastic region. The results suggest that dynamic shear modul are more promising than linear and nonlinear elastic modul in differentiation of various classes of abnormal tissues.

Characterizing phantom arteries with multi-channel laser ultrasonics and photo-acoustics

Multi-channel photo-acoustic and laser ultrasonic waves are used to sense the characteristics of proxies for healthy and diseased vessels. The acquisition system is non-contacting and non-invasive with a pulsed laser source and a laser vibrometer detector. As the wave signatures of our targets are typically low in amplitude, we exploit multi-channel acquisition and processing techniques. These are commonly used in seismology to improve the signal-to-noise ratio of data. We identify vessel proxies with a diameter on the order of 1 mm, at a depth of 18 mm. Variations in scattered and photo-acoustic signatures are related to differences in vessel wall properties and content. The methods described have the potential to improve imaging and better inform interventions for atherosclerotic vessels, such as the carotid artery.

Development of oil-in-gelatin phantoms for viscoelasticity measurement in ultrasound shear wave elastography

Because tissues consist of solid and fluid materials, their mechanical properties should be characterized in terms of both elasticity and viscosity. Although the elastic properties of tissue-mimicking phantoms have been extensively studied and well characterized in commercially available phantoms, their viscous properties have not been fully investigated. In this article, a set of 14 tissue-mimicking phantoms with different concentrations of gelatin and castor oil were fabricated and characterized in terms of acoustic and viscoelastic properties. The results indicate that adding castor oil to gelatin phantoms decreases shear modulus, but increases shear wave dispersion. For 3% gelatin phantoms containing 0%, 10%, 20% and 40% oil, the measured shear moduli are 2.01 ± 0.26, 1.68 ± 0.25, 1.10 ± 0.22 and 0.88 ± 0.17 kPa, and the Voigt-model coupled shear viscosities are 0.60 ± 0.11, 0.89 ± 0.07, 1.05 ± 0.11 and 1.06 ± 0.13 Pa·s, respectively. The results also confirm that increasing the gelatin concentration increases shear modulus. For phantoms containing 3%, 4%, 5%, 6% and 7% gelatin, the measured shear moduli are 2.01 ± 0.26, 3.10 ± 0.34, 4.18 ± 0.84, 8.05 ± 1.00 and 10.24 ± 1.80 kPa at 0% oil and 1.10 ± 0.22, 1.97 ± 0.20, 3.13 ± 0.63, 4.60 ± 0.60 and 8.43 ± 1.39 kPa at 20% oil, respectively. The phantom recipe developed in this study can be used in validating ultrasound shear wave elastography techniques for soft tissues.

Frequency dependence of laser ultrasonic SAW phase velocities measurements

Advances in the field of laser ultrasonics have opened up new possibilities in applications in many areas. This paper verifies the relationship between phase velocities of different materials, including hard solid and soft solid, and the frequency range of SAW signal. We propose a novel approach that utilizes a low coherence interferometer to detect the laser-induced surface acoustic waves (SAWs). A Nd:YAG focused laser line-source is applied to steel, iron, plastic plates and a 3.5% agar-agar phantom. The generated SAW signals are detected by a time domain low coherence interferometry system. SAW phase velocity dispersion curves were calculated, from which the elasticity of the specimens was evaluated. The relationship between frequency content and phase velocities was analyzed. We show that the experimental results agreed well with those of the theoretical expectations.

A nonlinear elasticity phantom containing spherical inclusions

The strain image contrast of some in vivo breast lesions changes with increasing applied load. This change is attributed to differences in the nonlinear elastic properties of the constituent tissues suggesting some potential to help classify breast diseases by their nonlinear elastic properties. A phantom with inclusions and long-term stability is desired to serve as a test bed for nonlinear elasticity imaging method development, testing, etc. This study reports a phantom designed to investigate nonlinear elastic properties with ultrasound elastographic techniques. The phantom contains four spherical inclusions and was manufactured from a mixture of gelatin, agar and oil. The phantom background and each of the inclusions have distinct Young's modulus and nonlinear mechanical behavior. This phantom was subjected to large deformations (up to 20%) while scanning with ultrasound, and changes in strain image contrast and contrast-to-noise ratio between inclusion and background, as a function of applied deformation, were investigated. The changes in contrast over a large deformation range predicted by the finite element analysis (FEA) were consistent with those experimentally observed. Therefore, the paper reports a procedure for making phantoms with predictable nonlinear behavior, based on independent measurements of the constituent materials, and shows that the resulting strain images (e.g., strain contrast) agree with that predicted with nonlinear FEA.

Evaluating elastic properties of heterogeneous soft tissue by surface acoustic waves detected by phase-sensitive optical coherence tomography

The combined use of surface acoustic wave (SAW) and phase-sensitive optical coherence tomography (PhS-OCT) is useful to evaluate the elasticity of layered biological tissues, such as normal skin. However, the pathological tissue is often originated locally, leading to the alternation of mechanical properties along both axial and lateral directions. We present a feasibility study on whether the SAW technique is sensitive to detect the alternation of mechanical property along the lateral direction within tissue, which is important for clinical utility of this technique to localize diseased tissue. Experiments are carried out on purposely designed tissue phantoms and ex vivo chicken breast samples, simulating the localized change of elasticity. A PhS-OCT system is employed not only to provide the ultra-high sensitive measurement of the generated surface waves on the tissue surface, but also to provide the real time imaging of the tissue to assist the elasticity evaluation of the heterogeneous tissue. The experimental results demonstrate that with PhS-OCT used as a pressure sensor, the SAW is highly sensitive to the elasticity change of the specimen in both vertical and lateral directions with a sensing depth of ∼5 mm with our current system setup, thus promising its useful clinical applications where the quantitative elasticity of localized skin diseases is needed to aid in diagnosis and treatment.

Quantitative imaging of nonlinear shear modulus by combining static elastography and shear wave elastography

The study of new tissue mechanical properties such as shear nonlinearity could lead to better tissue characterization and clinical diagnosis. This work proposes a method combining static elastography and shear wave elastography to derive the nonlinear shear modulus by applying the acoustoelasticity theory in quasi-incompressible soft solids. Results demonstrate that by applying a moderate static stress at the surface of the investigated medium, and by following the quantitative evolution of its shear modulus, it is possible to accurately and quantitatively recover the local Landau (A) coefficient characterizing the shear nonlinearity of soft tissues.

Recent results in nonlinear strain and modulus imaging

We report a summary of recent developments and current status of our team's efforts to image and quantify in vivo nonlinear strain and tissue mechanical properties. Our work is guided by a focus on applications to cancer diagnosis and treatment using clinical ultrasound imaging and quasi-static tissue deformations. We review our recent developments in displacement estimation from ultrasound image sequences. We discuss cross correlation approaches, regularized optimization approaches, guided search methods, multiscale methods, and hybrid methods. Current implementations can return results of high accuracy in both axial and lateral directions at several frames per second.We compare several strain estimators. Again we see a benefit from a regularized optimization approach. We then discuss both direct and iterative methods to reconstruct tissue mechanical property distributions from measured strain and displacement fields. We review the formulation, discretization, and algorithmic considerations that come into play when attempting to infer linear and nonlinear elastic properties from strain and displacement measurements. Finally we illustrate our progress with example applications in breast disease diagnosis and tumor ablation monitoring. Our current status shows that we have demonstrated quantitative determination of nonlinear parameters in phantoms and in vivo, in the context of 2D models and data. We look forward to incorporating 3D data from 2D transducer arrays to noninvasively create calibrated 3D quantitative maps of nonlinear elastic properties of breast tissues in vivo.

Elastic properties of soft tissue-mimicking phantoms assessed by combined use of laser ultrasonics and low coherence interferometry

Advances in the field of laser ultrasonics have opened up new possibilities in medical applications. This paper evaluates this technique as a method that would allow for rapid characterization of the elastic properties of soft biological tissue. In doing so, we propose a novel approach that utilizes a low coherence interferometer to detect the laser-induced surface acoustic waves (SAW) from the tissue-mimicking phantoms. A Nd:YAG focused laser line-source is applied to one- and two-layer tissue-mimicking agar-agar phantoms, and the generated SAW signals are detected by a time domain low coherence interferometry system. SAW phase velocity dispersion curves are calculated, from which the elasticity of the specimens is evaluated. We show that the experimental results agree well with those of the theoretical expectations. This study is the first report that a laser-generated SAW phase velocity dispersion technique is applied to soft materials. This technique may open a way for laser ultrasonics to detect the mechanical properties of soft tissues, such as skin.