INTRODUCTION
Helical TomoTherapy is a feasible solution of image-guided radiation therapy to fulfil highly conformal intensity-modulated radiotherapy (IMRT).Reference Mackie, Kapatoes, Ruchala, Lu, Wu, Olivera, Forrest, Tome, Welsh, Jeraj, Harari, Reckwerdt, Paliwal, Ritter, Keller, Fowler and Mehta1,Reference Welsh, Lock, Harari, Tomé, Fowler, Mackie, Ritter, Kapatoes, Forrest, Chappell, Paliwal and Mehta2 The presence of the integrated online megavoltage computed tomography (MVCT) unit results in innovative approaches to adaptive radiation therapy.Reference Welsh, Lock, Harari, Tomé, Fowler, Mackie, Ritter, Kapatoes, Forrest, Chappell, Paliwal and Mehta2 Adaptive radiation therapy is a concept/technique to modify radiation treatments in a closed-loop process, utilising a systematic feedback of measurements.Reference Yan, Vicini, Wong and Martinez3 MVCT, which would be acquired just before each treatment, enables verification of daily set-up and corrections for internal organ motion. These MVCTs also serve as inputs for the adaptive module of helical TomoTherapy, the Planned Adaptive software.
The Planned Adaptive software is capable of registering manually and automatically the MVCT images generated on the day of the treatment with the kilovoltage computed tomography (kVCT) images used for treatment planning.Reference Woodford, Yartsev, Dar, Bauman and Van Dyk4 With the image registration, the volume changes in the target and/or in the normal tissues during the treatment course of TomoTherapy could be evaluated.Reference Woodford, Yartsev, Dar, Bauman and Van Dyk4–Reference Siker, Tomé and Mehta6 Based on the registered kVCT-MVCT images, a verification dose distribution is calculated to indicate the actual dose delivered to the target and the adjacent normal tissues.Reference Han, Chen, Liu, Schultheiss and Wong7 Using the verification doses and sinograms from different treatment fractions, a summation dose is generated for the partial or the whole course of treatment. Based on the dose deviation of the summation dose from the planned prescription, the adaptive modification of the rest treatment fractions could subsequently be made.Reference Woodford, Yartsev, Dar, Bauman and Van Dyk4,Reference Siker, Tomé and Mehta6,Reference Ramsey, Langen, Kupelian, Scaperoth, Meeks, Mahan and Seibert8 Although adaptive planning is considered to be a promising way to continuously adjust the dose that would be delivered to the target and organs at risk, the clinical benefit in relation to the cost of extra time and effort of re-planning remains mostly unproven.
Biologically effective uniform dose () is able to report physical dose delivery, but more importantly, it facilitates the analysis of dose-response relations, which link the radiation dose to the clinical effect.Reference Mavroidis, Lind and Brahme9 Therefore, the concept of the has been used to assess the difference between planned and delivered IMRT dose distributions.Reference Mavroidis, Plataniotis, Górka and Lind10 Previous studies have utilised for plan comparison between helical TomoTherapy and multileaf collimator–based IMRT.Reference Mavroidis, Ferreira, Shi, Lind and Papanikolaou11,Reference Su, Mavroidis and Swanson12
With the ability to retrospectively review the delivered doses to the target and critical organs by using the Planned Adaptive software plus the capability of converting the physical doses to clinical effectiveness using the , the goal of this study is to predict the radiobiological effects of the adaptive TomoTherapy using the in addition to physical dose indices (such as the dose-volume histogram (DVH) and the mean dose). Furthermore, the clinical impact of the adaptive plans in TomoTherapy would be assessed by using the and the complication-free tumour control probability, P +.
MATERIALS AND METHODS
Study candidates and adaptive software module of TomoTherapy
Helical TomoTherapy was delivered using the Hi-Art TomoTherapy unit (TomoTherapy Inc., Madison, WI). Four lung cancer patients who underwent helical TomoTherapy were randomly selected retrospectively. The general information of these four patients is listed in Table 1. Patient 4 had the original treatment plan for 16 fractions followed by the second-phase treatment plan for the remaining 9 fractions. Daily MVCT scans were performed and registered with the planning kVCT images to correct the patient set-up and internal organ motions. Evaluation of the delivered dose distribution in each treatment fraction of the TomoTherapy was performed by using the Planned Adaptive software (TomoTherapy Inc., Madison, WI).
Before the assessment of the delivered dose distribution for each fraction, image registration between the MVCT and kVCT was performed using an automatic registration function in Planned Adaptive. Woodford et al.Reference Woodford, Yartsev and Van Dyk5 suggested that with the selection of full image fusion and fine resolution will result in a low residual error in most of the cases when modifying image registration of MVCT and kVCT for lung cancer patients treated with helical TomoTherapy.
Concept of biologically effective uniform dose ()
The uniform dose that causes the same tumour control probability or normal tissue complication rate as the actual dose distribution given to the patient was evaluated using the biologically effective uniform dose, (ref. 9,13). The general definition of the can be expressed as the equation below:
where the denotes the three-dimensional dose distribution.
The radiobiological model that was used to describe the dose-response relation of tumours and organs at risk was the linear-quadratic Poisson modelReference Källman, Agren and Brahme14,Reference Ågren Cronqvist, Kallman, Brahme and Baltas15:
where the P(D) is the probability to control the tumour or induce a certain injury to an organ that is irradiated uniformly with a dose D. Since this model takes into account the fractionation effects that are introduced by the irradiation schedule, d (equals to D/n) is the dose per fraction and n is the number of fractions. D 50 is the dose which gives a response probability of 50% and γ is the maximum normalised value of the dose-response gradient. Variables α and β are the fractionation parameters of the model and account for the early and late effects, respectively. The dose-response parameters of the target and organs at risk used in this study are listed in Table 2.
D 50 is the 50% response dose, γ is the maximum normalised value of the dose-response gradient and s is the relative seriality, which characterises the volume dependence of the organ.
is the biologically effective uniform dose, which is calculated based on the radiological characteristics of the target and it is associated with the clinical benefit. is the biologically effective uniform dose, which is calculated based on the radiological characteristics of normal tissues and it is associated with the radiation-induced injury.Reference Mavroidis, Lind and Brahme9
Statistical methods for results
For four lung cancer patients, the physical doses to the target and critical organs were evaluated through the DVH at the prescription dose level for all treatment plans. The physical dose distributions were also calculated through the Planned Adaptive software for every five treatment fractions of the whole treatment course to periodically monitor the delivered dose distributions and to compare the measured dose distributions with the planned ones.Reference Han, Chen, Liu, Schultheiss and Wong7 Based on the physical dose distributions, the DVHs and the dose-response parameters of the target and the critical organs of the treatment plans and actual delivered fractions, the , tumour control probability (P B), normal tissue injury rate (P I) and complication-free tumour control probability (P +) were determined for these four lung cancer patients.Reference Mavroidis, Lind and Brahme9 Biological evaluations also included the assessment of the optimal and P + for all the treatment plans to indicate the ideal complication-free tumour control probability. The optimal P + is the maximum of the probability distribution calculated based on the (ref. 9). Additionally, the ratio of the delivered and the planned physical dose, biological doses and P + were calculated to show the efficiency of the TomoTherapy delivery physically and biologically.
Among these four patients, patient 4 had a two-phase treatment. In the first-phase treatment, because of the obvious tumour shrinkage observed in the planning target volume through the daily MVCT scans and the high cumulative dose to the ipsilateral lung, the physician decided to reduce the sizes of the irradiating fields and prescribe a second-phase plan after 16 fractions of TomoTherapy. The second-phase treatment of the patient 4 in this situation could be deemed as adaptive TomoTherapy. In order to get enough sampling points for the two-phase TomoTherapy of patient 4, the treatment plan and the actual delivered fractions were sampled for every two to five fractions instead of every five fractions.
RESULTS
Figure 1 demonstrates the ratio between the delivered and the planned physical doses to the target, that ratio calculated using the biologically effective uniform doses () and the ratio obtained by dividing the delivered to the planned P + for these four patients. Tables 3–6 summarise the physical mean dose (), , , P +, P B and P I, which were calculated based on the treatment plans and the corresponding delivered fractions for patient 1 to patient 4. These tables also indicate the optimum (in the first column) of the biologically effective uniform dose and the complication-free tumour control rate for the treatment plan of each patient.
From Figure 1a, up to 9% under-dosage of the target is observed from the ratio of the delivered to the planned physical dose over a total of 30 treatment fractions (range from 0.908 to 0.912). When evaluating the ratio between the delivered divided by the planned , which was derived with the consideration of the biological characteristics of the target,Reference Mavroidis, Lind and Brahme9,Reference Mavroidis, Lind, Van Dijk, Koedooder, De Neve, De Wagter, Planskoy, Rosenwald, Proimos, Kappas, Claudia, Benassi, Chierego and Brahme13 the biologically effective uniform dose that was delivered to the target is lower up to about 6%. These inferior ratios (range from 0.937 to 0.939) of over 30 fractions result in up to 10% lower P + compared to the planned P + throughout the treatment course of patient 1. According to Table 3, the optimal is 89.00 Gy and the optimal is 18.00 Gy, resulting in a P + of 92.2%.
In Figure 1b, the results of patient 2 show that the ratio of the delivered dose range from 0.969 to 0.970 in comparison with the planned prescribed dose. This results in deviations in from 0.978 to 0.981 over 25 treatment fractions of this patient. As a result, up to 4% lower P + compared to that of the treatment plan is observed during the treatment course. According to Table 4, the optimal of the treatment plan for patient 2 is 67.95 Gy and the optimal is 21.25 Gy, resulting in a P + of 72.1%.
For patient 3, the ratios between the delivered and the planned dose to the target range from 0.903 to 0.915 in 20 treatment fractions, resulting in ratios between 0.928 and 0.937 of the delivered to the planned as shown in Figure 1c. Furthermore, according to Figure 1c, up to 7% lower P + is observed at the end of 20th fraction for patient 3. Table 5 demonstrates that the optimal of the treatment plan for patient 3 is 70.00 Gy and the optimal is 19.45 Gy, resulting in a P + of 79.9%.
Figure 1d shows the results of patient 4, who had a second-phase treatment after the 16th of 30 treatment fractions. In the first-phase treatment, the ratio of the physical dose range from 0.953 to 0.956, which results in the ratio of the to vary from 0.978 to 0.983. As a result, up to 4.5% lower P + was calculated in the first-phase delivery of patient 4. In the second-phase treatment, the ratio of the physical dose increased to the values of 0.972 and the ratio of the rose up to 0.989. Therefore, the ratio of the P + grew up to 0.977 in the second-phase delivery. In Table 6, the optimal of 62.00 Gy, the optimal of 24.25 Gy and the P + of 52.6% for the first-phase treatment, and the optimal of 74.00 Gy, the optimal of 20.75 Gy and the P + of 79.4% for the second-phase treatment are listed.
Figure 2 shows the changes in terms of P + within the first-phase and the second-phase treatments and between the two treatment phases of patient 4. The P + of the second-phase treatment is significantly higher by about 13.2% as compared to the average P + in the first-phase treatment. On the other hand, the variation among the treatment fractions in terms of P + in the first-phase treatment is small. Also, there was no considerable change in P + within the second-phase treatment. The substantial increased P + in the second-phase TomoTherapy in comparison with the first-phase treatment indicates the clinical benefits of adaptive TomoTherapy for this lung patient.
DISCUSSION
Helical TomoTherapy of four lung cancer patients has been evaluated comprehensively using the Planned Adaptive software as well as the biologically effective uniform dose. In addition to the dosimetric evaluation, which is based on DVH and dose statistics (e.g. mean dose) obtained from the Planned Adaptive module, calculation was also employed to quantify dose differences in terms of changes in the expected clinical outcome.Reference Mavroidis, Ferreira, Papanikolaou, Svensson, Kappas, Lind and Brahme16 Quantification of the variation in tumour control probability and normal tissue complication rate for the delivered fractions would be more useful in adaptive planning during radiation therapy. The results in Figure 1 show the ratio of the delivered to the planned dose distributions in terms of the physical dose, the biological dose and the P +. Over the whole course of treatment using TomoTherapy, the P + ratio curves revealed equal or more pronounced changes in comparison with the physical dose ratio curves for these four lung cancer patients. From Figure 1a,c, the clinical effect of different doses per fraction could be assessed by comparing results of patient 1 and patient 3. In Figure 1a,c, both curves of the physical dose ratios for patient 1 and patient 3, respectively, varied within a very similar range (in Figure 1a: 0.908–0.912; in Figure 1c: 0.903–0.915), indicating no significant difference (p = 0.5950). On the other hand, the P + ratios of patient 3, who had a larger dose-per-fraction, were significantly higher (p = 0.0019) than those of patient 1. The greater biological impact of the larger dose-per-fraction could be foundReference Song, Battista and Van Dyk17,Reference Song, Chiu, Bauman, Lock, Rodrigues, Ash, Lewis, Fenster, Battista and Van Dyk18 by using the P + evaluation, which was calculated with consideration of the dose-response parameters of the targets and organs at risk.
From the results in Figure 1 and Tables 3–6, the estimated doses delivered to the targets of the four lung cancer patients were different in trend. The same holds for the biologically effective uniform dose and complication-free tumour control rate over the whole course of treatment. Also, previous reports have proven that it is difficult to predict the volume changes in the target and treatment results based on the patient, treatment schedule or tumour characteristic.Reference Woodford, Yartsev, Dar, Bauman and Van Dyk4,Reference Siker, Tomé and Mehta6 Therefore, the individual evaluation of the dosimetric and clinical effects for each lung cancer patient treated with TomoTherapy using plus P + is necessary. Figure 1a shows the flat trends in the three curves of physical dose, biological dose and P + ratios after the 10th fraction (i.e., the 2nd week) of TomoTherapy for patient 1. The consistent suboptimal ratios over the whole treatment interval when compared to the treatment plan indicate the lack and, thus, the necessity of adaptive planning. The re-optimization of adaptive planning usually results in an improved tumour control probability combined with a limited normal tissue complication rate.Reference Yan, Vicini, Wong and Martinez3 From the results of patient 2 in Figure 1b, it is shown that the decreasing trend of the P + ratios for the first 15 fractions contributed to the lowest complication-free tumour control rate observed at 15th fraction. The dosimetric evaluation of the Planned Adaptive software and the assessment of the clinical effects using the and the P + identify the dose deviations and, as a result, the lower P +. With this information in the middle of the treatment, the physician could re-check the target contour, and the medical physicist could investigate the accuracy of dose delivery and the necessity of planning adaptation. For patient 3, the continuously decreasing ratios of the physical dose, the biological dose and the P + between delivered and planned dose distributions showed the suboptimal dose delivering of the treatment plan (cf. Figure 1c). With surveillance of physical doses and biological doses during TomoTherapy, the under-dosage of the target, resulting in lower tumour control, should be identified and corrected in time.
The results of patient 4 show the potential benefit of adaptive TomoTherapy as it is denoted by the significantly increased complication-free tumour control rate, resulting from the notably decreased normal tissue complication possibility. After the first-phase TomoTherapy (16 fractions), the physician decided to have a second-phase plan because of the considerably reduced tumour volume observed from the daily MVCTs and the high dose to the ipsilateral lung. From Figure 1d and Table 6, the P + ratios of the second-phase treatment (from 17th fraction to the end) increased by 1.2% as compared to the P + ratios of the first-phase treatment. When comparing the P + of the two treatment phases, P + of the second-phase TomoTherapy increased by 13.2% to the P + of the first-phase TomoTherapy. The reason for this significant increase in P + was mainly due to the notably lower P I (8.6%) of the treatment plan in the second-phase treatment. The results of patient 4 with adaptive TomoTherapy are in line with the findings of Woodford et al.Reference Woodford, Yartsev, Dar, Bauman and Van Dyk4 Woodford et al.Reference Woodford, Yartsev, Dar, Bauman and Van Dyk4 suggested that adaptive planning can yield significant improvements in cumulative doses to organs at risk if the gross tumour volume decreases considerably. More evaluations are required on the need for adaptive planning especially for the cases in which the organs at risk impose significant dose limitation.Reference Woodford, Yartsev, Dar, Bauman and Van Dyk4 In comparison to previous reports,Reference Woodford, Yartsev, Dar, Bauman and Van Dyk4,Reference Siker, Tomé and Mehta6 our results not only show the changes in tumour volume with respect to the dosimetric effects when using the Planned Adaptive software, but it also demonstrate the increased tumour control rate and, thus, the escalating clinical effectiveness of the adaptive planning during the treatment course of TomoTherapy by using the . Several studies have suggested the benefits of adaptive TomoTherapy by using the Planned Adaptive software to monitor the tumour volume changes and the margins around the tumour during the course of treatment.4,8,19 In this study, the use of the leads to a closer association of the DVH and dose statistics with tumour control or normal tissue injury. With more information concerning the clinical impact of the delivered treatment, it would be possible to obtain a plan that is better tailored to the individual patient. For an adaptive schedule to be clinically feasible with respect to the clinical workload, the physicians could adjust the plan and the target coverage based on the predicted changes in complication-free tumour control rate at certain fractions. Therefore, the time and resources spent in re-planning would be justified by the increasing tumour control rate and/or decreasing normal tissue complication possibility.
Another advantage of the evaluation for the treatment plans and delivered fractions is the reports of the optimal P + calculated from the optimal compensation between the P B and P I indices. With the data of the optimal prescription level based on the tumour control probability and the normal tissue injury rate, it would be easier to pick up the optimal treatment plans before the treatment and to monitor the difference in the P + between the actual delivered dose and the optimal amount. It is also important to assess the available therapeutic window in terms of tumour control and normal tissue injury probabilities when considering about the necessity of using the adaptive planning. From the results in Tables 3–6, comparing the optimum in the first columns to the rest of the columns, large gaps can be observed in the and the P + between the optimal situations and the treatment plans for all four lung cancer patients. Even larger differences in the and the P + were found between the optimums and the actually delivered radiation. With the pre-knowledge of the optimal level of the and the P +, adaptive planning should be performed to maximise the tumour control probability for a clinically acceptable normal tissue complication rate.Reference Yan, Vicini, Wong and Martinez3,Reference Mavroidis, Lind and Brahme9
The under-dosage of the target evaluated by using the Planned Adaptive software was as large as over 9% in both patient 1 and patient 3. Possible reasons for the notable dose deviations evaluated by the Planned Adaptive software are set-up error, organ motion and errors due to the image registration. Han et al.Reference Han, Chen, Liu, Schultheiss and Wong7 assessed the actual dose variation to the target and the critical organs of patients treated with TomoTherapy using the Planned Adaptive module. They concluded that with the daily set-up corrections using the MVCT image registration, the variation in the dose could be as large as 7.7% around the average dose. In our study, in order to minimise the errors generated during the image registrations that were done right before the treatment, automatic image registrations were performed again for each treatment fraction of each lung cancer patient. Woodford et al.Reference Woodford, Yartsev and Van Dyk5 suggested the optimal way of the MVCT registration setting for thoracic cases on helical TomoTherapy. Our study adopted their MVCT registration setting, registering the MVCT with the kVCT using either coarse or fine spacing with full image fusion technique and fine resolution selection. According to the conclusion of Woodford et al.,Reference Woodford, Yartsev and Van Dyk5 the residual errors with the applied MVCT registration setting would be small with negligible influence on the dose calculation. Deformable registration was not used in this study, but could be a valuable extension. With the deformable image registration, more accurate assignment of doses to all structures could be defined with more reliable DVHs. The biggest advantage it could endow is the ability to create adapted plans that compensate for under-dosage or over-dosage of the target or critical organs.Reference Woodford, Yartsev, Dar, Bauman and Van Dyk4 In this study, the actual adaptive plans for these lung cancer patients were not delivered. The closest one is the second-phase plan in the second-phase treatment of patient 4. The main purpose of this study is to evaluate the feasibility of using the calculation in addition to the Planned Adaptive software to predict the clinical effect of helical TomoTherapy for these lung cancer patients. Also, our study tends to demonstrate the potential clinical benefits of adaptive planning like the second-phase plan of patient 4 through the assessment. The actual protocol of doing adaptive planning using the Planned Adaptive software is out of the scope of this paper.
Another possible source of errors in the dose estimation using the Planned Adaptive software might be attributed to the merged CT image. Planned Adaptive software merged the MVCT and the kVCT images after image registration. This merged MVCT-kVCT image has the same slice spacing as the kVCT images and is created by inserting the registered MVCT in the proper location and filling the remaining slices with the kVCT (TomoTherapy Inc., Planned Adaptive Guide). For the dose calculation in the Planned Adaptive software, both kVCT and MVCT CT numbers would be converted to electron densities by using the CT number to electron density calibration curves, respectively. Although the uniformity and spatial resolutions of the MVCT images are comparable to that of the diagnostic CT images, the MVCT unit does not have the same performance characteristics as those of the diagnostic kVCT scanner. Investigations have been done on the stability of the CT number to electron density calibration curve for both MVCT and kVCT units.Reference Kilby, Sage and Rabett20,Reference Langen, Meeks, Poole, Wagner, Willoughby, Kupelian, Ruchala, Haimerl and Olivera21 Langen et al.Reference Langen, Meeks, Poole, Wagner, Willoughby, Kupelian, Ruchala, Haimerl and Olivera21 concluded that although the variation in the MVCT number is larger than that of the kVCT image, the resulting electron density difference and, thus, the dose deviation after converting electron density to dose are similar in magnitude for a 6 MV beam. Therefore, the possible dosimetric errors in the merged MVCT-kVCT images generated by the Planned Adaptive software should be within 2–3% (ref. 20,21).
CONCLUSION
The dosimetric impact and the clinical effect of the adaptive TomoTherapy have been evaluated using the Planned Adaptive software plus the . The knowledge of the tumour control and normal tissue injury over the partial or whole course of helical TomoTherapy would help the physician evaluate the necessity of adaptive planning. Based on our results for lung cancer patients treated with helical TomoTherapy, it is necessary to have the individual assessment with the in addition to the physical dose evaluated by the Planned Adaptive software. In contrast to the physical dose assessment, the better biological effect of the use of larger dose-per-fraction can only be observed in the increase of P + using the evaluation. Moreover, significant increase in the P + by 13.2% in the adaptive TomoTherapy plan in one of the lung cancer patients indicates the clinical benefit of the adaptive planning. In conclusion, the protocol of periodical evolution of the delivered dose using the may be set to achieve the end points of higher tumour control probability and/or lower normal tissue complication probability upon the adaptive radiation therapy.