Abstract
Background/Aim: This study pursued two goals: Firstly, to search for anatomical structures strongly correlating with dose deterioration, and secondly to investigate the effectiveness of image registration focusing on critical anatomy by comparing it with a conventional method. The aim was to achieve robust image registration to correct for anatomical changes during treatment. Patients and Methods: Twenty patients with head and neck cancer were enrolled, and 68 simulation computed tomography (CT) and rescan CT image sets were retrospectively analyzed. Forty volumetric-modulated arc therapy and intensity-modulated proton therapy plans were generated and recalculated according to the rescan CT to evaluate the dose effects of anatomical changes. Correlation coefficients were calculated for the relationships between the six-axis motion of the anatomy and the dose indices for the clinical target volume (CTV) and organs at risk. In the image registration, we compared a conventional method and target-based registration that limited the registration range to the CTV and vertebrae. Results: The CTV coverage and spinal cord dose were correlated with the position error associated with the pitch and vertical position of the vertebrae, and the parotid gland and oral cavity dose were strongly correlated with the position error associated with the roll of the clivus and mandible. The target registration improved CTV coverage and suppressed the increase in dose to organs at risk compared with conventional methods. Conclusion: Monitoring vertebral alignment, the assessment and correction of positioning errors associated with the clivus and mandible position errors are important to ensure the quality of daily treatment. Target-based registration may allow for more robust image registration.
- Head and neck cancer
- image guided radiotherapy (IGRT)
- Volumetric-modulated arc therapy (VMAT)
- intensity-modulated proton therapy (IMPT)
- image registration
- anatomical change
- positioning error
- patient positioning
In the case of head and neck cancer, intensity-modulated proton therapy (IMPT) can achieve better dose reduction to normal tissues than intensity-modulated radiotherapy (IMRT) using X-ray beams (1-3). However, when IMPT and IMRT are performed on the head and neck region, the dose distribution deteriorates because of anatomical changes and setup errors occurring over the treatment period (4-6). Moreover, these anatomical changes have a greater impact on dose delivery in proton beam therapy compared with photon beam therapy (7-11). Robust planning is therefore employed in IMPT to compensate for the effects of anatomical variations and setup errors (12-16). However, although the plan so produced may be robust in countering anatomical variations and positioning errors, the dose to organs at risk (OARs) is generally increased, and the plan cannot accommodate large changes beyond the compensation range. Adaptive radiation therapy (ART), which uses a rescan computed tomography (CT) acquired during the treatment period to modify the plan, has also been employed to update the treatment plan for the effects of anatomical variations (17, 18). Several studies using ART have been reported, and various strategies have been devised (19). One of the typical ART strategies is the online adaptive strategy, in which the dose distribution is updated in daily treatment and provides a more robust and superior dose distribution (18, 20, 21). However, online adaptive ART requires high-quality imaging in the treatment room, creation of the clinical target volume (CTV) and OAR contours during treatment, recalculation of the particle transport and optimization, dose distribution verification, and physician approval. This requires advanced technology and is resource-intense. In contrast, image-guided radiation therapy (IGRT) and patient repositioning are simpler adaptive treatments accounting for anatomical changes, which do not require complicated procedures and can be applied in most clinics. However, because IGRT and repositioning are performed on the basis of the therapist’s experience, the quality of the procedures depends on the operator. In this study, we searched for anatomical structures showing strong associations with deterioration in the dose distribution to maximize the improvement of the dose distribution through IGRT and repositioning. Furthermore, we performed image registration using regions of interest (ROIs) fitted to highly relevant anatomical structures (target-based registration) and evaluated the effectiveness of this target-based registration by comparing it with a conventional method.
Patients and Methods
Patients and imaging datasets. Twenty patients with head and neck cancer were enrolled (Table I), and 20 simulation CT and 48 rescan CT images were retrospectively analyzed. The patients were treated with IMRT or IMPT, and CT images were acquired 2-4 times during the treatment period. All patients were immobilized using a four-point thermoplastic mask covering the shoulder and a patient-specific pillow. CT datasets were acquired using an Aquilion One scanner (Canon Medical Systems, Tochigi, Japan). Rescan CT images of six patients were acquired once at the 20th treatment, and rescan CT images of 14 patients were acquired at the 10th, 20th, and 30th treatments. Patient information that could identify an individual was anonymized. Although patient-informed consent was not required, details of this study were published on the center’s homepage and patients were given the opportunity to refuse to participate in the study. The contents of this study, including the investigation procedure and the handling of patient information, were approved by the Institutional Review Board of the National Cancer Center Hospital East (IRB No. 2020-282).
Patient characteristics.
Treatment planning. Forty volume modulated arc therapy (VMAT) and IMPT plans based on the initial CT images of the 20 patients were created. VMAT plans were created using two arcs, and the IMPT plans were created using 2-4 fields. All treatments were planned using an Eclipse treatment planning system (Varian Medical Systems, Palo Alto, CA, USA) by a well-trained therapist and medical physicist. The CTV and OARs were contoured on the simulation CT images by a radiation oncologist, and the CTV and OARs on the rescan CT were defined by the medical physicist using a deformed ROI based on the contouring by the radiation oncologist. All treatment plans were generated with 35 fractions of 2 Gy [relative biological effectiveness (RBE)] to give a total dose of 70 GyRBE. IMPT plans were generated to include a robustness to range uncertainty of 3.5% and isocenter uncertainty of 3 mm, using multiple field optimization. All plans were designed in accordance with the clinical protocol: CTV (high-risk), >98% of the volume to receive >95% of 70 GyRBE; CTV (low-risk), >98% of the volume to receive >95% of 54 GyRBE; CTV, ≤2% to receive <110% of 70 GyRBE; maximum dose to spinal cord <45 GyRBE; maximum dose to brain stem <54 GyRBE; mean dose to parotid gland <26 GyRBE; dose to cochlea <40 GyRBE; mean dose to oral cavity <35 GyRBE; ≤2% of lower jawbone to receive <105% of 70 GyRBE.
Anatomical structures and corrections. The typical anatomical structures of the head and neck region were selected for analysis: The clivus, second cervical vertebra, fifth cervical vertebra, first thoracic vertebra, lower jawbone, sternum, hyoid bone, and both shoulders. The positioning errors were obtained by calculating the differences in the six-axis vectors [anterior–posterior (A-P), superior–inferior (S-I), left–right (L-R), yaw, pitch, and roll] for each anatomical structure between the initial and replanning CT images. The positioning error determined according to the imaging positions in both CT images was corrected by triaxial translation. Calculation of the position error was performed by placing a rectangular ROI on the anatomical structure and using rigid image registration to align the ROIs in the two image sets. MIM Maestro® software (MIM Software Inc., Cleveland, OH, USA) was used for the rigid image registration algorithm. The three-dimensional vector of the position error for translational displacements (A-P, S-I, R-L) was calculated using Equation 1 below.
(Eq. 1)
Correlation coefficients between anatomical structures and dose indices. Pearson’s correlation coefficients were calculated using a combination of the eight anatomical structures and their respective four correction factors (3D vector, yaw, pitch and roll) and the 12 dose indices (CTVhigh D98, CTVlow D98, CTVhigh D2, spinal cord D1cc, parotid gland D50, oral cavity D50, brain stem D1cc, lower jawbone D2, larynx D50, conformity index, maximum dose, minimum dose in CTV). Differences in dose index were calculated for VMAT and IMPT by subtracting the respective recalculated plans based on the rescan CT from the initial plans. Differences in the corrected amounts and dose indices were calculated for 48 sets of initial CT and rescan CT, and correlation coefficients between each anatomical structure and important dose indices on the dose distribution were obtained.
Automatic image registration. Automatic registration was performed using the box image registration function of MIM Maestro®, and was used to calculate the dose distribution after the six-axis correction performed at the treatment couch. The target ROIs for automatic registration were set as follows: i) The entire patient image (all-image registration); ii) limited to the range covering the PTV (PTV registration), iii). limited to the clivus-C6 range (spine registration). Automatic registration was performed by transforming the rescan CT to the initial planning CT using each target ROI, and the CT after six-axis correction was reconstructed as a new CT using MIM Maestro®. The slice thickness of the new reconstructed CT images was 2 mm, which was the same as the initial planning CT. The new reconstructed CT images were transferred to the Eclipse treatment planning system and the treatment plans were recalculated using the initial VMAT and IMPT plans. The dose distributions obtained from the three automatic registration methods (six degrees of freedom) and three-axis correction were compared with the initial plan dose index.
Statistical analysis. Pearson’s correlation coefficient was considered weak, intermediate, and strong at <0.4, 0.4-0.7, and >0.7, respectively. The average dose differences from the initial plan of VMAT or IMPT were assessed using pairwise comparisons with t-tests. Bonferroni corrections were applied for multiple comparisons. All the statistical analyses were two-sided and were performed using R v. 4.2.0; statistical significance was set at p<0.05.
Results
Position error of the anatomical structures according to elapsed treatment days. Figure 1 and Supplementary Table I show the position error of each of the anatomical structures on CT images acquired at 10, 20, and 30 fractions. In the R-L direction, the position error was small, and was not dependent on the treatment day. In the S-I direction, the positional errors of the hyoid bone and shoulder were relatively large. In the A-P direction, the vertebrae (C2, C5, TH1), sternum, and shoulder position errors were large, and those of the C2, C5, TH1 and sternum increased in association with elapsed treatment days from initial treatment. In the yaw and roll directions, position error was random, and was not dependent on the treatment day. In the pitch direction, the errors for C2, TH1, hyoid bone and sternum increased with the elapsed treatment days. especially those for C2, for which the average position error exceeded 2.5 degrees at 30 fractions.
The average position error of anatomical structures in translation and rotation from simulation computed tomography for 20 patients. In the vertical plane and pitch, increasing position error of spinal structures was seen with elapsed treatment.
Correlations between dose index and anatomical position errors. Supplementary Table II and Supplementary Table III show the correlation coefficients between anatomical position errors and dose indices for VMAT and IMPT planning. For the CTV of high-risk coverage, the correlation coefficients of C2-pitch, C5-3Dvector, and TH1-3Dvector were −0.46, −0.43, and −0.43, respectively for VMAT, and – 0.39, −0.41, and −0.35 for IMPT, indicating intermediate correlations. These anatomical indices also negatively correlated with the minimum CTV dose. The correlation coefficient between hyoid bone-pitch and the CTV for high-risk D98 was −0.44 for VMAT, while the correlation coefficients between hyoid bone-3D vector and the conformity index were 0.47 for VMAT and 0.38 for IMPT. For the OAR dose, the correlation coefficients between parotid gland D50 dose and C2-roll, lower jawbone-roll, and clivus-roll were 0.55, 0.74, and 0.74, respectively, for VMAT, and 0.72, 0.71, and 0.81 for IMPT, indicating strong correlations. The correlation coefficients between oral cavity dose D50 and C2-roll, lower jawbone-roll, and clivus-roll were 0.74, 0.53, and 0.65, respectively, for VMAT, but no strong correlation was found for IMPT. The correlation coefficients between spinal cord D1cc and the C5-3Dvector and TH1-3Dvector were 0.39 and 0.45, respectively, for VMAT, but there were weak correlations for IMPT. For the hyoid bone, sternum, and shoulder, weak correlations with dose indices was found for either VMAT or IMPT.
Image registration strategy and dose index. We evaluated three-axis image registration and three methods for six-axis automatic image registration with different registration ROIs. Figure 2 shows boxplots of the CTV high-risk D98, CTV high-risk D2, CTV low-risk D98, spinal cord D1cc, parotid D50, and oral cavity D50 in VMAT. In the CTV high-risk D2, the average increase was more than 1% with the three-axis correction, but less than 0.5% with all three methods for six-axis correction. The coverage of CTV high-risk D98 improved slightly with six-axis correction compared with three-axis correction. For spinal cord-D1cc, the average dose increased by 2 Gy, and the maximum average dose increased by more than 12 Gy with three-axis correction. Six-axis correction with spine registration showed the best performance in respect to the mean and standard deviation of the dose. For parotid D50, an average dose increase of 2 Gy and a maximum average dose of more than 25 Gy was found with the three-axis correction. Of the three six-axis methods, the six-axis correction with spine registration provided the lowest increase in the mean parotid dose. Figure 3 shows boxplots of CTV high-risk D98, CTV high-risk D2, CTV low-risk D98, spinal cord D1cc, parotid D50, and oral cavity D50 in IMPT. Coverage of the CTV high-risk D98 improved slightly with six-axis correction compared with three-axis correction. However, the other dose indices showed no significant differences between the six-axis correction and three-axis correction.
The average dose difference from the initial plan in volumetric-modulated arc therapy. Boxplots represent each image registration method focusing on different regions of interest (ROIs). Spine: Translation and rotation correction focusing on spinal ROIs, 6 degrees of freedom (DoF); PTV: 6 DoF focusing on planning target volume ROI; All: 6 DoF using all image information (without focusing on ROIs); 3DoF: translation correction only. The boxes represent the interquartile range and the horizontal lines the median values, crosses indicate the average values. Whisker lines include data within the 1.5-fold of the interquartile range; dots are outliers. CTV: Clinical target volume; D2/50/98: dose to 2%/50%/98% of the volume; D1cc: Dose delivered to 1 cm3 volume. *Significantly different.
The average dose difference from the initial plan in intensity-modulated proton therapy. Boxplots represent each image registration method focusing on different regions of interest (ROIs). Spine: Translation and rotation correction focusing on spinal ROIs, 6 degrees of freedom (DoF); PTV: 6 DoF focusing on planning target volume ROI; All: 6 DoF using all image information (without focusing on ROIs); 3DoF: translation correction only. The boxes represent the interquartile range and the horizontal lines the median values, crosses indicate the average values. Whisker lines include data within the 1.5-fold of the interquartile range, dots are outliers. CTV: Clinical target volume; D2/50/98: dose to 2%/50%/98% of the volume; D1cc: Dose delivered to 1 cm3 volume.
Discussion
In this study, setup errors and their effects on dose distribution were retrospectively investigated in 20 patients with head and neck cancer. The relationship between setup error and dose deterioration was previously investigated in patients with nasopharyngeal cancer (22), but the study did not include rotational motion as a setup error, and position error was evaluated only at the isocenter position. In contrast, we investigated the relationship between the six-axis movement of each anatomical structure and the dose distribution. In addition, we evaluated the effectiveness of target registration focusing on identified anatomical structures. In the relationship between position error and treatment day, the position errors of C2-pitch and C2, C5, and TH1 in the vertical direction increased with the elapsed treatment days. We believe this may be because the thickness of the posterior neck and back decreased due to weight loss during the treatment period, and the vertebrae were therefore displaced downward. C2-pitch, C5-3D vector, and TH1-3D vector were associated with decreased CTV coverage. Therefore, the positional errors of C2-pitch, C5-3D vector, and TH1-3D vector increased because of the decrease in thickness of the posterior neck and back, and changes in these anatomical structures may be factors that can predict a decrease in CTV coverage. Spine registration with an emphasis on the vertebrae improved the CTV coverage. We believe this was because the spine registration corrected the downward displacement of the vertebrae and brought the vertebral alignment at treatment closer to that of the initial CT. The C5-3D vector and TH1-3D vector correlated with the dose increase to the spinal cord, and the spine registration suppressed the dose increase to the spinal cord the most out of the three methods evaluated. The increase in spinal cord dose is caused by a shift of the high-dose area towards the vertebrae, which in turn is caused by displacement of the vertebral alignment due to the reduction in thickness of the posterior neck and back. We believe that the spine registration corrected the downward displacement of the vertebrae well, and thereby suppressed the increase in the spinal cord dose.
An increase in parotid gland dose strongly correlated with C2-roll, lower jawbone-roll, and clivus-roll. A high correlation between mandibular rotation and parotid gland dose is consistent with the report by Otsuka et al. (23). However, the position errors in the roll direction were less related to the elapsed treatment time. Position errors of the C2, mandible, and clivus in the roll direction mean that the orientation of the patient’s face was different from that at the time of the initial planning. Because no correlation with the elapsed treatment time was observed, these errors appear to have been caused by daily patient positioning. Change in a patient’s position resulting in rotation of the face with the body facing forward cannot be corrected by six-axis correction (24). Spine registration was able to minimize the increase in the mean dose to the parotid gland; since the image registration range was limited to the vertebrae (C1-6) in the spine registration, anatomical structures such as the clivus and C2 were emphasized in the six-axis correction, and we believe that the orientation of the patient’s face was therefore well corrected. In the VMAT plans, the spine registration and CTV registration (six-axis correction using limited ROI) contributed to improved CTV coverage, fewer hot-spots, and lower OAR dose. In comparison, in the IMPT plans, the range-limited six-axis correction was somewhat effective at improving the CTV coverage but did not contribute to reducing the CTV D2 or OAR dose. This may be related to the characteristics of the IMPT plans with robust planning and means that the effects of alignment distortion and changes in body shape that cannot be corrected by image registration are greater with IMPT plans than with VMAT plans, and that six-axis correction cannot correct them. Deterioration of the dose distribution due to anatomical changes was more severe in IMPT than in VMAT, consistent with the report of Góra et al. (17). Therefore, in head and neck IMPT, it is more important to correct vertebral alignment distortion by repositioning and to correct body thickness changes by ART.
The relationships between the position error of each anatomical structure and dose deterioration found in this study can be applied to optimization of the margin size. A previous study reported that anisotropic margins had lower margin size than isotropic margins (25). Therefore, it is possible that expanding the margin for the CTV around the vertebrae, clivus and mandible may help suppress deterioration of the CTV coverage. However, in the treatment of nasopharyngeal cancer, OARs such as the cochlea, brainstem and chiasm are in close proximity and margin expansion is often not possible, although there is a report that uncompromised IGRT is possible using an anisotropic margin generated with emphasis on the OARs (26). However, the CTV for head and neck cancer is generally broad and there are multiple OARs, which means that this method cannot be used. According to our results, the anisotropic margins at the initial planning may not be optimal because the positional variation of each anatomical feature increased with the number of treatment days and varied from patient to patient. Therefore, monitoring for vertebral alignment, and repositioning the patient and triggering the ART timing according to it are important to ensure the quality of dose distributions (27).
This study has several limitations. The number of cases examined was limited to 20 cases. The 10th, 20th, and 30th CT images included the effects of both positioning errors and changes in body shape, and it was not possible to separate these effects in the dose index evaluation. We used MIM Maestro® for automatic image registration, other automatic registration algorithms may provide different results.
Conclusion
We searched for anatomical structures that affect the dose indices in IMPT and VMAT. CTV coverage was intermediately correlated with positioning errors of vertebrae and OAR doses were strongly correlated with positioning errors of the C2, clivus, and mandible. In particular, the positional errors in the pitch and vertical directions of the vertebrae increased with the elapsed treatment days, and may be affected by changes in body shape. Six-axis correction improved the CTV coverage and OAR dose in VMAT.
Acknowledgements
The Authors wish to express sincere gratitude to Yoshihisa Muramatsu Ph.D. for supporting the research facilities and the environment for conducting this research. We thank Edanz (https://jp.edanz.com/ac) for editing a draft of this article.
Footnotes
Authors’ Contributions
Concept and design: KH, KT and SM. Treatment planning: KH. Data analysis and interpretation: KH and KT. Research management and supervision: MI and TS. All Authors read and approved the final article.
↵Supplementary Material
Supplementary Table I available at: https://www.dropbox.com/s/erfwoly8frzkxia/Supplementary%20table1.docx?dl=0
Supplementary Table II available at: https://www.dropbox.com/scl/fi/75jeexlcgi997blnxb0od/Supplementary-table2.docx?dl=0&rlkey=w319b9n6fz7myn0umc7mdifey
Supplementary Table III available at: https://www.dropbox.com/scl/fi/jmlas5objihsohzktw9qc/Supplementary-table3.docx?dl=0&rlkey=50pcpkbfu9poyus6cd5u21ai1
Conflicts of Interest
The Authors declare no conflicts of interest.
- Received January 17, 2023.
- Revision received February 5, 2023.
- Accepted February 9, 2023.
- Copyright © 2023 International Institute of Anticancer Research (Dr. George J. Delinasios), All rights reserved.









