Demographic baseline characteristics
This study enrolled 1362 pediatric patients with osteosarcoma based on inclusion and exclusion criteria.The mean age of the patients was 12.8 years. Among them, 1015 (74.5 %) were white and 210 (15.4%) were black. 754 (56.1%) were male and 598(43.9%) were female. Grade included 23 (1.69%) I, 44 (3.23%) II, 320 (23.5%) III, 615 (45.2%) IV, and 360 (26.4%) unknow. T stage included 497 (36.5%) T1, 817 (60.0%) T2, 46 (3.38%) T3, 2 (0.15%) T4. 1332 patients (97.8%) had stage N0. The mean tumor size was 105 mm. The primary site of osteosarcoma in 1237(90.8%) children was in the limbs. 1,292(94.9%) received chemotherapy and only 43(3.16%) received radiation. There were 79 (5.80%) patients without surgery, 111 (8.15%) patients with partial resection, 885 (65.0%)patients with radical excision, and 287 (21.1%) patients with amputation. Cancer-specific survival (CSS) was 0.73. The mean survival time was 64.6 months. By 2018, 71 percent of patients were alive. The clinicopathological data of all patients are shown in Table 1. Based on the baseline information, the training cohort and validation cohort showed no statistically significant differences (P > 0.05).
Table 1. Clinicopathological information descrition table of pediatric osteosarcoma group.
|
ALL
|
Training set
|
Validation set
|
p
|
|
N=1362
|
N=965
|
N=397
|
|
Age
|
12.8 (3.47)
|
12.8 (3.44)
|
12.8 (3.54)
|
0.954
|
Race
|
|
|
|
0.999
|
white
|
1015 (74.5%)
|
719 (74.5%)
|
296 (74.6%)
|
|
black
|
210 (15.4%)
|
149 (15.4%)
|
61 (15.4%)
|
|
other
|
137 (10.1%)
|
97 (10.1%)
|
40 (10.1%)
|
|
Sex
|
|
|
|
0.075
|
Male
|
764 (56.1%)
|
526 (54.5%)
|
238 (59.9%)
|
|
Female
|
598 (43.9%)
|
439 (45.5%)
|
159 (40.1%)
|
|
Grade
|
|
|
|
0.601
|
I
|
23 (1.69%)
|
16 (1.66%)
|
7 (1.76%)
|
|
II
|
44 (3.23%)
|
30 (3.11%)
|
14 (3.53%)
|
|
III
|
320 (23.5%)
|
230 (23.8%)
|
90 (22.7%)
|
|
IV
|
615 (45.2%)
|
424 (43.9%)
|
191 (48.1%)
|
|
Unknown
|
360 (26.4%)
|
265 (27.5%)
|
95 (23.9%)
|
|
T
|
|
|
|
0.911
|
T1
|
497 (36.5%)
|
347 (36.0%)
|
150 (37.8%)
|
|
T2
|
817 (60.0%)
|
583 (60.4%)
|
234 (58.9%)
|
|
T3
|
46 (3.38%)
|
33 (3.42%)
|
13 (3.27%)
|
|
T4
|
2 (0.15%)
|
2 (0.21%)
|
0 (0.00%)
|
|
N
|
|
|
|
0.362
|
N0
|
1332 (97.8%)
|
941 (97.5%)
|
391 (98.5%)
|
|
N1
|
30 (2.20%)
|
24 (2.49%)
|
6 (1.51%)
|
|
Tumor size
|
105 (61.3)
|
107 (65.1)
|
101 (50.7)
|
0.080
|
Primary site
|
|
|
|
0.220
|
Limbs
|
1237 (90.8%)
|
870 (90.2%)
|
367 (92.4%)
|
|
Axial
|
125 (9.18%)
|
95 (9.84%)
|
30 (7.56%)
|
|
Laterality
|
|
|
|
0.390
|
Left
|
682 (50.1%)
|
486 (50.4%)
|
196 (49.4%)
|
|
Right
|
621 (45.6%)
|
433 (44.9%)
|
188 (47.4%)
|
|
Not pairs
|
59 (4.33%)
|
46 (4.77%)
|
13 (3.27%)
|
|
Chemotherapy
|
|
|
|
0.767
|
No/Unknown
|
70 (5.14%)
|
48 (4.97%)
|
22 (5.54%)
|
|
Yes
|
1292 (94.9%)
|
917 (95.0%)
|
375 (94.5%)
|
|
Radiation
|
|
|
|
0.040
|
No/Unknown
|
1319 (96.8%)
|
928 (96.2%)
|
391 (98.5%)
|
|
Yes
|
43 (3.16%)
|
37 (3.83%)
|
6 (1.51%)
|
|
Surgery
|
|
|
|
0.082
|
No
|
79 (5.80%)
|
66 (6.84%)
|
13 (3.27%)
|
|
Partial resection
|
111 (8.15%)
|
76 (7.88%)
|
35 (8.82%)
|
|
Radical excision
|
885 (65.0%)
|
622 (64.5%)
|
263 (66.2%)
|
|
Amputation
|
287 (21.1%)
|
201 (20.8%)
|
86 (21.7%)
|
|
Mets at lung
|
|
|
|
0.603
|
No
|
1105 (81.1%)
|
779 (80.7%)
|
326 (82.1%)
|
|
Yes
|
257 (18.9%)
|
186 (19.3%)
|
71 (17.9%)
|
|
CSS
|
0.73 (0.44)
|
0.72 (0.45)
|
0.77 (0.42)
|
0.049
|
Survival months
|
64.6 (49.8)
|
63.3 (48.8)
|
68.0 (52.0)
|
0.119
|
Univariate and Multivariate Logistic Regression Analysis
Using univariate logistic regression model, seven important prognostic factors were screened, including age (hazard ratio (HR)=1.06, 95% CI 1.01-1.11, P < 0.05), laterality (HR=0.27, 95% CI 0.08-0.87, P < 0.05), T stage (HR=4.51, 95% CI 2.06-9.88, P < 0.05), N stage (HR=6.26, 95% CI 2.73-14.33, P < 0.05), surgery (HR=0.39, 95% CI 0.22-0.69, P < 0.05), radiotherapy (HR=6.07, 95% CI 3.1-11.88, P < 0.05), tumor Size (HR=1.001, 95% CI 1-1.01, P < 0.05). Then,multivariate logistic regression was used to analyze these factors.There were four independent prognostic factors, which were T stage(HR=2.194, 95% CI 1.457-3.305, P < 0.05), N stage (HR=5.288, 95% CI 2.155-12.975, P < 0.05), surgery (HR=0.366, 95% CI 0.195-0.688, P < 0.05) and radiotherapy (HR=5.155, 95% CI 2.417-10.996, P < 0.05).(Table 2) We can combine these risk factors into a nomogram to predict the risk of pulmonary metastasis of pediatric osteosarcoma.
Table 2. Univariate and multivariate analyses of CSS in training cohort.
|
Univariate
|
Multivariate
|
|
HR
|
95%CI
|
P
|
HR
|
95%CI
|
P
|
Age
|
1.06
|
1.01-1.11
|
0.03
|
1.048
|
0.995-1.103
|
0.076
|
Race
|
|
|
|
|
|
|
white
|
reference
|
|
|
|
|
|
black
|
1.09
|
0.7-1.7
|
0.7
|
|
|
|
other
|
1.2
|
0.71-2.01
|
0.5
|
|
|
|
Sex
|
|
|
|
|
|
|
Male
|
reference
|
|
|
|
|
|
Female
|
0.79
|
0.57-1.1
|
0.16
|
|
|
|
Grade
|
|
|
|
|
|
|
I
|
reference
|
|
|
|
|
|
II
|
1.67
|
0.16-17.47
|
0.67
|
|
|
|
III
|
3.25
|
0.42-25.33
|
0.26
|
|
|
|
IV
|
3.71
|
0.48-28.45
|
0.21
|
|
|
|
Unknown
|
4.11
|
0.53-31.78
|
0.18
|
|
|
|
T
|
|
|
|
|
|
|
T1
|
reference
|
|
|
|
|
|
T2
|
2.36
|
1.6-3.46
|
0
|
2.194
|
1.457-3.305
|
0
|
T3
|
4.51
|
2.06-9.88
|
0
|
3.888
|
1.666-9.072
|
0.002
|
T4
|
7.9
|
0.48-128.81
|
0.15
|
15.933
|
0.764-332.179
|
0.074
|
N
|
|
|
|
|
|
|
N0
|
reference
|
|
|
|
|
|
N1
|
6.26
|
2.73-14.33
|
0
|
5.288
|
2.155-12.975
|
0
|
Tumor size
|
1.001
|
1-1.01
|
0
|
|
|
|
Primary site
|
|
|
|
|
|
|
Limbs
|
reference
|
|
|
|
|
|
Axial
|
0.7
|
0.39-1.32
|
0.24
|
|
|
|
Laterality
|
|
|
|
|
|
|
Left
|
reference
|
|
|
|
|
|
Right
|
0.89
|
0.64-1.23
|
0.48
|
0.873
|
0.617-1.234
|
0.442
|
Not pairs
|
0.27
|
0.08-0.87
|
0.03
|
0.285
|
0.079-1.023
|
0.054
|
Chemotherapy
|
|
|
|
|
|
|
No/Unknown
|
reference
|
|
|
|
|
|
Yes
|
2.27
|
0.97-7.68
|
0.06
|
|
|
|
Radiation
|
|
|
|
|
|
|
No/Unknown
|
reference
|
|
|
|
|
|
Yes
|
6.07
|
3.1-11.88
|
0
|
5.155
|
2.417-10.996
|
0
|
Surgery
|
|
|
|
|
|
|
No
|
reference
|
|
|
|
|
|
Partial resection
|
0.18
|
2.73-11.33
|
0
|
0.211
|
0.088-0.505
|
0
|
Radical excision
|
0.22
|
0.08-0.41
|
0
|
0.243
|
0.137-0.432
|
0
|
Amputation
|
0.39
|
0.13-0.38
|
0
|
0.366
|
0.195-0.688
|
0.002
|
Construction and validation a novel nomogram
We constructed a novel Nomogram to predict the risk of pulmonary metastasis for pediatric osteosarcoma (Figure 2). T stage was the most significant risk factor for pulmonary metastasis, followed by N stage, radiotherapy and surgery. Verify the nomogram accuracy of the training and validation queues using calibration curves. Calibration curves showed a good degree of consistency between the two groups..(Figure 3). In the training cohorts and validation cohorts, the C-index was 0.699 (95% CI, 0.656-0.741) and 0.736(95% CI, 0.675-0.797), respectively, indicating good accuracy. The AUC of the training cohort and the validation cohort also suggested excellent predictive power, with 0.676 (95%CI, 0.632–0.719) and 0.587 (95%CI, 0.517–0.657), respectively(Figure 4).
Clinical Application of Nomogram
DCA shown good clinical application in training and validation cohorts(Figure 5). Nomogram has significant predictive advantages over traditional TN staging. There was a high-risk group (total score > 71.3) and a low-risk group (≤71.3 overall) in the training and validation cohorts. The K-M curve indicated that the 1-, 3-, and 5-year survival rates in the high-risk group were 0.936 (95% CI, 0.917–0.955),0.727 (95% CI, 0.692–0.764), and 0.655 (95% CI, 0.617–0.696), respectively. The 1-, 3- and 5-year survival rates of patients in the low-risk group were 0.964 (95% CI, 0.949–0.978), 0.847 (95%CI, 0.818–0.877), and 0.783 (95% CI, 0.749–0.819), respectively(Figure 6). Outcomes proved the high-risk group had a poor prognosis and was more prone to distant metastasis.We analyzed surgical procedures in the high-risk and low-risk groups.In the low-risk group, partial resection had the highest survival rate, followed by radical resection and amputation. Patients who did not have surgery had the lowest survival rates(Figure 7A). In the high-risk group, most children received radical resection, and patients undergoing radical resection possessed substantially higher survival rates than other groups. Patients who underwent surgery had significantly higher survival rates than those who did not (Figure 7B).