SEER-Based Survival Nomogram (1998-2015) Based on ‘Stage,Lymph Node Dissection,Tumor Size and Degree of Differentiation,and Therapies’ for Prognosis of Primary Pulmonary Sarcoma

Abstract

Objective: Primary pulmonary sarcoma (PPS) is very rare in terms of incidence, henceforth, the clinical evidence pertinent to the prognosis of PPS is limited. The aim of this study was to construct a nomogram for evaluating the overall survival (OS) of patients diagnosed with PPS based on the stage, lymph node dissection, tumor size and degree of differentiation, and therapies. Methods: A total of 515 patients diagnosed with PPS during the period of 1998 to 2015 were obtained from the surveillance, epidemiology, and end results database and randomly segregated into ‘training group’ and ‘validation group’ with a ratio of 7:3. Regression analysis was executed for the training group to obtain the independent factors influencing prognosis of PPS patients. A nomogram was constructed as per the results obtained through multivariate Cox regression analysis subsequently validated using C index, receiver operating characteristic (ROC) curve, and calibration curves. Results: Age, tumor size, histology type, lymph node surgery, summary stage and differentiation grade were independent factors affecting the prognosis. C index was 0.775 and 0.737 for both training group, and validation group, respectively. Areas under the ROC curve of 1-year, 3-year, and 5-year OS were 87.6 (95% CI: 83.8-91.3), 90.1 (95% CI: 86.2-94.0) and 90.6 (95% CI: 85.8-95.4), respectively, in training group. Area under the curve values of 1-year, 3-year, and 5-year OS in the validation group were 83.1 (95% CI: 75.8-90.5), 82.9 (95% CI: 73.2-92.7) and 87.0 (95% CI: 75.9-98.1), respectively. Based on the nomogram, patients were segregated into low-risk group and high-risk group (degree of risk: cutoff score 193). OS of low-risk group was significantly higher when compared to high-risk group (P < .001) in the training group and validation group. Radiotherapy was a risk factor for the low-risk group and adjuvant chemotherapy has not exhibited influence on OS pertinent to low-risk group. However, adjuvant radiotherapy or chemotherapy both significantly improved the prognosis of PPS patients (P < .001) in the high-risk group. Conclusion: Constructed nomogram could have a strong predictive ability with higher accuracy for the prognosis of patients with PPS. Patients at low risk could not benefit from adjuvant radiotherapy or chemotherapy, while the prognosis clearly improved in the high-risk populations treated with either radiotherapy or chemotherapy.

Keywords

primary pulmonary sarcoma (PPS)therapies nomograms prognosis histology type lymph node surgery summary stage and differentiation grade

Introduction

Primary pulmonary sarcoma (PPS) is a very rare and aggressive non-epithelial (mesenchymal) malignancy originating in the mesenchymal tissue of bronchial wall, blood vessels, and lung mesenchyme.¹ Incidence of PPS accounts 0.4% to 1.1% among all other types of lung malignant tumors.² PPS infiltrates into the lung parenchyma and spreads into the bronchi.³ However, PPS is more aggressive and exhibits a worse prognosis followed by 35% of 5-year overall survival (OS) when compared to non-small cell lung cancer (NSCLC).^4–6 Lungs are the most metastatic sites for soft tissue sarcomas, and the formation of extrapulmonary metastatic sarcomas tend to be more common in terms of occurrence than PPS.^7,8

Diagnosis of PPS is considered as an intricate process due to the difficulty in distinguishing PPS with primary and extrapulmonary metastatic sarcoma, pulmonary epithelial tumors, pulmonary sarcomatoid carcinoma, and malignant melanoma.⁹ Therefore, a patient's comprehensive medical history, and a more detailed examination of biopsy, or pathological samples are required to confirm the sarcoma of lung origin.¹⁰

Very limited patients cases with PPS have been reported in the literature so far due to the low, and sporadic incidence rate.^11,12 Generally, the prediction of survival rate of patients with PPS is based on the personal experience of doctors.² Till now, there is no standard predictive model for the risk and prognosis of PPS patients. Therefore, it is necessary to establish a prediction model to help clinicians for evaluating the prognosis of PPS patients more accurately.

Five-year OS rate pertinent to PPS is influenced by several factors including patient age, tumor size, surgery, metastasis, degree of histological malignancy, and therapeutic modalities.^2,11,12 Tumor prognosis can be accurately predicted using nomograms. Data pertinent to PPS patients procured from the registry of ‘Surveillance, Epidemiology, and End Results (SEER) program database’. Thus, nomogram is an effective predictive model for evaluating tumor outcomes by integrating multiple prognostic factors based on variable regression coefficients.¹³ In certain cancers, nomogram has been proved to have better predictive ability than traditional tumor node and metastasis (TNM) classification.¹⁴ In present work, a nomogram was constructed based on the SEER database to evaluate the survival probability and prognostic factors more accurately for PPS patients subsequently to prefer personalized therapeutic regimen. In addition, PPS patients were stratified according to the degree of risk based on nomogram, and the influence of adjuvant radiotherapy and chemotherapy on the prognosis of different groups was analyzed.

Methods

Data Sources

Data collection was performed from the SEER database; the data acquired from this database was approved by the local ethical committee for analysis. Hence, our study has not required any ethical approval statement. The SEER program registry has the authoritative registered malignant tumor data in the United States which has been collecting data from 18 cancer registries, and covering 30% of the U.S. population. We obtained the information of patients diagnosed with PPS during the period of 1998 to 2015 in SEER database and collected the following clinical data including, age, sex, race, marital status, ICD-O-3 histological type, location of primary tumors, number of primary tumors, degree of differentiation, tumor size, involved lymph nodes, laterality, number of positive lymph nodes, surgery, lymph node resection, summary stage, radiotherapy, chemotherapy, follow-up time, survival status, and cause of death.

Inclusion and Exclusion Criteria

Inclusion criteria was performed according to the factors such as (1) the diagnosis time, 1998 to 2015, (2) the primary sites of PPS in lung and bronchus (C34.0-C34.9), (3) ICD-O-3 histological type confirmed by pathology: epithelioid hemangioendothelioma, hemangiosarcoma, pleuropulmonary blastoma, chondroma, monophasic/biphasic synovial sarcoma, pulmonary artery/vein sarcoma, and (4) exact follow-up information.

Exclusion criteria were executed according to the factors including (1) non-pathologically positive data and (2) unknown survival time.

Statistical Methods

The patient's age diagnosed with PPS and tumor diameter were analyzed by x-title software to obtain the best cut-off value in our study. Later, they were grouped as per the obtained value. Baseline characteristics of PPS patients were described with the aid of descriptive statistics. Categorical variables were compared by Chi-squared test or Fisher's exact test using SPSS21.0 software. The results of univariate and multivariate Cox risk regression analysis were expressed in risk ratio (HR) and corresponding 95% confidence interval.

Based on the cox proportional risk model, a nomogram was constructed using R4.1.0. Nomogram performance was evaluated by the concordance index (C-index), receiver operating characteristic (ROC) curves and calibration curves. Kaplan–Meier method was used to calculate the survival probability and draw the survival curve. Log-rank test was used to compare the differences between groups. A two-sided P value < .05 was considered statistically significant.

Establishment of Nomogram

Data were randomly segregated into ‘training group’ and ‘validation group’ (7:3 ratio) using R studio software “sample” package. For the training group, we used univariate and multivariate Cox regression analysis to predict the independent risk factors pertinent to OS; subsequently HR and corresponding 95% confidence interval were calculated. Forward logistic regression analysis was performed to eliminate redundant variables. Forest map was drawn with the aid of “forest” package in R statistical computing software to further visualize results of multivariate Cox regression analysis. This study used “rms”, “foreign” and “survival” packages to obtain the nomogram model and to assess the prediction ability and accuracy of nomogram. P values < .05 indicate statistical significance.

Validation of Nomograph

The accuracy and reliability of nomogram was verified in the training group and validation group respectively. C-Index was determined in order to assess the ability of a model for predicting occurrence of the event. The value range of C index was 0.5 to 1. When the C index was equal to 0.5, then the model cannot be distinguished. A higher value of C index indicates the stronger ability of its model to predict the occurrence of events. In addition, we used the “risk regression” package to obtain the ROC curve of multivariate Cox regression to quantify the prediction ability of nomogram. A higher area under the curve (AUC) of ROC depicts a stronger prediction and judgment ability of the model. Bootstrap self sampling method was used internally to verify the effectiveness of nomogram-survival prediction model. The calibration curve was obtained by resampling the above data for 1000 times. It can be used to evaluate the accuracy of the model based on the difference between the probability of events predicted by the model and the observed outcome frequencies. A closer predicted calibration curve to the standard curve indicates higher accuracy of the model.

Risk Stratification Based on Nomogram

Total score related to the individual patient was determined using nomogram formula package. The optimal cutoff value was calculated using the “surv_cutpoint” function. We divided the training group and validation group into high-risk group and low-risk group according to the best cut-off value for determining nomogram accuracy. Subsequently, we performed Kaplan–Meier survival analysis for assessing the significance of survival differences between low-risk and high-risk groups; P < .05 was considered as statistically significant.

Effects of Adjuvant Radiotherapy and Chemotherapy in Different Risk Groups

The influence of adjuvant radiotherapy and chemotherapy in high risk and low risk groups were deciphered as per the OS analysis and prognosis.

Results

The Best Cutoff Values for Age and Tumor Size

The results of X-title software analysis showed that the best cut-off values for age were 50 and 80 years old, therefore, the patients were segregated into three age groups viz., ≤50 years old, 50 to 80 years, and ≥80 years. A significant difference in OS among the three age groups (P < .05) was observed. The optimum cut-off values of tumor diameter were 82 mm and 130 mm and they were segregated into three groups such as ≤ 82 mm, 82 to 130 mm, and ≥ 130 mm; significant differences were observed among three patient groups in terms of OS (P < .05).

Clinicopathological Features of Training Group and Validation Group

Patients diagnosed for PPS during the time period of 1998 to 2015 were obtained from SEER database; a total of 515 eligible patients were included in this study. Training group composed of 360 patients whereas validation group composed of 155 patients after random division as per the ratio of 7:3. Significant differences in ‘tumor laterality distribution’ and ‘surgery on other sites’ between the two groups (P < .05) was observed. The demographic distribution and pathological characteristics of two groups were described in Table 1.

Table 1.

Demographic Characteristics Distribution and Clinical Characteristics of Training Group and Validation Group.

Variables	Training Group (n = 360)	Validation Group (n = 155)	Total (n = 515)	P-Value
Age (years)				0.771
≤50	90	37	127
50-80	202	92	294
≥80	68	26	94
Sex				0.495
Male	216	88	304
Female	144	67	211
Race				0.695
Black	35	17	52
White	301	125	426
Other	24	13	37
Marital status				0.662
Single	66	31	97
Married	200	77	277
Divorced/separated/widowed	81	41	122
Others	13	6	19
Primary site				0.890
Main bronchus	12	3	15
Upper lobe	136	53	189
Middle lobe	23	10	33
Lower lobe	100	47	147
Overlapping lesion	11	6	17
Unknown	78	36	114
Spread status				0.043
Unilateral	338	152	490
Bilateral	22	3	25
Histological type				0.853
Pulmonary artery/vein sarcoma	219	97	316
Pleuropulmonary blastoma	28	13	41
Monophasic/biphasic synovial sarcoma	47	18	65
Hemangiosarcoma	62	25	87
Epithelioid hemangioendothelioma	2	0	2
Chondroma	2	2	4
Tumor size (mm)				0.284
≤82	167	69	236
82-130	60	24	84
≥130	20	16	36
Unknown	113	46	159
Involved lymph nodes				0.169
None	203	79	282
Yes	127	55	182
Other	30	21	51
Number of primary tumors				0.895
1	246	105	351
>1	114	50	164
Differentiation				0.683
Well/moderate	20	12	32
Poor	77	33	110
Undifferentiated	87	41	128
Unknown	176	69	245
Summary stage				0.269
Localized	78	35	113
Regional	90	27	117
Distant	165	82	247
Unknown	27	11	38
Surgery				0.335
None	192	82	274
Sub-lobectomy	63	19	82
Lobectomy	72	32	104
Pneumonectomy	28	18	46
Unknown	5	4	9
Lymph node resection				0.365
None	265	110	375
Yes	90	40	130
Unknown	5	5	10
Surgery on other sites				0.022
None	342	138	480
Yes	17	14	31
Unknown	1	3	4
Radiation				0.114
None/unknown	279	110	389
Yes	81	45	126
Chemotherapy				0.098
None/unknown	241	92	333
Yes	119	63	182
Positive lymph nodes				0.964
None	82	37	119
Yes	21	9	30
Unknown	257	109	366

Univariate and Multivariate cox Regression Analysis

Univariate analysis of cox proportional regression model showed that age, sex, marital status, ICD-O-3 histological type, degree of differentiation, tumor size, involved lymph nodes, laterality, summary stage, number of positive lymph nodes, surgery of primary site, lymph node resection, and chemotherapy could significantly affect the survival state (P < .05). Consequently, we further incorporated these factors analyzed through multivariate cox proportional regression analysis (P < .05). Age, tumor size, histological type, lymph node resection, summary stage and degree of differentiation were significantly correlated with OS of the patients diagnosed with PPS. Forest map described the results of multivariate cox regression (Figure 1) (Table 2).

Figure 1.

Forest map depicts the results obtained through multivariate Cox proportional regression pertinent to the training group analysis.

Table 2.

The Outcomes Observed Through the Univariate and Multivariate Cox Proportional Regression Analysis in Training Group.

Variables	Univariate Analysis		Multivariate Analysis
Variables	HR (95%CI)	P-Value	HR (95%CI)	P-Value
Age(years)
≤50	Reference		Reference
50-80	2.563(1.860-3.533)	0.000	1.589(1.110-2.274)	0.011
≥80	4.466(3.063-6.512)	0.000	3.181(2.083-4.858)	0.000
Sex
Male	Reference
Female	0.761(0.600-0.964)	0.024
Race
Black	Reference
White	0.986(0.662-1.469)	0.946
Other	1.332(0.752-2.361)	0.325
Marital status
Single	Reference
Married	2.521(1.724-3.685)	0.000
Divorced/separated/widowed	3.723(2.451-5.655)	0.000
Others	1.561(0.744-3.272)	0.239
Primary site
Main bronchus	Reference
Upper lobe	1.149(0.582-2.27)	0.689
Middle lobe	1.016(0.452-2.281)	0.970
Lower lobe	0.932(0.466-1.861)	0.841
Over lapping lesion	0.855(0.328-2.226)	0.748
unknown	1.943(0.969-3.895)	0.061
Spread status
Unilateral	Reference
Bilateral	2.216(1.428-3.439)	0.000
Histological type
Pulmonary artery/vein sarcoma	Reference		Reference
Pleuropulmonary blastoma	0.088(0.036-0.215)	0.000	0.141(0.053-0.374)	0.000
Monophasic/biphasic synovial sarcoma	0.356(0.24-0.527)	0.000	0.489(0.319-0.752)	0.001
Hemangiosarcoma	0.95(0.703-1.283)	0.737	1.032(0.739-1.442)	0.853
Epithelioid hemangioendothelioma	0.522(0.125-2.178)	0.372	1.329(0.264-6.689)	0.730
Chondroma	0.253(0.035-1.804)	0.170	0.555(0.076-4.070)	0.562
Tumor size (mm)
≤82	Reference		Reference
82-130	1.162(0.830-1.627)	0.381	1.195(0.840-1.698)	0.322
≥130	2.432(1.513-3.910)	0.000	2.586(1.584-4.221)	0.000
Unknown	1.449(1.112-1.889)	0.006	1.315(0.964-1.795)	0.084
Involved lymph nodes
None	Reference
Yes	1.928(1.509-2.464)	0.000
Other	1.731(1.147-2.613)	0.009
Number of primary tumors
1	Reference
＞1	0.975(0.763-1.247)	0.841
Differentiation
Well/Moderate	Reference		Reference
Poor	2.993(1.576-5.684)	0.001	3.217(1.553-6.666)	0.002
Undifferentiated	3.275(1.727-6.21)	0.000	3.493(1.697-7.190)	0.001
Unknown	2.197(1.182-4.086)	0.013	2.841(1.407-5.738)	0.004
Summary stage
Localized	Reference
Regional	2.570(1.745-3.786)	0.000	2.383(1.592-3.568)	0.000
Distant	4.291(3.000-6.137)	0.000	2.673(1.800-3.969)	0.000
Unknown	3.65(2.206-6.037)	0.000	1.990(1.091-3.627)	0.025
Surgery
None	Reference
Sub-lobectomy	0.274(0.190-0.395)	0.000
Lobectomy	0.280(0.201-0.390)	0.000
Pneumonectomy	0.325(0.205-0.515)	0.000
Unknown	0.410(0.151-1.115)	0.000
Lymph node resection
None	Reference		Reference
Yes	0.478(0.360-0.634)	0.000	0.592(0.429-0.816)	0.001
Unknown	0.863(0.321-2.321)	0.770	2.397(0.846-6.794)	0.100
Surgery on other sites
None	Reference
Yes	1.371(0.814-2.308)	0.235
Unknown	0.780(0.109-5.565)	0.804
Radiation
None/unknown	Reference
Yes	1.066(0.813-1.398)	0.645
Chemotherapy
None/unknown	Reference
Yes	0.715(0.556-0.919)	0.009
Positive lymph nodes
None	Reference
Yes	1.493(0.867-2.569)	0.148
Unknown	2.055(1.529-2.763)	0.000

Construction of Nomogram

Six important prognostic factors including age, tumor size, histological type, lymph node resection, summary stage and degree of differentiation were determined for the training group using multivariate analysis. A nomogram was built to predict the 1-year, 3-year and 5-year OS of PPS patients based on these factors (Figure 2).

Figure 2.

Nomogram pertinent to 1-year, 3-year, and 5-year overall survival (OS). PB: pleuropulmonary blastoma; EH: epithelioid hemangioendothelioma; HS: hemangiosarcoma; CS: chondroma; SS: monophasic/biphasic synovial sarcoma; PS: pulmonary artery/vein sarcoma. I: Well/Moderate, II: Poor, III: Undifferentiated.

Verification of Nomograph

The predictive power of nomogram was verified by the C-index and ROC curves. The C index of patient prognosis model in the training group and the validation group were 0.775 and 0.737 respectively. In the training group, the areas under the ROC curve of 1-year, 3-year and 5-year OS were 87.6 (95% CI：83.8-91.3), 90.1 (95% CI：86.2-94.0), 90.6 (95% CI：85.8-95.4). AUC values of 1-year, 3-year and 5-year OS in the validation group were 83.1 (95% CI: 75.8-90.5), 82.9 (95% CI: 73.2-92.7) and 87.0 (95% CI: 75.9-98.1) respectively (Figure 3a, 3b).

Figure 3.

Receiver operating characteristic curve (ROC) and AUC values related to 1-year, 3-year, and 5-year OS in (A) training group and (B) validation group.

C index observed for both high risk and low risk groups was consistent with the AUC values, and concluded that the model exhibited higher prediction ability. The calibration curves were constructed for both training group and validation group in order to ascertain the predictive ability of nomogram (Figure 4a, 4b). They were the same as the actual prediction curves, indicating that the 1-year, 3-year, and 5-year OS predicted by nomogram was consistent with actual results, which indicated the accuracy of the model was relatively high.

Figure 4.

Calibration curves depict 1-year, 3-year, and 5-year overall survival in (A) training group and (B) validation group.

Risk Stratification Based on Nomogram

A total nomogram score was determined for individual patients. Subsequently, the best cut-off score was found to be 193 points. Patients were segregated into high-risk group (≥193 points) and low-risk group (<193 points) based on the best cut-off value of prognostic risk score. Kaplan–Meier survival curve showed that the OS rate of the high-risk group in the training group was lesser when compared to low-risk group (P < .05). This was confirmed in the validation group in which the same cut-off value was used to distinguish high-risk group from low-risk group. These curves demonstrated that there was a strong positive correlation between low-risk scores and OS (P < .05), which further verified the discrimination of the model ( Figure 5a , 5b).

Figure 5.

Kaplan–Meier curves depict overall survival (OS) in the patients diagnosed with PPS and differentiated into high-risk and low risk groups, respectively: (A) the training group, (B) validation group.

Role of Adjuvant Radiotherapy and Chemotherapy

Good OS was observed in the high risk group who were treated with adjuvant radiotherapy and adjuvant chemotherapy. Adjuvant radiotherapy was accompanied by the poor OS but the adjuvant chemotherapy was independent of OS in the low risk group ( Figure 6 ).

Figure 6.

Kaplan–Meier curves related to overall survival (OS) in different risk groups received different treatment regimens: (A) low risk patients with or without radiotherapy, (B) low risk patients with or without chemotherapy, (C) high risk patients with or without radiotherapy, and (D) high risk patients with or without chemotherapy.

Discussion

PPS is derived from lung mesenchymal tissue and is one of the significant pulmonary malignant tumors.¹⁵ Currently, a minimal number of clinical studies are available related to PPS which mainly focus on the case reports and retrospective analyses. Heterogeneity of PPS is one of the major driving factors making oncologists or clinicians difficult to stratify the prognosis of patients with PPS. The current study was a large cohort study mainly based on the data related to 515 PPS patients in the SEER database. It is aimed to find out the independent prognostic factors affecting the OS of PPS patients and construct a nomogram that can individually evaluate the prognosis of PPS patients.

Our study reported that PPS was more common in the male patients aged 50 to 80 years, and its occurrence mostly could be observed in white individuals. Most of the tumors originating in the upper and lower lobes of the lung were unilateral with a diameter of ≤82 mm. For patients with known tumor differentiation, most tumors are poorly differentiated or undifferentiated. More than half of the patients did not receive surgery, lymph node dissection, adjuvant chemotherapy or radiotherapy. Analysis of the training group concluded that age, tumor size, summary stage and degree of differentiation were independent risk factors affecting the prognosis of PPS patients. However, lymph node dissection and histological types (synovial sarcoma, pleuropulmonary blastoma and chondroma) could improve the OS of PPS patients.

Six variables were selected and incorporated into the nomogram. C index, ROC curve and calibration curve were used to verify the performance and accuracy of nomogram. The results showed that nomogram exhibited good prediction ability and could be used as a reference standard for clinicians to assess the prognosis of patients with different PPS. In addition, we divided the patients into high-risk group and low-risk group, and the OS of the low-risk group was significantly higher than that of the high-risk group, which further verified the discriminative ability of the model. Adjuvant chemoradiotherapy was not an independent prognostic factor for PPS in this study, and there was no unified view on the effect of adjuvant chemoradiotherapy on PPS. Therefore, we performed a subgroup analysis of PPS patients. The results showed that radiation therapy was a risk factor for patients in the low-risk population, while adjuvant chemotherapy was not associated with OS. In high-risk groups, both radiotherapy and adjuvant chemotherapy could improve the prognosis of PPS patients.

Furthermore, most of the PPS patients were diagnosed during the age of 50 to 80 years, and age was an independent risk factor affecting the prognosis of PPS patients. White people were the highest in number in PPS patients, but ethnic differences were not related to the prognosis of PPS patients. These results were consistent with the previous literature.¹⁶ A total of 60% PPS patients were males rather than females, but the difference in the incidence rate between the gender groups was not as obvious as that of lung cancer.¹⁷ Gender was not associated with the prognosis of PPS in present study. However, men exhibited worse OS and PFS than women among the patients who underwent surgery.⁷

Malignant fibrous histiocytoma, synovial sarcoma and smooth muscle sarcoma are the three most common histological types of PPS.^18,19 However, the most common histological type of PPS among 515 patients was pulmonary artery/vein sarcoma, followed by angiosarcoma. This may be related to the polymorphism and heterogeneity of PPS in histology or the continuous changes of PPS classification standards in recent years. Although some studies have not typically shown the effect of the histological type of PPS on OS in the patients with PPS, our findings suggested that the patients with pathological types mainly related to synovial sarcoma and pleuropulmonary blastoma exhibited a better prognosis.^20,21

Many studies have shown that the tumor diameter could play a key role in evaluating the prognosis of PPS patients. Tumor diameter of 5 cm was an important factor in assessing the development stage of soft tissue sarcoma.²² Janssen et al¹⁷ divided PPS patients into two groups (diameter ≤ 4 cm and > 4 cm) according to the median tumor diameter (4 cm) of PPS patients. In this study,¹⁷ authors reported the 5-year survival rate for malignant fibrous histiocytoma, synovial sarcoma and leiomyosarcoma was 30%, 50%, and 60% respectively. According to this study, the survival rate of patients with tumor diameter ≤ 4 cm was higher. Nascimento et al²³ summarized the clinical data of 18 patients diagnosed with PPS and the patients with tumor diameter less than 3 cm exhibited a better prognosis. The results of our study showed that patients with tumor diameter greater than 8.2 cm reported a poor prognosis. Therefore, the tumor diameter could be used as one of the significant factors for determining the prognosis of PPS patients, and a large tumor diameter can be a risk factor for the prognosis of PPS patients. However, our study has limitations such as small sample size, single evaluation index and different selection of dividing points. Majority of the patients with PPS exhibited a high degree of histological malignancy; the high-grade tumors were considered as adverse factors affecting the prognosis of patients with PPS. This is consistent with previous findings that a higher tumor grade is associated with lower OS in the patients with PSS.^16,17,23

Staging enables a better evaluation of the prognosis pertinent to several cancers in order to prefer a suitable personalized treatment. However, there are very limited reports available specifically related to the TNM staging standard suitable for PPS. The staging systems for lung cancer, trunk and limb soft tissue sarcoma were often used as the staging standards of PPS in previous studies.^20,21,24,25 In the eighth edition of AJCC (American Joint Committee on cancer) cancer staging manual, a special staging system for STS of abdominal and thoracic internal organs was proposed for the first time, indicating that PPS could be adapted to this system. But Collaud et al² reported that the new staging system did not provide any value in predicting death or recurrence in their study cohort. Previous studies showed that the lung cancer staging system could accurately predict the survival and prognosis of patients with PPS and was better suitable for PPS.^19–21 However, the SEER database has not reported AJCC staging for PPS, and there was no information available on pathological lymph nodes and metastases. Therefore, we were unable to perform TNM staging analysis. According to the extent of tumor invasion, SEER database classified the tumors into localized, regional and distant metastases. Most patients were diagnosed with distant metastasis in our study, and resulted in the worse survival.

PPS is advanced to the stage IV through local invasion and blood metastasis with few lymph node metastases. But, Regnard et aj.¹⁹ described that 5 out of 20 resected patients reported with lymph node metastasis. Spraker et al¹⁶ delineated that the proportion of patients with lymph node metastasis in a large sample cohort of 365 patients was 16%. The rate of lymph node invasion was between 8% and 22.5% in other studies.^24,26,27 In the present study, 35% of patients reported positive lymph node. In the case of limb sarcomas, metastasis to lymph nodes only occurs in 2% to 6% of patients.^28–30 These findings may indicate that sarcomas are confined in the lungs spread in a different way than the sarcomas in other locations. PPS can attain lymph node metastasis at a higher rate than limb soft tissue sarcoma, therefore, the prognosis of PPS was reported to be worse. This proved the necessity of systematic lymph node dissection, which could not only optimize the stage, but also refine choosing the adjuvant therapy in order to improve the prognosis. Lymph node dissection was an important factor affecting OS in the patients with PPS, which was consistent with the results of this study. We can advise sentinel lymph node biopsy and radical lymphadenectomy for the PPS patients who were clinically positive lymph node.

The treatment strategy of PPS was similar to that of lung cancer. Complete resection with a clear edge could be considered as an important strategy of treatment.^24,25 Early studies described that the survival rate of patients who could undergo complete surgical resection to treat PPS was significantly higher than that of the patients with positive surgical margin or inoperable surgery.^17,23,31 These findings confirmed the results of improved survival in patients who can undergo resection. Gołota et al²⁵ examined the influence of surgical interventions on OS in the patients with PPS. According to these authors, when compared to the patients with non-anatomical surgery, the patients with anatomical complete resection exhibited lower 5-year OS (6% vs 83%), but were less prone to local recurrence. The high recurrence rate and increased lymph node metastasis rate in the patients undergoing non-anatomical complete resection reported that the metastasis ability of PPS could be similar to that of lung cancer. The scope and type of different surgical procedures affected the survival and recurrence of patients.

Surgical treatment can enhance clinical outcomes in the PPS patients, but the role of adjuvant radiotherapy or chemotherapy was still inconclusive due to the chemo/radio-resistance acquired by sarcoma.^8,25,26 Currently, there was no standard chemotherapy regimen for treating PPS, which was mostly carried out with reference to the therapeutic regimen given for treating NSCLC. The effect of chemotherapy on PPS was poor, consistent with the results of this study, and chemotherapy had no significant impact on survival rate (P > .05). There was no sufficient evidence-based data to prove that radiotherapy could improve the prognosis of PPS patients. Some studies reported that radiotherapy exhibited certain value for soft tissue sarcoma and could improve the local control rate, but there were also published data that did not support adjuvant radiotherapy.^32,33 In a large retrospective analysis based on the SEER database, Spraker et al¹⁶ described that the 5-year survival rates of PPS patients receiving surgery alone, radiotherapy alone, surgery and radiotherapy were 41%, 7%, and 25%, respectively. Yeo et al¹² reported the effect of ‘stereotactic ablation radiotherapy’ (SABR) for treating PPS. The results showed that SABR achieved excellent local control with minimal toxicity in both extrapulmonary metastatic sarcoma and PPS. This reflected the survival benefits obtained by using adjuvant radiotherapy to improve local area control. We speculate that different inclusion criteria for patients, different radiotherapy sites or methods in different institutions, and different incidence of radiotherapy complications may affect the outcome of radiotherapy.

In our study, neither radiotherapy nor chemotherapy was acting as an independent influencing factor of PPS. However, after stratification of PPS patients according to the degree of risk, we found that radiotherapy was a risk factor for low-risk population, but adjuvant chemotherapy has not conferred any influence on the patient's OS in low-risk group. In the high-risk population, the adjuvant radiotherapy or chemotherapy could significantly improve the prognosis of PPS patients. This is suggesting that we may choose the appropriate adjuvant treatment strategy based on the degree of risk in the patients. Unfortunately, we are unable to obtain specific chemotherapy drugs from SEER database, and the optimal chemotherapy needs further research. Risk stratification treatment can provide us with a new idea. According to different risk stratification, different treatment modalities can be adopted, including surgical tumor resection, radiotherapy and chemotherapy regimen. High risk group indicates that the tumor has worse biological behavior, which requires advanced surgical strategies and appropriate adjuvant treatment. Surgical treatment is the first choice for low-risk group, and adjuvant treatment is not recommended temporarily; further studies with larger groups of patients are needed to confirm these findings.

This was the first large-scale data obtained through the SEER database to construct nomogram for evaluating the prognosis of PPS patients. The nomogram contains more predictors and is more accurate than the traditional staging model when compared to the traditional TNM staging. Moreover, it is more convenient to evaluate the individual prognosis by using nomogram due to existing clinical data. There are also some limitations for our current study. First, it was a single center retrospective analysis, which may lead to selection bias. Second, this study was limited by the incomplete information of the SEER database and lacking covariate information such as smoking history, positive or negative surgical margin, clinical symptoms, laboratory examination, and specific radiotherapy technology. Third, the bootstrap self sampling method was used for internal verification, and there was a lack of external verification data.

Conclusion

This study concluded the prognostic factors affecting the OS of PPS patients, subsequently constructed and validated the nomogram; this nomogram could provide a practical, convenient and reliable tool for PPS patients to evaluate the prognosis consequently to choose personalized oncomedicine. Radiotherapy or chemotherapy was not recommended for the low-risk individuals, but the individuals at high risk could benefit from adjuvant chemotherapy or radiotherapy.

Footnotes

Abbreviations

Author Contributions

Hao Gu (HG), Ruixia Song (RS), Narasimha M. Beeraka (NMB), Tingxuan Li (TL), Di Zhao (DZ), Junqi Liu (JL), Ruitai Fan (RF) has conceptualized and collected the literature. HG, NMB, RS, TL performed literature search and developed the manuscript. NMB, JL, RF proof-read the manuscript. All authors have reviewed and approved the manuscript before submission.

Declaration of Conflicting Interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: The Natural Science Foundation of Henan Province, (grant number 202300410453).

Informed Consent

Informed consent is not required and not relevant to this study.

Ethics Approval

Consent for Publication

NA.

ORCID iDs

Junqi Liu

Ruitai Fan

Correction (May 2023):

Article updated to link authors who contributed equally with contribution statement.

References

Al Sbihi

Manasrah

Al Haj

Al Qasem

Appel

. A paraneoplastic leukemoid reaction in primary lung sarcoma. Cureus. 2021;13(5):e15047.

Collaud

Stork

Schildhaus

, et al. Multimodality treatment including surgery for primary pulmonary sarcoma: size does matter. J Surg Oncol. 2020;122(3):506‐514.

Suster

. Primary sarcomas of the lung. Semin Diagn Pathol. 1995;12(2):140‐157.

Cakir

Topal

Bayram

Tolunay

. Sarcomas: rare primary malignant tumors of the thorax. Diagn Interv Radiol. 2005;11(1):23‐27.

Etienne-Mastroianni

Falchero

Chalabreysse

, et al. Primary sarcomas of the lung: a clinicopathologic study of 12 cases. Lung Cancer. 2002;38(3):283‐289.

Berzenji

Van Schil

. Commentary: primary pulmonary sarcomas and pulmonary carcinosarcomas, challenging and enigmatic, but treatable. J Thorac Cardiovasc Surg. 2021;162(1):286‐287.

Golota

Osowiecka

Orlowski

. Primary pulmonary sarcoma—long-term treatment outcomes and prognostic factors. Kardiochir Torakochirurgia Pol. 2018;15(3):162‐169.

Sayan

Kankoc

Ozkan

Celik

Kurul

Tastepe

. Prognostic analysis of primary pulmonary malignant mesenchymal tumors treated surgically. J Chest Surg. 2021;54(5):356‐360.

Plestina

Librenjak

Marusic

Batelja Vuletic

Janevski

Jakopovic

. An extremely rare primary sarcoma of the lung with peritoneal and small bowel metastases: a case report. World J Surg Oncol. 2019;17(1):147.

10.

Beasley

Brambilla

Travis

. The 2004 World Health Organization classification of lung tumors. Semin Roentgenol. 2005;40(2):90‐97.

11.

Belbaraka

Ismaili

. Primary sarcoma of the lung: a very rare diagnosis and poor prognosis. Clin Cancer Investig J. 2014;3(2):176.

12.

Yeo

S-G

. Primary lung sarcoma treated with stereotactic ablative radiotherapy: a case report. Oncotargets Ther. 2017;Volume 10:3285‐3288.

13.

Liu

Xie

, et al. Development and validation of a SEER-based prognostic nomogram for cervical cancer patients below the age of 45 years. Bosn J Basic Med Sci. 2021;21(5):620‐631.

14.

Zuo

Zhang

Song

, et al. Survival nomogram for stage IB non-small-cell lung cancer patients, based on the SEER database and an external validation cohort. Ann Surg Oncol. 2021;28(7):3941‐3950.

15.

Huwer

Kalweit

Straub

Feindt

Volkmer

Gams

. Pulmonary carcinosarcoma: diagnostic problems and determinants of the prognosis. Eur J Cardiothorac Surg. 1996;10(6):403‐407.

16.

Spraker

Bair

Connell

Mahmood

Koshy

. An analysis of patient characteristics and clinical outcomes in primary pulmonary sarcoma. J Thorac Oncol. 2013;8(2):147‐151.

17.

Janssen

Mulder

Wagenaar

Elbers

van den Bosch

. Primary sarcoma of the lung: a clinical study with long-term follow-up. Ann Thorac Surg. 1994;58(4):1151‐1155.

18.

Keel

Bacha

Mark

Nielsen

Rosenberg

. Primary pulmonary sarcoma: a clinicopathologic study of 26 cases. Mod Pathol. 1999;12(12):1124‐1131.

19.

Regnard

Icard

Guibert

de Montpreville

Magdeleinat

Levasseur

. Prognostic factors and results after surgical treatment of primary sarcomas of the lung. Ann Thorac Surg. 1999;68(1):227‐231.

20.

Porte

Metois

Leroy

Conti

Gosselin

Wurtz

. Surgical treatment of primary sarcoma of the lung. Eur J Cardiothorac Surg. 2000;18(2):136‐142.

21.

Bacha

Wright

Grillo

, et al. Surgical treatment of primary pulmonary sarcomas. Eur J Cardiothorac Surg. 1999;15(4):456‐460.

22.

Coindre

. [New WHO classification of tumours of soft tissue and bone]. Ann Pathol. 2012;32(5 Suppl):S115‐S116.

23.

Nascimento

Unni

Bernatz

. Sarcomas of the lung. Mayo Clin Proc. 1982;57(6):355‐359.

24.

Petrov

Vlassov

Kalaydjiev

, et al. Primary pulmonary sarcomas and carcinosarcomas--postoperative results and comparative survival analysis. Eur J Cardiothorac Surg. 2003;23(4):461‐466.

25.

Golota

Osowiecka

Orlowski

. Primary pulmonary sarcoma—treatment outcomes depending on the different types of radical operation. Kardiochir Torakochirurgia Pol. 2019;16(1):1‐6.

26.

Yamada

Kaplan

Soltermann

, et al. Surgical outcomes and risk analysis of primary pulmonary sarcoma. Thorac Cardiovasc Surg. 2021;69(1):101‐108.

27.

Robinson

Babacan

Tanvetyanon

Henderson-Jackson

Bui

Druta

. Results of treating primary pulmonary sarcomas and pulmonary carcinosarcomas. J Thorac Cardiovasc Surg. 2021;162(1):274‐284.

28.

Fong

Coit

Woodruff

Brennan

. Lymph node metastasis from soft tissue sarcoma in adults. Analysis of data from a prospective database of 1772 sarcoma patients. Ann Surg. 1993;217(1):72‐77.

29.

Riad

Griffin

Liberman

, et al. Lymph node metastasis in soft tissue sarcoma in an extremity. Clin Orthop Relat Res. 2004;426:129‐134.

30.

Mazeron

Suit

. Lymph nodes as sites of metastases from sarcomas of soft tissue. Cancer. 1987;60(8):1800‐1808.

31.

Boffa

Rosen

Mallin

, et al. Using the National Cancer Database for outcomes research: a review. JAMA Oncol. 2017;3(12):1722‐1728.

32.

Cosgun

Tezel

Akyil

Kolbas

Sen

Tezel

. Primary pulmonary malignant fibrous histiocytoma. Turk Thorac J. 2017;18(2):54‐56.

33.

Zagars

Mullen

Pollack

. Malignant fibrous histiocytoma: outcome and prognostic factors following conservation surgery and radiotherapy. Int J Radiat Oncol Biol Phys. 1996;34(5):983‐994.