Feasibility and potential benefits of defining the internal gross tumor volume of hepatocellular carcinoma using contrast-enhanced 4D CT images obtained by deformable registration
© Xu et al.; licensee BioMed Central Ltd. 2014
Received: 21 April 2014
Accepted: 24 September 2014
Published: 16 October 2014
To study the feasibility and the potential benefits of defining the internal gross tumor volume (IGTV) of hepatocellular carcinoma (HCC) using contrast-enhanced 4D CT images obtained by combining arterial-phase (AP) contrast-enhanced (CE) 3D CT and non-contrast-enhanced (NCE) 4D CT images using deformable registration (DR).
Ten HCC patients who had received radiotherapy beforehand were selected for this study. The following CT simulation images were acquired sequentially: NCE 4D CT in free breathing, NCE 3D CT and APCE 3D CT in end-expiration breath holding. All 4D CT images were sorted into ten phases according to breath cycle (CT00 ~ CT90). Gross tumor volumes (GTVs) were contoured on all CT images and the IGTV-1 was obtained by merging the GTVs in each phase of 4D CT images. The GTV on the APCE 3D CT image was deformably registered to each 4D CT phase image according to liver shape using RayStationTM 22.214.171.124 version treatment planning system. The IGTV-DR was obtained by merging the GTVs after DR on the 4D CT images. Volume differences among the GTVs and between the IGTV-1 and the IGTV-DR were compared.
The edge of most lesions could be definitively identified using APCE 3D CT images compared to NCE 4D and 3D CT images. The GTV volume on APCE 3D CT images increased by an average of 34.79% (P < 0.05). There was no significant difference among the GTV volumes obtained using NCE 4D and 3D CT images (P > 0.05). The GTV volumes after DR on 4D CT different phase images increased by an average of 36.29% (P < 0.05), as was observed using the APCE 3D CT image (P > 0.05). Lastly, the volume of IGTV-DR increased by an average of 19.91% compared to that of IGTV-1 (P < 0.05).
NCE 4D CT imaging alone has the potential risk of missing a partial volume of the HCC. The combination of APCE 3D CT and NCE 4D CT images using the DR technique improved the accuracy of the definition of the IGTV in HCC.
The liver shape and position vary significantly during breathing. Indeed, it has been shown that liver motion can reach 0.5 ~ 5.0 cm during breathing -. Therefore, it is imperative to accurately determine the location of hepatocellular carcinoma (HCC) to improve the accuracy of radiotherapy . Indeed changes in liver position could adversely affect radiotherapy planning and treatment, which could lead to inaccurate definition of the volumes of tumor and normal tissues, errors in dosimetric calculations and potential increases in radiation toxicity -. Therefore, eliminating the position uncertainty derived from breathing motion could improve the accuracy of HCC radiotherapy administration and therefore the clinical outcome of cancer patients .
Breath motion management methods for precision radiotherapy of HCC include active breathing control (ABC), abdominal compression and four-dimensional computed tomography (4D CT) techniques. In essence, dynamic 4D CT images could reflect the movement and deformation regularity of liver and tumors by sorting the volumetric CT images according to breath cycle which are recorded and segmented applying a respiratory position management (RPM) system ,. Application of 4D CT in HCC radiotherapy has shown promising clinical outcomes. However, a contrast-enhanced (CE) 4D CT scan is lacking because of the difficulty of capturing the timing of contrast agent injection during the extended 4D CT image acquisition process . In addition, breath motions bring more artifacts to all phases of non-contrast-enhanced (NCE) 4D CT image reconstruction . It is well known that CE 3D CT image acquisition during the arterial phase yields a clear tumor edge for HCC . Therefore, improving the detection methods to determine the location of the HCC target would improve the accuracy of radiotherapy for HCC.
Medical image deformable registration (DR) technique might prove a feasible approach in the treatment of HCC patients. DR could achieve a point-to-point correspondence between two given images, namely the objective image and the reference image, by voxel deformation and translation. The desired correspondence would describe the location of each tissue element in the first image relative to the second image . DR has become the core technique of adaptive radiotherapy and recent studies have shown that it can both improve the accuracy of delineation of tumors and organs at risk (OARs) based on multi-modality imaging and assess the accumulation dose -.
In this study we investigated the feasibility and benefits of defining the individual internal gross target volume (IGTV) of HCC using CE 4D CT images obtained by combining APCE 3D CT and NCE 4D CT images using deformable registration (DR).
Materials and methods
Ten HCC patients (5 men, 5 women; ages 58 ~ 72) were studied at the Shandong Cancer Hospital following the guidelines of the ethics committee of the institution. The study were approved by the Research Ethics Board of the Shandong Cancer Hospital. All patients provided informed consent at enrollment and had their disease confirmed by radiography. None of the patients received transcatheter arterial chemoembolization (TACE). Of the ten tumors studied one was located in the left lobe of liver and the rest in the right hepatic lobe. Patients' Karnofsky performance status was higher than 80 prior to receiving initial radiotherapy.
CT acquisition was performed using a Brilliance Big Bore CT scanner (Philips Medical Systems, Highland Heights, OH) with patients in a supine position. Scans (NCE 4D CT scans in free breathing, NCE 3D and APCE 3D CT scans in end-expiration breath holding) were acquired sequentially. 4D CT images were sorted into 10 series of CT images (CT00 ~ CT90) according to the respiratory phase, with CT00 being defined as the end inspiration and CT50 the end expiration. All CT images were exported into the RayStation treatment planning system (TPS), version 126.96.36.199 (RaySearch Laboratories, Sweden) for determination of GTVs, IGTVs, liver volumes and DR.
GTVs and liver contours were delineated on all images. IGTV-1 was acquired by merging the 10 GTVs at all phases of 4D CT images. The GTV on APCE 3D CT image was deformably registered to each phase of 4D CT images and labeled as -DR. The IGTV-DR was obtained by merging the 10 GTVs-DR after registration. The volume differences amongst the GTVs and between IGTV-1 and IGTV-DR were compared. All GTVs and liver contours were delineated by the same radiation oncologist and reviewed by a second radiation oncologist for accuracy. To reduce inconsistency in the contours, standardized window settings were used for the liver (W/L 350/ 50).
DR analysis and evaluation
DR analysis was performed using MORFEUS, a finite element model (FEM)-based multi-organ deformable image registration method that has been developed by RayStation TPS . In this study, the APCE 3D CT image in end-expiration breath holding was selected as the target image and 10 phases of 4D CT images were selected as the reference images respectively. The target image was registered to different 4D CT phase images. The overall registration was performed in the CT region which included the size of the entire patient's external body. The local registration was processed by tracking the liver volume as the organ of interest. After registration, the GTV on APCE 3D CT image was mapped to all phases of 4D CT images and 10 GTVs-DR were generated. Brock et al. evaluated the accuracy of the algorithm and reported that the precision of the vessel bifurcation identification was less than 0.1 cm in the left-right (LR), anterior-posterior (AP) and superior-inferior (SI) directions, with the accuracy of the DR of the liver being lower than 0.2 cm in each direction .
The data were analyzed using the SPSS 17.0 software package (IBM, Chicago, IL). The GTVs obtained from the different scans before and after DR were compared using one way analysis of variance. The paired t- test was used for the group comparison of GTVs and IGTVs before and after DR. A P- value <0.05 was considered statistically significant.
Comparison of GTVs before DR
A comparison of GTVs obtained from 3D and 4D CT scans in different breathing phases before and after DR ( (_x ± S)
GTV before DR (cm3)
GTV after DR (cm3)
Comparison of GTVs after DR
Comparison of IGTVs before and after DR
The tumor target can move and deform significantly during respiration, making the accurate definition of liver lesions a critical step for high-precision radiotherapy ,. 4D CT imaging is increasingly being used in radiation oncology to account for the effect of breathing on organ and tumor position in the thorax and abdomen. With this technique, the patient can breathe naturally during the 4D image acquisition process, making the therapy more comfortable. Therefore, 4D CT scanning in HCC radiotherapy could improve the accuracy of IGTV definition, reducing the margins of planning target volumes, and sparing more normal liver tissue from radiation. In summary, the 4D CT scanning technique would improve the accuracy of dose delivery and potentially result in better clinical outcomes for patients ,.
Conventionally, arterial-phase enhancement and venous-phase "washout" images are used in liver cancer patients to diagnose HCC . In this study, we found that the tumor areas on the NCE 3D and 4D CT scans were not accurately seen due to a lack of image contrast, which makes the determination of tumor boundaries difficult and impairs the precise delivery of radiotherapy. In addition, we found that the volumes of GTVs obtained using enhanced 3D CT scans were on average higher than those obtained using non-enhanced CT scans, consistent with previous reports . Currently, CE 4D CT scanning is not fully developed . When the patients are administered complete TACE, the location and extent of the lesions on 4D CT images can be determined with the help of iodized oil. In contrast, when the patients receive incomplete TACE or initial treatment alone, there is a greater uncertainty in the definition of tumor volumes based on NCE 4D CT images alone . Therefore, if the GTV generated on an artifact-free APCE 3D CT image could be mapped to all phase images of 4D CT with DR, some of the drawbacks of target delineation could be overcome.
MORFEUS allows for the image registration of both single and multiple organs. The latter can be challenging because of the difficulty of modeling the interface between organs using splines, fluid-flow and optical flow. This DR algorithm based on FEM allows the organs to be explicitly deformed by assigning the material properties and determining the boundary conditions. For some organs, including the liver, the boundary representation is consistent between different imaging modalities although the internal image intensity representation may vary. The feasibility and accuracy of MORFEUS was validated on MR thoracic and abdominal images from healthy volunteers at inspiration and expiration. For the lung and liver deformation, the average accuracy as measured by tracking visible bifurcations in the LR, AP and SI directions was 0.19 cm, 0.28 cm and 0.17 cm, respectively, with an average accuracy vector magnitude of 0.44 cm . A FEM-based algorithm has been previously described for single organ DR. In a recent study, Ferrant et al. used this algorithm to register 3D intraoperative MR images of the brain and they reported that surface-based landmark discrepancies were reduced from 1.0 cm to less than 0.1 cm . Likewise, Liang et al. tested the algorithm on a simulated deformation of the rectum and found that the error in registration was reduced from 0.9 cm based on an initial boundary point correspondence to 0.16 cm after optimization . Zhang et al. also used this algorithm in a deformable lung model based on negative surface pressure . Together, these reports support the accuracy and reliability of organ registration using MORFEUS even with internal organ motion and deformation.
Based on the above, the actual geometric verification of the mapped GTVs on 4D CT images, such as the shape and position of the GTV volumes, had not been performed. Apparently, the conventional way to verify the geometric change of the registered volumes is to introduce additional images, such as MR images, as a reference. MR images can yield high soft-tissue resolution. However, they require a long acquisition time and are susceptible to respiration and gastrointestinal peristalsis making their acquisition challenging. Moreover, a certain fusion error could inevitably exist when registering the MR and CT images. 4D MR would be expected to enable direct measurement of the tumor motion without contrast.
In this study, DR was performed in different series of images with the same imaging modality. Specifically, a CE breath-holding CT image was obtained followed by a 4D CT image. Therefore, the CT image sets had the same coordinate system and patient position. Thus, the registration error would be acceptable for the target delineation and mapping analysis. The GTV on APCE 3D CT images was deformed to match all the respiratory phases of the 4D CT images. The IGTV obtained after DR contained more tumor motion information and further improved the delineation accuracy of the tumor. To ensure the registration accuracy of whole liver, all 3D CT scans were acquired with the patient holding breath so the liver contours correlated with those on the single-phase 4D CT images. Compared with the IGTV generated from NCE 4D CT images, the IGTV obtained after DR can significantly increase the accuracy of target definition.
Beddar et al. investigated enhanced 4D CT imaging of liver tumors using synchronized intravenous contrast agent injection. Before scanning, the patient's intravenous line was connected to an automatic injector. The timing of the enhancement phase was programmed during 4D CT image acquisition to obtain good visualization of the tumor. However, the method was invasive due to the intervention with iodine oil. In addition, the optimal time window was difficult to capture as the automatic program would not work well with all patients. In contrast, as the DR technique is incorporated into modern TPS, the technique used in our study could be easily applied.
In summary, here we combined conventional target delineation protocols with image post-processing techniques to improve the accuracy of the definition of the IGTV of HCC using 4D CT. We showed that this protocol was more accurate in the determination of tumor boundaries. In future studies, we will evaluate the physiological impact of 4D radiotherapy administration.
Here we show that contrast-enhanced 4D CT images, acquired by combining arterial-phase contrast-enhanced 3D CT and non-enhanced 4D CT images based on the DR technique could improve the definition of the IGTV of HCC and have the potential ability to improve the success of radiotherapy.
This research was supported by the National Natural Science Fund of China (No. 81301936 and 81272699).
- Balter JM, Dawson LA, Kazanjian S, McGinn C, Brock KK, Lawrence T, Ten Haken R: Determination of ventilatory liver movement via radiographic evaluation of diaphragm position. Int J Radiat Oncol Biol Phys 2001, 51: 267-270. 10.1016/S0360-3016(01)01649-2View ArticlePubMedGoogle Scholar
- Davies SC, Hill AL, Holmes RB, Halliwell M, Jackson PC: Ultrasound quantitation of respiratory organ motion in the upper abdomen. Br J Radiol 1994, 67: 1096-1102. 10.1259/0007-1285-67-803-1096View ArticlePubMedGoogle Scholar
- Suramo I, Paivansalo M, Myllyla V: Cranio-caudal movements of the liver, pancreas and kidneys in respiration. Acta Radiol Diagn (Stockh) 1984, 25: 129-131.Google Scholar
- Milano MT, Constine LS, Okunieff P: Normal tissue toxicity after small field hypofractionated stereotactic body radiation. Radiat Oncol 2008, 3: 36. 10.1186/1748-717X-3-36PubMed CentralView ArticlePubMedGoogle Scholar
- Kuo HC, Liu WS, Wu A, Mah D, Chuang KS, Hong L, Yaparpalvi R, Guha C, Kalnicki S: Biological impact of geometric uncertainties: what margin is needed for intra-hepatic tumors? Radiat Oncol 2010, 5: 48. 10.1186/1748-717X-5-48PubMed CentralView ArticlePubMedGoogle Scholar
- Balter JM, Lam KL, McGinn CJ, Lawrence TS, Ten Haken RK: Improvement of CT-based treatment-planning models of abdominal targets using static exhale imaging. Int J Radiat Oncol Biol Phys 1998, 41: 939-943. 10.1016/S0360-3016(98)00130-8View ArticlePubMedGoogle Scholar
- Ten Haken RK, Balter JM, Marsh LH, Robertson JM, Lawrence TS: Potential benefits of eliminating planning target volume expansions for patient breathing in the treatment of liver tumors. Int J Radiat Oncol Biol Phys 1997, 38: 613-617. 10.1016/S0360-3016(97)00009-6View ArticlePubMedGoogle Scholar
- Yang W, Fraass BA, Reznik R, Nissen N, Lo S, Jamil LH, Gupta K, Sandler H, Tuli R: Adequacy of inhale/exhale breathhold CT based ITV margins and image-guided registration for free-breathing pancreas and liver SBRT. Radiat Oncol 2014, 9: 11. 10.1186/1748-717X-9-11PubMed CentralView ArticlePubMedGoogle Scholar
- Hugo GD, Rosu M: Advances in 4D radiation therapy for managing respiration: part I - 4D imaging. Z Med Phys 2012, 22: 258-271. 10.1016/j.zemedi.2012.06.009PubMed CentralView ArticlePubMedGoogle Scholar
- Moorees J, Bezak E: Four dimensional CT imaging: a review of current technologies and modalities. Australas Phys Eng Sci Med 2012, 35: 9-23. 10.1007/s13246-012-0124-6View ArticlePubMedGoogle Scholar
- Baron RL: Understanding and optimizing use of contrast material for CT of the liver. AJR Am J Roentgenol 1994, 163: 323-331. 10.2214/ajr.163.2.8037023View ArticlePubMedGoogle Scholar
- Yamamoto T, Langner U, Loo BW Jr, Shen J, Keall PJ: Retrospective analysis of artifacts in four-dimensional CT images of 50 abdominal and thoracic radiotherapy patients. Int J Radiat Oncol Biol Phys 2008, 72: 1250-1258. 10.1016/j.ijrobp.2008.06.1937PubMed CentralView ArticlePubMedGoogle Scholar
- Choi BI, Han JK, Cho JM, Choi DS, Han MC, Lee HS, Kim CY: Characterization of focal hepatic tumors. Value of two-phase scanning with spiral computed tomography. Cancer 1995, 76: 2434-2442. 10.1002/1097-0142(19951215)76:12<2434::AID-CNCR2820761206>3.0.CO;2-GView ArticlePubMedGoogle Scholar
- Sarrut D: Deformable registration for image-guided radiation therapy. Z Med Phys 2006, 16: 285-297.View ArticlePubMedGoogle Scholar
- Wijesooriya K, Weiss E, Dill V, Dong L, Mohan R, Joshi S, Keall PJ: Quantifying the accuracy of automated structure segmentation in 4D CT images using a deformable image registration algorithm. Med Phys 2008, 35: 1251-1260. 10.1118/1.2839120PubMed CentralView ArticlePubMedGoogle Scholar
- van Dam IE, Van Sornsen de Koste JR, Hanna GG, Muirhead R, Slotman BJ, Senan S: Improving target delineation on 4-dimensional CT scans in stage I NSCLC using a deformable registration tool. Radiother Oncol 2010, 96: 67-72. 10.1016/j.radonc.2010.05.003View ArticlePubMedGoogle Scholar
- Huang TC, Liang JA, Dilling T, Wu TH, Zhang G: Four-dimensional dosimetry validation and study in lung radiotherapy using deformable image registration and Monte Carlo techniques. Radiat Oncol 2010, 5: 45. 10.1186/1748-717X-5-45PubMed CentralView ArticlePubMedGoogle Scholar
- Jung SH, Yoon SM, Park SH, Cho B, Park JW, Jung J, Park JH, Kim JH, Ahn SD: Four-dimensional dose evaluation using deformable image registration in radiotherapy for liver cancer. Med Phys 2013, 40: 011706. 10.1118/1.4769427View ArticlePubMedGoogle Scholar
- Brock KK, Sharpe MB, Dawson LA, Kim SM, Jaffray DA: Accuracy of finite element model-based multi-organ deformable image registration. Med Phys 2005, 32: 1647-1659. 10.1118/1.1915012View ArticlePubMedGoogle Scholar
- Hallman JL, Mori S, Sharp GC, Lu HM, Hong TS, Chen GT: A four-dimensional computed tomography analysis of multiorgan abdominal motion. Int J Radiat Oncol Biol Phys 2012, 83: 435-441. 10.1016/j.ijrobp.2011.06.1970View ArticlePubMedGoogle Scholar
- Habermehl D, Debus J, Ganten T, Ganten MK, Bauer J, Brecht IC, Brons S, Haberer T, Haertig M, Jäkel O, Parodi K, Welzel T, Combs SE: Hypofractionated carbon ion therapy delivered with scanned ion beams for patients with hepatocellular carcinoma - feasibility and clinical response. Radiat Oncol 2013, 8: 59. 10.1186/1748-717X-8-59PubMed CentralView ArticlePubMedGoogle Scholar
- Rietzel E, Chen GT, Choi NC, Willet CG: Four-dimensional image-based treatment planning: target volume segmentation and dose calculation in the presence of respiratory motion. Int J Radiat Oncol Biol Phys 2005, 61: 1535-1550. 10.1016/j.ijrobp.2004.11.037View ArticlePubMedGoogle Scholar
- Liu J, Wang JZ, Zhao JD, Xu ZY, Jiang GL: Use of combined maximum and minimum intensity projections to determine internal target volume in 4-dimensional CT scans for hepatic malignancies. Radiat Oncol 2012, 7: 11. 10.1186/1748-717X-7-11PubMed CentralView ArticlePubMedGoogle Scholar
- Befeler AS, Di Bisceglie AM: Hepatocellular carcinoma: diagnosis and treatment. Gastroenterology 2002, 122: 1609-1619. 10.1053/gast.2002.33411View ArticlePubMedGoogle Scholar
- Fujishima T, Yoshida H, Obi S, Shiina S, Kanda M, Tateishi R, Akamatsu M, Koike Y, Sato S, Teratani T, Kawabe T, Shiratori Y, Omata M: Analysis of factors influencing hepatocellular carcinoma detection: efficient use of computed tomography during arterial portography and during hepatic arteriography. J Gastroenterol 2005, 40: 266-273. 10.1007/s00535-004-1533-4View ArticlePubMedGoogle Scholar
- Beddar AS, Briere TM, Balter P, Pan T, Tolani N, Ng C, Szklaruk J, Krishnan S: 4D-CT imaging with synchronized intravenous contrast injection to improve delineation of liver tumors for treatment planning. Radiother Oncol 2008, 87: 445-448. 10.1016/j.radonc.2007.12.009View ArticlePubMedGoogle Scholar
- Oh D, Lim do H, Park HC, Paik SW, Koh KC, Lee JH, Choi MS, Yoo BC, Lim HK, Lee WJ, Rhim H, Shin SW, Park KB: Early three-dimensional conformal radiotherapy for patients with unresectable hepatocellular carcinoma after incomplete transcatheter arterial chemoembolization: a prospective evaluation of efficacy and toxicity. Am J Clin Oncol 2010, 33: 370-375. 10.1097/COC.0b013e3181b0c298View ArticlePubMedGoogle Scholar
- Ferrant M, Nabavi A, Macq B, Jolesz FA, Kikinis R, Warfield SK: Registration of 3-D intraoperative MR images of the brain using a finite-element biomechanical model. IEEE Trans Med Imaging 2001, 20: 1384-1397. 10.1109/42.974933View ArticlePubMedGoogle Scholar
- Liang J, Yan D: Reducing uncertainties in volumetric image based deformable organ registration. Med Phys 2003, 30: 2116-2122. 10.1118/1.1587631View ArticlePubMedGoogle Scholar
- Zhang T, Keller H, Jeraj R, Manon R, Welsh J, Patel R, Fenwick J, Mehta M, Mackie TR, Paliwal B: Breathing synchronized delivery-a new technique for radiation treatment of the targets with respiratory motion. Int J Radiat Oncol Biol Phys 2003, 57: S185-S186. 10.1016/S0360-3016(03)00980-5View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.