Abstract

Background. The primary purpose of this study was to determine the risk factors affecting overall survival (OS) in patients with fibrosarcoma after surgery and to develop a prognostic nomogram in these patients. Methods. Data were collected from the Surveillance, Epidemiology, and End Results database on 439 postoperative patients with fibrosarcoma who underwent surgical resection from 2004 to 2015. Independent risk factors were identified by performing Cox regression analysis on the training set, and based on this, a prognostic nomogram was created. The accuracy of the prognostic model in terms of survival was demonstrated by the area under the curve (AUC) of the receiver operating characteristic curves. In addition, the prediction consistency and clinical value of the nomogram were validated by calibration curves and decision curve analysis. Results. All included patients were divided into a training set (n = 308) and a validation set (n = 131). Based on univariate and multivariate analyses, we determined that age, race, grade, and historic stage were independent risk factors for overall survival after surgery in patients with fibrosarcoma. The AUC of the receiver operating characteristic curves demonstrated the high predictive accuracy of the prognostic nomogram, while the decision curve analysis revealed the high clinical application of the model. The calibration curves showed good agreement between predicted and observed survival rates. Conclusion. We developed a new nomogram to estimate 1-year, 3-year, and 5-year OS based on the independent risk factors. The model has good discriminatory performance and calibration ability for predicting the prognosis of patients with fibrosarcoma after surgery.

1. Introduction

Fibrosarcoma is a group of malignant soft tissue tumors. As solid tumors composed of fibroblasts, FS often involve the deep soft tissues of the extremities, trunk, head, and neck. Surgery is the primary treatment of fibrosarcoma, but the prognosis of its patients varies depending on the type of tumor [1].

Adult fibrosarcoma used to be the most common soft tissue sarcoma. However, with the increased awareness of soft tissue tumors and fibrosarcoma, the diagnosis of adult fibrosarcoma has been tightly defined so that its incidence may be as low as 1% of adult soft tissue sarcomas [2, 3]. Unlike the highly malignant adult fibrosarcoma, congenital or infantile fibrosarcoma is intermediate in 2020 WHO classification, with few metastases and a better prognosis [46]. Since neurotrophic tropomyosin receptor kinase (NTRK) fusion often occurs in infantile fibrosarcoma, the application of NTRK inhibitors such as larotrectinib is highly efficacious in patients [79]. Dermatofibrosarcoma protuberans (DFSP) is a low-grade, growing by infiltration but highly locally recurrent fibroblastic malignancy. DFSP most commonly occurs in middle-aged and young adults. Through surgery, radiation therapy, or the targeted drug imatinib mesylate, patients with DFSP have a good prognosis and a high survival rate [1012]. Surgical resection is the standard treatment for most patients, whether it is a typical or a particular type of fibrosarcoma. Therefore, the evaluation of postoperative prognostic indicators of fibrosarcoma will provide valuable guidance for the clinical treatment of fibrosarcoma.

Nomograms integrate multiple risk indicators and are often used to predict disease survival [1315]. Fibrosarcoma is primarily treated surgically but has a high recurrence rate after simple resection [2, 16]. Recently, Xiang et al. constructed a nomogram of 663 cases of fibrosarcoma and determined that age, sex, surgical use, tumor stage, pathologic grade, and tumor size may affect cancer-specific survival [15]. However, studies related to predicting overall survival (OS) after surgery in patients with fibrosarcoma have not been reported. This study aimed to develop a nomogram to identify the clinical and pathological factors associated with improved OS in fibrosarcoma patients after surgery using the Surveillance, Epidemiology, and End Results (SEER) database.

2. Methods

2.1. Inclusion and Exclusion Criteria of Patients

The SEER database is updated annually with the latest cancer information. We identified patients with fibrosarcoma from 2004 to 2015 by the SEER  Stat version 8.3.9. The inclusion criteria were as follows: (1) histology ICD-O-3 was limited to 8810/3, 8812/3, 8813/3, 8814/3, 8823/3, 8832/3, 8833/3, 9321/3, and 9330/3; (2) fibrosarcoma as the patient’s primary tumor. The exclusion criteria were as follows: (1) the detailed information lacks age, race, grade, primary site, tumor size, and marital status; (2) unknown historic stage and information related to treatment (radiotherapy and chemotherapy). The final dataset we screened had 439 patients with fibrosarcoma. All patients were randomly divided into a training set (70%) and a validation set (30%) for constructing and validating the nomogram.

2.2. Variable Declaration

The age was regrouped as <60 years old (young) and ≥60 years old (old). The race was categorized as white, black, and others. The primary site of the tumor was classified as the head and neck, trunk, and extremities. Tumor size was reclassified as <50 and ≥50 mm. Information on radiotherapy was divided into no radiotherapy, preoperative radiotherapy, and postoperative radiotherapy.

2.3. Statistical Analysis

The primary endpoint of this study, overall survival, was analyzed by Kaplan–Meier survival curves, and a log-rank test for significance was performed [14]. Prognostic factors with a value <0.05 in univariate Cox regression analysis were included in multivariate Cox regression analysis. In multivariate Cox regression analysis, values <0.05 were considered statistically significant, indicating that the variables examined were independent risk factors. In addition, hazard ratios and 95% confidence intervals were used to adjust prognostic variables.

Based on the above results, we plotted the nomogram against the training set by the rms package. The discrimination of the nomogram was demonstrated by the area under the curve (AUC) of the time-dependent receiver operating characteristic (ROC) curves. In addition, 1-year, 3-year, and 5-year calibration curves and decision curve analysis (DCA) were generated to validate the nomogram’s predictive consistency and clinical value. Through the risk assessment system established by X-tile software, the patients with fibrosarcoma after surgery were classified into high-risk, intermediate-risk, and low-risk groups, further demonstrating the application value of the predictive model. All statistics were analyzed by R software (version 4.0.3), while <0.05 (two-sided) was considered significant.

3. Results

3.1. Baseline Characteristics of the Patients

A total of 439 patients with fibrosarcoma after surgery from the SEER database were included according to our criteria. In addition, 308 patients were included in the training set, and 131 patients were included in the validation set. As given in Table 1, the proportion of male patients was higher than that of females (54.5% vs. 45.5%; ). The median age of the primary diagnosis was 52 (25–75%, 37–65). Most patients were white (67.9%), and 65.3% were <60 years old. 7.5% of the population had primary disease sites in the head and neck, 43.2% in the trunk, and 49.4% in the extremities. 35.1% of screened patients with postoperative fibrosarcoma were in stage I, 31.5% in stage II, 17.5% in stage III, and 15.9% in stage IV. 9.7% of the patients received chemotherapy, and 30.5% received radiation therapy, of which 3.6% were preoperative and 26.9% were postoperative. The majority of patients were localized (75.6%), and most tumors were ≥50 mm (55.5%).

3.2. Risk Factors in Postoperative Prognosis of Fibrosarcoma

To identify factors that may predict the occurrence of postoperative fibrosarcoma, univariate and multivariate analyses were performed on the training set. As given in Table 2, the age, race, grade, chemotherapy, tumor size, and historic stage were strongly associated with the OS of the patients with fibrosarcoma after surgery. Consistent with univariate analysis, Kaplan–Meier analysis also showed that clinical factors (age, race, grade, chemotherapy, tumor size, and historic stage) were significantly associated with OS (Figure 1). The multivariate Cox regression analysis confirmed that age, race, grade, and historic stage were independent risk factors (Table 2).

3.3. Development and Validation of Nomogram

We integrated multiple predictors to make a nomogram based on Cox regression analysis to express the interrelationship between the variables (Figure 2). The tumor grade of fibrosarcoma has the most significant impact on patients, and as tumor grade increases, patients have a worse prognosis. Populations of other races have a better prognosis than those of black or white races. Patients who are ≥60 years old or whose tumors have distant metastases have a worse prognosis.

Next, we examine the prognostic models for the training and validation sets, respectively, by plotting ROC curves (Figure 3). The results have demonstrated that the AUC of the training set is 0.827, 0.812, and 0.814 for 12 months, 36 months, and 60 months, respectively (Figure 3(a)), while the validation set is 0.834, 0.817, and 0.766, respectively (Figure 3(b)). The DCA can be used to evaluate the net benefit of nomogram-aided decision-making under different threshold probabilities to assess the benefit degree of patients and the clinical application value of the mode [17]. As shown in Figure 4, this nomogram model has good clinical application in predicting the prognosis of fibrosarcoma after surgery. In addition, the calibration curves for both the training and validation sets have identified a strong agreement between the nomogram prediction and realistic observation for 1-year, 3-year, and 5-year of OS (Figure 5).

3.4. Risk Assessment System

The total score for each patient was calculated, and all patients were divided into three groups using X-tile software, including high-risk (>231), intermediate-risk (205–231), and low-risk (<205) groups. Kaplan–Meier survival curves for each risk subgroup were plotted, and the results in both the training and validation sets showed differences in OS for patients with different risk levels (Figure 6). While patients with high-risk scores had the lowest survival rate, those with low-risk scores had the highest survival rate, indicating that this risk grouping system has a strong predictive value for the prognosis of patients with fibrosarcoma after surgery.

4. Discussion

SEER is a typical medical database that provides systematic evidence support and valuable first-hand information for clinicians’ practice and medical research [13, 15, 1820]. Nomogram is conducive to the promotion of personalized medicine and has been proposed to improve disease prediction [2123]. Wang et al. developed a prognostic nomogram model for small-cell lung cancer patients. They verified that the performance of nomogram was better than that of early models, including those using AJCC staging [24]. As a pictorial representation of a complex mathematical formula, nomograms respond to our exploration of comprehensive biological and clinical models [17, 22, 23].

Fibrosarcoma is a malignant tumor derived from mesenchymal cells, which is usually invasive and has a high postoperative recurrence rate [16]. Adult fibrosarcoma occurs most frequently in the trunk and limbs, followed by the head and neck, consistent with our results [25]. Tumor size of malignant fibrous neoplasms (MFN) of long bones >10 cm is a poor prognostic factor for OS and cancer-specific survival [26]. A report on primary intracranial fibrosarcoma suggested that large tumor volume (≥5 cm) and high Ki-67 index (≥30%) were independent risk factors for OS [27, 28]. Like most types of fibrosarcoma, renal fibrosarcoma is highly malignant and has a poor prognosis. Radical nephrectomy is the primary clinical treatment strategy [29, 30]. Mitotic activity and Ki-67 positive cells were identified as important factors in diagnosing ovarian fibrosarcoma [31, 32]. Less than 5% of fibrosarcoma originated in the genitourinary tract, and fibrosarcoma of the penis is even rarer. To our knowledge, there are few reports on the prognosis of patients with fibrosarcoma. Xiang et al. demonstrated that age, gender, operation, tumor stage, pathological grade, and tumor size might be independent risk factors affecting the survival of FS patients, among which age is the main factor affecting the prognosis [15]. Surgery is the standard treatment for fibrosarcoma at various sites, but the prognosis is abysmal [33]. A reported case of cavernous fibrosarcoma died ten months after the operation, although he received chemotherapy [34]. Patel et al. analyzed 51 sinonasal fibrosarcoma patients from the SEER database. They found that the disease-specific survival rate of patients treated with surgery was better than patients treated with primary radiotherapy alone [35, 36]. However, there has been no evaluation of postoperative prognosis in patients with fibrosarcoma.

This study aims to develop a clinically helpful nomogram that could predict the prognosis of patients with fibrosarcoma after surgery. We finally included 439 patients with postsurgical fibrosarcoma from the SEER database. A randomly divided training set was used to establish the prediction model, while the validation set was used to test the model’s accuracy. Through univariate and multivariate analyses, we determined that age, race, tumor grade, and historic stage were closely associated with prognosis (Table 2). We did not find that chemotherapy was an independent prognostic factor for OS, and the results were consistent with a recent study [26]. Based on this, a nomogram was created to facilitate clinical work (Figure 2). Validation is the process of testing models on different groups to obtain unbiased estimates of model performance and judge their applicability to these groups [17]. Knowledge of the identification and calibration of nomograms in this patient group will enable clinicians and patients to comprehensively evaluate the reliability and accuracy of nomograms [17, 37]. Validated by calibration curves, ROC curves, and DCA, our prognostic model showed good internal and external performance in predicting the prognosis of patients with fibrosarcoma after surgery (Figures 35).

This study may have limitations due to its retrospective nature. Although the nomogram has good performance, it ignores patients’ and doctors’ satisfaction and quality of life and lacks clinical utility [17, 38]. Therefore, all nomograms require more validation using large independent groups. Nomograms must be rigorously reviewed before they can be used in clinical decision-making. When we understand the performance and limitations of the predictive model, we can provide a better prognosis for patients.

5. Conclusion

In this study, we used the SEER database to analyze the prognosis data of patients with fibrosarcoma after surgery for the first time. Nomograms for estimating OS at 1, 3, and 5 years were established based on a large study cohort. The current model has good prediction ability for patient diagnosis, risk assessment, and clinical decision-making to help clinicians provide highly customized patient management in the future.

Abbreviations

NTRK:Neurotrophic tropomyosin receptor kinase
DFSP:Dermatofibrosarcoma protuberans
OS:Overall survival
SEER:Surveillance, Epidemiology, and End Results
AUC:Area under the curve
ROC:Receiver operating characteristic
DCA:Decision curve analysis.

Data Availability

The dataset from the SEER database that was generated and/or analyzed during the current study is available in the SEER dataset repository (https://seer.cancer.gov/).

Ethical Approval

The data were collected from the SEER database, and due to its retrospective nature, the study was exempted by the SEER database administrators.

Not applicable.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

YY C, CX C, and ZP Z conceived and designed the study. SJ C and BB Y performed the literature search. YY C and SJ H generated the figures and tables. YY Cand CX C analyzed the data and wrote the manuscript. ZP Z critically reviewed the manuscript. YY C and ZP Z supervised the research. Yuyuan Chen, Changxing Chi, and Dedian Chen contributed equally to this work. All authors have read and approved the manuscript.