Dosimetric Evaluation of Body Contour Changes to Target Volumes and Organs at Risk for Cervix and Head and Neck Radiotherapy Plans ()
1. Introduction
A radiotherapy (RT) treatment starts with the acquisition of computed tomography (CT) scan, which is used to plan an individualized treatment for the patient. Highly precise RT techniques enable delivered dose in accordance with planned dose to the clinical target volumes (CTVs) and organs and risk (OARs), based on the premise that the anatomy is unchanged since the planning stage [1]. However, it takes about 1 to 2 weeks from initial CT scan to the start of RT, and a course of intensity-modulated radiation therapy (IMRT) or volumetric modulated arc therapy (VMAT) for cervical cancer or nasopharyngeal carcinoma (NPC) is typically 5 to 7 weeks [2]. During this long time, the weight change, tumor shrinkage, and shoulder position variations possibly happen, so the delivered dose to CTVs and OARs are not the same as the planned dose owing to the body contour changes.
The abdomen of cervical cancer patients usually has a lot of adipose tissue deformed with weight fluctuation. The dosimetric effects of changes in body contour, e.g., due to weight change, have been studied for head and neck [3] and prostate cancer patients [4]. Lee et al. [5] reported the weight change during RT on the development of toxicity in patients with locally advanced cervical cancer treated with IMRT. Stauch et al. [6] and Sun et al. [7] reported dosimetric effects of the body contour changes for prostate and H&N cancer. Chow et al. [8] reported the dosimetric estimation on variations of patient size in prostate VMAT therapy. Astrid et al. [9] studied the dosimetric effects of changes in body contour for pancreatic cancer. The pieces of literature [3] - [9] showed that involuntary weight changes have a dosimetric effect on RT plans for abdominal neoplasms. Several studies show that H&N patients’ anatomy changes during the course of the treatment, and that this results in dosimetric changes from the original plans [1] [3] [7]. As a result, the patient’s body contour on the treatment day could differ from the CT scan compared with cone-beam CT (CBCT) images taken on the treatment day [10]. Figure 1 illustrated a cervical cancer patient’s pre-treatment and post-treatment image which is from
Figure 1. The body contour of pre- and post-treatment comparison.
CBCT registration, the anterior body contour shrank 1.55 cm in this example. So body contour changes induced perturbations in the dose distribution, although generally only for large changes. There were no studies found about the dosimetric evaluation of target volumes and OARs from different body contour directions for tumor sites in the abdomen and pelvis, including the cervix.
The purpose of this study is to quantify for RT of cervical cancer and NPC patients the impacts on target dose coverage and OARs sparing dose parameters as a result of changes in body contour in IMRT or VMAT plans to ensure treatment efficacy, and provide a prejudgment whether further re-assessment of the plan is needed for RT staff.
2. Materials and Methods
2.1. Patient Selection
This study included ten cervical cancer patients and ten NPC patients. A total of 20 patients were randomly selected and retrospectively analyzed according to the prescription dose based on the TNM staging. The basic information of 20 selected patients was shown in Table 1. RT planning images of the 20 patients were used for this dosimetric evaluation. Their mean age of cervical cancer patients and NPC were 46.8 ± 9.5 and 58.3 ± 12.4 years old, respectively. The patients were positioned supine and immobilized with thermoplastic fixation. For each of cervix and H&N patient, a CT image with a slice thickness of 3 mm was obtained by a CT simulator (Brilliance-16, Philips Medical Systems Inc., Cleveland, OH, USA).
2.2. Body Contour Changes
To simulate the weight’s variation during the RT, different body contours change from the original CT (o-CT) were created treatment planning system (TPS, Eclipse 13.6, Varian Medical Systems Inc., Palo Alto, CA, USA). The new CT (n-CT) was produced by o-CT through body contour shrinkage or expansion.
For cervical patients, to simulate the weight’s gain or loss, the patients’ external body contours were expanded or shrunk anteriorly, posteriorly, laterally and uniformly by 1 cm and 2 cm as shown in Figure 2. For NPC patients, to simulate the weight’s gain or loss, the patients’ external body contours were expanded or shrunk uniformly 0.5 cm and 1 cm in the facial area (from the level of pituitary fossa to 1 cm inferior of the most inferior slice of mandible), and the superior-inferior direction of shoulder expanded or shrunk 1 cm, 1.5 cm and 2 cm to simulate shoulder positional changes as shown in Figure 2(e) & Figure 2(f). The expanded tissue was assigned a CT of 0 HU (approximately equivalent to water). The body contour changed sizes for cervical cancer and NPC was chosen based on the previous works [7] [11] [12].
2.3. RT Plans
All scans were exported to Eclipse for target volumes and OARs delineation and
Table 1. Basic information of 20 selected patients.
Figure 2. Body contour change model in different directions. For cervical plans: (a) anterior, (b) posterior, (c) lateral, (d) uniform directions. For NPC plans: (e) shoulder inferior, (f) facial uniform directions.
treatment planning. The plan was delivered with 6-MV photon beams from a linear accelerator (Trilogy and Millennium 120 MLC, Varian Medical System, Palo Alto, CA, USA).
Ten cervical cancer patients had been treated with external-beam radiation therapy using IMRT (n = 5) and VMAT (n = 5). Of the 5 patients who underwent IMRT, 1 was planned to 46 Gy at 2.0 Gy/fraction, 4 was planned to 50 Gy at 2.0 Gy/fraction. Of the 5 patients who underwent VMAT, 2 were planned to 46 Gy, 3 was planned to 50 Gy at 2.0 Gy/fraction. The IMRT plans consisted of an 8-beam arrangement, with gantry angles of 0˚, 45˚, 90˚, 135˚, 179˚, 225˚, 270˚, and 315˚. The VMAT plans used two 360˚ arcs.
Ten NPC patients’ plans were also generated using IMRT (n = 5) and VMAT (n = 5). These 10 patients (p1 to p10) had already received RT with a simultaneous integrated boost to planning target volumes (PTV) including primary gross target volumes (PTV1: 66 - 70 Gy), positive lymph nodal regions (PTV2: 66 - 70 Gy), high-risk metastasis regions (PTV3: 59 - 63 Gy) and low-risk metastasis regions (PTV4: 50 - 54 Gy). The IMRT plans consisted of a 9-beam arrangement, with gantry angles of 11˚, 52˚, 93˚, 134˚, 175˚, 207˚, 248˚, 289˚, and 330˚. The VMAT plans used two 360˚ arcs.
The original plan (o-plan) in o-CT was copied to n-CT to recomputed (not reoptimized) to evaluate the effect of weight change on the dose distribution over numerous regions of interest. The recomputed plan was named a new plan (n-plan). The dosimetric parameters of PTV and OARs for body contour changes were statistically compared for IMRT vs VMAT plans, using a 2-sided Wilcoxon rank-sum test in SPSS (IBM Armonk, NY, USA). A p-value of 0.05 or less was considered statistically significant.
The dose to D95% (the dose that covers 95% of the volumes) variation of the target was defined as ΔD95% which could be expressed as follows:
(1)
where, D95%,n-plan and D95%,o-plan are the D95% of the n-plan and o-plan, respectively.
The maximal dose variation of OARs was defined as ΔDmax which could be expressed as follows:
(2)
where, Dmax,n-plan and Dmax,o-plan is the Dmax of the n-plan and o-plan, respectively.
3. Results
3.1. Target Volumes
Figure 3 shows the ΔD95% variation for ten cervical plans as body contour change in anterior (A), posterior (P), uniform (U) and lateral (L) direction. We have assigned a short name to the body contour change. + and − means contour expansion and shrinkage, respectively. For instance, A + 2 means the body contour expanded 2 cm in an anterior direction and A − 2 means body contour shrank
Figure 3. Variation in ΔD95% of PTV due to body contour change in A, P, U, L direction for cervical plans.
2cm in the anterior direction, similarly for other directions. The ΔD95% of target volumes increased near linearly as body contour shrank. The ΔD95% variations were fitted linearly. The mean slopes of the lines were −1.98% ± 0.1%, −1.21% ± 0.4%, −2.83% ± 0.5% and −1.3% ± 0.2% per cm for IMRT in A, P, U and L direction expanding, respectively. And the mean slopes of the lines were −2.16% ± 0.1%, −1.25% ± 0.6%, −2.87% ± 0.3% and −1.18% ± 0.3% per cm for VMAT. The difference between IMRT and VMAT was significant (p < 0.05) in A, U, L directions. For both IMRT and VMAT plans, the anterior region ΔD95% showed larger variations than do in the posterior region, this was mainly because the anterior body contour changed length is longer than the posterior body contour change. In A and U direction, the ΔD95% for VMAT was overall larger than that for IMRT, which is mainly because in VMAT plans the dose is delivered by two full arcs with multileaf collimator modulation and the IMRT plans the dose is delivered in an 8-beam arrangement. Whereas in L direction, the ΔD95% shows slightly larger than that for VMAT, probably because in IMRT plans, the gantry angle of 45˚, 90˚ and 135˚ had a greater influence than VMAT plans during the lateral body contour change. The difference between IMRT and VMAT was not significant (p > 0.05) in the P direction.
The results of ΔD95% variation as body contour change in superior(S), inferior (I), facial (F) direction for ten H&N plans are shown in Figure 4. The ΔD95% of PTV1 decreased near linearly as body contour expanded in F direction, and the body contour change of S, I directions almost did not influence PTV1. The mean slopes of the PTV1 ΔD95% lines were −3.88% ± 0.18% per cm for IMRT and −4.13% ± 0.25% per cm for VMAT in F direction expanding (P < 0.05). The ΔD95% of PTV2 has a similar influence to PTV1 whereas the PTV2 of p5, p8 covered to supraclavicular nodes, so the PTV2 of p5 and p8 was influenced by S and I directions. The mean slopes of the PTV2 ΔD95% lines were −3.41% ± 0.31% per cm for IMRT and −4.08% ± 0.22% per cm for VMAT in F direction expanding. The ΔD95% for VMAT was overall larger than that for IMRT (P < 0.05), which is also mainly because in VMAT plans the dose is delivered by two full arcs with multileaf collimator modulation and the IMRT plans the dose is delivered in a 9-beam arrangement. Whether the body contour change affects the PTV3 or not depends on the position of PTV3. As shown in Figure 4 PTV3, the ΔD95% loss was seen in the C7-T2 region in p1 and p2, and the ΔD95% loss was seen in the
p4, p6: without PTV2; p2, p3, p5, p9: without PTV4.
Figure 4. Variation in ΔD95% of PTV due to body contour change in S, I, F direction for NPC plans.
C5-C6 region in p6, p7, p8, so the shoulder body contour expansion had larger influence to PTV3, particularly those with the lower neck region. The shoulder body contour shrinkage had a relatively smaller influence on PTV3. The facial body contour expansion led to ΔD95% loss in the 10 plans, but the facial body contour shrinkage had both ΔD95% loss and gain. For both cervix and NPC patients, p1, p2 ∙∙∙ p5 are IMRT plans, and p6, p7 ∙∙∙ p10 are VMAT plans. For both IMRT and VMAT plans, shoulder body contour change in the superior-inferior directions resulted in the greatest loss of coverage for PTV4 and up to 4.45%/cm (p8), because the PTV4 is lower neck target (from C6 to T2 vertebrae). The facial body contour change did not show a large loss in PTV4 because PTV4 was far enough from the facial region.
3.2. OARs Sparing
In addition to decreasing the dose to the target volumes, body contour change also has the potential to increase the dose to OARs. The change in dose to the bladder, rectum, small bowel for cervical patients and brain stem and the spinal cord for NPC patients in each direction change was evaluated in Eclipse for both IMRT and VMAT plans.
For the ten selected cervical patients, the rectum, bladder and small bowel were investigated from the dose-volume histogram. Figures 5(a)-(d) show the ΔD50% changes for rectum in A, P, U, L directions, respectively. The ΔD50% were 1.86% ± 0.16%, 0.96% ± 0.31%, 2.83% ± 0.31%, 1.32% ± 0.23% per cm for IMRT and 1.49% ± 0.21%, 1.86% ± 0.59%, 2.40% ± 0.24%, 1.25% ± 0.16% per cm for VMAT in A, P, U and L directions shrinkage, respectively. The difference between IMRT and VMAT was significant (p < 0.05) in all directions. Only in the P direction, the ΔD50% for VMAT was larger than that for IMRT, which is mainly because the rectum is close to the posterior body contour and the range of beam path length affected by body contour change and the relative portion of the dose delivered in the corresponding direction. Figures 6(a)-(d) show the ΔD50% changes for bladder in A, P, U, L directions, respectively. The ΔD50% were 2.06% ± 0.48%, 1.24% ± 0.53%, 2.32% ± 0.41%, 1.12% ± 0.21% per cm for IMRT and 2.41% ± 0.29%, 1.14% ± 0.24%, 2.98% ± 0.32%, 1.31% ± 0.16% per cm for VMAT in A, P, U and L direction shrinkage, respectively(P < 0.05). The results of ΔD50% for bladder were opposite to the ΔD50% for rectum for IMRT and VMAT in A, P, U, L directions, which is mainly because the engorged bladder is enough larger than the rectum in volume which involved the target area. The dose is delivered when the gantry was moving for VMAT plans. So the bladder ΔD50% for VMAT shows overall larger than it in the rectum for IMRT except for P direction. Figure 7 shows the small bowel ΔDmax variation. The body contour expansion leads to Dmax decrease and shrinkage leads to dose increase for small bowel. The Dmax of small bowel increased by up to 391 cGy, 163 cGy, 439 cGy, 276 cGy for IMRT and 295 cGy, 378 cGy, 348 cGy, 188 cGy for VMAT in A-2, P-2, U-2, and L-2 direction, respectively.
Figure 5. Variation in ΔD50% of the rectum due to body contour change in A, P, U, L directions for cervical plans.
Figure 6. Variation in ΔD50% of the bladder due to body contour change in A, P, U, L directions for cervical plans.
For the ten selected NPC patients, the Dmax change to the brain stem and spinal cord was displayed in Figure 8. The ΔDmax of the brain stem of all of the examined dosimetric cases was up to 195 cGy for IMRT and 210 cGy for VMAT in F-1 direction, and the ΔDmax of the spinal cord increased by up to 209 cGy for IMRT and 628 cGy for VMAT in S-2 direction.
Figure 7. Variation in ΔDmax of small bowel due to body contour change in A, P, U, L directions for cervical plans.
Figure 8. Variation in ΔDmax of the brain stem and spinal cord due to body contour change in S, I, F directions for NPC plans.
4. Discussion
Accurate dose delivery to target volumes and organs at risk (OARs) is essential to ensure the radiotherapeutic effect and minimize the risk of normal tissue toxicity, whereas weight fluctuation frequently occurs during the whole radiotherapy course, which may cause the body contour changes to effect on treatment accuracy. Therefore, we evaluated the impacts of body contour changes to target volumes and OARs in the radiotherapy plans of cervical cancer and head and neck tumors. A course of RT is typically several weeks. During the long period from the CT scans to RT finish, the weight change for abdominal and H&N patients was reported in previous literature [13] [14] [15] [16]. There is evidence showing that weight change is correlated with external contour changes [17].
Booth et al. [18] reported that 68% of the 198 analyzed CBCT images from 19 prostate patients were in the range of 0 - 1 cm, 28% 1 - 2 cm, and 4% > 2 cm with deviations occurring mostly in the postero-lateral direction. Chow et al. [13] studied body contour shrunk by 0.5, 1, 1.5, 2 cm in anterior, left and right direction for five patients' IMRT and VMAT prostate plans. Sun et al. [7] reported that prostate patients who have body contour change less than 2 cm at a single side or less than 1 cm uniformly are unlikely to need further assessment. For H&N patients, Chen et al. [14] reported 25 NPC patients shrank the external contours with different margins (2, 3, and 5 mm). Our study illustrated the detailed dosimetry of PTV and OARs for cervical cancer and NPC in different body contour changes based on previous works [4] [7].
Currently, IMRT or VMAT are widely used for planning to treat patients with cervical cancer and NPC, which the dosimetric effect of body contour change in different directions is still unknown. Our contour-size change model is consistent with that used in a recent study by Sun et al. [7], in which the contour-size effect was assessed for VMAT only. Also, Sun et al. only evaluated the dose percentage and isodose line shift whereas our works focus on the dose variations of D95% for whole PTV and D50% or Dmax for OARs. Another discrepancy between our means with Sun et al. is that the different PTVs with simultaneous integrated boost for H&N patients were discussed in this paper. This study showed that the body contour shrinkage (expansion) caused the dose increase (decrease) to PTV or OARs and this finding matches the results from the study published by Pair et al. [4]. This paper showed the cervical patients’ D95% of target volumes decreased by up to 1.98%, 1.21%, 2.87%, 1.3% per cm for IMRT and 2.16%, 1.38%, 2.83% and 1.18% for VMAT in A, P, U and L direction expanding, respectively. Pair et al. [4] reported that the prostate patients' target mean dose decreased by 2.9% per cm for IMRT and by 3.6% for VMAT in U direction. The result of IMRT plans is larger than do VMAT plans in our research which was mainly because the IMRT Gantry angle arrangement is different. In Pair’s [4] study, the change of D10% for rectum and bladder were 2.8% to 3.5% which were close to our study that the D50% for rectum and bladder were increased by up to 2.98%.
For H&N plans, we studied the body contour changes in facial area and shoulder based on the previous study by Sun et al. [7] and Neubauer et al. [19]. Neubauer et al. [19] examined ten patients and 243 CTs, and found that 2% of shoulder shifts were greater than 1 cm. Noble et al. [1] measured lateral neck diameter which is 175 mm on the first day and 162 mm on the final treatment day. In our study, we chose shoulder changed 1 cm and 2 cm in S-I directions and 0.5 cm and 1 cm in the facial area. Chen et al. [14] shrank external contour with different margins (2, 3, and 5 mm) and reported that the D95% of PTV1 was increased by 1.9% to 2.9%, which was similar to our result 1.95% ± 0.58%/5mm. The PTV3/4 in Figure 4 shows relative irregular changes for D95% compared with PTV1/2. The PTV3/4 extended outside of the shrunk body contour, so the volumes within the new contour were affected by the build-up effect. Zhao et al. [20] found an increase in the maximum dose to the spinal cord and brainstem volumes of 560 cGy and 250 cGy, respectively, by comparing repeat CT imaging to the dose distribution on the original planning CT. Wang et al. [21] reported the NPC repeated CT scan after 18 fractions, the mean volume of the left and right parotid decreased 6.19 mL and 6.44 mL, and the center of C2 vertebral body slices contracted with the mean contraction of 8.2 mm, 9.4 mm, and 7.6 mm while the maximum dose to the brain stem and spinal cord increased by 0.08 to 6.51 Gy and 0.05 to 7.8 Gy. Our study showed that the ΔDmax of the brain stem of all of the examined dosimetric cases was up to 195 cGy for IMRT and 210 cGy for VMAT in F-1 direction, and the ΔDmax of the spinal cord increased by up to 209 cGy for IMRT and 628 cGy for VMAT in S-2 direction. The dosimetrists need to take the dosimetric changes into account during the RT plan design.
The factors which affect the dose to target volumes and OARs were anatomy and setup error, this paper studied dose discrepancy which resulted from the anatomy change. The limitation of this study was that the location, geometry, and size of the tumor, OARs may change during the RT period. Moreover, the patients' body contour change may not be as regular in real clinical status. These situations were not discussed in this study. The body contour changed methods in this work have been widely reported in previous literature [4] [7] [8] [13] [14]. It is a simplified way but good choice to predict such complicated dosimetric problems. Therefore, the RT staff could make a preliminary judgment of dosimetric parameters induced by body contour changes for cervical cancer and NPC patients based on the findings throughout this work.
5. Conclusion
The dosimetric evaluation of body contour changes to PTV and OARs for cervical cancer and NPC plans was studied. The body contour shrinkage or expansion affects the IMRT and VMAT dose delivery. The body contour changes may impact the dosimetry of the PTV and OARs to a different extent, depending on the directions and magnitude of the body contour changes. The RT staff could determine whether resimulation and replanning or not according to which body contour directions were changed.
Acknowledgements
This work was supported by National Natural Science Foundation of China (Grant No. 11575221), National Key Research and Development Program of China (2018YFA0404300).
NOTES
*Corresponding authors.