In vivo real-time cavitation imaging in moving organs

3D ultrasound localization microscopy of the nonhuman primate brain

BACKGROUND

Haemodynamic changes occur in stroke and neurodegenerative diseases. Developing imaging techniques allowing the in vivo visualisation and quantification of cerebral blood flow would help better understand the underlying mechanism of these cerebrovascular diseases.

METHODS

3D ultrasound localization microscopy (ULM) is a recently developed technology that can map the microvasculature of the brain at large depth and has been mainly used until now in rodents. In this study, we tested the feasibility of 3D ULM of the nonhuman primate (NHP) brain with a single 256-channel programmable ultrasound scanner.

FINDINGS

We achieved a highly resolved vascular map of the macaque brain at large depth (down to 3 cm) in presence of craniotomy and durectomy using an 8-MHz multiplexed matrix probe. We were able to distinguish vessels as small as 26.9 μm. We also demonstrated that transcranial imaging of the macaque brain at similar depth was feasible using a 3-MHz probe and achieved a resolution of 60 μm.

INTERPRETATION

This work paves the way to clinical applications of 3D ULM. In particular, transcranial 3D ULM in humans could become a tool for the non-invasive study and monitoring of the brain cerebrovascular changes occurring in neurological diseases.

FUNDING

This work was supported by the New Frontier in Research Fund (NFRFE-2022-00590), by the Canada Foundation for Innovation under grant 38095, by the Natural Sciences and Engineering Research Council of Canada (NSERC) under discovery grant RGPIN-2020-06786, by Brain Canada under grant PSG2019, and by the Canadian Institutes of Health Research (CIHR) under grant PJT-156047 and MPI-452530. Computing support was provided by the Digital Research Alliance of Canada.

Self-Sensing Cavitation Detection for Pulsed Cavitational Ultrasound Therapy

OBJECTIVES

Monitoring cavitation during ultrasound therapy is crucial for assessing the procedure safety and efficacy. This work aims to develop a self-sensing and low-complexity approach for robust cavitation detection in moving organs such as the heart.

METHODS

An analog-to-digital converter was connected onto one channel of the therapeutic transducer from a clinical system dedicated to cardiac therapy, allowing to record signals on a computer. Acquisition of successive echoes backscattered by the cavitation cloud on the therapeutic transducer was performed at a high repetition rate. Temporal variations of the backscattered echoes were analyzed with a Singular-Value Decomposition filter to discriminate signals associated to cavitation, based on its stochastic nature. Metrics were derived to classify the filtered backscattered echoes. Classification of raw backscattered echoes was also performed with a machine learning approach. The performances were evaluated on 155 in vitro acquisitions and 110 signals acquired in vivo during transthoracic cardiac ultrasound therapy on 3 swine.

RESULTS

Cavitation detection was achieved successfully in moving tissues with high signal to noise ratio in vitro (cSNR = 25±5) and in vivo (cSNR = 20±6) and outperformed conventional methods (cSNR = 11±6). Classification methods were validated with spectral analysis of hydrophone measurements. High accuracy was obtained using either the clutter filter-based method (accuracy of 1) or the neural network-based method (accuracy of 0.99).

CONCLUSION

Robust self-sensing cavitation detection was demonstrated to be possible with a clutter filter-based method and a machine learning approach.

SIGNIFICANCE

The self-sensing cavitation detection method enables robust, reliable and low complexity cavitation activity monitoring during ultrasound therapy.

Dynamic mode decomposition based Doppler monitoring of de novo cavitation induced by pulsed HIFU: an in vivo feasibility study

Pulsed high-intensity focused ultrasound (pHIFU) has the capability to induce de novo cavitation bubbles, offering potential applications for enhancing drug delivery and modulating tissue microenvironments. However, imaging and monitoring these cavitation bubbles during the treatment presents a challenge due to their transient nature immediately following pHIFU pulses. A planewave bubble Doppler technique demonstrated its potential, yet this Doppler technique used conventional clutter filter that was originally designed for blood flow imaging. Our recent study introduced a new approach employing dynamic mode decomposition (DMD) to address this in an ex vivo setting. This study demonstrates the feasibility of the application of DMD for in vivo Doppler monitoring of the cavitation bubbles in porcine liver and identifies the candidate monitoring metrics for pHIFU treatment. We propose a fully automated bubble mode identification method using k-means clustering and an image contrast-based algorithm, leading to the generation of DMD-filtered bubble images and corresponding Doppler power maps after each HIFU pulse. These power Doppler maps are then correlated with the extent of tissue damage determined by histological analysis. The results indicate that DMD-enhanced power Doppler map can effectively visualize the bubble distribution with high contrast, and the Doppler power level correlates with the severity of tissue damage by cavitation. Further, the temporal characteristics of the bubble modes, specifically the decay rates derived from DMD, provide information of the bubble dissolution rate, which are correlated with tissue damage level-slower rates imply more severe tissue damage.

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.

Comparative Study of Histotripsy Pulse Parameters Used to Inactivate Escherichia coli in Suspension

OBJECTIVE

Bacterial loads can be effectively reduced using cavitation-mediated focused ultrasound, or histotripsy. In this study, gram-negative bacteria (Escherichia coli) in suspension were used as model bacteria to evaluate the effectiveness of two regimens of histotripsy treatments: cavitation histotripsy (CH) and boiling histotripsy (BH).

METHODS

Ten-milliliter volumes of Escherichia coli were treated at different negative focal pressure amplitudes and over time periods up to 40 min. Cavitation activity was characterized with coaxial passive cavitation detection (PCD) and synchronized plane wave B-mode imaging.

RESULTS

CH treatments exhibited a threshold behavior that was consistent with PCD metrics of cavitation. Above the threshold, bacterial inactivation followed a monotonically increasing log-linear relationship that indicated an exponential inactivation rate. BH exhibited no threshold, but instead followed a different monotonically increasing inactivation rate. Inactivation rates were larger for BH at or below the CH threshold, and larger for CH substantially above the threshold. CH studies performed at different pulse lengths at the same duty cycle had similar inactivation rates, suggesting that at any given pressure amplitude, the "on time" was the most important variable for inactivating E. coli. The maximum inactivation was produced by CH at the highest pressure amplitudes used, leading to a log reduction >4.2 for a 40 min treatment.

CONCLUSION

The results of this study suggest that both CH and BH can be used to inactivate E. coli in suspension, with the optimal regimen depending on the attainable peak negative focal pressure at the target.

The fundamental mechanisms of the Korotkoff sounds generation

Blood pressure measurement is the most widely performed clinical exam to predict mortality risk. The gold standard for its noninvasive assessment is the auscultatory method, which relies on listening to the so-called "Korotkoff sounds" in a stethoscope placed at the outlet of a pneumatic arm cuff. However, more than a century after their discovery, the origin of these sounds is still debated, which implies a number of clinical limitations. We imaged the Korotkoff sound generation in vivo at thousands of images per second using ultrafast ultrasound. We showed with both experience and theory that Korotkoff sounds are paradoxically not sound waves emerging from the brachial artery but rather shear vibrations conveyed in surrounding tissues by the nonlinear pulse wave propagation. When these shear vibrations reached the stethoscope, they were synchronous, correlated, and comparable in intensity with the Korotkoff sounds. Understanding this mechanism could ultimately improve blood pressure measurement and provide additional understanding of arterial mechanical properties.

Spatially segmented SVD clutter filtering in cardiac blood flow imaging with diverging waves

Ultrafast ultrasound imaging enables the visualization of rapidly changing blood flow dynamics in the chambers of the heart. Singular value decomposition (SVD) filters outperform conventional high pass clutter rejection filters for ultrafast blood flow imaging of small and shallow fields of view (e.g., functional imaging of brain activity). However, implementing SVD filters can be challenging in cardiac imaging due to the complex spatially and temporally varying tissue characteristics. To address this challenge, we describe a method that involves excluding the proximal portion of the image (near the chest wall) and divides the reduced field of view into overlapped segments, within which tissue signals are expected to be spatially and temporally coherent. SVD filtering with automatic selection of cut-off singular vector orders to remove tissue and noise signals is implemented for each segment. Auto-thresholding is based on the coherence of spatial singular vectors, delineating tissue, blood, and noise subspaces within a spatial similarity matrix calculated for each segment. Filtered blood flow signals from the segments are reconstructed and then combined and Doppler processing is used to form a set of blood flow images. Preliminary experimental results suggest that the spatially segmented approach improves the separation of the tissue and blood subsets in the spatial similarity matrix so that automatic thresholding is significantly improved, and tissue clutter can then be rejected more effectively in cardiac ultrafast imaging, compared to using the full field of view. In the case studied, spatially segmented SVD improved the rate of correct automatic selection of thresholds from 78% to 98.7% for the investigated cases and improved the post-filter power of blood signals by an average of more than 10 dB during a cardiac cycle.

Frame rate effects and their compensation on super-resolution microvessel imaging using ultrasound localization microscopy

Ultrasound localization microscopy (ULM) breaks the diffraction limit and allows imaging microvasculature at micrometric resolution while preserving the penetration depth. Frame rate plays an important role for high-quality ULM imaging, but there is still a lack of review and investigation of the frame rate effects on ULM. This work aims to clarify how frame rate influences the performance of ULM, including the effects of microbubble detection, localization and tracking. The performance of ULM was evaluated using an in vivo rat brain dataset (15.6 MHz, 3 tilted plane waves (-5°, 0°, +5°), at a compounded frame rate of 1000 Hz) with different frame rates. Quantification methods, including Fourier ring correlation and saturation parameter, were applied to analyze the spatial resolution and reconstruction efficiency, respectively. In addition, effects on each crucial step in ULM processing were further analyzed. Results showed that when frame rates dropped from 1000 Hz to 250 Hz, the spatial resolution deteriorated from 9.9 μm to 15.0 μm. Applying a velocity constraint was able to improve the ULM performance, but inappropriate constraint may artificially result in high apparent resolution. For the dataset, compared with the results of 1000 Hz frame rate, the velocity was underestimated at 100 Hz with 47.18% difference and the saturation was reduced from 55.00% at 1000 Hz to 43.34% at 100 Hz. Analysis showed that inadequate frame rate generated unreliable microbubble detection, localization and tracking as well as incomplete track reconstruction, resulting in the deterioration in spatial resolution, the underestimation in velocity measurement and the decrease in saturation. Finally, a guidance of determining the frame rate requirement was discussed by considering the required spatial sampling points based on vessel morphology, clutter filtering method, tracking algorithm and acquisition time, which provides indications for future clinical application of ULM method.

A simplified and accelerated implementation of SVD for filtering ultrafast power Doppler images

BACKGROUND AND OBJECTIVE

Ultrafast Power Doppler (UPD) is a growing ultrasound modality for imaging and diagnosing microvasculature disease. A key element of UPD is using singular value decomposition (SVD) as a highly selective filter for tissue and electronic noise. However, two significant drawbacks of SVD are its computational burden and the complexity of its algorithms. These limitations hinder the development of fast and specific SVD algorithms for UPD imaging. This study introduces power SVD (pSVD), a simplified and accelerated algorithm for filtering tissue and noise in UPD images.

METHODS

pSVD exploits several mathematical properties of SVD specific to UPD images. In particular, pSVD allows the direct computation of blood-related SVD components from the temporal singular vectors. This feature simplifies the expression of SVD while significantly accelerating its computation. After detailing the theory behind pSVD, we evaluate its performances in several in vitro and in vivo experiments and compare it to SVD and randomized SVD (rSVD).

RESULTS

pSVD strongly decreases the running time of SVD (between 5 and 12 times in vivo) without impacting the quality of UPD images. Compared to rSVD, pSVD can be significantly faster (up to 3 times) or slightly slower but gives access to more estimators to isolate tissue subspaces.

CONCLUSION

pSVD is highly valuable for implementing UPD imaging in clinical ultrasound and provides a better understanding of SVD for ultrasound imaging in general.

Super-Resolution Ultrasound Localization Microscopy Using High-Frequency Ultrasound to Measure Ocular Perfusion Velocity in the Rat Eye

Imaging of the ocular vasculature can provide new insights into the pathophysiology of ocular diseases. This study proposes a novel high-frequency super-resolution ultrasound localization microscopy (SRULM) technique and evaluates its ability to measure in vivo perfusion changes in the rat eye at elevated intraocular pressure (IOP). A 38.4 MHz center frequency linear array transducer on a VisualSonics Vevo F2 imaging platform was used to collect high frame rate (1 kHz) radiofrequency data of the posterior rat eye following systemic microbubble contrast injection. Following clutter and spatiotemporal non-local means filtering, individual microbubbles were localized and tracked. The microbubble tracks were accumulated over 10,000 frames to generate vascular images quantifying perfusion velocity and direction. Experiments were performed using physiologic relevant controlled flow states for algorithm validation and subsequently performed in vivo on the rat eye at 10 mm Hg IOP increments from 10 to 60 mm Hg. The posterior vasculature of the rat eye, including the ophthalmic artery, long posterior ciliary arteries and their branches, central retinal artery and retinal arterioles and venules were successfully visualized, and velocities quantified at each IOP level. Significant reductions in arterial flow were measured as IOP was elevated. High-frequency SRULM can be used to visualize and quantify the perfusion velocity of the rat eye in both the retrobulbar and intraocular vasculature simultaneously. The ability to detect ocular perfusion changes throughout the depth of the eye may help elucidate the role ischemia has in the pathophysiology of ocular diseases such as glaucoma.

Barrier-breaking effects of ultrasonic cavitation for drug delivery and biomarker release

Recently, emerging evidence has demonstrated that cavitation actually creates important bidirectional channels on biological barriers for both intratumoral drug delivery and extratumoral biomarker release. To promote the barrier-breaking effects of cavitation for both therapy and diagnosis, we first reviewed recent technical advances of ultrasound and its contrast agents (microbubbles, nanodroplets, and gas-stabilizing nanoparticles) and then reported the newly-revealed cavitation physical details. In particular, we summarized five types of cellular responses of cavitation in breaking the plasma membrane (membrane retraction, sonoporation, endocytosis/exocytosis, blebbing and apoptosis) and compared the vascular cavitation effects of three different types of ultrasound contrast agents in breaking the blood-tumor barrier and tumor microenvironment. Moreover, we highlighted the current achievements of the barrier-breaking effects of cavitation in mediating drug delivery and biomarker release. We emphasized that the precise induction of a specific cavitation effect for barrier-breaking was still challenged by the complex combination of multiple acoustic and non-acoustic cavitation parameters. Therefore, we provided the cutting-edge in-situ cavitation imaging and feedback control methods and suggested the development of an international cavitation quantification standard for the clinical guidance of cavitation-mediated barrier-breaking effects.

Independent Component Analysis Filter for Small Vessel Contrast Imaging During Fast Tissue Motion

Suppressing tissue clutter is an essential step in blood flow estimation and visualization, even when using ultrasound contrast agents. Blind source separation (BSS)-based clutter filter for high-framerate ultrasound imaging has been reported to perform better in tissue clutter suppression than the conventional frequency-based wall filter and nonlinear contrast pulsing schemes. The most notable BSS technique, singular value decomposition (SVD) has shown compelling results in cases of slow tissue motion. However, its performance degrades when the tissue motion is faster than the blood flow speed, conditions that are likely to occur when imaging the small vessels, such as in the myocardium. Independent component analysis (ICA) is another BSS technique that has been implemented as a clutter filter in the spatiotemporal domain. Instead, we propose to implement ICA in the spatial domain where motion should have less impact. In this work, we propose a clutter filter with the combination of SVD and ICA to improve the contrast-to-background ratio (CBR) in cases where tissue velocity is significantly faster than the flow speed. In an in vitro study, the range of fast tissue motion velocity was 5-25 mm/s and the range of flow speed was 1-12 mm/s. Our results show that the combination of ICA and SVD yields 7-10 dB higher CBR than SVD alone, especially in the tissue high-velocity range. The improvement is crucial for cardiac imaging where relatively fast myocardial motions are expected.

Real-Time Visualization of a Focused Ultrasound Beam Using Ultrasonic Backscatter

Focused ultrasound (FUS) therapies induce therapeutic effects in localized tissues using either temperature elevations or mechanical stresses caused by an ultrasound wave. During an FUS therapy, it is crucial to continuously monitor the position of the FUS beam in order to correct for tissue motion and keep the focus within the target region. Toward the goal of achieving real-time monitoring for FUS therapies, we have developed a method for the real-time visualization of an FUS beam using ultrasonic backscatter. The intensity field of an FUS beam was reconstructed using backscatter from an FUS pulse received by an imaging array and then overlaid onto a B-mode image captured using the same imaging array. The FUS beam visualization allows one to monitor the position and extent of the FUS beam in the context of the surrounding medium. Variations in the scattering properties of the medium were corrected in the FUS beam reconstruction by normalizing based on the echogenicity of the coaligned B-mode image. On average, normalizing by echogenicity reduced the mean square error between FUS beam reconstructions in nonhomogeneous regions of a phantom and baseline homogeneous regions by 21.61. FUS beam visualizations were achieved, using a single diagnostic imaging array as both an FUS source and an imaging probe, in a tissue-mimicking phantom and a rat tumor in vivo with a frame rate of 25-30 frames/s.

Blind Source Separation for Clutter and Noise Suppression in Ultrasound Imaging: Review for Different Applications

Blind source separation (BSS) refers to a number of signal processing techniques that decompose a signal into several "source" signals. In recent years, BSS is increasingly employed for the suppression of clutter and noise in ultrasonic imaging. In particular, its ability to separate sources based on measures of independence rather than their temporal or spatial frequency content makes BSS a powerful filtering tool for data in which the desired and undesired signals overlap in the spectral domain. The purpose of this work was to review the existing BSS methods and their potential in ultrasound imaging. Furthermore, we tested and compared the effectiveness of these techniques in the field of contrast-ultrasound super-resolution, contrast quantification, and speckle tracking. For all applications, this was done in silico, in vitro, and in vivo. We found that the critical step in BSS filtering is the identification of components containing the desired signal and highlighted the value of a priori domain knowledge to define effective criteria for signal component selection.

Cavitation in soft matter

Cavitation is the sudden, unstable expansion of a void or bubble within a liquid or solid subjected to a negative hydrostatic stress. Cavitation rheology is a field emerging from the development of a suite of materials characterization, damage quantification, and therapeutic techniques that exploit the physical principles of cavitation. Cavitation rheology is inherently complex and broad in scope with wide-ranging applications in the biology, chemistry, materials, and mechanics communities. This perspective aims to drive collaboration among these communities and guide discussion by defining a common core of high-priority goals while highlighting emerging opportunities in the field of cavitation rheology. A brief overview of the mechanics and dynamics of cavitation in soft matter is presented. This overview is followed by a discussion of the overarching goals of cavitation rheology and an overview of common experimental techniques. The larger unmet needs and challenges of cavitation in soft matter are then presented alongside specific opportunities for researchers from different disciplines to contribute to the field.

Observation and modulation of the dissolution of histotripsy-induced bubble clouds with high-frame rate plane wave imaging

Focused ultrasound therapies are a noninvasive means to ablate tissue. Histotripsy utilizes short ultrasound pulses with sufficient tension to nucleate bubble clouds that impart lethal strain to the surrounding tissues. Tracking bubble cloud dissolution between the application of histotripsy pulses is critical to ensure treatment efficacy. In this study, plane wave B-mode imaging was employed to monitor bubble cloud motion and grayscale at frame rates up to 11.25 kHz. Minimal changes in the area or position of the bubble clouds were observed 50 ms post excitation. The bubble cloud grayscale was observed to decrease with the square root of time, indicating a diffusion-driven process. These results were qualitatively consistent with an analytic model of gas diffusion during the histotripsy process. Finally, the rate of bubble cloud dissolution was found to be dependent on the output of the imaging pulse, indicating an interaction between the bubble cloud and imaging parameters. Overall, these results highlight the utility of plane wave B-mode imaging for monitoring histotripsy bubble clouds.

SVD-Based Separation of Stable and Inertial Cavitation Signals Applied to Passive Cavitation Mapping During HIFU

Detection of inertial and stable cavitation is important for guiding high-intensity focused ultrasound (HIFU). Acoustic transducers can passively detect broadband noise from inertial cavitation and the scattering of HIFU harmonics from stable cavitation bubbles. Conventional approaches to cavitation noise diagnostics typically involve computing the Fourier transform of the time-domain noise signal, applying a custom comb filter to isolate the frequency components of interest, followed by an inverse Fourier transform. We present an alternative technique based on singular value decomposition (SVD) that efficiently separates the broadband emissions and HIFU harmonics. Spatiotemporally resolved cavitation detection was achieved using a 128-element, 5-MHz linear-array ultrasound imaging system operating in the receive mode at 15 frames/s. A 1.1-MHz transducer delivered HIFU to tissue-mimicking phantoms and excised liver tissue for a duration of 5 s. Beamformed radio frequency signals corresponding to each scan line in a frame were assembled into a matrix, and SVD was performed. Spectra of the singular vectors obtained from a tissue-mimicking gel phantom were analyzed by computing the peak ratio ( R ), defined as the ratio of the peak of its fifth-order polynomial fit and the maximum spectral peak. Singular vectors that produced an were classified as those representing stable cavitation, i.e., predominantly containing harmonics of HIFU. The projection of data onto this singular base reproduced stable cavitation signals. Similarly, singular vectors that produced an were classified as those predominantly containing broadband noise associated with inertial cavitation. These singular vectors were used to isolate the inertial cavitation signal. The R -value thresholds determined using gel data were then employed to analyze cavitation data obtained from bovine liver ex vivo. The SVD-based method faithfully reproduced the structural details in the spatiotemporal cavitation maps produced using the more cumbersome comb-filter approach with a maximum root-mean-squared error of 10%.

MRI-guided transurethral insonation of silica-shell phase-shift emulsions in the prostate with an advanced navigation platform

PURPOSE

In this study, the efficacy of transurethral prostate ablation in the presence of silica-shell ultrasound-triggered phase-shift emulsions (sUPEs) doped with MR contrast was evaluated. The influence of sUPEs on MR imaging assessment of the ablation zone was also investigated.

METHODS

sUPEs were doped with a magnetic resonance (MR) contrast agent, Gd2 O3 , to assess ultrasound transition. Injections of saline (sham), saline and sUPEs alone, and saline and sUPEs with Optison microbubbles were performed under guidance of a prototype interventional MRI navigation platform in a healthy canine prostate. Treatment arms were evaluated for differences in lesion size, T1  contrast, and temperature. In addition, non-perfused areas (NPAs) on dynamic contrast-enhanced (DCE) MRI, 55°C isotherms, and areas of 240 cumulative equivalent minutes at 43°C (CEM43 ) dose or greater computed from MR thermometry were measured and correlated with ablated areas indicated by histology.

RESULTS

For treatment arms including sUPEs, the computed correlation coefficients between the histological ablation zone and the NPA, 55°C isotherm, and 240 CEM43 area ranged from 0.96-0.99, 0.98-0.99, and 0.91-0.99, respectively. In the absence of sUPEs, the computed correlation coefficients between the histological ablation zone and the NPA, 55°C isotherm, and 240 CEM43 area were 0.69, 0.54, and 0.50, respectively. Across all treatment arms, the areas of thermal tissue damage and NPAs were not significantly different (P = 0.47). Areas denoted by 55°C isotherms and 240 CEM43 dose boundaries were significantly larger than the areas of thermal damage, again for all treatment arms (P = 0.009 and 0.003, respectively). No significant differences in lesion size, T1 contrast, or temperature were observed between any of the treatment arms (P > 0.0167). Lesions exhibiting thermal fixation on histological analysis were present in six of nine insonations involving sUPE injections and one of five insonations involving saline sham injections. Significantly larger areas (P = 0.002), higher temperatures (P = 0.004), and more frequent ring patterns of restricted diffusion on ex vivo diffusion-weighted imaging (P = 0.005) were apparent in lesions with thermal fixation.

CONCLUSIONS

T1 contrast suggesting sUPE transition was not evident in sUPE treatment arms. The use of MR imaging metrics to predict prostate ablation was not diminished by the presence of sUPEs. Lesions generated in the presence of sUPEs exhibited more frequent thermal fixation, though there were no significant changes in the ablation areas when comparing arms with and without sUPEs. Thermal fixation corresponded to some qualitative imaging features.

2D and 3D real-time passive cavitation imaging of pulsed cavitation ultrasound therapy in moving tissues

Pulsed cavitation ultrasound therapy (PCUT) is an effective non-invasive therapeutic approach in various medical indications that relies on the mechanical effects generated by cavitation bubbles. Even though limited by the poor contrast, conventional ultrasound B-Mode imaging has been widely used for the guidance and monitoring of the therapeutic procedure, allowing the visualization of the cavitation bubble cloud. However, the visualization of the bubble cloud is often limited in deep organs such as the liver and the heart and remains moreover completely subjective for the operator. Our goal is to develop a new imaging mode to better identify the cavitation cloud. Active and passive cavitation imaging methods have been developed but none of them has been able to locate the cavitation bubble created by PCUT in real-time and in moving organs. In this paper we propose a passive ultrasound imaging approach combined with a spatiotemporal singular value decomposition filter to detect and map the bubble cloud with high sensitivity and high contrast. In moving applications at a maximal motion speed of 10 mm s^-1, the contrast-to-noise ratio for passive cavitation imaging is up to 10 times higher than for active cavitation imaging, with a temporal resolution of about 100 ms. The mapping of the bubble cloud can be overlaid in real-time to the conventional B-Mode, which permits to locate the cavitation phenomena in relation to the anatomic image. Finally, we extend the technique to volumetric imaging and show its feasibility on moving phantoms.

Adaptive Spatiotemporal SVD Clutter Filtering for Ultrafast Doppler Imaging Using Similarity of Spatial Singular Vectors

Singular value decomposition of ultrafast imaging ultrasonic data sets has recently been shown to build a vector basis far more adapted to the discrimination of tissue and blood flow than the classical Fourier basis, improving by large factor clutter filtering and blood flow estimation. However, the question of optimally estimating the boundary between the tissue subspace and the blood flow subspace remained unanswered. Here, we introduce an efficient estimator for automatic thresholding of subspaces and compare it to an exhaustive list of thirteen estimators that could achieve this task based on the main characteristics of the singular components, namely the singular values, the temporal singular vectors, and the spatial singular vectors. The performance of those fourteen estimators was tested in vitro in a large set of controlled experimental conditions with different tissue motion and flow speeds on a phantom. The estimator based on the degree of resemblance of spatial singular vectors outperformed all others. Apart from solving the thresholding problem, the additional benefit with this estimator was its denoising capabilities, strongly increasing the contrast to noise ratio and lowering the noise floor by at least 5 dB. This confirms that, contrary to conventional clutter filtering techniques that are almost exclusively based on temporal characteristics, efficient clutter filtering of ultrafast Doppler imaging cannot overlook space. Finally, this estimator was applied in vivo on various organs (human brain, kidney, carotid, and thyroid) and showed efficient clutter filtering and noise suppression, improving largely the dynamic range of the obtained ultrafast power Doppler images.

Power cavitation-guided blood-brain barrier opening with focused ultrasound and microbubbles

Image-guided monitoring of microbubble-based focused ultrasound (FUS) therapies relies on the accurate localization of FUS-stimulated microbubble activity (i.e. acoustic cavitation). Passive cavitation imaging with ultrasound arrays can achieve this, but with insufficient spatial resolution. In this study, we address this limitation and perform high-resolution monitoring of acoustic cavitation-mediated blood-brain barrier (BBB) opening with a new technique called power cavitation imaging. By synchronizing the FUS transmit and passive receive acquisition, high-resolution passive cavitation imaging was achieved by using delay and sum beamforming with absolute time delays. Since the axial image resolution is now dependent on the duration of the received acoustic cavitation emission, short pulses of FUS were used to limit its duration. Image sets were acquired at high-frame rates for calculation of power cavitation images analogous to power Doppler imaging. Power cavitation imaging displays the mean intensity of acoustic cavitation over time and was correlated with areas of acoustic cavitation-induced BBB opening. Power cavitation-guided BBB opening with FUS could constitute a standalone system that may not require MRI guidance during the procedure. The same technique can be used for other acoustic cavitation-based FUS therapies, for both safety and guidance.

Post Hoc Analysis of Passive Cavitation Imaging for Classification of Histotripsy-Induced Liquefaction in Vitro

Histotripsy utilizes focused ultrasound to generate bubble clouds for transcutaneous tissue liquefaction. Bubble activity maps are under development to provide image guidance and monitor treatment progress. The aim of this paper was to investigate the feasibility of using plane wave B-mode and passive cavitation images to be used as binary classifiers of histotripsy-induced liquefaction. Prostate tissue phantoms were exposed to histotripsy pulses over a range of pulse durations (5- ) and peak negative pressures (12-23 MPa). Acoustic emissions were recorded during the insonation and beamformed to form passive cavitation images. Plane wave B-mode images were acquired following the insonation to detect the hyperechoic bubble cloud. Phantom samples were sectioned and stained to delineate the liquefaction zone. Correlation between passive cavitation and plane wave B-mode images and the liquefaction zone was assessed using receiver operating characteristic (ROC) curve analysis. Liquefaction of the phantom was observed for all the insonation conditions. The area under the ROC (0.94 versus 0.82), accuracy (0.90 versus 0.83), and sensitivity (0.81 versus 0.49) was greater for passive cavitation images relative to B-mode images ( ) along the azimuth of the liquefaction zone. The specificity was greater than 0.9 for both imaging modalities. These results demonstrate a stronger correlation between histotripsy-induced liquefaction and passive cavitation imaging compared with the plane wave B-mode imaging, albeit with limited passive cavitation image range resolution.

Gold nanoparticle nucleated cavitation for enhanced high intensity focused ultrasound therapy

High intensity focused ultrasound (HIFU) or focused ultrasound surgery is a non-invasive technique for the treatment of cancerous tissue, which is limited by difficulties in getting real-time feedback on treatment progress and long treatment durations. The formation and activity of acoustic cavitation, specifically inertial cavitation, during HIFU exposures has been demonstrated to enhance heating rates. However, without the introduction of external nuclei its formation an activity can be unpredictable, and potentially counter-productive. In this study, a combination of pulse laser illumination (839 nm), HIFU exposures (3.3 MHz) and plasmonic gold nanorods (AuNR) was demonstrated as a new approach for the guidance and enhancement of HIFU treatments. For imaging, short duration HIFU pulses (10 μs) demonstrated broadband acoustic emissions from AuNR nucleated cavitation with a signal-to-noise ranging from 5-35 dB for peak negative pressures between 1.19-3.19  ±  0.01 MPa. In the absence of either AuNR or laser illumination these emissions were either not present or lower in magnitude (e.g. 5 dB for 3.19 MPa). Continuous wave (CW) HIFU exposures for 15 s, were then used to generate thermal lesions for peak negative pressures from 0.2-2.71  ±  0.01 MPa at a fluence of 3.4 mJ [Formula: see text]. Inertial cavitation dose (ICD) was monitored during all CW exposures, where exposures combined with both laser illumination and AuNRs resulted in the highest level of detectable emissions. This parameter was integrated over the entire exposure to give a metric to compare with measured thermal lesion area, where it was found that a minimum total ICD of [Formula: see text] a.u. was correlated with the formation of thermal lesions in gel phantoms. Furthermore, lesion area (mm²) was increased for equivalent exposures without either AuNRs or laser illumination. Once combined with cancer targeting AuNRs this approach could allow for the future theranostic use of HIFU, such as providing the ability to identify and treat small multi-focal cancerous regions with minimal damage to surrounding healthy tissue.