Objective
In the context of developing boiling histotripsy (BH) as a potential clinical approach for non-invasive mechanical ablation of kidney tumors, the concept of BH dose (BHD) was quantitatively investigated in porcine and canine kidney models in vivo and ex vivo.

Methods
Volumetric lesions were produced in renal tissue using a 1.5-MHz 256-element HIFU-array with various pulsing protocols: pulse duration t_p = 1–10 ms, number of pulses per point ppp = 1–15. Two BHD metrics were evaluated: BHD1 = ppp, BHD2 = t_p × ppp. Quantitative assessment of lesion completeness was performed by their histological analysis and assignment of damage score to different renal compartments (i.e., cortex, medulla, and sinus). Shear wave elastography (SWE) was used to measure the Young's modulus of renal compartments in vivo vs ex vivo, and before vs after BH treatments.

Results
In vivo tissue required lower BH doses to achieve identical degree of fractionation as compared to ex vivo. Renal cortex (homogeneous, low in collagen) was equal or higher in stiffness than medulla (anisotropic, collagenous), 5.8–12.2 kPa vs 4.7–9.6 kPa, but required lower BH doses to be fully fractionated. Renal sinus (fatty, irregular, with abundant collagenous structures) was significantly softer ex vivo vs in vivo, 4.9–5.1 kPa vs 9.7–15.2 kPa, but was barely damaged in either case with any tested BH protocols. BHD1 was shown to be relevant for planning the treatment of renal cortex (sufficient BHD1 = 5 pulses in vivo and 10 pulses ex vivo), while none of the tested doses resulted in complete fractionation of medulla or sinus. Post-treatment SWE imaging revealed reduction of tissue stiffness ex vivo by 27–58%, increasing with the applied dose, and complete absence of shear waves within in vivo lesions, both indicative of tissue liquefaction.

Conclusion
The results imply that tissue resistance to mechanical fractionation, and hence required BH dose, are not solely determined by tissue stiffness but also depend on its composition and structural arrangement, as well as presence of perfusion. The SWE-derived reduction of tissue stiffness with increasing BH doses correlated with tissue damage score, indicating potential of SWE for post-treatment confirmation of BH lesion completeness.

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.

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.