Pilot Study on Boiling Histotripsy Treatment of Human Leiomyoma Ex Vivo.

OBJECTIVE: As an alternative to surgical excision and magnetic resonance-guided thermal high-intensity focused ultrasound ablation of uterine leiomyoma, this work was aimed at pilot feasibility demonstration of use of ultrasound-guided boiling histotripsy for non-invasive non-thermal fractionation of human uterine leiomyoma ex vivo. METHODS: A custom-made sector ultrasound transducer of 1.5-MHz operating frequency and nominal f-number F# = 0.75 was used to produce a volumetric lesion (two layers of 5 × 5 foci with a 1 mm step) in surgically resected human leiomyoma ex vivo. A sequence of 10 ms pulses (P+/P-/As = 157/-25/170 MPa in situ) with 1% duty cycle was delivered N = 30 times per focus under B-mode guidance. The treatment outcome was evaluated via B-mode imaging and histologically with hematoxylin and eosin and Masson's trichrome staining. RESULTS: The treatment was successfully performed in less than 30 min and resulted in formation of a rectangular lesion visualized on B-mode images during the sonication as an echogenic region, which sustained for about 10 min post-treatment. Histology revealed loss of cellular structure, necrotic debris and globules of degenerated collagen in the target volume surrounded by injured smooth muscle cells. CONCLUSION: The pilot experiment described here indicates that boiling histotripsy is feasible for non-invasive mechanical disintegration of human uterine leiomyoma ex vivo under B-mode guidance, encouraging further investigation and optimization of this potential clinical application of boiling histotripsy.

Fully nonlinear three-dimensional modeling of parametric interactions in the field of a dual-frequency acoustic array.

An algorithm is developed for fully nonlinear three-dimensional (3D) simulation of a difference-frequency acoustic beam resulting from the interaction of two high-intensity pump waves. Simulations are performed in the frequency domain based on the Khokhlov-Zabolotskaya-Kuznetsov equation. A spectrum filtering method is used to enable accurate solutions for the difference-frequency fields in strongly nonlinear beams and with a high downshift frequency ratio using only dozens of spectral components retained in the algorithm. As an example, the dual-frequency operation of an underwater multi-element ellipsoidal array is considered, and numerical solutions describing parametric interactions in the array field are analyzed. It is shown that difference-frequency beams are more symmetric in transverse directions compared with the pump beams. The most efficient parametric generation of difference-frequency beams corresponded to close and beyond shock-forming conditions. Axial pressure amplitude of the difference frequency was shown to grow first quadratically with the source pressure following the quasi-linear solution and then linearly once shocks start to develop. The percentage of the total power converted to the difference frequency from pump waves increased at high power outputs without saturation. Up to twofold increase in directivity angles of difference-frequency beams under shock-forming conditions was observed compared with quasi-linear conditions.

Elastic Properties of Aging Human Hematoma Model In Vitro and Its Susceptibility to Histotripsy Liquefaction.

OBJECTIVE: Tissue susceptibility to histotripsy disintegration has been reported to depend on its elastic properties. This work was aimed at investigation of histotripsy efficiency for liquefaction of human hematomas, depending on their stiffness and degree of retraction over time (0-10 d). METHODS: As an in vitro hematoma model, anticoagulated human blood samples (200 mL) were recalcified at different temperatures. In one set of samples, the shear modulus was measured by shear wave elastography during blood clotting at 10℃, 22℃ and 37℃, and then daily during further aging. The ultrastructure of the samples was analyzed daily with scanning electron microscopy (SEM). Another set of blood samples (50-200 mL) were recalcified at 37℃ for density and retraction measurements over aging and exposed to histotripsy at varying time points. Boiling histotripsy (2.5 ms pulses) and hybrid histotripsy (0.2 ms pulses) exposures (2 MHz, 1% dc, P+/P-/As = 182/-27/207 MPa in situ) were used to produce either individual cigar-shaped or volumetric (0.8-3 mL) lesions in samples incubated for 3 h, 5 d and 10 d. The obtained lesions were sized, then the lysate aspirated under B-mode guidance was analyzed ultrastructurally and diluted in distilled water for sizing of residual fragments. RESULTS: It was found that clotting time decreased from 113 to 25 min with the increase in blood temperature from 10℃ to 37℃. The shear modulus increased to 0.53 ± 0.17 kPa during clotting and remained constant within 8 d of incubation at 2℃. Sample volumes decreased by 57% because of retraction within 10 d. SEM revealed significant echinocytosis but unchanged ultrastructure of the fibrin meshwork. Liquefaction rate and lesion dimensions produced with the same histotripsy protocols correlated with the increase in the degree of retraction and were lower in retracted samples versus freshly clotted samples. More than 80% of residual fibrin fragments after histotripsy treatment were shorter than 150 µm; the maximum length was 208 µm, allowing for unobstructed aspiration of the lysate with most clinically used needles. CONCLUSION: The results indicate that hematoma susceptibility to histotripsy liquefaction is not entirely determined by its stiffness, and correlates with the retraction degree.

Histology-based quantification of boiling histotripsy outcomes via ResNet-18 network: Towards mechanical dose metrics.

This work was focused on the newly developed ultrasonic approach for non-invasive surgery - boiling histotripsy (BH) - recently proposed for mechanical ablation of tissues using pulsed high intensity focused ultrasound (HIFU). The BH lesion is known to depend in size and shape on exposure parameters and mechanical properties, structure and composition of tissue being treated. The aim of this work was to advance the concept of BH dose by investigating quantitative relationships between the parameters of the lesion, pulsing protocols, and targeted tissue properties. A HIFU focus of a 1.5 MHz 256-element array driven by power-enhanced Verasonics system was electronically steered along the grid within 12 × 4 × 12 mm volume to produce volumetric lesions in porcine liver (soft, with abundant collagenous structures) and bovine myocardium (stiff, homogenous cellular) ex vivo tissues with various pulsing protocols (1-10 ms pulses, 1-15 pulses per point). Quantification of the lesion size and completeness was performed through serial histological sectioning, and a computer vision approach using a combination of manual and automated detection of fully fractionated and residual tissue based on neural network ResNet-18 was developed. Histological sample fixation led to underestimation of BH ablation rate compared to the ultrasound-based estimations, and provided similar qualitative feedback as did gross inspection. This suggests that gross observation may be sufficient for qualitatively evaluating the BH treatment completeness. BH efficiency in liver tissue was shown to be insensitive to the changes in pulsing protocol within the tested parameter range, whereas in bovine myocardium the efficiency increased with either increasing pulse length or number of pulses per point or both. The results imply that one universal mechanical dose metric applicable to an arbitrary tissue type is unlikely to be established. The dose metric as a product of the BH pulse duration and the number of pulses per sonication point (BHD1) was shown to be more relevant for initial planning of fractionation of collagenous tissues. The dose metric as a number of pulses per point (BHD2) is more suitable for the treatment planning of softer targets primarily containing cellular tissue, allowing for significant acceleration of treatment using shorter pulses.