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.

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.

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.

Noninvasive acoustic manipulation of objects in a living body.

In certain medical applications, transmitting an ultrasound beam through the skin to manipulate a solid object within the human body would be beneficial. Such applications include, for example, controlling an ingestible camera or expelling a kidney stone. In this paper, ultrasound beams of specific shapes were designed by numerical modeling and produced using a phased array. These beams were shown to levitate and electronically steer solid objects (3-mm-diameter glass spheres), along preprogrammed paths, in a water bath, and in the urinary bladders of live pigs. Deviation from the intended path was on average <10%. No injury was found on the bladder wall or intervening tissue.

Robust and durable aberrative and absorptive phantom for therapeutic ultrasound applications.

Phase aberration induced by soft tissue inhomogeneities often complicates high-intensity focused ultrasound (HIFU) therapies by distorting the field and, previously, we designed and fabricated a bilayer gel phantom to reproducibly mimic that effect. A surface pattern containing size scales relevant to inhomogeneities of a porcine body wall was introduced between gel materials with fat- and muscle-like acoustic properties-ballistic and polyvinyl alcohol gels. Here, the phantom design was refined to achieve relevant values of ultrasound absorption and scattering and make it more robust, facilitating frequent handling and use in various experimental arrangements. The fidelity of the interfacial surface of the fabricated phantom to the design was confirmed by three-dimensional ultrasound imaging. The HIFU field distortions-displacement of the focus, enlargement of the focal region, and reduction of focal pressure-produced by the phantom were characterized using hydrophone measurements with a 1.5 MHz 256-element HIFU array and found to be similar to those induced by an ex vivo porcine body wall. A phase correction approach was used to mitigate the aberration effect on nonlinear focal waveforms and enable boiling histotripsy treatments through the phantom or body wall. The refined phantom represents a practical tool to explore HIFU therapy systems capabilities.

Holographic extraction of plane waves from an ultrasound beam for acoustic characterization of an absorbing layer of finite dimensions.

For the acoustic characterization of materials, a method is proposed for interpreting experiments with finite-sized transducers and test samples in terms of the idealized situation in which plane waves are transmitted through an infinite plane-parallel layer. The method uses acoustic holography, which experimentally provides complete knowledge of the wave field by recording pressure waveforms at points on a surface intersected by the acoustic beam. The measured hologram makes it possible to calculate the angular spectrum of the beam to decompose the field into a superposition of plane waves propagating in different directions. Because these waves cancel one another outside the beam, the idealized geometry of an infinite layer can be represented by a sample of finite size if its lateral dimensions exceed the width of the acoustic beam. The proposed method relies on holograms that represent the acoustic beam with and without the test sample in the transmission path. The method is described theoretically, and its capabilities are demonstrated experimentally for silicone rubber samples by measuring their frequency-dependent phase velocities and absorption coefficients in the megahertz frequency range.

Bilayer aberration-inducing gel phantom for high intensity focused ultrasound applications.

Aberrations induced by soft tissue inhomogeneities often complicate high-intensity focused ultrasound (HIFU) therapies. In this work, a bilayer phantom made from polyvinyl alcohol hydrogel and ballistic gel was built to mimic alternating layers of water-based and lipid tissues characteristic of an abdominal body wall and to reproducibly distort HIFU fields. The density, sound speed, and attenuation coefficient of each material were measured using a homogeneous gel layer. A surface with random topographical features was designed as an interface between gel layers using a 2D Fourier spectrum approach and replicating different spatial scales of tissue inhomogeneities. Distortion of the field of a 256-element 1.5 MHz HIFU array by the phantom was characterized through hydrophone measurements for linear and nonlinear beam focusing and compared to the corresponding distortion induced by an ex vivo porcine body wall of the same thickness. Both spatial shift and widening of the focal lobe were observed, as well as dramatic reduction in focal pressures caused by aberrations. The results suggest that the phantom produced levels of aberration that are similar to a real body wall and can serve as a research tool for studying HIFU effects as well as for developing algorithms for aberration correction.

Simulation of nonlinear trans-skull focusing and formation of shocks in brain using a fully populated ultrasound array with aberration correction.

Multi-element high-intensity focused ultrasound phased arrays in the shape of hemispheres are currently used in clinics for thermal lesioning in deep brain structures. Certain side effects of overheating non-targeted tissues and skull bones have been revealed. Here, an approach is developed to mitigate these effects. A specific design of a fully populated 256-element 1-MHz array shaped as a spherical segment (F-number, F# = 1) and filled by randomly distributed equal-area polygonal elements is proposed. Capability of the array to generate high-amplitude shock fronts at the focus is tested in simulations by combining three numerical algorithms for linear and nonlinear field modeling and aberration correction. The algorithms are based on the combination of the Rayleigh integral, a linear pseudo-spectral time domain Kelvin-Voigt model, and nonlinear Westervelt model to account for the effects of inhomogeneities, aberrations, reflections, absorption, nonlinearity, and shear waves in the skull. It is shown that the proposed array can generate nonlinear waveforms with shock amplitudes >60 MPa at the focus deep inside the brain without exceeding the existing technical limitation on the intensity of 40 W/cm2 at the array elements. Such shock amplitudes are sufficient for mechanical ablation of brain tissues using the boiling histotripsy approach and implementation of other shock-based therapies.

Irregular reflection of spark-generated shock pulses from a rigid surface: Mach-Zehnder interferometry measurements in air.

The irregular reflection of weak acoustic shock waves, known as the von Neumann reflection, has been observed experimentally and numerically for spherically diverging waves generated by an electric spark source. Two optical measurement methods are used: a Mach-Zehnder interferometer for measuring pressure waveforms and a Schlieren system for visualizing shock fronts. Pressure waveforms are reconstructed from the light phase difference measured by the interferometer using the inverse Abel transform. In numerical simulations, the axisymmetric Euler equations are solved using finite-difference time-domain methods and the spark source is modeled as an instantaneous energy injection with a Gaussian shape. Waveforms and reflection patterns obtained from the simulations are in good agreement with those measured by the interferometer and the Schlieren methods. The Mach stem formation is observed close to the surface for incident pressures within the range of 800 to 4000 Pa. Similarly, as for strong shocks generated by blasts, it is found that for spherical weak shocks the Mach stem length increases with distance following a parabolic law. This study confirms the occurrence of irregular reflections at acoustic pressure levels and demonstrates the benefits of the Mach-Zehnder interferometer method when microphone measurements cannot be applied.

Effect of surface roughness on nonlinear reflection of weak shock waves.

The authors have recently shown that irregular reflections of spark-generated pressure weak shocks from a smooth rigid surface can be studied using an optical interferometer [Karzova, Lechat, Ollivier, Dragna, Yuldashev, Khokhlova, and Blanc-Benon, J. Acoust. Soc. Am. 145(1), 26-35 (2019)]. The current study extends these results to the reflection from rough surfaces. A Mach-Zehnder interferometer is used to measure pressure waveforms. Simulations are based on the solution of axisymmetric Euler equations. It is shown that roughness causes a decrease of the Mach stem height and the appearance of oscillations in the pressure waveforms. Close to rough surfaces, the pressure was higher compared to the smooth surface.

Ultrasound-guided tissue fractionation by high intensity focused ultrasound in an in vivo porcine liver model.

The clinical use of high intensity focused ultrasound (HIFU) therapy for noninvasive tissue ablation has been recently gaining momentum. In HIFU, ultrasound energy from an extracorporeal source is focused within the body to ablate tissue at the focus while leaving the surrounding organs and tissues unaffected. Most HIFU therapies are designed to use heating effects resulting from the absorption of ultrasound by tissue to create a thermally coagulated treatment volume. Although this approach is often successful, it has its limitations, such as the heat sink effect caused by the presence of a large blood vessel near the treatment area or heating of the ribs in the transcostal applications. HIFU-induced bubbles provide an alternative means to destroy the target tissue by mechanical disruption or, at its extreme, local fractionation of tissue within the focal region. Here, we demonstrate the feasibility of a recently developed approach to HIFU-induced ultrasound-guided tissue fractionation in an in vivo pig model. In this approach, termed boiling histotripsy, a millimeter-sized boiling bubble is generated by ultrasound and further interacts with the ultrasound field to fractionate porcine liver tissue into subcellular debris without inducing further thermal effects. Tissue selectivity, demonstrated by boiling histotripsy, allows for the treatment of tissue immediately adjacent to major blood vessels and other connective tissue structures. Furthermore, boiling histotripsy would benefit the clinical applications, in which it is important to accelerate resorption or passage of the ablated tissue volume, diminish pressure on the surrounding organs that causes discomfort, or insert openings between tissues.

Statistics of peak overpressure and shock steepness for linear and nonlinear N-wave propagation in a kinematic turbulence.

Linear and nonlinear propagation of high amplitude acoustic pulses through a turbulent layer in air is investigated using a two-dimensional KZK-type (Khokhlov-Zabolotskaya-Kuznetsov) equation. Initial waves are symmetrical N-waves with shock fronts of finite width. A modified von Kármán spectrum model is used to generate random wind velocity fluctuations associated with the turbulence. Physical parameters in simulations correspond to previous laboratory scale experiments where N-waves with 1.4 cm wavelength propagated through a turbulence layer with the outer scale of about 16 cm. Mean value and standard deviation of peak overpressure and shock steepness, as well as cumulative probabilities to observe amplified peak overpressure and shock steepness, are analyzed. Nonlinear propagation effects are shown to enhance pressure level in random foci for moderate initial amplitudes of N-waves thus increasing the probability to observe highly peaked waveforms. Saturation of the pressure level is observed for stronger nonlinear effects. It is shown that in the linear propagation regime, the turbulence mainly leads to the smearing of shock fronts, thus decreasing the probability to observe high values of steepness, whereas nonlinear effects dramatically increase the probability to observe steep shocks.

Shock formation and nonlinear saturation effects in the ultrasound field of a diagnostic curvilinear probe.

Newer imaging and therapeutic ultrasound technologies may benefit from in situ pressure levels higher than conventional diagnostic ultrasound. One example is the recently developed use of ultrasonic radiation force to move kidney stones and residual fragments out of the urinary collecting system. A commercial diagnostic 2.3 MHz C5-2 array probe has been used to deliver the acoustic pushing pulses. The probe is a curvilinear array comprising 128 elements equally spaced along a convex cylindrical surface. The effectiveness of the treatment can be increased by using higher transducer output to provide a stronger pushing force; however nonlinear acoustic saturation can be a limiting factor. In this work nonlinear propagation effects were analyzed for the C5-2 transducer using a combined measurement and modeling approach. Simulations were based on the three-dimensional Westervelt equation with the boundary condition set to match low power measurements of the acoustic pressure field. Nonlinear focal waveforms simulated for different numbers of operating elements of the array at several output power levels were compared to fiber-optic hydrophone measurements and were found to be in good agreement. It was shown that saturation effects do limit the acoustic pressure in the focal region of a diagnostic imaging probe.

Histotripsy methods in mechanical disintegration of tissue: towards clinical applications.

In high intensity focused ultrasound (HIFU) therapy, an ultrasound beam is focused within the body to locally affect the targeted site without damaging intervening tissues. The most common HIFU regime is thermal ablation. Recently there has been increasing interest in generating purely mechanical lesions in tissue (histotripsy). This paper provides an overview of several studies on the development of histotripsy methods toward clinical applications. Two histotripsy approaches and examples of their applications are presented. In one approach, sequences of high-amplitude, short (microsecond-long), focused ultrasound pulses periodically produce dense, energetic bubble clouds that mechanically disintegrate tissue. In an alternative approach, longer (millisecond-long) pulses with shock fronts generate boiling bubbles and the interaction of shock fronts with the resulting vapour cavity causes tissue disintegration. Recent preclinical studies on histotripsy are reviewed for treating benign prostatic hyperplasia (BPH), liver and kidney tumours, kidney stone fragmentation, enhancing anti-tumour immune response, and tissue decellularisation for regenerative medicine applications. Potential clinical advantages of the histotripsy methods are discussed. Histotripsy methods can be used to mechanically ablate a wide variety of tissues, whilst selectivity sparing structures such as large vessels. Both ultrasound and MR imaging can be used for targeting and monitoring the treatment in real time. Although the two approaches utilise different mechanisms for tissue disintegration, both have many of the same advantages and offer a promising alternative method of non-invasive surgery.

Acoustic holography as a metrological tool for characterizing medical ultrasound sources and fields.

Acoustic holography is a powerful technique for characterizing ultrasound sources and the fields they radiate, with the ability to quantify source vibrations and reduce the number of required measurements. These capabilities are increasingly appealing for meeting measurement standards in medical ultrasound; however, associated uncertainties have not been investigated systematically. Here errors associated with holographic representations of a linear, continuous-wave ultrasound field are studied. To facilitate the analysis, error metrics are defined explicitly, and a detailed description of a holography formulation based on the Rayleigh integral is provided. Errors are evaluated both for simulations of a typical therapeutic ultrasound source and for physical experiments with three different ultrasound sources. Simulated experiments explore sampling errors introduced by the use of a finite number of measurements, geometric uncertainties in the actual positions of acquired measurements, and uncertainties in the properties of the propagation medium. Results demonstrate the theoretical feasibility of keeping errors less than about 1%. Typical errors in physical experiments were somewhat larger, on the order of a few percent; comparison with simulations provides specific guidelines for improving the experimental implementation to reduce these errors. Overall, results suggest that holography can be implemented successfully as a metrological tool with small, quantifiable errors.

Mach stem formation in reflection and focusing of weak shock acoustic pulses.

The aim of this study is to show the evidence of Mach stem formation for very weak shock waves with acoustic Mach numbers on the order of 10(-3) to 10(-2). Two representative cases are considered: reflection of shock pulses from a rigid surface and focusing of nonlinear acoustic beams. Reflection experiments are performed in air using spark-generated shock pulses. Shock fronts are visualized using a schlieren system. Both regular and irregular types of reflection are observed. Numerical simulations are performed to demonstrate the Mach stem formation in the focal region of periodic and pulsed nonlinear beams in water.

Characterization of spark-generated N-waves in air using an optical schlieren method.

Accurate measurement of high-amplitude, broadband shock pulses in air is an important part of laboratory-scale experiments in atmospheric acoustics. Although various methods have been developed, specific drawbacks still exist and need to be addressed. Here, a schlieren optical method was used to reconstruct the pressure signatures of nonlinear spherically diverging short acoustic pulses generated using an electric spark source (2.5 kPa, 33 µs at 10 cm from the source) in homogeneous air. A high-speed camera was used to capture light rays deflected by refractive index inhomogeneities, caused by the acoustic wave. Pressure waveforms were reconstructed from the light intensity patterns in the recorded images using an Abel-type inversion method. Absolute pressure levels were determined by analyzing at different propagation distances the duration of the compression phase of pulses, which changed due to nonlinear propagation effects. Numerical modeling base on the generalized Burgers equation was used to evaluate the smearing of the waveform caused by finite exposure time of the high-speed camera and corresponding limitations in resolution of the schlieren technique. The proposed method allows the study of the evolution of spark-generated shock waves in air starting from the very short distances from the spark, 30 mm, up to 600 mm.

Dynamic Mode Decomposition for Transient Cavitation Bubbles Imaging in Pulsed High-Intensity Focused Ultrasound Therapy.

Pulsed high-intensity focused ultrasound (pHIFU) can induce sparse de novo inertial cavitation without the introduction of exogenous contrast agents, promoting mild mechanical disruption in targeted tissue. Because the bubbles are small and rapidly dissolve after each HIFU pulse, mapping transient bubbles and obtaining real-time quantitative metrics correlated with tissue damage are challenging. Prior work introduced Bubble Doppler, an ultrafast power Doppler imaging method as a sensitive means to map cavitation bubbles. The main limitation of that method was its reliance on conventional wall filters used in Doppler imaging and its optimization for imaging blood flow rather than transient scatterers. This study explores Bubble Doppler enhancement using dynamic mode decomposition (DMD) of a matrix created from a Doppler ensemble for mapping and extracting the characteristics of transient cavitation bubbles. DMD was first tested in silico with a numerical dataset mimicking the spatiotemporal characteristics of backscattered signal from tissue and bubbles. The performance of DMD filter was compared to other widely used Doppler wall filter-singular value decomposition (SVD) and infinite impulse response (IIR) high-pass filter. DMD was then applied to an ex vivo tissue dataset where each HIFU pulse was immediately followed by a plane wave Doppler ensemble. In silico DMD outperformed SVD and IIR high-pass filter and ex vivo provided physically interpretable images of the modes associated with bubbles and their corresponding temporal decay rates. These DMD modes can be trackable over the duration of pHIFU treatment using k-means clustering method, resulting in quantitative indicators of treatment progression.

Dynamic Mode Decomposition for Transient Cavitation Bubbles Imaging in Pulsed High Intensity Focused Ultrasound Therapy.

Pulsed high-intensity focused ultrasound (pHIFU) can induce sparse de novo inertial cavitation without the introduction of exogenous contrast agents, promoting mild mechanical disruption in targeted tissue. Because the bubbles are small and rapidly dissolve after each HIFU pulse, mapping transient bubbles and obtaining real-time quantitative metrics correlated to tissue damage are challenging. Prior work introduced Bubble Doppler, an ultrafast power Doppler imaging method as a sensitive means to map cavitation bubbles. The main limitation of that method was its reliance on conventional wall filters used in Doppler imaging and optimized for imaging blood flow rather than transient scatterers. This study explores Bubble Doppler enhancement using dynamic mode decomposition (DMD) of a matrix created from a Doppler ensemble for mapping and extracting the characteristics of transient cavitation bubbles. DMD was first tested in silico with a numerical dataset mimicking the spatiotemporal characteristics of backscattered signal from tissue and bubbles. The performance of DMD filter was compared to other widely used Doppler wall filters - singular value decomposition (SVD) and infinite impulse response (IIR) highpass filter. DMD was then applied to an ex vivo tissue dataset where each HIFU pulse was immediately followed by a plane wave Doppler ensemble. In silico DMD outperformed SVD and IIR high pass filter and ex vivo provided physically interpretable images of the modes associated with bubbles and their corresponding temporal decay rates. These DMD modes can be trackable over the duration of pHIFU treatment using k-means clustering method, resulting in quantitative indicators of treatment progression.

Treatment Planning and Aberration Correction Algorithm for HIFU Ablation of Renal Tumors.

High-intensity focused ultrasound (HIFU) applications for thermal or mechanical ablation of renal tumors often encounter challenges due to significant beam aberration and refraction caused by oblique beam incidence, inhomogeneous tissue layers, and presence of gas and bones within the beam. These losses can be significantly mitigated through sonication geometry planning, patient positioning, and aberration correction using multielement phased arrays. Here, a sonication planning algorithm is introduced, which uses the simulations to select the optimal transducer position and evaluate the effect of aberrations and acoustic field quality at the target region after aberration correction. Optimization of transducer positioning is implemented using a graphical user interface (GUI) to visualize a segmented 3-D computed tomography (CT)-based acoustic model of the body and to select sonication geometry through a combination of manual and automated approaches. An HIFU array (1.5 MHz, 256 elements) and three renal cell carcinoma (RCC) cases with different tumor locations and patient body habitus were considered. After array positioning, the correction of aberrations was performed using a combination of backpropagation from the focus with an ordinary least squares (OLS) optimization of phases at the array elements. The forward propagation was simulated using a combination of the Rayleigh integral and k-space pseudospectral method (k-Wave toolbox). After correction, simulated HIFU fields showed tight focusing and up to threefold higher maximum pressure within the target region. The addition of OLS optimization to the aberration correction method yielded up to 30% higher maximum pressure compared to the conventional backpropagation and up to 250% higher maximum pressure compared to the ray-tracing method, particularly in strongly distorted cases.