Real time shear waves elastography monitoring of thermal ablation: in vivo evaluation in pig livers

Acoustic Cavitation Emissions Predict Near-complete/complete Histotripsy Treatment in Soft Tissues

OBJECTIVE

Histotripsy is a non-invasive acoustic ablation technique that leverages cavitation to impart mechanical damage to a viscoelastic medium, such as tissue. Although histotripsy bubbles and lesions can be imaged with a variety of modalities, reliable methods to predict tissue disruption across different tissue-types remain to be determined.

APPROACH

Several ex-vivo bovine tissues were ablated by intrinsic threshold histotripsy over a range of pulse-per-location acoustic doses. Acoustic Cavitation Emission (ACE) signals were captured following every other therapeutic pulse using transmit-receive capable histotripsy arrays. Final bubble lifespan, lifespan-slope, and percent-reduction were calculated and correlated against histologic necrosis score (0-5: 0=0% necrosis, 5=>95% necrosis) and residual structure score (0-4: 0=none present, 4=intact) to evaluate the ability of features from ACE-signals to predict histotripsy-induced damage. Further, optimal ACE-feature thresholds were determined for binary evaluation of whether a necrosis score equal or greater than 4 had been reached.

RESULTS

Measured lifespans increased and lifespan-slopes decreased with pulses per location (ppl) and eventually plateaued in all tissue types, in similar trends to those previously observed in tissue phantoms. Necrosis score increased and residual structure decreased with increasing acoustic dose. Bubble lifespan-slope and percent-reduction correlated well with necrosis score. Thresholds able to predict the necrosis score of 4 or greater in brain, liver, and kidney were calculated with high sensitivity/specificity (>80%). The necrosis score of 4 and 5 is expected to correspond to near-complete/complete ablation by histological evaluation.

CONCLUSION

Features measured from ACE-signals, particularly the lifespan-slope and percent reduction, were used to predict near-complete/complete ablation of large-volume histotripsy treatments in ex vivo bovine liver, kidney, and brain tissues with good accuracy. Tissue heterogeneities were observed to impact the histotripsy damage and corresponding ACE-signals, and thus the predication accuracy.

Image-Based Monitoring of Thermal Ablation

Thermal therapy is a commonly used local treatment technique in clinical practice. Monitoring the treatment process is essential for ensuring its success. In this review, we analyze recent image-based methods for thermal therapy monitoring, focusing particularly on their feasibility for synchronous or immediate postoperative monitoring. This includes thermography and other techniques that track the physical changes in tissue during thermal ablation. Potential directions and challenges for further clinical applications are also summarized.

Therapeutic ultrasound transducer technology and monitoring techniques: a review with clinical examples

The exponential growth of therapeutic ultrasound applications demonstrates the power of the technology to leverage the combinations of transducer technology and treatment monitoring techniques to effectively control the preferred bioeffect to elicit the desired clinical effect.Objective: This review provides an overview of the most commonly used bioeffects in therapeutic ultrasound and describes existing transducer technologies and monitoring techniques to ensure treatment safety and efficacy.Methods and materials: Literature reviews were conducted to identify key choices that essential in terms of transducer design, treatment parameters and procedure monitoring for therapeutic ultrasound applications. Effective combinations of these options are illustrated through descriptions of several clinical indications, including uterine fibroids, prostate disease, liver cancer, and brain cancer, that have been successful in leveraging therapeutic ultrasound to provide effective patient treatments.Results: Despite technological constraints, there are multiple ways to achieve a desired bioeffect with therapeutic ultrasound in a target tissue. Visualizations of the interplay of monitoring modality, bioeffect, and applied acoustic parameters are presented that demonstrate the interconnectedness of the field of therapeutic ultrasound. While the clinical indications explored in this review are at different points in the clinical evaluation path, based on the ever expanding research being conducted in preclinical realms, it is clear that additional clinical applications of therapeutic ultrasound that utilize a myriad of bioeffects will continue to grow and improve in the coming years.Conclusions: Therapeutic ultrasound will continue to improve in the next decades as the combination of transducer technology and treatment monitoring techniques will continue to evolve and be translated in clinical settings, leading to more personalized and efficient therapeutic ultrasound mediated therapies.

Alternating direction method of multipliers for displacement estimation in ultrasound strain elastography

BACKGROUND

Ultrasound strain imaging, which delineates mechanical properties to detect tissue abnormalities, involves estimating the time delay between two radio-frequency (RF) frames collected before and after tissue deformation. The existing regularized optimization-based time-delay estimation (TDE) techniques suffer from at least one of the following drawbacks: (1) The regularizer is not aligned with the tissue deformation physics due to taking only the first-order displacement derivative into account; (2) TheL 2 $L2$ -norm of the displacement derivatives, which oversmooths the estimated time-delay, is utilized as the regularizer; (3) The modulus function defined mathematically should be approximated by a smooth function to facilitate the optimization ofL 1 $L1$ -norm.

PURPOSE

Our purpose is to develop a novel TDE technique that resolves the aforementioned shortcomings of the existing algorithms.

METHODS

Herein, we propose employing the alternating direction method of multipliers (ADMM) for optimizing a novel cost function consisting ofL 2 $L2$ -norm data fidelity term andL 1 $L1$ -norm first- and second-order spatial continuity terms. ADMM empowers the proposed algorithm to use different techniques for optimizing different parts of the cost function and obtain high-contrast strain images with smooth backgrounds and sharp boundaries. We name our technique ADMM for totaL variaTion RegUlarIzation in ultrasound STrain imaging (ALTRUIST). ALTRUIST's efficacy is quantified using absolute error (AE), Structural SIMilarity (SSIM), signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR), and strain ratio (SR) with respect to GLUE, OVERWIND, andL 1 $L1$ -SOUL, three recently published energy-based techniques, and UMEN-Net, a state-of-the-art deep learning-based algorithm. Analysis of variance (ANOVA)-led multiple comparison tests and paired t $t$ -tests at5 % $5\%$ overall significance level were conducted to assess the statistical significance of our findings. The Bonferroni correction was taken into account in all statistical tests. Two simulated layer phantoms, three simulated resolution phantoms, one hard-inclusion simulated phantom, one multi-inclusion simulated phantom, one experimental breast phantom, and three in vivo liver cancer datasets have been used for validation experiments. We have published the ALTRUIST code at http://code.sonography.ai.

RESULTS

ALTRUIST substantially outperforms the four state-of-the-art benchmarks in all validation experiments, both qualitatively and quantitatively. ALTRUIST yields up to573 % ∗ ${573\%}^{*}$ ,41 % ∗ ${41\%}^{*}$ , and51 % ∗ ${51\%}^{*}$ SNR improvements and443 % ∗ ${443\%}^{*}$ ,53 % ∗ ${53\%}^{*}$ , and15 % ∗ ${15\%}^{*}$ CNR improvements overL 1 $L1$ -SOUL, its closest competitor, for simulated, phantom, and in vivo liver cancer datasets, respectively, where the asterisk (*) indicates statistical significance. In addition, ANOVA-led multiple comparison tests and paired t $t$ -tests indicate that ALTRUIST generally achieves statistically significant improvements over GLUE, UMEN-Net, OVERWIND, andL 1 $L1$ -SOUL in terms of AE, SSIM map, SNR, and CNR.

CONCLUSIONS

A novel ultrasonic displacement tracking algorithm named ALTRUIST has been developed. The principal novelty of ALTRUIST is incorporating ADMM for optimizing anL 1 $L1$ -norm regularization-based cost function. ALTRUIST exhibits promising performance in simulation, phantom, and in vivo experiments.

Reverse time migration and genetic algorithms Combined for Reconstruction in Transluminal Shear Wave Elastography: An In Silico Case Study

A new reconstruction approach that combines Reverse Time Migration (RTM) and Genetic Algorithms (GAs) is proposed for solving the inverse problem associated with transluminal shear wave elastography. The transurethral identification of the first thermal lesion generated by transrectal High Intensity Focused Ultrasound (HIFU) for the treatment of prostate cancer, was used to preliminarily test in silico the combined reconstruction method. The RTM method was optimised by comparing reconstruction images from several cross-correlation techniques, including a new proposed one, and different device configurations in terms of the number and arrangement of emitters and receivers of the conceptual transurethral probe. The best results were obtained for the new proposed cross-correlation method and a device configuration with 3 emitters and 32 receivers. The RTM reconstructions did not completely contour the shape of the HIFU lesion, however, as planned for the combined approach, the areas in the RTM images with high level of correlation were used to narrow down the search space in the GA-based technique. The GA-based technique was set to find the location of the HIFU lesion and the increment in stiffness and viscosity due to thermal damage. Overall, the combined approach achieves lower level of error in the reconstructed values, and in a shorter computational time, compared to the GA-based technique alone. The lowest errors were accomplished for the location of HIFU lesion, followed by the contrast ratio of stiffness between thermally treated tissue and non-treated normal tissue. The homologous ratio of viscosity obtained higher level of error. Further investigation considering diverse scenarios to be reconstructed and with experimental data is required to fully evaluate the feasibility of the combined approach.

Magnetic Resonance Acoustic Radiation Force Imaging (MR-ARFI) for the monitoring of High Intensity Focused Ultrasound (HIFU) ablation in anisotropic tissue

OBJECTIVE

We introduce a non-invasive MR-Acoustic Radiation Force Imaging (ARFI)-based elastography method that provides both the local shear modulus and temperature maps for the monitoring of High Intensity Focused Ultrasound (HIFU) therapy.

MATERIALS AND METHODS

To take tissue anisotropy into account, the local shear modulus μ is determined in selected radial directions around the focal spot by fitting the phase profiles to a linear viscoelastic model, including tissue-specific mechanical relaxation time τ. MR-ARFI was evaluated on a calibrated phantom, then applied to the monitoring of HIFU in a gel phantom, ex vivo and in vivo porcine muscle tissue, in parallel with MR-thermometry.

RESULTS

As expected, the shear modulus polar maps reflected the isotropy of phantoms and the anisotropy of muscle. In the HIFU monitoring experiments, both the shear modulus polar map and the thermometry map were updated with every pair of MR-ARFI phase images acquired with opposite MR-ARFI-encoding. The shear modulus was found to decrease (phantom and ex vivo) or increase (in vivo) during heating, before remaining steady during the cooling phase. The mechanical relaxation time, estimated pre- and post-HIFU, was found to vary in muscle tissue.

DISCUSSION

MR-ARFI allowed for monitoring of viscoelasticity changes around the HIFU focal spot even in anisotropic muscle tissue.

Simultaneous proton resonance frequency T1 - MR shear wave elastography for MR-guided focused ultrasound multiparametric treatment monitoring

PURPOSE

To develop an efficient MRI pulse sequence to simultaneously measure multiple parameters that have been shown to correlate with tissue nonviability following thermal therapies.

METHODS

A 3D segmented EPI pulse sequence was used to simultaneously measure proton resonance frequency shift (PRFS) MR thermometry (MRT), T1 relaxation time, and shear wave velocity induced by focused ultrasound (FUS) push pulses. Experiments were performed in tissue mimicking gelatin phantoms and ex vivo bovine liver. Using a carefully designed FUS triggering scheme, a heating duty cycle of approximately 65% was achieved by interleaving FUS ablation pulses with FUS push pulses to induce shear waves in the tissue.

RESULTS

In phantom studies, temperature increases measured with PRFS MRT and increases in T1 correlated with decreased shear wave velocity, consistent with material softening with increasing temperature. During ablation in ex vivo liver, temperature increase measured with PRFS MRT initially correlated with increasing T1 and decreasing shear wave velocity, and after tissue coagulation with decreasing T1 and increasing shear wave velocity. This is consistent with a previously described hysteresis in T1 versus PRFS curves and increased tissue stiffness with tissue coagulation.

CONCLUSION

An efficient approach for simultaneous and dynamic measurements of PRSF, T1 , and shear wave velocity during treatment is presented. This approach holds promise for providing co-registered dynamic measures of multiple parameters, which correlates to tissue nonviability during and following thermal therapies, such as FUS.

Monitoring MR-guided high intensity focused ultrasound therapy using transient supersonic shear wave MR-elastography

Objective.The aim of the paper is to propose an all-in-one method based on magnetic resonance-supersonic shear wave imaging (MR-SSI) and proton resonance frequency shift (PRFS) to monitor high intensity focused ultrasound (HIFU) thermal ablations.Approach.Mechanical properties have been shown to be related to tissue damage induced by thermal ablations. Monitoring elasticity in addition to temperature changes may help in ensuring the efficacy and the accuracy of HIFU therapies. For this purpose, an MR-SSI method has been developed where the ultrasonic transducer is used for both mechanical wave generation and thermal ablation. Transient quasi-planar shear waves are generated using the acoustic radiation force, and their propagation is monitored in motion-sensitized phase MR images. Using a single-shot gradient-echo echo-planar-imaging sequence, MR images can be acquired at a sufficiently high temporal resolution to provide an update of PRFS thermometry and MR-SSI elastography maps in real time.Main results.The proposed method was first validated on a calibrated elasticity phantom, in which both the possibility to detect inclusions with different stiffness and repeatability were demonstrated. The standard deviation between the 8 performed measurements was 2% on the background of the phantom and 11%, at most, on the inclusions. A second experiment consisted in performing a HIFU heating in a gelatin phantom. The temperature increase was estimated to be 9 °C and the shear modulus was found to decrease from 2.9 to 1.8 kPa, reflecting the gel softening around the HIFU focus, whereas it remained steady in non-heated areas.Significance.The proposed MR-SSI technique allows monitoring HIFU ablations using thermometry and elastography simultaneously, without the need for an additional external mechanical exciter such as those used in MR elastography.

Incorporating Gradient Similarity for Robust Time Delay Estimation in Ultrasound Elastography

Energy-based ultrasound elastography techniques minimize a regularized cost function consisting of data and continuity terms to obtain local displacement estimates based on the local time-delay estimation (TDE) between radio frequency (RF) frames. The data term associated with the existing techniques takes only the amplitude similarity into account and hence is not sufficiently robust to the outlier samples present in the RF frames under consideration. This drawback creates noticeable artifacts in the strain image. To resolve this issue, we propose to formulate the data function as a linear combination of the amplitude and gradient similarity constraints. We estimate the adaptive weight concerning each similarity term following an iterative scheme. Finally, we optimize the nonlinear cost function in an efficient manner to convert the problem to a sparse system of linear equations which are solved for millions of variables. We call our technique rGLUE: robust data term in GLobal Ultrasound Elastography. rGLUE has been validated using simulation, phantom, in vivo liver, and breast datasets. In all our experiments, rGLUE substantially outperforms the recent elastography methods both visually and quantitatively. For simulated, phantom, and in vivo datasets, respectively, rGLUE achieves 107%, 18%, and 23% improvements of signal-to-noise ratio (SNR) and 61%, 19%, and 25% improvements of contrast-to-noise ratio (CNR) over global ultrasound elastography (GLUE), a recently published elastography algorithm.

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.

Viscoelasticity Imaging of Biological Tissues and Single Cells Using Shear Wave Propagation

Changes in biomechanical properties of biological soft tissues are often associated with physiological dysfunctions. Since biological soft tissues are hydrated, viscoelasticity is likely suitable to represent its solid-like behavior using elasticity and fluid-like behavior using viscosity. Shear wave elastography is a non-invasive imaging technology invented for clinical applications that has shown promise to characterize various tissue viscoelasticity. It is based on measuring and analyzing velocities and attenuations of propagated shear waves. In this review, principles and technical developments of shear wave elastography for viscoelasticity characterization from organ to cellular levels are presented, and different imaging modalities used to track shear wave propagation are described. At a macroscopic scale, techniques for inducing shear waves using an external mechanical vibration, an acoustic radiation pressure or a Lorentz force are reviewed along with imaging approaches proposed to track shear wave propagation, namely ultrasound, magnetic resonance, optical, and photoacoustic means. Then, approaches for theoretical modeling and tracking of shear waves are detailed. Following it, some examples of applications to characterize the viscoelasticity of various organs are given. At a microscopic scale, a novel cellular shear wave elastography method using an external vibration and optical microscopy is illustrated. Finally, current limitations and future directions in shear wave elastography are presented.

Wave Propagation in a Fractional Viscoelastic Tissue Model: Application to Transluminal Procedures

In this article, a wave propagation model is presented as the first step in the development of a new type of transluminal procedure for performing elastography. Elastography is a medical imaging modality for mapping the elastic properties of soft tissue. The wave propagation model is based on a Kelvin Voigt Fractional Derivative (KVFD) viscoelastic wave equation, and is numerically solved using a Finite Difference Time Domain (FDTD) method. Fractional rheological models, such as the KVFD, are particularly well suited to model the viscoelastic response of soft tissue in elastography. The transluminal procedure is based on the transmission and detection of shear waves through the luminal wall. Shear waves travelling through the tissue are perturbed after encountering areas of altered elasticity. These perturbations carry information of medical interest that can be extracted by solving the inverse problem. Scattering from prostate tumours is used as an example application to test the model. In silico results demonstrate that shear waves are satisfactorily transmitted through the luminal wall and that echoes, coming from reflected energy at the edges of an area of altered elasticity, which are feasibly detectable by using the transluminal approach. The model here presented provides a useful tool to establish the feasibility of transluminal procedures based on wave propagation and its interaction with the mechanical properties of the tissue outside the lumen.

Magnetic resonance shear wave elastography using transient acoustic radiation force excitations and sinusoidal displacement encoding

A magnetic resonance (MR) shear wave elastography technique that uses transient acoustic radiation force impulses from a focused ultrasound (FUS) transducer and a sinusoidal-shaped MR displacement encoding strategy is presented. Using this encoding strategy, an analytic expression for calculating the shear wave speed in a heterogeneous medium was derived. Green's function-based simulations were used to evaluate the feasibility of calculating shear wave speed maps using the analytic expression. Accuracy of simulation technique was confirmed experimentally in a homogeneous gelatin phantom. The elastography measurement was compared to harmonic MR elastography in a homogeneous phantom experiment and the measured shear wave speed values differed by less than 14%. This new transient elastography approach was able to map the position and shape of inclusions sized from 8.5 to 14 mm in an inclusion phantom experiment. These preliminary results demonstrate the feasibility of using a straightforward analytic expression to generate shear wave speed maps from MR images where sinusoidal-shaped motion encoding gradients are used to encode the displacement-time history of a transiently propagating wave-packet. This new measurement technique may be particularly well suited for performing elastography before, during, and after MR-guided FUS therapies since the same device used for therapy is also used as an excitation source for elastography.

Study on the relationship between reduced scattering coefficient and Young's modulus of tumors in microwave ablation

BACKGROUND AND OBJECTIVE

In the clinical treatment of tumors using microwave ablation (MWA), although temperature can be used as an important reference index for evaluating the curative effect of ablation, it cannot fully reflect the biological activity status of tumor tissue during thermal ablation. Finding multi-parameter comprehensive evaluation factors to achieve real-time evaluation of therapeutic effects has become the key for precise ablation. More and more scholars use the reduced scattering coefficient (μs') and Young's modulus (E) to evaluate the treatment outcomes of MWA. However, the intrinsic relationship between these parameters is unclear. This paper aims to investigate the specific relationship between μs' and E during MWA.

MATERIAL AND METHODS

The MWA experiment was conducted on porcine liver in vitro, the two-parameter simultaneous acquisition system was designed to obtain the reduced scattering coefficient and Young's modulus of the liver tissue during MWA. The relationship between reduced scattering coefficient and Young's modulus was investigated.

RESULTS

It is found that the trend of change of μs' is very similar to E in the process of MWA, i.e. first increasing and then reaching a steady state, and in some experiments there are synchronous changes. Based on this, the quantitative relationship between E-μs'  is established, enabling the quantitative estimation of Young's modulus of liver tissue based on reduced scattering coefficient. The maximum absolute error is 29.37 kPa and the minimum absolute error is 0.88 kPa.

CONCLUSION

This study contributes to the further establishment of a multi-parameter MWA effectiveness evaluation model. It is also valuable for clinically evaluating the ablation outcomes of tumor in real time.

Photoacoustic Spectral Sensing Technique for Diagnosis of Biological Tissue Coagulation: In-Vitro Study

Thermal coagulation of abnormal tissues has evolved as a therapeutic technique for different diseases including cancer. Tissue heating beyond 55 °C causes coagulation that leads to cell death. Noninvasive diagnosis of thermally coagulated tissues is pragmatic for performing efficient therapy as well as reducing damage of surrounding healthy tissues. We propose a noninvasive, elasticity-based photoacoustic spectral sensing technique for differentiating normal and coagulated tissues. Photoacoustic diagnosis is performed for quantitative differentiation of normal and coagulated excised chicken liver and muscle tissues in vitro by characterizing a dominant frequency of photoacoustic frequency spectrum. Pronounced distinction in the spectral parameter (i.e., dominant frequency) was observed due to change in tissue elastic property. We confirmed nearly two-fold increase in dominant frequencies for the coagulated muscle and liver tissues as compared to the normal ones. A density increase caused by tissue coagulation is clearly reflected in the dominant frequency composition. Experimental results were consistent over five different sample sets, delineating the potential of proposed technique to diagnose biological tissue coagulation and thus monitor thermal coagulation therapy in clinical applications.

Simultaneous fat-referenced proton resonance frequency shift thermometry and MR elastography for the monitoring of thermal ablations

PURPOSE

Simultaneous fat-referenced proton resonance frequency shift (FRPRFS) thermometry combined with MR elastography (MRE) is proposed, to continuously monitor thermal ablations for all types of soft tissues, including fat-containing tissues. Fat-referenced proton resonance frequency shift thermometry makes it possible to measure temperature even in the water fraction of fat-containing tissues while enabling local field-drift correction. Magnetic resonance elastography allows measuring the mechanical properties of tissues that are related to tissue structural damage.

METHODS

A gradient-echo MR sequence framework was proposed that combines the need for multiple TE acquisitions for the water-fat separation of FRPRFS, and the need for multiple MRE phase offsets for elastogram reconstructions. Feasibility was first assessed in a fat-containing gelatin phantom undergoing moderate heating by a hot water circulation system. Subsequently, high intensity focused ultrasound heating was conducted in porcine muscle tissue ex vivo (N = 4; 2 samples, 2 locations/sample).

RESULTS

Both FRPRFS temperature maps and elastograms were updated every 4.1 seconds. In the gelatin phantom, FRPRFS was in good agreement with optical fiber thermometry (average difference 1.2 ± 1°C). In ex vivo high-intensity focused ultrasound experiments on muscle tissue, the shear modulus was found to decrease significantly by 34.3% ± 7.7% (experiment 1, sample 1), 17.9% ± 10.0% (experiment 2, sample 1), 55.1% ± 8.7% (experiment 3, sample 2), and 34.7% ± 8.4% (experiment 4, sample 2) as a result of temperature increase (ΔT = 22.5°C ± 4.2°C, 14.0°C ± 2.8°C, 14.7°C ± 3.7°C, and 14.5°C ± 3.0°C, respectively).

CONCLUSION

This study demonstrated the feasibility of monitoring thermal ablations with FRPRFS thermometry together with MRE, even in fat-containing tissues. The acquisition time is similar to non-FRPRFS thermometry combined with MRE.

Preliminary analysis of ultrasound elastography imaging-based thermometry on non-perfused ex vivo swine liver

AIMS

Real-time monitoring of tissue temperature during percutaneous tumor ablation improves treatment efficacy, leading clinicians in adjustment of treatment settings. This study aims at assessing feasibility of ultrasound thermometry during laser ablation of biological tissue using a specific ultrasound imaging techniques based on elastography acoustic radiation force impulse (ARFI).

METHODS

ARFI uses high-intensity focused ultrasound pulses to generate 'radiation force' in tissue; this provokes tissue displacements trackable using correlation-based ultrasound methods: the sensitivity of shear waves velocity is able to detect temperature changes. Experiments were carried out using a Nd:YAG laser (power: 5 W) in three non-perfused ex vivo pig livers. In each organ, a thermocouple was placed close to the applicator tip (distance range 1.5-2.5 cm) used to record a reference temperature. Positioning of laser applicator and thermocouple was eco-guided. The organ was scanned by an echography system equipped with ARFI; propagation velocity was measured in a region of interest of 1 × 0.5 cm located close to thermocouple, to investigate influence of tissue temperature on shear waves velocity.

RESULTS

Shear wave velocity has a very low sensitivity to temperature up to 55-60 °C, and in all cases, velocity is < 5 m s^-1; for temperature > 55-60 °C, velocity shows a steep increment. The system measures a value "over limit", meaning a velocity > 5 m s^-1.

CONCLUSIONS

Ultrasound thermometry during laser ablation of biological tissue based on elastography shows an abrupt output change at temperatures > 55-60 °C. This issue can have a relevant clinical impact, considering tumor necrosis when temperature crosses 55 °C to define the boundary of damaged volume.

Performance of Shear Wave Elastography in Delineating the Radiofrequency Ablation Boundary: An in Vivo experiment

This study was aimed at exploring the cutoff value of Young's modulus of ablated tissue and the optimal scale at which shear wave elastography (SWE) can delineate the ablation boundary. The livers of 30 rabbits were radiofrequency (RF) ablated, and ultrasonic imaging, including SWE and contrast-enhanced ultrasound (CEUS), was performed. The ablation boundary in the SWE image was located using CEUS, and the SWE parameters of the boundary were measured to calculate the cutoff value of Young's modulus. The cutoff value of the ablated tissue was 48-50 kPa 2 h to 28 d post-ablation. The regions of increased stiffness in SWE images at a scale of 0-50 kPa overlapped well with the non-enhanced regions of CEUS images in 88% of specimens. Therefore, elasticity values differed significantly between ablated and non-ablated tissues, and the cutoff value for Young's modulus differentiated these tissues. SWE delineated the ablation boundary well at the optimal SWE scale with respect to the cutoff value.

Monitoring Microwave Ablation Using Ultrasound Echo Decorrelation Imaging: An ex vivo Study

In this study, a microwave-induced ablation zone (thermal lesion) monitoring method based on ultrasound echo decorrelation imaging was proposed. A total of 15 cases of ex vivo porcine liver microwave ablation (MWA) experiments were carried out. Ultrasound radiofrequency (RF) signals at different times during MWA were acquired using a commercial clinical ultrasound scanner with a 7.5-MHz linear-array transducer. Instantaneous and cumulative echo decorrelation images of two adjacent frames of RF data were calculated. Polynomial approximation images were obtained on the basis of the thresholded cumulative echo decorrelation images. Experimental results showed that the instantaneous echo decorrelation images outperformed conventional B-mode images in monitoring microwave-induced thermal lesions. Using gross pathology measurements as the reference standard, the estimation of thermal lesions using the polynomial approximation images yielded an average accuracy of 88.60%. We concluded that instantaneous ultrasound echo decorrelation imaging is capable of monitoring the change of thermal lesions during MWA, and cumulative ultrasound echo decorrelation imaging and polynomial approximation imaging are feasible for quantitatively depicting thermal lesions.

Evaluation of liver lesions by use of shear wave elastography and computed tomography perfusion imaging after radiofrequency ablation in clinically normal dogs

OBJECTIVE To evaluate acute changes of the liver by use of shear wave elastography (SWE) and CT perfusion after radiofrequency ablation (RFA). ANIMALS 7 healthy Beagles. PROCEDURES RFA was performed on the liver (day 0). Stiffness of the ablation lesion, transitional zone, and normal parenchyma were evaluated by use of SWE, and blood flow, blood volume, and arterial liver perfusion of those regions were evaluated by use of CT perfusion on days 0 and 4. All RFA lesions were histologically examined on day 4. RESULTS Examination of the SWE color-coded map distinctly revealed stiffness of the liver tissue, which increased from the normal parenchyma to the transitional zone and then to the ablation zone. For CT perfusion, blood flow, blood volume, and arterial liver perfusion decreased from the transitional zone to the normal parenchyma and then to the ablation zone. Tissue stiffness and CT perfusion variables did not differ significantly between days 0 and 4. Histologic examination revealed central diffuse necrosis and peripheral hyperemia with infiltration of lymphoid cells and macrophages. CONCLUSIONS AND CLINICAL RELEVANCE Coagulation necrosis induced a loss of blood perfusion and caused tissue hardening (stiffness) in the ablation zone. Hyperemic and inflammatory changes of the transitional zone resulted in increased blood perfusion. Acute changes in stiffness and perfusion of liver tissue after RFA could be determined by use of SWE and CT perfusion. These results can be used to predict the clinical efficacy of RFA and to support further studies, including those involving hepatic neoplasia.

Stiffness distribution in the ablated zone after radiofrequency ablation for liver: An ex-vivo study with a tissue elastometer

OBJECTIVE

To investigate the stiffness distribution in the ablated zone after radiofrequency ablation (RFA), we used a device called tissue elastometer based on gross liver samples.

MATERIALS

AND METHODS: Twelve freshly excised porcine livers were subject to RFA under a same setup to form elliptic ablated samples. Each sample was cut open for gross examination, and then the surface of the section plane was sliced into one piece for Young's modulus test using the tissue elastometer. Five test points along the long- and short-axis on each piece were selected to evaluate stiffness distribution respectively. Among them, four points distributed equidistantly from center to boundary in the ablated zone and one was in the unablated zone.

RESULTS

In the ablated zone, we found the Young's moduli were significantly different among the four test points both in long- (F = 99.04, p <0.001) and short-axis (F = 79.47, p <0.001) directions. The Young's modulus showed a downtrend in each direction, and was linearly related to the distance from the center to the test point (for long axis, R2 = 0.968; for short axis, R2 = 0.984, both p <0.001). A more significant downtrend was observed in short-axis direction. The Young's moduli gained from the inner edge of ablated zone were comparable and significantly higher than those from the outer edge for both directions. The maximum value of 24.71kPa for Young's modulus was the appropriate threshold to ensure the tissues were necrotic completely.

CONCLUSION

The stiffness inside the ablated zone represented a radial distribution with downtrend, following a linear law. The stiffness at the inner edge of ablated zone is stable and significantly higher than that at the outer edge. The maximum value of 24.71 kPa close to the inner edge of Wz may be used as the standard of complete ablation.

Ultrasound single-phase CBE imaging for monitoring radiofrequency ablation

Radiofrequency (RF) ablation (RFA) is the most commonly used minimally invasive procedure for thermal ablation of liver tumors. Ultrasound not only provides real-time feedback of the electrode location for RFA guidance but also enables visualization of the tissue temperature. Changes in backscattered energy (CBE) have been widely applied to ultrasound temperature imaging for assessing thermal ablation. Pilot studies have revealed that significant shadowing features appear in CBE imaging and are caused by the electrode and RFA-induced gas bubbles. To resolve this problem, the current study proposed ultrasound single-phase CBE imaging based on positive CBE values. An in vitro model with tissue samples derived from the porcine tenderloin was used to validate the proposed method. During RFA with various electrode lengths, ultrasound scans of tissue samples were obtained using a clinical ultrasound scanner equipped with a convex array transducer of 3 MHz. Raw image data comprising 256 scan lines of backscattered RF signals were acquired for B-mode, conventional CBE, and single-phase CBE imaging by using the proposed algorithmic scheme. The ablation sizes estimated using CBE imaging and gross examinations were compared to calculate the correlation coefficient. The experimental results indicated that single-phase CBE imaging largely suppressed artificial CBE information in the shadowed region. Moreover, compared with conventional CBE imaging, single-phase CBE imaging provided a more accurate estimation of ablation sizes (the correlation coefficient was higher than 0.8).

Acoustic Radiation Force Based Ultrasound Elasticity Imaging for Biomedical Applications

Pathological changes in biological tissue are related to the changes in mechanical properties of biological tissue. Conventional medical screening tools such as ultrasound, magnetic resonance imaging or computed tomography have failed to produce the elastic properties of biological tissues directly. Ultrasound elasticity imaging (UEI) has been proposed as a promising imaging tool to map the elastic parameters of soft tissues for the clinical diagnosis of various diseases include prostate, liver, breast, and thyroid gland. Existing UEI-based approaches can be classified into three groups: internal physiologic excitation, external excitation, and acoustic radiation force (ARF) excitation methods. Among these methods, ARF has become one of the most popular techniques for the clinical diagnosis and treatment of disease. This paper provides comprehensive information on the recently developed ARF-based UEI techniques and instruments for biomedical applications. The mechanical properties of soft tissue, ARF and displacement estimation methods, working principle and implementation instruments for each ARF-based UEI method are discussed.

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.

FEM-based elasticity reconstruction using ultrasound for imaging tissue ablation

PURPOSE

Success of ablation treatment depends on the accurate placement of the target ablation focus and the complete destruction of the pathological tissue. Thus, monitoring the formation, location, and size of the ablated lesion is essential. As ablated tissue gets stiffer, an option for ablation monitoring is ultrasound elastography, for imaging the tissue mechanical properties. Reconstruction of elasticity distribution can be achieved by solving an inverse problem from observed displacements, based on a deformable tissue model, commonly discretized by the finite element method (FEM). However, available reconstruction techniques are prone to noise and may achieve suboptimal accuracy.

METHODS

We propose a novel inverse problem formulation and elasticity reconstruction method, in which both the elasticity parameters and the model displacements are estimated as independent parameters of an unconstrained optimization problem. Total variation regularization of spatial elasticity distribution is introduced in this formulation, providing robustness to noise.

RESULTS

Our approach was compared to state of the art direct and iterative harmonic elastography techniques. We employed numerical simulation studies using various noise and inclusion contrasts, given multiple excitation frequencies. Compared to alternatives, our method leads to a decrease in RMSE of up to 50% and an increase in CNR of up to 11 dB in numerical simulations. The methods were also compared on an ex vivo bovine liver sample that was locally subjected to ablation, for which improved lesion delineation was obtained with our proposed method. Our method takes [Formula: see text] for [Formula: see text] reconstruction grid.

CONCLUSIONS

We present a novel FEM problem formulation that improves reconstruction accuracy and inclusion delineation compared to currently available techniques.

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.

Nanoparticle-based hyperthermia, a local treatment modulating the tumor extracellular matrix

The structural complexity and physical properties of the tumor microenvironment negatively affect the penetration and efficiency of conventional anticancer drugs. While previously underestimated, the tumor microenvironment now becomes a potential target for cancer treatment. This microenvironment can be modulated either systemically by pharmacological means, or locally, through physical effects mediated by certain nanoparticles. Some of them, such as magnetic, plasmonic or carbon-based nanoparticles, can generate heat on demand in a spatially and temporally controlled manner. In addition, the nanoparticles can be either activated by light or magnetic stimuli. The impact of the resulting local heating can be observed on the ultrastructural level, as it strongly affects the organization of collagen fibers, and on the macroscopic level, since the thermal damages alter the mechanical properties of the tumor. Nanoparticle-based hyperthermia thus improves the effect of conventional anticancer drugs, as it allows their better penetration through the altered extracellular matrix. Here we suggest the use of nanoparticle-generated hyperthermia, obtained after magnetic or light activation, as an adjuvant treatment to prime the tumor microenvironment and improve the efficacy of chemotherapy.

Tumor Stiffening, a Key Determinant of Tumor Progression, is Reversed by Nanomaterial-Induced Photothermal Therapy

Tumor stiffening, stemming from aberrant production and organization of extracellular matrix (ECM), has been considered a predictive marker of tumor malignancy, non-invasively assessed by ultrasound shear wave elastography (SWE). Being more than a passive marker, tumor stiffening restricts the delivery of diagnostic and therapeutic agents to the tumor and per se could modulate cellular mechano-signaling, tissue inflammation and tumor progression. Current strategies to modify the tumor extracellular matrix are based on ECM-targeting chemical agents but also showed deleterious systemic effects.

On-demand excitable nanomaterials have shown their ability to perturb the tumor microenvironment in a spatiotemporal-controlled manner and synergistically with chemotherapy. Here, we investigated the evolution of tumor stiffness as well as tumor integrity and progression, under the effect of mild hyperthermia and thermal ablation generated by light-exposed multi-walled carbon nanotubes (MWCNTs) in an epidermoid carcinoma mouse xenograft. SWE was used for real-time mapping of the tumor stiffness, both during the two near infrared irradiation sessions and over the days after the treatment. We observed a transient and reversible stiffening of the tumor tissue during laser irradiation, which was lowered at the second session of mild hyperthermia or photoablation. In contrast, over the days following photothermal treatment, the treated tumors exhibited a significant softening together with volume reduction, whereas non-treated growing tumors showed an increase of tumor rigidity. The organization of the collagen matrix and the distribution of CNTs revealed a spatio-temporal correlation between the presence of nanoheaters and the damages on collagen and cells. This study highlights nanohyperthermia as a promising adjuvant strategy to reverse tumor stiffening and normalize the mechanical tumor environment.

Comparison of three dimensional strain volume reconstructions using SOUPR and wobbler based acquisitions: A phantom study

PURPOSE

Ultrasound strain imaging is a relatively low cost and portable modality for monitoring percutaneous thermal ablation of liver neoplasms. However, a 3D strain volume reconstruction from existing 2D strain images is necessary to fully delineate the thermal dose distribution. Tissue mimicking (TM) phantom experiments were performed to validate a novel volume reconstruction algorithm referred to as sheaf of ultrasound planes reconstruction (SOUPR), based on a series of 2D rotational imaging planes.

METHODS

Reconstruction using SOUPR was formulated as an optimization problem with constraints on data consistency with 2D strain images and data smoothness of the volume data. Reconstructed ablation inclusion dimensions, volume, and elastographic signal to noise ratio (SNRe) and contrast to noise ratio (CNRe) were compared with conventional 3D ultrasound strain imaging based on interpolating a series of quasiparallel 2D strain images with a wobbler transducer.

RESULTS

Volume estimates of the phantom inclusion were in a similar range for both acquisition approaches. SNRe and CNRe obtained with SOUPR were significantly higher on the order of 250% and 166%, respectively. The mean error of the inclusion dimension reconstructed with a wobbler transducer was on the order of 10.4%, 3.5%, and 19.0% along the X, Y, and Z axes, respectively, while the error with SOUPR was on the order of 2.6%, 2.8%, and 9.6%. A qualitative comparison of SOUPR and wobbler reconstruction was also performed using a thermally ablated region created in ex vivo bovine liver tissue.

CONCLUSIONS

The authors have demonstrated using experimental evaluations with a TM phantom that the reconstruction results obtained with SOUPR were superior when compared with a conventional wobbler transducer in terms of inclusion shape preservation and detectability.

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.

Use of Shear Wave Ultrasound Vibrometry for Detection of Simulated Esophageal Malignancy in Ex Vivo Porcine Esophagi

Esophageal cancer is a malignant neoplasm with poor outcomes. Determination of local disease progression is a major determining factor in treatment modality, radiation dose, radiation field and subsequent surgical therapy. Discrimination of true tumor extent is difficult given the similarity of soft tissues of the malignancy compared to non-malignant tissues using current imaging modalities. A possible method to discriminate between these tissues may be to exploit mechanical properties to diagnostic advantage, as malignant tissues tend to be stiffer relative to normal adjacent tissue. Shear waves propagate faster in stiffer tissues relative to softer tissues. This may be measured by using ultrasound based shear wave vibrometry. In this method, acoustic radiation force is used to create a shear wave in the tissue of interest and ultrafast ultrasound imaging is used to track the propagating wave to measure the wave velocity and estimate the shear moduli. In this study we created simulated malignant lesions (1.5 cm length) using radiofrequency ablation in ex vivo esophageal samples with varied progression (partial thickness n = 4, and full thickness n = 5) and used normal regions of the same esophageal specimen as controls. Shear wave vibrometry was used to measure shear wave group velocity and shear wave phase velocity in the ex vivo specimens. These values were used to estimate shear moduli using an elastic shear wave model and elastic and viscoelastic Lamb wave models. Our results show that the group and phase velocities increase due to both full and mucosal ablation, and that discrimination may be provided by higher order analysis using viscoelastic Lamb wave fitting. This technique may have application for determination of extent of early esophageal malignancy and warrants further investigation using in vivo approaches to determine performance compared to current imaging modalities.

Monitoring Radiofrequency Ablation Using Ultrasound Envelope Statistics and Shear Wave Elastography in the Periablation Period: An In Vitro Feasibility Study

Radiofrequency ablation (RFA) is a minimally invasive method for treating tumors. Shear wave elastography (SWE) has been widely applied in evaluating tissue stiffness and final ablation size after RFA. However, the usefulness of periablation SWE imaging in assessing RFA remains unclear. Therefore, this study investigated the correlation between periablation SWE imaging and final ablation size. An in vitro porcine liver model was used for experimental validation (n = 36). During RFA with a power of 50 W, SWE images were collected using a clinical ultrasound system. To evaluate the effects of tissue temperature and gas bubbles during RFA, changes in the ablation temperature were recorded, and image echo patterns were measured using B-mode and ultrasound statistical parametric images. After RFA, the gross pathology of each tissue sample was compared with the region of change in the corresponding periablation SWE image. The experimental results showed that the tissue temperature at the ablation site varied between 70°C and 100°C. Hyperechoic regions and changes were observed in the echo amplitude distribution induced by gas bubbles. Under this condition, the confounding effects (including the temperature increase, tissue stiffness increase, and presence of gas bubbles) resulted in artifacts in the periablation SWE images, and the corresponding region correlated with the estimated final ablation size obtained from the gross pathology (r = 0.8). The findings confirm the feasibility of using periablation SWE imaging in assessing RFA.

Relative Elastic Modulus Imaging Using Sector Ultrasound Data for Abdominal Applications: An Evaluation of Strategies and Feasibility

We reconstruct the elastic modulus distribution for one tissue mimicking (TM) phantom and two in vivo biopsy-confirmed liver tumors using curvilinear ultrasound echo data. Spatial distribution of the relative elastic modulus values is determined by solving an inverse problem within a region of interest (ROI). This inverse problem solution requires knowledge of the ultrasonically measured displacement field in a uniform rectilinear grid to ensure that the resolution on the resultant relative elastic modulus elastogram will be uniform over the entire ROI. Taking advantage of a new speckle tracking algorithm, two different displacement tracking strategies are investigated: 1) sector-shaped ultrasound data were converted to ultrasound data on a rectilinear grid prior to speckle tracking and 2) axial and lateral displacements directly obtained from sector-shaped data were converted to vertical and horizontal displacements on a rectilinear grid after speckle tracking. Compared with strain elastography (SE), TM phantom results show that relative elastic modulus imaging (REMI) using Strategy 2 provided higher contrast-to-noise ratios (>300% and 25% increases compared with SE and REMI using Strategy 1, respectively). Furthermore, in phantoms, REMI using Strategy 2 more accurately (a 1.3% difference to shear wave elastography measurements) estimated the elastic contrast ratio between the target and the background, compared with both SE (20%-25%) and REMI using Strategy 1 (4.1%). It was also observed that relative modulus elastograms were more consistent with anatomical structures visualized on corresponding B-mode images for the two in vivo liver cases. Overall, we conclude that applying REMI is feasible for abdominal organs such as the liver. Strategy 2 offered improved and consistent results for the data investigated.

Diagnostic Accuracy of 2D-Shear Wave Elastography for Liver Fibrosis Severity: A Meta-Analysis

Purpose

To evaluate the accuracy of shear wave elastography (SWE) in the quantitative diagnosis of liver fibrosis severity.

Methods

The published literatures were systematically retrieved from PubMed, Embase, Web of science and Scopus up to May 13^th, 2016. Included studies reported the pooled sensitivity, specificity, positive and negative predictive values, as well as the diagnostic odds ratio of SWE in populations with liver fibrosis. A bivariate mixed-effects regression model was used, which was estimated by the I² statistics. The quality of articles was evaluated by quality assessment of diagnostic accuracy studies (QUADAS).

Results

Thirteen articles including 2303 patients were qualified for the study. The pooled sensitivity and specificity of SWE for the diagnosis of liver fibrosis are as follows: ≥F1 0.76 (p<0.001, 95% CI, 0.71–0.81, I² = 75.33%), 0.92 (p<0.001, 95% CI, 0.80–0.97, I² = 79.36%); ≥F2 0.84 (p = 0.35, 95% CI, 0.81–0.86, I² = 9.55%), 0.83 (p<0.001, 95% CI, 0.77–0.88, I² = 86.56%); ≥F3 0.89 (p = 0.56, 95% CI, 0.86–0.92, I² = 0%), 0.86 (p<0.001, 95% CI, 0.82–0.90, I² = 75.73%); F4 0.89 (p = 0.24, 95% CI, 0.84–0.92, I² = 20.56%), 0.88 (p<0.001, 95% CI, 0.84–0.92, I² = 82.75%), respectively. Sensitivity analysis showed no significant changes if any one of the studies was excluded. Publication bias was not detected in this meta-analysis.

Conclusions

Our study suggests that SWE is a helpful method to appraise liver fibrosis severity. Future studies that validate these findings would be appropriate.

Comparison of real-time contrast-enhanced ultrasonography and standard ultrasonography in liver cancer microwave ablation

Primary liver cancer has a high incidence and high mortality rates, and currently the only viable option is surgery, although there are a number of difficulties related to this method. The aim of the present study was to investigate the potential advantages of the real-time contrast-enhanced ultrasonography (CEUS) for microwave ablation of primary liver cancer. One hundred patients with primary liver cancer were included in the study. The patients were divided into the ordinary ultrasonography and the CEUS groups. For the ordinary ultrasonography group, the ordinary ultrasonography-guided microwave ablation method was used, while microwave ablation under the guidance of CEUS was conducted for the CEUS group. The size of lesions and clearness of the tumor boundary prior to surgery in the two groups were compared. Additionally, postoperative complications and the survival rate were monitored. Lesion boundary areas measured by CEUS were significantly larger than those measured with ordinary ultrasonography. The incidence rate of postoperative pain, fever, intra-abdominal hemorrhage and infection and other complications in the ordinary ultrasonography group were significantly higher than that in the CEUS group. The tumor recurrence rate in the CEUS group was significantly lower than that in the ordinary ultrasonography group. Seventy-two percent of patients in the CEUS group showed no progress, compared to 48% of in the ordinary ultrasonography group. The progress-free survival rate in the CEUS group after 6 months was significantly higher than that in the ordinary ultrasonography group. Disease-free survival time in the CEUS group was considerably longer than the control group. In conclusion, the guidance of real-time CEUS on the primary liver cancer microwave ablation treatment can achieve good intra-operative results. It offers a real-time guidance effect, improves survival time and reduces the incidence of complications.

Tumor characterization and treatment monitoring of postsurgical human breast specimens using harmonic motion imaging (HMI)

BACKGROUND

High-intensity focused ultrasound (HIFU) is a noninvasive technique used in the treatment of early-stage breast cancer and benign tumors. To facilitate its translation to the clinic, there is a need for a simple, cost-effective device that can reliably monitor HIFU treatment. We have developed harmonic motion imaging (HMI), which can be used seamlessly in conjunction with HIFU for tumor ablation monitoring, namely harmonic motion imaging for focused ultrasound (HMIFU). The overall objective of this study was to develop an all ultrasound-based system for real-time imaging and ablation monitoring in the human breast in vivo.

METHODS

HMI was performed in 36 specimens (19 normal, 15 invasive ductal carcinomas, and 2 fibroadenomas) immediately after surgical removal. The specimens were securely embedded in a tissue-mimicking agar gel matrix and submerged in degassed phosphate-buffered saline to mimic in vivo environment. The HMI setup consisted of a HIFU transducer confocally aligned with an imaging transducer to induce an oscillatory radiation force and estimate the resulting displacement.

RESULTS

3D HMI displacement maps were reconstructed to represent the relative tissue stiffness in 3D. The average peak-to-peak displacement was found to be significantly different (p = 0.003) between normal breast tissue and invasive ductal carcinoma. There were also significant differences before and after HMIFU ablation in both the normal (53.84 % decrease) and invasive ductal carcinoma (44.69 % decrease) specimens.

CONCLUSIONS

HMI can be used to map and differentiate relative stiffness in postsurgical normal and pathological breast tissues. HMIFU can also successfully monitor thermal ablations in normal and pathological human breast specimens. This HMI technique may lead to a new clinical tool for breast tumor imaging and HIFU treatment monitoring.

Utilizing confocal laser endomicroscopy for evaluating the adequacy of laparoscopic liver ablation

Background

Laparoscopic liver ablation therapy can be used for the treatment of primary and secondary liver malignancy. The increased incidence of cancer recurrence associated with this approach, has been attributed to the inability of monitoring the extent of ablated liver tissue.

Methods

The feasibility of assessing liver ablation with probe‐based confocal laser endomicroscopy (CLE) was studied in a porcine model of laparoscopic microwave liver ablation. Following the intravenous injection of the fluorophores fluorescein and indocyanine green, CLE images were recorded at 488 nm and 660 nm wavelength and compared to liver histology. Statistical analysis was performed to assess if fluorescence intensity change can predict the presence of ablated liver tissue.

Results

CLE imaging of fluorescein at 488 nm provided good visualization of the hepatic microvasculature; whereas, CLE imaging of indocyanine green at 660 nm enabled detailed visualization of hepatic sinusoid architecture and interlobular septations. Fluorescence intensity as measured in relative fluorescence units was found to be 75–100% lower in ablated compared to healthy liver regions. General linear mixed modeling and ROC analysis found the decrease in fluorescence to be statistically significant.

Conclusion

Laparoscopic, dual wavelength CLE imaging using two different fluorophores enables clinically useful visualization of multiple liver tissue compartments, in greater detail than is possible at a single wavelength. CLE imaging may provide valuable intraoperative information on the extent of laparoscopic liver ablation. Lasers Surg. Med. 48:299–310, 2016. © 2015 The Authors. Lasers in Surgery and Medicine Published by Wiley Periodicals, Inc.

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.

Cardiac shear-wave elastography using a transesophageal transducer: application to the mapping of thermal lesions in ultrasound transesophageal cardiac ablation

Heart rhythm disorders, such as atrial fibrillation or ventricular tachycardia can be treated by catheter-based thermal ablation. However, clinically available systems based on radio-frequency or cryothermal ablation suffer from limited energy penetration and the lack of lesion's extent monitoring. An ultrasound-guided transesophageal device has recently successfully been used to perform High-Intensity Focused Ultrasound (HIFU) ablation in targeted regions of the heart in vivo. In this study we investigate the feasibility of a dual therapy and imaging approach on the same transesophageal device. We demonstrate in vivo that quantitative cardiac shear-wave elastography (SWE) can be performed with the device and we show on ex vivo samples that transesophageal SWE can map the extent of the HIFU lesions. First, SWE was validated with the transesophageal endoscope in one sheep in vivo. The stiffness of normal atrial and ventricular tissues has been assessed during the cardiac cycle (n = 11) and mapped (n = 7). Second, HIFU ablation has been performed with the therapy-imaging transesophageal device in ex vivo chicken breast samples (n  =  3), then atrial (left, n = 2) and ventricular (left n = 1, right n = 1) porcine heart tissues. SWE provided stiffness maps of the tissues before and after ablation. Areas of the lesions were obtained by tissue color change with gross pathology and compared to SWE. During the cardiac cycle stiffness varied from 0.5   ±   0.1 kPa to 6.0   ±   0.3 kPa in the atrium and from 1.3   ±   0.3 kPa to 13.5   ±   9.1 kPa in the ventricles. The thermal lesions were visible on all SWE maps performed after ablation. Shear modulus of the ablated zones increased to 16.3   ±   5.5 kPa (versus 4.4   ±   1.6 kPa before ablation) in the chicken breast, to 30.3   ±   10.3 kPa (versus 12.2   ±   4.3 kPa) in the atria and to 73.8   ±   13.9 kPa (versus 21.2   ±   3.3 kPa) in the ventricles. On gross pathology, the size of the lesions ranged from 0.1 to 1.5 cm(2) in the imaging plane area. Elasticity-estimated depths and widths of the lesions differed respectively with a median of 0.2 mm (first quartile Q1:  -0.8 mm; third quartile Q3: 2.6 mm) for a mean squared error (MSE) of 5.1 mm(2) and a median of 0.2 mm (Q1:  -2.7 mm; Q3: 2.7 mm) for a MSE of 11.1 mm(2) from gross pathology. We have demonstrated the feasibility of the HIFU thermal ablation monitoring using a dual therapy and imaging transesophageal device. The combination of HIFU, ultrasound imaging and SWE on the same transesophageal system could lead to a new clinical device for a safer and controlled treatment of a wide variety of cardiac arrhythmias.

Quantitative evaluation of atrial radio frequency ablation using intracardiac shear-wave elastography

PURPOSE

Radio frequency catheter ablation (RFCA) is a well-established clinical procedure for the treatment of atrial fibrillation (AF) but suffers from a low single-procedure success rate. Recurrence of AF is most likely attributable to discontinuous or nontransmural ablation lesions. Yet, despite this urgent clinical need, there is no clinically available imaging modality that can reliably map the lesion transmural extent in real time. In this study, the authors demonstrated the feasibility of shear-wave elastography (SWE) to map quantitatively the stiffness of RFCA-induced thermal lesions in cardiac tissues in vitro and in vivo using an intracardiac transducer array.

METHODS

SWE was first validated in ex vivo porcine ventricular samples (N = 5). Both B-mode imaging and SWE were performed on normal cardiac tissue before and after RFCA. Areas of the lesions were determined by tissue color change with gross pathology and compared against the SWE stiffness maps. SWE was then performed in vivo in three sheep (N = 3). First, the stiffness of normal atrial tissues was assessed quantitatively as well as its variation during the cardiac cycle. SWE was then performed in atrial tissue after RFCA.

RESULTS

A large increase in stiffness was observed in ablated ex vivo regions (average shear modulus across samples in normal tissue: 22 ± 5 kPa, average shear-wave speed (ct): 4.5 ± 0.4 m s(-1) and in determined ablated zones: 99 ± 17 kPa, average ct: 9.0 ± 0.5 m s(-1) for a mean shear modulus increase ratio of 4.5 ± 0.9). In vivo, a threefold increase of the shear modulus was measured in the ablated regions, and the lesion extension was clearly visible on the stiffness maps.

CONCLUSIONS

By its quantitative and real-time capabilities, Intracardiac SWE is a promising intraoperative imaging technique for the evaluation of thermal ablation during RFCA.

The role of intraoperative ultrasonography in the diagnosis and management of focal hepatic lesions

Intraoperative ultrasonography (IOUS) has been widely utilized in hepatic surgery both as a diagnostic technique and in the course of treatment. Since IOUS involves direct-contact imaging of the target organ, it can provide high spatial resolution without interference from the surrounding structures. Therefore, IOUS may improve the detection, characterization, localization, and local staging of hepatic tumors. IOUS is also a real-time imaging modality capable of providing interactive information and valuable guidance in a range of procedures. Recently, contrast-enhanced IOUS, IOUS elastography, and IOUS-guided hepatic surgery have attracted increasing interest and are expected to lead to the broader implementation of IOUS. Herein, we review the various applications of IOUS in the diagnosis and management of focal hepatic lesions.

Thermometry and ablation monitoring with ultrasound

In this review we present the current status of ultrasound thermometry and ablation monitoring, with emphasis on the diverse approaches published in the literature and with an eye on which methods are closest to clinical reality. It is hoped that this review will serve as a guide to the expansion of sonographic methods for treatment monitoring and thermometry since the last brief review in 2007.

Monitoring radiofrequency ablation using real-time ultrasound Nakagami imaging combined with frequency and temporal compounding techniques

Gas bubbles induced during the radiofrequency ablation (RFA) of tissues can affect the detection of ablation zones (necrosis zone or thermal lesion) during ultrasound elastography. To resolve this problem, our previous study proposed ultrasound Nakagami imaging for detecting thermal-induced bubble formation to evaluate ablation zones. To prepare for future applications, this study (i) created a novel algorithmic scheme based on the frequency and temporal compounding of Nakagami imaging for enhanced ablation zone visualization, (ii) integrated the proposed algorithm into a clinical scanner to develop a real-time Nakagami imaging system for monitoring RFA, and (iii) investigated the applicability of Nakagami imaging to various types of tissues. The performance of the real-time Nakagami imaging system in visualizing RFA-induced ablation zones was validated by measuring porcine liver (n = 18) and muscle tissues (n = 6). The experimental results showed that the proposed algorithm can operate on a standard clinical ultrasound scanner to monitor RFA in real time. The Nakagami imaging system effectively monitors RFA-induced ablation zones in liver tissues. However, because tissue properties differ, the system cannot visualize ablation zones in muscle fibers. In the future, real-time Nakagami imaging should be focused on the RFA of the liver and is suggested as an alternative monitoring tool when advanced elastography is unavailable or substantial bubbles exist in the ablation zone.

A survey of ultrasound elastography approaches to percutaneous ablation monitoring

Percutaneous thermal ablation has been widely used as a minimally invasive treatment for tumors. Treatment monitoring is essential for preventing complications while ensuring treatment efficacy. Mechanical testing measurements on tissue reveal that tissue stiffness increases with temperature and ablation duration. Different types of imaging methods can be used to monitor ablation procedures, including temperature or thermal strain imaging, strain imaging, modulus imaging, and shear modulus imaging. Ultrasound elastography demonstrates the potential to become the primary imaging modality for monitoring percutaneous ablation. This review briefly presented the state-of-the-art ultrasound elastography approaches for monitoring radiofrequency ablation and microwave ablation. These techniques were divided into four groups: quasi-static elastography, acoustic radiation force elastography, sonoelastography, and applicator motion elastography. Their advantages and limitations were compared and discussed. Future developments were proposed with respect to heat-induced bubbles, tissue inhomogeneities, respiratory motion, three-dimensional monitoring, multi-parametric monitoring, real-time monitoring, experimental data center for percutaneous ablation, and microwave ablation monitoring.