Skip to main content

Imaging biomarkers of contrast-enhanced computed tomography predict survival in oesophageal cancer after definitive concurrent chemoradiotherapy



This study aimed to evaluate the predictive potential of contrast-enhanced computed tomography (CT)-based imaging biomarkers (IBMs) for the treatment outcomes of patients with oesophageal squamous cell carcinoma (OSCC) after definitive concurrent chemoradiotherapy (CCRT).


Altogether, 154 patients with OSCC who underwent definitive CCRT were included in this retrospective study. All patients were randomised to the training cohort (n = 99) or the validation cohort (n = 55). Pre-treatment contrast-enhanced CT scans were obtained for all patients and used for the extraction of IBMs. An IBM score, was constructed by using the least absolute shrinkage and selection operator with Cox regression analysis, which was equal to the log-partial hazard of the Cox model in the training cohort and tested in the validation cohort. IBM nomograms were built based on IBM scores for individualised survival estimation. Finally, a decision curve analysis was performed to estimate the clinical usefulness of the nomograms.


Altogether, 96 IBMs were extracted from each contrast-enhanced CT scan. IBM scores were constructed from 11 CT-based IBMs for overall survival (OS) and 8 IBMs for progression-free survival (PFS), using the LASSO-Cox regression method in the training cohort. Multivariate analysis revealed that IBM score was an independent prognostic factor correlated with OS and PFS. In the training cohort, the C-indices of IBM scores were 0.734 (95% CI 0.664–0.804) and 0.658 (95% CI 0.587–0.729) for OS and PFS, respectively. In the validation cohort, C-indices were 0.672 (95% CI 0.578–0.766) and 0.666 (95% CI 0.574–0.758) for OS and PFS, respectively. Kaplan–Meier survival analysis showed a significant difference between risk subgroups in the training and validation cohorts. Decision curve analysis confirmed the clinical usefulness of the IBM score.


The IBM score based on pre-treatment contrast-enhanced CT could predict the OS and PFS for patients with OSCC after definitive CCRT. Further multicentre studies with larger sample sizes are warranted.


Oesophageal cancer (OC) is one of the most common cancers globally and in 2018 its incidence and number of cancer-related deaths ranked seventh and sixth, respectively [1]. The management of OC typically involves multidisciplinary therapy, including definitive concurrent chemoradiotherapy (CCRT), which is the main standard treatment for oesophageal squamous cell carcinoma (OSCC) for medically unresectable tumours; and is also an option for resectable tumours. However, the outcomes of CCRT among these patients are still disappointing, with 3-year overall survival (OS) rates of 23–44.7% [2,3,4,5]. More than 50% of patients in the RTOG 85-01 trial and INT 0123 trial experienced locoregional disease progression [3, 5]. Patients with higher mortality risk following CCRT may benefit from more intensive primary treatment (e.g., planned radical surgery after CCRT), adjuvant therapy (e.g. chemotherapy), or more frequent follow-up. The application of these strategies requires the identification of patients with high mortality risk prospectively to achieve personalised management. Thus, to improve the overall survival of patients with OC after CCRT, it is crucial to predict the mortality risk of each individual patient.

Prediction of outcomes among patients with OC after CCRT remains an unmet clinical need. One of the most commonly used methodologies for prognostic evaluation in the clinic is the TNM staging system, which stratifies patients into different stages according to their tumour burden. Although the clinical stage system provides important insights for evaluating outcomes of patients with different stages, its role in survival prediction among patients with the same disease stage is non-significant. Indeed, previous studies have shown that the clinical stage system fails to predict heterogeneous outcomes of patients with locally advanced disease following CCRT [6,7,8]. A variety of other clinical factors and biomarkers have also been assessed for their prognostic potential [9,10,11]. Yet none of these factors have been widely used for the clinical stratification of patients and decision-making, as each of them has weaknesses and limitations.

By using current imaging techniques, quantitative imaging biomarkers (IBMs) could become an interesting way of assessing multiple cancer diagnosis and prognosis. These so-called radiomics could extract relevant information from commonly available images with a high throughput [12,13,14]. Previous studies have reported on the potential information to be gleaned from computed tomography (CT) IBMs in OC, and were able to assess and predict histopathological characteristics, treatment response, or survival outcome among patients with OC to some extent [15,16,17,18,19]. CT scans play an important role in the radiation treatment of OC, including diagnosis, staging, treatment planning, quality control, and follow-up. Non-contrast enhanced CT-based IBMs have been shown to be correlated with patient outcomes for a number of cancer types, including OSCC [18, 20]. However, the most commonly available imaging modalities for patients after undergoing definitive CCRT were not non-contrast enhanced CT scans but contrast-enhanced CT scans, which were performed during treatment planning. A previous study suggested that post-treatment IBMs extracted from contrast-enhanced CT images might have a correlation with OS in patients with OC who received definitive CCRT [16]. Although the sample size of this study was small and only included 26 cases of squamous cell carcinoma (SCC), it provided encouragement for the pursuit of further studies. The most common pathological type of OC in China is SCC, and radiotherapy is administered with a total dose of 60–66 Gy [21]. This range is much higher than the dose used in the standard treatment of OC via conventional fractionated radiotherapy (50.4 Gy). Oesophageal oedema is a common acute adverse event after definitive CCRT. There are fewer residual lesions that could be used for objective analysis or evaluation after definitive CCRT. Further, some patients with a complete response did not have residual lesions. It is unclear whether contrast-enhanced CT images obtained before treatment could serve as a feasible source for radiomics analysis in OSCC. Therefore, stronger evidence is needed in support of the implications for survival outcomes and the reliability of the methodology.

In this study, based on pre-treatment contrast-enhanced CT images, we sought to develop and validate an IBM score to predict OS and progression-free survival (PFS) for patients with OSCC and assess its value for individual OS and PFS estimation.



The protocol for this retrospective study was obtained from the local ethics and institutional review board. Approval and the need for informed consent had been waived. This study included patients with OC who underwent definitive CCRT at AAA between September 2009 and August 2015. The inclusion criteria were: (1) pathological diagnosis of OSCC; (2) primary tumour located in the cervical, upper thoracic, or middle thoracic oesophagus; and (3) contrast-enhanced CT scan findings, which were used in treatment planning before definitive CCRT. The exclusion criteria were: (1) patients who only received radiotherapy or chemotherapy; (2) prior surgery or administration of chest radiotherapy or chemotherapy. As shown in Fig. 1, the final study population consisted of 154 patients. All patients received intensity modulated radiation therapy (IMRT) combined with chemotherapy. Of these, 78 patients with OSCC were from a phase II prospective clinical study, using simultaneous modulated accelerated radiotherapy (SMART) combined with chemotherapy [22]. The 154 patients were randomly assigned into a training cohort (n = 99) and validation cohort (n = 55).

Fig. 1

Flowchart of inclusion in the present study. OSCC, oesophageal squamous cell carcinoma; CCRT, concurrent chemoradiotherapy; CT, computed tomography

All patients underwent simulation CT scans for treatment planning. Seventy-eight patients underwent SMART, followed by radiation therapy with a prescribed dose of 66 Gy/30F, 5 days/week. Other patients underwent radiation therapy with a prescribed dose of 64/32F, 5 days/week. Most patients (90.9%) received concurrent chemotherapy based on the cisplatin and 5-fluorouracil (PF) regimen. The intensity of concurrent chemotherapy was relatively reduced among patients of advanced age or with poor performance status. Data regarding clinical characteristics of patients were collected in both cohorts, including age, sex, clinical stage, and tumour location. Dose-volume information for the primary tumour was collected from the radiotherapy planning system. Further details are shown in Table 1.

Table 1 Clinical characteristics of 154 patients with OSCC after definitive CCRT

Contrast-enhanced CT image acquisition

The CT scans of all patients were acquired (Philips Brilliance CT Big Bore Oncology Configuration, Cleveland, OH, USA; voxel size: 1.0 × 1.0 × 3.0 mm3 for 79 patients and 1.0 × 1.0 × 5.0 mm3 for 72 patients; convolution kernel: Philips Healthcare’s B), using a scanning voltage of 120 kVp with a slice thickness of 3–5 mm after an intravenous injection of 75 ml of 300 mg/mL iodinated contrast agent at a rate of 1.8–2 mL/sec with a pump injector (Medrad Stellant; Bayer, Beijing, China). The CT images were transmitted to the radiation therapy planning system (Eclipse Planning System version 10.0) via the DICOM 3.0 port.

Region of interest (ROI) delineation and IBMs extracted

Pre-treatment contrast-enhanced CT scan images of patients were exported for analysis. The primary tumour was delineated by experienced radiation oncologists on the mediastinal window of the planning CT scan. IBMs were extracted by internal programming software using MATLAB R2016a (Mathworks, Natick, USA) and its toolbox. From the contrast-enhanced CT images of each patient, 96 IBMs were extracted, including the following types: (1) 24 CT intensity IBMs, describing the distribution of voxel parameter values in the volume of interest, such as the min, max and skewness of the primary tumour intensity; (2) 20 geometric IBMs that calculated the size and shape of the volume of interest, such as sphericity, volume, surface and long axis length; and (3) 52 texture IBMs, that described the difference in voxel density distribution of the three-dimensional contoured structure and consisted of four different matrices: grey level co-occurrence (GLCM) [23], grey level run-length (GLRLM) [24], neighbourhood grey-tone difference (NGTDM) [25], and grey level size-zone (GLSZM) matrices [26]. More details on the algorithms for IBM extraction and application have been discussed in previous studies [14, 27].

Pre-selection Method and IBM score building

Because high correlations between most of the IBM variables were expected, in order to reduce the statistical probability of multi-collinearity, three rules were implemented to pre-select IBM variables for further analysis. First, IBM variables were assessed in the univariable analysis; variables with a p value less than 0.25 were used for the next analysis. Second, from highly correlated pairs of IBMs (i.e. the Pearson correlation coefficient r ≥ 0.8) variables with the higher p value in the Cox univariable analysis were omitted. Third, we performed the least absolute shrinkage and selection operator (LASSO) for the Cox regression model to select the most useful prognostic IBM variables from the potential predictors.[28]. The multiple-IBM-based scores (defined as the IBM scores), which were equal to the log-partial hazard of the Cox model, were calculated for each patient to reflect the risk of mortality or tumour progression and variance inflation factor (VIF) used to evaluate the collinearity among these final IBMs.

IBM score performance and validation

As patients with OSCC were assigned into two cohorts, the performances of the IBM score were evaluated by the concordance indices (C-indices), respectively. The potential correlation of the IBM score with the OS and PFS for both the training and validation cohorts was assessed by using Kaplan–Meier survival curve analyses. Time-dependent receiver operating characteristic (ROC) curves were plotted for both the training and validation cohorts in term of OS and PFS. 95% confidence intervals were used as the confidence level on the ROC curves in this study. The optimal cut-off values of the ROC curves were determined using the Youden Indices (YIs) in the training cohort and patients were divided into high- and low IBM score subgroups; the thresholds of which were stratified by the maximum YIs. The same cut-off values were then applied to the validation cohort. Multivariable Cox proportional hazards analysis was used to assess the IBM score as an independent predictor by integrating clinical risk factors. In the training cohort, nomograms based on the IBM score were developed to assess individual patient-level probability estimates for the median survival time and 1-year, 3-year, and 5-year OS or PFS rates according to each patient’s unique combination of baseline characteristics. To estimate the clinical utility of the IBM nomograms, decision curve analysis (DCA) was used to quantify the net benefits at different threshold probabilities in both cohorts.


The survival estimates mainly assessed in this study were OS and PFS. OS was defined as the time from the beginning of radiation therapy to death due to any cause or the last day of clinical follow-up, while PFS was defined as the time from the beginning of radiation therapy to first relapse at any site or death from any cause, whichever occurred first, or the last day of clinical follow-up.

Statistical analysis

The clinical features of the patients in the two cohorts were compared using the independent t-test or chi-squared test, with a statistical significance level of 0.05 for 2-tailed test. All statistical analyses were performed using R version 3.6.0 (The R Foundation for Statistical Computing, Vienna, Austria) and SPSS version 23.0 (IBM Corp, Armonk, NY, USA). The LASSO algorithm was implemented using the glmnet package in the R environment [29]. The ROC and Kaplan–Meier curves were plotted using the pROC and survminer packages, respectively, in the R environment. Nomograms were constructed using the rms and survival packages in the R environment. The DCA curves were created using the rmda package in the R environment.


Baseline clinical results

The clinical factors for the training and validation cohorts are listed in Table 1. Age was tested by independent t-test, and other clinical factors were tested by chi-squared test. No significant differences in patients’ clinical characteristics were found between the two cohorts. Of the 154 patients included in the study, 119 (77.3%) were men, and the median (interquartile range, IQR) age of all patients was 60 years (55–64.25 years). In the training cohort, the median (IQR) survival time for OS and PFS were 43 (19–59) and 36 (10–55) months, respectively. In the validation cohort, the median (IQR) survival time for OS and PFS were 37 (22–55) and 36 (12–54) months, respectively.

IBM selection results

In the univariable analysis, 46 IBM variables were used and 18 IBMs remained after comparing the inter-variable correlations (Additional file 1: Table S1). 11 potential predictors with non-zero coefficients were selected in the LASSO Cox regression model. We plotted the partial likelihood deviance versus log (λ), where λ is the tuning parameter (Fig. 2). A dotted vertical line was drawn at log (λ) = –2.643, which corresponded to the best value λ = 0.071. The optimal tuning parameter resulted in 11 non-zero coefficients. With their corresponding coefficients in the LASSO-Cox model, the calculation formulas of IBM score for OS (Formula 1) was constructed as:

$$\begin{aligned} & {\text{IBM}}\,{\text{score}} = 0.019 \times {\text{Range}} - 3.231 \times {\text{Q}}975 - 0.260 \times {\text{Sphericity}} \\ & \quad + 0.093 \times {\text{Major\_axis}}\_{\text{length}} + 0.035 \times {\text{Maximum}}\_{\text{Probability}}\_{\text{GLCM}} \\ & \quad - 0.304 \, \times {\text{Sum}}\_{\text{of}}\_{\text{Square}}\_{\text{Variance}}\_{\text{GLCM}} - 0.014 \times {\text{Coarseness}}\_{\text{NGTDM}} \\ & \quad - 0.020 \times {\text{Contrast}}\_{\text{NGTDM}} + 0.020 \times {\text{Busyness}}\_{\text{NGTDM}} + 0.670 \\ & \quad \times {\text{Small}}\_{\text{Zone}}\_{\text{Emphasis}}\_{\text{GLSZM}}- 0.064 \times {\text{Zone}}\_{\text{percentage}}\_{\text{GLSZM}} + 2.5. \\ \end{aligned}$$
Fig. 2

IBM selection using the LASSO-Cox regression model. a Ten-fold cross-validation was performed to select IBMs using LASSO method. At the optimal values of the tuning parameter (λ), the partial likelihood deviance curve was plotted versus log (λ). The dotted vertical line was set with minimum criteria, a λ value of 0.071127 with log (λ) of − 2.643291, where 11 IBMs were selected. b The coefficient profiles of 18 IBMs in the LASSO. The coefficient profile plot was generated against the log (λ) sequence. With ten-fold cross-validation, the dotted vertical line showed the non-zero coefficients selected, where eleven IBMs were selected. LASSO, least absolute shrinkage and selection operator; IBMs, imaging biomarkers

The constant value was 2.5, which was used to obtain IBM scores > 0 from the calculation formula. The VIFs of the 11 IBMs were acceptable, ranging from 1.150–3.403, indicating no collinearity problems (Additional file 1: Table S2).

The same analysis was used to select the IBMs which were associated with PFS in the training cohort. 32 IBM variables were preselected in the univariable analysis and 12 IBMs remained after comparing the inter-variable correlations (Additional file 1: Table S1). A dotted vertical line was drawn at log (λ) = − 2.702, which corresponded to the best value λ = 0.067 (Additional file 2: Figure S1a and Figure S1b). The optimal tuning parameter resulted in 8 non-zero coefficients. With their corresponding coefficients in the LASSO Cox model, the calculation formulas of IBM score for PFS (Formula 2) was constructed as:

$$\begin{aligned} & {\text{IBM}}\,{\text{score}} = - 0.663 \times {\text{Q}}75 - 0.234 \times {\text{Q}}975 - 0.246 \\ & \quad \times {\text{Volume}}\_{\text{Density}}- 0.01 \times {\text{Sphericity}} + 0.122 \\ & \quad \times {\text{Major}}\_{\text{axis}}\_{\text{length}} - 0.021 \times {\text{Contrast}}\_{\text{NGTDM}} \\ & \quad + 0.468\times {\text{Small}}\_{\text{Zone}}\_{\text{Emphasis}}\_{\text{GLSZM}} - 0.069\\& \quad \times {\text{Zone}}\_{\text{percentage}}\_{\text{GLSZM}} - 1.4. \end{aligned}$$

The VIFs of the 8 IBMs were acceptable, ranging from 1.266 to 4.524, indicating no collinearity problems (Additional file 1: Table S2).

Performance of IBM score

In the training cohort, we evaluated the predictive accuracy of the IBM score using ROCs analysis at different time points of follow-up. As shown in Fig. 3a, b, the area under the curves (AUCs) of the IBM score were 0.856 (95% CI 0.756–0.955, p < 0.001) and 0.779 (95% CI 0.663–0.895, p < 0.001) in terms of 5-year OS and 5-year PFS, respectively; and the optimal cut-off values of YI were 1.012 and 0.688, respectively. According to the maximum YI, the optimal cut-off values generated by ROC curves were 1.012 for 5-year OS and 0.688 for 5-year PFS. Patients were then stratified into high-risk or low-risk subgroups. In the training cohort, the 5-year OS and 5-year PFS were 85.0% and 74.1% respectively for the low-risk subgroup and 35.9% and 35.2% respectively for the high-risk subgroup (hazard ratios [HRs]:6.003 (95% CI 2.646–13.618) and 3.416 (95% CI 1.698–6.873), respectively; all p < 0.001, log-rank test; Fig. 4a, b). We then tested the same analyses using the ROC and Kaplan–Meier analysis, and similar results were observed in the validation cohort. As shown in Fig. 3c, d, the AUCs of the IBM score were 0.867 (95% CI 0.726–1.000, p = 0.001) and 0.852 (95% CI 0.713–0.990, p = 0.002) for the 5-year OS and 5-year PFS, respectively. Patients were then stratified into high-risk or low-risk subgroups. In the validation cohort the 5-year OS and 5-year PFS were 67.9% and 66.0% respectively for the low-risk subgroup; and 30.8% and 35.9% respectively for the high-risk subgroup (HR 2.957 (95% CI 1.104–7.919) and 3.051 (95% CI 1.324–7.034), respectively; all p < 0.05; Fig. 4c, d). Patients with OSCC with lower IBM scores were more likely to obtain a survival benefit from definitive CCRT. Those with high IBM scores had significantly poorer OS and PFS according to univariable Cox regression analysis (Additional file 1: Table S3). Multivariable Cox regression analysis for clinical factors and IBM score also revealed that the IBM score remained a powerful and independent predictive factor for OS and PFS in both training and validation cohorts (Table 2).

Fig. 3

Time-dependent ROC curves for the performance of IBM scores in the training cohort and validation cohort. a Time-dependent OS ROC for training cohort; b Time-dependent PFS ROC for training cohort; c Time-dependent OS ROC for validation cohort; d Time-dependent PFS ROC for validation cohort. AUCs for 1, 3, 5 years were used for survival prediction. ROC, receiver operation characteristic; AUC, area under the curve; CI, confidence interval; OS, overall survival; PFS, progression-free survival; IBM, image biomarker

Fig. 4

Kaplan–Meier survival analysis of overall survival and progression-free survival according to the optimum cut-offs of IBM scores. a Kaplan–Meier survival analysis of OS for the training cohort; b Kaplan–Meier survival analysis of PFS for the training cohort; c Kaplan–Meier survival analysis of OS for the validation cohort; d Kaplan–Meier survival analysis of PFS for validation cohort. We calculated the p value using the log-rank test. HR, hazard ratio; CI, confidence interval; IBM, imaging biomarker; OS, overall survival; PFS, progression-free survival

Table 2 Multivariable association of IBM score, clinical factors with OS and PFS in the training and validation cohort (likelihood Ratio: Backward stepwise)

Clinical benefit of IBM score

Clinical stage was associated with OS and PFS when using the Cox univariable analysis; however, it was not identified as a predictive independent factor for OS or PFS using multivariable analysis in the training cohort. IBM nomograms using only IBM scores for OS and PFS were constructed (Fig. 5a, b). In the training cohort, the C-indices of the models were 0.734 (95% CI 0.664–0.804) and 0.658 (95% CI 0.587–0.729) for OS and PFS, respectively. Similar results were observed in the validation cohort; the C-indices were 0.672 (95% CI 0.578–0.766) and 0.666 (95% CI 0.574–0.758) for OS and PFS, respectively. The C-index values showed that the IBM nomograms had good prognostic performance in both training and validation cohorts.

Fig. 5

IBM nomograms for OS (a) and PFS (b). The constructed IBM nomograms were used to estimate OS and PFS for individual OSCC patients. IBM, imaging biomarker; OS, overall survival; PFS, progression-free survival; GLCM, grey level co-occurrence matrices; NGTDM, neighbourhood grey-tone difference matrices; GLSZM, grey level size-zone matrices

The decision curve analysis showed that IBM score had higher overall net benefits than clinical stage, within a major range of reasonable threshold probability (Fig. 6). Compared to clinical stage, the IBM score demonstrated better discrimination capability in both training and validation cohorts.

Fig. 6

Decision curve analysis of PFS and OS were compared between the IBM score and clinical stage in the training and validation cohort, respectively. a DCA of OS for the training cohort; b DCA of PFS for the training cohort; c DCA of OS for the validation cohort; d DCA of PFS for the validation cohort. The y-axis represents the net benefit. The x-axis represents the threshold probability. Across the full range of threshold probabilities, the horizontal black line indicated that no patient chose to undergo follow-up, and the green line indicated that all patients underwent follow-up. The red line represented IBM score. The blue line represented clinical stage. Compared to the clinical stage, the IBM score had the higher net benefit. DCA, decision curve analysis; IBM, imaging biomarker; OS, overall survival; PFS, progression-free survival


This study showed that IBMs from contrast-enhanced CT images might allow prediction of OS and PFS for OSCC patients. The IBM score was revealed to be an independent prognostic factor for OSCC patients. Patients were successfully stratified into low-risk and high-risk subgroups by the IBM score, with significant differences in OS and PFS. IBM nomograms demonstrated better discrimination capability compared to traditional clinical stage, indicating the clinical value of the IBM score for individualised OS and PFS estimation.

The new IBM scores demonstrated significant associations with the OS and PFS of patients with OSCC. For geometric IBMs, sphericity, volume-density and the major axis length quantified the sphericity and size of tumours. Previous studies revealed that these IBMs basically represented tumour volume, which were significantly associated with treatment outcomes [30, 31]. In our study, the discrimination performances of the IBM nomograms were decreased when volume-related IBMs were omitted from the IBM score (C-index for the radiomics nomogram: OS, 0.672 (95% CI, 0.588–0.757); PFS, 0.629 (95% CI, 0.545–0.713) in the training cohort). These volume-related IBMs can promote the objective evaluation of subtle changes within tumours and provide clues to lesion invasiveness and growth-patterns [30, 32]. Range and Q975 were obtained from the histogram of voxel intensities and represented the heterogeneity of voxel intensities within the ROI [27].

The higher value texture IBMs, including maximum probability and sum of square variance, indicated the greater distribution variability of grey-level intensity values in the image [33, 34]. Coarseness, contrast, and busyness were all textural IBMs derived from NGTDM. Coarseness was used to quantify the granularity of the VOI of the tumour. Our study showed that higher value of busyness, lower values of coarseness or of contrast, might all be associated with poorer OS. The predictive and prognostic value of these pre-treatment IBMs had been previously demonstrated in several types of cancer [35,36,37]. Tixier et al. reported that coarseness from NGTDM was a strong predictor of treatment response for patients with OC following definitive CCRT [35]. However, IBMs derived from NGTDM were influenced by the reconstruction settings; therefore multicentre trials are still needed to standardise these IBMs [38]. Small zone emphasis measures the distribution of small size zones and small dependencies, while zone percentage assesses the distribution of large zones of the same intensity, and not of small groups of pixels or segments in any given direction [26, 35]. These texture IBMs containing spatial information among voxels could strongly reflect intra-tumour heterogeneity which was highly relevant to poor prognosis [12]. In order to correlate the multiple IBMs with the pathophysiological basis of tumours in an intuitive method, we constructed the multi-feature IBM score, which provided novel oncological biomarkers for obtaining phenotypic information, potentially assisting clinicians in formulating management strategies.

Current guidelines recommend definitive CCRT as a standard component for locally advanced OSCC therapies. However, several studies suggested that certain subgroups of patient failed to benefit from the present definitive CCRT strategies [6, 9]. Therefore, accurately distinguishing the risk subgroups of OSCC patients will help improve the current prognostic system and guide towards more personalised treatment. A few studies have focused on the correlations between radiomics analysis and treatment outcomes evaluation. Zhai et al. [30]. found that heterogeneous IBMs on CT images were significantly correlated with OS and helped improve the performance of clinical factors for OS among head and neck cancer patients. Mule et al. [39] investigated contrast-enhanced CT outcomes that might help predict survival in patients with advanced hepatocellular carcinoma treated with sorafenib. In the present study, our findings indicated that patients with OSCC with higher IBM scores had a greater likelihood of worse survival rates and failure to respond to CCRT. High-risk patients with OSCC identified in the present studies may benefit from more effective approaches to improve survival outcomes [40, 41]. Thus, the IBM score may serve as a prognostic tool for OSCC patients after definitive CCRT.

TNM staging is the most useful tool to stratify OSCC patients into different stages according to their tumour burden. However, its role in survival prediction among OSCC patients with the same clinical stage was non-significant. To develop an individualised easy-to-use tool for clinicians, we attempted to construct nomograms based on the IBM score to predict the prognosis of individual patients. These IBM nomograms could be used to predict the median survival time, and the probability of 1-year, 3-year and 5-year OS and PFS for individual OSCC patients. The nomograms performed well, with significant C-indices, and demonstrated good discrimination and clinical utility in both the training and validation cohorts. The decision curve analysis indicated that the IBM score was superior to the clinical stage, within a major range of reasonable threshold probability. Notably, time-dependent ROC curves showed that the IBM score did not have a good predictive performance for survival within 1 year. It was unclear why discrepancies remained for the training and validation cohorts. One possible explanation was the small sample size, retrospective nature of our study, and model-fitting differences. Further analysis for OSCC patients is needed to establish this.

For OC patients, contrast-enhanced CT scan is the main imaging procedure performed in conventional clinical practice [42]. It has been reported that IBMs extracted from contrast-enhanced CT images might be correlated with the spatial variability in microvessel density [43]. However, in standard CT images, IBMs might be associated with variability in tissue densities due to spatially variable fibrosis, cell density, and necrosis [13]. Badic at el. suggested that IBMs extracted from standard CT and contrast-enhanced CT images could provide complementary prognostic information from both approaches [44]. In view of the wide availability of contrast-enhanced CT scans among patients undergoing definitive radiotherapy, our study provides an important basis for conducting large-scale and multicentre research. It is important to note that quality assurance of contrast-enhanced CT scans will have a critical impact on radiomics based on these images. Furthermore, verification is needed on whether IBMs extracted from contrast-enhanced CT images could provide prognostic information for patients with oesophageal adenocarcinoma.

The limitations of our retrospective design include several aspects that were insufficient for the model [45]. This was a retrospective, single-centre study, involving a relatively small sample size. This could be addressed more thoroughly in future by using a larger sample size with multicentre validation cohorts to acquire high-level evidence for survival outcomes. Compared to the IBM score, the clinical factors used in this study demonstrated poor discrimination ability in predicting OS and PFS; other potential prognostic biomarkers could be incorporated into our IBM nomograms. A combination of multiple biomarkers and IBMs may improve the capability of predicting OS and PFS among patients with OSCC undergoing definitive CCRT.


We demonstrated that IBMs extracted from contrast-enhanced CT images could effectively predict survival among OSCC patients. The IBM score might serve as a non-invasive predictive tool to guide individualised treatment decisions. Further studies with a larger sample size and multicentre validation are required.

Availability of data and materials

All data supporting the conclusions of this article is available upon request from the corresponding author.



Oesophageal squamous cell carcinoma;


Computed tomography


Imaging biomarkers


Concurrent chemoradiotherapy


American Joint Committee on Cancer (AJCC) staging system (version 6.0th)


Squamous cell carcinoma


Cisplatin and 5-fluorouracil


Intensity modulated radiation therapy


Simultaneous modulated accelerated radiotherapy


Grey level co-occurrence matrices


Grey level run-length matrices


Neighbourhood grey-tone difference matrices


Grey level size-zone matrices


Least absolute shrinkage and selection operator


Variance inflation factor


Youden Index


Decision curve analysis


Receiver operation characteristic


Area under the curve


Confidence interval


Overall survival


Progression-free survival


  1. 1.

    Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA A Cancer J Clin. 2018;68(6):394–424.

    Article  Google Scholar 

  2. 2.

    Kato K, Muro K, Minashi K, Ohtsu A, Ishikura S, Boku N, Takiuchi H, Komatsu Y, Miyata Y, Fukuda H. Phase II study of chemoradiotherapy with 5-fluorouracil and cisplatin for Stage II-III esophageal squamous cell carcinoma: JCOG trial (JCOG 9906). Int J Radiat Oncol Biol Phys. 2011;81(3):684–90.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  3. 3.

    Minsky BD, Pajak TF, Ginsberg RJ, Pisansky TM, Martenson J, Komaki R, Okawara G, Rosenthal SA, Kelsen DP. INT 0123 (Radiation Therapy Oncology Group 94–05) phase III trial of combined-modality therapy for esophageal cancer: high-dose versus standard-dose radiation therapy. J Clin Oncol. 2002;20(5):1167–74.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  4. 4.

    Ohtsu A, Boku N, Muro K, Chin K, Muto M, Yoshida S, Satake M, Ishikura S, Ogino T, Miyata Y, et al. Definitive chemoradiotherapy for T4 and/or M1 lymph node squamous cell carcinoma of the esophagus. J Clin Oncol. 1999;17(9):2915–21.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  5. 5.

    Cooper JS, Guo MD, Herskovic A, Macdonald JS, Martenson JA Jr, Al-Sarraf M, Byhardt R, Russell AH, Beitler JJ, Spencer S, et al. Chemoradiotherapy of locally advanced esophageal cancer: long-term follow-up of a prospective randomized trial (RTOG 85–01). Radiation Therapy Oncology Group. JAMA. 1999;281(17):1623–7.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  6. 6.

    Gertler R, Stein HJ, Langer R, Nettelmann M, Schuster T, Hoefler H, Siewert JR, Feith M. Long-term outcome of 2920 patients with cancers of the esophagus and esophagogastric junction: evaluation of the New Union Internationale Contre le Cancer/American Joint Cancer Committee staging system. Ann Surg. 2011;253(4):689–98.

    PubMed  Article  PubMed Central  Google Scholar 

  7. 7.

    Hamai Y, Hihara J, Emi M, Furukawa T, Ibuki Y, Yamakita I, Kurokawa T, Okada M. Effects of neoadjuvant chemoradiotherapy on pathological TNM stage and their prognostic significance for surgically-treated esophageal squamous cell carcinoma. Anticancer Res. 2017;37(10):5639–46.

    CAS  PubMed  PubMed Central  Google Scholar 

  8. 8.

    Oweira H, Schmidt J, Mehrabi A, Kulaksiz H, Schneider P, Schob O, Giryes A, Abdel-Rahman O. Validation of the eighth clinical American Joint Committee on Cancer stage grouping for esophageal cancer. Future Oncol (London, England). 2018;14(1):65–75.

    CAS  Article  Google Scholar 

  9. 9.

    Chen CZ, Chen JZ, Li DR, Lin ZX, Zhou MZ, Li DS, Chen ZJ. Long-term outcomes and prognostic factors for patients with esophageal cancer following radiotherapy. World J Gastroenterol. 2013;19(10):1639–44.

    PubMed  PubMed Central  Article  Google Scholar 

  10. 10.

    Chen M, Huang J, Zhu Z, Zhang J, Li K. Systematic review and meta-analysis of tumor biomarkers in predicting prognosis in esophageal cancer. BMC Cancer. 2013;13:539.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  11. 11.

    Qiao Y, Chen C, Yue J, Yu Z. Tumor marker index based on preoperative SCC and CYFRA 21–1 is a significant prognostic factor for patients with resectable esophageal squamous cell carcinoma. Cancer Biomark Sect A Disease Mark. 2019;25(3):243–50.

    CAS  Google Scholar 

  12. 12.

    Lambin P, Rios-Velazquez E, Leijenaar R, Carvalho S, van Stiphout RG, Granton P, Zegers CM, Gillies R, Boellard R, Dekker A, et al. Radiomics: extracting more information from medical images using advanced feature analysis. Eur J Cancer (Oxford, England: 1990). 2012;48(4):441–6.

    Article  Google Scholar 

  13. 13.

    Aerts HJ, Velazquez ER, Leijenaar RT, Parmar C, Grossmann P, Carvalho S, Bussink J, Monshouwer R, Haibe-Kains B, Rietveld D, et al. Decoding tumour phenotype by noninvasive imaging using a quantitative radiomics approach. Nat Commun. 2014;5:4006.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  14. 14.

    Gillies RJ, Kinahan PE, Hricak H. Radiomics: images are more than pictures, they are data. Radiology. 2016;278(2):563–77.

    PubMed  Article  PubMed Central  Google Scholar 

  15. 15.

    Ganeshan B, Skogen K, Pressney I, Coutroubis D, Miles K. Tumour heterogeneity in oesophageal cancer assessed by CT texture analysis: preliminary evidence of an association with tumour metabolism, stage, and survival. Clin Radiol. 2012;67(2):157–64.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  16. 16.

    Yip C, Landau D, Kozarski R, Ganeshan B, Thomas R, Michaelidou A, Goh V. Primary esophageal cancer: heterogeneity as potential prognostic biomarker in patients treated with definitive chemotherapy and radiation therapy. Radiology. 2014;270(1):141–8.

    PubMed  Article  PubMed Central  Google Scholar 

  17. 17.

    Hou Z, Ren W, Li S, Liu J, Sun Y, Yan J, Wan S. Radiomic analysis in contrast-enhanced CT: predict treatment response to chemoradiotherapy in esophageal carcinoma. Oncotarget. 2017;8(61):104444–54.

    PubMed  PubMed Central  Article  Google Scholar 

  18. 18.

    Larue R, Klaassen R, Jochems A, Leijenaar RTH, Hulshof M, van Berge Henegouwen MI, Schreurs WMJ, Sosef MN, van Elmpt W, van Laarhoven HWM, et al. Pre-treatment CT radiomics to predict 3-year overall survival following chemoradiotherapy of esophageal cancer. Acta oncologica (Stockholm, Sweden). 2018;57(11):1475–81.

    CAS  Article  Google Scholar 

  19. 19.

    Cheng L, Wu L, Chen S, Ye W, Liu Z, Liang C. CT-based radiomics analysis for evaluating the differentiation degree of esophageal squamous carcinoma. Zhong nan da xue xue bao Yi xue ban = J Cent South Univ Medi Sci. 2019;44(3):251–6.

    Google Scholar 

  20. 20.

    Xie C, Yang P, Zhang X, Xu L, Wang X, Li X, Zhang L, Xie R, Yang L, Jing Z, et al. Sub-region based radiomics analysis for survival prediction in oesophageal tumours treated by definitive concurrent chemoradiotherapy. EBioMedicine. 2019;44:289–97.

    PubMed  PubMed Central  Article  Google Scholar 

  21. 21.

    Ge X, Yang X, Lu X, Wen W, Zhen F, Ye H, Zhu H, Cao Y, Zhang S, Cheng H, et al. Long-term clinical outcome of intensity-modulated radiation therapy for locally advanced esophageal squamous cell carcinoma. Tumori. 2015;101(2):168–73.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  22. 22.

    Chen J, Guo H, Zhai T, Chang D, Chen Z, Huang R, Zhang W, Lin K, Guo L, Zhou M, et al. Radiation dose escalation by simultaneous modulated accelerated radiotherapy combined with chemotherapy for esophageal cancer: a phase II study. Oncotarget. 2016;7(16):22711–9.

    PubMed  PubMed Central  Article  Google Scholar 

  23. 23.

    Haralick RM, Shanmugam K, Dinstein I. Textural features for image classification. IEEE Trans Syst Man Cybernet. 1973;SMC-3(6):610–21.

    Article  Google Scholar 

  24. 24.

    Galloway M. Texture analysis using gray level run lengths. Comput Graph Image Proc. 1974;4:172–99.

    Article  Google Scholar 

  25. 25.

    Sun C, Wee WG. Neighboring gray level dependence matrix for texture classification. Comput Vis Graph Image Process. 1983;23(3):341–52.

    Article  Google Scholar 

  26. 26.

    Thibault G, Fertil B, Navarro C, Pereira S, Cau P, Levy N, Sequeira J, Mari JL. Texture indexes and gray level size zone matrix application to cell nuclei classification. In: 10th international conference on pattern recognition and information processing. 2009.

  27. 27.

    Lambin P, Leijenaar RTH, Deist TM, Peerlings J, de Jong EEC, van Timmeren J, Sanduleanu S, Larue R, Even AJG, Jochems A, et al. Radiomics: the bridge between medical imaging and personalized medicine. Nat Rev Clin Oncol. 2017;14(12):749–62.

    PubMed  Article  PubMed Central  Google Scholar 

  28. 28.

    Hall M. Correlation-based feature selection for discrete and numeric class machine learning. In: Proceedings of the seventeenth international conference on machine learning. 2000. pp. 359–366.

  29. 29.

    Friedman J, Hastie T, Tibshirani R. Regularization paths for generalized linear models via coordinate descent. J Stat Softw. 2010;33(1):1–22.

    PubMed  PubMed Central  Article  Google Scholar 

  30. 30.

    Zhai TT, van Dijk LV, Huang BT, Lin ZX, Ribeiro CO, Brouwer CL, Oosting SF, Halmos GB, Witjes MJH, Langendijk JA, et al. Improving the prediction of overall survival for head and neck cancer patients using image biomarkers in combination with clinical parameters. Radiother Oncol J Eur Soc Therapeut Radiol Oncol. 2017;124(2):256–62.

    Article  Google Scholar 

  31. 31.

    Bogowicz M, Jochems A, Deist TM, Tanadini-Lang S, Huang SH, Chan B, Waldron JN, Bratman S, O’Sullivan B, Riesterer O, et al. Privacy-preserving distributed learning of radiomics to predict overall survival and HPV status in head and neck cancer. Sci Rep. 2020;10(1):4542.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  32. 32.

    Chen LL, Nolan ME, Silverstein MJ, Mihm MC Jr, Sober AJ, Tanabe KK, Smith BL, Younger J, Michaelson JS. The impact of primary tumor size, lymph node status, and other prognostic factors on the risk of cancer death. Cancer. 2009;115(21):5071–83.

    PubMed  Article  PubMed Central  Google Scholar 

  33. 33.

    van Griethuysen JJM, Fedorov A, Parmar C, Hosny A, Aucoin N, Narayan V, Beets-Tan RGH, Fillion-Robin JC, Pieper S, Aerts H. Computational radiomics system to decode the radiographic phenotype. Can Res. 2017;77(21):e104–7.

    Article  CAS  Google Scholar 

  34. 34.

    Lee G, Lee HY, Park H, Schiebler ML, van Beek EJR, Ohno Y, Seo JB, Leung A. Radiomics and its emerging role in lung cancer research, imaging biomarkers and clinical management: State of the art. Eur J Radiol. 2017;86:297–307.

    PubMed  Article  PubMed Central  Google Scholar 

  35. 35.

    Tixier F, Le Rest CC, Hatt M, Albarghach N, Pradier O, Metges JP, Corcos L, Visvikis D. Intratumor heterogeneity characterized by textural features on baseline 18F-FDG PET images predicts response to concomitant radiochemotherapy in esophageal cancer. J Nuclear Med. 2011;52(3):369–78.

    Article  Google Scholar 

  36. 36.

    Lovinfosse P, Polus M, Van Daele D, Martinive P, Daenen F, Hatt M, Visvikis D, Koopmansch B, Lambert F, Coimbra C, et al. FDG PET/CT radiomics for predicting the outcome of locally advanced rectal cancer. Eur J Nuclear Med Mol Imaging. 2018;45(3):365–75.

    Article  Google Scholar 

  37. 37.

    Oh JS, Kang BC, Roh JL, Kim JS, Cho KJ, Lee SW, Kim SB, Choi SH, Nam SY, Kim SY. Intratumor textural heterogeneity on pretreatment (18)F-FDG PET images predicts response and survival after chemoradiotherapy for hypopharyngeal cancer. Ann Surg Oncol. 2015;22(8):2746–54.

    PubMed  Article  PubMed Central  Google Scholar 

  38. 38.

    Yan J, Chu-Shern JL, Loi HY, Khor LK, Sinha AK, Quek ST, Tham IW, Townsend D. Impact of image reconstruction settings on texture features in 18F-FDG PET. J Nuclear Med. 2015;56(11):1667–73.

    CAS  Article  Google Scholar 

  39. 39.

    Mule S, Thiefin G, Costentin C, Durot C, Rahmouni A, Luciani A, Hoeffel C. Advanced hepatocellular carcinoma: pretreatment contrast-enhanced ct texture parameters as predictive biomarkers of survival in patients treated with sorafenib. Radiology. 2018;288(2):445–55.

    PubMed  Article  PubMed Central  Google Scholar 

  40. 40.

    Janmaat VT, Steyerberg EW, van der Gaast A, Mathijssen RH, Bruno MJ, Peppelenbosch MP, Kuipers EJ, Spaander MC. Palliative chemotherapy and targeted therapies for esophageal and gastroesophageal junction cancer. Cochrane Database Syst Rev. 2017;11:CD004063.

    PubMed  PubMed Central  Google Scholar 

  41. 41.

    Bolm L, Kasmann L, Paysen A, Karapetis C, Rades D, Wellner UF, Keck T, Watson DI, Hummel R, Hussey DJ. Multimodal anti-tumor approaches combined with immunotherapy to overcome tumor resistance in esophageal and gastric cancer. Anticancer Res. 2018;38(6):3231–42.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  42. 42.

    Hayes T, Smyth E, Riddell A, Allum W. Staging in esophageal and gastric cancers. Hematol Oncol Clin North Am. 2017;31(3):427–40.

    PubMed  Article  PubMed Central  Google Scholar 

  43. 43.

    Win T, Miles KA, Janes SM, Ganeshan B, Shastry M, Endozo R, Meagher M, Shortman RI, Wan S, Kayani I, et al. Tumor heterogeneity and permeability as measured on the CT component of PET/CT predict survival in patients with non-small cell lung cancer. Clin Cancer Res. 2013;19(13):3591–9.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  44. 44.

    Badic B, Desseroit MC, Hatt M, Visvikis D. Potential complementary value of noncontrast and contrast enhanced CT radiomics in colorectal cancers. Academic Radiol. 2019;26(4):469–79.

    Article  Google Scholar 

  45. 45.

    Moons KG, Altman DG, Reitsma JB, Ioannidis JP, Macaskill P, Steyerberg EW, Vickers AJ, Ransohoff DF, Collins GS. Transparent Reporting of a multivariable prediction model for Individual Prognosis or Diagnosis (TRIPOD): explanation and elaboration. Ann Intern Med. 2015;162(1):W1-73.

    PubMed  Article  PubMed Central  Google Scholar 

Download references


Not applicable.


This study was supported by the following grants: Science and Technology Planning Project of Guangdong Province (2016ZC0166 and 2017ZC0279); Shantou University Medical College Clinical Research Enhancement Initiative (Grant Number: N0201424); Collaborative and Creative Center, Molecular Diagnosis and Personalized Medicine, Shantou University, Guangdong, China (grant number: 2501220302); Special Fund for Guangdong Science and Technology Innovation Strategy (Grant Number: 180918114960704); Project category of Medical and Health, Shantou, Guangdong Province, China (Grant Number: 190816095262212); Science and Technology Special Fund of Guangdong Province, China (190829105556145); Strategic and Special Fund for Science and Technology Innovation of Guangdong Province of China (180918114960704).

Author information




Conception and design of the study: CBZ and CZC. Acquisition of data: CBZ and TTZ. Analysis and interpretation of the data: CBZ, TTZ and CZC. All authors participated in clinical data acquisition. Writing and revision of the manuscript: CBZ, TTZ, JZC and CZC. All authors have read and approved the final manuscript.

Corresponding author

Correspondence to Chuangzhen Chen.

Ethics declarations

Ethics approval and consent to participate

The Cancer Hospital of Shantou University Medical College institutional ethics committee approved this retrospective study (IRB No. 2019036), and the need for informed consent has been waived. This study complies with the standards of the Declaration of Helsinki and current ethics guidelines.

Consent for publication

Not applicable.

Competing interests

The authors have no conflicts of interest to declare.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1:

Supplementary Tables.

Additional file 2:

Figure S1.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Zeng, C., Zhai, T., Chen, J. et al. Imaging biomarkers of contrast-enhanced computed tomography predict survival in oesophageal cancer after definitive concurrent chemoradiotherapy. Radiat Oncol 16, 8 (2021).

Download citation


  • Oesophageal cancer
  • Imaging biomarker
  • Computed tomography
  • Prognosis