|Year : 2011 | Volume
| Issue : 2 | Page : 78-84
Feasibility of deformation-independent tumor-tracking radiotherapy during respiration
Seonkyu Kim1, Myonggeun Yoon1, Dong Ho Shin1, Dongwook Kim1, Sangyeob Lee1, Se Byeong Lee1, Sung Yong Park1, Sang Hyuk Song2
1 Proton Therapy Center, National Cancer Center, Goyang, Korea
2 Seoul National University Hospital, Seoul, Korea
|Date of Submission||03-Nov-2010|
|Date of Decision||21-Dec-2010|
|Date of Acceptance||14-Jan-2011|
|Date of Web Publication||18-Apr-2011|
Proton Therapy Center, National Cancer Center, 809 Madu 1-dong, Ilsandong-gu, Goyang, 411-769
Source of Support: None, Conflict of Interest: None
| Abstract|| |
To evaluate the feasibility of tumor-tracking radiotherapy that does not consider tumor deformation during respiration. Four-dimensional computed tomography (4D-CT) data, which considers 10 phases of the respiration cycle, were acquired in 4 patients with lung cancer and 4 patients with liver cancer. Initial treatment plans were established at the end of the inhalation phase (phase 1). As a simulation of deformation-free tumor-tracking radiotherapy, the beam center of the initial plan was moved to the tumor center for all other phases, and the tumor shape acquired from phase 1 was used for all 10 phases. The feasibility of this method was analyzed based on assessment of equivalent uniform dose (EUD), homogeneity index (HI) and coverage index (COV). In photon radiation treatment, movement-induced dose reduction was not particularly significant, with 0.5%, 17.3% and 2.8% average variation in EUD, HI and COV, respectively. In proton radiation treatment, movement-induced dose reduction was more significant, with 0.3%, 40.5% and 2.2% average variation in EUD, HI and COV, respectively. Proton treatment is more sensitive to tumor movement than is photon treatment, and that it is reasonable to disregard tumor deformation during photon therapy employing tumor-tracking radiotherapy.
Keywords: Deformation, real-time tumor tracking, respiration, 4D-CT
|How to cite this article:|
Kim S, Yoon M, Shin DH, Kim D, Lee S, Lee SB, Park SY, Song SH. Feasibility of deformation-independent tumor-tracking radiotherapy during respiration. J Med Phys 2011;36:78-84
|How to cite this URL:|
Kim S, Yoon M, Shin DH, Kim D, Lee S, Lee SB, Park SY, Song SH. Feasibility of deformation-independent tumor-tracking radiotherapy during respiration. J Med Phys [serial online] 2011 [cited 2019 Apr 25];36:78-84. Available from: http://www.jmp.org.in/text.asp?2011/36/2/78/79691
| Introduction|| |
The objective of radiation therapy is to administer radiation to a tumor, with minimization of radiation received by adjacent normal tissues. Treatment planning is used to accurately target a tumor and to optimize the administration of radiation, which is typically three-dimensional conformal radiation therapy (3D-CRT) or intensity-modulated radiation therapy (IMRT). 3D-CRT and IMRT require high accuracy, because a multileaf collimator (MLC) must be shaped to match the tumor. , The shape of MLC is based on stationary medical images, such as those provided by computed tomography (CT), magnetic resonance imaging (MRI) and positron emission tomography (PET). Patient respiration or movement can change the location and shape of a tumor in a lung or liver. To compensate for these changes, the tumor margin is generally considered during treatment planning. For example, liver cancer has a setup margin of at least 1.5 cm during respiration, and is considered in planning for the target margin. 
Various techniques have been developed to solve the intra-fractional motion of tumors, including the motion-encompassing method, holding of breath, respiratory gating, and real-time tumor-tracking. The motion-encompassing method requires an estimation of the position and range of tumor motion during respiration, as provided by 4D-CT. ,, However, if patient respiration is irregular, tumor size may be overestimated, and the radiation dose may be unnecessarily high. In the breath-holding method, the patient holds his/ her breath, and the position of the patient is maintained in all fractions. In the forced-shallow breathing method, respiration movements are limited by a physical plate positioned on the abdominal region. , However, this procedure may result in significant patient discomfort and poor reproducibility.
Several research centers are studying the use of respiratory gating to deal with respiratory motion during radiotherapy of thoracic and abdominal tumors.  This procedure is essentially the same as 3-D conformal therapy. In particular, imaging and treatment are synchronized with the patient's respiratory cycle, thereby reducing the margin of difference between the clinical target volume (CTV) and the planning target volume (PTV).  Respiratory-gating methods are categorized as internal or external, depending on the use of surrogates. During the respiratory cycle, the position and width of the gate are obtained by monitoring the respiratory motion of the patient using an external signal and a fiducial marker. External gating uses a surrogate marker (typically on the patient's abdomen), and internal gating employs an implanted marker. A drawback of this method is that marker and tumor motion may be different. In addition, the relationship between tumor motion and the surrogate signal is unstable and can vary over time.
Real-time tumor-tracking systems provide a method for reducing the effect of respiration-induced target motion. Such tumor-tracking using an MLC (multileaf collimator) has been introduced in previous studies; it allows easy adjustment of the size and shape of the irradiating photon beam. , The simplest approach is to move the collimator. However, in tumor-tracking radiotherapy, it is essential to know the exact tumor location and shape. The location of a tumor can be tracked using various methods, but real-time tracking of tumor shape cannot be achieved by the 2D imaging systems that are conventionally used in treatment rooms.
There have been many studies of real-time tumor-tracking, but little is known about the effect of radiation dose when changes in tumor shape are ignored during tumor-tracking radiotherapy. In the present study, we evaluated the feasibility of tumor-tracking radiotherapy used without consideration of tumor deformation.
| Materials and Methods|| |
Patient characteristics and treatment planning
Four lung cancer patients and 4 liver cancer patients were analyzed at our institution using 4D-CT, which divided the respiratory cycle into 10 phases. Phase 1 corresponds to the end of an inhalation phase (0%). [Figure 1] shows a schematic of our tumor-tracking simulation. This figure shows that the beam center of the initial plan was moved to the tumor center in other phases, using tumor shape (considered to be constant) acquired from phase 1. The targets were defined in accordance with the report of the International Commission on Radiation Units and Measurements (ICRU 50). In particular, gross tumor volume (GTV) encompassed all detectable tumors and lymph nodes as observed on CT scans. Planning target volume (PTV) included the GTV plus a 10-15-mm margin. To be consistent for tumor contouring, the same clinician manually drew all the tumor volumes in the 10 phases of CT.
|Figure 1: Schematic picture for tumor-tracking simulation. The beam center is moved to the tumor center based on the location of moving tumor|
Click here to view
[Table 1] lists the clinically relevant patient characteristics and treatment schemes. All treatment plans for lung and liver cancer patients used 3 to 5 beams for IMRT and 2 to 3 beams for Proton beam therapy (PBT). The proximal, distal and transverse margins were 2, 2 and 10 mm, respectively; and the border smoothing and smearing margins were set at 0 and 3 mm, respectively. We used relatively small margins in planning, because various uncertainties had already been included in our PTV.
Adequate assessment of the homogeneity of target volume is important for calculating the dose volume histogram (DVH). Losses in tumor-control probability will occur for cold-spot-induced dose inhomogeneity, so it is essential to properly evaluate the homogeneity of the target volume. The homogeneity index (HI) is defined as 
whereD2 and D98 represent the doses to 2% and 98% of the target volume, respectively, and D p is the prescribed dose. D98 , which may be considered the minimum dose, is the dose received by 98% of the target volume; D 2 , which may be considered the maximum dose, is the dose received by 2% of the target volume. A lower HI indicates a more homogeneous target dose. The normalized HI (nHI) is the ratio of the HI of the shifted target to the HI of the original target.
In addition to providing information on the homogeneity of radiation doses, DVHs can also be used to assess target coverage based on the coverage index (COV), defined as the percentage of tumor volume that received the prescribed dose. Ideally, tumor DVH would be a step function, with 100% of the target receiving the exact prescribed dose. However, actual DVH curves are not step functions, because of constraints imposed by tumor volume and other organs at risk (OAR). The conditions for clinically acceptable target volume coverage include (i) no more than 20% of any planning target volume (PTV) will receive > 110% of the prescribed dose; (ii) the prescribed dose is the isodose that encompasses at least 95% of the PTV; and (iii) no more than 1% of any PTV will receive < 93% of the prescribed dose.  The last two conditions indicate that coverage indices at D p and at 93% of D p should be more than 95% and 99%, respectively. The normalized COV (nCOV) is the ratio of the COV for the shifted target to the COV of the original target.
Equivalent uniform dose
When evaluating the dose homogeneity of a treatment plan, it is necessary to consider the radiobiological impact of dose inhomogeneity on target volume. The equivalent uniform dose (EUD) is one parameter that describes the relationship between homogeneity and radiobiological effect. The EUD is defined as the biologically equivalent dose that if given uniformly, would lead to the same reduction in tumor volume as the actual inhomogeneous dose distribution.  It is calculated as
The reference dose (Dref) and surviving fraction (SF 2 ) are typically set at 2 Gy and 0.5 Gy, respectively. Summation is performed over all bins of the DVH with a volume element (Vi) and a clonogenic cell density (r i), both of which are uniform. EUD is calculated as the percentage of the prescribed dose, and normalized EUD (nEUD) is the ratio of the EUD for the shifted target to the EUD of the original target.
| Results|| |
[Figure 2] shows an example of a radiotherapy plan for treatment of lung cancer using photon and proton therapies. The dose map for lung cancer proton therapy [Figure 2]c and d is clearly better than that for photon therapy [Figure 2]a and b. Similarly, the treatment plans for liver cancer show better dosimetric behavior of proton therapy compared with photon therapy [Figure 3]. All treatment plans were established with a CT set during phase 1, and satisfy conventional dosimetric factors (dose conformity, dose homogeneity and dose coverage).
|Figure 2: An example of a radiotherapy plan for treatment of lung cancer using photon and proton therapies. (a) Axial view, (b) coronal view in photon treatment; and (c) axial view, (d) coronal view in proton treatment are shown|
Click here to view
|Figure 3: An example of a radiotherapy plan for treatment of liver cancer using photon and proton therapies. (a) Axial view, (b) coronal view in photon treatment; and (c) axial view, (d) coronal view in proton treatment are shown|
Click here to view
[Figure 4] shows an example of 10 coronal images of the 10 phases of respiration in 4D-CT. This figure shows the relationship between tumor shape and location with respect to respiration phase. Although it is not particularly clear, tumor shape definitely changes throughout the 10 phases, with a change in position of up to 2-3 cm. This suggests that dosimetric factors may vary because of tumor motion during radiation therapy.
|Figure 4: Coronal slices depicting the peak-inhale phase (phase 1, top left), peak-exhale phases (phases 5-6, bottom left and top right), peak-inhale phase (phase 10, bottom right) and a few others of a 4-D computed tomography (4D CT) scan. The contour represents the tumor volume as delineated by the physician on the planning CT|
Click here to view
[Table 2] shows variation in tumor size during the 10 respiratory phases. Average tumor size varied by up to ~8%, suggesting that volume changes during respiration may need to be considered. Tumor size is largest in phases 5 to 7.
|Table 2: Planning target volume (relative to phase 1) during 10 phases of respiration|
Click here to view
[Figure 5] shows an example of dosimetric variation of DVH for photon and proton treatment of 1 lung cancer patient. DVHs were assessed when tumor-tracking radiotherapy was performed, assuming constant tumor shape (acquired from phase 1). In photon therapy [Figure 5]a, DVH changed very little, suggesting that tumor-tracking radiotherapy without consideration of tumor deformation may be effective in photon therapy [Figure 5]a. DVH values for proton therapy [Figure 5]b also show little variation from the phase-1 DVH, suggesting that this method may also be suitable for proton therapy. [Figure 6] shows the DVH for phase 1 and the average DVH of phases 2 to 10 for photon and proton treatments of lung cancer patient seen in [Figure 5]. Clearly, the DVH of phase 1 and the average DVH of phases 2 to 10 are similar, suggesting that tumor-tracking radiotherapy may be successful if lung tumor deformation is ignored.
|Figure 5: An example of dosimetric variation of DVH for (a) photon treatment and (b) proton treatment of 1 lung cancer patient when tumor-tracking radiotherapy was performed, assuming constant tumor shape (acquired from phase 1)|
Click here to view
|Figure 6: DVH for phase 1 and the average DVH of phases 2 to 10 for (a) photon treatment and (b) proton treatment of lung cancer patient seen in Fig. 5, which are similar, indicating that the same level of DVH can be acquired based on the proposed method|
Click here to view
[Figure 7] shows an example of dosimetric variation of DVH for photon and proton treatments of 1 liver cancer patient. The results are similar to those for lung cancer [Figure 5], indicating that tumor deformation during respiration can also be ignored for radiotherapy of liver cancer. [Figure 8] shows that the DVH of phase 1 and the average DVH of phases 2 to 10 are similar for liver cancer, as was the case with lung cancer [Figure 6].
|Figure 7: An example of dosimetric variation of DVH for (a) photon treatment and (b) proton treatment of 1 liver cancer patient when tumor-tracking radiotherapy was performed, assuming constant tumor shape (acquired from phase 1)|
Click here to view
|Figure 8: DVH for phase 1 and the average DVH of phases 2 to 10 for (a) photon treatment and (b) proton treatment of liver cancer patient seen in Figure 7, which are similar, indicating that the same level of DVH can be acquired based on the proposed method|
Click here to view
[Table 3] shows the averaged dosimetric changes during the 10 phases of respiration for lung and liver cancers together when tumor deformation is ignored. For photon therapy, the greatest change in minimum dose was 23.8% (during phase 6), but the average dose varied by only about 0.5%. For proton therapy, the greatest change in minimum dose was 38.6% (during phase 6), and average dose varied by about 1.0%.
|Table 3: Radiation dose for the 10 phases of respiration, relative to phase 1; and the averaged data both for lung and liver cancers in photon and proton treatments|
Click here to view
Finally, we calculated 3 dosimetric factors, nEUD, nHI, and nCOV [Table 4], to assess the role played by tumor movement during photon and proton therapies. In photon therapy, the average variation of nEUD, nHI and nCOV was 0.5%, 17.3% and 2.8%, respectively. In proton therapy, the average variation of nEUD, nHI and nCOV was 0.3%, 40.5% and 2.2%, respectively.
|Table 4: Relative dosimetric factors for 10 phases compared to the factor for phase 1 in photon and proton treatments|
Click here to view
| Discussion|| |
In the present study, we examined the dosimetric effect of the radiation treatment plan during tumor-tracking radiotherapy, when tumor deformation is ignored. There have been numerous previous studies of tumor deformation during radiotherapy. In particular, Chhatkuli et al. simulated the deformation of lungs and lung tumors during inspiration using a mesh-free simulation technique termed the moving particle semi-implicit method (MPS). The cited authors modeled deformation by considering lung tissues to be homogenous, isotropic and visco-elastic. The regional deformation of lung tumors with a superior-inferior (SI), right-left (RL) and anterior-posterior (AP) orientation was compared with experimental CT data taken at the end of inspiration. The authors compared their numerical results with experimental data obtained by tracking the movement of gold fiducial markers. Their results showed that deformation varied from less than 5 mm in the upper region to over 20 mm in the lower segment. For the tumor, however, both the experimental and numerical predictions showed that there was no volumetric change in tumor shape by the end of inspiration.
Thus, our experimental results with lung and liver cancer patients, along with previous simulation data on lung cancer, indicate that tumor volume and shape do not have significance throughout inspiration. However, both experimental results and numerical simulations indicate that the center point of a lung tumor was 5 mm lower (along the SI direction) at the end of inspiration. This suggests that tumor-center-tracking radiotherapy which assumes a constant tumor shape may be feasible and supports the proposal that tumor deformation during respiration may be disregarded in radiation therapy.
| Conclusion|| |
A previous simulation study showed it was possible to establish a radiation treatment plan by tracking the movement of the tumor center, ignoring changes in tumor shape during respiration. Our experimental results with 4 lung cancer and 4 liver cancer patients suggest that proton therapy is more sensitive to tumor movement than is photon therapy, and that tumor deformation may be safely disregarded in tumor-tracking radiotherapy that employs photon therapy.
| References|| |
|1.||Saw CB, Ayyangar KM, Enke CA. MIMiC-based IMRT-Part I. Med Dosm 2001;26:111-1. |
|2.||Saw CB, Komanduri KM, Enke. CA. MLC-based IMRT-Part II. Med Dosm 2001;26:111-2. |
|3.||Wagman R, Yorke E, Ford E, Giraud P, Mageras G, Minsky B, et al.Respiratory gating for liver tumors: Use in dose escalation. Int J Radiat Oncol Biol Phys 2003;55:659-8. |
|4.||Underberg RW, Lagerwaard FT, Slotman BJ, Cuijpers JP, Senan S. Benefit of respiration-gated stereotactic radiotherapy for stage I lung cancer: An analysis of 4DCT datasets. Int J Radiat Oncol Biol Phys 2005;63:253-60. |
|5.||Vedam SS, Keall PJ, Kini VR, Mostafavi H, Shukla HP, Mohan R. Acquiring a four-dimensional computed tomography dataset using an external respiratory signal. Phys Med Biol 2003;48:45-62. |
|6.||Low DA, Nystrom M, Kalinin E, Parikh P, Dempsey JF, Bradley JD, et al. A method for the reconstruction of 4-dimensional synchronized CT scans acquired during free breathing. Med Phys 2003;30:1254-63. |
|7.||Lax I, Blomgren H, Naslund I, Svanstrom R. Stereotactic radiotherapy of malignancies in the abdomen.Methodological aspects. Acta Oncol 1994;33:677-83. |
|8.||Blomgren H, Lax I, Naslund I, Svanstrom R. Stereotactic high dose fraction radiation therapy of extracranial tumors using an accelerator. Acta Oncol 1995;34:861-70. |
|9.||Vedam SS, Kini VR, Keall PJ, Ramakrishnan V, Mostafavi H, Moha R. Quantifying the predictability of diaphragm motion during respiration with a noninvasive external marker. Med Phys2003;30:505-13. |
|10.||Keall PJ, Gig S, James M. The management of respiratory motion in radiation oncology report of AAPM Task Group 76. Med Phys 2006;33:3874-900 . |
|11.||Keall PJ, Kini VR, Vedam SS, Mohan R. Motion adaptive X-ray therapy: A feasibility study. Phys Med Biol 2001;46:1-10 . |
|12.||Neicu T, Shirato H, Seppenwoolde y, Jiang SB. Average Tumor Trajectory for Lung Patients Phys Med Biol 2003;48:587-98. |
|13.||Wu Q, Mohan R, Morris M, Lauve A, Schmidt-Ullrich R. Simultaneous integrated boost intensity-modulated radiotherapy for locally advanced head-and-neck squamous cell carcinomas. In: Dosimetric results. Int J Radiat Oncol Biol Phys 2003;56:573-85. |
|14.||Radiation Therapy Oncology Group RTOG 0022: Phase I/II Study of Conformal and Intensity Modulated Irradiation for Oropharyngeal Cancer - CLOSED. In: Eisbruch A, Chao KS, Garden AS, editors. Philadelphia: Radiation Therapy Oncology Group; 2001. |
|15.||Niemierko A. Reporting and analyzing dose distributions: A concept of equivalent uniform dose. Med Phys 1997;24:103-10. |
|16.||Chhatkuli S, Kojima H, Koshizuka S, Uesaka M. Meshfree Simulation of Lung Deformation during Inspiration. Proceedings of the World Congress on Engineering and Computer ScienceVol 2. USA: WCECS; 2009. |
[Figure 1], [Figure 2], [Figure 3], [Figure 4], [Figure 5], [Figure 6], [Figure 7], [Figure 8]
[Table 1], [Table 2], [Table 3], [Table 4]