Bubble-Induced Color Doppler Feedback for Histotripsy Tissue Fractionation

Acoustic Cavitation Emissions Predict Near-complete/complete Histotripsy Treatment in Soft Tissues

OBJECTIVE

Histotripsy is a non-invasive acoustic ablation technique that leverages cavitation to impart mechanical damage to a viscoelastic medium, such as tissue. Although histotripsy bubbles and lesions can be imaged with a variety of modalities, reliable methods to predict tissue disruption across different tissue-types remain to be determined.

APPROACH

Several ex-vivo bovine tissues were ablated by intrinsic threshold histotripsy over a range of pulse-per-location acoustic doses. Acoustic Cavitation Emission (ACE) signals were captured following every other therapeutic pulse using transmit-receive capable histotripsy arrays. Final bubble lifespan, lifespan-slope, and percent-reduction were calculated and correlated against histologic necrosis score (0-5: 0=0% necrosis, 5=>95% necrosis) and residual structure score (0-4: 0=none present, 4=intact) to evaluate the ability of features from ACE-signals to predict histotripsy-induced damage. Further, optimal ACE-feature thresholds were determined for binary evaluation of whether a necrosis score equal or greater than 4 had been reached.

RESULTS

Measured lifespans increased and lifespan-slopes decreased with pulses per location (ppl) and eventually plateaued in all tissue types, in similar trends to those previously observed in tissue phantoms. Necrosis score increased and residual structure decreased with increasing acoustic dose. Bubble lifespan-slope and percent-reduction correlated well with necrosis score. Thresholds able to predict the necrosis score of 4 or greater in brain, liver, and kidney were calculated with high sensitivity/specificity (>80%). The necrosis score of 4 and 5 is expected to correspond to near-complete/complete ablation by histological evaluation.

CONCLUSION

Features measured from ACE-signals, particularly the lifespan-slope and percent reduction, were used to predict near-complete/complete ablation of large-volume histotripsy treatments in ex vivo bovine liver, kidney, and brain tissues with good accuracy. Tissue heterogeneities were observed to impact the histotripsy damage and corresponding ACE-signals, and thus the predication accuracy.

Overview of Therapeutic Ultrasound Applications and Safety Considerations: 2024 Update

A 2012 review of therapeutic ultrasound was published to educate researchers and physicians on potential applications and concerns for unintended bioeffects (doi: 10.7863/jum.2012.31.4.623). This review serves as an update to the parent article, highlighting advances in therapeutic ultrasound over the past 12 years. In addition to general mechanisms for bioeffects produced by therapeutic ultrasound, current applications, and the pre-clinical and clinical stages are outlined. An overview is provided for image guidance methods to monitor and assess treatment progress. Finally, other topics relevant for the translation of therapeutic ultrasound are discussed, including computational modeling, tissue-mimicking phantoms, and quality assurance protocols.

A comparative study of histotripsy parameters for the treatment of fibrotic ex-vivo human benign prostatic hyperplasia tissue

Histotripsy is a noninvasive focused ultrasound therapy that mechanically fractionates tissue to create well-defined lesions. In a previous clinical pilot trial to treat benign prostatic hyperplasia (BPH), histotripsy did not result in consistent objective improvements in symptoms, potentially because of the fibrotic and mechanically tough nature of this tissue. In this study, we aimed to identify the dosage required to homogenize BPH tissue by different histotripsy modalities, including boiling histotripsy (BH) and cavitation histotripsy (CH). A method for histotripsy lesion quantification via entropy (HLQE) analysis was developed and utilized to quantify lesion area of the respective treatments. These data were correlated to changes in mechanical stiffness measured by ultrasound shear-wave elastography before and after treatment with each parameter set and dose. Time points corresponding to histologically observed complete lesions were qualitatively evaluated and quantitatively measured. For the BH treatment, complete lesions occurred with >  = 30 s treatment time, with a corresponding maximum reduction in stiffness of -90.9 ± 7.2(s.d.)%. High pulse repetition frequency (PRF) CH achieved a similar reduction to that of BH at 288 s (-91.6 ± 6.0(s.d.)%), and low-PRF CH achieved a (-82.1 ± 5.1(s.d.)%) reduction in stiffness at dose >  = 144 s. Receiver operating characteristic curve analysis showed that a >  ~ 75% reduction in stiffness positively correlated with complete lesions observed histologically, and can provide an alternative metric to track treatment progression.

A Comparative Study of Histotripsy Parameters for the Treatment of Fibrotic ex-vivo Human Benign Prostatic Hyperplasia Tissue

Histotripsy is a noninvasive focused ultrasound therapy that mechanically fractionates tissue to create well-defined lesions. In a previous clinical pilot trial to treat benign prostatic hyperplasia (BPH), histotripsy did not result in consistent objective improvements in symptoms, potentially because of the fibrotic and mechanically tough nature of this tissue. In this study, we aimed to identify the dosage required to homogenize BPH tissue by different histotripsy modalities, including boiling histotripsy (BH) and cavitation histotripsy (CH). A method for histotripsy lesion quantification via entropy (HLQE) analysis was developed and utilized to quantify lesion area of the respective treatments. These data were correlated to changes in mechanical stiffness measured by ultrasound shear-wave elastography before and after treatment with each parameter set and dose. Time points corresponding to histologically observed complete lesions were qualitatively evaluated and quantitatively measured. For the BH treatment, complete lesions occurred with >=30s treatment time, with a corresponding maximum reduction in stiffness of -90.9±7.2(s.d.)%. High pulse repetition frequency (PRF) CH achieved a similar reduction to that of BH at 288s (-91.6±6.0(s.d.)%), and low-PRF CH achieved a (-82.1±5.1(s.d.)%) reduction in stiffness at dose >=144s. Receiver operating characteristic curve analysis showed that a >~75% reduction in stiffness positively correlated with complete lesions observed histologically, and can provide an alternative metric to track treatment progression.

Contrast-Enhanced Imaging of Histotripsy Bubble Clouds Using Chirp-Coded Excitation and Volterra Filtering

Histotripsy is a focused ultrasound therapy that ablates tissue via bubble cloud activity. Real-time ultrasound image guidance is used to ensure safe and effective treatment. Plane-wave imaging enables tracking of histotripsy bubble clouds at a high frame rate but lacks adequate contrast. Furthermore, bubble cloud hyperechogenicity is reduced in abdominal targets, making the development of contrast-specific sequences for deep-seated targets an active area of research. Chirp-coded subharmonic imaging was reported previously to enhance histotripsy bubble cloud detection by a modest 4-6 dB compared to the conventional sequence. Incorporating additional steps into the signal processing pipeline could enhance bubble cloud detection and tracking. In this study, we evaluated the feasibility of combining chirp-coded subharmonic imaging with Volterra filtering for enhancing bubble cloud detection in vitro. Chirped imaging pulses were used to track bubble clouds generated in scattering phantoms at a 1-kHz frame rate. Fundamental and subharmonic matched filters were applied to the received radio frequency signals, followed by a tuned Volterra filter to extract bubble-specific signatures. For subharmonic imaging, the application of the quadratic Volterra filter improved the contrast-to-tissue ratio from 5.18 ± 1.29 to 10.90 ± 3.76 dB, relative to the application of the subharmonic matched filter. These findings demonstrate the utility of the Volterra filter for histotripsy image guidance.

The histotripsy spectrum: differences and similarities in techniques and instrumentation

Since its inception about two decades ago, histotripsy - a non-thermal mechanical tissue ablation technique - has evolved into a spectrum of methods, each with distinct potentiating physical mechanisms: intrinsic threshold histotripsy, shock-scattering histotripsy, hybrid histotripsy, and boiling histotripsy. All methods utilize short, high-amplitude pulses of focused ultrasound delivered at a low duty cycle, and all involve excitation of violent bubble activity and acoustic streaming at the focus to fractionate tissue down to the subcellular level. The main differences are in pulse duration, which spans microseconds to milliseconds, and ultrasound waveform shape and corresponding peak acoustic pressures required to achieve the desired type of bubble activity. In addition, most types of histotripsy rely on the presence of high-amplitude shocks that develop in the pressure profile at the focus due to nonlinear propagation effects. Those requirements, in turn, dictate aspects of the instrument design, both in terms of driving electronics, transducer dimensions and intensity limitations at surface, shape (primarily, the F-number) and frequency. The combination of the optimized instrumentation and the bio-effects from bubble activity and streaming on different tissues, lead to target clinical applications for each histotripsy method. Here, the differences and similarities in the physical mechanisms and resulting bioeffects of each method are reviewed and tied to optimal instrumentation and clinical applications.

Quantitative Assessment of Boiling Histotripsy Progression Based on Color Doppler Measurements

Boiling histotripsy (BH) is a mechanical tissue liquefaction method that uses sequences of millisecond-long high intensity focused ultrasound (HIFU) pulses with shock fronts. The BH treatment generates bubbles that move within the sonicated volume due to acoustic radiation force. Since the velocity of the bubbles and tissue debris is expected to depend on the lesion size and liquefaction completeness, it could provide a quantitative metric of the treatment progression. In this study, the motion of bubble remnants and tissue debris immediately following BH pulses was investigated using high-pulse repetition frequency (PRF) plane-wave color Doppler ultrasound in ex vivo myocardium tissue. A 256-element 1.5 MHz spiral HIFU array with a coaxially integrated ultrasound imaging probe (ATL P4-2) produced 10 ms BH pulses to form volumetric lesions with electronic beam steering. Prior to performing volumetric BH treatments, the motion of intact myocardium tissue and anticoagulated bovine blood following isolated BH pulses was assessed as two limiting cases. In the liquid blood the velocity of BH-induced streaming at the focus reached over 200 cm/s, whereas the intact tissue was observed to move toward the HIFU array consistent with elastic rebound of tissue. Over the course of volumetric BH treatments tissue motion at the focus locations was dependent on the axial size of the forming lesion relative to the corresponding size of the HIFU focal area. For axially small lesions, the maximum velocity after the BH pulse was directed toward the HIFU transducer and monotonically increased over time from about 20-100 cm/s as liquefaction progressed, then saturated when tissue was fully liquefied. For larger lesions obtained by merging multiple smaller lesions in the axial direction, the high-speed streaming away from the HIFU transducer was observed at the point of full liquefaction. Based on these observations, the maximum directional velocity and its location along the HIFU propagation axis were proposed and evaluated as candidate metrics of BH treatment completeness.

Histotripsy for the Treatment of Uterine Leiomyomas: A Feasibility Study in Ex Vivo Uterine Fibroids

Uterine fibroids (leiomyomas), the most common benign tumors in women of reproductive age, are a frequent cause of abnormal vaginal bleeding and other reproductive complaints among women. This study investigates the feasibility of using histotripsy, a non-invasive, non-thermal focused ultrasound ablation method, to ablate uterine fibroids. Human fibroid samples (n = 16) were harvested after hysterectomy or myomectomy procedures at Carilion Memorial Hospital. Histotripsy was applied to ex vivo fibroids in two sets of experiments using a 700-kHz clinical transducer to apply multicycle histotripsy pulses and a prototype 500-kHz transducer to apply single-cycle histotripsy pulses. Ultrasound imaging was used for real-time treatment monitoring, and post-treatment ablation was quantified histologically using hematoxylin and eosin and Masson trichrome stains. Results revealed that multicycle histotripsy generated diffuse cavitation in targeted fibroids, with minimal cellular ablative changes after treatment with 2000 pulses/point. Single-cycle pulsing generated well-confined bubble clouds with evidence of early coagulative necrosis on histological assessment in samples treated with 2000 pulses/point, near-complete ablation in samples treated with 4000 pulses/point and complete tissue destruction in samples treated with 10,000 pulses/point. This study illustrates that histotripsy is capable of fibroid ablation under certain pulsing parameters and warrants further investigation as an improved non-invasive ablation method for the treatment of leiomyomas.

Estimating the mechanical energy of histotripsy bubble clouds with high frame rate imaging

Mechanical ablation with the focused ultrasound therapy histotripsy relies on the generation and action of bubble clouds. Despite its critical role for ablation, quantitative metrics of bubble activity to gauge treatment outcomes are still lacking. Here, plane wave imaging was used to track the dissolution of bubble clouds following initiation with the histotripsy pulse. Information about the rate of change in pixel intensity was coupled with an analytic diffusion model to estimate bubble size. Accuracy of the hybrid measurement/model was assessed by comparing the predicted and measured dissolution time of the bubble cloud. Good agreement was found between predictions and measurements of bubble cloud dissolution times in agarose phantoms and murine subcutaneous SCC VII tumors. The analytic diffusion model was extended to compute the maximum bubble size as well as energy imparted to the tissue due to bubble expansion. Regions within tumors predicted to have undergone strong bubble expansion were collocated with ablation. Further, the dissolution time was found to correlate with acoustic emissions generated by the bubble cloud during histotripsy insonation. Overall, these results indicate a combination of modeling and high frame rate imaging may provide means to quantify mechanical energy imparted to the tissue due to bubble expansion for histotripsy.

Doppler Passive Acoustic Mapping

In therapeutic ultrasound using microbubbles, it is essential to drive the microbubbles into the correct type of activity and the correct location to produce the desired biological response. Although passive acoustic mapping (PAM) is capable of locating where microbubble activities are generated, it is well known that microbubbles rapidly move within the ultrasound beam. We propose a technique that can image microbubble movement by estimating their velocities within the focal volume. Microbubbles embedded within a wall-less channel of a tissue-mimicking material were sonicated using 1-MHz focused ultrasound. The acoustic emissions generated by the microbubbles were captured with a linear array (L7-4). PAM with robust Capon beamforming was used to localize the microbubble acoustic emissions. We spectrally analyzed the time trace of each position and isolated the higher harmonics. Microbubble velocity maps were constructed from the position-dependent Doppler shifts at different time points during sonication. Microbubbles moved primarily away from the transducer at velocities on the order of 1 m/s due to primary acoustic radiation forces, producing a time-dependent velocity distribution. We detected microbubble motion both away and toward the receiving array, revealing the influence of acoustic radiation forces and fluid motion due to the ultrasound exposure. High-speed optical images confirmed the acoustically measured microbubble velocities. Doppler PAM enables passive estimation of microbubble motion and may be useful in therapeutic applications, such as drug delivery across the blood-brain barrier, sonoporation, sonothrombolysis, and drug release.

Sonobactericide: An Emerging Treatment Strategy for Bacterial Infections

Ultrasound has been developed as both a diagnostic tool and a potent promoter of beneficial bio-effects for the treatment of chronic bacterial infections. Bacterial infections, especially those involving biofilm on implants, indwelling catheters and heart valves, affect millions of people each year, and many deaths occur as a consequence. Exposure of microbubbles or droplets to ultrasound can directly affect bacteria and enhance the efficacy of antibiotics or other therapeutics, which we have termed sonobactericide. This review summarizes investigations that have provided evidence for ultrasound-activated microbubble or droplet treatment of bacteria and biofilm. In particular, we review the types of bacteria and therapeutics used for treatment and the in vitro and pre-clinical experimental setups employed in sonobactericide research. Mechanisms for ultrasound enhancement of sonobactericide, with a special emphasis on acoustic cavitation and radiation force, are reviewed, and the potential for clinical translation is discussed.

For Whom the Bubble Grows: Physical Principles of Bubble Nucleation and Dynamics in Histotripsy Ultrasound Therapy

Histotripsy is a focused ultrasound therapy for non-invasive tissue ablation. Unlike thermally ablative forms of therapeutic ultrasound, histotripsy relies on the mechanical action of bubble clouds for tissue destruction. Although acoustic bubble activity is often characterized as chaotic, the short-duration histotripsy pulses produce a unique and consistent type of cavitation for tissue destruction. In this review, the action of histotripsy-induced bubbles is discussed. Sources of bubble nuclei are reviewed, and bubble activity over the course of single and multiple pulses is outlined. Recent innovations in terms of novel acoustic excitations, exogenous nuclei for targeted ablation and histotripsy-enhanced drug delivery and image guidance metrics are discussed. Finally, gaps in knowledge of the histotripsy process are highlighted, along with suggested means to expedite widespread clinical utilization of histotripsy.

MRI-guided transurethral insonation of silica-shell phase-shift emulsions in the prostate with an advanced navigation platform

PURPOSE

In this study, the efficacy of transurethral prostate ablation in the presence of silica-shell ultrasound-triggered phase-shift emulsions (sUPEs) doped with MR contrast was evaluated. The influence of sUPEs on MR imaging assessment of the ablation zone was also investigated.

METHODS

sUPEs were doped with a magnetic resonance (MR) contrast agent, Gd2 O3 , to assess ultrasound transition. Injections of saline (sham), saline and sUPEs alone, and saline and sUPEs with Optison microbubbles were performed under guidance of a prototype interventional MRI navigation platform in a healthy canine prostate. Treatment arms were evaluated for differences in lesion size, T1  contrast, and temperature. In addition, non-perfused areas (NPAs) on dynamic contrast-enhanced (DCE) MRI, 55°C isotherms, and areas of 240 cumulative equivalent minutes at 43°C (CEM43 ) dose or greater computed from MR thermometry were measured and correlated with ablated areas indicated by histology.

RESULTS

For treatment arms including sUPEs, the computed correlation coefficients between the histological ablation zone and the NPA, 55°C isotherm, and 240 CEM43 area ranged from 0.96-0.99, 0.98-0.99, and 0.91-0.99, respectively. In the absence of sUPEs, the computed correlation coefficients between the histological ablation zone and the NPA, 55°C isotherm, and 240 CEM43 area were 0.69, 0.54, and 0.50, respectively. Across all treatment arms, the areas of thermal tissue damage and NPAs were not significantly different (P = 0.47). Areas denoted by 55°C isotherms and 240 CEM43 dose boundaries were significantly larger than the areas of thermal damage, again for all treatment arms (P = 0.009 and 0.003, respectively). No significant differences in lesion size, T1 contrast, or temperature were observed between any of the treatment arms (P > 0.0167). Lesions exhibiting thermal fixation on histological analysis were present in six of nine insonations involving sUPE injections and one of five insonations involving saline sham injections. Significantly larger areas (P = 0.002), higher temperatures (P = 0.004), and more frequent ring patterns of restricted diffusion on ex vivo diffusion-weighted imaging (P = 0.005) were apparent in lesions with thermal fixation.

CONCLUSIONS

T1 contrast suggesting sUPE transition was not evident in sUPE treatment arms. The use of MR imaging metrics to predict prostate ablation was not diminished by the presence of sUPEs. Lesions generated in the presence of sUPEs exhibited more frequent thermal fixation, though there were no significant changes in the ablation areas when comparing arms with and without sUPEs. Thermal fixation corresponded to some qualitative imaging features.

Using the cavitation collapse time to indicate the extent of histotripsy-induced tissue fractionation

Histotripsy is an ultrasonic tissue ablation method based on acoustic cavitation. It has been shown that cavitation dynamics change depending on the mechanical properties of the host medium. During histotripsy treatment, the target-tissue is gradually fractionated and eventually liquefied to acellular homogenate. In this study, the change in the collapse time (t _col) of the cavitation bubble cloud over the course of histotripsy treatment is investigated as an indicator for progression of the tissue fractionation process throughout treatment. A 500 kHz histotripsy transducer is used to generate single-location lesions within tissue-mimicking agar phantoms of varying stiffness levels as well as ex vivo bovine liver samples. Cavitation collapse signals are acquired with broadband hydrophones, and cavitation is imaged optically using a high-speed camera in transparent tissue-mimicking phantoms. The high-speed-camera-acquired measurements of t _col validate the acoustic hydrophone measurements. Increases in t _col are observed both with decreasing phantom stiffness and throughout histotripsy treatment with increasing number of pulses applied. The increasing trend of t _col throughout the histotripsy treatment correlates well with the progression of lesion formation generated in tissue-mimicking phantoms (R ²  =  0.87). Finally, the increasing trend of t _col over the histotripsy treatment is validated in ex vivo bovine liver.

Bubble-Induced Color Doppler Feedback Correlates with Histotripsy-Induced Destruction of Structural Components in Liver Tissue

Bubble-induced color Doppler (BCD) is a histotripsy-therapy monitoring technique that uses Doppler ultrasound to track the motion of residual cavitation nuclei that persist after the collapse of the histotripsy bubble cloud. In this study, BCD is used to monitor tissue fractionation during histotripsy tissue therapy, and the BCD signal is correlated with the destruction of structural and non-structural components identified histologically to further understand how BCD monitors the extent of treatment. A 500-kHz, 112-element phased histotripsy array is used to generate approximately 6- × 6- × 7-mm lesions within ex vivo bovine liver tissue by scanning more than 219 locations with 30-1000 pulses per location. A 128-element L7-4 imaging probe is used to acquire BCD signals during all treatments. The BCD signal is then quantitatively analyzed using the time-to-peak rebound velocity (tprv) metric. Using the Pearson correlation coefficient, the tprv is compared with histologic analytics of lesions generated by various numbers of pulses using a significance level of 0.001. Histologic analytics in this study include viable cell count, reticulin-stained type III collagen area and trichrome-stained type I collagen area. It is found that the tprv metric has a statistically significant correlation with the change in reticulin-stained type III collagen area with a Pearson correlation coefficient of -0.94 (p <0.001), indicating that changes in BCD are more likely because of destruction of the structural components of tissue.

Post Hoc Analysis of Passive Cavitation Imaging for Classification of Histotripsy-Induced Liquefaction in Vitro

Histotripsy utilizes focused ultrasound to generate bubble clouds for transcutaneous tissue liquefaction. Bubble activity maps are under development to provide image guidance and monitor treatment progress. The aim of this paper was to investigate the feasibility of using plane wave B-mode and passive cavitation images to be used as binary classifiers of histotripsy-induced liquefaction. Prostate tissue phantoms were exposed to histotripsy pulses over a range of pulse durations (5- ) and peak negative pressures (12-23 MPa). Acoustic emissions were recorded during the insonation and beamformed to form passive cavitation images. Plane wave B-mode images were acquired following the insonation to detect the hyperechoic bubble cloud. Phantom samples were sectioned and stained to delineate the liquefaction zone. Correlation between passive cavitation and plane wave B-mode images and the liquefaction zone was assessed using receiver operating characteristic (ROC) curve analysis. Liquefaction of the phantom was observed for all the insonation conditions. The area under the ROC (0.94 versus 0.82), accuracy (0.90 versus 0.83), and sensitivity (0.81 versus 0.49) was greater for passive cavitation images relative to B-mode images ( ) along the azimuth of the liquefaction zone. The specificity was greater than 0.9 for both imaging modalities. These results demonstrate a stronger correlation between histotripsy-induced liquefaction and passive cavitation imaging compared with the plane wave B-mode imaging, albeit with limited passive cavitation image range resolution.

Non-invasive, Rapid Ablation of Tissue Volume Using Histotripsy

Histotripsy is a non-invasive, non-thermal ablation technique that uses high-amplitude, focused ultrasound pulses to fractionate tissue via acoustic cavitation. The goal of this study was to illustrate the potential of histotripsy with electronic focal steering to achieve rapid ablation of a tissue volume at a rate matching or exceeding those of current clinical techniques (∼1-2 mL/min). Treatment parameters were established in tissue-mimicking phantoms and applied to ex vivo tissue. Six-microsecond pulses were delivered by a 250-kHz array. The focus was electrically steered to 1000 locations at a pulse repetition frequency of 200 Hz (0.12% duty cycle). Magnetic resonance imaging and histology of the treated tissue revealed a distinct region of necrosis in all samples. Mean lesion volume was 35.6 ± 4.3 mL, generated at 0.9-3.3 mL/min, a speed faster than that of any current ablation method for a large volume. These results suggest that histotripsy has the potential to achieve non-invasive, rapid, homogeneous ablation of a tissue volume.

Effect of pulse repetition frequency of high-intensity focused ultrasound on in vitro thrombolysis

Vascular occlusion by the thrombi is the main reason for ischemic stroke and deep vein thrombosis. High-intensity focused ultrasound (HIFU) and histotripsy or microtripsy pulses can effectively dissolve the blood clot with no use of thrombolytic agent and ultrasound contrast agent (microbubbles). In this study, HIFU bursts at the same duty cycle (2%) but varied pulse repetition frequency (PRF) from 1Hz to 1000Hz were delivered to in vitro porcine blood clot for 30s. Thrombolysis efficiency initially increases slightly with the PRF, 86.4±10.3%, 89.9±11.9, and 92.9±12.8% at the PRF of 1Hz, 10Hz, and 100Hz, respectively, without significant difference (p>0.05), but then drops dramatically to 37.9±6.9% at the PRF of 1000Hz (p<0.05). The particle size in the supernatant of dissolution is 547.1±129.5nm, which suggests the disruption of thrombi into the subcellular level. Thrombi motion during HIFU exposure shows violent motion and significant curling at the low PRF, rotation about its axis with occasional curling at the moderate PRF, and localized vibration at the high PRF due to the generation of acoustic radiation force and streaming. Quantitative analysis of recorded motion shows the axial displacement decreases with the PRF of delivered HIFU bursts, from 3.9±1.5mm at 1Hz to 0.7±0.4mm at 1000Hz. Bubble cavitation during HIFU exposure to the blood clot was also monitored. The increase of PRF led to the increase of inertial cavitation but the decrease of stable cavitation. In summary, the PRF of delivered HIFU bursts at the same output energy has a significant effect on the thrombi motion, bubble cavitation activities, and subsequently thrombolysis efficiencies.