Proposal of a Pretreatment Nomogram for Predicting Local Recurrence after Intensity-Modulated Radiation Therapy in T4 Nasopharyngeal Carcinoma: A Retrospective Review of 415 Chinese Patients.

PURPOSE: Local relapse-free survival (LRFS) differs widely among patients with T4 category nasopharyngeal carcinoma (NPC). We aimed to build a nomogram incorporating clinicopathological information to predict LRFS in T4 NPC after definitive intensity-modulated radiation therapy (IMRT).
MATERIALS AND METHODS: Retrospective study of 415 Chinese patients with non-metastatic T4 NPC treated with definitive IMRT with or without chemotherapy at our cancer center between October 2009 and September 2013. The nomogram for LRFS at 3 and 5 years was generated based on multivariate Cox proportional hazards regression, and validated using bootstrap resampling, assessing discriminative performance using the concordance index (C-index) and determining calibration ability via calibration curves.
RESULTS: Five-year LRFS was 88.8%. We identified and incorporated four independent prognostic factors for LRFS: ethmoid sinus invasion, primary gross tumor volume, age, and pretreatment body mass index. The C-index of the nomogram for local recurrence was 0.732 (95% confidence interval, 0.726 to 0.738), indicating excellent predictive accuracy. The calibration curve revealed excellent agreement between nomogram-predicted and observed LRFS probabilities. Risk subgroups based on total point score cutoff values enabled effective discrimination of LRFS.
CONCLUSION: This pretreatment nomogram enables clinicians to accurately predict LRFS in T4 NPC after definitive IMRT, and could help to facilitate personalized patient counselling and treatment strategies.

Introduction

Nasopharyngeal carcinoma (NPC), a malignant tumor that originates in the nasopharyngeal epithelium, is endemic in southern China and Southeast Asia where the age-standardized incidence ranges from 20 to 50 per 100,000 person-years [1]. Radiotherapy is the primary treatment modality for nondisseminated NPC as a result of its complex anatomical location and high radiosensitivity [2,3]. The introduction of intensity-modulated radiation therapy (IMRT) has greatly improved locoregional control in NPC [4]. Recently, Wu et al. [5] reported the 10-year local relapse-free survival (LRFS) rates for patients with T1, T2, and T3 NPC receiving IMRT were satisfactory (94.2%, 92.5%, and 91.4%, respectively; p > 0.05); however, the 10-year LRFS rate for T4 NPC remains poor (79.3%, p < 0.05 vs. T1-T3 NPC). Recurrent NPC remains a challenging clinical problem and traditional management does little to limit local control and, more seriously, frequently induces severe, life-threatening late complications [6]. Thus, it is necessary to develop accurate prognostic models to identify patients with T4 NPC at high risk of local recurrence as early as possible; this may provide an opportunity establish more intensive treatment and follow-up strategies for patients with T4 NPC at high risk of local recurrence.

The American Joint Committee on Cancer (AJCC)/Union for International Cancer Control (UICC) TNM staging system is the most widely used prognostic tool for predicting the risk of recurrence and survival outcomes. Patients with non-metastatic NPC are stratified according to tumor size, invasion and lymph node involvement. However, significant variations in local recurrence after definitive IMRT are observed within each individual T category [5]. This indicates additional, more relevant variables should be integrated to improve prognostication of clinical outcomes for patients with T4 NPC. Previous research demonstrated a number of independent prognostic factors such as such as age, World Health Organization (WHO) histologic type, primary gross tumor volume (GTVp), and pretreatment lactate dehydrogenase (LDH), can significantly predict local recurrence [7-9]. Therefore, a comprehensive tool that estimates individual risk of recurrence by incorporating TNM stage and clinicopathological information, including host factors, treatment strategies, pathological factors, and anatomic factors, could represent a valuable decision-making tool for clinicians.

Nomograms incorporate a variety of important prognostic factors and have been accepted as useful predictive statistical models for quantifying individual risk in several cancers [10,11]. Nomograms have been demonstrated to have a better ability to predict prognosis than TNM stage in many types of cancer [10-12]. Hence, the present study aimed to establish an effective prognostic pretreatment nomogram to enable individualized prediction of the risk of local recurrence in patients with T4 NPC treated with IMRT.

Materials and Methods

1. Patient selection and staging evaluation

A total of 415 patients with NPC treated at Sun Yat-sen University Cancer Center (Guangzhou, China) between October 2009 and September 2013 were enrolled in this retrospective study. The eligibility criteria were as follows: (1) newly diagnosed, pathologically proven, previously untreated NPC; (2) no evidence of distant metastasis; (3) stage T4N0-3M0 NPC; (4) treated with IMRT and completed planned radiotherapy with good radiation treatment compliance; and (5) complete clinicopathologic and treatment data available.

Pretreatment staging workup included a complete patient history, physical and neurological examinations, hematology and biochemistry profiles, fiberoptic nasopharyngoscopy, magnetic resonance imaging (MRI) of the nasopharynx and neck, chest radiography, abdominal ultrasonography, and single positron emission tomography computed tomography (PET-CT). Each patient was restaged according to the eighth edition of the AJCC/UICC staging system. Two senior physicians experienced in head and neck cancer (L.C. and W.F.L.) separately evaluated all scans and all disagreements were resolved by consensus.

2. T classification criteria and stage grouping according to the eighth edition of the AJCC/UICC staging system and proposed staging system

According to the eighth edition of the AJCC/UICC staging system, tumors confined to the nasopharynx or with extension to the oropharynx and/or nasal cavity without parapharyngeal involvement are classified as T1 disease; tumors with extension to the parapharyngeal space and/or adjacent soft tissue involvement (medial peterygoid, lateral pterygoid, and prevertebral muscle) are classified as T2 disease; tumors with infiltration of bony structures at the skull base, cervical vertebra, pterygoid structures and/or paranasal sinus are classified as T3 disease; tumors with intracranial extension, involvement of the cranial nerves, hypopharynx or orbit or extensive soft tissue involvement (beyond the lateral surface of the lateral pterygoid muscle, parotid gland) are classified as T4 disease.

As T4 lesions could have coincidental lesions that should belong to T1-3 categories, we investigated the association between such coincidental lesions and local recurrence in patients with T4 disease.

3. Treatment strategies

During the study period, institutional guidelines recommended concurrent chemoradiotherapy with or without neoadjuvant/adjuvant chemotherapy for stage IVA-IVB. All patients were treated with fractionated IMRT 5 days/wk. Target volumes were contoured according to the institutional treatment protocol [13], in agreement with International Commission on Radiation Units and Measurements Reports 50 and 62 [14]. The prescribed radiation doses were: total of 66-72 Gy at 2.12-2.43 Gy/fraction to the planning target volume (PTV) of the GTVp, 64-70 Gy/28-33 fractions to the involved nodal gross tumor volume PTV, 60-63 Gy/28-33 fractions to the PTV of CTV-1 (high-risk clinical target volume), and 54-56 Gy/28-33 fractions to the PTV of CTV-2 (low-risk clinical target volumes). A total of 192/415 patients (46.3%) completed radiotherapy on schedule and 223/415 patients (53.7%) had prolonged radiotherapy schedules.

Concurrent chemotherapy was 30-40 mg/m2 cisplatin administered weekly or 80-100 mg/m2 cisplatin administered on weeks 1, 4, and 7 of radiotherapy, beginning first day of IMRT. Neoadjuvant chemotherapy was two or three cycles of cisplatin (75 mg/m2) with docetaxel (75 mg/m2); cisplatin (80 mg/m2) with 5-fluorouracil (1,000 mg/m2); or cisplatin (60 mg/m2) with 5-fluorouracil (600 mg/m2) and docetaxel (60 mg/m2) every 3 weeks.

4. Follow-up and end-points

After completing IMRT, patients returned for follow-up appointments at least every 3 months in the first three years and every 6 months thereafter (or until death). Follow-up was measured from first day of therapy to day of last examination or death. Patients who achieved a complete response but developed recurrence at the primary site after completing IMRT were defined as suffering local recurrence. Conventional examinations during follow-up were similar to the pretreatment assessment. Local relapse was diagnosed via MRI of the nasopharynx, biopsy, or both; PET-CT imaging was used when necessary. The primary end-point was LRFS, calculated from first day of therapy to date of local failure. Patients whose attendance at recent follow-up examinations was not recorded in their medical records were followed up by telephone.

5. Statistical analysis

This study aimed to establish an effective prognostic nomogram for individualized prediction of local recurrence and survival. SPSS ver. 22.0 (IBM Corp., Armonk, NY) was used to identify independent prognostic factors for LRFS. Forty-three candidate predictors were included in univariable analysis: age, sex, WHO pathologic type, N category, family history of NPC, smoking, pretreatment hemoglobin, body mass index (BMI), pretreatment viral capsid antigen IgA, pretreatment early antigen IgA, pretreatment Epstein-Barr virus DNA, pretreatment LDH; cumulative dose to the GTVp; chemotherapy; pretreatment white blood cell (WBC) count; pretreatment neutrophil (Neu) count; pretreatment Neu/WBC ratio; pretreatment C-reactive protein; GTVp; prolonged radiation treatment; nasal cavity invasion; oropharyngeal invasion; parapharyngeal space invasion; medial pterygoid invasion; lateral pterygoid invasion; prevertebral muscle invasion; cervical vertebrae invasion; maxillary sinus invasion; frontal sinus invasion; sphenoid sinus invasion; ethmoidal sinus invasion; pterygopalatine fossa invasion; foramen lacerum invasion; foramen ovale invasion; clivus invasion; petrous apex invasion; great wing of sphenoid bone invasion; base of sphenoid bone invasion; cranial nerve invasion; intracranial invasion; orbit invasion; hypopharynx invasion and extensive soft tissue involvement (beyond the lateral surface of the lateral pterygoid muscle, parotid gland).

LRFS curves were plotted using the Kaplan-Meier method. Univariate and multivariate analysis were performed via Cox proportional hazard regression. Variables achieving significance (p < 0.05) were included in the Cox regression model with backward elimination for multivariate analyses.

A nomogram based on the results of the multivariate analysis was formulated to provide visualized risk prediction using the rms package of R ver. 3.2.5 (https://www.r-project.org/). A final model was selected using a backward stepdown process that incorporated Akaike’s information criterion [15]. Nomograms were generated for the probabilities of local recurrence at 36 and 60 months.

The predictive accuracy of the nomogram was assessed by comparing discrimination and calibration. Predictive accuracy was assessed over 1,000 bootstrap replicates [10,11]. The discriminatory performance of the nomograms was measured using the concordance index (C-index) by comparing nomogram-predicted versus observed Kaplan-Meier estimates of survival probability, a variable equivalent to the area under curve of receiver operating characteristic curves for censored data. Maximum C-index (1.0) indicates a perfect match between predicted and actual probabilities, while 0.5 indicates the random chance of the model correctly predicting outcomes. Calibration of the nomograms for 3- and 5-year LRFS was performed using a calibration curve incorporating model-predicted and actual observed probabilities. Biases of predictive performance were first estimated using bootstrap sampling, then corrected with apparent measures to produce bias-corrected values to avoid overfitting [12]. Patients were classified into three risk subgroups for local recurrence based on the scores calculated by the nomograms using X-tile software [16]. Two-tailed p-values < 0.05 were considered significant.

6. Ethical statement

This study was approved by our institutional review board (approve number: YB2017-032). As this was a retrospective analysis of routine clinical data, a waiver of the requirement for individual informed consent was granted by the institutional ethics committee.

Results

1. Baseline characteristics, patterns of treatment failure, and survival

Median follow-up for the whole cohort was 59.8 months (range, 3.3 to 84.4 months); median age was 47 years (range, 17 to 78 years). All 415 patients had T4N0-3M0 (stage IV) disease. Detailed clinicopathological characteristics and treatment factors are summarized in Table 1. A total of 45/415 patients (10.8%) developed local recurrence by last follow-up. The latency from the end of treatment to first recurrence ranged from 3.4 months to 85.6 months. Three- and 5-year LRFS were 91.6% and 88.8%, respectively.

Table 1.

Clinical characteristics and treatment parameters of 415 patients with T4 nasopharyngeal carcinoma and association with local recurrence-free survival in univariate analysis

Characteristic	No. of patients (%) (n=415)	p-value[a)]
Age (yr)
≤ 47	216 (52.0)	0.011
> 47	199 (48.0)
Sex
Male	328 (79)	0.659
Female	87 (21)
WHO pathology
Type I	5 (1.2)	0.168
Type II	15 (3.6)
Type III	395 (95.2)
N category[b)]
N0	46 (11.1)	0.589
N1	262 (63.1)
N2	70 (16.9)
N3	37 (8.9)
Family history of NPC
Yes	38 (9.2)	0.628
No	377 (90.8)
Smoking
Yes	180 (43.4)	0.985
No	235 (56.6)
Hb (g/L)
< 113	22 (5.3)	0.402
113-151	289 (69.6)
> 151	104 (25.1)
BMI (kg/m2)
< 22.2	204 (49.2)	0.003
≥ 22.2	211 (50.8)
VCA-IgA
< 1:80	65 (15.7)	0.354
1:80-1:320	255 (61.4)
≥ 1:640	95 (22.9)
EA-IgA
< 1:10	103 (24.8)	0.181
1:10-1:20	139 (33.5)
≥ 1: 40	173 (41.7)
EBV-DNA (copy/mL)
< 11,600	206 (49.6)	0.245
≥ 11,600	209 (50.4)
LDH (U/L)
< 245	370 (89.2)	0.564
≥ 245	39 (9.4)
Cumulative dose to GTVp (Gy)
< 68	195 (47.0)	0.430
≥ 68	219 (52.8)
Chemotherapy
None	7 (1.7)	0.066
Concurrent	104 (25.1)
Neoadjuvant±concurrent	304 (73.3)
WBC (×109/L)
< 4	14 (3.4)	0.135
4-0	345 (83.1)
≥ 10	56 (13.5)
Neu (×109/L)
< 2	17 (4.1)	0.048
2-7	334 (80.5)
≥ 7	64 (15.4)
Neu/WBC
< 0.7	289 (69.6)	0.819
≥ 0.7	126 (30.4)
CRP (g/mL)
< 1	92 (22.2)	0.336
1-3	127 (30.6)
≥ 3	196 (47.2)
GTVp (cm3)
< 65.7	207 (49.9)	0.002
≥ 65.7	208 (50.1)
Prolonged radiation treatment, total days
No	192 (46.3)	0.644
Yes	223 (53.7)

WHO, World Health Organization; NPC, nasopharyngeal carcinoma; Hb, hemoglobin; BMI, body mass index; VCA-IgA, viral capsid antigen IgA; EA-IgA, pretreatment early antigen IgA; EBV, Epstein-Barr virus; LDH, pretreatment lactate dehydrogenase; GTVp, primary gross tumor volume; WBC, white blood cells; Neu, neutrophils; CRP, C-reaction protein.

p-values were calculated using univariate Cox regression analysis,

According to the 8th edition of the American Joint Committee on Cancer/Union for International Cancer Control staging system.

2. Independent prognostic factors for LRFS

Eleven factors including age, BMI, Neu, GTVp, orbit invasion, foramen lacerum invasion, maxillary sinus invasion, intracranial invasion, ethmoidal sinus invasion, petrous apex invasion, and invasion of the great wing of the sphenoid bone were significantly associated with local recurrence in univariate analysis (Tables 1 and 2).

Table 2.

Lesions in T1-T4 category structures among 415 patients with T4 nasopharyngeal carcinoma and association with local recurrence-free survival in univariate analysis

Structure	No. of patients (%) (n=415)	p-value[a)]
Coincidental lesions in T1 category structures
Nasal cavity invasion
No	188 (45.3)	0.397
Unilateral	145 (34.9)
Bilateral	82 (19.8)
Oropharyngeal invasion
No	343 (82.7)	0.296
Unilateral	60 (14.5)
Bilateral	12 (2.9)
Coincidental lesions in T2 category structures
Parapharyngeal space invasion
No	50 (12.0)	0.638
Unilateral	252 (60.8)
Bilateral	113 (27.2)
Medial pterygoid invasion
No	152 (36.6)	0.385
Unilateral	235 (56.6)
Bilateral	28 (6.7)
Lateral pterygoid invasion
No	286 (68.9)	0.499
Unilateral	124 (29.9)
Bilateral	5 (1.2)
Prevertebral muscle
No	83 (20.0)	0.122
Unilateral	189 (45.5)
Bilateral	143 (34.5)
Coincidental lesions in T3 categories structures
Cervival vertebrae invasion
No	400 (96.4)	0.410
Unilateral	13 (3.1)
Bilateral	2 (0.5)
Maxillary sinus invasion
No	351 (84.6)	0.007
Unilateral	54 (13.0)
Bilateral	10 (2.4)
Frontal sinus invasion
No	414 (99.8)	0.825
Unilateral	1 (0.2)
Bilateral	0
Sphenoid sinus invasion
No	196 (47.2)	0.076
Unilateral	218 (52.5)
Bilateral	1 (0.2)
Ethmoidal sinus invasion
No	339 (81.7)	< 0.001
Unilateral	43 (10.4)
Bilateral	33 (8.0)
Pterygopalatine fossa invasion
No	194 (46.7)	0.066
Unilateral	180 (43.4)
Bilateral	41 (9.9)
Foramen lacerum invasion
No	73 (17.6)	0.013
Unilateral	264 (63.6)
Bilateral	78 (18.8)
Foramen ovale invasion
No	169 (40.7)	0.069
Unilateral	226 (54.5)
Bilateral	20 (4.8)
Clivus invasion
No	36 (8.7)	0.246
Unilateral	378 (91.1)
Bilateral	1 (0.2)
Petrous apex invasion
No	43 (10.4)	0.004
Unilateral	255 (61.4)
Bilateral	117 (28.2)
Great wing of sphenoid bone invasion
No	175 (42.2)	0.041
Unilateral	209 (50.4)
Bilateral	31 (7.5)
Base of sphenoid bone invasion
No	28 (6.7)	0.095
Unilateral	384 (92.6)
Bilateral	3 (0.7)
Coincidental lesions in T4 categories structures
Cranial nerve invasion
Yes	137 (33.0)	0.457
No	278 (67.0)
Intracranial invasion
No	74 (17.8)	0.032
Unilateral	303 (73.0)
Bilateral	38 (9.2)
Orbit invasion
No	356 (85.8)	0.019
Unilateral	54 (13.0)
Bilateral	5 (1.2)
Hypopharynx invsaion
No	393 (94.7)	0.159
Unilateral	16 (3.9)
Bilateral	6 (1.4)
Extensive soft tissue invasion (beyond lateral surface of lateral pterygoid muscle, parotid gland)
No	346 (83.4)	0.649
Unilateral	67 (16.1)
Bilateral	2 (0.5)

p-values were calculated using univariate Cox regression analysis.

The significant 11 factors from univariate analysis were included in Cox multivariate analysis of LRFS to adjust for various prognostic factors. Ethmoidal sinus invasion, GTVp, age, and pretreatment BMI were independent prognostic factors for LRFS in multivariate analysis. Bilateral ethmoidal sinus invasion, large GTVp, old age, and low BMI were significantly associated with poorer LRFS (Table 3).

Table 3.

Summary of multivariate analysis of independent prognostic factors for LRFS in T4 nasopharyngeal carcinoma

Characteristic	LRFS
	p-value[a)]	HR	95% CI
Ethmoidal sinus invasion
No	Reference
Unilateral	0.087	1.999	0.905-4.415
Bilateral	0.010	2.798	1.285-6.090
GTVp (cm3)	0.015	2.277	1.170-4.430
< 65.7
≥ 65.7
Age (yr)	0.010	2.243	1.215-4.142
≤ 47
> 47
BMI (kg/m2)	0.020	0.459	0.239-0.883
< 22.2
≥ 22.2

LRFS, local relapse-free survival; HR, hazard ratio; CI, confidence interval; GTVp, primary gross tumor volume; BMI, body mass index.

p-values were calculated using the Cox proportional hazards model.

3. Development and internal validation of a nomogram for local recurrence

The prognostic nomogram that incorporates all significant independent factors for LRFS is presented in Fig. 1. The values for each of the model covariates were mapped to points on a scale axis ranging from 0 to 100. By adding up the total score for all variables and locating this value on the total point scale, the probability of local recurrence can be determined by drawing a vertical line to the total score. The nomogram indicates ethmoid sinus involvement has the most significant effect on LRFS, followed by GTVp, age and BMI. In Fig. 2, the x-axes denote the predicted probabilities of LRFS from the nomograms, the y-axes denote actual probabilities of LRFS calculated using the Kaplan-Meier method, and the 45-degree line represents the ideal reference line where nomogram-predicted and observed LRFS probabilities are exactly the same. The calibration plots indicated excellent agreement between nomogram-predicted LRFS and actual 3- and 5-year LRFS probabilities. The c-index of the nomogram for local recurrence was 0.732 (95% confidence interval [CI], 0.726 to 0.738), demonstrating the nomogram represents a reasonable model for predicting local recurrence.

Fig. 1.

Nomogram for local recurrence in non-metastatic T4 nasopharyngeal carcinoma after radical intensity-modulated radiation therapy. BMI, body mass index; GTVp, primary gross tumor volume.

Fig. 2.

Calibration curves for predicting local recurrence at 3 years (A) and 5 years (B). Nomogram-predicted probability of local recurrence is plotted on the x-axis; actual observed local recurrence-free probability is plotted on the y-axis.

4. Stratification of the risk of local recurrence for patients with T4 NPC

Based on the cutoff values of the total point scores for local recurrence determined using X-tile software, patients were classified into a low-risk group (score < 123), intermediate-risk group (123-202), and high-risk group (203-325). In summary, a total of 193 (46.5%), 152 (36.5%), and 70 (16.9%) patients were assigned to the low-, intermediate-, and highrisk groups, respectively. Three-year LRFS for the low-, intermediate-, and high-risk groups was 97.9%, 90.1%, and 76.2%, respectively; and 5-year LRFS was 96.7%, 86.2% and 70.8%, respectively. Kaplan-Meier survival analysis demonstrated this risk stratification effectively discriminated the risk of LRFS in patients with T4 NPC treated with IMRT (p < 0.001) (Fig. 3).

Fig. 3.

Kaplan-Meier local recurrence-free survival curves for patients with T4 nasopharyngeal carcinoma based on risk group stratification. Patients were classified into a low-risk group (score < 123), intermediate-risk group (123- 202), and high-risk group (203-325) using total nomogram point score cutoff values determined using X-tile software.

Discussion

T category is not the sole risk factor for local recurrence in NPC: the LRFS rates of patients with the same T category can vary widely [5]. Thus, T category is not capable of accurately predicting local recurrence in patients with NPC treated with IMRT. Furthermore, while previous research demonstrated no significant differences in LRFS between patients with T1, T2, and T3 NPC; patients with T4 disease have a significantly poorer local control rate [5]. Establishing a prognostic nomogram for each T category may help to provide more individualized and accurate estimation of the risk of local recurrence following IMRT; this would be especially significant for patients with T4 NPC who have the highest risk of local recurrence. More importantly, analysis of the risk of local relapse in T4 NPC is crucial for formulation of appropriate treatment approaches; a nomogram may enable clinicians to develop more individualized multidisciplinary treatment and follow-up strategies.

To the best of our knowledge, this is the first study to develop a prognostic nomogram for LRFS based on widelyassessed clinicopathological factors in patients with nonmetastatic T4 NPC treated with definitive IMRT. The multi-variable Cox analyses demonstrated ethmoidal sinus invasion, GTVp, age and BMI were independent prognostic factors for LRFS in non-metastatic T4 NPC. The resulting nonograms demonstrated excellent discriminative ability (0.732; 95% CI, 0.726 to 0.738) with excellent agreement between nomogram-predicted LRFS and actual LRFS probabilities, as confirmed by the calibration curve. In addition, the nomogram-generated scores enabled patients with T4 NPC to be further stratified into three distinct risk groups for LRFS, and as expected, the proposed risk groups significantly discriminated the risk of LRFS within T4 NPC, especially for the high-risk group. Thus, this nomogram may represent a clinically useful tool to supplement the traditional TNM staging system.

The nomogram revealed ethmoid sinus involvement had the most significant association with LRFS, followed by GTVp, age, and BMI. Consistent with numerous previous studies, GTV-p was a significant prognostic factor for LRFS in T4 NPC after IMRT [6,7,17,18]. A large GTVp was associated with poorer LRFS; the following factors may explain this phenomenon. First, a large GTVp, which reflects a high tumor burden, often indicates tumor hypoxia. Tumor hypoxia plays a key role in the development of radio- and chemoresistance, and promotes disease progression such as local recurrence [19-22]. Secondly, a large GTVp indicates the primary tumor is closer to critical normal organs at risk, which negatively affects local radiation dose escalation to the involved areas.

However, Sakata et al. [23] reported that deep primary tumor infiltration was a more important adverse prognostic factor than a large GTVp. we found both ethmoid sinus invasion and a large GTVp were adverse independent prognostic factors for local recurrence in T4 NPC, with ethmoid sinus invasion having a greater contribution than GTVp to local recurrence. The ethmoid sinus is located distant to the nasopharynx and is a low-risk site of invasion in NPC. After lesions in various anatomical structures belonging to the T1-T4 categories were incorporated into the analysis, we found that only ethmoid sinus invasion was an adverse independent prognostic factor for local recurrence in T4 NPC. In fact, deep infiltration of the ethmoid sinus is a marker of severe extension; therefore, ethmoid sinus invasion is inherently associated a larger GTVp.

BMI and age were also significant prognostic factors for local recurrence in T4 NPC; older age (> 47-years-old) and lower BMI (< 22.2 kg/m2) were associated with poor LRFS. Consistent with these results, Cheng et al. [9] showed older age (> 40-years-old) adversely affected locoregional control in NPC after definitive conventional radiotherapy. Moreover, Tang et al. [24] demonstrated low pretreatment BMI was associated with local and systemic recurrence. It is plausible that lower BMI and older age may be associated with poorer immune system function and reduced tolerance to intensive treatment, which may ultimately result poorer survival.

The nomogram established in the present study has some limitations. First, this was a retrospective analysis of medical records; therefore, selection bias may exist. However, the relatively large sample size should reduce the influence of such bias. Secondly, the nomogram was not subject to external verification due to the lack of an external validation cohort. However, we conducted internal validation via 1,000 bootstrap replicates to avoid over-fitting of the data. In summary, the nomogram developed in this study may enable clinicians to make accurate predictions of LRFS for patients with T4 NPC after definitive IMRT. We hope the nomograms could be applied to aid patient counseling or guide individual treatment, which may help to further improve local control. The prognostic value offered by this pretreatment nomogram for predicting local recurrence could enable more aggressive treatment protocols to be delivered to patients at high-risk of local recurrence. Concurrent chemoradiotherapy (CCRT) is now recommended as a standard treatment for T4N0-3M0 NPC. However, induction chemotherapy has been shown to reduce the risk of locoregional recurrence [25]. Thus, induction chemotherapy combined with CCRT may be able to achieve an improvement in local control for patients with T4 NPC, although further evidence is needed to prove this hypothesis. Additional studies are required to further validate this nomogram in other patient cohorts.