Determining appropriate timing of adaptive radiation therapy for nasopharyngeal carcinoma during intensity-modulated radiation therapy

Background To determine appropriate timing of an adaptive radiation therapy (ART) replan by evaluating anatomic and dosimetric changes of target volumes and organs at risk (OARs) during intensity-modulated radiation therapy (IMRT) for nasopharyngeal carcinoma (NPC). Methods Nineteen NPC patients were recruited. Each patient had repeat computed tomography (CT) scans after each five fractions and at treatment completion. Automatic re-contouring the targets and OARs by using deformable registration algorithm was conducted through CT-CT fusion. Anatomic changes were assessed by comparing the initial CT and repeated CT. Hybrid plans with re-contouring were generated and the dose-volume histograms (DVH) of the hybrid plan and the original plan were compared. Results Progressive volume reductions in gross target volume for primary disease (GTVnx), gross target volume for involved lymph nodes (GTVnd), and parotids were observed over time. Comparing with the original plan, each hybrid plan had no significant difference in homogeneity index (HI) for all the targets. Some parameters for planning target volumes for primary disease and high-risk clinical target volume (PTVnx and PTV1, respectively) improved significantly, notably starting from the 10th fraction. These parameters included mean dose (Dmean), dose to 95 % of the volume (D95), percentage of the volume receiving 95 % of the prescription dose (V95), and conformity index (CI) for PTVnx, and Dmean, D95, and CI for PTV1. The dosimetric parameters for PTVnd remained the same in general except for D95 and V95 which had significant improvement at specific time points; whereas for PTV2, similar trend of dosimetric changes was also observed. Dose to some OARs increased significantly at some time points. Conclusions There were significant anatomic and dosimetric changes in the targets and OARs. The target dose coverage in the hybrid plans did not get worse, but overdose occurred in some critical structures. Significant dosimetric changes should be considered as a trigger point at which ART replanning is indicated. D95/V95/CI for PTV2, Dmax for the brain stem, spinal cord, right eyeball and left lens, and Dmean/V30 for the parotids and glottis were taken into account for predicting the need for ART. Two replans at the 5th and 15th fractions were suggested.


Background
Radiation therapy (RT) is considered a main treatment approach in the management of nasopharyngeal carcinoma (NPC). Intensity-modulated radiation therapy (IMRT), formally introduced in clinical practice since early 1990s, has now replaced conventional radiation therapy (CRT) and three-dimensional radiation therapy (3D-CRT) and become a standard treatment for NPC. This technique provides adequate target coverage while maintaining steep dose gradients at the border between the target and adjacent normal tissues, through which dose escalation for the targets may be achieved without delivering excessive dose to the organs at risk (OARs) [1]. Numerous studies have demonstrated that IMRT results in a decreased incidence of radiation-induced side effects and an improvement in quality of life (QOL) for NPC patients, comparing with CRT or 3D-CRT [2][3][4][5][6].
However, NPC patients may experience significant anatomical changes throughout the entire treatment course. These changes include the shrinkage of the primary disease and metastatic lymph nodes, external contour because of significant weight loss, and displacement/size of the normal structures. As a consequence of the anatomic changes, significant difference between the actual delivered dose and the initial planned dose would be anticipated, which may ultimately result in underdose to the targets and/or overdose to the critical structures. Wang et al. [7] compared a repeated CT after treatment at a dose of 40 Gy with the initial planning CT for 20 patients with locally advanced NPC, and found that the dose coverage of the targets remained unchanged; whereas the dose delivered to the parotid glands and spinal cord increased significantly. Excessive irradiation to OARs increased the risk of incidence of radiation-induced complications such as xerostomia and myelitis, which may have a great negative impact on QOL [8][9][10][11][12].
It is suggested that the initial planning based on pretreatment condition may not truly reflect the dosimetric variations during the course of IMRT. Thus adaptive radiation therapy (ART), a plan modification and implementation according to tumor response and anatomic changes of normal structures, becomes particularly important. In recent years, many studies have focused on ART for NPC patients, however, the optimal timing and frequency of ART remains unanswered [13][14][15][16][17][18][19]. In this prospective study, we aimed to determine appropriate timing to perform an ART replan by evaluating anatomic and dosimetric changes of the target volumes and OARs during IMRT for NPC.

Eligibility criteria
Patients with histologically proven NPC and treated with curative IMRT were enrolled into this prospective study. Inclusion criteria were as follows: age ≥18 years; Eastern Cooperative Oncology Group (ECOG) Performance Status 0-2; stages I-IVb according to the 2010 AJCC Staging System. Patients diagnosed with, or treated for other malignances, or treated with non-IMRT techniques were excluded in the study. Written informed consent was obtained for all patients. The study was approved by the Institutional Review Board (IRB) .

Immobilization and simulation
All patients were immobilized in a supine position with the head in a neutral position with a tailored thermoplastic mask covering the head, neck, and shoulders. Intravenous contrast-enhanced CT using 2 mm slice from the vertex to the manubriosternal joint was performed for planning. The CT data were imported to the a treatment planning system (Pinnacle 3 , version 9.2).

Delineation of the targets and OARs
The target delineation was in accordance with the International Commission on Radiation Units and Measurements Reports 50 and 62. Briefly, the primary gross volume (GTVnx) and the involved cervical lymphadenopathy (GTVnd) included all known gross disease as determined by the imaging, clinical, and endoscopic findings. The high-risk clinical target volume (CTV1) was defined as GTVnx plus 5-mm margin and entire nasopharyngeal mucosa plus 5-mm submucosal volume. The low-risk clinical target volume (CTV2) covered CTV1, entire nasopharynx, parapharyngeal space, pterygopalatine fossa, posterior third of the nasal cavity and maxillary sinuses, inferior sphenoid sinus, posterior ethmoid sinus, skull base, and anterior half of the clivus. CTV2 also covered elective neck nodal regions, including bilateral retropharyngeal lymph nodes and ipsilateral levels II, III, and Va for nodenegative neck, or full length of ipsilateral neck for nodepositive neck. Level Ib was not routinely irradiated unless there was confirmed lymphadenopathy in level Ib, or large metastatic node size (≥3 cm)/extracapsular spread was present in level IIa. PTVnx, PTVnd, PTV1, and PTV2 were generated by adding 5-mm margin to GTVnx, GTVnd, CTV1, and CTV2, respectively. Care was taken to ensure at least 5-mm gap was present between the PTVs and the skin. The contoured critical structures included the brain stem, chiasm, optic nerves, spinal cord, eyes, lens, parotid glands, oral cavity, larynx, mandible, and temporomandibular joints.

Treatment design and delivery
The plans were designed and optimized using the Pinnacle inverse planning system. The prescribed radiation dose was 69.76 Gy at 2.18 Gy per fraction delivered to the PTVnx and PTVnd, and 60.8 Gy at 1.9 Gy per fraction delivered to the PTV1. The PTV2 was treated to 54.4 Gy at 1.7 Gy per fractions. All patients were treated once daily, five fractions weekly. Dose constrains to the critical structures were within the tolerance according to the RTOG 0225 protocol, and efforts were made to meet the criteria as closely as possible. IMRT was delivered via seven fixedgantry angles with an Elekta Synergy Linear Accelerator (Elekta Ltd.) with step-and-shoot treatment techniques.
Acquiring and processing CT data during IMRT Repeat CT scans were acquired for each patient with the same mask and isocenter as the initial simulation CT scan after each five fractions and at the end of the full treatment course using a 24-slice CT scanner (Somatom Sensation Open, Siemens Medical Solutions, Erlangen, Germany). The datasets were denoted as CT-1, CT-2, CT-3, CT-4, CT-5, and CT-6, respectively. Each new CT dataset was registered with the initial planning CT dataset through VoxAlign Deformation Engine provided by the MiM Maestro software. Auto-propagating the planning contours on the new CT was conducted and manual modification was performed if needed. Changes in the volume of GTVnx, GTVnd, and parotid and shift of the parotid centroid were calculated by comparing the new CT and the planning CT.
A hybrid IMRT plan was generated by superimposing the initial treatment plan (Plan-0) to each repeated new CT image. After recalculation of the dose distribution on the new CT images, the dose delivered by the hybrid plans to the redelineated target volumes and OARs were recorded and analyzed (adaptive replanning was not suitable for CT-6 since the CT-6 images were acquired at the end of IMRT. There were only 5 hybrid plans, namely Plan-1, Plan-2, Plan-3, Plan-4, and Plan-5).

Statistical analysis
The Kolmogorov-Smirnov test was used to test data for normality. Mean ± standard deviation was used for data with normal distribution; whereas median (interquartile range) was used for data with skewed distribution. A paired sample T-test or Wilcoxon rank-sum test was chosen based on the data types. The correlation between weight loss and parameters associated with anatomic changes was estimated using Pearson's correlation. Analysis of variance (ANOVA) was performed to determine significant dosimetric changes at any time point (trigger point) of the entire treatment course. A probability value less than 0.05 was considered significant. Analyses were performed using SPSS 19.0 software (SPSS, Inc., Chicago, IL, USA).

Patient characteristics
Between August 2012 and December 2013, a total of 19 patients diagnosed with undifferentiated non-keratinizing NPC were enrolled into this study. There were 11 men and 8 women with median age of 46 years (range, 22-70). Stage distributions according to the 2010 AJCC Staging System were as follows: stages II, 3 patients; stage III, 13 patients; and stage IV, 3 patients. Concurrent platinumbased chemotherapy was given to 17 patients, and no chemotherapy to 2 patients. The characteristics of the patient cohort are listed in Table 1.

Changes in volume of GTVnx and GTVnd
Steady volume reduction in GTVnx and GTVnd was observed over time. Compared with the baseline, the volumes of GTVnx and GTVnd at the end of IMRT were decreased significantly by 65.6 ± 13.3 % and 72.7 ± 13.3 % (p = 0.001 and 0.046), corresponding to a reduction rate of 17.3 ± 10.4 %/week and 22.2 ± 15.0 %/week, respectively.

Changes in volume and displacement of the parotid
At the treatment completion, the volume of the left parotid decreased by 38.0 ± 15.3 %, corresponding to a reduction rate of 7.9 ± 9.6 %/week; whereas for the right parotid, the volume reduction was 39.2 ± 14.7 %, corresponding to a reduction rate of 7.8 ± 11.3 %/week. The most significant volume reduction occurred at the third week, with its value of 15.8 ± 6.8 % (p = 0.016) and 14.4 ± 7.6 % (p = 0.003) for the left and right parotids, respectively.
The displacements of the centroid for the left parotid were −0.27 ± 0.25 cm, 0.04 ± 0.12 cm, and 0.14 ± 0.20 cm in the right-left (RL), anterior-posterior (AP), and superior-inferior (SI) directions, respectively. Significant shifts were noticed in the RL and AP directions (p = 0.001 and 0.018). While for the right parotid, the displacements of its centroid were 0.31 ± 0.19 cm, 0.01 ± 0.21 cm, and 0.03 ± 0.20 cm in the RL, AP, and SI directions, respectively. Significant shift occurred in the AP direction (p < 0.001). The centroid displacement for the parotids is illustrated in Fig. 1.

Dosimetric changes in the target volumes
There was no significant difference in each homogeneity index (HI) for all the targets between the original plan and the hybrid plan. Some parameters for PTVnx and PTV1 improved significantly as treatment continued, notably starting from the 10th fraction. These parameters included mean dose (Dmean), dose to 95 % of the volume (D95), percentage of the volume receiving 95 % of the prescription dose (V95), and conformity index (CI) for PTVnx, and Dmean, D95, and CI for PTV1. Other parameters including percentage of the volume receiving >110 % of the prescription dose (V110) for PTVnx, and V95 and V110 for PTV1 changed less or not significantly. The dosimetric parameters for PTVnd remained the same in general except for D95 which had significant improvement at the 15th and 25th fraction, and for V95 which had significant improvement starting from the 10th fraction to 25th fraction; whereas for PTV2, a similar trend of dosimetric changes was also observed (Table 3).

Dosimetric changes in OARs
Dmean for the brain stem and spinal cord increased significantly over the IMRT course, occurring as early as at the 15th and 5th fractions, respectively. Dmean and V30 for the left and right parotids increased significantly, beginning at the fifteenth and twentieth fraction, respectively. For the glottis, right eye, and left lens, significant increase in Dmax or Dmean was found in some time point (s) during the treatment. Table 4 illustrates the dosimetric changes in some OARs. No significant differences were found in dosimetric parameters between the initial treatment plan and the hybrid plan in other OARs including the temporal lobes, optic nerves, optic chiasm, left eye, right lens, cochleae, mandible, temporomandibular joints, esophagus, oral cavity, supraglottis, and subglottis (not shown in Table 4).

Timing of adaptive treatment replanning
It is reasonable to postulate that any anatomic changes during the treatment course may not certainly result in marked dosimetric changes. So ART may not be necessary for a patient who has no significant dosimetric changes at some time point even if statistical anatomic changes exist at the same time. Based on this concept, we considered a trigger point at which significant dosimetric variation for a specific parameter was present as an indicator for ART replanning. First of all, dosimetric parameters in Plan-0 were compared with those in the following hybrid plans separately until statistical differences were found. For example, if significant difference was found in the target dose coverage or overdose to the OARs between Plan-0 and Plan-1, adaptive replanning would be initiated after the 5th fraction (here the 5th fraction was considered a trigger point). Then parameters in Plan-1 were compared with those in the following hybrid plans until another trigger point was found, and so on. In this way, the trigger points for each selected parameter throughout the treatment course could be identified. Table 5 illustrates the suggested trigger points at which adaptive replanning would be designed. However, for practical consideration, two replans at the 5th and 15th fractions were proposed. This was largely because most of the trigger points occurred at these time points.

Discussion
Similar to other head-and-neck cancers, patients with NPC disease also experience marked changes during IMRT. These changes include shrinkage of the primary disease and involved lymph nodes, weight loss, and geometric/volumetric variations of the OARs. Fung et al. [14] acquired daily megavoltage CT (MVCT) images and registered them to the corresponding planning CT images for 30 NPC patients treated with helical tomotherapy. At the end of treatment, there was a volume reduction of 35.70 ± 20.06 % for the posterolateral wall of nasopharynx (P-NP) which was a surrogate for the primary disease, corresponding to a daily reduction of 0.99 ± 0.55 %. In a study by Cheng et al. [20], the mean shrinkages of the nasopharyngeal disease and metastatic lymph nodes were 9.1 and 16.2 % at 30 Gy, and 13.1 and 28.7 % at 50 Gy, respectively. Barker et al. [19] found that the GTV for head-and-neck cancer treated with definitive external beam RT decreased throughout the course of fractionated RT, at a median rate of 0.2 cm 3 . On the last day of treatment, the median volume loss of the initial GTV was 69.5 % (range, 9.9-91.9 %). This was consistent with our findings. In the present study, progressive tumor regression for GTVnx and GTVnd was observed over time. The volumes of GTVnx and GTVnd at the end of IMRT were decreased significantly by 65.6 ± 13.3 % and 72.7 ± 13.3 % (p = 0.001 and 0.046), corresponding to reduction rates of 17.3 ± 10.4 %/week and 22.2 ± 15.0 %/week, respectively.
Weight loss is a common event for head-and-neck cancer patients when RT is given with or without concurrent chemotherapy. Ng et al. [21] found that at the end of RT, 82 % of NPC patients had significant weight loss and were in negative energy balance, which persisted for more than 6 months. Although NPC patients who received IMRT lost less weight than those who received conventional RT, weight loss in patients treated with IMRT should not be underestimated. In a study by Qiu et al. [22], weight loss of 5.81 ± 2.34 kg was observed. Cheng et al. [20] found that most NPC patients who were treated with IMRT experienced significant weight loss during the RT course. Compared with the baseline, mean weight loss at 30 and 50 Gy were 5.4 and 9.3 %, respectively. In the present study, we found an average weight loss of 1.5 ± 2.2 % per week. Besides, there was a significant weight loss at each time point from the first week to the end of treatment course, indicating that the patients experienced progressive weight loss. Patients in our study lost less body weight than those in Cheng's study. This may be attributed to the discrepancy in the patients' baseline physical status, severity of complications caused by treatment itself, and support care between the two studies. Significant weight loss during RT may have an negative impact on treatment outcome, including lower quality-of-life (QOL) scores, poorer treatment compliance, prolonged recovery time and hospital stay, and poorer prognosis [23][24][25][26].
The parotids are the most frequently changed organs in NPC patients throughout the treatment course, both in the volume and in the displacement. Wang et al. [7] showed a volume reduction of 14.7 % for the left parotid and 18.2 % for the right parotid during IMRT for locally advanced NPC by comparing the initial planning CT and a repeat CT obtained after treatment at a dose of 40 Gy. Lu et al. [16] found a volume reduction of 35.1 ± 20.0 % and 24.6 ± 11.9 % for the left and right parotids, respectively, after 25 fractions. Similar reduction trends for the parotid could also be reflected by another study in which the mean percentage volume loss was 47.54 ± 14.27 % at the treatment completion, and the mean loss rate was estimated to be 1.35 ± 0.39 %/day. In addition, the mean center of mass (COM) of the parotid shifted progressively towards the medial and superior aspects during treatment. Also the mean medial and superior displacement for both sides of the parotid were 0.34 ± 0.27 cm and 0.24 ± 0.39 cm, respectively (both P < 0.001) [14]. In our study, the volume of both sides of the parotid at the treatment completion decreased significantly, comparing with the baseline. We found that the most significant  volume reduction occurred at the 3rd week. Apart from volumetric variation, the displacement of the centroid for the parotids also changed. Significant shifts were found in the AP direction for the right parotid, and in the RL and AP directions for the left parotid. It has been reported that shrinkage of the parotid and shift of its centroid are correlated with weight loss. Barker et al. [19] found that in head-and-neck cancer, the medial shift in parotid COM correlated highly with weight loss during fractionated RT course. Similar findings were also shown in our study. Weight loss significantly correlated with the shrinkage of parotids and their centroid shift towards the medial direction, but no correlation was observed between the weight loss and the volumetric changes of GTVnx and GTVnd. The parotids may be exposed to a high dose region due to their medial shift or volume reduction during the IMRT course. Significant weight loss may be used as a surrogate for a large medial parotid shift or larger parotid volume loss, implying that modified replanning may be needed to avoid excessive irradiation to the parotid.
Significant anatomic changes during the treatment course in NPC patients may affect the dose distribution. The initial plan based on pretreatment CT images that are only snapshots of the patient's anatomy at a static time point do not precisely reflect the actual dose distribution during fractionated IMRT. Fung et al. [13] created two new adaptive plans (PII-ART and PIII-ART) for 10 NPC patients undergoing Hi-Art Tomotherapy based on up-to-date CT images and contours and used these plans for treatment in phase two (PII; after 25th fraction) and phase three (PIII; after 35th fraction),  respectively. Two hybrid plans (PII-NART and PIII-NART) were generated using the original contours pasted on the PII-and PIII-CT sets by CT-CT fusion. Dosimetric comparisons were made between the NART plans and the corresponding ART plans. They found that without replanning, the doses to D95 for all the target volumes were increased with better dose uniformity, whereas the OARs received higher doses compared with the corresponding ART plans. The total dose to D1 for the brain stem and spinal cord significantly increased by 7.87 ± 7.26 and 10.69 ± 6.72 %, respectively. There were also significantly increased maximum doses to the optic chiasm and pituitary gland, and significantly increased mean doses to both sides of the parotid. Cheng et al. [20] performed CT and MRI scans at 30 Gy and 50 Gy intervals for 19 NPC patients treated with IMRT. When comparing the initial plan with the hybrid plans which were generated by superimposing the initial plan to the repeat CT images, they found that the hybrid plans demonstrated significantly higher dose to most of target volumes with greater dose inhomogeneity, higher maximum doses to the spinal cord and brainstem, and higher median doses to the parotid. Wang et al. [7] noticed that after IMRT at the 40 Gy dose point, dose coverage of all the targets remained unchanged, whereas the dose delivered to the parotid glands and spinal cord increased significantly. The results were in line with ours in this study. Findings mentioned above indicated that although the actually delivered doses to the targets may be higher or at least not worse than the planned ones, doses to OARs, particularly the parotids, spinal cord, and brain stem, however, significantly increased. Determining an appropriate time point at which ART is intervened in a timely manner is critical to ensure that the planned dose to the targets and OARs can be delivered faithfully throughout the entire IMRT course. Workflow for ART includes re-simulating, re-contouring, and treatment re-designing, which is a time-consuming and heavy workload process. Thus daily ART is not always applicable in many cancer centers. Many reports in recent years suggested that mid-to-late phase of the treatment course was appropriate timing for ART. Lu et al. [16] recommended that ART should be initiated after 25 fractions of IMRT to ensure adaptive doses to the targets and critical normal tissues. Cheng et al. [20] showed that doses to some OARs such as the brain stem, spinal cord and parotids were significantly increased at 30 Gy during the course of IMRT, while the target coverage remained adequate. So replanning was suggested at this time point. In another study by Fung et al., significant anatomic changes were found at the 9th, 19th, and 29th fractions. The authors thus recommended ART replanning at these three time points. Contrary to Fung's study, we selected parameters only related with dose distributions as the endpoint to determine whether a replan was needed. As explained in the previous section, significant anatomic changes may not certainly result in remarkable changes in dosimetric effects. So only dosimetric changes were chosen as determinant for ART replanning in this study. Bearing this concept in mind, we identified the trigger points at which ART replanning would be initiated. However, for practical consideration, two replans at the 5th and 15th fractions were proposed. This was largely because most of the trigger points occurred at these time points, as shown in Table 5.
It should be noted that all patients in the present study were treated with the original plan. The repeat CT scans during IMRT were only used for research, and no ART replanning was prepared for patient treatment. Whether ART replanning can transfer into dosimetric or clinical benefits remains unknown. To answer these questions, a further study is underway to investigate the impact of ART on dosimetric parameters and clinical outcomes.

Conclusions
In summary, there were significant anatomic and dosimetric changes in the targets and OARs during the course of IMRT. Weight loss significantly correlated with the shrinkage of parotids and their centroid shift towards the medial direction. The target dose coverage in the hybrid plans did not get worse, comparing with the original plan, but overdose occurred in some critical structures. Significant dosimetric changes should be considered as a trigger point at which ART replanning would be initiated. Based on the dosimetric analysis, D95/V95/CI for PTV2, Dmax for the brain stem, spinal cord, right eyeball and left lens, and Dmean/V30 for the parotids and glottis were taken into account for predicting the need for ART. Two replans at the 5th and 15th fractions were suggested for performing ART.