Harmonic motion imaging of human breast masses: an in vivo clinical feasibility

Harmonic Motion Imaging-Guided Focused Ultrasound Ablation: Comparison of Three Focused Ultrasound Interference Filtering Methods

OBJECTIVE

Harmonic motion imaging (HMI) is an acoustic radiation force-based elasticity imaging technique, which can be used to monitor changes in tissue mechanical properties caused by focused ultrasound (FUS)-induced thermal ablation. In conventional HMI, the amplitude-modulated FUS sequence and imaging pulse are transmitted simultaneously. With this method, the high-amplitude FUS signal must be separated from the imaging data for tissue displacement estimation. Frequency domain notch filtering and bandpass filtering were previously used to filter FUS interference from imaging data. However, FUS interference becomes more pronounced at the increased intensities and durations necessary for thermal ablation, rendering frequency domain filtering challenging. In this study, three methods for FUS interference filtering during HMI-guided FUS ablation (HMIgFUS) were compared. The methods were notch filtering; interleaved HMIgFUS, with interleaved FUS and imaging pulses; and FUS-net, a convolutional neural network-based U-Net autoencoder developed by our group for FUS interference filtering in radiofrequency data. FUS-net was applied here for the first time for the purpose of ablation displacement monitoring.

METHODS

The three filtering methods were tested during 20 HMIgFUS acquisitions in an ex vivo canine liver using a range of peak positive pressures from 11 to 18 MPa and durations ranging from 60 to 180 seconds. The B-mode mean squared errors (MSEs) and displacement amplitude contrast-to-noise ratios (CNRs) of the three methods were calculated and compared.

RESULTS

The interleaved method for HMIgFUS was found to be significantly robust in avoiding FUS interference in all tested cases for FUS ablation monitoring, especially cases with high FUS pressure and long durations, as opposed to traditional notch filtering and FUS-net filtering. CNRs obtained from displacement amplitude maps in the interleaved data set were significantly higher in all cases than those obtained from the notch filtered and FUS-net data sets. There was not a significant trend in displacement CNR between the FUS-net and notch filtered data sets. However, B-mode MSE was found to be significantly higher when comparing the FUS-net and interleaved data sets as opposed to the notch filtered and interleaved data sets, suggesting further potential of FUS-net as an FUS interference filtering method.

CONCLUSION

These findings indicate the robustness of interleaved HMIgFUS in avoiding FUS interference during HMIgFUS monitoring and the advantages, limitations and future potential of FUS-net and notch filtering.

Optimization of the Tracking Beam Sequence in Harmonic Motion Imaging

Harmonic motion imaging (HMI) is an ultrasound elastography technique that estimates the viscoelastic properties of tissues by inducing localized oscillatory motion using focused ultrasound (FUS). The resulting displacement, assumed to be inversely proportional to the tissue local stiffness, is estimated using an imaging array based on RF speckle tracking. In conventional HMI, this is accomplished with plane-wave (PW) imaging, which inherently suffers from low lateral resolution. Coherent PW compounding (PWC) leverages spatial and temporal resolution using synthetic focusing in transmit. In this study, we introduced focused imaging with parallel tracking in HMI and compared parallel tracking of various transmit F-numbers (F/2.6, 3, 4, and 5) qualitatively and quantitatively with PW and PWC imaging at various compounded angle ranges (6°, 12°, and 18°). An in silico model of a 56-kPa spherical inclusion (diameter: 3.6 mm) embedded in a 5.3-kPa background and a 5.3-kPa elastic phantom with cylindrical inclusions (Young's moduli: 22-56 kPa, diameters: 2.0-8.6 mm) were imaged to assess different tracking beam sequences. Speckle biasing in displacement estimation associated with parallel tracking was also investigated and concluded to be negligible in HMI. Parallel tracking in receive (Rx) resulted in 2%-7% and 8%-12% increase compared to PW imaging ( ) in HMI contrast and contrast-to-noise ratio in silico and phantoms. Focused imaging with parallel tracking in Rx was concluded to be most robust among PW and PWC imaging for displacement estimation, and its preclinical feasibility was demonstrated in postsurgical human cancerous breast tissue specimens and in vivo murine models of breast cancer.

The impact of amplitude modulation frequency in harmonic motion imaging on inclusion characterization

OBJECTIVE

Ultrasound elasticity imaging techniques aim to provide a non-invasive characterization of tissue mechanical properties to detect pathological changes and monitor disease progression. Harmonic motion imaging (HMI) is an ultrasound-based elasticity imaging technique that utilizes an oscillatory acoustic radiation force to induce localized displacements and estimate relative tissue stiffness. Previous studies have applied a low amplitude modulation (AM) frequency of 25 or 50 Hz in HMI to assess the mechanical properties of different tissue types. In this study, we investigate the dependence of AM frequency in HMI and whether the frequency can be adjusted based on the size and mechanical properties of the underlying medium for enhanced image contrast and inclusion detection.

METHODS

A tissue-mimicking phantom with embedded inclusions at different sizes and stiffnesses was imaged within a range of AM frequencies from 25 to 250 Hz at 25-Hz step size.

DISCUSSION

The AM frequency at which the maximum contrast and CNR are achieved depends on the size and stiffness of the inclusions. A general trend shows that contrast and CNR peak at higher frequencies for smaller inclusions. In addition, for some inclusions with the same size but different stiffnesses, the optimized AM frequency increases with the stiffness of the inclusion. Nevertheless, there is a shift between the frequencies at which the contrast peaks and those with maximum CNR. Finally, in agreement with the phantom findings, imaging an ex-vivo human specimen with a 2.7-cm breast tumor at a range of AM frequencies showed that the highest contrast and CNR are achieved at the AM frequency of 50 Hz.

CONCLUSION

These findings indicate that the AM frequency can be optimized in different applications of HMI, especially in the clinic, for improved detection and characterization of tumors with different geometries and mechanical properties.

An Efficient and Multi-Focal Focused Ultrasound Technique for Harmonic Motion Imaging

Harmonic motion imaging (HMI) is an ultrasound-based elasticity imaging technique that utilizes oscillatory acoustic radiation force to estimate the mechanical properties of tissues, as well as monitor high-intensity focused ultrasound (HIFU) treatment. Conventionally, in HMI, a focused ultrasound (FUS) transducer generates oscillatory tissue displacements, and an imaging transducer acquires channel data for displacement estimation, with each transducer being driven with a separate system. The fixed position of the FUS focal spot requires mechanical translation of the transducers, which can be a time-consuming and challenging procedure. In this study, we developed and characterized a new HMI system with a multi-element FUS transducer with the capability of electronic focal steering of ±5 mm and ±2 mm from the geometric focus in the axial and lateral directions, respectively. A pulse sequence was developed to drive both the FUS and imaging transducers using a single ultrasound data acquisition (DAQ) system. The setup was validated on a tissue-mimicking phantom with embedded inclusions. Integrating beam steering with the mechanical translation of the transducers resulted in a consistent high contrast-to-noise ratio (CNR) for the inclusions with Young's moduli of 22 and 44 kPa within a 5-kPa background while the data acquisition speed is increased by 4.5-5.2-fold compared to the case when only mechanical movements were applied. The feasibility of simultaneous generation of multiple foci and tracking the induced displacements is demonstrated in phantoms for applications where imaging or treatment of a larger region is needed. Moreover, preliminary feasibility is shown in a human subject with a breast tumor, where the mean HMI displacement within the tumor was about 4 times lower than that within perilesional tissues. The proposed HMI system facilitates data acquisition in terms of flexibility and speed and can be potentially used in the clinic for breast cancer imaging and treatment.

Imaging of Single Transducer-Harmonic Motion Imaging-Derived Displacements at Several Oscillation Frequencies Simultaneously

Mapping of mechanical properties, dependent on the frequency of motion, is relevant in diagnosis, monitoring treatment response, or intra-operative surgical resection planning. While shear wave speeds at different frequencies have been described elsewhere, the effect of frequency on the "on-axis" acoustic radiation force (ARF)-induced displacement has not been previously investigated. Instead of generating single transducer-harmonic motion imaging (ST-HMI)-derived peak-to-peak displacement (P2PD) image at a particular frequency, a novel multi-frequency excitation pulse is proposed to generate P2PD images at 100-1000 Hz simultaneously. The performance of the proposed excitation pulse is compared with the ARFI by imaging 16 different inclusions (Young's moduli of 6, 9, 36, 70 kPa and diameters of 1.6, 2.5, 6.5, and 10.4 mm) embedded in an 18 kPa background. Depending on inclusion size and stiffness, the maximum CNR and contrast were achieved at different frequencies and were always higher than ARFI. The frequency, at which maximum CNR and contrast were achieved, increased with stiffness for fixed inclusion's size and decreased with size for fixed stiffness. In vivo feasibility is tested by imaging a 4T1 breast cancer mouse tumor on Day 6, 12, and 19 post-injection of tumor cells. Similar to phantoms, the CNR of ST-HMI images was higher than ARFI and increased with frequency for the tumor on Day 6. Besides, P2PD at 100-1000 Hz indicated that the tumor became stiffer with respect to the neighboring non-cancerous tissue over time. These results indicate the importance of using a multi-frequency excitation pulse to simultaneously generate displacement at multiple frequencies to better delineate inclusions or tumors.

Feasibility of Harmonic Motion Imaging Using a Single Transducer: In Vivo Imaging of Breast Cancer in a Mouse Model and Human Subjects

Harmonic motion imaging (HMI) interrogates the mechanical properties of tissues by simultaneously generating and tracking harmonic oscillation using focused ultrasound and imaging transducers, respectively. Instead of using two transducers, the objective of this work is to develop a single transducer HMI (ST-HMI) to both generate and track harmonic motion at "on-axis" to the force for facilitating data acquisition. In ST-HMI, the amplitude-modulated force was generated by modulating excitation pulse duration and tracking of motion was performed by transmitting tracking pulses interleaved between excitation pulses. The feasibility of ST-HMI was performed by imaging two elastic phantoms with three inclusions (N = 6) and comparing it with acoustic radiation force impulse (ARFI) imaging, in vivo longitudinal monitoring of 4T1, orthotropic breast cancer mice (N = 4), and patients (N = 3) with breast masses in vivo. Six inclusions with Young's moduli of 8, 10, 15, 20, 40, and 60 kPa were embedded in a 5 kPa background. The ST-HMI-derived peak-to-peak displacement (P2PD) successfully detected all inclusions with [Formula: see text] of the linear regression between the P2PD ratio of background to inclusion versus Young's moduli ratio of inclusion to background. The contrasts of 10 and 15 kPa inclusions were higher in ST-HMI than ARFI-derived images. In the mouse study, the median P2PD ratio of tumor to non-cancerous tissues was 3.0, 5.1, 6.1, and 7.7 at 1, 2, 3, and 4 weeks post-injection of the tumor cells, respectively. In the clinical study, ST-HMI detected breast masses including fibroadenoma, pseudo angiomatous stromal hyperplasia, and invasive ductal carcinoma with a P2PD ratio of 1.37, 1.61, and 1.78, respectively. These results indicate that ST-HMI can assess the mechanical properties of tissues via generation and tracking of harmonic motion "on-axis" to the ARF. This study is the first step towards translating ST-HMI in clinics.

An analytical model of full-field displacement and strain induced by amplitude-modulated focused ultrasound in harmonic motion imaging

The majority of disease processes involves changes in the micro-structure of the affected tissue, which can translate to changes in the mechanical properties of the corresponding tissue. Harmonic motion imaging (HMI) is an elasticity imaging technique that allows the study of the mechanical parameters of tissue by detecting the tissue response by a harmonic motion field, which is generated by oscillatory acoustic radiation force (ARF). HMI has been demonstrated in tumor detection and characterization as well as monitoring of ablation procedures. In this study, an analytical HMI model is demonstrated and compared with a finite element model (FEM), allowing rapid and accurate computation of the displacement, strain, and shear wave velocity (SWV) at any location in a homogeneous linear elastic material. Average absolute differences between the analytical model and the FEM were respectively 1.2 % for the displacements and 0.5 % for the strains for 41940 force voxels at 0.22 seconds per displacement evaluation. A convergence study showed that the average difference could be further decreased to 1.0 % and 0.15 % for the displacements and strains, respectively, if force resolution is increased. SWV fields, as calculated with the FEM and the analytical model, have regional differences in velocities up to 0.57 m/s with an average absolute difference of 0.11±0.07 m/s, primarily due to imperfections in the non-reflecting FEM boundary conditions. The apparent SWV differed from the commonly used plane-wave approximation by up to 1.2 m/s due to near and intermediate field effects. Maximum displacement amplitudes for a model with an inclusion stabilize within 10 % of the homogeneous model at an inclusion radius of 10 mm while the maximum strain reacts faster, stabilizing at an inclusion radius of 3 mm. In conclusion, an analytical model for HMI stiffness estimation is presented in this paper. The analytical model has advantages over FEM as the full-field displacements do not need to be calculated to evaluate the model at a single measurement point. This advantage, together with the computational speed, makes the analytical model useful for real-time imaging applications. However, the analytical model was found to have restrictive assumptions on tissue homogeneity and infinite dimensions, while the FEM approaches were shown adaptable to variable geometry and non-homogeneous properties.

Synchronous temperature variation monitoring during ultrasound imaging and/or treatment pulse application: a phantom study

Ultrasound attenuation through soft tissues can produce an acoustic radiation force (ARF) and heating. The ARF-induced displacements and temperature evaluations can reveal tissue properties and provide insights into focused ultrasound (FUS) bio-effects. In this study, we describe an interleaving pulse sequence tested in a tissue-mimicking phantom that alternates FUS and plane-wave imaging pulses at a 1 kHz frame rate. The FUS is amplitude modulated, enabling the simultaneous evaluation of tissue-mimicking phantom displacement using harmonic motion imaging (HMI) and temperature rise using thermal strain imaging (TSI). The parameters were varied with a spatial peak temporal average acoustic intensity (I _spta ) ranging from 1.5 to 311 W.cm^-2, mechanical index (MI) from 0.43 to 4.0, and total energy (E) from 0.24 to 83 J.cm^-2. The HMI and TSI processing could estimate displacement and temperature independently for temperatures below 1.80°C and displacements up to ~117 μm (I _spta <311 W.cm^-2, MI<4.0, and E<83 J.cm^-2) indicated by a steady-state tissue-mimicking phantom displacement throughout the sonication and a comparable temperature estimation with simulations in the absence of tissue-mimicking phantom motion. The TSI estimations presented a mean error of ±0.03°C versus thermocouple estimations with a mean error of ±0.24°C. The results presented herein indicate that HMI can operate at diagnostic-temperature levels (i.e., <1°C) even when exceeding diagnostic acoustic intensity levels (720 mW.cm^-2 < I _spta < 207 W.cm^-2). In addition, the combined HMI and TSI can potentially be used for simultaneous evaluation of safety during tissue elasticity imaging as well as FUS mechanism involved in novel ultrasound applications such as ultrasound neuromodulation and tumor ablation.