1
|
Mislati R, Uccello TP, Lin Z, Iliza KT, Toussaint KC, Gerber SA, Doyley MM. Shear wave elastography can stratify rectal cancer response to short-course radiation therapy. Sci Rep 2023; 13:16149. [PMID: 37752156 PMCID: PMC10522682 DOI: 10.1038/s41598-023-43383-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2023] [Accepted: 09/22/2023] [Indexed: 09/28/2023] Open
Abstract
Rectal cancer is a deadly disease typically treated using neoadjuvant chemoradiotherapy followed by total mesorectal excision surgery. To reduce the occurrence of mesorectal excision surgery for patients whose tumors regress from the neoadjuvant therapy alone, conventional imaging, such as computed tomography (CT) or magnetic resonance imaging (MRI), is used to assess tumor response to neoadjuvant therapy before surgery. In this work, we hypothesize that shear wave elastography offers valuable insights into tumor response to short-course radiation therapy (SCRT)-information that could help distinguish radiation-responsive from radiation-non-responsive tumors and shed light on changes in the tumor microenvironment that may affect radiation response. To test this hypothesis, we performed elastographic imaging on murine rectal tumors (n = 32) on days 6, 10, 12, 16, 18, 20, 23, and 25 post-tumor cell injection. The study revealed that radiation-responsive and non-radiation-responsive tumors had different mechanical properties. Specifically, radiation-non-responsive tumors showed significantly higher shear wave speed SWS (p < 0.01) than radiation-responsive tumors 11 days after SCRT. Furthermore, there was a significant difference in shear wave attenuation (SWA) (p < 0.01) in radiation-non-responsive tumors 16 days after SCRT compared to SWA measured just one day after SCRT. These results demonstrate the potential of shear wave elastography to provide valuable insights into tumor response to SCRT and aid in exploring the underlying biology that drives tumors' responses to radiation.
Collapse
Affiliation(s)
- Reem Mislati
- Department of Electrical and Computer Engineering, University of Rochester, Rochester, NY, USA
| | - Taylor P Uccello
- Department of Microbiology and Immunology, University of Rochester, Rochester, NY, USA
| | - Zixi Lin
- School of Engineering, Brown University, Providence, RI, USA
| | - Katia T Iliza
- Department of Biomedical Engineering, University of Rochester, Rochester, NY, USA
| | | | - Scott A Gerber
- Department of Surgery, University of Rochester, Rochester, NY, USA
| | - Marvin M Doyley
- Department of Electrical and Computer Engineering, University of Rochester, Rochester, NY, USA.
| |
Collapse
|
2
|
Wang Y, Ono S, Johnson MP, Larson NB, Lynch T, Urban MW. Evaluating Variability of Commercial Liver Fibrosis Elastography Phantoms. ULTRASOUND IN MEDICINE & BIOLOGY 2023; 49:1018-1030. [PMID: 36690519 DOI: 10.1016/j.ultrasmedbio.2022.12.017] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/01/2022] [Revised: 12/12/2022] [Accepted: 12/22/2022] [Indexed: 06/17/2023]
Abstract
OBJECTIVE Liver fibrosis has been found to increase the mechanical stiffness of the liver. To mimic different stages of liver fibrosis, commercially available phantoms (Model 039, CIRS, Inc.) have been produced for clinical quality assurance and research purposes. The purpose of this study was to investigate the mechanical property variability of the phantoms in two lots of CIRS Model 039 phantoms. METHODS Each lot consisted of phantoms of four stiffness types, and there were 8-10 phantoms of each type. Shear wave elastography measurements were conducted on each phantom at 10 different angles. Group velocity measurements and phase velocity curves were calculated for every SWE acquisition. Multilevel functional principal component analysis (MFPCA) was performed on phase velocity data, which decomposes each phase velocity curve into the sum of eigenfunctions of two levels. The variance of the component scores of levels 1 and 2 were used to represent inter-phantom and intra-phantom variability, respectively. The 95% confidence intervals of phase velocity in a phantom type were calculated to reflect curve variability. DISCUSSION The standard deviations of the group velocity for phantoms of any type were less than 0.04 and 0.02 m/s for lots 1 and 2, respectively. For both lots, in every type, the phase velocity curves of most individual phantoms fall within the 95% confidence interval. CONCLUSION MFPCA is an effective tool for analyzing the inter- and intra-phantom variability of phase velocity curves. Given the known variability of a fully tested lot, estimation of the variability of a new lot can be performed with a reduced number of phantoms tested.
Collapse
Affiliation(s)
- Yuqi Wang
- Department of Radiology, Mayo Clinic, Rochester, MN, USA.
| | | | - Matthew P Johnson
- Department of Quantitative Health Sciences, Division of Clinical Trials and Biostatistics, Mayo Clinic, Rochester, MN, USA
| | - Nicholas B Larson
- Department of Quantitative Health Sciences, Division of Clinical Trials and Biostatistics, Mayo Clinic, Rochester, MN, USA
| | | | - Matthew W Urban
- Department of Radiology, Mayo Clinic, Rochester, MN, USA; Department of Physiology and Biomedical Engineering, Mayo Clinic, Rochester, MN, USA
| |
Collapse
|
3
|
Xiao Y, Jin J, Yuan Y, Zhao Y, Li D. A New Estimation Scheme for Improving the Performance of Shear Wave Elasticity Imaging. ULTRASOUND IN MEDICINE & BIOLOGY 2023; 49:289-308. [PMID: 36283938 DOI: 10.1016/j.ultrasmedbio.2022.09.003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/04/2022] [Revised: 08/09/2022] [Accepted: 09/04/2022] [Indexed: 06/16/2023]
Abstract
Shear wave velocity (SWV) reconstruction based on time-of-flight (TOF) is widely adopted to realize shear wave elasticity imaging (SWEI). It typically breaks down the reconstruction of a SWV image into many kernels and treats them independently. We hypothesized that information exchange among kernels improves the performance of SWEI. Therefore, we propose the approach of iterative re-weighted least squares based on inter-kernel communication (IKC-IRLS). We also hypothesized that time-to-peak (TTP) is superior to cross-correlation (CC) in visualizing small targets because TTP uses higher shear wave frequencies than CC. To examine the hypotheses, IKC-IRLS was combined with TTP data and compared with four established methods. The five methods were tested by imaging several small-size stiff targets (2.5, 4.0 and 6.4 mm in diameter) using different kernel sizes in the simulation and real experiments. The results indicate that the IKC-IRLS approach can mitigate speckle noise and is robust to TTP outliers. Consequently, the proposed method achieves the highest contrast-to-noise ratio and the lowest mean absolute percentage error of target in almost all tested cases.
Collapse
Affiliation(s)
- Yang Xiao
- Department of Control Science and Engineering, Harbin Institute of Technology, Harbin, Heilongjiang Province, China
| | - Jing Jin
- Department of Control Science and Engineering, Harbin Institute of Technology, Harbin, Heilongjiang Province, China.
| | - Yu Yuan
- Department of Control Science and Engineering, Harbin Institute of Technology, Harbin, Heilongjiang Province, China
| | - Yue Zhao
- Department of Control Science and Engineering, Harbin Institute of Technology, Harbin, Heilongjiang Province, China
| | - Dandan Li
- Department of Control Science and Engineering, Harbin Institute of Technology, Harbin, Heilongjiang Province, China
| |
Collapse
|
4
|
Goswami S, Ahmed R, Feng F, Khan S, Doyley MM, McAleavey SA. Imaging the Local Nonlinear Viscoelastic Properties of Soft Tissues: Initial Validation and Expected Benefits. IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL 2022; 69:975-987. [PMID: 34986096 PMCID: PMC9815723 DOI: 10.1109/tuffc.2021.3140203] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Imaging tissue mechanical properties has shown promise in noninvasive assessment of numerous pathologies. Researchers have successfully measured many linear tissue mechanical properties in laboratory and clinical settings. Currently, multiple complex mechanical effects such as frequency-dependence, anisotropy, and nonlinearity are being investigated separately. However, a concurrent assessment of these complex effects may enable more complete characterization of tissue biomechanics and offer improved diagnostic sensitivity. In this work, we report for the first time a method to map the frequency-dependent nonlinear parameters of soft tissues on a local scale. We recently developed a nonlinear elastography model that combines strain measurements from arbitrary tissue compression with radiation-force-based broadband shear wave speed (WS) measurements. Here, we extended this model to incorporate local measurements of frequency-dependent shear modulus. This combined approach provides a local frequency-dependent nonlinear parameter that can be obtained with arbitrary, clinically realizable tissue compression. Initial assessments using simulations and phantoms validate the accuracy of this approach. We also observed improved contrast in nonlinearity parameter at higher frequencies. Results from ex-vivo liver experiments show 32, 25, 34, and 38 dB higher contrast in elastograms than traditional linear elasticity, elastic nonlinearity, viscosity, and strain imaging methods, respectively. A lesion, artificially created by injection of glutaraldehyde into a liver specimen, showed a 59% increase in the frequency-dependent nonlinear parameter and a 17% increase in contrast ratio.
Collapse
|
5
|
Romero SE, Naemi R, Flores G, Allan D, Ormachea J, Gutierrez E, Casado FL, Castaneda B. Plantar Soft Tissue Characterization Using Reverberant Shear Wave Elastography: A Proof-of-Concept Study. ULTRASOUND IN MEDICINE & BIOLOGY 2022; 48:35-46. [PMID: 34702642 DOI: 10.1016/j.ultrasmedbio.2021.09.011] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/23/2021] [Revised: 09/10/2021] [Accepted: 09/13/2021] [Indexed: 06/13/2023]
Abstract
Plantar soft tissue stiffness provides relevant information on biomechanical characteristics of the foot. Therefore, appropriate monitoring of foot elasticity could be useful for diagnosis, treatment or health care of people with complex pathologies such as a diabetic foot. In this work, the reliability of reverberant shear wave elastography (RSWE) applied to plantar soft tissue was investigated. Shear wave speed (SWS) measurements were estimated at the plantar soft tissue at the first metatarsal head, the third metatarsal head and the heel from both feet in five healthy volunteers. Experiments were repeated for a test-retest analysis with and without the use of gel pad using a mechanical excitation frequency range between 400 and 600 Hz. Statistical analysis was performed to evaluate the reliability of the SWS estimations. In addition, the results were compared against those obtained with a commercially available shear wave-based elastography technique, supersonic imaging (SSI). The results indicate a low coefficient of variation for test-retest experiments with gel pad (median: 5.59%) and without gel pad (median: 5.83%). Additionally, the values of the SWS measurements increase at higher frequencies (median values: 2.11 m/s at 400 Hz, 2.16 m/s at 450 Hz, 2.24 m/s at 500 Hz, 2.21 m/s at 550 Hz and 2.31 m/s at 600 Hz), consistent with previous reports at lower frequencies. The SWSs at the plantar soft tissue at the first metatarsal head, third metatarsal head and heel were found be significantly (p<0.05) different, with median values of 2.42, 2.16 and 2.03 m/s, respectively which indicates the ability of the method to differentiate between shear wave speeds at different anatomical locations. The results indicated better elastographic signal-to-noise ratios with RSWE compared to SSI because of the artifacts presented in the SWS generation. These preliminary results indicate that the RSWE approach can be used to estimate the plantar soft tissue elasticity, which may have great potential to better evaluate changes in biomechanical characteristics of the foot.
Collapse
Affiliation(s)
- Stefano E Romero
- Laboratorio de Imagenes Medicas, Pontificia Universidad Catolica del Peru, San Miguel, Lima, Peru.
| | - Roozbeh Naemi
- Centre for Biomechanics and Rehabilitation Technologies, School of Health Science and Wellbeing, Staffordshire University, Stoke-on-Trent, United Kingdom
| | - Gilmer Flores
- Laboratorio de Imagenes Medicas, Pontificia Universidad Catolica del Peru, San Miguel, Lima, Peru
| | - David Allan
- Centre for Biomechanics and Rehabilitation Technologies, School of Health Science and Wellbeing, Staffordshire University, Stoke-on-Trent, United Kingdom
| | - Juvenal Ormachea
- Laboratorio de Imagenes Medicas, Pontificia Universidad Catolica del Peru, San Miguel, Lima, Peru; Department of Electrical and Computer Engineering, University of Rochester, Rochester, New York, USA
| | - Evelyn Gutierrez
- Laboratorio de Imagenes Medicas, Pontificia Universidad Catolica del Peru, San Miguel, Lima, Peru
| | - Fanny L Casado
- Instituto de Ciencias Omicas y Biotecnologia Aplicada, Pontificia Universidad Catolica del Peru, San Miguel, Lima, Peru
| | - Benjamin Castaneda
- Laboratorio de Imagenes Medicas, Pontificia Universidad Catolica del Peru, San Miguel, Lima, Peru
| |
Collapse
|
6
|
Shear Wave Velocity Estimation Using the Real-Time Curve Tracing Method in Ultrasound Elastography. APPLIED SCIENCES-BASEL 2021. [DOI: 10.3390/app11052095] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
The estimation of shear wave velocity is very important in ultrasonic shear wave elasticity imaging (SWEI). Since the stability and accuracy of ultrasonic testing equipment have been greatly improved, in order to further improve the accuracy of shear wave velocity estimation and increase the quality of shear wave elasticity maps, we propose a novel real-time curve tracing (RTCT) technique to accurately reconstruct the motion trace of shear wave fronts. Based on the curve fitting of each frame shear wave, the propagation velocity of two-dimensional shear waves can be estimated. In this paper, shear wave velocity estimation and shear wave image reconstruction are implemented for homogeneous regions and stiff spherical inclusion regions with different elasticity, respectively. The experimental result shows that the proposed shear wave velocity estimation method based on the real-time curve tracing method has advantages in accuracy and anti-noise performance. Moreover, by eliminating artifacts of shear wave videos, the velocity map acquired can restore the shape of inclusions better. The real-time curve tracing method can provide a new idea for the estimation of shear wave velocity and elastic parameters.
Collapse
|
7
|
Ahmed S, Kamal U, Hasan MK. DSWE-Net: A deep learning approach for shear wave elastography and lesion segmentation using single push acoustic radiation force. ULTRASONICS 2021; 110:106283. [PMID: 33166787 DOI: 10.1016/j.ultras.2020.106283] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/24/2019] [Revised: 10/02/2020] [Accepted: 10/10/2020] [Indexed: 06/11/2023]
Abstract
Ultrasound-based non-invasive elasticity imaging modalities have received significant consideration for tissue characterization over the last few years. Though substantial advances have been made, the conventional Shear Wave Elastography (SWE) methods still suffer from poor image quality in regions far from the push location, particularly those which rely on single focused ultrasound push beam to generate shear waves. In this study, we propose DSWE-Net, a novel deep learning-based approach that is able to construct Young's modulus maps from ultrasonically tracked tissue velocity data resulting from a single acoustic radiation force (ARF) push. The proposed network employs a 3D convolutional encoder, followed by a recurrent block consisting of several Convolutional Long Short-Term Memory (ConvLSTM) layers to extract high-level spatio-temporal features from different time-frames of the input velocity data. Finally, a pair of coupled 2D convolutional decoder blocks reconstructs the modulus image and additionally performs inclusion segmentation by generating a binary mask. We also propose a multi-task learning loss function for end-to-end training of the network with 1260 data samples obtained from a simulation environment which include both bi-level and multi-level phantom structures. The performance of the proposed network is evaluated on 140 synthetic test data and the results are compared both qualitatively and quantitatively with that of the current state of the art method, Local Phase Velocity Based Imaging (LPVI). With an average SSIM of 0.90, RMSE of 0.10 and 20.69 dB PSNR, DSWE-Net performs much better on the imaging task compared to LPVI. Our method also achieves an average IoU score of 0.81 for the segmentation task which makes it suitable for localizing inclusions as well. In this initial study, we also show that our method gains an overall improvement of 0.09 in SSIM, 4.81 dB in PSNR, 2.02 dB in CNR, and 0.09 in RMSE over LPVI on a completely unseen set of CIRS tissue mimicking phantom data. This proves its better generalization capability and shows its potential for use in real-world clinical practice.
Collapse
Affiliation(s)
- Shahed Ahmed
- Department of Electrical and Electronic Engineering, Bangladesh University of Engineering and Technology, Dhaka 1205, Bangladesh
| | - Uday Kamal
- Department of Electrical and Electronic Engineering, Bangladesh University of Engineering and Technology, Dhaka 1205, Bangladesh
| | - Md Kamrul Hasan
- Department of Electrical and Electronic Engineering, Bangladesh University of Engineering and Technology, Dhaka 1205, Bangladesh.
| |
Collapse
|
8
|
Baek J, Ahmed R, Ye J, Gerber SA, Parker KJ, Doyley MM. H-scan, Shear Wave and Bioluminescent Assessment of the Progression of Pancreatic Cancer Metastases in the Liver. ULTRASOUND IN MEDICINE & BIOLOGY 2020; 46:3369-3378. [PMID: 32907773 PMCID: PMC9066934 DOI: 10.1016/j.ultrasmedbio.2020.08.006] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/24/2020] [Revised: 08/04/2020] [Accepted: 08/05/2020] [Indexed: 05/24/2023]
Abstract
The non-invasive quantification of tumor burden and the response to therapies remain an important objective for imaging modalities. To characterize the performance of two newly optimized ultrasound-based analyses, we applied shear wave and H-scan scattering analyses to repeated trans-abdominal ultrasound scans of a murine model of metastatic pancreatic cancer. In addition, bioluminescence measurements were obtained as an alternative reference. The tumor metastases grow aggressively and result in death at approximately 4 wk if untreated, but longer for those treated with chemotherapy. We found that our three imaging methods (shear wave speed, H-scan, bioluminescence) trended toward increasing output measures with time during tumor growth, and these measures were delayed for the group receiving chemotherapy. The relative sensitivity of H-scan tracked closely with bioluminescence measurements, particularly in the early to mid-stages of tumor growth. The correlation between H-scan and bioluminescence was found to be strong, with a Spearman's rank correlation coefficient greater than 0.7 across the entire series. These preliminary results suggest that non-invasive ultrasound imaging analyses are capable of tracking the response of tumor models to therapeutic agents.
Collapse
Affiliation(s)
- Jihye Baek
- Department of Electrical and Computer Engineering, University of Rochester, Rochester, New York, USA
| | - Rifat Ahmed
- Department of Electrical and Computer Engineering, University of Rochester, Rochester, New York, USA
| | - Jian Ye
- Department of Surgery, University of Rochester Medical Center, Rochester, New York, USA
| | - Scott A Gerber
- Department of Surgery, University of Rochester Medical Center, Rochester, New York, USA
| | - Kevin J Parker
- Department of Electrical and Computer Engineering, University of Rochester, Rochester, New York, USA.
| | - Marvin M Doyley
- Department of Electrical and Computer Engineering, University of Rochester, Rochester, New York, USA
| |
Collapse
|
9
|
Goswami S, Ahmed R, Khan S, Doyley MM, McAleavey SA. Shear Induced Non-Linear Elasticity Imaging: Elastography for Compound Deformations. IEEE TRANSACTIONS ON MEDICAL IMAGING 2020; 39:3559-3570. [PMID: 32746104 PMCID: PMC8527856 DOI: 10.1109/tmi.2020.2999439] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/12/2023]
Abstract
The goal of non-linear ultrasound elastography is to characterize tissue mechanical properties under finite deformations. Existing methods produce high contrast non-linear elastograms under conditions of pure uni-axial compression, but exhibit bias errors of 10-50% when the applied deformation deviates from the uni-axial condition. Since freehand transducer motion generally does not produce pure uniaxial compression, a motion-agnostic non-linearity estimator is desirable for clinical translation. Here we derive an expression for measurement of the Non-Linear Shear Modulus (NLSM) of tissue subject to combined shear and axial deformations. This method gives consistent nonlinear elasticity estimates irrespective of the type of applied deformation, with a reduced bias in NLSM values to 6-13%. The method combines quasi-static strain imaging with Single-Track Location-Shear Wave Elastography (STL-SWEI) to generate local estimates of axial strain, shear strain, and Shear Wave Speed (SWS). These local values were registered and non-linear elastograms reconstructed with a novel nonlinear shear modulus estimation scheme for general deformations. Results on tissue mimicking phantoms were validated with mechanical measurements and multiphysics simulations for all deformation types with an error in NLSM of 6-13%. Quantitative performance metrics with the new compound-motion tracking strategy reveal a 10-15 dB improvement in Signal-to-Noise Ratio (SNR) for simple shear versus pure compressive deformation for NLSM elastograms of homogeneous phantoms. Similarly, the Contrast-to-Noise Ratio (CNR) of NLSM elastograms of inclusion phantoms improved by 25-30% for simple shear over pure uni-axial compression. Our results show that high fidelity NLSM estimates may be obtained at ~30% lower strain under conditions of shear deformation as opposed axial compression. The reduction in strain required could reduce sonographer effort and improve scan safety.
Collapse
|
10
|
Ahmed R, Doyley MM. Parallel Receive Beamforming Improves the Performance of Focused Transmit-Based Single-Track Location Shear Wave Elastography. IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL 2020; 67:2057-2068. [PMID: 32746171 PMCID: PMC7590368 DOI: 10.1109/tuffc.2020.2998979] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/05/2023]
Abstract
Single-track location shear wave elastography (STL-SWEI) is robust against speckle-induced noise in shear wave speed (SWS) estimates; however, it is not immune to other incoherent sources of noise (such as electronic noise) that increases the variance in SWS estimates. Although estimation averaging enabled by parallel receive beamforming adequately suppresses these noise sources, these beamforming techniques often rely on broad transmit beams (plane or diverging). While broad beam approaches, such as plane-wave imaging, are becoming ubiquitous in research ultrasound systems, clinical systems usually employ focused transmit beams due to compatibility with hardware beamforming and deeper penetration. Consequently, improving the noise robustness of focused transmit-based STL-SWEI may enable easier translation to clinical scenarios. In this article, we experimentally evaluated the performance of parallel beamforming for STL-SWEI using fixed or multiple transmit focus. By imaging tissue-mimicking phantoms, we found that parallel beamforming improved the focal zone elastographic signal-to-noise ratio (SNRe) by 40.9%. For a receive line spacing equivalent to transducer pitch, averaging estimates from three parallel lines produced peak SNRe at the focal zone (25 mm), while, at the shallower regions (< 20 mm), a larger number of parallel lines (>7) were needed. Increasing the beamforming line density by a factor of 8 increased the focal zone SNRe only by 13.2%. When SWS quantification was desirable at a fixed depth (such as within the push focal depth), using a deeper tracking focal zone enabled higher parallel line count and improved the peak SNRe by 33%. The multifocusing strategy produced a lower SNRe than the single-focus configurations. For a fixed tracking focal zone, a depth-dependent averaging based on the simulated transmit intensity adequately accounted for the transmit beamwidth. The results in this work demonstrated that STL-SWEI can be implemented using focused transmit beams with robust noise-suppression capability.
Collapse
|
11
|
Wiseman LM, Urban MW, McGough RJ. A parametric evaluation of shear wave speeds estimated with time-of-flight calculations in viscoelastic media. THE JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA 2020; 148:1349. [PMID: 33003848 PMCID: PMC7482672 DOI: 10.1121/10.0001813] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/19/2019] [Revised: 07/30/2020] [Accepted: 08/11/2020] [Indexed: 06/11/2023]
Abstract
Shear wave elasticity imaging (SWEI) uses an acoustic radiation force to generate shear waves, and then soft tissue mechanical properties are obtained by analyzing the shear wave data. In SWEI, the shear wave speed is often estimated with time-of-flight (TOF) calculations. To characterize the errors produced by TOF calculations, three-dimensional (3D) simulated shear waves are described by time-domain Green's functions for a Kelvin-Voigt model evaluated for multiple combinations of the shear elasticity and the shear viscosity. Estimated shear wave speeds are obtained from cross correlations and time-to-peak (TTP) calculations applied to shear wave particle velocities and shear wave particle displacements. The results obtained from these 3D shear wave simulations indicate that TTP calculations applied to shear wave particle displacements yield effective estimates of the shear wave speed if noise is absent, but cross correlations applied to shear wave particle displacements are more robust when the effects of noise and shear viscosity are included. The results also show that shear wave speeds estimated with TTP methods and cross correlations using shear wave particle velocities are more sensitive to increases in shear viscosity and noise, which suggests that superior estimates of the shear wave speed are obtained from noiseless or noisy shear wave particle displacements.
Collapse
Affiliation(s)
- Luke M Wiseman
- Department of Electrical and Computer Engineering, Michigan State University, East Lansing, Michigan 48824, USA
| | - Matthew W Urban
- Department of Radiology, Mayo Clinic, Rochester, Minnesota 55905, USA
| | - Robert J McGough
- Department of Electrical and Computer Engineering, Michigan State University, East Lansing, Michigan 48824, USA
| |
Collapse
|
12
|
Wiseman LM, Urban MW, McGough RJ. A parametric evaluation of shear wave speeds estimated with time-of-flight calculations in viscoelastic media. THE JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA 2020; 148:1349. [PMID: 33003848 PMCID: PMC7482672 DOI: 10.1121/10.0001813#suppl] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/19/2019] [Revised: 07/30/2020] [Accepted: 08/11/2020] [Indexed: 06/11/2023]
Abstract
Shear wave elasticity imaging (SWEI) uses an acoustic radiation force to generate shear waves, and then soft tissue mechanical properties are obtained by analyzing the shear wave data. In SWEI, the shear wave speed is often estimated with time-of-flight (TOF) calculations. To characterize the errors produced by TOF calculations, three-dimensional (3D) simulated shear waves are described by time-domain Green's functions for a Kelvin-Voigt model evaluated for multiple combinations of the shear elasticity and the shear viscosity. Estimated shear wave speeds are obtained from cross correlations and time-to-peak (TTP) calculations applied to shear wave particle velocities and shear wave particle displacements. The results obtained from these 3D shear wave simulations indicate that TTP calculations applied to shear wave particle displacements yield effective estimates of the shear wave speed if noise is absent, but cross correlations applied to shear wave particle displacements are more robust when the effects of noise and shear viscosity are included. The results also show that shear wave speeds estimated with TTP methods and cross correlations using shear wave particle velocities are more sensitive to increases in shear viscosity and noise, which suggests that superior estimates of the shear wave speed are obtained from noiseless or noisy shear wave particle displacements.
Collapse
Affiliation(s)
- Luke M Wiseman
- Department of Electrical and Computer Engineering, Michigan State University, East Lansing, Michigan 48824, USA
| | - Matthew W Urban
- Department of Radiology, Mayo Clinic, Rochester, Minnesota 55905, USA
| | - Robert J McGough
- Department of Electrical and Computer Engineering, Michigan State University, East Lansing, Michigan 48824, USA
| |
Collapse
|
13
|
Ahmed R, Ye J, Gerber SA, Linehan DC, Doyley MM. Preclinical Imaging Using Single Track Location Shear Wave Elastography: Monitoring the Progression of Murine Pancreatic Tumor Liver Metastasis In Vivo. IEEE TRANSACTIONS ON MEDICAL IMAGING 2020; 39:2426-2439. [PMID: 32012006 PMCID: PMC7329602 DOI: 10.1109/tmi.2020.2971422] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/12/2023]
Abstract
Recently, researchers have discovered the direct impact of the tumor mechanical environment on the growth, drug uptake and prognosis of tumors. While estimating the mechanical parameters (solid stress, fluid pressure, stiffness) can aid in the treatment planning and monitoring, most of these parameters cannot be quantified noninvasively. Shear wave elastography (SWE) has shown promise as a means of noninvasively measuring the stiffness of soft tissue. However, stiffness is still not a recognized imaging biomarker. While SWE has been shown to be capable of measuring tumor stiffness in humans, much important research is done in small animal preclinical models, where tumors are often too small for the resolution of traditional SWE tools. Single-track location SWE (STL-SWE) has previously been shown to overcome the fundamental resolution limit of SWE imposed by ultrasound speckle, which may make it suitable for preclinical imaging. Using STL-SWE, in this work, we demonstrate, for the first time, that the stiffness changes occurring inside metastatic murine pancreatic tumors can be monitored over long time scales (up to 9 weeks). To prevent the respiration motion from degrading the STL-SWE estimates, we developed a real-time software-based respiration gating scheme that we implemented on a Verasonics ultrasound imaging system. By imaging the liver of three healthy mice and performing correlation analysis, we confirmed that the respiration-gated STL-SWE data was free from motion corruption. By performing coregistered power-doppler imaging, we found that the local variability in liver shear wave speed (SWS) measurements increased from 5.4% to 9.9% due to blood flow. We performed a longitudinal study using a murine model of pancreatic cancer liver metastasis to assess the temporal changes (over nine weeks) in SWS in two groups: a controlled group receiving no treatment (n=8), and an experimental group (n=6) treated with Gemcitabine, a chemotherapy agent. We independently evaluated tumor burden using bioluminescence imaging (BLI). The initial and endpoint SWS measurements were statistically different (p<0.05). Additionally, when the liver SWS exceeded 2.5 ± 0.3 and 2.73 ± 0.34 m/s in untreated and treated mice, respectively, the death of the mice was imminent within approximately 10 days. The time taken for the SWS to exceed the thresholds was 17 days (on average) longer in Gemcitabine treated mice compared to the untreated ones. The survival statistics corroborated the effectiveness of Gemcitabine. Spearman correlation analysis revealed a monotonic relationship between SWE measurements (SWS) and BLI measurements (radiance) for tumors whose radiance exceeded 1×107 photons/s/cm2/sr. Longitudinal measurements on the liver of four healthy mice revealed a maximum coefficient of variation of 11.4%. The results of this investigation demonstrate that with appropriate gating, researchers can use STL-SWE for small animal imaging and perform longitudinal studies using preclinical cancer models.
Collapse
|
14
|
Wang H, Mills B, Mislati R, Ahmed R, Gerber SA, Linehan D, Doyley MM. Shear Wave Elastography Can Differentiate between Radiation-Responsive and Non-responsive Pancreatic Tumors: An ex Vivo Study with Murine Models. ULTRASOUND IN MEDICINE & BIOLOGY 2020; 46:393-404. [PMID: 31727378 PMCID: PMC7060930 DOI: 10.1016/j.ultrasmedbio.2019.10.005] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/21/2019] [Revised: 10/01/2019] [Accepted: 10/09/2019] [Indexed: 05/04/2023]
Abstract
Neither contrast-enhanced computed tomography nor magnetic resonance imaging can monitor changes in the pancreatic ductal adenocarcinoma microenvironment during therapy. We hypothesized that shear wave elastography could overcome this limitation. To test this hypothesis, we measured the shear modulus of two groups of murine pancreatic tumors (KCKO, n = 30; PAN02, n = 30) treated with stereotactic body radiation therapy (SBRT). The mean shear modulus of KCKO tumors was 7.651 kPa higher than that of PAN02 tumors (p < 0.001). SBRT reduced the shear modulus in KCKO tumors by 8.914 kPa (p < 0.001). No significant difference in the shear modulus of SBRT-treated PAN02 tumors was observed. Additionally, necrotic and collagen densities were reduced only in the SBRT-treated KCKO tumors. Shear modulus was dependent on collagen distribution and histological texture parameters (i.e., entropy and fractional dimension). Shear wave elastography imaging differentiates between SBRT-responsive (KCKO) and non-responsive (PAN02) tumors.
Collapse
Affiliation(s)
- Hexuan Wang
- Department of Electrical & Computer Engineering, University of Rochester, Rochester, New York, USA
| | - Bradley Mills
- Department of Surgery, University of Rochester Medical Center, Rochester, New York, USA
| | - Reem Mislati
- Department of Electrical & Computer Engineering, University of Rochester, Rochester, New York, USA
| | - Rifat Ahmed
- Department of Electrical & Computer Engineering, University of Rochester, Rochester, New York, USA
| | - Scott A Gerber
- Department of Surgery, University of Rochester Medical Center, Rochester, New York, USA
| | - David Linehan
- Department of Surgery, University of Rochester Medical Center, Rochester, New York, USA
| | - Marvin M Doyley
- Department of Electrical & Computer Engineering, University of Rochester, Rochester, New York, USA.
| |
Collapse
|
15
|
Goswami S, Ahmed R, Doyley MM, McAleavey SA. Nonlinear Shear Modulus Estimation With Bi-Axial Motion Registered Local Strain. IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL 2019; 66:1292-1303. [PMID: 31150340 PMCID: PMC6684490 DOI: 10.1109/tuffc.2019.2919600] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/10/2023]
Abstract
Nonlinear elasticity imaging provides additional information about tissue behavior that is potentially diagnostic and avoids errors inherent in applying a linear elastic model to tissue under large strains. Nonlinear elasticity imaging is challenging to perform due to the large deformations required to obtain sufficient tissue strain to elicit nonlinear behavior. This work uses a method of axial and lateral displacement tracking to estimate local axial strain with simultaneous measurement of shear modulus at multiple compression levels. By following the change in apparent shear modulus and the stress deduced from the strain maps, we are able to accurately quantify nonlinear shear modulus (NLSM). We have validated our technique with a mechanical NLSM measurement system. Our results demonstrate that 2-D tracking provides more consistent NLSM estimates than those obtained by 1-D (axial) tracking alone, especially where lateral motion is significant. The elastographic contrast-to-noise ratio in heterogeneous phantoms was 12.5%-60% higher using our method than that of 1-D tracking. Our method is less susceptible to mechanical variations, with deviations in mean elastic values of 2%-4% versus 5%-37% for 1-D tracking.
Collapse
|
16
|
Ahmed R, Doyley MM. Distributing Synthetic Focusing Over Multiple Push-Detect Events Enhances Shear Wave Elasticity Imaging Performance. IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL 2019; 66:1170-1184. [PMID: 30990427 PMCID: PMC6701192 DOI: 10.1109/tuffc.2019.2911036] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/04/2023]
Abstract
Plane wave (PW) imaging is a commonly used method for tracking waves during shear wave elasticity imaging (SWEI), but its unfocused transmission beam reduces tracking accuracy and precision. Coherent compounding minimizes this problem, but SWEI's stringent frame rate requirement and the coarse pitch of most clinical transducers limit its effectiveness. Synthetic aperture imaging (SAI) is an alternate ultrasound imaging approach with a much tighter focus than PW imaging, but its lower transmission power has deterred researchers from using SAI in SWEI. Hadamard-encoded multielement SAI can overcome this limitation. However, only a limited number of subapertures (3-5) can be transmitted in a single push-detect event. We have developed methods to distribute more subapertures or more compounding angles over multiple push-detect events. In this paper, we report the results of experiments conducted on phantoms to assess SWEI's performance when using Hadamard-encoded distributed-multielement synthetic aperture (HDMSA) imaging or distributed plane wave compounding (DPWC) to track shear waves. Tracking shear waves with HDMSA improved the elastographic signal-to-noise ratio (SNRe) by 61.6%-89.5% depending on the phantom employed. Similarly, DPWC tracking improved SNRe by 56.2%-93.3% for the same group of phantoms. Compared to focused ultrasound tracking (at the focus), SNRe improved by 28.6% and 32.5% when tracking shear waves with HDMSA and DPWC, respectively. Long acquisitions could introduce decoding errors that decrease the performance when performing HDMSA tracking within the clinical setting. Nevertheless, the results of studies performed on the bicep muscle of three healthy volunteers demonstrate that for stationary organs, tracking shear waves with HDMSA yielded repeatable elastograms that offer better elastographic performance than those produced with current tracking methods.
Collapse
|
17
|
Wang H, Mislati R, Ahmed R, Vincent P, Nwabunwanne SF, Gunn JR, Pogue BW, Doyley MM. Elastography Can Map the Local Inverse Relationship between Shear Modulus and Drug Delivery within the Pancreatic Ductal Adenocarcinoma Microenvironment. Clin Cancer Res 2018; 25:2136-2143. [PMID: 30352906 DOI: 10.1158/1078-0432.ccr-18-2684] [Citation(s) in RCA: 38] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2018] [Revised: 10/05/2018] [Accepted: 10/19/2018] [Indexed: 12/18/2022]
Abstract
PURPOSE High tissue pressure prevents chemotherapeutics from reaching the core of pancreatic tumors. Therefore, targeted therapies have been developed to reduce this pressure. While point probes have shown the effectiveness of these pressure-reducing therapies via single-location estimates, ultrasound elastography is now widely available as an imaging technique to provide real-time spatial maps of shear modulus (tissue stiffness). However, the relationship between shear modulus and the underlying tumor microenvironmental causes of high tissue pressure has not been investigated. In this work, elastography was used to investigate how shear modulus influences drug delivery in situ, and how it correlates with collagen density, hyaluronic acid content, and patent vessel density-features of the tumor microenvironment known to influence tissue pressure. EXPERIMENTAL DESIGN Intravenous injection of verteporfin, an approved human fluorescent drug, was used in two pancreatic cancer xenograft models [AsPC-1 (n = 25) and BxPC-3 (n = 25)]. RESULTS Fluorescence intensity was higher in AsPC-1 tumors than in BxPC-3 tumors (P < 0.0001). Comparing drug uptake images and shear wave elastographic images with histologic images revealed that: (i) drug delivery and shear modulus were inversely related, (ii) shear modulus increased linearly with increasing collagen density, and (iii) shear modulus was marginally correlated with the local assessment of hyaluronic acid content. CONCLUSIONS These results demonstrate that elastography could guide targeted therapy and/or identify patients with highly elevated tissue pressure.See related commentary by Nia et al., p. 2024.
Collapse
Affiliation(s)
- Hexuan Wang
- Department of Electrical and Computer Engineering, University of Rochester, Rochester, New York
| | - Reem Mislati
- Department of Electrical and Computer Engineering, University of Rochester, Rochester, New York
| | - Rifat Ahmed
- Department of Electrical and Computer Engineering, University of Rochester, Rochester, New York
| | - Phuong Vincent
- Thayer School of Engineering, Dartmouth College, Hanover, New Hampshire
| | | | - Jason R Gunn
- Thayer School of Engineering, Dartmouth College, Hanover, New Hampshire
| | - Brian W Pogue
- Thayer School of Engineering, Dartmouth College, Hanover, New Hampshire
| | - Marvin M Doyley
- Department of Electrical and Computer Engineering, University of Rochester, Rochester, New York.
| |
Collapse
|