Combined ultrasonic thermal ablation with interleaved ARFI image monitoring using a single diagnostic curvilinear array: a feasibility study

Technical note: Evaluation of the acoustic radiation force imaging for predicting HIFU focus with in vitro and ex vivo experiments

BACKGROUND

High-intensity focused ultrasound (HIFU) is currently used for the treatment of various diseases, but it still lacks a reliable technique in the preoperative stage to accurately place its "energy blade" onto diseased targets. Acoustic radiation force imaging (ARFI) was recently introduced to tackle this issue, but its applicability and limitations were not clear.

PURPOSE

The aim of this study was to evaluate the performance of ARFI method in prediction of HIFU focal location at the preoperative stage.

METHODS

A point spread function (PSF) localization method, which was borrowed from the ultrasound super resolution field, was used to validate the core autocorrelation-based motion estimation algorithm in the ARFI procedure. Accuracy of the ARFI method for estimating the HIFU focus were tested with in vitro and ex vivo experiments with a clinically equivalent HIFU system. Comparisons were made between the estimated focal locations and those of the damaged area after the testing objects were cut open.

RESULTS

Results showed that the PSF localization was able to serve as a validating method for motion detection only when the tissue displacement was large. With the ARFI method, location of the HIFU focus could be accurately predicted by a 2D motion map preoperatively, and the axial spatial errors were less than 0.5 mm. However, the derived 2D motion maps can only be valuable when the acoustic stimulation in ARFI were strong enough, which was probably due to the fact that the HIFU focal locations were at large depths and the ultrasound imaging signal had low signal to noise ratio.

CONCLUSION

The ARFI method was indeed an accurate technique for preoperatively predicting HIFU focus in vitro and ex vivo. If clinical applications were to be considered, particularly in deep tissues, efforts might need to be made to improve ability of the ultrasound motion estimation technique.

Methods of monitoring thermal ablation of soft tissue tumors - A comprehensive review

Thermal ablation is a form of hyperthermia in which oncologic control can be achieved by briefly inducing elevated temperatures, typically in the range 50-80°C, within a target tissue. Ablation modalities include high intensity focused ultrasound, radiofrequency ablation, microwave ablation, and laser interstitial thermal therapy which are all capable of generating confined zones of tissue destruction, resulting in fewer complications than conventional cancer therapies. Oncologic control is contingent upon achieving predefined coagulation zones; therefore, intraoperative assessment of treatment progress is highly desirable. Consequently, there is a growing interest in the development of ablation monitoring modalities. The first section of this review presents the mechanism of action and common applications of the primary ablation modalities. The following section outlines the state-of-the-art in thermal dosimetry which includes interstitial thermal probes and radiologic imaging. Both the physical mechanism of measurement and clinical or pre-clinical performance are discussed for each ablation modality. Thermal dosimetry must be coupled with a thermal damage model as outlined in Section 4. These models estimate cell death based on temperature-time history and are inherently tissue specific. In the absence of a reliable thermal model, the utility of thermal monitoring is greatly reduced. The final section of this review paper covers technologies that have been developed to directly assess tissue conditions. These approaches include visualization of non-perfused tissue with contrast-enhanced imaging, assessment of tissue mechanical properties using ultrasound and magnetic resonance elastography, and finally interrogation of tissue optical properties with interstitial probes. In summary, monitoring thermal ablation is critical for consistent clinical success and many promising technologies are under development but an optimal solution has yet to achieve widespread adoption.

HIFU-induced changes in optical scattering and absorption of tissue over nine orders of thermal dose

The optical properties of tissue change during thermal ablation. Multi-modal methods such as acousto-optic (AO) and photo-acoustic (PA) imaging may provide a real-time, direct measure of lesion formation. Baseline changes in optical properties have been previously measured over limited ranges of thermal dose for tissues exposed to a temperature-controlled water bath, however, there is scant data for optical properties of lesions created by HIFU. In this work, the optical scattering and absorption coefficients from 400-1300 nm of excised chicken breast exposed to HIFU were measured using an integrating sphere spectrophotometric technique. HIFU-induced spatiotemporal temperature elevations were measured using an infrared camera and used to calculate the thermal dose delivered to a localized region of tissue. Results obtained over a range of thermal dose spanning 9 orders of magnitude show that the reduced scattering coefficient increases for HIFU exposures exceeding a threshold thermal dose of CEM₄₃  =  600  ±  81 cumulative equivalent minutes. HIFU-induced thermal damage results in changes in scattering over all optical wavelengths, with a 2.5-fold increase for thermal lesions exceeding 70 °C. The tissue absorption coefficient was also found to increase for thermally lesioned tissue, however, the magnitude was strongly dependent on the optical wavelength and there was substantial sample-to-sample variability, such that the existence of a threshold thermal dose could not be determined. Therapeutic windows, where the optical penetration depth is expected to be greatest, were identified in the near infrared regime centered near 900 nm and 1100 nm. These data motivate further research to improve the real-time AO and PA sensing of lesion formation during HIFU therapy as an alternative to thermometry.

Detection of gaps between high-intensity focused ultrasound (HIFU)-induced lesions using transient axial shear strain elastograms

PURPOSE

High-intensity focused ultrasound (HIFU) is becoming an effective and noninvasive treatment modality for cancer and solid tumors. In order to avoid the cancer relapse and guarantee the success of ablation, there should be no gaps left among all HIFU-generated lesions. However, there are few imaging approaches available for detecting the HIFU lesion gaps in real time during ablation.

METHODS

Transient axial shear strain elastograms (ASSEs) were proposed and evaluated both numerically and experimentally to detect the lesion gaps immediately after the cessation of therapeutic HIFU exposure. Acoustic intensity and subsequent acoustic radiation force were first calculated by solving the nonlinear Khokhlov-Zabolotskaya-Kuznetzov (KZK) equation. Motion of being- and already-treated lesions during and after HIFU exposure was simulated using the transient dynamic analysis module of finite element method (FEM). The corresponding B-mode sonography of tissue-mimicking phantom with two HIFU lesions inside was simulated by FIELD II, and then axial strain elastograms (ASEs) under static compression and transient ASSEs were reconstructed. An ultrasound imaging probe was integrated with the HIFU transducer and used to obtain radio frequency (RF) echo signals at high frame rate using plane wave imaging (PWI). The resulting strains were mapped using the correlation-based method and block search strategy.

RESULTS

Acoustic radiation force from the therapeutic HIFU burst is sufficiently strong to produce significant displacement. As a result, large and highly localized axial shear strain appears in the gap zone between two HIFU-generated lesions and then disappears after sufficient HIFU ablation (no gap between them). Such capability of detecting the lesion gap is validated at the varied acoustic radiation force density, gap width, and the size of the lesion. In contrast, conventional ASEs using the static compression cannot distinguish whether a gap exists between lesions. Static ASEs and transient ASSEs reconstructed using both high-speed photography and sonography in the gel phantom show the same conclusion as that in the simulation. Ex vivo tissue experiments further confirmed that the presence of large axial shear strain in the gap zone. The ratios of axial shear strain in the porcine kidney and liver samples had statistical differences for two HIFU-generated lesions without and with a gap (P < 0.05).

CONCLUSIONS

Large axial shear strain induced by the acoustic radiation force from therapeutic HIFU burst only appears between two HIFU-generated lesions with a gap between them. Transient ASSEs reconstructed immediately after the cession of HIFU exposure can easily, reliably, and sensitively detect the gap between produced lesions, which would provide real-time feedback to enhance the success of HIFU ablation.

Delineation of Post-Procedure Ablation Regions with Electrode Displacement Elastography with a Comparison to Acoustic Radiation Force Impulse Imaging

We compared a quasi-static ultrasound elastography technique, referred to as electrode displacement elastography (EDE), with acoustic radiation force impulse imaging (ARFI) for monitoring microwave ablation (MWA) procedures on patients diagnosed with liver neoplasms. Forty-nine patients recruited to this study underwent EDE and ARFI with a Siemens Acuson S2000 system after an MWA procedure. On the basis of visualization results from two observers, the ablated region in ARFI images was recognizable on 20 patients on average in conjunction with B-mode imaging, whereas delineable ablation boundaries could be generated on 4 patients on average. With EDE, the ablated region was delineable on 40 patients on average, with less imaging depth dependence. Study of tissue-mimicking phantoms revealed that the ablation region dimensions measured on EDE and ARFI images were within 8%, whereas the image contrast and contrast-to-noise ratio with EDE was two to three times higher than that obtained with ARFI. This study indicated that EDE provided improved monitoring results for minimally invasive MWA in clinical procedures for liver cancer and metastases.

Post-Procedure Evaluation of Microwave Ablations of Hepatocellular Carcinomas Using Electrode Displacement Elastography

Microwave ablation has been used clinically as an alternative to surgical resection. However, lack of real-time imaging to assess treated regions may compromise treatment outcomes. We previously introduced electrode displacement elastography (EDE) for strain imaging and verified its feasibility in vivo on porcine animal models. In this study, we evaluated EDE on 44 patients diagnosed with hepatocellular carcinoma, treated using microwave ablation. The ablated region was identified on EDE images for 40 of the 44 patients. Ablation areas averaged 13.38 ± 4.99 cm² on EDE, compared with 7.61 ± 3.21 cm² on B-mode imaging. Contrast and contrast-to-noise ratios obtained with EDE were 232% and 98%, respectively, significantly higher than values measured on B-mode images (p < 0.001). This study indicates that EDE is feasible in patients and provides improved visualization of the ablation zone compared with B-mode ultrasound.

Real-time Monitoring of High Intensity Focused Ultrasound (HIFU) Ablation of In Vitro Canine Livers Using Harmonic Motion Imaging for Focused Ultrasound (HMIFU)

Harmonic Motion Imaging for Focused Ultrasound (HMIFU) is a technique that can perform and monitor high-intensity focused ultrasound (HIFU) ablation. An oscillatory motion is generated at the focus of a 93-element and 4.5 MHz center frequency HIFU transducer by applying a 25 Hz amplitude-modulated signal using a function generator. A 64-element and 2.5 MHz imaging transducer with 68kPa peak pressure is confocally placed at the center of the HIFU transducer to acquire the radio-frequency (RF) channel data. In this protocol, real-time monitoring of thermal ablation using HIFU with an acoustic power of 7 W on canine livers in vitro is described. HIFU treatment is applied on the tissue during 2 min and the ablated region is imaged in real-time using diverging or plane wave imaging up to 1,000 frames/second. The matrix of RF channel data is multiplied by a sparse matrix for image reconstruction. The reconstructed field of view is of 90° for diverging wave and 20 mm for plane wave imaging and the data are sampled at 80 MHz. The reconstruction is performed on a Graphical Processing Unit (GPU) in order to image in real-time at a 4.5 display frame rate. 1-D normalized cross-correlation of the reconstructed RF data is used to estimate axial displacements in the focal region. The magnitude of the peak-to-peak displacement at the focal depth decreases during the thermal ablation which denotes stiffening of the tissue due to the formation of a lesion. The displacement signal-to-noise ratio (SNRd) at the focal area for plane wave was 1.4 times higher than for diverging wave showing that plane wave imaging appears to produce better displacement maps quality for HMIFU than diverging wave imaging.

Using passive cavitation images to classify high-intensity focused ultrasound lesions

Passive cavitation imaging provides spatially resolved monitoring of cavitation emissions. However, the diffraction limit of a linear imaging array results in relatively poor range resolution. Poor range resolution has limited prior analyses of the spatial specificity and sensitivity of passive cavitation imaging in predicting thermal lesion formation. In this study, this limitation is overcome by orienting a linear array orthogonal to the high-intensity focused ultrasound propagation direction and performing passive imaging. Fourteen lesions were formed in ex vivo bovine liver samples as a result of 1.1-MHz continuous-wave ultrasound exposure. The lesions were classified as focal, "tadpole" or pre-focal based on their shape and location. Passive cavitation images were beamformed from emissions at the fundamental, harmonic, ultraharmonic and inharmonic frequencies with an established algorithm. Using the area under a receiver operating characteristic curve (AUROC), fundamental, harmonic and ultraharmonic emissions were found to be significant predictors of lesion formation for all lesion types. For both harmonic and ultraharmonic emissions, pre-focal lesions were classified most successfully (AUROC values of 0.87 and 0.88, respectively), followed by tadpole lesions (AUROC values of 0.77 and 0.64, respectively) and focal lesions (AUROC values of 0.65 and 0.60, respectively).

Harmonic motion imaging for abdominal tumor detection and high-intensity focused ultrasound ablation monitoring: an in vivo feasibility study in a transgenic mouse model of pancreatic cancer

Harmonic motion imaging (HMI) is a radiationforce- based elasticity imaging technique that tracks oscillatory tissue displacements induced by sinusoidal ultrasonic radiation force to assess the resulting oscillatory displacement denoting the underlying tissue stiffness. The objective of this study was to evaluate the feasibility of HMI in pancreatic tumor detection and high-intensity focused ultrasound (HIFU) treatment monitoring. The HMI system consisted of a focused ultrasound transducer, which generated sinusoidal radiation force to induce oscillatory tissue motion at 50 Hz, and a diagnostic ultrasound transducer, which detected the axial tissue displacements based on acquired radio-frequency signals using a 1-D cross-correlation algorithm. For pancreatic tumor detection, HMI images were generated for pancreatic tumors in transgenic mice and normal pancreases in wild-type mice. The obtained HMI images showed a high contrast between normal and malignant pancreases with an average peak-to-peak HMI displacement ratio of 3.2. Histological analysis showed that no tissue damage was associated with HMI when it was used for the sole purpose of elasticity imaging. For pancreatic tumor ablation monitoring, the focused ultrasound transducer was operated at a higher acoustic power and longer pulse length than that used in tumor detection to simultaneously induce HIFU thermal ablation and oscillatory tissue displacements, allowing HMI monitoring without interrupting tumor ablation. HMI monitoring of HIFU ablation found significant decreases in the peak-to-peak HMI displacements before and after HIFU ablation with a reduction rate ranging from 15.8% to 57.0%. The formation of thermal lesions after HIFU exposure was confirmed by histological analysis. This study demonstrated the feasibility of HMI in abdominal tumor detection and HIFU ablation monitoring.

High-intensity focused ultrasound monitoring using harmonic motion imaging for focused ultrasound (HMIFU) under boiling or slow denaturation conditions

Harmonic motion imaging for focused ultrasound (HMIFU) is a recently developed high-intensity focused ultrasound (HIFU) treatment monitoring method that utilizes an amplitude-modulated therapeutic ultrasound beam to induce an oscillatory radiation force at the HIFU focus and estimates the focal tissue displacement to monitor the HIFU thermal treatment. In this study, the performance of HMIFU under acoustic, thermal, and mechanical effects was investigated. The performance of HMIFU was assessed in ex vivo canine liver specimens (n = 13) under slow denaturation or boiling regimes. A passive cavitation detector (PCD) was used to assess the acoustic cavitation activity, and a bare-wire thermocouple was used to monitor the focal temperature change. During lesioning with slow denaturation, high quality displacements (correlation coefficient above 0.97) were observed under minimum cavitation noise, indicating the tissue initial-softening-then- stiffening property change. During HIFU with boiling, HMIFU monitored a consistent change in lesion-to-background displacement contrast (0.46 ± 0.37) despite the presence of strong cavitation noise due to boiling during lesion formation. Therefore, HMIFU effectively monitored softening-then-stiffening during lesioning under slow denaturation, and detected lesioning under boiling with a distinct change in displacement contrast under boiling in the presence of cavitation. In conclusion, HMIFU was shown under both boiling and slow denaturation regimes to be effective in HIFU monitoring and lesioning identification without being significantly affected by cavitation noise.

Sparse matrix beamforming and image reconstruction for 2-D HIFU monitoring using harmonic motion imaging for focused ultrasound (HMIFU) with in vitro validation

Harmonic motion imaging for focused ultrasound (HMIFU) utilizes an amplitude-modulated HIFU beam to induce a localized focal oscillatory motion simultaneously estimated. The objective of this study is to develop and show the feasibility of a novel fast beamforming algorithm for image reconstruction using GPU-based sparse-matrix operation with real-time feedback. In this study, the algorithm was implemented onto a fully integrated, clinically relevant HMIFU system. A single divergent transmit beam was used while fast beamforming was implemented using a GPU-based delay-and-sum method and a sparse-matrix operation. Axial HMI displacements were then estimated from the RF signals using a 1-D normalized cross-correlation method and streamed to a graphic user interface with frame rates up to 15 Hz, a 100-fold increase compared to conventional CPU-based processing. The real-time feedback rate does not require interrupting the HIFU treatment. Results in phantom experiments showed reproducible HMI images and monitoring of 22 in vitro HIFU treatments using the new 2-D system demonstrated reproducible displacement imaging, and monitoring of 22 in vitro HIFU treatments using the new 2-D system showed a consistent average focal displacement decrease of 46.7 ±14.6% during lesion formation. Complementary focal temperature monitoring also indicated an average rate of displacement increase and decrease with focal temperature at 0.84±1.15%/(°)C, and 2.03±0.93%/(°)C , respectively. These results reinforce the HMIFU capability of estimating and monitoring stiffness related changes in real time. Current ongoing studies include clinical translation of the presented system for monitoring of HIFU treatment for breast and pancreatic tumor applications.

Multi-parametric monitoring and assessment of high-intensity focused ultrasound (HIFU) boiling by harmonic motion imaging for focused ultrasound (HMIFU): an ex vivo feasibility study

Harmonic motion imaging for focused ultrasound (HMIFU) is a recently developed high-intensity focused ultrasound (HIFU) treatment monitoring method with feasibilities demonstrated in vitro and in vivo. Here, a multi-parametric study is performed to investigate both elastic and acoustics-independent viscoelastic tissue changes using the Harmonic Motion Imaging (HMI) displacement, axial compressive strain and change in relative phase shift during high energy HIFU treatment with tissue boiling. Forty three (n = 43) thermal lesions were formed in ex vivo canine liver specimens (n = 28). Two-dimensional (2D) transverse HMI displacement maps were also obtained before and after lesion formation. The same method was repeated in 10 s, 20 s and 30 s HIFU durations at three different acoustic powers of 8, 10, and 11 W, which were selected and verified as treatment parameters capable of inducing boiling using both thermocouple and passive cavitation detection (PCD) measurements. Although a steady decrease in the displacement, compressive strain, and relative change in the focal phase shift (Δϕ) were obtained in numerous cases, indicating an overall increase in relative stiffness, the study outcomes also showed that during boiling, a reverse lesion-to-background displacement contrast was detected, indicating potential change in tissue absorption, geometrical change and/or, mechanical gelatification or pulverization. Following treatment, corresponding 2D HMI displacement images of the thermal lesions also mapped consistent discrepancy in the lesion-to-background displacement contrast. Despite the expectedly chaotic changes in acoustic properties with boiling, the relative change in phase shift showed a consistent decrease, indicating its robustness to monitor biomechanical properties independent of the acoustic property changes throughout the HIFU treatment. In addition, the 2D HMI displacement images confirmed and indicated the increase in the thermal lesion size with treatment duration, which was validated against pathology. In conclusion, multi-parametric HMIFU was shown capable of monitoring and mapping tissue viscoelastic response changes during and after HIFU boiling, some of which were independent of the acoustic parameter changes.

Flexible integration of high-imaging-resolution and high-power arrays for ultrasound-induced thermal strain imaging (US-TSI)

Ultrasound-induced thermal strain imaging (USTSI) for carotid artery plaque detection requires both high imaging resolution (<100 μm) and sufficient US-induced heating to elevate the tissue temperature (~1°C to 3°C within 1 to 3 cardiac cycles) to produce a noticeable change in sound speed in the targeted tissues. Because the optimization of both imaging and heating in a monolithic array design is particularly expensive and inflexible, a new integrated approach is presented which utilizes independent ultrasound arrays to meet the requirements for this particular application. This work demonstrates a new approach in dual-array construction. A 3-D printed manifold was built to support both a high-resolution 20 MHz commercial imaging array and 6 custom heating elements operating in the 3.5 to 4 MHz range. For the application of US-TSI in carotid plaque characterization, the tissue target site is 20 to 30 mm deep, with a typical target volume of 2 mm (elevation) × 8 mm (azimuthal) × 5 mm (depth). The custom heating array performance was fully characterized for two design variants (flat and spherical apertures), and can easily deliver 30 W of total acoustic power to produce intensities greater than 15 W/cm(2) in the tissue target region.

Ultrasonic evaluation of microwave-induced thermal lesions based on wavelet analysis of mean scatterer spacing

The microwave ablation has become an important manner for tumor treatment. In this paper, we proposed a new method for evaluation of microwave-induced thermal lesions using the wavelet analysis of the mean scatterer spacing (MSS). First, the ultrasonic radiofrequency (RF) data of normal and coagulated porcine liver tissues was collected through the temperature-controlled water bath heating experiments. The convex array ultrasound probe with a center frequency of 3.5 MHz was used. Second, the wavelet analysis was used to compute the MSS of normal and coagulated porcine liver tissues, respectively. Finally, the microwave-induced thermal lesions were detected based on the differences in the MSS between normal tissues and thermal lesions. Eighteen cases of microwave ablation experiments and 20 cases of water bath heating experiments were conducted on fresh porcine liver samples. The MSS of normal porcine liver tissues was 1.15±0.12 mm, and the MSS of coagulated porcine liver tissues was 0.93±0.07 mm. Six cases of thermal lesions were compared between the MSS-detected area and the caliper-measured area, and the MSS-detected area had an error of (13.55±5.29) %. The experimental results indicated that the proposed method could be used in preliminary detection and evaluation of microwave-induced thermal lesions.