Assessment of histotripsy-induced liquefaction with diagnostic ultrasound and magnetic resonance imaging in vitro and ex vivo

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.

Enhancing boiling histotripsy efficacy with a tandem pulse sequence: Immediate hypoechoic Sonograhy and expanded lesion size

Boiling histotripsy (BH) has proven effective in noninvasively disintegrating various soft tissues through cavitation effects. Although liquefied tissue appears as hypoechoic in sonography, the transition of BH-induced bubbles from hyperechoic to hypoechoic in the focal region typically requires several minutes. To facilitate rapid clinical assessment, a tandem pulse sequence of high-power BH pulses (with an acoustic power of 1484 W, a pulse duration of 10 ms, and a pulse repetition frequency of 1 Hz) followed by low-power long ultrasound pulses (with an acoustic power of 240 W, a pulse duration of 100 ms, and a pulse repetition frequency of 1 Hz) was introduced to expedite bubble clearance, resulting in an immediate hypoechoic presentation in sonography. This method was evaluated through high-speed photography, red blood cell (RBC) phantom, and ex vivo tissue experiments. High-speed photography experiments captured the enhanced bubble clearance induced by the low-power long pulses, validating our hypothesis. In RBC phantom experiments, conventional BH sequences yielded hypoechic patterns after 4.39 ± 0.84 min, whereas the tandem pulse sequences achieved hypoechic appearance instantaneously post-treatment (p < 0.05). Moreover, the tandem pulse sequences increased the erosion area in the RBC layer by 7.8 folds, from 2.36 ± 0.88 mm2 to 18.43 ± 5.15 mm2 (p < 0.05), at the equivalent energy output. Ex vivo bovine liver experiments mirrored these findings, with hypoechoic appearance at > 10 min vs. 0 min (p < 0.05) and liquefied areas of 33.78 ± 3.28 mm2 vs. 66.52 ± 11.24 mm2 (p < 0.05), respectively. In summary, our results suggest that the strategic modulation of cavitation activities not only accelerates the immediate hypoechoic appearance in sonography but also enlarges the area of BH-induced disintegration. The tandem pulse sequence strategy presents a promising avenue for enhancing the efficacy and efficiency of BH treatment in clinical applications.

Ultrasonic Nakagami imaging for automatically positioning and identifying the treated lesion induced by histotripsy

Histotripsy has been proposed as a non-invasive surgical procedure for clinical use that liquefies the tissue into acellular debris by utilizing the mechanical mechanism of bubbles. Accurate and reliable imaging guidance is essential for successful clinical histotripsy implementation. Nakagami imaging is a promising method to evaluate the microstructural change induced by high intensity focused ultrasound. However, practically, it is difficult for the Nakagami imaging to distinguish the treated lesion induced by histotripsy from the surrounding normal biological tissues. In this study, we introduce the use of noise-assisted correlation algorithm (NCA) in Nakagami images as a solution to suppress the background normal tissue and identify the treated lesion induced by histotripsy. Experiments are conducted on fresh porcine liver ex vivo by cavitation-cloud histotripsy. Results show that the contrast-to-noise ratio between the treated lesion and surrounding tissue corresponding to the Nakagami image after NCA and original Nakagami image is 3.434 and 0.505, respectively. The optimal artificial noise level is 1-fold of the background normal tissue amplitude, and the corresponding optimal threshold of correlation coefficient should be between 0.6 and 0.8 in the application of NCA. Therefore, the use of NCA in Nakagami image can suppress the background normal tissues without affecting the information of treated lesion for an appropriate artificial noise level and threshold used in the NCA. Moreover, the Nakagami images after the application of the NCA can also be used for automatically distinguishing and measuring the tissue fractionation accurately using binarization. The proposed Nakagami images overlaid on the B-mode images can provide a promising method for positioning and visualizing the treated lesion to achieve precise histotripsy treatment.

Real-Time Passive Acoustic Mapping With Enhanced Spatial Resolution in Neuronavigation-Guided Focused Ultrasound for Blood-Brain Barrier Opening

OBJECTIVE

Passive acoustic mapping (PAM) provides the spatial information of acoustic energy emitted from microbubbles during focused ultrasound (FUS), which can be used for safety and efficacy monitoring of blood-brain barrier (BBB) opening. In our previous work with a neuronavigation-guided FUS system, only part of the cavitation signal could be monitored in real time due to the computational burden although full-burst analysis is required to detect transient and stochastic cavitation activity. In addition, the spatial resolution of PAM can be limited for a small-aperture receiving array transducer. For full-burst real-time PAM with enhanced resolution, we developed a parallel processing scheme for coherence-factor-based PAM (CF-PAM) and implemented it onto the neuronavigation-guided FUS system using a co-axial phased-array imaging transducer.

METHODS

Simulation and in-vitro human skull studies were conducted for the performance evaluation of the proposed method in terms of spatial resolution and processing speed. We also carried out real-time cavitation mapping during BBB opening in non-human primates (NHPs).

RESULTS

CF-PAM with the proposed processing scheme provided better resolution than that of traditional time-exposure-acoustics PAM with a higher processing speed than that of eigenspace-based robust Capon beamformer, which facilitated the full-burst PAM with the integration time of 10 ms at a rate of 2 Hz. In vivo feasibility of PAM with the co-axial imaging transducer was also demonstrated in two NHPs, showing the advantages of using real-time B-mode and full-burst PAM for accurate targeting and safe treatment monitoring.

SIGNIFICANCE

This full-burst PAM with enhanced resolution will facilitate the clinical translation of online cavitation monitoring for safe and efficient BBB opening.

In silico assessment of histotripsy-induced changes in catheter-directed thrombolytic delivery

Introduction: For venous thrombosis patients, catheter-directed thrombolytic therapy is the standard-of-care to recanalize the occluded vessel. Limitations with thrombolytic drugs make the development of adjuvant treatments an active area of research. One potential adjuvant is histotripsy, a focused ultrasound therapy that lyses red blood cells within thrombus via the spontaneous generation of bubbles. Histotripsy has also been shown to improve the efficacy of thrombolytic drugs, though the precise mechanism of enhancement has not been elucidated. In this study, in silico calculations were performed to determine the contribution of histotripsy-induced changes in thrombus diffusivity to alter catheter-directed therapy. Methods: An established and validated Monte Carlo calculation was used to predict the extent of histotripsy bubble activity. The distribution of thrombolytic drug was computed with a finite-difference time domain (FDTD) solution of the perfusion-diffusion equation. The FDTD calculation included changes in thrombus diffusivity based on outcomes of the Monte Carlo calculation. Fibrin degradation was determined using the known reaction rate of thrombolytic drug. Results: In the absence of histotripsy, thrombolytic delivery was restricted in close proximity to the catheter. Thrombolytic perfused throughout the focal region for calculations that included the effects of histotripsy, resulting in an increased degree of fibrinolysis. Discussion: These results were consistent with the outcomes of in vitro studies, suggesting histotripsy-induced changes in the thrombus diffusivity are a primary mechanism for enhancement of thrombolytic drugs.

An X-Ray C-Arm Guided Automatic Targeting System for Histotripsy

OBJECTIVE

Histotripsy is an emerging noninvasive, nonionizing and nonthermal focal cancer therapy that is highly precise and can create a treatment zone of virtually any size and shape. Current histotripsy systems rely on ultrasound imaging to target lesions. However, deep or isoechoic targets obstructed by bowel gas or bone can often not be treated safely using ultrasound imaging alone. This work presents an alternative x-ray C-arm based targeting approach and a fully automated robotic targeting system.

METHODS

The approach uses conventional cone beam CT (CBCT) images to localize the target lesion and 2D fluoroscopy to determine the 3D position and orientation of the histotripsy transducer relative to the C-arm. The proposed pose estimation uses a digital model and deep learning-based feature segmentation to estimate the transducer focal point relative to the CBCT coordinate system. Additionally, the integrated robotic arm was calibrated to the C-arm by estimating the transducer pose for four preprogrammed transducer orientations and positions. The calibrated system can then automatically position the transducer such that the focal point aligns with any target selected in a CBCT image.

RESULTS

The accuracy of the proposed targeting approach was evaluated in phantom studies, where the selected target location was compared to the center of the spherical ablation zones in post-treatment CBCTs. The mean and standard deviation of the Euclidean distance was 1.4 ±0.5 mm. The mean absolute error of the predicted treatment radius was 0.5 ±0.5 mm.

CONCLUSION

CBCT-based histotripsy targeting enables accurate and fully automated treatment without ultrasound guidance.

SIGNIFICANCE

The proposed approach could considerably decrease operator dependency and enable treatment of tumors not visible under ultrasound.

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.

Clot Degradation Under the Action of Histotripsy Bubble Activity and a Lytic Drug

Deep vein thrombosis is a major source of morbidity worldwide. For critical obstructions, catheter-directed thrombolytics are the frontline therapy to achieve vessel recanalization. Techniques that aid lytic therapy are under development to improve treatment efficacy and reduce procedure-related complications. Histotripsy is one such adjuvant under development that relies on focused ultrasound for in situ nucleation of bubble clouds. Prior studies have demonstrated synergistic effects for clot dissolution when histotripsy is combined with lytic therapy. The success of this combination approach is hypothesized to promote thrombolytic efficacy via two mechanisms: erythrocyte fractionation (hemolysis) and increased lytic activity (fibrinolysis). In this study, the contributions of hemolysis and fibrinolysis to clot degradation under histotripsy and a lytic were quantified with measurements of hemoglobin and D-dimer, respectively. A linear regression analysis was used to determine the relationship between hemoglobin, D-dimer, and the overall treatment efficacy (clot mass loss). A similar analysis was conducted to gauge the role of bubble activity, which was assessed with passive cavitation imaging, on hemolysis and fibrinolysis. Tabulation of these data demonstrated hemolysis and fibrinolysis contributed equally to clot mass loss. Furthermore, bubble cloud activity promoted the generation of hemoglobin and D-dimer in equal proportion. These studies indicate a multifactorial process for clot degradation under the action of histotripsy and a lytic therapy.

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.

Transcostal Histotripsy Ablation in an In Vivo Acute Hepatic Porcine Model

PURPOSE

To determine whether histotripsy can create human-scale transcostal ablations in porcine liver without causing severe thermal wall injuries along the beam path.

MATERIALS AND METHODS

Histotripsy was applied to the liver using a preclinical prototype robotic system through a transcostal window in six female swine. A 3.0 cm spherical ablation zone was prescribed. Duration of treatment (75 min) was longer than a prior subcostal treatment study (24 min, 15 s) to minimize beam path heating. Animals then underwent contrast-enhanced MRI, necropsy, and histopathology. Images and tissue were analyzed for ablation zone size, shape, completeness of necrosis, and off-target effects.

RESULTS

Ablation zones demonstrated complete necrosis with no viable tissue remaining in 6/6 animals by histopathology. Ablation zone volume was close to prescribed (13.8 ± 1.8 cm³ vs. prescribed 14.1 cm³). Edema was noted in the body wall overlying the ablation on T2 MRI in 5/5 (one animal did not receive MRI), though there was no gross or histologic evidence of injury to the chest wall at necropsy. At gross inspection, lung discoloration in the right lower lobe was present in 5/6 animals (mean size: 1 × 2 × 4 cm) with alveolar hemorrhage, preservation of blood vessels and bronchioles, and minor injuries to pneumocytes noted at histology.

CONCLUSION

Transcostal hepatic histotripsy ablation appears feasible, effective, and no severe injuries were identified in an acute porcine model when prolonged cooling time is added to minimize body wall heating.

Assessment of Collaborative Robot (Cobot)-Assisted Histotripsy for Venous Clot Ablation

OBJECTIVE

The application of bubble-based ablation with the focus ultrasound therapy histotripsy is gaining traction for the treatment of venous thrombosis, among other pathologies. For extensive clot burden, the histotripsy source must be translated to ensure uniform bubble activity throughout the vascular obstruction. The purpose of this study was to evaluate the targeting accuracy of a histotripsy system comprised of a focused source, ultrasound image guidance, and a collaborative robot (cobot) positioner. The system was designed with a primary emphasis for treating deep vein thrombosis.

METHODS

Studies to test treatment planning and targeting bubble activity with the histotripsy-cobot system were conducted in an in vitro clot model. A tissue-mimicking phantom was also targeted with the system, and the predicted and actual areas of liquefaction were compared to gauge the spatial accuracy of ablation.

RESULTS

The system provided submillimeter accuracy for both tracking along an intended path (within 0.6 mm of a model vessel) and targeting bubble activity within the venous clot model (0.7 mm from the center of the clot). Good correlation was observed between the planned and actual liquefaction locations in the tissue phantom, with an average Dice similarity coefficient of 77.8%, and average Hausdorff distance of 1.6 mm.

CONCLUSION

Cobots provide an effective means to apply histotripsy pulses over a treatment volume, with the ablation precision contingent on the quality of image guidance.

SIGNIFICANCE

Overall, these results demonstrate cobots can be used to guide histotripsy ablation for targets that extend beyond the natural focus of the transducer.

In vitro assessment of stiffness-dependent histotripsy bubble cloud activity in gel phantoms and blood clots

As a bubble-based ablative therapy, the efficacy of histotripsy has been demonstrated in healthy or acutely diseased models. Chronic conditions associated with stiff tissues may require additional bubble activity prior to histotripsy liquefaction. In this study, histotripsy pulses were generated in agarose phantoms of Young's moduli ranging from 12.3 to 142 kPa, and in vitro clot models with mild and strong platelet-activated retraction. Bubble cloud emissions were tracked with passive cavitation imaging, and the threshold acoustic power associated with phantom liquefaction was extracted with receiver operator characteristic analysis. The power of histotripsy-generated emissions and the degree of liquefaction were tabulated for both clot models. For the agarose phantoms, the acoustic power associated with liquefaction increased with Young's modulus. When grouped based on agarose concentration, only two arms displayed a significant difference in the liquefaction threshold acoustic power (22.1 kPa versus 142 kPa Young's modulus). The bubble cloud dynamics tracked with passive cavitation imaging indicated no strong changes in the bubble dynamics based on the phantom stiffness. For identical histotripsy exposure, the power of acoustic emissions and degree of clot lysis did not vary based on the clot model. Overall, these results indicate that a fixed threshold acoustic power mapped with passive cavitation imaging can be utilized for predicting histotripsy liquefaction over a wide range of tissue stiffness.