Abstract

Objective. To establish and validate a nomogram to predict the overall survival (OS) and cancer-specific survival (CSS) in patients with nasopharyngeal carcinoma (NPC) receiving radiotherapy by integrating multiple independent prognostic factors. Materials and Methods. Data from 5663 patients with NPC who received definite radiotherapy between 2004 and 2018 were included and divided into training and validation cohorts. Univariate and multivariate Cox regression analyses were performed to determine the independent prognostic factors of patients with NPC after radiotherapy. Thereafter, the predictive accuracy of the nomogram model was evaluated. Results. Age, race, marital status, pathological type, tumor size, T stage, N stage, M stage, American Joint Committee on Cancer stage, and chemotherapy were independent factors affecting the prognosis of patients with NPC receiving radiotherapy. Nomograms with a concordance index of 0.726 (95% confidence interval (CI): 0.675–0.777) and 0.732 (95% CI: 0.680–0.785) were able to predict OS and CSS, respectively. The area under the curve showed excellent predictive performance. Additionally, the calibration curve indicated that the predicted survival rate was consistent with the actual survival rate, and the decision curve indicated its clinical value. The established risk stratification system was able to accurately stratify patients receiving radiotherapy for NPC into three risk subgroups with significant differences in prognosis (). Conclusions. The constructed nomogram had good prognostic performance and could be used as an effective tool to evaluate the prognosis of patients with NPC after radiotherapy. This nomogram could be further used to guide clinical decisions and personalized treatment plans.

1. Introduction

Nasopharyngeal carcinoma (NPC) is one of the most common head and neck malignancies, which occurs mainly in Guangdong, Hainan, Hunan, and other regions in China [1, 2]. NPC is sensitive to radiotherapy, which is the first choice for radical treatment, as alone it can achieve a good prognosis for early lesions [3, 4]. However, the initial symptoms of NPC are often unnoticed, and regional lymph nodes and/or blood vessel metastasis can occur easily. The majority of NPC patients are already in the locally advanced stage when they are first admitted to a hospital, and the efficacy of radiotherapy alone is not favorable. The 5-year survival rate for NPC patients has been reported to range between 50% and 60%, with the main reasons for treatment failure being local recurrence and distant metastasis [57]. With the development of imaging and intensity-modulated radiotherapy technology and the application of radiochemical combined therapies, the efficacy of NPC treatment has significantly improved, with the 5-year local control and survival rates being approximately 90% and 80%, respectively [8]. The factors influencing the prognosis of patients with NPC after radiotherapy mainly include the general health of patients, disease factors, and treatment factors [5, 7]. Individual prognostic prediction is important for NPC patients undergoing radiotherapy. Traditional staging systems, such as the American Joint Committee on Cancer (AJCC) staging classification [9], are generally applicable to patient groups rather than to specific individuals. Clinical guidelines for identifying clinical risk are usually based on the tumor node metastases (TNM) staging at diagnosis. Additionally, TNM staging only considers tumor size, lymph node, and distant metastases, regardless of age, sex, tumor grade, treatment, and other prognostic factors [10, 11]. Therefore, there is an urgent need to develop a more accurate analysis tool to predict the individual prognosis of NPC patients after radiotherapy and to accurately stratify patients with different risk profiles [12, 13].

In recent years, various studies have used nomograms to construct prediction models as a tool for oncology and medical prognosis evaluation. Nomograms can generate individual predictions by integrating multiple prognostic and determining variables for personalized patient identification and stratification [13, 14]. The Surveillance, Epidemiology, and End Results (SEER) database is the largest publicly available tumor registry database in the United States (US), which includes demographic, clinicopathological, survival, and clinical data of patients with various types of tumors, covering approximately 30% of the US population [15, 16]. The SEER database collates data from 18 cancer registry databases [15]. The aim of this study was to use the SEER database to screen patients with NPC and randomly divide them into training and validation cohorts to construct a new risk stratification system based on the individual prediction scores of the nomogram. Additionally, it is aimed at predicting overall survival (OS) and cancer-specific survival (CSS) and at evaluating the performance of the statistical prediction model. In doing so, we hoped to identify high-risk patients and guide clinical treatment decisions for NPC.

2. Materials and Methods

2.1. Data Sources and Study Populations

This was a retrospective study where information on patients diagnosed with NPC between 2004 and 2018 was obtained from the SEER database using the SEERStat software (version 8.3.9, http://www.seer.cancer.gov/seerstat). The inclusion criteria were as follows: (1) patients were diagnosed between 2004 and 2018; (2) their tumors were staged according to the 6th, 7th, and 8th versions of AJCC (2004–2009, 2010–2017, and 2018, respectively) and were available in the SEER database; (3) they had newly diagnosed NPC without a history of surgery; (4) they had received definite radiotherapy; and (5) the International Classification of Diseases for Oncology, 3rd Edition (ICD-O-3), was used to classify tumors according to histological types determined by the World Health Organization (WHO) classification scheme. Squamous cell carcinoma (SCC) (ICD-O-3 codes 8070 and 8071) represents the histological subtypes of keratinizing squamous cell carcinoma (KSCC), while nonkeratinizing carcinoma (NKC) (ICD-O-3 codes 8072 and 8073) represents the histological subtypes of differentiated NKC. Undifferentiated, anaplastic, and lymphoepithelial carcinomas (ICD-O-3 codes 8020, 8021, and 8082) represent undifferentiated NKC histologic subtypes. Carcinoma that was not specified (ICD-O-3 code 8010) was assigned to form an NOS histologic subtype. For such cases, the causes of death and survival time were known.

Patients were excluded if (1) the diagnostic method was unknown; (2) their race and/or marital status was unknown; (3) their cause of death was unknown; (4) they had an unknown or incomplete TNM staging and/or an unknown metastasis; (5) their radiotherapy status was unknown; (6) their available follow-up time was less than 1 month; and (7) they received nonbeam radiation. Beam radiotherapy was defined as external beam radiation using gamma rays generated from a cobalt-60 source or high-energy photons generated from a linear accelerator. Nonbeam radiation was defined as any type of brachytherapy, unsealed source, intraoperative radiation not otherwise specified, or radiation not otherwise specified. This group of patients was excluded to reduce potential selection bias caused by different radiotherapy modalities.

2.2. Patient Data Collection

Data from 8599 patients with NPC was collected from the SEER database. The data processing flowchart is shown in Figure 1. Overall, 5663 patients with NPC were included in this study. All patients were randomly assigned to training () and validation () cohorts using a 7 : 3 ratio. The variables included for each patient were age, sex, race, year of diagnosis, marital status, pathological type, histological grade, tumor size, TNM stage, AJCC stage, chemotherapy record, diagnostic confirmation record, cause-specific death classification, follow-up time, and survival status. The main endpoints for this study were OS and CSS. OS was defined as the time between the date of NPC diagnosis and the date of death from any cause or the last follow-up date. CSS was defined as the time from the date of diagnosis to death attributable to NPC.

2.3. Statistical Methods

Statistical analysis was conducted using SPSS version 25.0 and R version 4.1.2. The R function “createDataPartition” was used to divide the patients into training and validation cohorts using a 7 : 3 ratio to ensure randomization. The training cohort was used to filter variables and build prediction models. Variables that were significant in the univariate Cox regression model were included in the multivariate Cox regression model, and multivariate proportional hazard models (hazard ratios and 95% confidence intervals (CI)) were used to identify independent prognostic factors associated with OS and CSS. The prognostic nomogram was constructed based on the multivariate analysis of the training cohort, and it was used to predict the 3- and 5-year OS and CSS by representing the total points for each variable [17]. Individual scores were read based on each clinicopathological indicator for each patient, and the scores were added together to obtain a total score. The lower the total score, the higher the 3- and 5-year survival rates. The accuracy of the prognostic model was evaluated by calculating the concordance index (C-index) and the receiver operating characteristic (ROC) curve. The time-dependent area under the curve (AUC) was not only used to evaluate the predictive power of the nomogram but was also used to compare the predictive ability of the AJCC stage. The calibration curve was used to compare the consistency of the predicted OS and CSS with the actual survival rate, and the decision curve was used to evaluate the clinical applicability of this nomogram by quantifying the net improvement benefit at different threshold probabilities. The risk stratification system was established based on the scores of each patient in the training cohort, and the optimal cut-off values of continuous variables were determined using the X-tile software. The patients were divided into three subgroups, namely, low-, medium-, and high-risk groups. Each risk subgroup represented a different prognosis. The Kaplan-Meier survival analysis was used to generate survival curves for the different risk groups, and log-rank tests were used for comparison between the groups. was considered statistically significant.

3. Results

3.1. General Information

In total, 5663 screened patients with NPC had a median follow-up of 44 months. Overall, 3967 and 1696 patients were randomly divided into the training and validation cohorts, of whom 1556 and 673 died, respectively. Clinical and demographic information for the two groups is shown in Table 1; there were no significant differences between the two groups ().

3.2. Independent Prognostic Factors in the Training Cohort

Univariate Cox regression analysis showed that there were significant differences in OS and CSS for the evaluated prognostic factors () with the exception of sex (Table 2). In the multivariate Cox regression analysis of OS and CSS, age, race, marital status, pathological type, tumor size, T stage, N stage, M stage, AJCC stage, and chemotherapy were independent factors affecting OS and CSS in NPC patients undergoing radiotherapy (all ) (Table 3).

3.3. Construction and Validation of a Nomogram for Predicting the Total Survival Rate

Based on independent prognostic factors determined in the multivariate Cox regression analysis, a nomogram was constructed to predict 3- and 5-year OS and CSS for the training cohort (Figure 2). Age, race, and M stage had greater efficacy in the construction of the prediction model. The scores of each variable were added together to obtain 3- and 5-year OS and CSS. The C-index for OS was higher than that of the AJCC staging system (0.726 vs. 0.603), and that of the CSS was higher than that of OS (0.732 vs. 0.726). The AUC of the 3- and 5-year survival rates in the OS prediction model of the training cohort were 0.743 and 0.738, respectively, and those in the CSS prediction model were 0.751 and 0.734, respectively. The curve in Figure 3 shows that the AUC value for the OS/CSS nomogram was higher than that of the AJCC stage and separate T, N, and M stages, and the predictive ability of the nomogram was better. Additionally, the calibration curves show that the nomogram’s prediction of the 3- and 5-year OS and CSS was in agreement with the actual survival rates (Figure 4). The decision curve for the nomogram and AJCC staging prediction model is shown in Figure 5. The net benefit was higher than that of the AJCC model in OS and CSS prediction at 3 and 5 years.

The Kaplan-Meier survival curve for the prognostic risk stratification system found that OS and CSS in the three subgroups were accurately separated by the system, and there was a significant difference in survival among the three groups () (Figure 6).

4. Discussion

Radiotherapy is still the cornerstone of NPC treatment [18, 19]. In the past two decades, there have been significant changes in the standard of care for NPC, including the use of the AJCC/Union for International Cancer Control (UICC) staging system [20, 21], and the use of magnetic resonance imaging (MRI) as standard radiologic staging for local disease. Cisplatin chemotherapy is used in combination with radiotherapy for patients with locally advanced disease, and intensity-modulated radiotherapy (IMRT) is used as standard radiotherapy [22].

The AJCC/UICC TNM staging system is the most widely accepted universal language for describing tumors, is the most common staging system for NPC, and is an important basis for guiding clinical activity in judging prognosis [2325]. This staging system is continually updated, with the latest edition (8th edition) relying on improved local control of IMRT, as well as more detailed MRI preprocessing imaging, which results in a more consistent definition of T2 (with T4) and N3 diseases [26]. However, AJCC stage is not the only factor affecting the prognosis of NPC, and the survival time of NPC patients with the same stage varies. A growing number of studies have incorporated clinicopathologic features and clinical treatment information, including demographic statistical variables, Epstein-Barr virus- (EBV-) DNA, hematological markers of inflammation, and gross tumor volume, to further improve the prognostic accuracy [27].

NPC is sensitive to radiotherapy and chemotherapy. The combination of chemotherapy and radiotherapy is considered a breakthrough in the treatment of locally advanced NPC. A number of randomized trials have shown relatively consistent results, with chemotherapy plus radiotherapy having better survival correlates [28, 29]. Au et al. demonstrated that the use of chemotherapy (neoadjuvant plus concurrent, concurrent, and adjuvant plus concurrent) was an independent favorable prognostic factor for OS and PFS compared with IMRT alone, confirming the benefit of chemotherapy for IMRT [22]. The presence of plasma EBV DNA in patients with NPC is associated with apoptosis of tumor tissue and has the same polymorphism as the primary site, and plasma EBV DNA is strongly associated with tumor grade. Consequently, it has become the most accurate biomarker for screening, diagnosis, monitoring response to treatment, monitoring recurrence, and prediction of NPC in endemic areas [30, 31]. Two meta-analyses have also shown that pretreatment plasma EBV DNA titers are an important prognostic factor in NPC [32, 33]. In 1978, the WHO classified NPC histologically into three categories: squamous cell carcinoma (WHO-I), NKC (WHO-II), and undifferentiated carcinoma (WHO-III). This WHO classification was revised in 1991, and NPCs were classified as KSCC and NKC. NKC was then further subdivided into differentiated and undifferentiated. The third edition of the WHO classification of NPC added basal cell squamous cell carcinoma as one of the histological subtypes in 2003 [34]. Histological subtypes of NPC have a distinct geographic and ethnic distribution in both endemic and nonendemic areas of NPC, and in endemic areas, such as southern China, more than 95% of NPC patients have NKC; however, KSCC is more common in nonendemic areas [22, 27]. Pan et al. showed that histological subtypes determine long-term survival outcomes in NPC patients. Long-term survival outcomes are worse in patients with KSCC, whereas survival is higher in patients with NKC [35].

In this study, age, race, marital status, pathologic type, tumor size, T stage, N stage, M stage, AJCC stage, and chemotherapy were identified as independent prognostic variables for OS. The nomogram is based on multivariate regression analysis, which calculates the individual end event forecast value using the individual score and end event occurrence probability function transformation. Previous studies have shown that complex regression equations can be transformed into visual graphs to make the results of the prediction model more readable and to facilitate the evaluation of patients [36, 37]. Traditionally, AJCC staging, which is closely related to OS, has been the first choice for the prediction of NPC prognosis. However, the prognosis of patients with the same staging can differ. This heterogeneity occurs because AJCC staging does not consider other factors such as age, pathology, adjuvant therapy, marital status, and/or race [38, 39]. In predicting the prognosis of NPC, the nomogram reduces the diversity due to different treatments and sociodemographic and treatment statuses. However, there is still a lack of a complete and effective nomogram prediction model for NPC patients undergoing radiotherapy, and only a few studies have applied it to visual prediction models [40, 41]. Moreover, the stability and generality of the alignment chart may be reduced as a result of the smaller sample size.

Our study developed a 3- and 5-year OS and CSS nomogram for NPC patients based on 5663 samples obtained from the SEER database. Among them, M staging, age, and race contributed the most substantially to the OS nomogram, while M staging, age, and AJCC staging contributed the most substantially to the CSS nomogram. Although we used different versions of AJCC staging system, it was helpful to analyze the results and showed a good comparability. By comparing the C-index of the nomogram and AJCC stages, the C-index of the nomogram for OS and CSS prediction was found to be higher than that of AJCC in the training and validation cohorts. These findings suggest that this model has a good recognition ability, and the C-index is better than the AJCC staging prediction model. The correction curve showed that the predicted probability of the nomogram was in agreement with the actual probability calculated using the Kaplan-Meier analysis, and the 3- and 5-year survival rates were in good agreement. These results highlight that this nomogram could be used in clinical practice with less deviation and better accuracy. Decision curve analysis combines benefits and harms to measure the net benefits of an approach or predictive model. Compared with the ROC curve, the decision curve considers the clinical application value, which is an important indicator to evaluate whether patients can benefit from it. However, only a few papers have applied it to the survival prediction for NPC. In our study, we calculated not only the net benefit of the nomogram but also evaluated the AJCC staging prediction model and found that the net benefit of the nomogram was higher than that of the AJCC staging prediction model, indicating that this nomogram is feasible for clinical use.

Additionally, we established a prognostic nomogram risk stratification system to accurately classify patients into three risk subgroups and identify high-risk subgroups that may require more intensive treatment. This provided guidance for patient counseling and clinical risk management. For high-risk patients, more psychological and emotional care may also be needed. Additionally, they may be encouraged to participate in clinical trials evaluating new drugs, such as immunotherapy and targeted drugs. The use of this model for stratification helped to reduce heterogeneity among the different treatment groups and thus should be considered for future use.

This study has some limitations which merit mentioning. First, it was a retrospective study, in which patient selection was based on rigorous inclusion and exclusion criteria; thus, potential selection bias may have occurred. Second, the SEER database does not specify procedures, operators, and other factors, which may induce bias due to the different experience levels of the operators and pathologists. Third, because of the limited number of variables collected in the SEER database, some important variables were not available, such as Epstein-Barr virus (EBV) infection, complications, detailed chemotherapy regimens, chemotherapy course and duration, AJCC integration, and staging information based on MRI images. Fourth, SEER does not mention tumor location (primary or metastatic) or dose and mode of radiotherapy, and it does not specify the specific radiotherapy technique, making the study of the survival benefit of patients treated with radiotherapy somewhat limited. Finally, in the absence of external validation, analysis of US population-based data is not necessarily generalizable to the global population, especially to patients with NPC in endemic areas. Therefore, to improve the predictive value of the model, it will be important to consider integrating multicenter Chinese and foreign databases in the future.

5. Conclusions

Based on clinical features and treatment-related information evaluated, our study established a useful nomogram for the prediction of 3- and 5-year OS and CSS in patients with NPC receiving radiotherapy. The results showed that our nomogram was superior to the AJCC staging prediction model and was successful in predicting individual survival. Further external validation using prospective data should be conducted to verify these findings. Nonetheless, this nomogram could be used in the future to assist clinicians in devising the best treatment plans for patients with NPC.

Data Availability

If the original data is reasonably required, the permission to obtain the data can be applied to Professor Zhiru Li.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Authors’ Contributions

Zhiru Li and Chao Li contributed equally to this work.

Acknowledgments

We would like to acknowledge all the participants for their participation and their contribution to this study. We also acknowledge that SEER database provided the data for this study. This study was supported by grants from the National Natural Science Foundation of China (Nos. 71964003, 81460460, 81760542, and 82160467), the Natural Science Foundation of Guangxi Zhuang Autonomous Region (No. 2018JJA141048), the Research Foundation of the Science and Technology Department of Guangxi Province, China (grant Nos. 2016GXNSFAA380252, 2018AB61001, and 2014GXNSFBA118114), the Research Foundation of the Health Department of Guangxi Province, China (No. S2018087), the Guangxi Medical University Training Program for Distinguished Young Scholars (2017), the Medical Excellence Award funded by the Creative Research Development Grant from the First Affiliated Hospital of Guangxi Medical University (2016), the Guangxi Medical High-Level Talents Training Program, and the Central Government Guide Local Science and Technology Development Projects (ZY18057006).