|
|
||||||||
1 Department of Medical Physics, The University of Wisconsin-Madison, 1530
Medical Sciences Center, 1300 University Ave., Madison, WI 53706.
2 Department of Electrical Engineering, The University of Wisconsin-Madison,
Madison, WI 53706.
3 Department of Radiology, The University of Wisconsin-Madison, Madison, WI
53706.
Received November 4, 2002;
accepted after revision March 18, 2003.
Supported in part by grants from The University of Wisconsin-Madison and by
grant R01 CA39224 from the National Institutes of Health.
Abstract
|
|
|---|
MATERIALS AND METHODS. Radiofrequency ablation was performed in vitro on liver tissue specimens encased in gelatin phantoms. The imaging plane for elastography was perpendicular to the axis of the radiofrequency electrode so that the ablated region was around the center of the plane. To obtain three-dimensional visualization of thermal lesions, we reconstructed the lesions from multiple elastograms by linearly translating the elastographic scanning plane. Pathology photographs were obtained in the same image plane used for elastography by slicing through the gelatin and tissue phantom using external markers. We used digitized gross pathology photographs obtained at a specified slice thickness to compute the areas and volumes of the lesions. These measurements were then compared to the measurements obtained from the elastograms.
RESULTS. In a sample of 40 thermal lesions, we obtained a correlation between in vitro elastographic and pathologic measurements of r = 0.9371 (p < 0.00001) for area estimates and r = 0.979 (p < 0.00001) for volume estimates.
CONCLUSION. We found excellent correlation between the measurements of the dimensions, areas, and volumes of thermal lesions that were based on elastographic images and the measurements that were based on digitized pathologic images. When compared with digitized pathologic measurements, elastographic measurements showed a tendency to slighly underestimate both the areas and volumes of lesions. Nevertheless, elastography is a reliable technique for delineating thermal lesions resulting from radiofrequency ablation.
|
|
|---|
Elastography has recently been shown to be useful in monitoring the impact of ablative therapies on tumors [17, 18, 19, 20]. Techniques such as radiofrequency ablation produce protein denaturation because of the elevated tissue temperatures; the elevated temperatures, in turn, induce elevation of the elastic modulus of soft tissue. Temperatures greater than 42°C (107.6°F) are considered potentially lethal, depending on the duration of application [21], and temperatures greater than 60°C (140°F) are associated with uniform tissue necrosis [22]. On elastography, radiofrequency ablation lesions exhibit extremely high contrast compared with the normal, untreated background liver tissue [17, 18]. Early work in other laboratories [17, 18, 19, 20] and our own preliminary results [17, 18] suggest that sonographic elastography may be substantially superior to conventional sonography for monitoring radiofrequency ablation lesions.
Conventional transabdominal sonography has been used extensively for guidance and monitoring of radiofrequency ablation of hepatic tumors [23, 24]. However, the zone of necrosis produced by radiofrequency ablation is not easily visualized on sonography because of the intrinsically low contrast between normal and ablated liver tissue and because of artifacts resulting from gas bubble formation and from the generator [23, 25]. Microbubble sonography contrast agents have also been investigated as aids in revealing the zone of necrosis generated after radiofrequency ablation [26]. Because these agents are used to differentiate between perfused and nonperfused tissue, the use of contrast sonography may help in the detection of residual or viable nonablated tumor regions.
MRIs of ablated regions reveal altered signals on both T1- and T2-weighted images [23, 25, 27, 28], with the treated areas devoid of gadolinium enhancement.
Radiologic and pathologic correlations in both experimental and clinical studies have shown that CT and MRI can aid radiologists in predicting the region of coagulation within 2-3 mm [29]. Cha et al. [30] compared the usefulness of CT and sonography as methods of monitoring radiofrequency ablation therapy. These researchers found that on CT, radiofrequency ablation lesions consistently showed lower attenuation than the surrounding liver. On sonography, echogenicity of the lesions varied, with approximately 59% of the lesions being hypoechoic; 25%, hyperechoic; and about 16%, isoechoic [30]. Both unenhanced and contrast-enhanced CT offered better conspicuity and edge detection of radiofrequency ablation lesions and fewer artifacts than did sonography. Unenhanced CT had the closest correlation to pathologic size of the lesions (r = 0.74), followed by contrast-enhanced CT (r = 0.72), and sonography (r = 0.56).
The delineation of the boundary of an ablated region is key to making accurate judgments concerning treated versus untreated regions. Our goal in the work described here was to understand more thoroughly how the boundaries of thermal lesions viewed on elastograms correspond to actual ablated tissue boundaries observed histologically. We compared the widths, heights, areas, and volumes of thermal lesions depicted on elastograms with the same parameters of the lesions measured at gross pathology. We performed our analysis under in vitro conditions to evaluate the feasibility of using elastography for monitoring the thermal lesions generated with radiofrequency ablation.
|
|
|---|
|
An electrosurgical device (model 1500, RITA Medical Systems, Mountain View, CA) was used for radiofrequency ablation procedures. The radiofrequency ablation electrode in this device consists of a 15-gauge shaft through which multiple sharp tines, each with a diameter of 0.053 cm (0.021 inches; 25 gauge), can be deployed. Fully extended, the tines are in an "umbrella" configuration at 45° intervals. The last 1 cm of the electrode tip and each tine constitute the electrically active surface. We inserted the electrode into the liver lobe and deployed the tines, taking care to keep them within the liver parenchyma. Radiofrequency ablation of the targeted tissue was conducted for 10 min by heating the tissue to set temperatures of 70°, 80°, 90°, or 100°C (158°, 176°, 194°, 212°F), using a 150-W power level.
Freshly excised porcine livers obtained from an adjoining laboratory were used for the analysis. In our initial set of experiments, the thermal lesions were created after normal sections of liver tissue were encased in gelatin. However, ablation of liver tissue inside the gelatin blocks caused some melting of the gelatin, which may have changed the stiffness properties. To avoid this occurrence, we performed radiofrequency ablation on whole liver samples before encasing them in gelatin blocks. Ablation was performed in three to five sections of each liver sample by inserting the radiofrequency ablation electrode about 2-4 cm deep and then switching on the power. After the lesion had cooled, the lesion and the surrounding tissue were encased in a gelatin block for subsequent imaging.
Gelatin phantoms were prepared using 200-bloom calfskin gelatin at a concentration of 15.4 g/L of distilled water. The gelatin powder was mixed with distilled water and cooked in a double-heated water bath to a temperature of 80°C (194°F). After the gelatin solution clarified, glass beads (mean diameter, 8 µm) were mixed with warm distilled water and added to the molten gelatin at a concentration of 1 g of beads per liter of gelatin to provide a source of scattered echo signals during sonography. The molten gelatin was then poured into plexiglass molds. First, a 2.5-cm thick layer of molten gelatin with a temperature of 29°C (82.4°F) was poured. After this layer congealed, we placed a section of liver with a radiofrequency ablation lesion on the gelatin layer. The mold was then filled with molten gel (temperature, 29°C or 82.4°F), and the entire phantom was allowed to congeal for 1 hr. The molds containing the gelatin blocks were refrigerated at 10°C (50°F) for 2 hr. Samples were then removed from the molds for elastography.
We used a real-time scanner (SSD 2000, Aloka, Tokyo, Japan) for elastographic imaging of the thermal lesion, with a 40-mm-wide 5.0-MHz linear array transducer and a 70% bandwidth. In addition, we obtained elastograms using a 3.5-MHz sector transducer to assess the use of lower frequency probes. The raw sonographic radiofrequency echo signals from complete B-mode frames were digitized using a 12-bit data acquisition board (Gage, Montreal, QC, Canada) at a sampling rate of 50 MHz. The apparatus includes a stepper motor-controlled motion-compression system and a compression plate that also supports the transducer. Both the transducer and compression plate were placed in direct contact with the gelatin block. A personal computer was used to control the operation of the system. Elastograms were generated using compressive increments of 0.5%. Signal processing was performed off-line. We computed tissue displacements with time-domain cross-correlation analysis using a 3-mm window length setting, with a 50% overlap between data segments. Elastograms were then obtained using a gradient operator on the estimates of tissue displacement.
|
|
|---|
|
|
|
|
|
|
|
A three-dimensional (3D) visualization of a thermal lesion was generated from parallel, multislice elastograms acquired at 1-mm increments (Fig. 3D). Segmentation of the two-dimensional (2D) elastograms to obtain the lesion surface was performed by setting a single threshold (a strain value 90% below the maximum strain surrounding the thermal lesion) to separate the background strains from the lower strains in the thermal lesion. The 3D representation allows computation of the lesion volume.
Thermal lesion dimensions, areas, and volumes were estimated for 40 thermal lesions to analyze correspondence between elastographic and visual pathologic depictions. Ten sets of in vitro thermal lesions (a set consists of four lesions generated at temperatures of 70°, 80°, 90°, and 100°C or (158°, 176°, 194°, and 212°F) were studied using elastography and compared with corresponding digitized photographs of gross pathology specimens that were sliced along the central sonographic and elastographic scanning planes. Fiducial markers were inscribed on these planes of the gelatin phantom to align both data sets.
Scatterplots of the lesion dimensions (width and height) are presented in Figures 4 and 5, respectively. As these scatterplots show, we found that the lesion widths obtained using elastography corresponded closely to the width of the thermal lesion on the digitized photographs of pathology specimens (r = 0.8693). The heights of the lesions on elastograms had a slightly better correspondence to measurements obtained from the pathologic specimens (r = 0.8742), implying a better fit between the data sets. This result is probably explained by the fact that lesion heights were measured in the axial direction, which provides the best sonographic and elastographic resolution. We obtained the areas of the thermal lesion from the width and height measurements so that we could compute the area of the bounding ellipse (Fig. 6), where the r value obtained was 0.9126. However, obtaining areas in this manner does not take into account the errors that may be caused by the shift in the angle that the thermal lesion makes along the principal axes.
|
|
|
Better estimates of the lesion areas were obtained by drawing contours of the lesion shape (including the angle produced by the thermal lesion with the principal axes) and then computing the area of the bounding ellipse. The areas of lesions obtained from the pathology specimens were measured by drawing contours (yellow outline in Figure 3C) that contained only the white zone (region of complete necrosis) of the thermal lesion from the optically scanned images of the gross pathology specimens. The area calculations for elastograms were obtained by drawing a contour (green outline in Figure 3B) that contained only the apparent thermal lesion (i.e., inside the black or stiffer region). Correlation between estimates of areas obtained from elastography and those obtained from pathology specimens improved (r = 0.9371) using this method (Fig. 7). The areas obtained using elastography slightly underestimated the lesion areas. Compared with the best correlation (r = 0.74) reported by Cha et al. [30] between pathologic and unenhanced CT measurements of lesion sizes, our results were significantly better.
|
Volume estimates for the lesions in gross pathology specimens (Fig. 8) were calculated from a sum of the product of the areas of the thermal lesions multiplied by the thickness of the thermal lesions section, sliced in 1-mm increments. The elastographic lesion volumes were calculated as the sum of the product of the areas obtained from 2D elastograms multiplied by the scanning plane separation of 1 mm. The r value between elastographic and pathologic volumes was 0.979. As we found in the estimates of thermal lesion areas, we discovered an apparent underestimation of the volumes using elastography. The underestimation of the lesion areas and volume using elastography may be an accurate reflection of the underlying stiffness changes, but we plan to investigate this result more extensively.
|
|
|
|---|
In our statistical results, the r value between elastographic lesion areas and pathologic lesion areas was approximately 0.93. The comparison with elastographic and pathologic volumes provided higher correlations (r = 0.97). These findings indicate that elastography is potentially a viable technique for delineating thermal lesions.
As a result of this finding, clinicians may be able to assess the success of the procedure or the need to perform additional ablation immediately, as opposed to waiting for a few days to obtain CT scans. In contrast to conventional sonography, 3D surface and volume rendering of thermal lesions generated from elastographic images is made possible by the marked contrast in stiffness between the lesion and the surrounding normal tissue. Imaging of the thermal lesion in 3D also allows computation of the lesion volume, size, and position within the treated region, parameters that indicate the success or failure of the procedure. Although the diameters of thermal lesions imaged in our study were 1-3 cm, the elastograms obtained using the sector probe and the 3.5-MHz center frequency showed that larger and deeper lesions can also be monitored with elastography.
For our analysis, we used in vitro samples in which the controlled quasistatic compressions required for elastography can be easily applied. The feasibility of using in vivo elastography to image radiofrequency-ablated lesions in the liver has already been shown [17, 18]. The use of the radiofrequency ablation electrode as the compressor-displacement device reduces the lateral slippage or nonaxial motion that may occur with externally applied compressions or imaging during the respiratory cycle [17, 18]. This technique provides controlled and reproducible compressions of the liver for in vivo elastography. A complete statistical analysis in a study using in vivo thermal lesions would confirm the usefulness of elastography in monitoring radiofrequency ablation procedures and will be the subject of future research.
Acknowledgments
We thank Larry Whitesell and Jennifer Buck for providing liver tissue
samples used for the experiments.
|
|
|---|
This article has been cited by other articles:
![]() |
O. Kolokythas, T. Gauthier, A. T. Fernandez, H. Xie, B. A. Timm, C. Cuevas, M. K. Dighe, L. M. Mitsumori, M. F. Bruce, D. A. Herzka, et al. Ultrasound-Based Elastography: A Novel Approach to Assess Radio Frequency Ablation of Liver Masses Performed With Expandable Ablation Probes: A Feasibility Study J. Ultrasound Med., June 1, 2008; 27(6): 935 - 946. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Mulier, Y. Ni, L. Frich, F. Burdio, A. L. Denys, J.-F. De Wispelaere, B. Dupas, N. Habib, M. Hoey, M. C. Jansen, et al. Experimental and Clinical Radiofrequency Ablation: Proposal for Standardized Description of Coagulation Size and Geometry Ann. Surg. Oncol., April 1, 2007; 14(4): 1381 - 1396. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Wilson, Q. Chen, J. A. Zagzebski, T. Varghese, and L. VanMiddlesworth Initial clinical experience imaging scatterer size and strain in thyroid nodules. J. Ultrasound Med., August 1, 2006; 25(8): 1021 - 1029. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Varghese, U. Techavipoo, J. A. Zagzebski, and F. T. Lee Jr Impact of Gas Bubbles Generated During Interstitial Ablation on Elastographic Depiction of In Vitro Thermal Lesions J. Ultrasound Med., April 1, 2004; 23(4): 535 - 544. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |