Research Paper Volume 13, Issue 14 pp 18645—18657

A nomogram for predicting late radiation-induced xerostomia among locoregionally advanced nasopharyngeal carcinoma in intensity modulated radiation therapy era

Kaixuan Yang1,2, , Wenji Xie1, , Xiangbin Zhang1, , Yu Wang3, , Arthur Shou4, , Qiang Wang1, , Jiangfang Tian1, , Jiangping Yang1, , Guangjun Li1, ,

  • 1 Department of Radiation Oncology, Cancer Center, West China Hospital, Sichuan University, Chengdu 610041, Sichuan, China
  • 2 Department of Radiation Oncology, West China Second University Hospital and Key Laboratory of Obstetrics and Gynecologic and Pediatric Diseases and Birth Defects of Ministry of Education, West China Second University Hospital, Sichuan University, Chengdu 610041, Sichuan, China
  • 3 West China School of Medicine, Sichuan University, Chengdu 610041, Sichuan, China
  • 4 School of Basic Medical Sciences and Forensic Medicine, Sichuan University, Chengdu 610041, Sichuan, China

Received: December 11, 2020       Accepted: June 29, 2021       Published: July 19, 2021      

https://doi.org/10.18632/aging.203308
How to Cite

Copyright: © 2021 Yang et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

Background: Dry mouth sensation cannot be improved completely even though parotids are spared correctly. Our purpose is to develop a nomogram to predict the moderate-to-severe late radiation xerostomia for patients with locoregionally advanced nasopharyngeal carcinoma (LA-NPC) in intensity modulated radiation therapy (IMRT) / volumetric modulated arc radiotherapy (VMAT) era.

Methods: A dataset of 311 patients was retrospectively collected between January 2010 and February 2013. The binary logistic regression was to estimate each factor’s prognostic value for development of moderate-to-severe patient-reported xerostomia at least 2 years (Xer2y) after completion of radiotherapy. Therefore, we can develop a nomogram according to binary logistic regression coefficients. This novel model was validated by bootstrapping analyses.

Results: Contralateral Parotid mean dose (coMD<24.4Gy), VMAT (yes), and platinum-based concurrent chemoradiotherapy (no) were significantly related to patient-reported xerostomia at least 2 years (Xer2y) (all p < 0.001), and were included in the nomogram. Receiver operating characteristic (ROC) analysis revealed AUC (area under the ROC curve) with the value of 0.811 (0.710-0.912) of the nomogram, which was significantly higher than coMD 0.698 (0.560-0.840) from QUANTEC2010 (p<0.001). Calibration plots illustrated that the predicted Xer2y was close to the actual observation, and decision curve analyses (DCA) indicated valid positive net benefits.

Conclusion: We developed a feasible nomogram to predict patient-rated Xer2y based on comprehensive individual data in patients with LA-NPC in the real world. The proposed model is able to facilitate the development of treatment plan and quality of life improvement.

Introduction

During the past decades, more than 80% of nasopharyngeal carcinoma (NPC) patients have prolonged survival after radical chemoradiotherapy. Concerns have arisen with late treatment-related side effects [1], such as xerostomia. It is the sensation of dryness resulting from salivary gland dysfunction or a variation in salivary structure. Decreased and/or thickened saliva is a common feature in xerostomia, which is usually related to oral health, speech, swallowing, and altered taste. Xerostomia with the Radiation Therapy Oncology Group (RTOG) grade 3/4 usually aggravates fatigue, sleeping domains and emotional functioning on quality of life (QoL) scales, which shows the diversified features of xerostomia [2, 3].

The QUANTEC (Quantitative Analyses of Normal Tissue Effects in the Clinic) Group proposed feasible protocols to prevent radiation-induced complications. For instance, severe xerostomia, defined as 25% of baseline in a stimulated salivary flow, could be relieved when at least one of parotids is spared with Dmean<20 Gy or both glands with Dmean<25 Gy [4]. However, the recommendation depends only on the dosimetric factors of the parotid glands. Even though parotids hypofunction contributes radiation-induced xerostomia, it is not the only prognostic factor [5]. Besides, the reduction in salivary function usually occurs in one week after the initiation of radiotherapy and continues thereafter. It takes approximately 2 years to recover after radiotherapy in most cases, which was confirmed by several patient self-reported questionnaires [6]. However, the majority of current studies focused on xerostomia less than 24 months [7, 8]. Actually, it is more reliable and reasonable to focus on 24-months monitor to avoid confounding from the gradual recovery of parotid function. Last but not least, methods used in salivary function measurements were unclear and fluctuated with 20–30% standard deviations for whole mouth evaluation [9, 10].

Recently, the patient-reported outcome has become the crucial step for normal tissue toxicity assessment and treatment tailoring. Compared with patient-reported events, xerostomia symptoms monitored by observer can underestimate the actual dryness symptoms [9, 11] Therefore, patient-reported outcomes should be reasonable research endpoints. It is essential to evaluate the patient-rated xerostomia two years after radiotherapy (Xer2y) based on large volumes of comprehensive clinical data.

Intensity modulated radiation therapy (IMRT), including IMRT with static ports and volumetric modulated arc radiotherapy (VMAT) IMRT, have been the most widely used forms of radiotherapy modalities for NPC. The previous planning studies have confirmed that VMAT improve parotid sparing without compromised target coverage compared to IMRT in terms of dosimetric parameters. However, it is unclear whether the improved parotid sparing can be translated into clinical benefits [12].

Therefore, it is meaningful to explore a general principle for preventing late xerostomia in contemporary technology. In this research, we developed and validated a novel nomogram to predict Xer2y based on a dataset of 311 individual data among patients with NPC in IMRT era.

Materials and Methods

Patients

A dataset of 311 patients with stage III/IVa NPC (AJCC/UICC 8th edition) in our center from January 2010 to February 2013 was used in this research. We predefined the inclusion criteria as follows: 1) histologically detection confirmed by biopsy; 2) treatment with curative IMRT or VMAT, either alone or in combination with platinum-based chemotherapy; 3) without previous radiotherapy, surgery, and/or chemotherapy; 4) without previous malignancies; 5) without moderate-to-severe dry month before treatment, because we focus on radiation-induced xerostomia; 6) no acute (within 3 months after radiotherapy) xerostomia patients censored in 2 years are regarded as no Xer2y as a reason of parotid function recovery in most case [13]. The exclusion criteria were: 1) censored within 3 months after treatment; 2) patients with xerostomia censored within two years; This study has been authorized by the Research Ethics Board of our institution.

Definition of dose-volume histograms (DVHs) and clinical factors

The following factors were used: tumor size, platinum-based dosage of chemotherapy, and DVH factors, et al (Table 1). The contralateral dose was defined as a larger proportion side of the parotid volume outside the PTV. Vd% was described as the parotid volume exposed in d Gy. Dv% was the minimal dose delivered into v% of the parotids volume. Treatment plans were restored, and DVHs parameters were extracted through our in-house script.

Table 1. Univariate and multivariate analyses for Xer2y.

Without xerostomiaWith xerostomiaUnivariate analysisMultivariate analysis
(n = 233)(n = 78)P-valueOR95%CIp value
Age(y)47.4±10.947.6±9.60.8451.1450.834-1.570.403
19-3010 (3.2)4 (1.3)
31-4049 (15.8)13 (4.2)
41-5087 (28)31 (10)
51-6058 (18.6)23 (7.4)
61-7428 (9)7 (2.3)
BMI22 (20.4-24.2)22.2 (20.2-24.5)0.697
Sex0.838
female/male60 (19.3)/173 (55.7)21 (6.8)/57 (18.3)
Family history of cancer0.314
No/Yes194(64.0)/39(12.5)61 (19.6)/17 (5.5)
Smoking0.252
No/Yes122(39.2)/111(35.7)35(11.3)/43(13.8)
Alcohol0.445
No/Yes166(53.4)/67(21.5)52(16.7)/26(8.4)
Histology0.847
Type11 (0.3)0(0.3)
Type2215 (69.1)74 (23.8)
Type312 (3.9)3 (10.0)
Type45 (1.6)1 (0.3)
GTVnx(cm3)53.0 (35.1-78.6)55.5 (34.9-83.7)0.448
GTVnd(cm3)15(8.5-29.6)13.5 (8.4-24)0.548
T category0.701
T2/T3/T47(2.3)/133(42.8)/93(29.9)1(0.3)/46(14.8)/31(10)
N category0.368
N1/N2/N320(6.4)/136(43.7)/77(24.8)4(1.3)/52(16.7)/22(7.1)
8th stage0.657
III/ IVA89(28.6)/144(46.3)32(10.3)/46(14.8)
PGTVnx (cGy)6844(6668.1-6998)6838(6514-7006.7)0.449
PGTVnd (cGy)6815.2(6700.4-6924)6762.6(6644.4-6800.4)0.150
VMAT<0.0010.0310.004-0.236<0.001
No/Yes182(58.5)/51(16.4)76(24.4)/2(0.6)
IGRT0.0161.2790.621-2.6340.505
No/Yes179(57.6)/54(17.4)57(15.8)/21(9.3)
RDD47(45-50)48(45-51)0.1071.0290.991-1.0680.135
Co
MD (cGy)3618(3345.2-3894.6)3831(3440.7-4302.2)<0.0011.0021.001-1.003<0.001
V15 (%)0.90(0.87-0.92)0.90(0.85-0.92)0.626
V20 (%)0.82(0.80-0.86)0.83(0.81-0.85)0.393
V30 (%)0.54(0.48-0.60)0.53(0.47-0.61)0.828
V45 (%)0.36(0.30-0.43)0.37(0.31-0.43)0.683
D503612(3286-3874)3570(3320-3777)0.772
Ip
MD (cGy)3700.4(3400-4085)3631(3299-3906)0.962
V15 (%)0.93(0.90-0.95)0.92(0.88-0.95)0.762
V20 (%)0.83(0.81-0.88)0.85(0.81-0.88)0.335
V30 (%)0.58(0.51-0.65)0.56(0.48-0.67)0.427
V45 (%)0.39(0.31-0.46)0.38(0.32-0.44)0.498
D503780(3522-4181)3831(3440-4387)0.012
Cetuximab0.658
No/Yes199(64)/34(10.9)65(20.9)/13(4.2)
Chemotherapy
IC0.2721.6150.597-4.3720.345
No/Yes7(2.3)/71(22.8)32(10.3)/201(64.6)
CCRT<0.0014.602.20-9.596<0.001
No/Yes13(4.2)/65(20.9)97(31.2)/136(43.7)
AC0.9190.7240.38-1.3780.325
No/Yes42(13.5)/36(11.6)127(40.8)/106(34.1)
IC-CCD140(67-186.3)142.6(66-204.6)0.156
CCRT-CCD75.4(0-151.3)90.6(68.4-179.3)0.011
AC-CCD0(0-114.7)0(0-143.2)0.166
IC*CCRT*AC0.92
RDD, radiotherapy duration days; OR, odds ratio; Co, Contralateral Parotid; 95% CI, 95% confidence interval; IP, Ipsilateral Parotid; MD, mean dose; CCD, cumulative cisplatin dose. Type1, Keratinizing SqCC, Type2, Non-Keratinizing Differentiated, Type3, Non-Keratinizing, Undifferentiated, Type4, Other/Unspecified.

Treatment

Gross tumor volume of nasopharynx (GTVnx) and gross tumor volume of cervical lymph node (GTVnd) were defined as visible tumour and the positive lymph nodes, respectively. Clinical target volume (CTV)-1 contained nasopharynx primary tumour with an additional a 5–10mm margin (2–3mm posteriorly adjacent to the spinal cord or brainstem). CTV-2 contained CTV-1 with the selective neck IB to V area and subclinical sites. GTVnx/nd, CTV-1 and CTV-2 are prescribed to 69.96/73.92, 59.4, and 54 Gy, respectively, with 33-fraction (2.12 or 2.24 Gy per fraction) scheme for 6-7 weeks using 6-MV photons. Over 95% of the prescribed doses are acceptable. IMRT in this research included less than 9 fixed-field beam angles and adopted the step-and-shoot technique, whereas VMAT included a large number of beam directions from the arc trajectory and supports the simultaneous variation in gantry rotation and dose delivery (Supplementary Materials). The regimens of induction chemotherapy contained TPF (docetaxel 60 mg/m2 IV on day 1, cisplatin, 75 mg/m2 IV on day 1 or within 3 days, 5-FU 600 mg/m2 IV on days 1 to 5), which was repeated every 3 weeks for 2–3 cycles. Concurrent chemoradiotherapy consisted of each of 2 regimes: Cisplatin 80 mg/m2 IV every 3 weeks; Cisplatin 30–40 mg/m2 IV weekly. For patients who received adjuvant chemotherapy, PF (cisplatin, 80 mg/m2 IV on day 1, 5-Fu 800 mg/m2/d continuously IV on day 1–5) or TPF regimen was repeated every 3 weeks for 1–4 cycles. Chemotherapy dose can be adjusted according to hematological or non-hematological toxicity.

Endpoints

As was shown in the European Organization for Research and Treatment of Cancer (EORTC) head and neck cancer module (QLQ-H&N35) questionnaire (Supplementary Materials), dry mouth item and sticky saliva item were utilized to quantify patient self-reported xerostomia in a 4-point Likert scale (i.e. none, a bit, quite a bit, to a lot). Each item score represents the degree of xerostomia. We defined moderate-to-severe patient-reported xerostomia at least 2 years (Xer2y) after completion of radiotherapy as the endpoint [7]. It corresponded with “quite a bit” to “a lot” on the 4-point scale.

Statistical analysis

To illustrate the proposed predictive nomogram to models the relationship between a set of predictors and a binary xerostomia response variable, we first conducted an univariate logistic regression analysis to evaluate the Xer2y-prediction ability of each factor and interaction of chemotherapy. For those factors with p < 0.15 in the previous step, we then assessed them in stepwise binary logistic regression which was adjusted by age, Induction Chemotherapy (IC) and Adjuvant Chemotherapy (AC). OR were calculated with the logistic regression mode, Finally, we incorporated significant factors into nomogram according to binary logistic regression using the rms package in R. In particular, the Spearman rank correlation analysis was applied before the stepwise binary logistic regression to decrease the degree of multicollinearity.

The proposed model was validated internally by 1000 bootstrap resamples. We utilized the receiver operating characteristic (ROC) analysis to compare the sensitivity and specificity, and to find the optimal cutoff value. In addition, we used the calibration curve to compare the actual Xer2y against the prediction probability. Decision curve analysis (DCA) illustrated the clinical use of our model by calculating net benefits of the continuous threshold probabilities, which is iterated by putting into the true positives and waving the false positives.

All computations were conducted using SPSS (IBM 22.0) and R software (version 3.5.2). P < 0.05 was recognized as statistically significant.

Results

Patient characteristics

We summarized patient characteristics in Table 2. In general, 78 (25.1%) patients were diagnosed with Xer2y. As shown in Figure 1, The median OS time for the entire cohort was 49 months, and the 5-year OS rate was 80.3% (95% CI 77.4%-83.2%) with a median follow-up of 49 months (ranging from 6 to 74 months).

Table 2. Patient characteristics (n = 311).

CharacteristicsPatients n=311
median (range)/
No. of patients (%)
Age(years), Mean (SD)47.5 (36.9-58.1)
BMI, Median (IQR)22 (20.4-24.5)
Gender
Female/male81 (26)/230 (74)
Family history of cancer
No/yes255 (82.0)/56 (18.8)
Cigarette smoking
No/Yes157 (50.5)/154 (49.5)
Alcohol
No/Yes218 (70.1)/93 (30.0)
Histology
Keratinizing SqCC1 (0.3)
Non-Keratinizing, Differentiated289 (92.9)
Non-Keratinizing, Undifferentiated15 (4.8)
Other/Unspecified6 (1.9)
GTVnx(cm3), Median (IQR)55.2 (35.1-78)
GTVnd(cm3), Median (IQR)15 (8.5-28.5)
T category
T2/T3/T48 (2.6)/179 (57.6)/124 (39.9)
N category
N1/N2/N324 (7.7)/188 (60.5)/99 (31.8)
8th UICC/AJCC stage
III/ IVA121 (38.9)/190 (61.1)
PGTVnx (cGy), Median (IQR)6842.4 (6631.7-6999)
PGTVnd (cGy), Median (IQR)6807 (6694-6900)
IMRT/VMAT258 (83)/53 (17)
IGRT
No/Yes236 (75.9)/75 (24.1)
Radiotherapy duration days48 (45-50)
Cetuximab
No/Yes264 (84.9)/47 (15.1)
Chemoradiotherapy
IC/CCRT/AC272 (87.5)/201 (64.6)/142 (45.7)
BMI, body mass index; IQR, interquartile range; SD, standard deviation; IC, Induction Chemotherapy; CCRT, Concurrent Chemoradiation Radiotherapy; AC, Adjuvant Chemotherapy; IGRT, image-guided radiotherapy; PGTV, planning gross target volume; VMAT, volumetric modulated arc radiotherapy.
Kaplan–Meier overall survival curves for all 311 patients.

Figure 1. Kaplan–Meier overall survival curves for all 311 patients.

Univariate and multivariate logistic regression analyses

In univariate logistic regression analysis, VMAT, image-guided radiotherapy (IGRT), radiotherapy duration days (RDD), coMD, ipsilateral Parotid D50 (ipD50), platinum-based concurrent chemoradiotherapy (CCRT), and CCRT cumulative cisplatin dose (CCD) were statistically associated with Xer2y in Table 1 (p < 0.15). According to Spearman’s analyses, the significant collinearity between CCRT and CCRT-CCD (r = 0.781, p < 0.001). coMD and ipD50 (r = 0.96, p < 0.001) were revealed in Table 3. We included CCRT and coMD instead of CCRT-CCD and ipD50 in the binary logistic regression by the use of a published selection criteria [14]. Age, IC and AC were also considered as adjusted factors. Therefore, age, IC, AC, VMAT, IGRT, RDD, coMD and CCRT were incorporated into the binary logistic regression: coMD (OR: 1.002, 95%CI: 1.001-1.003, p <0.001), VMAT (OR:0.031, 95%CI: 0.004-0.236, p <0.001), and CCRT (OR:4.60, 95%CI: 2.20-9.596, p <0.001) were predictors of Xer2y (Table 1). Since the majority of eligible patients did not receive cetuximab, for 264 patients who had not received the target therapy, 65 (25%) experienced Xer2y. After a subgroup binary logistic regression in these specific patients, as was noted previously, coMD (OR: 1.002, 95%CI: 1.001-1.003, p <0.001), VMAT (OR: 0.031, 95%CI: 0.004-0.242, p =0.001), and CCRT (OR: 5.389,95%CI: 2.375-12.23, p <0.001) were still significantly related with Xer2y.

Table 3. Critical values (p < 0.15 in univariate test) for Spearman’s rank correlation analyses.

VMATIGRTRDDcoMDipD50CCRTCC-CCD
VMAT10.0520.0740.046-0.0390.0490.101
IGRT0.05210.0020.0690.054-0.012-0.055
RDD0.0740.0021-0.049-0.075-0.018-0.003
coMD0.0460.069-0.04910.96***0.0940.044
ipD50-0.0390.054-0.0750.96***10.0840.042
CCRT0.049-0.012-0.0180.0940.08410.781***
CC-CCD0.101-0.055-0.0030.0440.0420.781***1
*, P < 0.05; **, P < 0.01; ***, P < 0.001.

Nomogram development and validation

According to coefficients of the multivariate logistic regression analyses, a nomogram (namely QQmodel) was developed visually to predict Xer2y (Figure 2). The QQmodel was internally validated using the bootstrapping analyses. As the ROC curves shown in Figure 3A, the QQmodel gained an AUC of 0.811 (95%CI: 0.71–0.912), which was significantly superior (p<0.001) to coMD (0.698, 95%CI: 0.560–0.840). The cut-off value (dose–volume constraints) of coMD was 24.4 Gy. The proposed nomogram also showed promising prediction efficiency (sensitivity: 56%, specificity: 96%). Calibration plots of QQmodel performed well as the predicted Xer2y and the actual observation gained an agreeable consistency (Figure 3B). DCA illustrated that QQmodel achieved much more positive net benefits than coMD for the majority of the threshold probabilities, which demonstrated an encouraging clinical effect of the QQmodel (Figure 3C).

Nomogram of Xer2y occurrence prediction. For each individual patient, the value of three variables (VMAT, CCRT and coMD) are translated into points by projecting them into the upper-most line (point scale), respectively. Summing the points of the three variables and projecting the total points value downward to the bottom-most line can determine the probability of this patient to have Xer2y (patient-reported xerostomia at least 2 years after radiotherapy). VMAT (volumetric modulated arc radiotherapy). CCRT (platinum-based concurrent chemoradiotherapy). coMD (Contralateral Parotid mean dose).

Figure 2. Nomogram of Xer2y occurrence prediction. For each individual patient, the value of three variables (VMAT, CCRT and coMD) are translated into points by projecting them into the upper-most line (point scale), respectively. Summing the points of the three variables and projecting the total points value downward to the bottom-most line can determine the probability of this patient to have Xer2y (patient-reported xerostomia at least 2 years after radiotherapy). VMAT (volumetric modulated arc radiotherapy). CCRT (platinum-based concurrent chemoradiotherapy). coMD (Contralateral Parotid mean dose).

(A) ROC curves comparing the sensitivity and specificity of the coMD and nomogram (B) Calibration plots of the nomogram for Xer2y incidence prediction. X-axis indicates the predicted probabilities with Xer2y while y-axis shows the actual events. The ideal prediction will correlate when slope equals to 1 (the black broken line in the figure). (C) Decision curves of two risk models for Xer2y incident prediction. The horizontal axis represents the risk threshold while the vertical-axis denotes standardized net benefit. The solid black line is the net benefit when no patients have Xer2y while the dash gray line suggests the net benefits where patients have Xer2y at a certain risk threshold. The red and blue curves imply the results of the Xer2y on the basis of coMD and nomogram, respectively. AUCROC (area under the receiver operating characteristic curve). 95% CI (95% confidence interval).

Figure 3. (A) ROC curves comparing the sensitivity and specificity of the coMD and nomogram (B) Calibration plots of the nomogram for Xer2y incidence prediction. X-axis indicates the predicted probabilities with Xer2y while y-axis shows the actual events. The ideal prediction will correlate when slope equals to 1 (the black broken line in the figure). (C) Decision curves of two risk models for Xer2y incident prediction. The horizontal axis represents the risk threshold while the vertical-axis denotes standardized net benefit. The solid black line is the net benefit when no patients have Xer2y while the dash gray line suggests the net benefits where patients have Xer2y at a certain risk threshold. The red and blue curves imply the results of the Xer2y on the basis of coMD and nomogram, respectively. AUCROC (area under the receiver operating characteristic curve). 95% CI (95% confidence interval).

Discussion

It is widely acknowledged that dysfunction of the parotids is one of the most problematic late effects in head and neck cancer patients [5]. We established and validated the first nomogram to predict the occurrence of moderate-to-severe Xer2y in IMRT/VMAT era. Results of the present study indicated that coMD, VMAT and CCRT are independent prognosticators of moderate-to-severe Xer2y from comprehensive clinical and dosimetric factors among patients with LANPC. We proved our model is superior to coMD from QUANTEC2010 [4] criteria.

The predictive ability of coMD (<24.4 Gy) for Xer2y was confirmed in our model. After QUANTEC2010 published, Recent research found that the contralateral mean parotid gland dose and baseline xerostomia are the top two important predictors to patients-reported QoL for dry mouth and 51.6% cases studies suffered from moderate-to-severe Xer6m in that study [15, 16]. Given that some xerostomia will steadily restore in around 2 years, 25.1% of Xer2y in the present research will be biologically plausible. Parotid glands mean dose is a significant risk factor for patient-reported moderate-to-severe xerostomia for post treatment within 6, 12, 18 and 24 months [7]. Furthermore, xerostomia was closely related to parotid glands mean dose in the actual delivered dose research [5, 17]. These evidences suggest that coMD could play a crucial role in radiation xerostomia.

In the present study, the dosimetric parameters of parotid gland V15, V20, V30, V45 and D50 were not significant predictors of Xer2y. There are several existing explanations. Firstly, V15, V30, V45 were assessed by observer-based grading [18]. Relationship between patient-reported events and observer-based grading have been verified to be inconsistent [9, 11]. Secondly, V20 is from the delivery dose rather than planning dose [17].

Compared with IMRT, we found that VMAT was better at preventing Xer2y. This finding supported the previous work [19]. Furthermore, we revealed that VMAT conserves salivary function when coMD is established. A possible explanation is that other salivary glands, such as the oral cavity accessory glands also play a key role in the sensation of chronic dry mouth apart from the parotid gland [20, 21]. Since IMRT is characterized by typically arranging less than 9 fixed-field beam angles, VMAT instead has numerous beam directions from the arc trajectory. Thanks to there being no beam modulation by the MLCs, VMAT achieves less high-dose areas, more homogeneous dose distribution, and better spare to the salivary glands while treatment beams are sequenced one after another in IMRT treatment plan [22]. Some researchers found that the isodose distribution in IMRT approaches to VMAT as increasing numbers of beams are utilized [23]. Moreover, VMAT has more beam entry angles and supports the simultaneous variation in MLC leaf positions, gantry rotation speed and dose rate, resulting in great reductions in monitor units and treatment delivery time per fraction by transforming IMRT into VMAT for head and neck cancer [24]. Lower MU has its advantages such as the potential decrease in total body dose due to leakage and scattering from MLCs [25, 26]. Extra low-dose areas may have potential protective effects on other oral cavity accessory glands [12]. Reducing treatment time increases patient comfort and decreases intra-fractional movement [27]. However, other studies mentioned that VMAT is not superior to IMRT in terms of dry month [28]. This inconsistency probably results from planning objectives, dose calculation algorithm and patient characteristics. Therefore, conditions should be declared prior to the clinical use.

Our research implied that CCRT, which was more powerful than CCRT-CCD, was an independent prognosticator for moderate-to-severe Xer2y. This conclusion confirmed several previous works [29]. Nonetheless, other published studies failed to reveal increased xerostomia in CCRT [30], because their datasets were moderate-to-severe xerostomia patients within 2 years after radiotherapy and gradual recovery of parotid glands function might attribute to unmeasured confounding. Recently, alternative concurrent chemotherapy regimens with mild toxicity for patients with LA-NPC have been reported [31]. According to our knowledge, there is no published data from randomized control trial to address the role of CCRT with IMRT versus IMRT alone for LA-NPC.

Validation is required to highlight the contribution of the proposed nomogram. External validation is usually a gold standard method. However, because of the modest sample size, even researchers in some eminent institutes preferred carrying out an internal validation. As one of the internal validation, the bootstrapping analysis, such as the calibration curve, is advantageous especially for a relatively small sample size [32]. However, conditions should be declared before the clinical use of model.

This work used the patient self-reported moderate-to-severe xerostomia at least 2 years (Xer2y) after completion of radiotherapy quantified by the EORTC QLQ-H&N35 questionnaire as the endpoints. Quality of life (QoL) reflects the awareness of the consequences of the disease and the burden that the disease imposes on the individual's daily function. This is performed by patients’ own sense. EORTC QLQ-H&N35 was one of the most widely used tools to reveal the shift of QoL [6]. As one of significant problems in QoL, Xerostomia (common feature of decreased and/or thickened saliva) has been successfully quantified by the Xerostomia item in EORTC QLQ-H&N35 [7]. Compared with xerostomia symptoms assessed by physicians, patient self-reported events were proved to be more reasonable [9]. However, other items associated with xerostomia-related QoL, such as global health status, swallowing, social eating and social contact, should be further collected and investigated in the future.

The AUC of coMD was 0.698(0.560-0.840), which was similar to former dosimetric study [15]. Such common predictive ability elucidates that salivary flow preservation due to accurately spared parotids fails to improve dry mouth sensation or of quality of life effectively [33]. In this research, the AUC of integration nomogram of coMD, VMAT and CCRT was increased to 0.811(0.710-0.912). Although parotids and submandibular glands are proposed to be spared in QUANTEC2010 and the Radiation Therapy Oncology Group 1016 trial in 2011, only 90% of bilateral parotids were outlined and submandibular glands were not delineated in 56% out of treatment plans in public hospitals [34], not to say other salivary glands. Therefore, our results have great significance for clinical work.

Some limitations of the present study are addressed here. Firstly, the research did not take position error and anatomical variation into consideration. Given the retrospective nature of this study design. our Cone-beam Computed Tomography data was insufficient for analysis. Further study should be conducted to assess their role in the nomogram. Secondly, the endpoint in this research was late xerostomia 2 years after radiotherapy. Although it helped to avoid interference with parotid function recovery, it might result in selective bias from unmeasured and unknown confounders as a reason of the inclusion / exclusion criteria. Thirdly, our participant samples were in a modest size and were taken from an endemic area in China. Fourthly, it is difficult to completely avoid effort and alertness of bias in planning optimization, although clinical plans with the same goals [35] are validated by a senior Medical Dosimetrist in our department. Finally, the reality of the nomogram would be more persuasive if an external validation is executed.

In conclusion, we established a clinically feasible nomogram involving coMD (<24.4 Gy), VMAT (yes), and CCRT (no) that could quantify the risk for patient-rated Xer2y among LA-NPC based on large volumes of comprehensive clinical data in IMRT/VMAT era. This novel predictive model reduced intervention from parotid function recovery, and will be valuable for the development of treatment plan and the improvement of quality of life.

Author Contributions

YK and XW carried out data analyses and drafted the manuscript. ZX and WQ collected the treatment planning data, TJ and YJ did the follow up; AS and YK performed the statistical analysis. LG, YK. WY designed, coordinated, and supervised the study and critically reviewed and discussed the manuscript. All authors have read and approved the final version of the manuscript.

Acknowledgments

The authors would like to thank Dr. Xu F from the Cancer Center at our institution for his fruitful discussions for the study. QQmodel is dedicated to Ms. Yunqi Zhong.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Funding

This study is funded by the National Natural Science Foundation of China (Grant No. 81472807).

References

  • 1. Brodin NP, Kabarriti R, Garg MK, Guha C, Tomé WA. Systematic Review of Normal Tissue Complication Models Relevant to Standard Fractionation Radiation Therapy of the Head and Neck Region Published After the QUANTEC Reports. Int J Radiat Oncol Biol Phys. 2018; 100:391–407. https://doi.org/10.1016/j.ijrobp.2017.09.041 [PubMed]
  • 2. Oba MK, Innocentini LM, Viani G, Ricz HM, de Carvalho Reis T, Ferrari TC, de Macedo LD. Evaluation of the correlation between side effects to oral mucosa, salivary glands, and general health status with quality of life during intensity-modulated radiotherapy for head and neck cancer. Support Care Cancer. 2021; 29:127–34. https://doi.org/10.1007/s00520-020-05454-5 [PubMed]
  • 3. Hawkins PG, Lee JY, Mao Y, Li P, Green M, Worden FP, Swiecicki PL, Mierzwa ML, Spector ME, Schipper MJ, Eisbruch A. Sparing all salivary glands with IMRT for head and neck cancer: Longitudinal study of patient-reported xerostomia and head-and-neck quality of life. Radiother Oncol. 2018; 126:68–74. https://doi.org/10.1016/j.radonc.2017.08.002 [PubMed]
  • 4. Deasy JO, Moiseenko V, Marks L, Chao KS, Nam J, Eisbruch A. Radiotherapy dose-volume effects on salivary gland function. Int J Radiat Oncol Biol Phys. 2010; 76:S58–63. https://doi.org/10.1016/j.ijrobp.2009.06.090 [PubMed]
  • 5. Teng F, Fan W, Luo Y, Xu S, Gong H, Ge R, Zhang X, Wang X, Ma L. A Risk Prediction Model by LASSO for Radiation-Induced Xerostomia in Patients With Nasopharyngeal Carcinoma Treated With Comprehensive Salivary Gland-Sparing Helical Tomotherapy Technique. Front Oncol. 2021; 11:633556. https://doi.org/10.3389/fonc.2021.633556 [PubMed]
  • 6. Bian X, Song T, Wu S. Outcomes of xerostomia-related quality of life for nasopharyngeal carcinoma treated by IMRT: based on the EORTC QLQ-C30 and H&N35 questionnaires. Expert Rev Anticancer Ther. 2015; 15:109–19. https://doi.org/10.1586/14737140.2015.961427 [PubMed]
  • 7. Beetz I, Steenbakkers RJ, Chouvalova O, Leemans CR, Doornaert P, van der Laan BF, Christianen ME, Vissink A, Bijl HP, van Luijk P, Langendijk JA. The QUANTEC criteria for parotid gland dose and their efficacy to prevent moderate to severe patient-rated xerostomia. Acta Oncol. 2014; 53:597–604. https://doi.org/10.3109/0284186X.2013.831186 [PubMed]
  • 8. Kaae JK, Stenfeldt L, Hyrup B, Brink C, Eriksen JG. A randomized phase III trial for alleviating radiation-induced xerostomia with chewing gum. Radiother Oncol. 2020; 142:72–78. https://doi.org/10.1016/j.radonc.2019.09.013 [PubMed]
  • 9. Meirovitz A, Murdoch-Kinch CA, Schipper M, Pan C, Eisbruch A. Grading xerostomia by physicians or by patients after intensity-modulated radiotherapy of head-and-neck cancer. Int J Radiat Oncol Biol Phys. 2006; 66:445–53. https://doi.org/10.1016/j.ijrobp.2006.05.002 [PubMed]
  • 10. Ni X, Tian T, Chen D, Liu L, Li X, Li F, Liang F, Zhao L. Acupuncture for Radiation-Induced Xerostomia in Cancer Patients: A Systematic Review and Meta-Analysis. Integr Cancer Ther. 2020; 19:1534735420980825. https://doi.org/10.1177/1534735420980825 [PubMed]
  • 11. van der Laan HP, Van den Bosch L, Schuit E, Steenbakkers RJ, van der Schaaf A, Langendijk JA. Impact of radiation-induced toxicities on quality of life of patients treated for head and neck cancer. Radiother Oncol. 2021; 160:47–53. https://doi.org/10.1016/j.radonc.2021.04.011 [PubMed]
  • 12. Vanetti E, Clivio A, Nicolini G, Fogliata A, Ghosh-Laskar S, Agarwal JP, Upreti RR, Budrukkar A, Murthy V, Deshpande DD, Shrivastava SK, Dinshaw KA, Cozzi L. Volumetric modulated arc radiotherapy for carcinomas of the oro-pharynx, hypo-pharynx and larynx: a treatment planning comparison with fixed field IMRT. Radiother Oncol. 2009; 92:111–17. https://doi.org/10.1016/j.radonc.2008.12.008 [PubMed]
  • 13. Blanco AI, Chao KS, El Naqa I, Franklin GE, Zakarian K, Vicic M, Deasy JO. Dose-volume modeling of salivary function in patients with head-and-neck cancer receiving radiotherapy. Int J Radiat Oncol Biol Phys. 2005; 62:1055–69. https://doi.org/10.1016/j.ijrobp.2004.12.076 [PubMed]
  • 14. Tang X, Li Y, Tian X, Zhou X, Wang Y, Huang M, Ren L, Zhou L, Xue J, Ding Z, Zhu J, Xu Y, Peng F, et al. Predicting severe acute radiation pneumonitis in patients with non-small cell lung cancer receiving postoperative radiotherapy: Development and internal validation of a nomogram based on the clinical and dose-volume histogram parameters. Radiother Oncol. 2019; 132:197–203. https://doi.org/10.1016/j.radonc.2018.10.016 [PubMed]
  • 15. Beetz I, Schilstra C, van der Schaaf A, van den Heuvel ER, Doornaert P, van Luijk P, Vissink A, van der Laan BF, Leemans CR, Bijl HP, Christianen ME, Steenbakkers RJ, Langendijk JA. NTCP models for patient-rated xerostomia and sticky saliva after treatment with intensity modulated radiotherapy for head and neck cancer: the role of dosimetric and clinical factors. Radiother Oncol. 2012; 105:101–06. https://doi.org/10.1016/j.radonc.2012.03.004 [PubMed]
  • 16. Wilkie JR, Mierzwa ML, Casper KA, Mayo CS, Schipper MJ, Eisbruch A, Worden FP, El Naqa I, Viglianti BL, Rosen BS. Predicting late radiation-induced xerostomia with parotid gland PET biomarkers and dose metrics. Radiother Oncol. 2020; 148:30–37. https://doi.org/10.1016/j.radonc.2020.03.037 [PubMed]
  • 17. Lou J, Huang P, Ma C, Zheng Y, Chen J, Liang Y, Li H, Yin Y, Liu D, Yu G, Li D. Parotid gland radiation dose-xerostomia relationships based on actual delivered dose for nasopharyngeal carcinoma. J Appl Clin Med Phys. 2018; 19:251–60. https://doi.org/10.1002/acm2.12327 [PubMed]
  • 18. Pan XB, Liu Y, Li L, Qu S, Chen L, Liang SX, Chen KH, Liang ZG, Zhu XD. Prognostic nomogram of xerostomia for patients with nasopharyngeal carcinoma after intensity-modulated radiotherapy. Aging (Albany NY). 2020; 12:1857–66. https://doi.org/10.18632/aging.102717 [PubMed]
  • 19. Cui T, Ward MC, Joshi NP, Woody NM, Murray EJ, Potter J, Dorfmeyer AA, Greskovich JF Jr, Koyfman SA, Xia P. Correlation between plan quality improvements and reduced acute dysphagia and xerostomia in the definitive treatment of oropharyngeal squamous cell carcinoma. Head Neck. 2019; 41:1096–103. https://doi.org/10.1002/hed.25594 [PubMed]
  • 20. Cao J, Zhang X, Jiang B, Chen J, Wang X, Wang L, Sahoo N, Zhu XR, Ye R, Blanchard P, Garden AS, Fuller CD, Gunn GB, Frank SJ. Intensity-modulated proton therapy for oropharyngeal cancer reduces rates of late xerostomia. Radiother Oncol. 2021; 160:32–39. https://doi.org/10.1016/j.radonc.2021.03.036 [PubMed]
  • 21. Teng F, Fan W, Luo Y, Ju Z, Gong H, Ge R, Tong F, Zhang X, Ma L. Reducing Xerostomia by Comprehensive Protection of Salivary Glands in Intensity-Modulated Radiation Therapy with Helical Tomotherapy Technique for Head-and-Neck Cancer Patients: A Prospective Observational Study. Biomed Res Int. 2019; 2019:2401743. https://doi.org/10.1155/2019/2401743 [PubMed]
  • 22. Pigorsch SU, Kampfer S, Oechsner M, Mayinger MC, Mozes P, Devecka M, Kessel KK, Combs SE, Wilkens JJ. Report on planning comparison of VMAT, IMRT and helical tomotherapy for the ESCALOX-trial pre-study. Radiat Oncol. 2020; 15:253. https://doi.org/10.1186/s13014-020-01693-2 [PubMed]
  • 23. Quan EM, Li X, Li Y, Wang X, Kudchadker RJ, Johnson JL, Kuban DA, Lee AK, Zhang X. A comprehensive comparison of IMRT and VMAT plan quality for prostate cancer treatment. Int J Radiat Oncol Biol Phys. 2012; 83:1169–78. https://doi.org/10.1016/j.ijrobp.2011.09.015 [PubMed]
  • 24. Fan J, Xing L, Ma M, Hu W, Yang Y. Verification of the machine delivery parameters of a treatment plan via deep learning. Phys Med Biol. 2020; 65:195007. https://doi.org/10.1088/1361-6560/aba165 [PubMed]
  • 25. Hauri P, Schneider U. Whole-body dose equivalent including neutrons is similar for 6 MV and 15 MV IMRT, VMAT, and 3D conformal radiotherapy. J Appl Clin Med Phys. 2019; 20:56–70. https://doi.org/10.1002/acm2.12543 [PubMed]
  • 26. Hall EJ. Intensity-modulated radiation therapy, protons, and the risk of second cancers. Int J Radiat Oncol Biol Phys. 2006; 65:1–7. https://doi.org/10.1016/j.ijrobp.2006.01.027 [PubMed]
  • 27. Mangesius J, Seppi T, Weigel R, Arnold CR, Vasiljevic D, Goebel G, Lukas P, Ganswindt U, Nevinny-Stickel M. Intrafractional 6D head movement increases with time of mask fixation during stereotactic intracranial RT-sessions. Radiat Oncol. 2019; 14:231. https://doi.org/10.1186/s13014-019-1425-7 [PubMed]
  • 28. Oliver M, McConnell D, Romani M, McAllister A, Pearce A, Andronowski A, Wang X, Leszczynski K. Evaluation of the trade-offs encountered in planning and treating locally advanced head and neck cancer: intensity-modulated radiation therapy vs dual-arc volumetric-modulated arc therapy. Br J Radiol. 2012; 85:1539–45. https://doi.org/10.1259/bjr/26344684 [PubMed]
  • 29. Hey J, Setz J, Gerlach R, Vordermark D, Gernhardt CR, Kuhnt T. Effect of Cisplatin on parotid gland function in concomitant radiochemotherapy. Int J Radiat Oncol Biol Phys. 2009; 75:1475–80. https://doi.org/10.1016/j.ijrobp.2008.12.071 [PubMed]
  • 30. Li XY, Chen QY, Sun XS, Liu SL, Yan JJ, Guo SS, Liu LT, Xie HJ, Tang QN, Liang YJ, Wen YF, Guo L, Mo HY, et al. Ten-year outcomes of survival and toxicity for a phase III randomised trial of concurrent chemoradiotherapy versus radiotherapy alone in stage II nasopharyngeal carcinoma. Eur J Cancer. 2019; 110:24–31. https://doi.org/10.1016/j.ejca.2018.10.020 [PubMed]
  • 31. Tang LQ, Chen DP, Guo L, Mo HY, Huang Y, Guo SS, Qi B, Tang QN, Wang P, Li XY, Li JB, Liu Q, Gao YH, et al. Concurrent chemoradiotherapy with nedaplatin versus cisplatin in stage II-IVB nasopharyngeal carcinoma: an open-label, non-inferiority, randomised phase 3 trial. Lancet Oncol. 2018; 19:461–73. https://doi.org/10.1016/S1470-2045(18)30104-9 [PubMed]
  • 32. Steyerberg EW, Harrell FE Jr. Prediction models need appropriate internal, internal-external, and external validation. J Clin Epidemiol. 2016; 69:245–47. https://doi.org/10.1016/j.jclinepi.2015.04.005 [PubMed]
  • 33. Eisbruch A. Reducing xerostomia by IMRT: what may, and may not, be achieved. J Clin Oncol. 2017; 25:4863–64. https://doi.org/10.1200/JCO.2007.13.4874 [PubMed]
  • 34. Vallard A, Guy JB, Mengue Ndong S, Vial N, Rivoirard R, Auberdiac P, Méry B, Langrand-Escure J, Espenel S, Moncharmont C, Ben Mrad M, Diao P, Goyet D, Magné N. Intensity-modulated radiotherapy or volumetric-modulated arc therapy in patients with head and neck cancer: Focus on salivary glands dosimetry. Head Neck. 2016; 38:1028–34. https://doi.org/10.1002/hed.24398 [PubMed]
  • 35. Castriconi R, Cattaneo GM, Mangili P, Esposito P, Broggi S, Cozzarini C, Deantoni C, Fodor A, Di Muzio NG, Vecchio AD, Fiorino C. Clinical Implementation of Knowledge-Based Automatic Plan Optimization for Helical Tomotherapy. Pract Radiat Oncol. 2021; 11:e236–44. https://doi.org/10.1016/j.prro.2020.09.012 [PubMed]