Skip to main content
  • Research article
  • Open access
  • Published:

Baseline MRI-based radiomics model assisted predicting disease progression in nasopharyngeal carcinoma patients with complete response after treatment



Accurate pretreatment prediction for disease progression of nasopharyngeal carcinoma is key to intensify therapeutic strategies to high-risk individuals. Our aim was to evaluate the value of baseline MRI-based radiomics machine-learning models in predicting the disease progression in nasopharyngeal carcinoma patients who achieved complete response after treatment.


In this retrospective study, 171 patients with pathologically confirmed nasopharyngeal carcinoma were included. Using hold-out cross validation scheme (7:3), relevant radiomic features were selected with the least absolute shrinkage and selection operator method based on baseline T2-weighted fat suppression and contrast-enhanced T1-weighted images in the training cohort. After Pearson’s correlation analysis of selected radiomic features, multivariate logistic regression analysis was applied to radiomic features and clinical characteristics selection. Logistic regression analysis and support vector machine classifier were utilized to build the predictive model respectively. The predictive accuracy of the model was evaluated by ROC analysis along with sensitivity, specificity and AUC calculated in the validation cohort.


A prediction model using logistic regression analysis comprising 4 radiomics features (HGLZE_T2H, HGLZE_T1, LDLGLE_T1, and GLNU_T1) and 5 clinical features (histology, T stage, N stage, smoking history, and age) showed the best performance with an AUC of 0.75 in the training cohort (95% CI: 0.66–0.83) and 0.77 in the validation cohort (95% CI: 0.64–0.90). The nine independent impact factors were entered into the nomogram. The calibration curves for probability of 3-year disease progression showed good agreement. The features of this prediction model showed satisfactory clinical utility with decision curve analysis.


A radiomics model derived from pretreatment MR showed good performance for predicting disease progression in nasopharyngeal carcinoma and may help to improve clinical decision making.


Nasopharyngeal carcinoma (NPC), the most common tumor in the nasopharynx, prevails in east and southeast Asia, with an incidence of up to an age-standardized rate of 3.0 per 100,000 in China [1]. Intensity modulated radiotherapy is now the recommended standard treatment for non-metastatic NPC [2]. Despite proven effectiveness of radiotherapy and chemotherapy, treatment failure due to locoregional recurrence or distant metastasis occurs in nearly 10–15% of patients during the first 2 years after tumor remission [3]. It is vital to identify those NPC patients with high risk of disease progression after remission, so as to individualize treatment plan and thus better manage NPC.

Tumor-node-metastasis (TNM) staging system has been widely used to predict prognosis of patients with NPC. However, NPCs with same TNM stage can have completely different responses to chemoradiotherapy and prognosis [4], part of which may be attributed to the fact that TNM staging system mainly reflects the relationship between tumor and surrounding anatomical structures and ignores intra-tumor characteristics, including tumor morphology itself, morphological heterogeneity, etc. [5].

Conventional radiological imaging modalities, such as CT, MRI and PET-CT are widely used in the evaluation of primary tumor extension and tumor staging of NPC. These multi-modality images contain more than just anatomical information about the primary NPC lesion and its surrounding structures, and they also can reflect intra-tumor characteristics, including morphological heterogeneity that is difficult to be quantified using conventionally visual-based reading strategy [6]. As a result, without sufficient data mining, the prognosis value of these medical images was not fully utilized. Radiomics, as an emerging field of medical imaging analysis, provides a brand-new way to non-invasively quantify high-dimensional imaging features that used to be nearly impossible to be quantified by our bare eyes, and thus could provide valuable evidence to improve decision-making in cancer management [7]. Multiple studies have demonstrated the robustness of developing a radiomic signature as a prognostic tool for various cancers [8, 9], including NPC [10]. However, these early reports either mainly focused on advanced NPC [11], or only predicted local-regional recurrence without metastasis in some studies [12, 13], or just focused on distant metastasis [5]. The same as radiomics, deep learning (DL) has also become one of the most important artificial intelligence (AI) tools, whose application in NPC have been gradually increasing since 2017 [14]. Jing [15] established an end-to-end multi-modality deep survival network (MDSN) to predict the risk of disease progression of NPC patients, but the best performance was a C-index of 0.651. A deep convolution model based on ResNet was established to predict the distant metastasis-free survival of locoregionally advanced NPC patients [16]. In this study, there was no significant difference between deep learning and radiomic signatures when they were used independently. The optimal AUC of the multiple models combined with the clinical features was 0.808. Therefore, it is necessary to establish and evaluate more comprehensive radiomics models based on baseline MRI of NPC in different clinical stages to predict the risk of disease progression of NPC to guide individualized treatment.

Herein, the purpose of this study was to evaluate the value of multiple radiomics models based on multi-parametric MR images at baseline in predicting disease progression probability of non-metastatic NPC after tumor remission. We present the following article in accordance with the STARD reporting checklist.


Patient selection and clinical characteristics

This retrospective study was approved by the institutional review board and the requirement for written informed consent was waived. Between January 2012 and July 2016, a total of 501 patients with NPC confirmed by pathological examinations were identified from the Picture Archiving and Communication System (PACS) of our institution. We only included those patients: 1) who underwent nasopharynx-neck MRI within 2 weeks before any type of antitumor treatment; 2) with no apparent artifacts of any type on MR images that may affect imaging analysis; 3) without evidence suggesting distant metastasis in the baseline assessment; 4) who achieved a complete response after initial treatment; 5) who were then followed up for at least 36 months after complete response. A total of 171 consecutive NPC patients who met the inclusion criteria were included in this study (Fig. 1). During the 36 months after complete response to treatment, NPC patients who developed local-regional recurrence or distant metastasis were assigned to disease progression group, otherwise to non-disease progression group.

Fig. 1
figure 1

Flow-chart for patient inclusion

From the medical records, clinical characteristics were collected including age, gender, histopathology, T-stage, N-stage, overall clinical stage, treatment regimen, pretreatment anti-EBV immunoglobulin A antibody (including early antigen immunoglobulin A [EA-IgA] and anti-viral capsid antigen immunoglobulin A [VCA-IgA]) and regular smoking history. Two senior radiologists (D.L. and Y.Z., with 31 and 17 years of diagnostic radiology experience respectively) performed tumor staging based on nasopharyngoscopy result and medical images according to the TNM classification (7th edition American Joint Committee on Cancer/Union for International Cancer Control) [17], with any disagreements resolved through discussion.

Treatment regimen and follow-up

All patients underwent an individualized treatment based on NCCN guideline [18]. The total radiation doses ranged from 66 to 75 Gy (mean, 70 Gy). All patients received treatment covering the primary tumor, cervical adenopathy, and adjacent at-risk tissues. Patients with stage I tumors received radiotherapy alone, while stage II to IVA patients usually were recommended to have neoadjuvant or adjuvant chemotherapy and/or concurrent chemotherapy with radiotherapy.

Three months after the completion of treatment, treatment response was evaluated for each patient with complete response indicated by MRI and confirmed by nasopharyngoscopy with biopsy. Patients with complete response after the completion of treatment were then followed up at least every 3 months during the first 2 years and every 6 months starting from the third year ever since the initiation of treatment. The routine clinical follow-up work-ups including physical examination, nasopharyngoscopy, one or multiple of imaging modalities, such as CT, MRI, PET/CT and whole-body bone scintigraphy were performed for each patient. Since over 90% of metastases occurred in the first 3 years after the initiation of treatment for NPCs, the cut-off follow-up time was arbitrarily set as 36 months, during which disease progression was defined as the primary clinical endpoint. Disease progression included local-regional recurrence and distant metastasis. Local recurrence for NPC was defined as recurrence at the primary site following complete response after treatment completion, which was confirmed by biopsy and/or MRI of the nasopharynx [12]. Regional recurrences were confirmed using fine needle aspiration biopsy when clinical examination of the neck and MRI indicated progressive cervical adenopathy. Distant metastases were diagnosed based on clinical symptoms, physical examination and imaging findings on nasopharyngeal-neck MRI, thoracic CT, abdominal sonography and whole-body bone scintigraphy, PET-CT, as appropriate.

MRI scan acquisition and tumor segmentation

All of the MR images were acquired using 3.0 T MR scanners (GE Discovery MR 750, General Electric Medical Systems, US) with an 8-channel head and neck phase array coil. For feature extraction, we used axial fast spin echo (FSE) T2-weighted with fat saturation (T2/FS-w) Digital Imaging and Communications in Medicine (DICOM) images and axial fast spoiled gradient-echo (FSPGR) contrast-enhanced T1-weighted (CET1-w) DICOM images that had been archived using PACS. Axial FSPGR CE-T1w imaging was obtained 50–60s after intravenous bolus injection of gadopentetate dimeglumine (Magnevist, Bayer, Leverkusen, Germany) at a dosage of 0.2 ml/kg of body weight at 2.5 ml/second with a power injector. The acquisition parameters were as follows: axial T2-weighted spin-echo images (repetition time [TR]/echo time [TE]: 8472/85 ms, field of view [FOV] = 24 × 24 cm, number of excitations [NEX] = 2.0, slice thickness = 4 mm, spacing between slices = 1.0 mm) and axial contrast-enhanced T1-weighted gradient-echo images (TR/TE: 315/2.9 ms, FOV = 24 × 24 cm, NEX = 1.0, slice thickness = 4 mm, spacing between slices = 1.0 mm).

The VOIs were based on the primary NPC lesions on the T2/FS-w and CET1-w images separately. Tumor segmentation was manually performed by a radiologist (D.B. with 3 years of experience in head and neck imaging) on the radiomics cloud platform V2.1.2 (Huiying Medical Technology Co., Ltd.), which were then reviewed by a senior radiologist (D.L.) who is specialized in head and neck imaging with 31 years of experience. For each lesion, the VOIs were drawn to cover the whole NPC tumor on each consecutive slice, with necrosis, hemorrhage, and cystic areas excluded (Fig. 2; Additional Fig. 1). The interobserver variability in manual segmentation between different readers was assessed. Two radiologists (D.B. and L.H.; with 4 and 6 years of experience, respectively) independently manually segmented the VOIs from the MR images of 30 patients randomly selected from the whole sample at a ratio of one-to-one between the two groups. The Dice similarity coefficient (DSC) for absolute agreement was calculated between segmented VOIs from each patient. According to the guidelines, Dice index of < 0.6, 0.6–0.8, 0.8–1.0, and 1.0 indicates inadequate, good, very good, and ideal consistency, respectively [19]. Then, intraclass correlation coefficient (ICC) for absolute agreement was calculated between radiomics features extracted from the VOIs of 30 randomly selected cases. An ICC above 0.75 was indicative of good agreement [20].

Fig. 2
figure 2

An example of manual segmentation in a 56-year-old male patient with NPC. The segmented tumor is within the red contour in one slice of oblique axial T2WI/FS sequence (A) and the red contour in one slice of axial CE-T1WI sequence (B)

Features extraction and dimensionality reduction

On the radiomics cloud platform V2.1.2 (Huiying Medical Technology Co., Ltd., the Python software, version 3.7.1) a total of 1409 features classified into four different feature classes were extracted, including intensity statistics features, shape- and size-based features, textural features, and filter and wavelet features. To ensure spatial consistency in texture analysis across the images, all images normalized by centering to the mean standard deviation, resampled to voxel size of 1 × 1 × 1 mm3 using B-Spline interpolation, and gray-level discretized by a fixed bin width of 25 in the histogram. Most feature extraction methods conform to the Imaging Biomarker Standardization Initiative (IBSI) standard [21]. With hold-out cross validation scheme (training set: validation set = 7:3, the number of seed = 1234) being used, the least absolute shrinkage and selection operator (LASSO) regression algorithm was used to reduce irrelevant features in the training set.

Feature selection methodology

To further reduce the redundancy of radiomic features in the training set, we performed Pearson correlation analysis with one of the paired features having correlation coefficient of more than 0.9 removed [22]. The remaining radiomic features and all clinical variables form a stable and reliable primary dataset. Then, the multivariate logistic regression analysis was performed for feature selection from the dataset with P < 0.05 indicating independent significant variable. Radiomics and clinical features that were independently significantly associated with disease progression were selected for the final classification modeling on the training set.

Construction and validation of radiomics-based model

Logistic regression analysis and support vector machine (SVM) classifier were applied for model building respectively. On the basis of the feature selection results, we constructed a classification model called a prediction of disease progression MRI-based (PDPM) model (in R with the glm package). An overview of the radiomics analysis process is presented in Fig. 3.

Fig. 3
figure 3

Workflow of the radiomic analysis in the current study

We drafted a nomogram integrating the variables included in the best PDPM model, in R with the rms package, which could visualize the weight of different variables in combined model. The area under the receiver operating characteristic curve (AUC) and its 95% confidence interval (CI) was calculated, which tested the model’s ability to discriminate patients with or without disease progression [23]. Nomogram was applied to reduce the statistical predictive model into a single numerical estimate of the probability of disease progression [24]. Calibration curves, which indicated the calibration ability of the nomogram, were assessed graphically by plotting the actual observed survival rates and the nomogram-predicted survival rates via a bootstrap method with 1000-iteration resampling.

In addition, decision curves analysis (DCA), in R with the rmda package, was introduced to estimate clinical utility of PDPM models based on the threshold probability (probability that triggers a medical intervention by a clinician or patient, equating to the probability at which the harm of a false-positive intervention exceeds the harm of a false-negative non-intervention) [25].

Statistical analysis

The Shapiro-Wilk test was used to determine whether the distribution of all the variables was normal, and Levene’s test was used for identifying the homogeneity of variance. The clinical characteristics of the disease progression group and non-disease progression group were compared using an independent samples t-test or the Mann–Whitney U test, Fisher’s exact test or the chi-squared (χ2) test, as appropriate. Diagnostic performance of the models was assessed by using receiver operating characteristic (ROC) (in R with the pROC package) curve quantified by the AUC, sensitivity, and specificity. The AUCs were compared across models by using the DeLong’s test. The Dice similarity coefficient (DSC) was used to test the consistency between the two radiologists on the manually segmented VOIs of 30 randomly selected samples.

The formula for Dice is as follows:

$$ Dice\left(A,B\right)=2\times \frac{\left|A\cap B\right|}{\left|A\right|+\left|B\right|} $$

A and B represents the segmented VOIs based on the same NPC lesion by two radiologists, respectively. A ∩ B indicates intersection between the manually segmented VOI A and VOI B.

All statistical analyses were performed using SPSS software (version 21.0; IBM Corp., Armonk, NY) and R software (version 4.0.2). P < 0.05 was considered to indicate a statistically significant difference.


Clinical characteristics

A total of 171 patients with pathologically confirmed NPC was included in our study. Baseline clinicopathologic characteristics are summarized in Table 1. The results of correlation analysis showed that age (P = 0.04), gender (P = 0.01) and smoking history (P = 0.02) were significantly associated with disease progression but not for other clinical variables (P > 0.05). During 3 years of follow-up, 76 patients who had developed local-regional recurrence or distant metastasis were assigned to disease progression group with the median follow-up duration of 14.2 months (range, 4–35.5 months), while 95 patients who had no local-regional recurrence or distant metastasis were assigned to non-disease progression group.

Table 1 Baseline clinical characteristics of the patients in the disease progression group and non-recurrence/non-disease progression group

Interobserver variability between readers

Two readers had very good agreement in the manual segmentation with Dice value of 0.85 ± 0.04 (range 0.74–0.91) for T2WI/FS sequences and 0.84 ± 0.05 (range 0.73–0.90) for CE-T1WI sequences. The inter-reader ICC between the two radiologists ranged was 0.93(range 0.83–0.97). These results indicated a favorable inter-observer reproducibility for feature extraction.

Radiomic feature extraction and features selection

A total of 1409 radiomic features were extracted from axial T2/FS-w and CET1-w images on the training set, respectively. After the dimensionality reduction using LASSO regression algorithm, 13 radiomic features were selected. We performed Pearson correlation analysis with the 13 radiomic features to prevent overfitting. None of the 13 radiomic features was found to be significant highly correlated (all r < 0.90) (Fig. 4). A detailed description list of these 13 radiomic features is presented in Table 2. The 13 radiomic features and 10 clinical variables form a stable and reliable primary dataset. No significant differences in clinical and radiomic features of patients were observed between the training and validation cohorts, as shown in Table 3 (all P > 0.05). We included the primary dataset in multivariate logistic regression analysis, HGLZE_T2H (P = 0.03), HGLZE_T1 (P = 0.01), LDLGLE_T1 (P = 0.01), GLNU_T1 (P = 0.03), histology (P = 0.02), T stage (P = 0.03), N stage (P = 0.02), and smoking history (P = 0.003) were found to be significantly associated with disease progression. Based on the result of univariate analysis of clinical characteristics and clinical consensus, we also included age as one of the clinical variables in the final modeling. Finally, we constructed 4 variable sets according to the results of feature selection for subsequent PDPM model establishment, 1) all 13 radiomic features selected and 10 clinical variables (Model 1); 2) 5 clinical variables only (Model 2), determined by multivariate logistic regression analysis and clinical consensus; 3) 4 radiomic features only (Model 3), determined by multivariate logistic regression analysis; 4) 4 radiomic features and 5 clinical variables (Model 4), determined by multivariate logistic regression analysis and clinical consensus.

Fig. 4
figure 4

Pearson correlation coefficient of the 13 significant features

Table 2 List of 13 radiomic features parameters
Table 3 The clinical and radiomic characteristics of patients in the training and validation cohorts

Development and performance of models

Four types of PDPM models were established with logistic regression (Model L1, L2, L3, L4) and support vector machine (SVM) (Model S1, S2, S3, S4) were trained with the four variable sets to predict 3-year disease progression of patients with NPC. Prognostic performances of different PDPM models are shown in Table 4. The best performance was found in the model trained with both 4 radiomic features and 5 clinical features (Model L4), with an AUC of 0.75 in the training cohort (95% CI: 0.66–0.83) and 0.77 in the validation cohort (95% CI: 0.64–0.90). The ROC curve plot of the eight models in validation cohort are shown in Fig. 5. The Model L4 [AUC = 0.77 (95% CI: 0.64–0.90)] exhibited significantly better prediction performance than Model L1 [AUC = 0.66 (95% CI: 0.51–0.81), P = 0.03], Model S2 [AUC = 0.61 (95% CI: 0.45–0.76), P = 0.02] and Model S4 [AUC = 0.67 (95% CI: 0.52–0.82), P = 0.04] in the validation cohorts, respectively.

Table 4 Predictive performance of the PDPM models in predicting the disease progression in the cohorts
Fig. 5
figure 5

The receiver operating characteristic (ROC) curve plots of the eight prediction of disease progression MRI-based (PDPM) models in the training (A, B, C, D) and validation cohort (E, F, G, H)

Model explanation with nomogram and DCA

To provide the clinician with a quantitative method for prediction the risk of disease progression in NPC patients, we constructed a nomogram integrating 4 radiomic features and 5 clinic features included in Model L4 (Fig. 6A). Calibration estimates how close the nomogram estimated risk is to the observed risk, depicted by a calibration plot (Fig. 6B and C). Good calibration was observed for the probability of disease progression in the training cohort and validation cohort. Calibration varies with nomogram calculated probabilities. For instance, the nomogram is more accurate at predicting a disease progression of 50% than 80% (Fig. 6B; note how at a disease progression of 50%, the red line overlaps the black dotted line indicating near perfect calibration however at a recurrence of 80%, the red line and black dotted line do not overlap).

Fig. 6
figure 6

(A): A nomogram incorporated 4 radiomic features and 5 clinical features of Model L4. (B, C): Calibration curves of nomogram developed in the training and validation cohorts

The decision curves analysis (DCA) was used to estimate the clinical utility of Model L4, which integrated with 4 radiomic features and 5 clinical variables, based on the threshold probability. The threshold probability is used to derive the net benefit (defined as the fraction of true-positives subtracted by the fraction of false-positives weighted by the relative harm of a false-positive and false-negative result, Fig. 7). Graphical analysis of the net benefit against the threshold probability yields a decision analysis curve, which can then be used to assess the net benefit of nomogram-assisted decisions at different threshold probabilities, compared to the net benefit of decisions made with the assumption that either all or no patient has the outcome of interest (Fig. 7). When a patient’s threshold probability was within the range of 5 to 75%, the application of Model L4 added more net benefit than the “treat none” or “treat all” strategies, and also added net benefits of “clinical variables only”. Furthermore, we plotted clinical impact curve with variables in Model L4 and only 5 clinical variables separately (Fig. 7B and Fig. 7C). The results revealed that the predictive value of PDPM model was increased by incorporating both clinical variables and radiomic features.

Fig. 7
figure 7

Decision curve analysis for Model L4 for the prediction of disease progression in patients with nasopharyngeal carcinoma (A). Clinical impact curves with variables in Model L4 and only 5 clinical variables in the training and validation cohorts (B, C)


In this study, we extracted radiomic features from pretreatment baseline MRI examinations, and constructed a machine learning model (PDPM model) by combining radiomic features and clinical variables to predict the disease progression of NPC patients after initial treatment. The PDPM Model L4 based on 4 radiomic and 5 clinical features showed good predictive performance, and its prognostic performance was significantly better than other models. We further presented and validated a radiomics nomogram that integrated variables in Model L4 for prediction of disease progression, which achieved satisfactory discrimination performance. The predictive value of PDPM model was increased by incorporating both clinical variables and radiomic features.

Several previous studies have made great efforts to evaluate the relationship between radiomics and clinical outcomes in NPC. Compared with previous studies of NPC on radiomics and disease progression in NPC, which focused only on risk of distant metastasis or local recurrence [5, 13], or was limited to NPC in T4 stage [12], our study included a wide range from early-stage NPC to locally advanced, but not distant metastatic NPC. Compared with the studies that only contoured the largest cross-section of the lesion [12], we obtained 3D VOIs in continuous multiple slices of the entire tumor, which may better contain the heterogeneity of the tumor. There is likely to be a discrepancy between radiomics features those obtained from the largest primary tumor cross-sectional area and the whole tumor. A study included 55 patients with primary colorectal cancer to determine whether there was a difference between contrast enhanced CT texture features from the largest cross-sectional area versus the whole tumor [26]. It demonstrated that whole tumor analysis may provide a more representative evaluation of tumor heterogeneity. The study of Yang also indicated that although 2D and 3D radiomics signature both had favorable prognosis, 3D signature had a better performance [27]. Although a clear superiority of 3D vs. 2D radiomics model performance remains to be scientifically proven, 3D analysis, by virtue of including the entire tumor volume, may better reflect overall tumor heterogeneity. The Dice similarity coefficient (DSC), which measures the degree of agreement between the three-dimensional consistency of tumor segmentation made by two or more readers [28], was calculated for absolute agreement between segmented VOIs in our study. In our study, patients were mainly from non-endemic areas of NPC, while in other studies [5, 29], cohort were mainly composed of homogeneous population from non-endemic regions. Moreover, our study adopted DCA to explain the best predictive model, making our results more favourable for clinical application.

We developed 8 models with two machine learning methods for predicting the 3-year disease progression of NPC. Based on the variable group 4, the Model L4 established by logistic regression algorithm could best discriminate patients who had developed disease progression. Interestingly, a significant difference was found in the prediction performance of Model L4 and Model L1. The 4 radiomic features and 5 clinical variables in Model L4 were selected from 13 radiomic features and 10 clinical variables by multivariate analysis. This result suggested that redundant features exist in a simple combination group [5].

The analysis identified “T2-w_wavelet-HLL_glszm_High Gray Level Zone Emphasis”, “CET1-w_wavelet-HLL_glszm_High Gray Level Zone Emphasis”, “CET1-w_wavelet-LHH_glszm_Gray Level Non-Uniformity”, and “CET1-w_wavelet-LHL_gldm_Large Dependence Low Gray Level Emphasis” as the most important factors of 13 radiomic features in predicting 3-year disease progression. The GLSZM analyzes the distance between groups of voxels with similar grey-levels by counting the number of groups of linked voxels [30]. The GLDM describes pair-wise arrangement of pixels with the set absolute difference in gray-level, in a given direction and distance, and used to highlight local heterogeneity information [31]. GLSZM and GLDM are regional texture features whose ability in differentiating patients with different prognosis has been already observed in various tumors [32, 33]. In line with those studies, the findings here suggest that the heterogeneity of intra-tumoral MRI distribution is an useful biomarker in predicting treatment outcome for patients with NPC. Besides, age, T stage, N stage, histology and smoking history were also found to be clear indicators of risk in the PDPM model, and these findings were consistent with conventional prior knowledge.

Features that were associated with disease progression would make a better and more interpretable model. To provide a clinically applicable method for individual prediction of 3-year disease progression, we further generated a radiomics nomogram that integrated 4 radiomic features and 5 clinical variables in PDPM Model L4. Our results showed that the radiomic nomogram provides better predictive performance compared with clinical nomogram (which only consisted of 5 clinical variables) and full features nomogram (which consisted of 4 radiomic features and 5 clinical variables). This finding was also supported by the results of the ROC analysis. The performance of the nomogram was further verified in the training cohort and validation cohort. The parameters of the nomogram can be easily acquired. Age, T stage, N stage, histology and smoking history were conventional predictive factors and could be obtained at the time of consultation. Moreover, the radiomic features could be obtained via common algorithms. In summary, our study showed that the nomogram can serve as either a scoring system or a visualization tool for disease progression prediction in patients with NPC, thus aiding physicians in rapid evaluation of the metastatic risk in the clinic.

There were some limitations to our study. First, this was a retrospective study with a relatively small sample and patients were included from a single-center. Due to our strict inclusion criteria, such patients with apparent artifacts on MR images were excluded from the study. The possible impact of the exclusion on the results was unknown, but this caused the limitations of the radiomics model established in our study in clinical application. A prospective study with large cohort will be needed to further confirm the conclusions of the current study. Second, it should be noted that we did not investigate the association between the plasma Epstein-Barr virus DNA copy number and disease progression in this study. Although this factor has been confirmed as a potentially useful predictor for the prognosis of NPC, this information was not available for most of our cohort. Third, in our study, we excluded areas of necrosis, hemorrhage, and cystic areas from VOIs, unlike some other studies of NPC. The impact of this practice on model performance and reproducibility might be an interesting topic for future research.


In conclusion, we proposed a noninvasive radiomics model that incorporates the radiomics features and clinical variables to facilitate pretreatment prediction of disease progression in patients with NPC. Our radiomics model may facilitate clinical decision making and potentially improve survival outcomes in selected patients, which required further studies to explore the generalized utility of our model and translate it into clinical practice.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.



area under the curve


logistic regression


nasopharyngeal carcinoma


receiver operating characteristic


support vector machine




  1. Chen YP, Chan ATC, Le QT, Blanchard P, Sun Y, Ma J. Nasopharyngeal carcinoma. Lancet. 2019;394(10192):64–80.

    Article  PubMed  Google Scholar 

  2. Li S, Wang K, Hou Z, Yang J, Ren W, Gao S, et al. Use of radiomics combined with machine learning method in the recurrence patterns after intensity-modulated radiotherapy for nasopharyngeal carcinoma: a preliminary study. Front Oncol. 2018;8:648.

    Article  PubMed  PubMed Central  Google Scholar 

  3. You R, Zou X, Wang SL, Jiang R, Tang LQ, Zhang WD, et al. New surgical staging system for patients with recurrent nasopharyngeal carcinoma based on the AJCC/UICC rTNM classification system. Eur J Cancer. 2015;51(13):1771–9.

    Article  PubMed  Google Scholar 

  4. Zhang L, Huang Y, Hong S, Yang Y, Yu G, Jia J, et al. Gemcitabine plus cisplatin versus fluorouracil plus cisplatin in recurrent or metastatic nasopharyngeal carcinoma: a multicentre, randomised, open-label, phase 3 trial. Lancet. 2016;388(10054):1883–92.

    Article  CAS  PubMed  Google Scholar 

  5. Zhang L, Dong D, Li H, Tian J, Ouyang F, Mo X, et al. Development and validation of a magnetic resonance imaging-based model for the prediction of distant metastasis before initial treatment of nasopharyngeal carcinoma: a retrospective cohort study. EBioMedicine. 2019;40:327–35.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Lambin P, Leijenaar RTH, Deist TM, Peerlings J, de Jong EEC, van Timmeren J, et al. Radiomics: the bridge between medical imaging and personalized medicine. Nat Rev Clin Oncol. 2017;14(12):749–62.

    Article  PubMed  Google Scholar 

  7. Gillies RJ, Kinahan PE, Hricak H. Radiomics: images are more than pictures, they are data. Radiology. 2016;278(2):563–77.

    Article  PubMed  Google Scholar 

  8. Bogowicz M, Riesterer O, Ikenberg K, Stieb S, Moch H, Studer G, et al. Computed tomography radiomics predicts HPV status and local tumor control after definitive radiochemotherapy in head and neck squamous cell carcinoma. Int J Radiat Oncol Biol Phys. 2017;99(4):921–8.

    Article  PubMed  Google Scholar 

  9. Wu S, Zheng J, Li Y, Wu Z, Shi S, Huang M, et al. Development and validation of an MRI-based radiomics signature for the preoperative prediction of lymph node metastasis in bladder cancer. EBioMedicine. 2018;34:76–84.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Raghavan Nair JK, Vallières M, Mascarella MA, el Sabbagh N, Duchatellier CF, Zeitouni A, et al. Magnetic resonance imaging texture analysis predicts recurrence in patients with nasopharyngeal carcinoma. Can Assoc Radiol J. 2019;70(4):394–402.

    Article  PubMed  Google Scholar 

  11. Zhang B, He X, Ouyang F, Gu D, Dong Y, Zhang L, et al. Radiomic machine-learning classifiers for prognostic biomarkers of advanced nasopharyngeal carcinoma. Cancer Lett. 2017;403:21–7.

    Article  CAS  PubMed  Google Scholar 

  12. Zhang LL, Huang MY, Li Y, Liang JH, Gao TS, Deng B, et al. Pretreatment MRI radiomics analysis allows for reliable prediction of local recurrence in non-metastatic T4 nasopharyngeal carcinoma. EBioMedicine. 2019;42:270–80.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Zhang L, Zhou H, Gu D, Tian J, Zhang B, Dong D, et al. Radiomic nomogram: pretreatment evaluation of local recurrence in nasopharyngeal carcinoma based on MR imaging. J Cancer. 2019;10(18):4217–25.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Li S, Deng YQ, Zhu ZL, Hua HL, Tao ZZ. A Comprehensive Review on Radiomics and Deep Learning for Nasopharyngeal Carcinoma Imaging. Diagnostics (Basel). 2021;24(11):1523.

    Article  Google Scholar 

  15. Jing B, Deng Y, Zhang T, Hou D, Li B, Qiang M, et al. Deep learning for risk prediction in patients with nasopharyngeal carcinoma using multi-parametric MRIs. Comput Methods Prog Biomed. 2020;197:105684.

    Article  Google Scholar 

  16. Zhang L, Wu X, Liu J, Zhang B, Mo X, Chen Q, et al. MRI-based deep-learning model for distant metastasis-free survival in Locoregionally advanced nasopharyngeal carcinoma. J Magn Reson Imaging. 2021;53(1):167–78.

    Article  CAS  PubMed  Google Scholar 

  17. Edge SB, Compton CC. The American joint committee on Cancer: the 7th edition of the AJCC cancer staging manual and the future of TNM. Ann Surg Oncol. 2010;17(6):1471–4.

    Article  PubMed  Google Scholar 

  18. Pfister DG, Spencer S, Brizel DM, Burtness B, Busse PM, Caudell JJ, et al. National Comprehensive Cancer Network. Head and neck cancers, version 2.2014. Clinical practice guidelines in oncology. J Natl Compr Cancer Netw. 2014;12(10):1454–87.

    Article  Google Scholar 

  19. Duane F, Aznar MC, Bartlett F, Cutter DJ, Darby SC, Jagsi R, et al. A cardiac contouring atlas for radiotherapy. Radiother Oncol. 2017;122(3):416–22.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Koo TK, Li MY. A guideline of selecting and reporting intraclass correlation coefficients for reliability research. J Chiropr Med. 2016;15:155–63.

    Article  Google Scholar 

  21. Zwanenburg A, Vallières M, Abdalah MA, Aerts HJWL, Andrearczyk V, Apte A, et al. The image biomarker standardization initiative: standardized quantitative radiomics for high-throughput image-based phenotyping. Radiology. 2020;295(2):328–38.

    Article  PubMed  Google Scholar 

  22. Schober P, Boer C, Schwarte LA. Correlation coefficients: appropriate use and interpretation. Anesth Analg. 2018;126(5):1763–8.

    Article  PubMed  Google Scholar 

  23. Steyerberg EW, Vickers AJ, Cook NR, Gerds T, Gonen M, Obuchowski N, et al. Assessing the performance of prediction models: a framework for traditional and novel measures. Epidemiology. 2010;21(1):128–38.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Balachandran VP, Gonen M, Smith JJ, DeMatteo RP. Nomograms in oncology: more than meets the eye. Lancet Oncol. 2015;16(4):e173–80.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Vickers AJ, Elkin EB. Decision curve analysis: a novel method for evaluating prediction models. Med Decis Mak. 2006;26(6):565–74.

    Article  Google Scholar 

  26. Ng F, Kozarski R, Ganeshan B, Goh V. Assessment of tumor heterogeneity by CT texture analysis: can the largest cross-sectional area be used as an alternative to whole tumor analysis. Eur J Radiol. 2013;82(2):342–8.

    Article  PubMed  Google Scholar 

  27. Yang L, Yang J, Zhou X, Huang L, Zhao W, Wang T, et al. Development of a radiomics nomogram based on the 2D and 3D CT features to predict the survival of non-small cell lung cancer patients. Eur Radiol. 2019;29(5):2196–206.

    Article  PubMed  Google Scholar 

  28. Polanski WH, Zolal A, Sitoci-Ficici KH, Hiepe P, Schackert G, Sobottka SB. Comparison of automatic segmentation algorithms for the subthalamic nucleus. Stereotact Funct Neurosurg. 2020;98(4):256–62.

    Article  PubMed  Google Scholar 

  29. Wan Y, Tian L, Zhang G, Xin H, Li H, Dong A, et al. The value of detailed MR imaging report of primary tumor and lymph nodes on prognostic nomograms for nasopharyngeal carcinoma after intensity-modulated radiotherapy. Radiother Oncol. 2019;131:35–44.

    Article  PubMed  Google Scholar 

  30. Bourbonne V, Fournier G, Vallières M, et al. External validation of an MRI-derived radiomics model to predict biochemical recurrence after surgery for high-risk prostate cancer. Cancers (Basel). 2020;12:814.

    Article  CAS  Google Scholar 

  31. Fan TW, Malhi H, Varghese B, et al. Computed tomography-based texture analysis of bladder cancer: differentiating urothelial carcinoma from micropapillary carcinoma. Abdom Radiol. (NY). 2019;44:201–8.

    Article  Google Scholar 

  32. Wang H, Xu X, Zhang X, Liu Y, Ouyang L, du P, et al. Elaboration of a multisequence MRI-based radiomics signature for the preoperative prediction of the muscle-invasive status of bladder cancer: a double-center study. Eur Radiol. 2020;30(9):4816–27.

    Article  PubMed  Google Scholar 

  33. Shiri I, Maleki H, Hajianfar G, Abdollahi H, Ashrafinia S, Hatt M, et al. Next-generation radiogenomics sequencing for prediction of EGFR and KRAS mutation status in NSCLC patients using multimodal imaging and machine learning algorithms. Mol Imaging Biol. 2020;22(4):1132–48.

    Article  CAS  PubMed  Google Scholar 

Download references


Not applicable.


This work was supported by the Non-profit Central Research Institute Fund of Chinese Academy of Medical Sciences (2019XK320073).

Author information

Authors and Affiliations



Dan Bao, Yanfeng Zhao, and Dehong Luo: conception and design; Lin Li, Xinming Zhao, Yanfeng Zhao, and Dehong Luo: administrative support; Dan Bao, Yayuan Geng, Lin Li, and Xinming Zhao: provision of study materials or patients; Dan Bao and Dehong Luo: Collection and assembly of data; Dan Bao, Zhou Liu, Yayuan Geng, Haijun Xu, Ya Zhang, Lei Hu, Yanfeng Zhao, and Dehong Luo: data analysis and interpretation. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Dehong Luo.

Ethics declarations

Ethics approval and consent to participate

The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013). The study was approved by Ethics Committee of National Cancer Center/Cancer Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College and individual consent for this retrospective analysis was waived (21/278–2949).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Additional Figure 1.

An example of manual segmentation with exclusion for necrosis/cystic area (yellow arrow). A 45-year-old male patient with NPC. The segmented tumor is within the red contour in one slice of oblique axial T2WI/FS sequence (a) and the red contour in one slice of axial CE-T1WI sequence (b).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Bao, D., Liu, Z., Geng, Y. et al. Baseline MRI-based radiomics model assisted predicting disease progression in nasopharyngeal carcinoma patients with complete response after treatment. Cancer Imaging 22, 10 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: