Mechanical damage induced by the appearance of rectified bubble growth in a viscoelastic medium during boiling histotripsy exposure

Histotripsy: A Method for Mechanical Tissue Ablation with Ultrasound

Histotripsy is a relatively new therapeutic ultrasound technology to mechanically liquefy tissue into subcellular debris using high-amplitude focused ultrasound pulses. In contrast to conventional high-intensity focused ultrasound thermal therapy, histotripsy has specific clinical advantages: the capacity for real-time monitoring using ultrasound imaging, diminished heat sink effects resulting in lesions with sharp margins, effective removal of the treated tissue, a tissue-selective feature to preserve crucial structures, and immunostimulation. The technology is being evaluated in small and large animal models for treating cancer, thrombosis, hematomas, abscesses, and biofilms; enhancing tumor-specific immune response; and neurological applications. Histotripsy has been recently approved by the US Food and Drug Administration to treat liver tumors, with clinical trials undertaken for benign prostatic hyperplasia and renal tumors. This review outlines the physical principles of various types of histotripsy; presents major parameters of the technology and corresponding hardware and software, imaging methods, and bioeffects; and discusses the most promising preclinical and clinical applications.

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.

Noninvasive mechanical destruction of liver tissue and tissue decellularisation by pressure-modulated shockwave histotripsy

Introduction

Boiling histotripsy (BH) is a promising High Intensity Focused Ultrasound (HIFU) technique that can be used to mechanically fractionate solid tumours at the HIFU focus noninvasively, promoting tumour immunity. Because of the occurrence of shock scattering phenomenon during BH process, the treatment accuracy of BH is, however, somewhat limited. To induce more localised and selective tissue destruction, the concept of pressure modulation has recently been proposed in our previous in vitro tissue phantom study. The aim of the present study was therefore to investigate whether this newly developed histotripsy approach termed pressure-modulated shockwave histotripsy (PSH) can be used to induce localised mechanical tissue fractionation in in vivo animal model.

Methods

In the present study, 8 Sprague Dawley rats underwent the PSH treatment and were sacrificed immediately after the exposure for morphological and histological analyses (paraffin embedded liver tissue sections were stained with H&E and MT). Partially exteriorised rat's left lateral liver lobe in vivo was exposed to a 2.0 MHz HIFU transducer with peak positive (P ₊) and negative (P _-) pressures of 89.1 MPa and -14.6 MPa, a pulse length of 5 to 34 ms, a pressure modulation time at 4 ms where P ₊ and P _- decreased to 29.9 MPa and - 9.6 MPa, a pulse repetition frequency of 1 Hz, a duty cycle of 1% and number of pulses of 1 to 20. Three lesions were produced on each animal. For comparison, the same exposure condition but no pressure modulation was also employed to create a number of lesions in the liver.

Results and Discussion

Experimental results showed that a partial mechanical destruction of liver tissue in the form of an oval in the absence of thermal damage was clearly observed at the HIFU focus after the PSH exposure. With a single pulse length of 7 ms, a PSH lesion created in the liver was measured to be a length of 1.04 ± 0.04 mm and a width of 0.87 ± 0.21 mm which was 2.37 times in length (p = 0.027) and 1.35 times in width (p = 0.1295) smaller than a lesion produced by no pressure modulation approach (e.g., BH). It was also observed that the length of a PSH lesion gradually grew towards the opposite direction to the HIFU source along the axial direction with the PSH pulse length, eventually leading to the generation of an elongated lesion in the liver. In addition, our experimental results demonstrated the feasibility of inducing partial decellularisation effect where liver tissue was partially destructed with intact extracellular matrix (i.e., intact fibrillar collagen) with the shortest PSH pulse length. Taken together, these results suggest that PSH could be used to induce a highly localised tissue fractionation with a desired degree of mechanical damage from complete tissue fractionation to tissue decellularisation through controlling the dynamics of boiling bubbles without inducing the shock scattering effect.

Sonothrombolysis with an acoustic net-assisted boiling histotripsy: A proof-of-concept study

Whilst sonothrombolysis is a promising and noninvasive ultrasound technique for treating blood clots, bleeding caused by thrombolytic agents used for dissolving clots and potential obstruction of blood flow by detached clots (i.e., embolus) are the major limitations of the current approach. In the present study, a new sonothrombolysis method is proposed for treating embolus without the use of thrombolytic drugs. Our proposed method involves (a) generating a spatially localised acoustic radiation force in a blood vessel against the blood flow to trap moving blood clots (i.e., generation of an acoustic net), (b) producing acoustic cavitation to mechanically destroy the trapped embolus, and (c) acoustically monitoring the trapping and mechanical fractionation processes. Three different ultrasound transducers with different purposes were employed in the proposed method: (1) 1-MHz dual focused ultrasound (dFUS) transducers for capturing moving blood clots, (2) a 2-MHz High Intensity Focused Ultrasound (HIFU) source for fractionating blood clots and (3) a passive acoustic emission detector with broad bandwidth (10 kHz to 20 MHz) for receiving and analysing acoustic waves scattered from a trapped embolus and acoustic cavitation. To demonstrate the feasibility of the proposed method, in vitro experiments with an optically transparent blood vessel-mimicking phantom filled with a blood mimicking fluid and a blood clot (1.2 to 5 mm in diameter) were performed with varying the dFUS and HIFU exposure conditions under various flow conditions (from 1.77 to 6.19 cm/s). A high-speed camera was used to observe the production of acoustic fields, acoustic cavitation formation and blood clot fragmentation within a blood vessel by the proposed method. Numerical simulations of acoustic and temperature fields generated under a given exposure condition were also conducted to further interpret experimental results on the proposed sonothrombolysis. Our results clearly showed that fringe pattern-like acoustic pressure fields (fringe width of 1 mm) produced in a blood vessel by the dFUS captured an embolus (1.2 to 5 mm in diameter) at the flow velocity up to 6.19 cm/s. This was likely to be due to the greater magnitude of the dFUS-induced acoustic radiation force exerted on an embolus in the opposite direction to the flow in a blood vessel than that of the drag force produced by the flow. The acoustically trapped embolus was then mechanically destructed into small pieces of debris (18 to 60 μm sized residual fragments) by the HIFU-induced strong cavitation without damaging the blood vessel walls. We also observed that acoustic emissions emitted from a blood clot captured by the dFUS and cavitation produced by the HIFU were clearly distinguished in the frequency domain. Taken together, these results can suggest that our proposed sonothrombolysis method could be used as a promising tool for treating thrombosis and embolism through capturing and destroying blood clots effectively.

Soft Tissue Aberration Correction for Histotripsy Using Acoustic Emissions From Cavitation Cloud Nucleation and Collapse

OBJECTIVE

Phase aberration from soft tissue limits the efficacy of histotripsy, a therapeutic ultrasound technique based on acoustic cavitation. Previous work has shown that the acoustic emissions from cavitation can serve as "point sources" for aberration correction (AC). This study compared the efficacy of soft tissue AC for histotripsy using acoustic cavitation emissions (ACE) from bubble cloud nucleation and collapse.

METHODS

A 750-kHz, receive-capable histotripsy array was pulsed to generate cavitation in ex vivo porcine liver through an intervening abdominal wall. Received ACE signals were used to determine the arrival time differences to the focus and compute corrective delays. Corrections from single pulses and from the median of multiple pulses were tested.

DISCUSSION

On average, ACE AC obtained 96% ± 3% of the pressure amplitude obtained by hydrophone-based correction (compared with 71% ± 5% without AC). Both nucleation- and collapse-based corrections obtained >96% of the hydrophone-corrected pressure when using medians of ≥10 pulses. When using single-pulse corrections, nucleation obtained a range of 49%-99% of the hydrophone-corrected pressure, while collapse obtained 95%-99%.

CONCLUSION

The results suggest that (i) ACE AC can recover nearly all pressure amplitude lost owing to soft tissue aberration and that (ii) the collapse signal permits robust AC using a small number of pulses.

Effects of medium viscoelasticity on bubble collapse strength of interacting polydisperse bubbles

Due to its physical and/or chemical effects, acoustic cavitation plays a crucial role in various emerging applications ranging from advanced materials to biomedicine. The cavitation bubbles usually undergo oscillatory dynamics and violent collapse within a viscoelastic medium, which are closely related to the cavitation-associated effects. However, the role of medium viscoelasticity on the cavitation dynamics has received little attention, especially for the bubble collapse strength during multi-bubble cavitation with the complex interactions between size polydisperse bubbles. In this study, modified Gilmore equations accounting for inter-bubble interactions were coupled with the Zener viscoelastic model to simulate the dynamics of multi-bubble cavitation in viscoelastic media. Results showed that the cavitation dynamics (e.g., acoustic resonant response, nonlinear oscillation behavior and bubble collapse strength) of differently-sized bubbles depend differently on the medium viscoelasticity and each bubble is affected by its neighboring bubbles to a different degree. More specifically, increasing medium viscosity drastically dampens the bubble dynamics and weakens the bubble collapse strength, while medium elasticity mainly affects the bubble resonance at which the bubble collapse strength is maximum. Differently-sized bubbles can achieve resonances and even subharmonic resonances at high driving acoustic pressures as the elasticity changes to certain values, and the resonance frequency of each bubble increases with the elasticity increasing. For the interactions between the size polydisperse bubbles, it indicated that the largest bubble generally has a dominant effect on the dynamics of smaller ones while in turn it is almost unaffected, exhibiting a pattern of destructive and constructive interactions. This study provides a valuable insight into the acoustic cavitation dynamics of multiple interacting polydisperse bubbles in viscoelastic media, which may offer a potential of controlling the medium viscoelasticity to appropriately manipulate the dynamics of multi-bubble cavitation for achieving proper cavitation effects according to the desired application.

Ultrasound Histotripsy on a Viable Perfused Whole Porcine Liver: Histological Aspects, Including 3D Reconstruction of the Histotripsy Site

Non-invasive therapeutic-focused ultrasound (US) can be used for the mechanical dissociation of tissue and is described as histotripsy. We have performed US histotripsy in viable perfused ex vivo porcine livers as a step in the development of a novel approach to hepatocyte cell transplantation. The histotripsy nidus was created with a 2 MHz single-element focused US transducer, producing 50 pulses of 10 ms duration, with peak positive and negative pressure values of P+ = 77.7 MPa and P− = –13.7 MPaat focus, respectively, and a duty cycle of 1%. Here, we present the histological analysis, including 3D reconstruction of histotripsy sites. Five whole porcine livers were retrieved fresh from the abattoir using human transplant retrieval and cold static preservation techniques and were then perfused using an organ preservation circuit. Whilst under perfusion, histotripsy was performed to randomly selected sites on the live. Fifteen lesional sites were formalin-fixed and paraffin-embedded. Sections were stained with Haematoxylin and Eosin and picro-Sirius red, and they were also stained for reticulin. Additionally, two lesion sites were used for 3D reconstruction. The core of the typical lesion consisted of eosinophilic material associated with reticulin loss, collagen damage including loss of birefringence to fibrous septa, and perilesional portal tracts, including large portal vein branches, but intact peri-lesional hepatic plates. The 3D reconstruction of two histotripsy sites was successful and confirmed the feasibility of this approach to investigate the effects of histotripsy on tissue in detail.

Liver Ultrasound Histotripsy: Novel Analysis of the Histotripsy Site Cell Constituents with Implications for Histotripsy Application in Cell Transplantation and Cancer Therapy

Introduction: Allogenic hepatocyte transplantation is an attractive alternative to whole-organ transplantation, particularly for the treatment of metabolic disorders and acute liver failure. However, the shortage of human donor organs for cell isolation, the low cell yield from decellularisation regimes, and low engraftment rates from portal administration of donor cells have restricted its clinical application. Using ultrasound histotripsy to provide a nidus in the liver for direct cell transplantation offers a new approach to overcoming key limitations in current cell therapy. We have analysed the liver cavity constituents to assess their potential as a site for cell delivery and implantation. Methods: Using human organ retrieval techniques, pig livers were collected from the abattoir and transported in ice-cold storage to the laboratory. Following 2 h of cold storage, the livers were flushed with organ preservation solution and placed on an organ perfusion circuit to maintain viability. Organs were perfused with Soltran™ organ preservation solution via the portal vein at a temperature of 24-30 °C. The perfusion circuit was oxygenated through equilibration with room air. Perfused livers (n=5) were subjected to ultrasound histotripsy, producing a total of 130 lesions. Lesions were generated by applying 50 pulses at 1 Hz pulse repetition frequency and 1% duty cycle using a single element 2 MHz bowl-shaped transducer (Sonic Concepts, H-148). Following histotripsy, a focal liver lesion was produced, which had a liquid centre. The fluid from each lesion was aspirated and cultured in medium (RPMI) at 37 °C in an incubator. Cell cultures were analysed at 1 and 7 days for cell viability and a live-dead assay was performed. The histotripsy sites were excised following aspiration and H&E staining was used to characterise the liver lesions. Cell morphology was determined by histology. Results: Histotripsy created a subcapsular lesion (~5 mm below the liver capsule; size ranging from 3 to 5 mm), which contained a suspension of cells. On average, 61×10⁴ cells per mL were isolated. Hepatocytes were present in the aspirate, were viable at 24 h post isolation and remained viable in culture for up to 1 week, as determined by phalloidin/DAPI cell viability stains. Cultures up to 21 days revealed metabolically active live hepatocyte. Live-dead assays confirmed hepatocyte viability at 1 week (Day 1: 12% to Day 7: 45% live cells; p < 0.0001), which retained metabolic activity and morphology, confirmed on assay and microscopy. Cell Titre-Glo^TM showed a peak metabolic activity at 1 week (average luminescence 24.6 RLU; p < 0.0001) post-culture compared with the control (culture medium alone), reduced to 1/3 of peak level (7.85 RLU) by day 21. Conclusions: Histotripsy of the liver allows isolation and culture of hepatocytes with a high rate of viability after 1 week in culture. Reproducing these findings using human livers may lead to wide clinical applications in cell therapy.

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.

Mechanical damage thresholds for hematomas near gas-containing bodies in pulsed HIFU fields

Objective. Boiling histotripsy (BH) is a novel high intensity focused ultrasound (HIFU) application currently being developed for non-invasive mechanical fractionation of soft tissues and large hematomas. In the context of development of BH treatment planning approaches for ablating targets adjacent to gas-containing organs, this study aimed at investigation of the ultrasound pressure thresholds of atomization-induced damage to the tissue-air interface and correlation of the danger zone dimensions with spatial structure of nonlinear HIFU field parameters.Approach. A flat interface with air of freshly clotted bovine blood was used as anex vivomodel due to its homogenous structure and higher susceptibility to ultrasound-induced mechanical damage compared to soft tissues. Three 1.5 MHz transducers of differentF-numbers (0.77, 1 and 1.5) were focused at various distances before or beyond a flat clot surface, and a BH exposure was delivered either at constant, high-amplitude output level, or at gradually increasing level until a visible damage to the clot surface occurred. The HIFU pressure field parameters at the clot surface were determined through a combination of hydrophone measurements in water, forward wave propagation simulation using 'HIFU beam' software and an image source method to account for the wave reflection from the clot surface and formation of a standing wave. The iso-levels of peak negative pressure in the resulting HIFU field were correlated to the outlines of surface erosion to identify the danger zone around the BH focus.Main results. The outline of the danger zone was shown to differ from that of a typical BH lesion produced in a volume of clot material. In the prefocal area, the zone was confined within the 4 MPa contour of the incident peak-to-peak pressure; within the main focal lobe it was determined by the maximum BH lesion width, and in the postfocal area-by the transverse size of the focal lobe and position of the first postfocal pressure axial null.Significance. The incident HIFU pressure-based danger zone boundaries were outlined around the BH focus and can be superimposed onto in-treatment ultrasound image to avoid damage to adjacent gas-containing bodies.

Investigation of the long-term healing response of the liver to boiling histotripsy treatment in vivo

Boiling histotripsy (BH) is a promising High-Intensity Focused Ultrasound technique that can be employed to mechanically fractionate solid tumours. Whilst studies have shown the feasibility of BH to destroy liver cancer, no study has reported on the healing process of BH-treated liver tissue. We therefore extensively investigated the evolution of the healing response of liver to BH in order to provide an insight into the healing mechanisms. In the present study, 14 Sprague Dawley rats underwent the BH treatment and were sacrificed on days 0, 3, 7, 14, and 28 for morphological, histological, serological and qPCR analyses. The area of the treated region was 1.44 cm² (1.2 cm × 1.2 cm). A well-defined BH lesion filled with coagulated blood formed on day 0. A week after the treatment, fibroblast activation was induced at the treatment site, leading to the formation of extracellular matrix structure (ECM). The ECM was then disrupted for 7 to 28 days. Regenerated normal hepatocytes and newly formed blood vessels were found within the BH region with the absence of hepatic fibrosis. No significant morphological, histological and genetic changes around the BH lesion occurred. These results suggest that BH could be a safe and promising therapeutic tool for treating solid tumours without inducing any significant adverse effect such as the formation of liver fibrosis.

Mechanisms of nuclei growth in ultrasound bubble nucleation

This paper interrogates the intersections between bubble dynamics and classical nucleation theory (CNT) towards constructing a model that describes intermediary nucleation events between the extrema of cavitation and boiling. We employ Zeldovich's hydrodynamic approach to obtain a description of bubble nuclei that grow simultaneously via hydrodynamic excitation by the acoustic field and vapour transport. By quantifying the relative dominance of both mechanisms, it is then possible to discern the extent to which viscosity, inertia, surface tension and vapour transport shape the growth of bubble nuclei through non-dimensional numbers that naturally arise within the theory. The first non-dimensional number Φ₁²/Φ₂ is analogous to the Laplace number, representing the balance between surface tension and inertial constraints to viscous effects. The second non-dimensional number δ represents how enthalpy transport into the bubble can reduce nucleation rates by cooling the surrounding liquid. This formulation adds to the current understanding of ultrasound bubble nucleation by accounting for bubble dynamics during nucleation, quantifying the physical distinctions between "boiling" and "cavitation" bubbles through non-dimensional parameters, and outlining the characteristic timescales of nucleation according to the growth mechanism of bubbles throughout the histotripsy temperature range. We observed in our simulations that viscous effects control the process of ultrasound nucleation in water-like media throughout the 0-120 °C temperature range, although this dominance decreases with increasing temperatures. Enthalpy transport was found to reduce nucleation rates for increasing temperatures. This effect becomes significant at temperatures above 30 °C and favours the creation of fewer nuclei that are larger in size. Conversely, negligible enthalpy transport at lower temperatures can enable the nucleation of dense clusters of small nuclei, such as cavitation clouds. We find that nuclei growth as modelled by the Rayleigh-Plesset equation occurs over shorter timescales than as modelled by vapour-dominated growth. This suggests that the first stage of bubble nuclei growth is hydrodynamic, and vapour transport effects can only be observed over longer timescales. Finally, we propose that this framework can be used for comparison between different experiments in bubble nucleation, towards standardisation and dosimetry of protocols.

An ultrasonically actuated needle promotes the transport of nanoparticles and fluids

Non-invasive therapeutic ultrasound (US) methods, such as high-intensity focused ultrasound (HIFU), have limited access to tissue targets shadowed by bones or presence of gas. This study demonstrates that an ultrasonically actuated medical needle can be used to translate nanoparticles and fluids under the action of nonlinear phenomena, potentially overcoming some limitations of HIFU. A simulation study was first conducted to study the delivery of a tracer with an ultrasonically actuated needle (33 kHz) inside a porous medium acting as a model for soft tissue. The model was then validated experimentally in different concentrations of agarose gel showing a close match with the experimental results, when diluted soot nanoparticles (diameter < 150 nm) were employed as delivered entity. An additional simulation study demonstrated a threefold increase in the volume covered by the delivered agent in liver under a constant injection rate, when compared to without US. This method, if developed to its full potential, could serve as a cost effective way to improve safety and efficacy of drug therapies by maximizing the concentration of delivered entities within, e.g., a small lesion, while minimizing exposure outside the lesion.

An ultrasonically actuated fine-needle creates cavitation in bovine liver

Ultrasonic cavitation is being used in medical applications as a way to influence matter, such as tissue or drug vehicles, on a micro-scale. Oscillating or collapsing cavitation bubbles provide transient mechanical force fields, which can, e.g., fractionate soft tissue or even disintegrate solid objects, such as calculi. Our recent study demonstrates that an ultrasonically actuated medical needle can create cavitation phenomena inside water. However, the presence and behavior of cavitation and related bioeffects in diagnostic and therapeutic applications with ultrasonically actuated needles are not known. Using simulations, we demonstrate numerically and experimentally the cavitation phenomena near ultrasonically actuated needles. We define the cavitation onset within a liver tissue model with different total acoustic power levels. We directly visualize and quantitatively characterize cavitation events generated by the ultrasonic needle in thin fresh bovine liver sections enabled by high-speed imaging. On a qualitative basis, the numerical and experimental results show a close resemblance in threshold and spatial distribution of cavitation. These findings are crucial for developing new methods and technologies employing ultrasonically actuated fine needles, such as ultrasound-enhanced fine-needle biopsy, drug delivery, and histotripsy.

Enhancing cavitation dynamics and its mechanical effects with dual-frequency ultrasound

Objective. Acoustic cavitation and its mechanical effects (e.g. stress and strain) play a primary role in ultrasound applications. Introducing encapsulated microbubbles as cavitation nuclei and utilizing dual-frequency ultrasound excitation are highly effective approaches to reduce cavitation thresholds and enhance cavitation effects. However, the cavitation dynamics of encapsulated microbubbles and the resultant stress/strain in viscoelastic tissues under dual-frequency excitation are poorly understood, especially for the enhancement effects caused by a dual-frequency approach. The goal of this study was to numerically investigate the dynamics of a lipid-coated microbubble and the spatiotemporal distributions of the stress and strain under dual-frequency excitation. Approach. The Gilmore–Zener bubble model was coupled with a shell model for the nonlinear changes of both shell elasticity and viscosity to accurately simulate the cavitation dynamics of lipid-coated microbubbles in viscoelastic tissues. Then, the spatiotemporal evolutions of the cavitation-induced stress and strain in the surrounding tissues were characterized quantitatively. Finally, the influences of some paramount parameters were examined to optimize the outcomes. Main results. We demonstrated that the cavitation dynamics and associated stress/strain were prominently enhanced by a dual-frequency excitation, highlighting positive correlations between the maximum bubble expansion and the maximum stress/strain. Moreover, the results showed that the dual-frequency ultrasound with smaller differences in its frequencies and pressure amplitudes could enhance the bubble oscillations and stress/strain more efficiently, whereas the phase difference manifested small influences under these conditions. Additionally, the dual-frequency approach seemed to show a stronger enhancement effect with the shell/tissue viscoelasticity increasing to a certain extent. Significance. This study might contribute to optimizing the dual-frequency operation in terms of cavitation dynamics and its mechanical effects for high-efficient ultrasound applications.

High-spatial-resolution, instantaneous passive cavitation imaging with temporal resolution in histotripsy: a simulation study

Purpose

In histotripsy, a shock wave is transmitted, and the resulting inertial bubble cavitation that disrupts tissue is used for treatment. Therefore, it is necessary to detect when cavitation occurs and track the position of cavitation occurrence using a new passive cavitation (PC) imaging method.

Methods

An integrated PC image, which is constructed by collecting the focused signals at all times, does not provide information on when cavitation occurs and has poor spatial resolution. To solve this problem, we constructed instantaneous PC images by applying delay and sum beamforming at instantaneous time instants. By calculating instantaneous PC images at all data acquisition times, the proposed method can detect cavitation when it occurs by using the property that when signals from the cavitation are focused, their amplitude becomes large, and it can obtain a high-resolution PC image by masking out side lobes in the vicinity of cavitation.

Results

Ultrasound image simulation confirmed that the proposed method has higher resolution than conventional integrated PC imaging and showed that it can determine the position and time of cavitation occurrence as well as the signal strength.

Conclusion

Since the proposed novel PC imaging method can detect each cavitation separately when the incidence of cavitations is low, it can be used to monitor the treatment process of shock wave therapy and histotripsy, in which cavitation is an important mechanism of treatment.

Nonlinear dynamics and acoustic emissions of interacting cavitation bubbles in viscoelastic tissues

The cavitation-mediated bioeffects are primarily associated with the dynamic behaviors of bubbles in viscoelastic tissues, which involves complex interactions of cavitation bubbles with surrounding bubbles and tissues. The radial and translational motions, as well as the resultant acoustic emissions of two interacting cavitation bubbles in viscoelastic tissues were numerically investigated. Due to the bubble-bubble interactions, a remarkable suppression effect on the small bubble, whereas a slight enhancement effect on the large one were observed within the acoustic exposure parameters and the initial radii of the bubbles examined in this paper. Moreover, as the initial distance between bubbles increases, the strong suppression effect is reduced gradually and it could effectively enhance the nonlinear dynamics of bubbles, exactly as the bifurcation diagrams exhibit a similar mode of successive period doubling to chaos. Correspondingly, the resultant acoustic emissions present a progressive evolution of harmonics, subharmonics, ultraharmonics and broadband components in the frequency spectra. In addition, with the elasticity and/or viscosity of the surrounding medium increasing, both the nonlinear dynamics and translational motions of bubbles were reduced prominently. This study provides a comprehensive insight into the nonlinear behaviors and acoustic emissions of two interacting cavitation bubbles in viscoelastic media, it may contribute to optimizing and monitoring the cavitation-mediated biomedical applications.

Cavitation Dynamics and Inertial Cavitation Threshold of Lipid Coated Microbubbles in Viscoelastic Media with Bubble-Bubble Interactions

Encapsulated microbubbles combined with ultrasound have been widely utilized in various biomedical applications; however, the bubble dynamics in viscoelastic medium have not been completely understood. It involves complex interactions of coated microbubbles with ultrasound, nearby microbubbles and surrounding medium. Here, a comprehensive model capable of simulating the complex bubble dynamics was developed via taking the nonlinear viscoelastic behaviors of the shells, the bubble–bubble interactions and the viscoelasticity of the surrounding medium into account simultaneously. For two interacting lipid-coated bubbles with different initial radii in viscoelastic media, it exemplified that the encapsulating shell, the inter-bubble interactions and the medium viscoelasticity would noticeably suppress bubble oscillations. The inter-bubble interactions exerted a much stronger suppressing effect on the small bubble within the parameters examined in this paper, which might result from a larger radiated pressure acting on the small bubble due to the inter-bubble interactions. The lipid shells make the microbubbles exhibit two typical asymmetric dynamic behaviors (i.e., compression or expansion dominated oscillations), which are determined by the initial surface tension of the bubbles. Accordingly, the inertial cavitation threshold decreases as the initial surface tension increases, but increases as the shell elasticity and viscosity increases. Moreover, with the distance between bubbles decreasing and/or the initial radius of the large bubble increasing, the oscillations of the small bubble decrease and the inertial cavitation threshold increases gradually due to the stronger suppression effects caused by the enhanced bubble–bubble interactions. Additionally, increasing the elasticity and/or viscosity of the surrounding medium would also dampen bubble oscillations and result in a significant increase in the inertial cavitation threshold. This study may contribute to both encapsulated microbubble-associated ultrasound diagnostic and emerging therapeutic applications.

Ultrastructural Analysis of Volumetric Histotripsy Bio-effects in Large Human Hematomas

Large-volume soft tissue hematomas are a serious clinical problem, which, if untreated, can have severe consequences. Current treatments are associated with significant pain and discomfort. It has been reported that in an in vitro bovine hematoma model, pulsed high-intensity focused ultrasound (HIFU) ablation, termed histotripsy, can be used to rapidly and non-invasively liquefy the hematoma through localized bubble activity, enabling fine-needle aspiration. The goals of this study were to evaluate the efficiency and speed of volumetric histotripsy liquefaction using a large in vitro human hematoma model. Large human hematoma phantoms (85 cc) were formed by recalcifying blood anticoagulated with citrate phosphate dextrose/saline-adenine-glucose-mannitol solution. Typical boiling histotripsy pulses (10 or 2 ms) or hybrid histotripsy pulses using higher-amplitude and shorter pulses (0.4 ms) were delivered at 1% duty cycle while continuously translating the HIFU focus location. Histotripsy exposures were performed under ultrasound guidance with a 1.5-MHz transducer (8-cm aperture, F# = 0.75). The volume of liquefied lesions was determined by ultrasound imaging and gross inspection. Untreated hematoma samples and samples of the liquefied lesions aspirated using a fine needle were analyzed cytologically and ultrastructurally with scanning electron microscopy. All exposures resulted in uniform liquid-filled voids with sharp edges; liquefaction speed was higher for exposures with shorter pulses and higher shock amplitudes at the focus (up to 0.32, 0.68 and 2.62 mL/min for 10-, 2- and 0.4-ms pulses, respectively). Cytological and ultrastructural observations revealed completely homogenized blood cells and fibrin fragments in the lysate. Most of the fibrin fragments were less than 20 μm in length, but a number of fragments were up to 150 μm. The lysate with residual debris of that size would potentially be amenable to fine-needle aspiration without risk for needle clogging in clinical implementation.

Modeling the Physics of Bubble Nucleation in Histotripsy

This work aims to establish a theoretical framework for the modeling of bubble nucleation in histotripsy. A phenomenological version of the classical nucleation theory was parametrized with histotripsy experimental data, fitting a temperature-dependent activity factor that harmonizes theoretical predictions and experimental data for bubble nucleation at both high and low temperatures. Simulations of histotripsy pressure and temperature fields are then used in order to understand spatial and temporal properties of bubble nucleation at varying sonication conditions. This modeling framework offers a thermodynamic understanding on the role of the ultrasound frequency, waveforms, peak focal pressures, and duty cycle on patterns of ultrasound-induced bubble nucleation. It was found that at temperatures lower than 50 °C, nucleation rates are more appreciable at very large negative pressures such as -30 MPa. For focal peak-negative pressures of -15 MPa, characteristic of boiling histotripsy, nucleation rates grow by 20 orders of magnitude in the temperature interval 60 °C-100 °C.

Control of the dynamics of a boiling vapour bubble using pressure-modulated high intensity focused ultrasound without the shock scattering effect: A first proof-of-concept study

Boiling histotripsy is a promising High-Intensity Focused Ultrasound (HIFU) technique that can be used to induce mechanical tissue fractionation at the HIFU focus via cavitation. Two different types of cavitation produced during boiling histotripsy exposure can contribute towards mechanical tissue destruction: (1) a boiling vapour bubble at the HIFU focus and (2) cavitation clouds in between the boiling bubble and the HIFU source. Control of the extent and degree of mechanical damage produced by boiling histotripsy is necessary when treating a solid tumour adjacent to normal tissue or major blood vessels. This is, however, difficult to achieve with boiling histotripsy due to the stochastic formation of the shock scattering-induced inertial cavitation clouds. In the present study, a new histotripsy method termed pressure-modulated shockwave histotripsy is proposed as an alternative to or in addition to boiling histotripsy without inducing the shock scattering effect. The proposed concept is (a) to generate a boiling vapour bubble via localised shockwave heating and (b) subsequently control its extent and lifetime through manipulating peak pressure magnitudes and a HIFU pulse length. To demonstrate the feasibility of the proposed method, bubble dynamics induced at the HIFU focus in an optically transparent liver tissue phantom were investigated using a high speed camera and a passive cavitation detection systems under a single 10, 50 or 100 ms-long 2, 3.5 or 5 MHz pressure-modulated HIFU pulse with varying peak positive and negative pressure amplitudes from 5 to 89 MPa and -3.7 to -14.6 MPa at the focus. Furthermore, a numerical simulation of 2D nonlinear wave propagation with the presence of a boiling bubble at the focus of a HIFU field was conducted by numerically solving the generalised Westervelt equation. The high speed camera experimental results showed that, with the proposed pressure-modulated shockwave histotripsy, boiling bubbles generated by shockwave heating merged together, forming a larger bubble (of the order of a few hundred micron) at the HIFU focus. This coalesced boiling bubble then persisted and maintained within the HIFU focal zone until the end of the exposure (10, 50, or 100 ms). Furthermore, and most importantly, no violent cavitation clouds which typically appear in boiling histotripsy occurred during the proposed histotripsy excitation (i.e. no shock scattering effect). This was likely because that the peak negative pressure magnitude of the backscattered acoustic field by the boiling bubble was below the cavitation cloud intrinsic threshold. The size of the coalesced boiling bubble gradually increased with the peak pressure magnitudes. In addition, with the proposed method, an oval shaped lesion with a length of 0.6 mm and a width of 0.1 mm appeared at the HIFU focus in the tissue phantom, whereas a larger lesion in the form of a tadpole (length: 2.7 mm, width: 0.3 mm) was produced by boiling histotripsy. Taken together, these results suggest that the proposed pressure-modulated shockwave histotripsy could potentially be used to induce a more spatially localised tissue destruction with a desired degree of mechanical damage through controlling the size and lifetime of a boiling bubble without the shock scattering effect.

The interaction of shockwaves with a vapour bubble in boiling histotripsy: The shock scattering effect

Boiling histotripsy is a High Intensity Focused Ultrasound (HIFU) technique which uses a number of short pulses with high acoustic pressures at the HIFU focus to induce mechanical tissue fractionation. In boiling histotripsy, two different types of acoustic cavitation contribute towards mechanical tissue destruction: a boiling vapour bubble and cavitation clouds. An understanding of the mechanisms underpinning these phenomena and their dynamics is therefore paramount to predicting and controlling the overall size of a lesion produced for a given boiling histotripsy exposure condition. A number of studies have shown the effects of shockwave heating in generating a boiling bubble at the HIFU focus and have studied its dynamics under boiling histotripsy insonation. However, not much is known about the subsequent production of cavitation clouds that form between the HIFU transducer and the boiling bubble. The main objective of the present study is to examine what causes this bubble cluster formation after the generation of a boiling vapour bubble. A numerical simulation of 2D nonlinear wave propagation with the presence of a bubble at the focus of a HIFU field was performed using the k-Wave MATLAB toolbox for time domain ultrasound simulations, which numerically solves the generalised Westervelt equation. The numerical results clearly demonstrate the appearance of the constructive interference of a backscattered shockwave by a bubble with incoming incident shockwaves. This interaction (i.e., the reflected and inverted peak positive phase from the bubble with the incoming incident rarefactional phase) can eventually induce a greater peak negative pressure field compared to that without the bubble at the HIFU focus. In addition, the backscattered peak negative pressure magnitude gradually increased from 17.4 MPa to 31.6 MPa when increasing the bubble size from 0.2 mm to 1.5 mm. The latter value is above the intrinsic cavitation threshold of -28 MPa in soft tissue. Our results suggest that the formation of a cavitation cloud in boiling histotripsy is a threshold effect which primarily depends (a) the size and location of a boiling bubble, and (b) the sum of the incident field and that scattered by a bubble.

Histotripsy: the first noninvasive, non-ionizing, non-thermal ablation technique based on ultrasound

Histotripsy is the first noninvasive, non-ionizing, and non-thermal ablation technology guided by real-time imaging. Using focused ultrasound delivered from outside the body, histotripsy mechanically destroys tissue through cavitation, rendering the target into acellular debris. The material in the histotripsy ablation zone is absorbed by the body within 1-2 months, leaving a minimal remnant scar. Histotripsy has also been shown to stimulate an immune response and induce abscopal effects in animal models, which may have positive implications for future cancer treatment. Histotripsy has been investigated for a wide range of applications in preclinical studies, including the treatment of cancer, neurological diseases, and cardiovascular diseases. Three human clinical trials have been undertaken using histotripsy for the treatment of benign prostatic hyperplasia, liver cancer, and calcified valve stenosis. This review provides a comprehensive overview of histotripsy covering the origin, mechanism, bioeffects, parameters, instruments, and the latest results on preclinical and human studies.

Effect of Stiffness of Large Extravascular Hematomas on Their Susceptibility to Boiling Histotripsy Liquefaction in Vitro

Large intra-abdominal, retroperitoneal and intramuscular hematomas are common consequences of sharp and blunt trauma and post-surgical bleeds, and often threaten organ failure, compartment syndrome or spontaneous infection. Current therapy options include surgical evacuation and placement of indwelling drains that are not effective because of the viscosity of the organized hematoma. We have previously reported the feasibility of using boiling histotripsy (BH)-a pulsed high-intensity focused ultrasound method-for liquefaction of large volumes of freshly coagulated blood and subsequent fine-needle aspiration. The goal of this work was to evaluate the changes in stiffness of large coagulated blood volumes with aging and retraction in vitro, and to correlate these changes with the size of the BH void and, therefore, the susceptibility of the material to BH liquefaction. Large-volume (55-200 mL) whole-blood clots were fabricated in plastic molds from human and bovine blood, either by natural clotting or by recalcification of anticoagulated blood, with or without addition of thrombin. Retraction of the clots was achieved by incubation for 3 h, 3 d or 8 d. The shear modulus of the samples was measured with a custom-built indentometer and shear wave elasticity (SWE) imaging. Sizes of single liquefied lesions produced with a 1.5-MHz high-intensity focused ultrasound transducer within a 30-s standard BH exposure served as the metric for susceptibility of clot material to this treatment. Neither the shear moduli of naturally clotted human samples (0.52 ± 0.08 kPa), nor their degree of retraction (ratio of expelled fluid to original volume 50%-58%) depended on the length of incubation within 0-8 d, and were significantly lower than those of bovine samples (2.85 ± 0.17 kPa, retraction 5%-38%). In clots made from anticoagulated bovine blood, the variation of calcium chloride concentration within 5-40 mmol/L did not change the stiffness, whereas lower concentrations and the addition of thrombin resulted in significantly softer clots, similar to naturally clotted human samples. Within the achievable shear modulus range (0.4-1.6 kPa), the width of the BH-liquefied lesion was more affected by the changes in stiffness than the length of the lesion. In all cases, however, the lesions were larger compared with any soft tissue liquefied with the same BH parameters, indicating higher susceptibility of hematomas to BH damage. These results suggest that clotted bovine blood with added thrombin is an acceptable in vitro model of both acute and chronic human hematomas for assessing the efficiency of BH liquefaction strategies.

Pilot in vivo studies on transcutaneous boiling histotripsy in porcine liver and kidney

Boiling histotripsy (BH) is a High Intensity Focused Ultrasound (HIFU) method for precise mechanical disintegration of target tissue using millisecond-long pulses containing shocks. BH treatments with real-time ultrasound (US) guidance allowed by BH-generated bubbles were previously demonstrated ex vivo and in vivo in exposed porcine liver and small animals. Here, the feasibility of US-guided transabdominal and partially transcostal BH ablation of kidney and liver in an acute in vivo swine model was evaluated for 6 animals. BH parameters were: 1.5 MHz frequency, 5–30 pulses of 1–10 ms duration per focus, 1% duty cycle, peak acoustic powers 0.9–3.8 kW, sonication foci spaced 1–1.5 mm apart in a rectangular grid with 5–15 mm linear dimensions. In kidneys, well-demarcated volumetric BH lesions were generated without respiratory gating and renal medulla and collecting system were more resistant to BH than cortex. The treatment was accelerated 10-fold by using shorter BH pulses of larger peak power without affecting the quality of tissue fractionation. In liver, respiratory motion and aberrations from subcutaneous fat affected the treatment but increasing the peak power provided successful lesion generation. These data indicate BH is a promising technology for transabdominal and transcostal mechanical ablation of tumors in kidney and liver.

Modeling tissue-selective cavitation damage

The destructive growth and collapse of cavitation bubbles are used for therapeutic purposes in focused ultrasound procedures and can contribute to tissue damage in traumatic injuries. Histotripsy is a focused ultrasound procedure that relies on controlled cavitation to homogenize soft tissue. Experimental studies of histotripsy cavitation have shown that the extent of ablation in different tissues depends on tissue mechanical properties and waveform parameters. Variable tissue susceptibility to the large stresses, strains, and strain rates developed by cavitation bubbles has been suggested as a basis for localized liver tumor treatments that spare large vessels and bile ducts. However, field quantities developed within microns of cavitation bubbles are too localized and transient to measure in experiments. Previous numerical studies have attempted to circumvent this challenge but made limited use of realistic tissue property data. In this study, numerical simulations are used to calculate stress, strain, and strain rate fields produced by bubble oscillation under histotripsy forcing in a variety of tissues with literature-sourced viscoelastic and acoustic properties. Strain field calculations are then used to predict a theoretical damage radius using tissue ultimate strain data. Simulation results support the hypothesis that differential tissue responses could be used to design tissue-selective treatments. Results agree with studies correlating tissue ultimate fractional strain with resistance to histotripsy ablation and are also consistent with experiments demonstrating smaller lesion size under exposure to higher frequency waveforms. Methods presented in this study provide an approach for modeling tissue-selective cavitation damage in general.

Boiling Histotripsy-induced Partial Mechanical Ablation Modulates Tumour Microenvironment by Promoting Immunogenic Cell Death of Cancers

Boiling histotripsy is a promising non-invasive High-Intensity Focused Ultrasound (HIFU) technique that employs HIFU mechanical effects to fractionate solid tumours without causing any significant thermal damage. It has been suggested that boiling histotripsy may induce a strong immune response due to the absence of denatured antigenic protein at the HIFU focus. The underlying immunological mechanisms of this technique are, however, poorly understood. In this study, we demonstrated the feasibility of using boiling histotripsy to mechanically fractionate human breast adenocarcinoma cells (MDA-MB-231) and the potential immunological effects induced by boiling histotripsy, for the first time. Our results showed that mechanical stresses produced by boiling histotripsy promote immunogenic cell death of cancer cells via TNF-induced necrosis signaling pathway. This immunogenic cell death significantly increases secretions of damage-associated molecular patterns (CRT, HSP70, HMGB-1), pro-inflammatory cytokines (IFN-γ, IL-1α, IL-1β, IL-18) and chemokines (IL-8) which are related to M1 macrophage activation. Furthermore, the levels of these signaling proteins increase with the degree of mechanical damage induced by boiling histotripsy. Together, the results presented can suggest that boiling histotripsy could be a potential therapeutic approach for not only mechanically destroying solid tumours (e.g., breast cancer) but also promoting immunogenic cell death via TNF-induced necrosis to trigger antitumour immunity.