Patient selection
Based on previous reports in the literature12,14, we assumed that the mean dose to the oral cavity would be 20 Gy lower with a standard deviation of 12 Gy in the adaptive IMPT (A-IMPT) compared to the adaptive IMXT (A-IMXT). To determine the differences in the mean dose to the oral cavity, the number of patients required to assure 80% power to detect the difference at a two-sided significance level of 5% was calculated by paired t test16. The sample size was calculated to be six by web tool and SAS®version9.4(Cary, NC).
There were nine patients with NPC who were all treated with A-IMXT (step and shoot approach) and concomitant chemotherapy between 2015 May and 2016 July at Hokkaido University Hospital. Two patients were excluded because they received induction chemotherapy before the IMXT with concomitant chemotherapy. One patient was excluded because of treatment by step and shoot IMXT in the initial part but by the volumetric arc therapy technique (VMAT) in the subsequent adaptive part. Finally, the remaining six patients were included in the dosimetric comparison. The characteristics of these six patients are shown in Table 1. The TNM classification and the clinical stage in Table 1 was according to the AJCC TNM classification 7th edition.
Table 1
Patient number | Age | Gender | T stage* | N stage* | M stage | Clinical stage |
1 | 58 | Male | T1 | N2 | M0 | Ⅲ |
2 | 66 | Female | T1 | N2 | M0 | Ⅲ |
3 | 60 | Male | T1 | N1 | M0 | Ⅱ |
4 | 67 | Male | T3 | N1 | M0 | Ⅲ |
5 | 58 | Male | T1 | N2 | M0 | Ⅲ |
6 | 69 | Male | T3 | N1 | M0 | Ⅲ |
*UICC 7th edition
Target Volume, OAR, And Dose Prescription
In the initial plan, Gross tumor volume (GTV), consisted of both a primary lesion (GTV primary) and metastatic lymph nodes (GTV node). These were contoured using fusion images of computed tomography (CT) and magnetic resonance imaging (MRI). The CTV primary and CTV node were basically made by adding a 1.0 cm margin to the GTV and an additional margin for high risk patients accounting for microscopic extensions to the GTV primary and GTV node, respectively. Clinical target volumes in the initial plan (CTV1) consisted of the CTV primary, the CTV node, and the prophylactic regional lymph node. The prophylactic lymph nodes were contoured according to the international consensus guideline17. Planning target volume for the initial plan, PTV1, was produced by adding an additional 3.0 mm as the set-up margin to the CTV1.
To make the adaptive plan, patients were rescanned by CT and MRI at the time when they were given 24 Gy/12fr-28Gy/14fr. The clinical target volume in the adaptive plan (CTV2) which consisted of the CTV2 primary and the CTV2 node was contoured using the rescanned CT. The planning target volume for the adaptive plan (PTV2) was produced by adding an additional 3.0 mm as the set-up margin to the CTV2.
The brain stem, optic nerve, optic chiasma, oral cavity, thyroid gland, and supraglottic larynx, glottis larynx and superior pharyngeal constrictor muscle, middle pharyngeal constrictor muscle, inferior constrictor muscle and cricopharyngeal muscle were contoured for OARs. Basically, a 3 mm Planning organ at risk volume (PRV) margin were added to the brain stem, optic nerve, optic chiasm, and spinal cord.
A-IMXT Planning
Pinnacle3 v9.0 (Phillips, Medical Systems, WI) was used in the inverse treatment planning of IMXT using 7 portals and step and shoot techniques.
In the actual treatment, 46 Gy in 23 fractions (46 Gy/23Fr) was prescribed using the initial plan and 24 Gy/12Fr was prescribed using the adaptive plan.
Dose constraints for OARs in the inverse planning are shown in Supplementary table S1 (Spinal cord PRV ≤ 50 Gy, Brainstem PRV ≤ 60 Gy, Optic nerves PRV ≤ 60 Gy, Optic chiasm PRV ≤ 60 Gy). The plans were optimized to ensure 100% of the prescription dose to cover 95% of the volume of PTV1 and PTV2.
Plans to irradiate 70 Gy to the PTV1 were made based on the initial CT (initial plan) and used in the plans to irradiate 70 Gy to PTV2 on the second CT (adaptive plan). The dose to the OARs and CTV in the initial plan was estimated by dividing the dose to the OARs and CTV multiplied by 46/70. The dose to the OARs and CTV in the adaptive plan was estimated by dividing the dose to the OARs and CTV multiplied 24/70.
A-IMPT Planning
The VQA version 3.077 (Hitachi co ltd., Tokyo) was used in the inverse treatment planning of the IMPT using 3 portals at 80, 160, and 200 degrees. The dose distribution was calculated using the beam profile of the spot scanning dedicated proton beam therapy system, PROBEAT-RT, (Hitachi, co ltd., Tokyo).
Relative biological effectiveness (RBE) was assumed to be 1.1 throughout the irradiated volume and we used 1GyE, which is the physical dose times 1.1, as the unit of absorbed dose in the treatment planning. In the actual treatment, 46GyE in 23 fractions (46GyE/23Fr) was prescribed using the initial plan and 24GyE/12Fr was prescribed using the adaptive plan.
Robust optimization was used to prescribe at least 70GyE at D99% of CTV by shifting the CT images by 3 mm and assuming an uncertainty in the range of 3.5%.
The dose constraints for OARs in the inverse planning are shown in Table 2 (Spinal cord PRV ≤ 50GyE, Brainstem PRV ≤ 60GyE, Optic nerves PRV ≤ 60GyE, and Optic chiasm PRV ≤ 60GyE). Plans were optimized to ensure 100% of the prescription dose to cover 99% of the volume of CTV1and CTV2.
For the initial CT, plans to irradiate 70GyE to CTV1 were made (initial plan) and for the second CT, plans to irradiate 70GyE to CTV2 were made (adaptive plan).The dose to the OARs and CTV in the initial plan was estimated by dividing the dose to the OARs and CTV multiplied by 46/70. The dose to the OARs and CTV in the adaptive plan was estimated just by dividing the dose to the OARs and CTV multiplied by 24/70.
Table 2
Mean dosimetric parameters, Conformity index and Homogeneity index of the target of the A-IMRT and A-IMPT treatment plans
| A-IMRT (Mean ± SD) | A-IMPT (Mean ± SD) | p value |
V70Gy (BODY) | 412.2 ± 149.3 | 348.8 ± 149.4 | 0.0013 |
V70Gy (CTV70) | 216.3 ± 85.2 | 218.2 ± 88.6 | 0.31 |
D5% (CTV70) | 73.9 ± 0.62 | 73.5 ± 1.43 | 0.007 |
D95% (CTV70) | 71.3 ± 0.64 | 71.6 ± 0.98 | 0.58 |
Conformity Index | 1.94 ± 0.23 | 1.59 ± 0.12 | 0.0068 |
Homogeneity Index | 1.04 ± 0.009 | 1.03 ± 0.008 | 0.073 |
Abbreviations: V70Gy, volume receiving 70 Gy; D5%, dose to the 5% of the volume; D95%, dose to the 95% of the volume.
Image Registration, Dose Accumulation, And Evaluation (Figure 1)
The way of adaptive treatment planning and dose accumulation in A-IMXT and A-IMPT has been based on the method previously described by Goras et al.18 and is illustrated in Fig. 1. We used the MIM version 6.92 (EuroMeditec, Tokyo) for the deformable image registration of the initial and adaptive plans. The CT images in the adaptive plan were deformed to be registered into the CT images in the initial plan. The CTV2 was deformed and registered into the CT used in the initial plan. The deformed CTV2 was denoted as CTV70. Dose volume statistics were carried out by adding the dose in the initial plan to the dose in adaptive plan on the CT images in the initial plan.
Statistical Comparison:
The conformity index (CI) was obtained by dividing the volume in the body which received 70 Gy or more by the volume in the CTV70 which received 70 Gy or more in A-IMXT and 70GyE or more in A-IMPT. The homogeneity index (HI) was obtained by dividing the dose which was received by 5% of CTV70, D5%, by the dose which was received by 95% of CTV70, D95%.
Dose volume statistics for the OARs were all analyzed and a comparison of A-IMXT and A-IMPT was made. The paired t-test was used to compare the mean dose (Dmean) for the OARs. We contoured the oral cavity, pharyngeal constrictor muscle (PCM), and larynx according to the Charlotte et al. method19.