Patient selection and characteristics
This study was approved by the ethics board of Shandong Cancer Hospital and Institute and informed consent has been obtained from the participants involved. Fifty-one patients with pathologically proven thoracic EC who had undergone preoperative or definitive concurrent chemoradiotherapy between May 2015 and June 2017 at Shandong Cancer Hospital and Institute were enrolled. Among the selected patients, there were seventeen cases each of upper, middle, and lower EC. One patient with lower EC was excluded due to the lack of PET-CT imaging data. All patients underwent a diagnostic imaging examination that included an endoscopy/EUS, esophagography, and FDG-PET/CT before receiving chemoradiotherapy. The average time for acquiring the diagnostic images was within the two-week period before chemoradiotherapy. Table 1 presents the patient characteristics.
Multimodal imaging
Endoscopy/EUS examination: All patients underwent diagnostic endoscopy examinations using an electronic gastroscope (Olympus GIF-Q260J) before treatment. Seven patients did not undergo EUS examinations due to esophageal stenosis. The ultrasonic probe (Olympus EVIS EUS EU-ME2) was inserted into the patient’s esophagus along the track of the biopsy forceps to detect the depth of tumor infiltration in the esophageal wall and the extent of proximal and distal tumor infiltration. The distances from the proximal and distal ends of the tumor to the incisors were recorded.
Esophagography (X-ray) image acquisition: Esophagography was performed before treatment using a digital radiography machine (Siemens Luminos dRT Max). All barium examinations were performed under fasting conditions, followed by a standard protocol (drinking 200 ml of diluted barium, in the upright, supine, and prone positions, with and without the gas powder).
PET/CT image acquisition: The PET-CT scan was performed within the two-week period prior to the planning CT scan as a part of the routine diagnostic management for EC. An 18F-FDG PET/CT scan of the chest was performed with an integrated PET/CT system (Philips Gemini TF Big Bore). The PET images were reconstructed with the CT-derived attenuation correction using an ordered subset expectation maximization algorithm with post-reconstruction Gaussian filtering, with a full width at half maximum of 5 mm.
Planning CT image acquisition: During the simulation, all patients were immobilized using a thermoplastic mask in the supine position with the arms placed along the side of the body. Each patient underwent an enhanced planning CT scan of the thoracic region on a 16-slice CT scanner (Philips Brilliance Bores CT) under free-breathing conditions. The planning CT images were reconstructed using a thickness of 3 mm and subsequently transferred to an Eclipse treatment planning system (Varian Eclipse 11).
Target volume delineation
A treatment planning system (Eclipse; Varian Medical Systems, Inc., Palo Alto, CA, USA) was used to contour the GTVs of the primary EC. The visualization parameter for delineation included the mediastinal window set to +40/400 HU. Before contouring, some clinical information such as the physical examination, pathological findings, and diagnostic CT image data were made available to the observers, while they were blind to the diagnostic endoscopy/EUS, esophagography, and FDG-PET/CT data. If the positive lymph nodes could not be separated from the primary tumor visually, they were delineated together with the primary tumor.
Five radiation oncologists (observers), who were blind to the diagnostic endoscopy/EUS, esophagography, and FDG-PET/CT patient data, were asked to independently delineate the GTVs with reference to different combinations of the multimodal images, including planning CT only (GTVC), CT combined with endoscopy/EUS (GTVCE), CT combined with endoscopy/EUS and esophagogram (X-ray) (GTVCEX), and CT combined with endoscopy/EUS, esophagogram, and FDG-PET/CT (GTVCEXP) (Figure 1). All observers were blind to the contours delineated by the other oncologists and their own former/previous contours. Observers 1 and 2 with clinical experience within five years were regarded as junior observers, while observers 3, 4, and 5 with more than ten years of clinical experience were regarded as senior observers. All contours were delineated in about two years. A delay of at least two months existed between each contouring of the tumor to eliminating a recall of the previous contouring for observers 1, 2, 3, and 5. The time interval for observer 4 was only one month, as the former observer 4 dropped out of the delineation process due to parturition.
Inter- / intra-observer variability analysis
Inter-/intra-observer variability in the volume, longitudinal length, generalized conformity index (CIgen), and position of the GTVs was assessed. The intra-observer variability can be generally regarded as the variability of the same observer when re-contouring a single case. However, in this study, it is defined as the variability of the contours on the four multimodal imaging/image combinations for one observer [23].
The mean volume and longitudinal length of the GTVs based on different multimodal imaging combinations for different observers were calculated. The inter-observer variability in the volume and longitudinal length on different multimodal imaging, combinations and the intra-observer variability for different observers were measured. The tumor length was measured using CT, endoscopy/EUS (43 cases), esophagography, and FDG-PET/CT, with the difference between the tumor length and corresponding longitudinal length of the GTVs subsequently evaluated.
CIgen was defined as the ratio of the common volume to encompassing volume [13,24]. The generalized CI (CIgen) was used to assess the overall consistency of all volume combinations delineated by different observers on the same imaging-modality combination, and that delineated by the same observer on different imaging-modality combinations. The formula is given by [13,25]:
CIgen is a good parameter for revealing the difference in the volumes delineated based on the size, shape, and location [10,23]. The use of CIgen tends to decrease the bias in the number of delineations [13]. The lower is the CIgen value for the same imaging-modality combination, the greater is the inter-observer variability. Similarly, a lower CIgen for the same observer suggests a greater intra-observer variability.
In addition, the x (right-left), y (anterior-posterior), and z (superior-inferior) axes of the center of mass (COM) of the volume were measured. The centroid shifts between the different volumes were then obtained. Finally, the three dimensional (3D) centroid shifts were calculated using the followed equation [24,26]:
Statistical analysis
Statistical analysis was performed using the SPSS software package (SPSS 25.0). All the data had an approximately normal distribution. The one-way ANOVA test was applied to detect the inter-/intra-observer variability in the volume, longitudinal length, CIgen, and position of the GTVs among different observers and different multimodal imaging combinations. The paired t-test was used to compare the volume, longitudinal length, CIgen, and position of the GTVs between two observers or two multimodal imaging combinations. A P<0.05 was considered significant.