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.

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.

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.

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. Expected final online publication date for the Annual Review of Biomedical Engineering, Volume 26 is May 2024. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

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.

Development of an Automated Ultrasound Signal Indicator of Lung Interstitial Syndrome.

OBJECTIVES: The number and distribution of lung ultrasound (LUS) imaging artifacts termed B-lines correlate with the presence of acute lung disease such as infection, acute respiratory distress syndrome (ARDS), and pulmonary edema. Detection and interpretation of B-lines require dedicated training and is machine and operator-dependent. The goal of this study was to identify radio frequency (RF) signal features associated with B-lines in a cohort of patients with cardiogenic pulmonary edema. A quantitative signal indicator could then be used in a single-element, non-imaging, wearable, automated lung ultrasound sensor (LUSS) for continuous hands-free monitoring of lung fluid. METHODS: In this prospective study a 10-zone LUS exam was performed in 16 participants, including 12 patients admitted with acute cardiogenic pulmonary edema (mean age 60 ± 12 years) and 4 healthy controls (mean age 44 ± 21). Overall,160 individual LUS video clips were recorded. The LUS exams were performed with a phased array probe driven by an open-platform ultrasound system with simultaneous RF signal collection. RF data were analyzed offline for candidate B-line indicators based on signal amplitude, temporal variability, and frequency spectrum; blinded independent review of LUS images for the presence or absence of B-lines served as ground truth. Predictive performance of the signal indicators was determined with receiving operator characteristic (ROC) analysis with k-fold cross-validation. RESULTS: Two RF signal features-temporal variability of signal amplitude at large depths and at the pleural line-were strongly associated with B-line presence. The sensitivity and specificity of a combinatorial indicator were 93.2 and 58.5%, respectively, with cross-validated area under the ROC curve (AUC) of 0.91 (95% CI = 0.80-0.94). CONCLUSION: A combinatorial signal indicator for use with single-element non-imaging LUSS was developed to facilitate continuous monitoring of lung fluid in patients with respiratory illness.

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.

Enhancement of Boiling Histotripsy by Steering the Focus Axially During the Pulse Delivery.

Boiling histotripsy (BH) is a pulsed high-intensity focused ultrasound (HIFU) method relying on the generation of high-amplitude shocks at the focus, localized enhanced shock-wave heating, and bubble activity driven by shocks to induce tissue liquefaction. BH uses sequences of 1-20 ms long pulses with shock fronts of over 60 MPa amplitude, initiates boiling at the focus of the HIFU transducer within each pulse, and the remainder shocks of the pulse then interact with the boiling vapor cavities. One effect of this interaction is the creation of a prefocal bubble cloud due to reflection of shocks from the initially generated mm-sized cavities: the shocks are inverted when reflected from a pressure-release cavity wall resulting in sufficient negative pressure to reach intrinsic cavitation threshold in front of the cavity. Secondary clouds then form due to shock-wave scattering from the first one. Formation of such prefocal bubble clouds has been known as one of the mechanisms of tissue liquefaction in BH. Here, a methodology is proposed to enlarge the axial dimension of this bubble cloud by steering the HIFU focus toward the transducer after the initiation of boiling until the end of each BH pulse and thus to accelerate treatment. A BH system comprising a 1.5 MHz 256-element phased array connected to a Verasonics V1 system was used. High-speed photography of BH sonications in transparent gels was performed to observe the extension of the bubble cloud resulting from shock reflections and scattering. Volumetric BH lesions were then generated in ex vivo tissue using the proposed approach. Results showed up to almost threefold increase of the tissue ablation rate with axial focus steering during the BH pulse delivery compared to standard BH.

Dual-Mode 1-D Linear Ultrasound Array for Image-Guided Drug Delivery Enhancement Without Ultrasound Contrast Agents.

Pulsed high-intensity focused ultrasound (pHIFU) uses nonlinearly distorted millisecond-long ultrasound pulses of moderate intensity to induce inertial cavitation in tissue without administration of contrast agents. The resulting mechanical disruption permeabilizes the tissue and enhances the diffusion of systemically administered drugs. This is especially beneficial for tissues with poor perfusion such as pancreatic tumors. Here, we characterize the performance of a dual-mode ultrasound array designed for image-guided pHIFU therapies in producing inertial cavitation and ultrasound imaging. The 64-element linear array (1.071 MHz, an aperture of 14.8×51.2 mm, and a pitch of 0.8 mm) with an elevational focal length of 50 mm was driven by the Verasonics V-1 ultrasound system with extended burst option. The attainable focal pressures and electronic steering range in linear and nonlinear operating regimes (relevant to pHIFU treatments) were characterized through hydrophone measurements, acoustic holography, and numerical simulations. The steering range at ±10% from the nominal focal pressure was found to be ±6 mm axially and ±11 mm azimuthally. Focal waveforms with shock fronts of up to 45 MPa and peak negative pressures up to 9 MPa were achieved at focusing distances of 38-75 mm from the array. Cavitation behaviors induced by isolated 1-ms pHIFU pulses in optically transparent agarose gel phantoms were observed by high-speed photography across a range of excitation amplitudes and focal distances. For all focusing configurations, the appearance of sparse, stationary cavitation bubbles occurred at the same P- threshold of 2 MPa. As the output level increased, a qualitative change in cavitation behavior occurred, to pairs and sets of proliferating bubbles. The pressure P- at which this transition was observed corresponded to substantial nonlinear distortion and shock formation in the focal region and was thus dependent on the focal distance of the beam ranging within 3-4 MPa for azimuthal F -numbers of 0.74-1.5. The array was capable of B-mode imaging at 1.5 MHz of centimeter-sized targets in phantoms and in vivo pig tissues at depths of 3-7 cm, relevant to pHIFU applications in abdominal targets.

Comparative Characterization of Nonlinear Ultrasound Fields Generated by Sonalleve V1 and V2 MR-HIFU Systems.

A Sonalleve magnetic resonance-guided high-intensity focused ultrasound (MR-HIFU) clinical system (Profound Medical, Mississauga, ON, Canada) has been shown to generate nonlinear ultrasound fields with shocks up to 100 MPa at the focus as required for HIFU applications such as boiling histotripsy of hepatic and renal tumors. The Sonalleve system has two versions V1 and V2 of the therapeutic array, with differences in focusing angle, focus depth, arrangement of elements, and the size of a central opening that is twice larger in the V2 system compared to the V1. The goal of this study was to compare the performance of the V1 and V2 transducers for generating high-amplitude shock-wave fields and to reveal the impact of different array geometries on shock amplitudes at the focus. Nonlinear modeling of the field in water using boundary conditions reconstructed from holography measurements shows that at the same power output, the V2 array generates 10-15-MPa lower shock amplitudes at the focus. Consequently, substantially higher power levels are required for the V2 system to reach the same shock-wave exposure conditions in histotripsy-type treatments. Although this difference is mainly caused by the smaller focusing angle of the V2 array, the larger central opening of the V2 array has a nontrivial impact. By excluding coherently interacting weakly focused waves coming from the central part of the source, the presence of the central opening results in a somewhat higher effective focusing angle and thus higher shock amplitudes at the focus. Axisymmetric equivalent source models were constructed for both arrays, and the importance of including the central opening was demonstrated. These models can be used in the "HIFU beam" software for simulating nonlinear fields of the Sonalleve V1 and V2 systems in water and flat-layered biological tissues.

Chronic effects of pulsed high intensity focused ultrasound aided delivery of gemcitabine in a mouse model of pancreatic cancer.

Pulsed high intensity focused ultrasound (pHIFU) is a non-invasive method that allows to permeabilize pancreatic tumors through inertial cavitation and thereby increase the concentration of systemically administered drug. In this study the tolerability of weekly pHIFU-aided administrations of gemcitabine (gem) and their influence on tumor progression and immune microenvironment were investigated in genetically engineered KrasLSL.G12D/þ; p53R172H/þ; PdxCretg/þ (KPC) mouse model of spontaneously occurring pancreatic tumors. KPC mice were enrolled in the study when the tumor size reached 4-6 mm and treated once a week with either ultrasound-guided pHIFU (1.5 MHz transducer, 1 ms pulses, 1% duty cycle, peak negative pressure 16.5 MPa) followed by administration of gem (n = 9), gem only (n = 5) or no treatment (n = 8). Tumor progression was followed by ultrasound imaging until the study endpoint (tumor size reaching 1 cm), whereupon the excised tumors were analyzed by histology, immunohistochemistry (IHC) and gene expression profiling (Nanostring PanCancer Immune Profiling panel). The pHIFU + gem treatments were well tolerated; the pHIFU-treated region of the tumor turned hypoechoic immediately following treatment in all mice, and this effect persisted throughout the observation period (2-5 weeks) and corresponded to areas of cell death, according to histology and IHC. Enhanced labeling by Granzyme-B was observed within and adjacent to the pHIFU treated area, but not in the non-treated tumor tissue; no difference in CD8 + staining was observed between the treatment groups. Gene expression analysis showed that the pHIFU + gem combination treatment lead to significant downregulation of 162 genes related to immunosuppression, tumorigenesis, and chemoresistance vs gem only treatment.

Initial Assessment of Boiling Histotripsy for Mechanical Ablation of Ex Vivo Human Prostate Tissue.

Boiling histotripsy (BH) is a focused ultrasound technology that uses millisecond-long pulses with shock fronts to induce mechanical tissue ablation. The pulsing scheme and mechanisms of BH differ from those of cavitation cloud histotripsy, which was previously developed for benign prostatic hyperplasia. The goal of the work described here was to evaluate the feasibility of using BH to ablate fresh ex vivo human prostate tissue as a proof of principle for developing BH for prostate applications. Fresh human prostate samples (N = 24) were obtained via rapid autopsy (<24 h after death, institutional review board exempt). Samples were analyzed using shear wave elastography to ensure that mechanical properties of autopsy tissue were clinically representative. Samples were exposed to BH using 10- or 1-ms pulses with 1% duty cycle under real-time B-mode and Doppler imaging. Volumetric lesions were created by sonicating 1-4 rectangular planes spaced 1 mm apart, containing a grid of foci spaced 1-2 mm apart. Tissue then was evaluated grossly and histologically, and the lesion content was analyzed using transmission electron microscopy and scanning electron microscopy. Observed shear wave elastography characterization of ex vivo prostate tissue (37.9 ± 22.2 kPa) was within the typical range observed clinically. During BH, hyperechoic regions were visualized at the focus on B-mode, and BH-induced bubbles were also detected using power Doppler. As treatment progressed, hypoechoic regions of tissue appeared, suggesting successful tissue fractionation. BH treatment was twofold faster using shorter pulses (1 ms vs. 10 ms). Histological analysis revealed lesions containing completely homogenized cell debris, consistent with histotripsy-induced mechanical ablation. It was therefore determined that BH is feasible in fresh ex vivo human prostate tissue producing desired mechanical ablation. The study supports further work aimed at translating BH technology as a clinical option for prostate ablation.

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.

In Vivo Aberration Correction for Transcutaneous HIFU Therapy Using a Multielement Array.

One of the challenges of transcutaneous high-intensity focused ultrasound (HIFU) therapies, especially ones relying heavily on shock formation, such as boiling histotripsy (BH), is the loss of focusing from aberration induced by the heterogeneities of the body wall. Here, a methodology to execute aberration correction in vivo is proposed. A custom BH system consisting of a 1.5-MHz phased array of 256 elements connected to a Verasonics V1 system is used in pulse/echo mode on a porcine model under general anesthesia. Estimation of the time shifts needed to correct for aberration in the liver and kidney is done by maximizing the value of the coherence factor on the acquired backscattered signals. As this process requires multiple pulse/echo sequences on a moving target to converge to a solution, tracking is also implemented to ensure that the same target is used between each iteration. The method was validated by comparing the acoustic power needed to generate a boiling bubble at one target with aberration correction and at another target within a 5-mm radius without aberration correction. Results show that the aberration correction effectively lowers the acoustic power required to reach boiling by up to 45%, confirming that it indeed restored formation of the nonlinear shock front at the focus.

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.

Development of Tough Hydrogel Phantoms to Mimic Fibrous Tissue for Focused Ultrasound Therapies.

Tissue-mimicking gels provide a cost-effective medium to optimize histotripsy treatment parameters with immediate feedback. Agarose and polyacrylamide gels are often used to evaluate treatment outcomes as they mimic the acoustic properties and stiffness of a variety of soft tissues, but they do not exhibit high toughness, a characteristic of fibrous connective tissue. To mimic pathologic fibrous tissue found in benign prostate hyperplasia (BPH) and other diseases that are potentially treatable with histotripsy, an optically transparent hydrogel with high toughness was developed that is a hybrid of polyacrylamide and alginate. The stiffness was established using shear wave elastography (SWE) and indentometry techniques and was found to be representative of human BPH ex vivo prostate tissue. Different phantom compositions and excised ex vivo BPH tissue samples were treated with a 700-kHz histotripsy transducer at different pulse repetition frequencies. Post-treatment, the hybrid gels and the tissue samples exhibited differential reduction in stiffness as measured by SWE. On B-mode ultrasound, partially treated areas were present as hyperechoic zones and fully liquified areas as hypoechoic zones. Phase contrast microscopy of the gel samples revealed liquefaction in regions consistent with the target lesion dimensions and correlated to findings identified in tissue samples via histology. The dose required to achieve liquefaction in the hybrid gel was similar to what has been observed in ex vivo tissue and greater than that of agarose of comparable or higher Young's modulus by a factor >10. These results indicate that the developed hydrogels closely mimic elasticities found in BPH prostate ex vivo tissue and have a similar response to histotripsy treatment, thus making them a useful cost-effective alternative for developing and evaluating different treatment protocols.

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.

Inertial Cavitation Behaviors Induced by Nonlinear Focused Ultrasound Pulses.

Inertial cavitation induced by pulsed high-intensity focused ultrasound (pHIFU) has previously been shown to successfully permeabilize tumor tissue and enhance chemotherapeutic drug uptake. In addition to HIFU frequency, peak rarefactional pressure ( p- ), and pulse duration, the threshold for cavitation-induced bioeffects has recently been correlated with asymmetric distortion caused by nonlinear propagation, diffraction and formation of shocks in the focal waveform, and therefore with the transducer F -number. To connect previously observed bioeffects with bubble dynamics and their attendant physical mechanisms, the dependence of inertial cavitation behavior on shock formation was investigated in transparent agarose gel phantoms using high-speed photography and passive cavitation detection (PCD). Agarose phantoms with concentrations ranging from 1.5% to 5% were exposed to 1-ms pulses using three transducers of the same aperture but different focal distances ( F -numbers of 0.77, 1.02, and 1.52). Pulses had central frequencies of 1, 1.5, or 1.9 MHz and a range of p- at the focus varying within 1-18 MPa. Three distinct categories of bubble behavior were observed as the acoustic power increased: stationary near-spherical oscillation of individual bubbles, proliferation of multiple bubbles along the pHIFU beam axis, and fanned-out proliferation toward the transducer. Proliferating bubbles were only observed under strongly nonlinear or shock-forming conditions regardless of frequency, and only where the bubbles reached a certain threshold size range. In stiffer gels with higher agarose concentrations, the same pattern of cavitation behavior was observed, but the dimensions of proliferating clouds were smaller. These observations suggest mechanisms that may be involved in bubble proliferation: enhanced growth of bubbles under shock-forming conditions, subsequent shock scattering from the gel-bubble interface, causing an increase in the repetitive tension created by the acoustic wave, and the appearance of a new growing bubble in the proximal direction. Different behaviors corresponded to specific spectral characteristics in the PCD signals: broadband noise in all cases, narrow peaks of backscattered harmonics in the case of stationary bubbles, and broadened, shifted harmonic peaks in the case of proliferating bubbles. The shift in harmonic peaks can be interpreted as a Doppler shift from targets moving at speeds of up to 2 m/s, which correspond to the observed bubble proliferation speeds.

Partial Respiratory Motion Compensation for Abdominal Extracorporeal Boiling Histotripsy Treatments With a Robotic Arm.

Extracorporeal boiling histotripsy (BH), a noninvasive method for mechanical tissue disintegration, is getting closer to clinical applications. However, the motion of the targeted organs, mostly resulting from the respiratory motion, reduces the efficiency of the treatment. Here, a practical and affordable unidirectional respiratory motion compensation method for BH is proposed and evaluated in ex vivo tissues. The BH transducer is fixed on a robotic arm following the motion of the skin, which is tracked using an inline ultrasound imaging probe. In order to compensate for system lags and obtain a more accurate compensation, an autoregressive motion prediction model is implemented. BH pulse gating is also implemented to ensure targeting accuracy. The system is then evaluated with ex vivo BH treatments of tissue samples undergoing motion simulating breathing with the movement of amplitudes between 5 and 10 mm, the frequency between 16 and 18 breaths/min, and a maximum speed of 14.2 mm/s. Results show a reduction of at least 89% of the value of the targeting error during treatment while only increasing the treatment time by no more than 1%. The lesions obtained by treating with the motion compensation were close in size and affected area to the no-motion case, whereas lesions obtained without the compensation were often incomplete and had larger affected areas. This approach to motion compensation could benefit extracorporeal BH and other histotripsy methods in clinical translation.