 Research
 Open Access
 Published:
Analysis of interfraction setup errors and organ motion by daily kilovoltage cone beam computed tomography in intensity modulated radiotherapy of prostate cancer
Radiation Oncology volume 7, Article number: 56 (2012)
Abstract
Background
Intensitymodulated radiotherapy (IMRT) enables a better conformality to the target while sparing the surrounding normal tissues and potentially allows to increase the dose to the target, if this is precisely and accurately determined. The goal of this work is to determine interfraction setup errors and prostate motion in IMRT for localized prostate cancer, guided by daily kilovoltage cone beam computed tomography (kVCBCT).
Methods
Systematic and random components of the shifts were retrospectively evaluated by comparing two matching modalities (automatic bone and manual softtissue) between each of the 641 daily kVCBCTs (18 patients) and the planning kVCT. A simulated Adaptive Radiation Therapy (ART) protocol using the average of the first 5 kVCBCTs was tested by nonparametric bootstrapping procedure.
Results
Shifts were < 1 mm in leftright (LR) and in superoinferior (SI) direction. In anteroposterior (AP) direction systematic prostate motion (2.7 ± 0.7 mm) gave the major contribution to the variability of results; the averages of the absolute total shifts were significantly larger in anterior (6.3 ± 0.2 mm) than in posterior (3.9 mm ± 0.2 mm) direction. The ART protocol would reduce margins in LR, SI and anterior but not in posterior direction.
Conclusions
The online softtissue correction based on daily kVCBCT during IMRT of prostate cancer is fast and efficient. The large random movements of prostate respect to bony anatomy, especially in the AP direction, where anisotropic margins are needed, suggest that daily kVCBCT is at the present time preferable for high dose and high gradients IMRT prostate treatments.
Background
Intensitymodulated radiation therapy (IMRT) for prostate cancer enables creating a steep dose gradient between prostate and rectum, allowing in principle higher doses to the target and high cure rates while reducing late rectal toxicity [1, 2]. Different filling conditions of bladder and rectum can significantly influence the interfraction position of the prostate during IMRT [3, 4], with consequent modifications of dose distribution in the target and adjacent organs. This may impair local control, with increased risk of late sequelae. To reduce the extra margin needed to allow for prostate motion, an accurate localization of the prostate position at the time of treatment is therefore needed.
Kilovoltage cone beam computed tomography (kVCBCT) is one method to assess and correct for interfraction prostate localization immediately before treatment [5]. It enables direct visualization of softtissue targets and organs at risk, and with a flatpanel imager may combine volumetric and radiographic/fluoroscopic imaging using the same device. However, daily imageguidance affects the traditional treatment workflow and can increase the patient's radiation dose [6]. Several studies have been published on direct visualization of the prostate observed with different modalities of image guidance, including implanted fiducial markers [7], ultrasound [8], kV or MV computed tomography [9–11] and electromagnetic responders [12]. However, studies are scarce on both setup and organ motion using daily kVCBCT.
In our current clinical practice at Bellaria Hospital of Bologna, daily target localization by kVCBCT (CBCT in the following) has been routinely performed since November 2008. The aim of this study is to determine the interfraction setup error and the prostate motion relative to the bony anatomy assessed with CBCT in patients treated for prostate cancer with radical radiotherapy. Results of daily corrections of 18 patients are presented. Moreover, an offline Adaptive Radiation Therapy (ART) imageguided correction protocol is retrospectively simulated to determine whether the number of images performed could be reduced.
Methods
Patients
Eighteen prostate cancer patients treated between January and September 2009 with IMRT were retrospectively studied. Patients were staged according to the National Comprehensive Cancer Network risk class assessment [13]. Each of these patients underwent a 3mm slice spacing planning CT (PCT) scan on a LightSpeed CT scanner (General Electric, UK). All patients were instructed to empty the rectum and fill the bladder, drinking 500 ml of water 30 min before the PCT and therapy. Planning target volumes (PTVs) were obtained by using threedimensional (3D) automatic expansions of clinical target volumes (CTVs), applying 8mm in leftright (LR) and anterior (A), 5mm in posterior (P) and 8mm in superoinferior (SI) directions. Intermediate and highrisk patients were prescribed 76 Gy to the CTV1 (the prostate and the first proximal third of the seminal vesicles) and 50 Gy to the CTV2 (the rest of seminal vesicles) in 38 fractions with Simultaneous Integrated Boost. Lowrisk patients were prescribed 74 Gy in 37 fractions soley to the prostate. Written informed consent was obtained from the patients for publication of this report.
Imageguidance procedure and treatment
Prostate position was assessed before each fraction through CBCT image guidance (XVI, Elekta, Crawley, UK). Patients were positioned first by using lasers and skin marks, then a CBCT was acquired. An average of 36 CBCT (2938) per patient were obtained and used in this analysis (641 CBCTs).
To quantify the setup error we first performed a pelvic bones match ("Bmatch") of the CBCT on the PCT scan based on a fully automated 3D chamfer algorithm [14]. Then, the matching was manually adjusted by a RT (technologist) to overlap the prostate on the PCT and CBCT scans through a greyvalue based softtissue matching ("Tmatch"), assessing the total interfraction error (setup + organ motion). Registrations were based on a rigidbody approach. The sodetermined corrections were automatically applied to a robotic table with 6 degrees of freedom (Hexapod Evo, ElektaMedical Intelligence, Crawley, UK) before the patient treatment (Linac Synergy S, Elekta, Crawley, UK). Values of interfraction setup and total positioning displacements were registered for the three principal axes, in leftright (LR) (X), superoinferior (SI) (Y), and anteroposterior (AP) (Z), and for the three rotation angles (pitch, roll and yaw). The setup translation components X_{B}, Y_{B}, Z_{B} refer to the Bmatch, and the total positioning translation components X_{T}, Y_{T}, Z_{T} refer to the Tmatch, respectively.
Motion of the prostate relative to the bony anatomy was defined as the difference between Tmatch and Bmatch. Therefore, the prostate interfraction movement vector relative to bony anatomy will be identified by the differences X_{T}X_{B}, Y_{T}Y_{B}, Z_{T}Z_{B}. Positive values for X, Y and Z shifts indicate a left, superior and anterior displacements of the isocenter.
Statistical analysis
Statistical analysis was performed using Statistical Package STATA 9.0 (StataCorp LP). For each patient, the average deviation μ_{i} and the standard deviation σ_{i} for Bmatch, Tmatch and for organ motion have been calculated in each direction. For the entire population of patients, the mean value of μ_{i} values (M, group systematic error), the root mean square of the patient σ_{i} (σ, random error), and the standard deviation of the patient means (Σ) were computed. The reported error for the means is one standard error. Normality test of each single variable distribution was carried out using ShapiroWilk test [15]. In order to test whether patients behaved differently for X, Y, and Z shifts, different statistics were carried out. Using an Ftest, Analysis of Variance (ANOVA) was performed to verify that the patient means were statistically different. As ANOVA is based on the strong assumption of normality of the underlining random variable, for sake of robustness a nonparametric test (KruskalWallis) was performed to compare medians of each patient (pvalue was determined using a χ^{2} distribution). The Levene test was carried out to test the equality of variances of each patient through an Ftest. All tests were twosided.
Simulation of an offline protocol: bootstrap analysis of the first 5 CBCT
To test the goodness of a first fivedays average estimation for correction of each patient instead of a more costly daily online procedure, an ART offline protocol was simulated and compared with a nocorrection protocol based only on skin marks alignment and with our daily online correction protocol. Displacements on X, Y, and Z were defined as the differences between daily correction shifts (total positioning translation components) and the average of the first five days patient's shifts. Those values were taken in absolute terms and cumulative distributions of errors were built. Cutoff margins ranging from 1 mm to 8 mm were considered, and the percentages of the observations outside each cutoff level were computed. For each of those cutoffs, confidence intervals were obtained using a nonparametric bootstrap procedure as in Efron, 1979 [16]. Nonparametric bootstrapping is a statistical method that can be implemented by constructing a large number of simulated samples from the original dataset, each of which is obtained by random sampling empirical observations from the original dataset. As a final results, instead of having a single sample, one is able to get a large number of samples, each having the same statistical properties of the previous one as they are derived from it. Meanwhile, the great advantage of bootstrap is to derive estimates of standard errors and confidence intervals for complex estimators of parameters of the distribution, how in the case of the present study.
Interobserver variability study
Tmatch method implies a manual intervention of the operator. To evaluate the interobserver variability, four RTs usually performing the CBCT/PCT matching were asked to match 9 randomly chosen CBCT acquisitions of 3 patients included in this study.
Results
Evaluation of systematic and random errors
The systematic rotations of the bones, as well as the rotations of the organ relative to the bones around the three axes were negligible (less than 1.2 degree). For that reason, further analysis has been performed on translations only.
Interobserver variability in Tmatch method was found to be within 0.9 mm (1 SD) in leftright (LR) (X) direction, and within 1.5 mm and 1.8 mm (1 SD) in superoinferior (SI) (Y) and anteroposterior (AP) (Z), respectively. Mean values μ_{i} and standard deviations σ_{i} of total positioning, setup and prostate motion translational shifts registered for all treatment fractions in each direction for the 18 patients are reported in Table 1 along with systematic (M) and random errors (σ) for each variable, and standard deviation of the μ_{i} values (Σ). In X and Y directions, M values are both on the order of 1 mm or less. The σ and Σ values in X_{T} and X_{B} are very similar, indicating a major contribution from bone misalignment, while systematic and random errors calculated for internal organ motion seem to be small in this direction. In Y direction, the setup and the internal organ motion components seem to contribute equally to the total σ and Σ values. It is worth noting the larger values for M, σ and Σ for the total shift in the Z direction (Z_{T}) and for the internal organ motion (Z_{T}Z_{B}), in comparison with those of the X and Y directions. The basically larger value of M for Z_{T} than for Z_{B} seems to indicate a significant systematic internal organ motion relative to the bone structure in that direction. Also values of σ and Σ are higher than for X and Y, with a slightly higher contribution from setup errors to the total positioning errors than from internal organ motion. These data are in agreement with setup errors previously reported for prostate cancer patients [10, 11], except for the higher interfraction prostate motion along the Z direction. Parametric and nonparametric tests show (Table 1) that all the differences are significant. The patients have significantly different means (F test, pvalue < 1%), medians (KW, pvalue < 1%) and standard deviations (Levene, pvalue < 5%) for all variables.
Figure 1 (a,b,c) shows the histograms of the shifts in the three directions, relative to the PCT scan, registered in the two modalities of matching (X_{B}, X_{T}; Y_{B}, Y_{T}; Z_{B}, Z_{T}) and the scatterplots of the couples of corresponding shifts, for the aggregate set of data (all patients, all fractions). The histograms allow us a rough evaluation of the shift distributions and to compare the two matching modalities. The null hypothesis of normality of the distribution is rejected in all cases (pvalue < 1%), due to the high values of the kurtosis, in turn due to the substantial tails on both sides of the distributions. On the other hand, all the distributions are symmetric (skewness close to 0), allowing one to consider the mean and the standard deviation of each set of the data as a good parameter to represent each distribution. The correlation matrix computed on all data between the three translation shifts acquired in the two matching modalities shows that the highest correlation coefficients are observed for the couples of shifts along the same axis (Figure 1).
Absolute shifts and evaluation of margins
The absolute values of the shifts determined by the Tmatch (X_{T}, Y_{T}, Z_{T}) allow us to evaluate the percentage of fractions for which the absolute values would exceed the margin by different protocols. Figure 2 shows the results for X, Y, Z.
In X and Y directions the percentage of treatments that would be outside the margins of 8 mm is less than 5%. By assuming, for example, a margin of 5 mm, about 20% of the treatments would result outside the margin along X and Y axes, in both positive and negative directions.
Along Z two cases have been considered: Z_{T} > 0 and Z_{T} < 0 (anterior direction, 446 observations, mean of the absolute values 6.3 ± 0.2 mm; posterior direction, 192 observations mean of the absolute values 3.9 ± 0.2 mm; statistically different p < < 1%). Only for three observations the value was zero and were not considered. With a margin of 8 mm in the anterior direction, about 30% of the treatments would result outside the margin. With margins of 5 mm this percentage would grow to about 55%. In the posterior the choice of 5 mm margin would have left outside about 30% of the treatments. Even an hypothetical margin of 8 mm (too deep in rectum to be applicable) would still leave about 10% of the treatments outside the margin.
Simulation of an offline protocol: bootstrap analysis of the first 5 CBCTs
The solid lines in Figure 2 simulate the results of an offline protocol with 5 CBCTs at the beginning of the treatment. In X and Y, by assuming for example, 5 mm margin, the percentage of fractions outside the margins would be reduced from 20% to less than 10% if the mean of the first 5 CBCTs was assumed for each patient and the corresponding average shifts applied for the following fractions. Along the positive direction of Z (Z_{T} > 0), by assuming a margin of 8 mm, the percentage would be reduced from 30% to about 5%. In the posterior direction, a 5 mm margin still would leave about 30% of the treatments outside. There would be no improvement in the Z negative direction by using the average of the first 5 CBCTs.
Discussion
The extensive statistical analysis of the corrections to make to the position of the patient in the imageguided radiation therapy (IGRT) of prostate cancer after daily kilovoltage CBCT allowed us to get information on the correctness of our protocol, on the organ movement relative to bony anatomy, on the large variability of patients and on the use of an ART protocol.
The systematic and random errors of setup and organ motion along LR (X) and SI (Y) for each patient are of the order of a few mm. This demonstrates the correctness of our positioning procedure. Both parametric and nonparametric statistical tests show that patients behave differently from each others.
In AP (Z) direction, it is the organ motion that gives the major contribution to the variability, both intra and interpatient. This axis is the more critical, probably due to the variability in patient's rectal distension, although we attempted to treat each patient with an empty rectum. Our findings are in agreement with other reported studies showing that variable rectal filling is the major cause of movement of prostate and seminal vesicles [3, 4].
Table 1 confirms that along X axis, the major contribution to the total positioning error is due to bone misalignment, while along Y and Z bone and organ motion shifts are of the same order, with a larger organ motion for Z than for X and Y.
The statistics of the absolute values of Z_{T} for positive and negative values confirm the utility of different margins in anteroposterior directions, with the anterior one larger than the posterior one. This is probably one of the major findings of this study.
More than 95% of the prostate displacements in X and Y direction would be within the 8 mmmargins currently used in our center, even in absence of any correction, as shown in Figure 2 (dasheddotted curves).
Conversely, in Z direction, only 70% of the displacements were within the currently used margins. As a consequence, our findings seem to indicate that, despite setting a daily patient alignment on skin marks, large displacements could occur. Facilities that use this traditional protocol without imageguidance should address strategies to reduce the internal organ motion such as regularly emptying the rectum [11, 17] or applying other modalities [18] before the treatment in order to avoid missing the target and/or increased toxicity to organs at risk.
Several ART protocols have been tested to try to reduce the costly daily image guidance procedure, limiting the systematic error Σ but not the random error, σ [19–21].
In this study we compared the percentages of treatments that would exceed the margins if no corrections are made with an ART offline protocol using the average of the first 5 CBCTs as an offset for all the following fractions of each patient. In X and Y direction, by assuming for example, 5 mm margin, the percentage of fractions outside the margins would be reduced from 20% to less then 10%. Moreover, the curves show that, by applying this offline protocol, a 6 mmmargin would probably be adequate to compensate for the majority of interfraction total displacements.
Again, the Z direction is the most critical. The dasheddotted lined is outside the area that represent the region contained inside the 5% and 95% confidence levels for X, Y and positiveZ directions. That means that the improvement from 30% to 5% that we would get with this offline correction protocol is significant, and not obtained by chance. For the negativeZ direction both solid and dotteddashed lines are contained in the 5%95% confidence level area. That means that they differ only by chance, and we can consider that there is no improvement in the Z negative direction by using the average of the first 5 CBCTs. This is another major finding of this study.
In the near future, at our centre, an hypofractionated scheme will be introduced for IMRT prostate cancer treatments. As seen in this report, setup and organ motion errors may be relevant, and, if not corrected, they could have a significant detrimental impact on the treatment. For this reason, the actual online daily correction protocol has been maintained for the near future.
The limitations of this study include the fact that uncertainties such as intrafraction organ motion, contouring delineation of CTVs, intraobserver variability and the intrinsic error caused by the rigid shift approximation were not considered. Further studies will address these issues at our institution.
Conclusions
The online softtissue imageguidance correction based on daily kilovoltage CBCT during intensitymodulated radiotherapy of prostate cancer patients is a fast and efficient procedure. Because of the large random errors that characterize the prostate movement with respect to the bony anatomy, especially in the anteroposterior (Z) direction, we showed that an offline protocol based on a few CBCTs could improve the treatment margins in leftright (X), superoinferior (Y) and positive Z direction, but not in negative Z direction with respect to no correction. At present it seems that the daily CBCT is preferable for high dose and high gradients intensitymodulated radiotherapy treatments. Further studies are needed in order to optimize the number of first CBCT scans to be used for an offline ART procedure.
Abbreviations
 IMRT:

Intensitymodulated radiotherapy
 kVCBCT:

CBCT: Kilovoltage cone beam computed tomography
 ART:

Adaptive radiation therapy
 LR:

X: Leftright
 SI:

Y: Superoinferior
 AP:

Z: Anteroposterior
 PCT:

Planning CT
 CTV:

Clinical target volume
 Tmatch:

Greyvalue based softtissue matching
 Bmatch:

Pelvic bones match
 X_{B}:

Y_{B}, Z_{B}: Translation components by the Bmatch
 X_{T}:

Y_{T}, Z_{T}: Translation components by the Tmatch.
References
 1.
Zelefsky MJ, Chan H, Hunt M, Yamada Y, Shippy AM, Amols H: Longterm outcome of high dose intensity modulated radiation therapy for patients with clinically localized prostate cancer. J Urol 2006, 176: 14151419. 10.1016/j.juro.2006.06.002
 2.
JereczekFossa BA, Orecchia R: Evidencebased radiation oncology: definitive, adjuvant and salvage radiotherapy for nonmetastatic prostate cancer. Radiother Oncol 2007, 84: 194215.
 3.
Guckenberger M, Pohl F, Baier K, Meyer J, Vordermark D, Flentje M: Adverse effect of a distended rectum in intensitymodulated radiotherapy (IMRT) treatment planning of prostate cancer. Radiother Oncol 2006, 79: 5963. 10.1016/j.radonc.2006.03.004
 4.
Pinkawa M, Siluschek J, Gagel B, Demirel C, Asadpour B, Holy R, Eble MJ: Influence of the initial rectal distension on posterior margins in primary and postoperative radiotherapy for prostate cancer. Radiother Oncol 2006, 81: 284290. 10.1016/j.radonc.2006.10.028
 5.
Jaffray DA, Siewerdsen JH: Conebeam computed tomography with flatpanel imager: initial performance characterization. Med Phys 2000, 27: 13111323. 10.1118/1.599009
 6.
Ding GX, Duggan DM, Coffey CW: Accurate patient dosimetry of kilovoltage conebeam CT in radiation therapy. Med Phys 2008, 35: 11351144. 10.1118/1.2839096
 7.
Crook JM, Raymond Y, Salhani D, Yang H, Esche B: Prostate motion during standard radiotherapy as assessed by fiducial markers. Radiother Oncol 1995, 37: 3542. 10.1016/01678140(95)01613L
 8.
Reddy NMS, Nori D, Sartin W, Maiorano S, Modena J, Mazur A, Osian A, Sood B, Ravi A, Sampath S, Lange CS: Influence of volumes of prostate, rectum, and bladder on treatment planning CT on interfraction prostate shifts during ultrasound imageguided IMRT. Med Phys 2009, 36: 56045611. 10.1118/1.3260840
 9.
Perks J, Turnbull H, Liu T, Purdy J, Valicenti R: Vector analysis of prostate patient setup with imageguided radiation therapy via kV cone beam computed tomography. Int J Radiat Oncol Biol Phys 2011, 79: 915919. 10.1016/j.ijrobp.2010.04.008
 10.
Bylund KC, Bayouth JE, Smith MC, Hass AC, Bhatia SK, Buatti JM: Analysis of interfraction prostate motion using megavoltage cone beam computed tomography. Int J Radiat Oncol Biol Phys 2008, 72: 949956. 10.1016/j.ijrobp.2008.07.002
 11.
Fiorino C, Di Muzio N, Broggi S, Cozzarini C, Maggiulli E, Alongi F, Valdagni R, Fazio F, Calandrino R: Evidence of limited motion of the prostate by carefully empting the rectum as assessed by daily MVCT image guidance with helical tomotherapy. Int J Radiat Oncol Biol Phys 2008, 71: 611617. 10.1016/j.ijrobp.2008.01.048
 12.
Shah AP, Kupelian PA, Willoughby TR, Langen KM, Meeks SL: An evaluation of intrafraction motion of the prostate in the prone and supine positions using electromagnetic tracking. Radiother Oncol 2011, 99: 3743. 10.1016/j.radonc.2011.02.012
 13.
Crawford ED: Use of algorithms as determinants for individual patient decision making: National Comprehensive Cancer Network versus artificial neural networks. Urology 2003, 62: 1319. 10.1016/j.urology.2003.10.008
 14.
Borgefors G: Hierarchical chamfer matching: a parametric edge matching algorithm. IEEE Transactions Pattern Analysis and Machine Intelligency 1998,10(6):849865.
 15.
Shapiro SS, Wilk MB: An analysis of variance test for normality (complete samples). Biometrika 1965,52(34):591611. 10.1093/biomet/52.34.591
 16.
Efron B: Bootstrap methods: another look to the Jackknife. Ann Statist 1979, 7: 126. 10.1214/aos/1176344552
 17.
Stasi M, Munoz F, Fiorino C, Pasquino M, Baiotto B, Marini P, Malinverni G, Valdagni R, Gabriele P: Emptying the rectum before treatment delivery limits the variations of rectal dosevolume parameters during 3DCRT of prostate cancer. Radiother Oncol 2006, 80: 363370. 10.1016/j.radonc.2006.08.007
 18.
Wang CW, Chong FC, Lai MK, Pu YS, Wu JK, Cheng JC: Setup errors due to endorectal balloon positioning in intensity modulated radiation therapy for prostate cancer. Radiother Oncol 2007, 84: 177184. 10.1016/j.radonc.2007.06.009
 19.
Hoogeman MS, van Herk M, de Bois J, Lebesque JV: Strategies to reduce the systematic error due to tumor and rectum motion in radiotherapy of prostate cancer. Radiother Oncol 2005, 74: 177185. 10.1016/j.radonc.2004.12.010
 20.
Baum C, Birkner M, Alber M, Paulsen F, Nüsslin F: Dosimetric consequences of the application of offline setup error correction protocols and a hullvolume definition strategy for intensity modulated radiotherapy of prostate cancer. Radiother Oncol 2005, 76: 3542.
 21.
Nijkamp J, Pos FJ, Nuver TT, de Jong R, Remeijer P, Sonke JJ, Lebesque JV: Adaptive radiotherapy for prostate cancer using kilovoltage conebeam computed tomography: first clinical results. Int J Radiat Oncol Biol Phys 2008, 70: 7582. 10.1016/j.ijrobp.2007.05.046
Acknowledgements
The authors would like to thank Mr. Cristiano Cavicchi, Mr. Mattia Baratozzi and Mr. Dario Capanni of Elekta Company for technical support in data collection.
Author information
Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
MP participated in the design of the study, carried out the treatment plannings, participated in data collection and interpretation, participated in drafting and final revising of the manuscript; SM performed the statistical analysis, including the bootstrap procedure for the ART protocol, participated in data interpretation, in drafting and final revising of the manuscript; PF participated in data interpretation and in drafting and final revising of the manuscript; CC participated in data collection and interpretation, and in drafting the manuscript; CDE participated in the design of the study and in data collection; GPF conceived the study and participated in its design and coordination. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
About this article
Cite this article
Palombarini, M., Mengoli, S., Fantazzini, P. et al. Analysis of interfraction setup errors and organ motion by daily kilovoltage cone beam computed tomography in intensity modulated radiotherapy of prostate cancer. Radiat Oncol 7, 56 (2012). https://doi.org/10.1186/1748717X756
Received:
Accepted:
Published:
Keywords
 Prostate
 Interfraction
 ConeBeam CT
 Imageguided radiotherapy
 Organ motion
 Adaptive radiation therapy
 Bootstrap statistical analysis