Patient characteristics
We identified 252 patients with pathologically confirmed inoperable ESCC treated with dCCRT from 2010 to 2019 in our study, with complete information of GTVp and GTVnd. The demographic characteristics of the ESCC patients are shown in Table 1. The age of the patients ranged from 30 to 82 years with a median age of 60 years. Median GTVp and GTVnd were 38.77 cm3 (interquartile range (IQR), 24.00-58.02 cm3) and 4.34 cm3 (IQR, 0.18-11.38 cm3), respectively. And the median ratio between the two (GTVnd/GTVp) was 0.11 with an IQR of 0.01-0.30.
Table 1. Characteristics of the patients with esophageal squamous cell carcinoma treated with definitive concurrent chemoradiotherapy (n=252)
Patient characteristics
|
No of patients (%)
|
Age
|
|
Median
|
60
|
Range
|
30-82
|
Gender
|
|
Male
|
220(87.3)
|
Female
|
32(12.7)
|
KPS
|
|
≥90
|
173(68.7)
|
<90
|
79(31.3)
|
Weight loss
|
|
Yes
|
89(35.3)
|
No
|
163(64.7)
|
Pain of chest and back
|
|
Yes
|
51(20.2)
|
No
|
201(79.8)
|
Tumor location
|
|
Cervical
|
22(8.7)
|
Upper thoracic
|
78(31.0)
|
Middle thoracic
|
93(36.9)
|
Lower thoracic
|
59(23.4)
|
Clinical T stage, 8th
|
|
T1
|
9(3.6)
|
T2
|
16(6.3)
|
T3
|
139(55.2)
|
T4
|
88(34.9)
|
Clinical N stage, 8th
|
|
N0
|
59(23.4)
|
N1
|
86(34.1)
|
N2
|
77(30.6)
|
N3
|
30(11.9)
|
Clinical TNM stage, 8th
|
|
I
|
7(2.8)
|
II
|
42(16.7)
|
III
|
102(40.5)
|
IV
|
101(40.1)
|
GTVp (cm3)
|
|
Median
|
38.77
|
IQR
|
24.00-58.02
|
GTVnd (cm3)
|
|
Median
|
4.34
|
IQR
|
0.18-11.38
|
GTVnd/GTVp
|
|
Median
|
0.11
|
IQR
|
0.01-0.30
|
Tumor length, cm
|
|
≤6
|
126(50.0)
|
>6
|
126(50.0)
|
Radiation dose, Gy
|
|
<54
|
75(29.8)
|
≥54
|
177(70.2)
|
Consolidation chemotherapy
|
|
Yes
|
106(42.1)
|
No
|
146(57.9)
|
Treatment outcomes
Within the median duration of follow-up for the whole patients of 38 months (ranging from 1 to 89 months), we identified 117 (46.4%) patients developing locoregional recurrence (LRR), 62 (24.6%) DM and 132(52.4%) death. The five-year OS, PFS, DMFS, LRFS and RRFS rates for the entire cohort were 30.6%, 28.4%, 58.0%, 46.9%, and 66.7%, respectively.
Prognostic value of volumetric parameters on survivals
GTVp, GTVnd and GTVnd/GTVp were segregated into two groups according to their median values. Subsequently, we investigated the prognostic role of the volumetric parameters using univariate and multivariate Cox proportional hazard regression. In univariate analysis (Table 2), GTVnd ≥4.34 cm3 and GTVnd/GTVp ≥0.11 were associated with poorer OS (both P<0.001), PFS (P=0.001; P=0.014; respectively), DMFS (both P<0.001) and RRFS (P=0.023; P=0.031; respectively). Then the variables with statistical significance in univariate analysis(P<0.05) were included in the multivariate analysis and P values were calculated using the Cox regression forward-LR model. Multivariable analyses validated the independent prognostic role of the GTVnd in OS (HR=1.949, P<0.001), PFS (HR=1.425, P=0.048) and DMFS (HR=2.548, P=0.001). In addition, the result demonstrated advanced clinical T stage (P=0.002) and male (P=0.046) were independently associated with shorter OS, shown as Table 3.
Table 2 Univariate analysis of prognostic factors in patients with ESCC (n=252)
Variable
|
OS
P
|
PFS
P
|
DMFS
P
|
RRFS
P
|
LRFS
P
|
Age (>60 years vs. ≤60 years)
|
0.796
|
0.923
|
0.412
|
0.116
|
0.640
|
Gender
|
0.034
|
0.005
|
0.025
|
0.016
|
0.084
|
KPS (≥90 vs. <90)
|
0.343
|
0.043
|
0.048
|
0.144
|
0.105
|
Weight loss
|
0.122
|
0.637
|
0.510
|
0.516
|
0.628
|
Pain of chest and back
|
0.973
|
0.665
|
0.167
|
0.602
|
0.783
|
Tumor location
|
0.010
|
<0.001
|
0.620
|
0.015
|
0.059
|
Clinical T stage
|
0.001
|
0.011
|
0.184
|
0.324
|
0.103
|
Clinical N stage
|
0.001
|
0.341
|
0.152
|
0.639
|
0.923
|
GTVp (≥38.77 vs.<38.77cm3)
|
0.167
|
0.496
|
0.945
|
0.475
|
0.480
|
GTVnd (≥4.34 vs.<4.34cm3)
|
< 0.001
|
0.001
|
< 0.001
|
0.023
|
0.464
|
GTVnd/GTVp (≥0.11 vs.<0.11)
|
< 0.001
|
0.014
|
< 0.001
|
0.031
|
0.887
|
Tumor length (>6cm vs. ≤6cm)
|
0.048
|
0.290
|
0.847
|
0.259
|
0.342
|
Radiation dose (≥54Gy vs. <54Gy)
|
0.242
|
0.127
|
0.060
|
0.093
|
0.110
|
Consolidation chemotherapy
|
0.633
|
0.189
|
0.423
|
0.110
|
0.478
|
Table 3 Multivariate cox proportional hazard regression analysis of prognostic factors in patients with ESCC (n=252)
|
HR (95%CI)
|
P
|
OS
|
|
|
T stage
|
|
0.002
|
T stage (T2 vs. T1)
|
2.492(0.290-21.405)
|
0.405
|
T stage (T3 vs. T1)
|
3.341(0.472-24.925)
|
0.223
|
T stage (T4 vs. T1)
|
6.134(0.844-44.550)
|
0.073
|
GTVnd (≥4.34 vs.<4.34cm3)
|
1.949(1.353-2.808)
|
< 0.001
|
Gender (male vs. female)
|
1.939(1.013-3.712)
|
0.046
|
PFS
|
|
|
Gender (male vs. female)
|
2.373(1.242-4.535)
|
0.009
|
GTVnd (≥4.34 vs.<4.34cm3)
|
1.425(1.003-2.024)
|
0.048
|
Tumor location
|
|
0.001
|
Tumor location (UT vs. Cervical)
|
0.274(0.146-0.514)
|
< 0.001
|
Tumor location (MT vs. Cervical)
|
0.543(0.310-0.950)
|
0.032
|
Tumor location (LT vs. Cervical)
|
0.475(0.261-0.8)
|
0.015
|
DMFS
|
|
|
GTVnd (≥4.34 vs.<4.34cm3)
|
2.548(1.491-4.355)
|
0.001
|
To confirm the optimal GTVnd cutoff value, we used the RCS with 3 knots and OS, PFS, DMFS, and RRFS as endpoint events. The cutoff value determined by RCS was 4.35 cm3 in this analysis (Fig.1). Multivariable analysis also showed that GTVnd≥4.35 cm3 was independent and significant negative prognostic factors for OS, PFS and DMFS.
Construction of Risk Grouping Using GTVnd by RPA model
Considering the prognostic value of GTVnd, we then performed RPA algorithm including T stage, GTVnd/GTVp, and gender to develop a new staging. The significant RPA-derived splits were only the T stage and GTVnd (Fig.2). The RPA model divided the 252 ESCC patients into the following three groups: high risk (T1-4 GTVnd≥4.35, n=126, III stage) , intermediate risk (T4 GTVnd<4.35,n=38,II stage), and low risk(T1-3GTVnd<4.35, n=88, I stage).
Clinical and treatment characteristics of the new risk grouping is shown in Table 4. The differences in age, gender, KPS, weight loss, pain of chest and back, tumor length, radiation dose, induction chemotherapy and consolidation chemotherapy were not statistically significant between the two groups (P > 0.05). Compared with those in the high-risk group, patients in the low-risk and intermediate-risk groups whose tumor were more located in the cervical and upper thoracic (48.9%, 52.6% vs.29.3%, P=0.033). And the proportion of N0-1 was higher in low-risk and intermediate-risk groups than that in the high-risk group (87.5%, 81.6% vs. 44.8%, P<0.001).
Table 4. Baseline characteristics of patients stratified by new risk stratification
Characteristic
|
Low-risk group
(n=88)
|
Intermediate-risk group
(n=38)
|
High-risk group
(n=126)
|
P
|
Age
|
|
|
|
0.751
|
≤60
|
44(50.0%)
|
19(50.0%)
|
69(54.8%)
|
|
>60
|
44(50.0%)
|
19(50.0%)
|
57(45.2%)
|
|
Gender
|
|
|
|
0.165
|
Male
|
76(86.4%)
|
30(78.9%)
|
114(90.5%)
|
|
Female
|
12(13.6%)
|
8(21.1%)
|
12(9.5%)
|
|
KPS
|
|
|
|
0.083
|
≥90
|
53(60.2%)
|
26(68.4%)
|
94(74.6%)
|
|
<90
|
35(39.8%)
|
12(31.6%)
|
32(25.4%)
|
|
Weight loss
|
|
|
|
0.127
|
Yes
|
27(30.7%)
|
10(26.3%)
|
52(41.3%)
|
|
No
|
61(69.3%)
|
28(73.7%)
|
74(58.7%)
|
|
Pain of chest and back
|
|
|
|
0.441
|
Yes
|
14(15.9%)
|
8(21.1%)
|
29(23.0%)
|
|
No
|
74(84.1%)
|
30(78.9%)
|
97(77.0%)
|
|
Tumor location
|
|
|
|
0.033
|
Cervical
|
7(8.0%)
|
4(10.5%)
|
11(8.7%)
|
|
Upper thoracic
|
36(40.9%)
|
16(42.1%)
|
26(20.6%)
|
|
Middle thoracic
|
28(31.8%)
|
12(31.6%)
|
53(42.1%)
|
|
Lower thoracic
|
17(19.3%)
|
6(15.8%)
|
36(28.6%)
|
|
Clinical T stage, 8th
|
|
|
|
<0.001
|
T1
|
7(8.0%)
|
0(0.0%)
|
2(1.6%)
|
|
T2
|
12(13.6%)
|
0(0.0%)
|
4(3.2%)
|
|
T3
|
69(78.4%)
|
0(0.0%)
|
70(55.6%)
|
|
T4
|
0(0.0%)
|
38(100.0%)
|
50(39.7%)
|
|
Clinical N stage, 8th
|
|
|
|
<0.001
|
N0
|
37(42.0%)
|
21(55.3%)
|
1(0.8%)
|
|
N1
|
40(45.5%)
|
10(26.3%)
|
36(28.6%)
|
|
N2
|
9(10.2%)
|
3(7.9%)
|
65(51.6%)
|
|
N3
|
2(2.3%)
|
4(10.5%)
|
24(19.0%)
|
|
GTVnd
|
|
|
|
<0.001
|
≥4.35
|
0(0.0%)
|
0(0.0%)
|
126(100.0%)
|
|
<4.35
|
88(100.0%)
|
38(100.0%)
|
0(0.0%)
|
|
Tumor length, cm
|
|
|
|
0.211
|
≤6
|
49(55.7%)
|
21(55.3%)
|
56(44.4%)
|
|
>6
|
39(44.3%)
|
17(44.7%)
|
70(55.6%)
|
|
Radiation dose, Gy
|
|
|
|
0.397
|
<54
|
26(29.5%)
|
8(21.1%)
|
41(32.5%)
|
|
≥54
|
62(70.5%)
|
30(78.9%)
|
85(67.5%)
|
|
Consolidation chemotherapy
|
|
|
|
0.174
|
Yes
|
33(37.5%)
|
21(55.3%)
|
52(41.3%)
|
|
No
|
55(62.5%)
|
17(44.7%)
|
74(58.7%)
|
|
The prognostic significance of the risk group
Then, we performed the Kaplan-Meier analysis to compare OS, PFS, DMFS, RRFS, and LRFS between the three groups derived by RPA (Fig.3). We found highly significant differences in OS among the three groups (P < 0.001; Fig. 3A), with corresponding 5-year OS rates of 17.0 % for high-risk group, 26.3% for intermediate-risk group, and 54.0% for low-risk group. The 5-year PFS rates of high-risk group, intermediate-risk group and low-risk group were 16.6%, 27.6% and 39.8%, respectively. By the log–rank test, there were significant differences in PFS among the three groups (P=0.002; Fig.3B). And the 5-year DMFS rates of 37.8% for patients with high-risk group and 72.0% for those with low-risk group showed significant differences (P =0.001; Fig.3C). Though a prognostic analysis demonstrated the 5-year RRFS rate in the high-risk group (59.1%) was lower than that of the intermediate-risk group (66.0%) and low-risk group (76.6%), the difference was not significant (P=0.063; Fig.3D). No significant difference in LRFS rate was observed among the three groups (P=0.194; Fig.3E).
Comparison of the performance of the new risk grouping and TNM staging system
The performance between the new risk grouping (TGTVndM) and TNM staging systems assessed by linear trend χ2 , likelihood ratio χ2, and the AIC is presented in Table 5. The TGTVndM staging system demonstrated higher linear trend χ2 (26.38 versus 25.77), high likelihood ratio χ2(24.39 versus 20.69), and lower AIC (1255.07 versus 1260.06) compared with the TNM stage, indicating the optimum prognostic stratification in predicting the survival of ESCC patients treated with dCCRT.
Table 5 Comparison of the prognostic performance of the TGTVndM and TNM staging systems
Classification
|
Subgroups
|
Figure
|
Linear trend χ2
|
Likelihood ratio χ2
|
AIC
|
N stage
|
N0, N1, N2, N3
|
SA
|
5.06
|
11.68
|
1271.29
|
GTVnd stage
|
GTVnd<4.35cm3
GTVnd≥4.35cm3
|
SB
|
16.23
|
18.48
|
1260.23
|
TNM stage
|
I, II, III, IV
|
SC
|
25.77
|
20.69
|
1260.06
|
TGTVndM stage
|
I, II, III
|
3A
|
26.38
|
24.39
|
1255.07
|