- Research article
- Open Access
- Open Peer Review
Normal tissue complication probability model parameter estimation for xerostomia in head and neck cancer patients based on scintigraphy and quality of life assessments
BMC Cancervolume 12, Article number: 567 (2012)
With advances in modern radiotherapy (RT), many patients with head and neck (HN) cancer can be effectively cured. However, xerostomia is a common complication in patients after RT for HN cancer. The purpose of this study was to use the Lyman–Kutcher–Burman (LKB) model to derive parameters for the normal tissue complication probability (NTCP) for xerostomia based on scintigraphy assessments and quality of life (QoL) questionnaires. We performed validation tests of the Quantitative Analysis of Normal Tissue Effects in the Clinic (QUANTEC) guidelines against prospectively collected QoL and salivary scintigraphic data.
Thirty-one patients with HN cancer were enrolled. Salivary excretion factors (SEFs) measured by scintigraphy and QoL data from self-reported questionnaires were used for NTCP modeling to describe the incidence of grade 3+ xerostomia. The NTCP parameters estimated from the QoL and SEF datasets were compared. Model performance was assessed using Pearson’s chi-squared test, Nagelkerke’s R2, the area under the receiver operating characteristic curve, and the Hosmer–Lemeshow test. The negative predictive value (NPV) was checked for the rate of correctly predicting the lack of incidence. Pearson’s chi-squared test was used to test the goodness of fit and association.
Using the LKB NTCP model and assuming n=1, the dose for uniform irradiation of the whole or partial volume of the parotid gland that results in 50% probability of a complication (TD50) and the slope of the dose–response curve (m) were determined from the QoL and SEF datasets, respectively. The NTCP-fitted parameters for local disease were TD50=43.6 Gy and m=0.18 with the SEF data, and TD50=44.1 Gy and m=0.11 with the QoL data. The rate of grade 3+ xerostomia for treatment plans meeting the QUANTEC guidelines was specifically predicted, with a NPV of 100%, using either the QoL or SEF dataset.
Our study shows the agreement between the NTCP parameter modeling based on SEF and QoL data, which gave a NPV of 100% with each dataset, and the QUANTEC guidelines, thus validating the cut-off values of 20 and 25 Gy. Based on these results, we believe that the QUANTEC 25/20-Gy spared-gland mean-dose guidelines are clinically useful for avoiding xerostomia in the HN cohort.
Head and neck (HN) cancer is a leading cause of cancer mortality in Taiwan, and radiotherapy (RT) plays an important role in its treatment. Xerostomia is a common complication after RT for HN [1–5]. Severe xerostomia is defined as long-term salivary dysfunction compared with the pre-RT function, based on the Late Effects of Normal Tissues–Subjective, Objective, Management, Analytic (LENT–SOMA) criteria [6–8].
Whole-mouth salivary function has been shown to be related to quality of life (QoL) [9, 10] and has been used to compare different treatment strategies in clinical oncology trials. Kakoei et al.  have shown that the decrease in saliva and xerostomia resulting from RT can negatively affect QoL for patients who undergo RT. Several prospective studies conducted over the past decade have reported the RT dose constraints to allow preservation of parotid gland function based on salivary flow measurements or salivary gland scintigraphy [10, 12–14].
In the present prospective study, we longitudinally observed parotid gland function using salivary scintigraphy to measure the salivary excretion factor (SEF) in patients receiving intensity-modulated radiotherapy (IMRT). Moreover, a self-reported QoL questionnaire (QLQ-C30) and a xerostomia-specific questionnaire (QLQ-H&N35) were completed by patients before RT and periodically after therapy to assess the interrelationships with salivary function. The normal tissue complication probability (NTCP) model proposed by Lyman has been used to determine the dose for uniform irradiation of the whole or partial volume of the parotid gland that results in 50% probability of a complication (TD50) in patients with local disease [15, 16]. NTCP-fitted parameters for patients with local disease were investigated using both datasets.
Between August 2007 and June 2008, 65 HN cancer patients who had undergone primary or postoperative RT for various malignancies were initially included in the study. Patients who suffered from Sјögren’s syndrome or any other medical cause of xerostomia were excluded. The use of any medication known to affect salivary gland function was prohibited. After eliminating patients because of missed appointments, refusal, and organizational problems, it was possible to objectively evaluate parotid gland function using scintigraphy and QoL questionnaires after RT initiation in 31 patients. The present prospective study enrolled these 31 HN cancer patients who received primary (n=15) or post-operative RT (n=16) with IMRT at Chiayi Chang Gung Memorial Hospital of the Chang Gung Medical Foundation. Nineteen patients received concurrent chemotherapy: 18 received five to seven courses of weekly cisplatin (40 mg/m2 CDDP), and one received two courses of a PF regimen (80 mg/m2 CDDP on day 1 + 800 mg/m2 5-FU on days 1–5, every 21 days). Five of these patients received additional adjuvant chemotherapy with a PF regimen for two to three courses (n=4) or a TEF regimen (60 mg/m2 taxol on day 1 + 20 mg/m2 CDDP on day 1 + 800 mg/m2 5-FU on days 1–2) for one course (n=1).
Patients with successful salivary flow scintigraphy imaging and full completion of QoL questionnaires before and during 1 year after treatment were included. No data were missing for these 31 patients. This study was approved by the institutional review board of the hospital (IRB-95-1430B).
Patients were immobilized from head to shoulders using a commercially available thermoplastic mask and/or an individually customized bite block. Computed tomography (CT) images (2.5-mm slice thickness) acquired from the top of the vertex to the level of the carina, containing 512 × 512 pixels in each slice, were examined. Both parotid glands were delineated by a radiation oncologist. We used the Pinnacle treatment planning system to perform inverse planning and dose optimization. For each patient, IMRT plans with five or seven coplanar portals were created. Dose distributions were calculated, and separate dose-volume histograms (DVHs) were generated for each parotid gland, enabling each gland to be analyzed separately. IMRT treatment mode was used in a sequential method .
IMRT was delivered by a computer-controlled and auto-sequencing segmented or dynamic multileaf collimator of a linear accelerator (Varian Clinac 21 EX or Elekta Precise), with the aim of sparing the parotid glands (predominantly contralateral side) while treating the primary targets and lymph nodes at risk. The prescribed doses were 67.4 to 70.8 Gy (mean dose, 69.8 Gy) to the macroscopic tumor planning target volume (PTV1), 54.8 to 70.8 Gy (mean dose, 62.0 Gy) to the resected tumor bed planning target volume (PTV2), and 46.8 Gy to the subclinical disease planning target volume (PTV3), delivered at 1.8 to 2 Gy per fraction.
Based on the Radiation Therapy Oncology Group studies 0615, and 0225 , the planning objectives for PTVs were a minimum dose to >95% of the target, with no more than 5% of any PTV1 receiving ≥110% of the prescribed dose. The structural constraints used were a parotid gland mean dose of ≤26 Gy or V30Gy ≤50%; for the oral cavity excluding the PTV, the mean dose must be ≤40 Gy. The mean DVH values for the parotid gland were calculated for each patient. All data are based on mean DVHs obtained from Pinnacle3® using a bin-size resolution of 0.01 Gy. The dose calculation resolution was 2.5 mm for all IMRT plans.
Salivary gland scintigraphy
All patients received salivary scintigraphy. Stimulated whole-mouth saliva was collected before RT and at various time intervals; for this analysis, the 1-year follow-up time point was used. Scintigraphy was performed after 4 h of fasting. After the patient received an intravenous injection of 10 mCi of 99mTc pertechnetate, sequential images of the left and right anterior views of the head and neck were acquired at 1 min/frame for 30 min. Major salivary gland function was represented by saliva excretion after sialogogue stimulation with acidic material. The salivary excretion factor (SEF) was determined as the maximal excretion activity per gland as a function of the maximal uptake .
Parotid gland function measured as the SEF by salivary scintigraphy was evaluated before RT and at 1 and 2 years after RT. All patients received scintigraphy 1 year after RT, whereas only 25 patients (25/31, 81%) were examined 2 years after RT. Scintigraphy was not performed for six patients because of tumor recurrence (n=2) or patient refusal (n=4). The excretion response was analyzed per patient and subsequently per individual gland. The primary end point was set as the salivary flow ≤45% of the pre-RT value , which is equivalent to grade 3+ xerostomia based on the dry mouth subscales of LENT-SOMA criteria (subjective: xerostomia, analytic: salivary flow), where grade 1 is 76–95% of pre-RT salivary flow; grade 2, 51–75%; grade 3, 26–50%; and grade 4, 0–25% [7, 8].
NTCP data fitting
All DVH data for each patient were transferred to MATLAB (version R2009b), and the analysis, including 95% confidence intervals, was performed with SPSS for Windows (version 17.0; SPSS, Chicago, IL) using the same dataset and selected variables. The data were fit to the Lyman-Kutcher-Burman (LKB) NTCP model [15, 16]. The model quantitatively assesses the effects of both the radiation dose and the volume of the gland irradiated on the probability of radiation-induced changes in parotid gland function. Three parameters are represented in the sigmoidal dose–response curve: n, m, and TD50. The parameter n accounts for the volume effect of an organ: n was set to 1 in this study. The parameter m describes the slope of the dose–response curve, where decreasing m indicates increasing steepness of the slope. The TD50 is the dose for uniform irradiation of the whole or partial volume resulting in 50% probability of a complication. The NTCP is calculated from the equivalent uniform dose (EUD), assuming a sigmoidal (integrated normal distribution) relationship between the complication and EUD :
The EUD is defined as the uniform dose that would lead to the same level of tumor-cell killing as a non-uniform dose. Recently, the EUD concept has also been applied in normal tissues to evaluate the harm of a non-uniform dose distribution with the same result as a specific uniform dose. The formula for EUD is as follows:
where N is the number of voxels of the organ; D i is the dose of the i-th voxel; v i is the volume of the i-th voxel; and n is a parameter reflecting the biological properties of the organ related to its serial (0 < n << 1) or parallel structure (n of approximately 1). When n=1, the EUD is equal to the mean dose, as described previously [15, 20, 21]. The simplified LKB model represents the integral used in the Lyman formula as an exponential of a second-degree polynomial of dose.
The EORTC questionnaire was chosen for this research because it is one of the most widely implemented questionnaires, with more than 10 years of research invested to develop an integrated, modular approach. Moreover, it has been used in international clinical trials, and the Taiwan Chinese version is easily completed by our patients. The traditional Chinese version of the EORTC QLQ-H&N35 questionnaire, obtained from the Quality of Life Unit, EORTC Data Center in Brussels, Belgium [22, 23], was used for a prospective QoL survey. The primary endpoint (grade 3+ xerostomia) was defined as moderate to severe xerostomia 1 year after the completion of RT based on the QLQ-HN35 questionnaires. To ensure that xerostomia was induced primarily by the radiation treatment, patients with moderate to severe xerostomia at baseline were excluded from the analysis. All scales pertaining to the EORTC QLQ-H&N35 ranged from 0 to 100. A high score for a functional or global QoL scale represents a relatively high/healthy level of functioning or global QoL, whereas a high score for a symptom scale represents the presence of a symptom or problems [24–27]. All patients completed all questionnaires at three time points (before RT and at 3 and 12 months after RT), but only 17 patients (54.8%) completed the questionnaires at 2 years after RT. Thus, only the 1-year data were analyzed in the present study.
The Quantitative Analysis of Normal Tissue Effects in the Clinic (QUANTEC) guidelines are a recent concerted effort by the RT community to review and summarize normal tissue toxicity, which may suggest dose-volume treatment planning guidelines and likely reduce the rates of side effects. QUANTEC guidelines to limit the probability of severe xerostomia recommend that at least one parotid gland should receive a mean dose of ≤20 Gy or both parotid glands should receive a mean dose of ≤25 Gy [28, 29]. We performed a validation test of these guidelines using prospectively collected QoL and salivary scintigraphic datasets.
Spearman’s correlation was used to check the correlation between parotid gland excretion recovery at 1 year after RT and the mean parotid gland dose. To identify the patient- or treatment-related factors associated with parotid function recovery, we statistically analyzed age, tumor site, parotid mean dose, surgery, and chemotherapy using Spearman’s correlation and univariate and multivariate analysis. The mean scores and standard deviations of the QoL scales were calculated according to the EORTC QLQ scoring manual. To validate the QUANTEC constraints at the cut-off points of 20 Gy and 25 Gy, the negative predictive value (NPV) was checked for the rate of correctly predicting the lack of xerostomia. Positive predictive value (PPV) is the proportion of patients with xerostomia who are correctly diagnosed. The equations used for the PPV and NPV are listed in equations 4 and 5. Pearson’s chi-squared test was used to test goodness of fit and associations for both the SEF and QoL data. Overall performance was measured by Nagelkerke’s R2, which quantifies the amount of variation explained by the model. Model performance was also evaluated using measures for discriminative ability, including the area under the receiver operating characteristic curve (AUC). The Hosmer–Lemeshow test was used for calibration processing to test the goodness of fit for the hypothesis that the model and observed outcomes were in agreement . Values of p < 0.05 indicated statistical significance. All analyses were performed using SPSS 17.0.
The demographic and tumor characteristics of the study population are listed in Table 1. After a median follow-up of 46.8 months (range, 34.9–62.3 months), 30 (97%) of the 31 patients were still alive. One patient had died of other disease (lung cancer). Among the 30 surviving patients, 25 were still disease free, and five patients had distant metastasis. All scintigraphic examinations and QoL assessments were performed during disease-free periods. Dosimetry analysis showed that the ipsilateral parotid gland received a dose ranging from 26.9 to 74.8 Gy (mean, 51.7 Gy), and the contralateral lobe received 7.6 to 57.6 Gy (mean, 36.7 Gy).
The SEF values before RT showed a normal distribution (mean, 48.1% ± 18.2%). Parotid gland output varied considerably, from 19.2 to 72.8%. The relationships between parotid gland excretion recovery at 1 year and patient age, tumor site, parotid gland mean dose, surgery, and chemotherapy were analyzed (Table 2). Only parotid gland mean dose was significantly correlated with recovery of parotid gland function at 1 year (r=−0.807; r2=0.651; p < 0.001).
Figure 1 shows the observed SEF data and the fitted dose–response curves (LKB NTCP model) for the incidence of xerostomia at 1 year after completion of RT. The local fitted parameters are TD50=43.6 Gy (CI: 41.3-45.9 Gy) and m=0.18 (CI: 0.17-0.19). Dashed lines show the 95% confidence interval for the model fit to the SEF dataset (solid line). The incidences of grade 3+ xerostomia at 1 year were ~1% and ~2% for the recommended cut-off points of 20 and 25 Gy, respectively. Figure 2 presents the observed QoL data and fitted dose–response curves for the incidence of xerostomia at 1 year after RT. The fitted parameters are TD50=44.1 Gy (CI: 41.7-46.5 Gy) and m=0.11 (CI: 0.10-0.12). The curves for the two datasets are compared in Figure 3.
The positive and negative predictive values are listed on Table 3, the cut-off points were set for the QUANTEC guidelines and our developed NTCP models. The 20- and 25-Gy QUANTEC guidelines are also applied to the SEF and QoL data in Figures 4 and 5, respectively. The incidences of xerostomia for treatment plans meeting the QUANTEC guidelines occur precisely when the spared parotid mean dose is less than the 20- or 25-Gy cut-off values, giving a NPV of 100% with each dataset. As seen, the rate of xerostomia for plans meeting our developed NTCP models are low, for 43.6 Gy cut-off point, resulting in NPV's of 92% for the SEF data and 95% for the QoL data, and for 44.1 Gy cut-off point, resulting in NPV's of 88.9% for the SEF data and 95.5% for the QoL data, respectively.
Pearson’s chi-squared test demonstrated that the SEF and QoL data gave similar results (p = 0.241). Based on Nagelkerke’s R2, the overall NTCP model performance was similar between the patient-rated xerostomia-related QoL and the measured SEF values (Table 4). Furthermore, the discrimination based on the AUC was almost equal between the two datasets, and the Hosmer–Lemeshow test showed no significant disagreement between the results determined from each.
Parotid gland excretion recovery at 1 year and mean parotid gland dose were strongly correlated, based on Spearman’s correlation, and mean parotid gland dose was the only significant predictive factor for xerostomia. This finding differs from the report by Beetz et al. [30, 31], which proposed that multiple factors are likely to have separate impacts on xerostomia, although there may be no racial differences in the parotid gland response to irradiation.
Whole-mouth salivary function has been shown to be related to QoL determined by questionnaires . With LKB NTCP modeling in the present study, the TD50 for xerostomia 1-year after RT was 43.6 Gy in the SEF analysis and 44.1 Gy in the QoL analysis. Although these values are higher than that reported by Moiseenko et al. (32.4 Gy) , they are similar to the TD50 reported by Dijkema et al. (39.9 Gy), who analyzed the combined and updated results from two institutions . Deasy et al. suggested that the wide range of reported TD50 values (28.4 to 52 Gy) may result from differences in dose distribution, salivary measurement methods, segmentation, intragland sensitivity, and/or patient geographical location .
Xerostomia-specific questionnaires are reliable and valid for measuring patient-reported xerostomia . In the present study, QoL data were shown to be as valid as SEF values for NTCP parameter modeling. Based on Pearson’s chi-squared test, SEF and QoL data gave similar results. Furthermore, both the Nagelkerke’s R2, which describes overall performance, and the AUC demonstrated that both datasets produced similar results, and the Hosmer–Lemeshow test showed no significant disagreement between the results determined from the SEF and QoL data. No significant difference was noted regarding dose distributions to the parotid glands.
For the IMRT planning goal, the mean dose to each parotid gland should be as low as possible while providing the desired clinical target volume coverage . In our analysis, the incidence of grade 3+ xerostomia at 1 year was only ~1% or ~2% for the QUANTEC-recommended cut-off points of 20 Gy or 25 Gy, respectively. Hence, the severe xerostomia would usually be avoided when at least one parotid gland is spared to a mean dose ≤20 Gy or when both glands have been spared to a mean dose ≤25 Gy . A lower parotid mean dose also results in better QoL for patients .
Potential limitations of the present study include the low number of patients with xerostomia toxicity. Although the SEF values before RT were normally distributed, confirmation in a larger sample is needed to validate the NTCP model. Grade selection for the endpoint is another potential limitation in the present study. Choosing a lower grade of xerostomia may provide more valuable dose constraints for preventing complications, as even grade 2 xerostomia significantly diminishes QoL for patients . We used a previously definition for moderate to severe xerostomia based on the QLQ-HN35 questionnaire [2, 30, 31]. However, to our knowledge, no direct evidence exists to clarify this definition or to determine whether it is similar to the grade 3+ xerostomia definition by the subscales of LENT-SOMA criteria. Here, we showed the similar NTCP mapping results between moderate-to-severe xerostomia and LENT-SOMA subscales grade 3+ xerostomia. The practical implications of our results are validation of the use of a QoL form (EORTC QLQ-H&N35) as a surrogate for whole mouth salivary function, and also an important validation of previously proposed QAUNTEC treatment planning constraints to avoid xerostomia. Therefore, follow the QUANTEC guidelines have benefits to result in generic QoL improvement. Preserving more gland function should be pursued as a planning goal when consistent with adequate target dose coverage. Further researches will investigate as to whether new radiation techniques or different study cohorts (combining multiple institutional or cooperative group data sets) could be further validated this finding.
Our study shows the agreement between the NTCP parameter modeling based on SEF and QoL data and the QUANTEC guidelines, thus validating the QUANTEC cut-off values of 20 and 25 Gy. Based on these results, we believe that the QUANTEC 25/20-Gy spared-gland mean-dose guidelines are clinically useful for avoiding xerostomia in the HN cohort.
Bär W, Schwarz M, Alber M, Bos LJ, Mijnheer BJ, Rasch C, Schneider C, Nüsslin F, Damen EMF: A comparison of forward and inverse treatment planning for intensity-modulated radiotherapy of head and neck cancer. Radiother Oncol. 2003, 69: 251-258. 10.1016/j.radonc.2003.08.002.
Beetz I, Schilstra C, Burlage FR, Koken PW, Doornaert P, Bijl HP, Chouvalova O, Leemans CR, de Bock GH, Christianen MEMC: Development of NTCP models for head and neck cancer patients treated with three-dimensional conformal radiotherapy for xerostomia and sticky saliva: The role of dosimetric and clinical factors. Radiother Oncol. 2012, 105: 86-93. 10.1016/j.radonc.2011.05.010.
Lee TF, Chao PJ, Ting HM, Lo SH, Wang YW, Tuan CC, Fang FM, Su TJ: Comparative analysis of SmartArc-based dual arc volumetric-modulated arc radiotherapy (VMAT) versus intensity-modulated radiotherapy (IMRT) for nasopharyngeal carcinoma. J Appl Clin Med Phys. 2011, 12: 158-174.
Lee TF, Ting HM, Chao PJ, Fang FM: Dual Arc Volumetric-modulated Arc Radiotherapy (VMAT) of Nasopharyngeal Carcinomas: A Simultaneous Integrated Boost Treatment Plan Comparison with Intensity-modulated Radiotherapies and Single Arc VMAT. Clini Oncol. 2011, 24: 196-207.
Lee T-F, Fang F-M, Chao P-J, Su TJ, Wang LK, Leung SW: Dosimetric comparisons of helical tomotherapy and step-and-shoot intensity-modulated radiotherapy in nasopharyngeal carcinoma. Radiother Oncol. 2008, 89: 89-96. 10.1016/j.radonc.2008.05.010.
Pow EHN, Kwong DLW, McMillan AS, Wong MCM, Sham JST, Leung LHT, Leung WK: Xerostomia and quality of life after intensity-modulated radiotherapy vs. conventional radiotherapy for early-stage nasopharyngeal carcinoma: Initial report on a randomized controlled clinical trial. Int J Radiat Oncol Biol Phys. 2006, 66: 981-991. 10.1016/j.ijrobp.2006.06.013.
Pavi J, Denekamp J, Letschert J: LENT-SOMA scales for all anatomic sites. Int J Radiat Oncol Biol Phys. 1995, 31: 1049-1091.
Pavy JJ, Denekamp J, Letschert J, Littbrand B, Mornex F, Bernier J, Gonzales-Gonzales D, Horiot JC, Bolla M, Bartelink H: Late effects toxicity scoring: The soma scale. Int J Radiat Oncol Biol Phys. 1995, 31: 1043-1047. 10.1016/0360-3016(95)00059-8.
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-453. 10.1016/j.ijrobp.2006.05.002.
Eisbruch A, Kim HM, Terrell JE, Marsh LH, Dawson LA, Ship JA: Xerostomia and its predictors following parotid-sparing irradiation of head-and-neck cancer. Int J Radiat Oncol Biol Phys. 2001, 50: 695-704. 10.1016/S0360-3016(01)01512-7.
Kakoei S, Haghdoost AA, Rad M, Mohammadalizadeh S, Pourdamghan N, Nakhaei MMB: Xerostomia after Radiotherapy and its Effect on Quality of Life in Head and Neck Cancer Patients. Arch Iran Med. 2012, 15: 214-218.
Roesink JM, Schipper M, Busschers W, Raaijmakers CPJ, Terhaard CHJ: A comparison of mean parotid gland dose with measures of parotid gland function after radiotherapy for head-and-neck cancer: Implications for future trials. Int J Radiat Oncol Biol Phys. 2005, 63: 1006-1009. 10.1016/j.ijrobp.2005.04.023.
Blanco AI, Chao KSC, 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-1069. 10.1016/j.ijrobp.2004.12.076.
Saarilahti K, Kouri M, Collan J, Kangasmäki A, Atula T, Joensuu H, Tenhunen M: Sparing of the submandibular glands by intensity modulated radiotherapy in the treatment of head and neck cancer. Radiother Oncol. 2006, 78: 270-275. 10.1016/j.radonc.2006.02.017.
Niemierko A, Goitein M: Calculation of normal tissue complication probability and dose-volume histogram reduction schemes for tissues with a critical element architecture. Radiother Oncol. 1991, 20: 166-176. 10.1016/0167-8140(91)90093-V.
Semenenko V, Li X: Lyman–Kutcher–Burman NTCP model parameters for radiation pneumonitis and xerostomia based on combined analysis of published clinical data. Phys Med Biol. 2008, 53: 737-10.1088/0031-9155/53/3/014.
Lee N, Harris J, Garden A, Straube W, Bosch W, Morrison W, Quivey J, Thorstadt W, Jones C, Ang K: Phase II multi-institutional study of IMRT±chemotherapy for nasopharyngeal carcinoma (RTOG 0225): Preliminary results. Int J Radiat Oncol Biol Phys. 2007, 69: S13-S14.
Roesink JM, Moerland MA, Hoekstra A, Rijk PPV, Terhaard CHJ: Scintigraphic assessment of early and late parotid gland function after radiotherapy for head-and-neck cancer: a prospective study of dose–volume response relationships. Int J Radiat Oncol Biol Phys. 2004, 58: 1451-1460. 10.1016/j.ijrobp.2003.09.021.
Seppenwoolde Y, Lebesque JV, De Jaeger K, Belderbos JSA, Boersma LJ, Schilstra C, Henning GT, Hayman JA, Martel MK, Ten Haken RK: Comparing different NTCP models that predict the incidence of radiation pneumonitis. Int J Radiat Oncol Biol Phys. 2003, 55: 724-735. 10.1016/S0360-3016(02)03986-X.
Wu Q, Mohan R, Niemierko A, Schmidt-Ullrich R: Optimization of intensity-modulated radiotherapy plans based on the equivalent uniform dose. Int J Radiat Oncol Biol Phys. 2002, 52: 224-235. 10.1016/S0360-3016(01)02585-8.
Gulliford SL, Partridge M, Sydes MR, Webb S, Evans PM, Dearnaley DP: Parameters for the Lyman Kutcher Burman (LKB) model of Normal Tissue Complication Probability (NTCP) for specific rectal complications observed in clinical practise. Radiother Oncol. 2012, 102: 347-351. 10.1016/j.radonc.2011.10.022.
Bjordal K, Hammerlid E, Ahlner-Elmqvist M, de Graeff A, Boysen M, Evensen JF, Biörklund A, de Leeuw JRJ, Fayers PM, Jannert M: Quality of life in head and neck cancer patients: validation of the European Organization for Research and Treatment of Cancer Quality of Life Questionnaire-H&N35. J Clin Oncol. 1999, 17: 1008-1008.
Aaronson NK, Ahmedzai S, Bergman B, Bullinger M, Cull A, Duez NJ, Filiberti A, Flechtner H, Fleishman SB, de Haes JCJM: The European Organization for Research and Treatment of Cancer QLQ-C30: a quality-of-life instrument for use in international clinical trials in oncology. J Natio Cancer Instit. 1993, 85: 365-376. 10.1093/jnci/85.5.365.
Fang FM, Tsai WL, Lee TF, Liao K-C, Chen H-C, Hsu H-C: Multivariate analysis of quality of life outcome for nasopharyngeal carcinoma patients after treatment. Radiother Oncol. 2010, 97: 263-269. 10.1016/j.radonc.2010.05.022.
Leung SW, Lee TF, Chien CY, Chao PJ, Tsai WL, Fang FM: Health-related Quality of life in 640 head and neck cancer survivors after radiotherapy using EORTC QLQ-C30 and QLQ-H&N35 questionnaires. BMC cancer. 2011, 11: 128-10.1186/1471-2407-11-128.
Fang FM, Tsai WL, Chien CY, Chen HC, Hsu HC, Huang TL, Lee TF, Huang HY, Lee CH: Pretreatment quality of life as a predictor of distant metastasis and survival for patients with nasopharyngeal carcinoma. J Clini Oncol. 2010, 28: 4384-4389. 10.1200/JCO.2010.28.8324.
Huang TL, Tsai WL, Chien CY, Lee TF, Fang FM: Quality of life for head and neck cancer patients treated by combined modality therapy: the therapeutic benefit of technological advances in radiotherapy. Qual Life Res. 2010, 19: 1243-1254. 10.1007/s11136-010-9688-3.
Moiseenko V, Wu J, Hovan A, Saleh Z, Apte A, Deasy JO, Harrow S, Rabuka C, Muggli A, Thompson A: Treatment planning constraints to avoid xerostomia in head-and-neck radiotherapy: an independent test of QUANTEC criteria using a prospectively collected dataset. Int J Radiat Oncol Biol Phys. 2012, 82: 1108-1114. 10.1016/j.ijrobp.2011.04.020.
Deasy JO, Moiseenko V, Marks L, Chao KSC, Nam J, Eisbruch A: Radiotherapy dose–volume effects on salivary gland function. Int J Radiat Oncol Biol Phys. 2010, 76: S58-S63. 10.1016/j.ijrobp.2009.06.090.
Beetz I, Schilstra C, van Luijk P, Christianen MEMC, Doornaert P, Bijl HP, Chouvalova O, van den Heuvel ER, Steenbakkers RJHM, Langendijk JA: External validation of three dimensional conformal radiotherapy based NTCP models for patient-rated xerostomia and sticky saliva among patients treated with intensity modulated radiotherapy. Radiother Oncol. 2012, 105: 101-106. 10.1016/j.radonc.2012.03.004.
Beetz I, Schilstra C, van der Schaaf A, van den Heuvel ER, Doornaert P, van Luijk P, Vissink A, van der Laan BFAM, Leemans CR, Bijl HP: 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: 94-100. 10.1016/j.radonc.2011.11.006.
Dijkema T, Raaijmakers CPJ, Ten Haken RK, Roesink JM, Braam PM, Houweling AC, Moerland MA, Eisbruch A, Terhaard CHJ: Parotid gland function after radiotherapy: the combined Michigan and Utrecht experience. Int J Radiat Oncol Biol Phys. 2010, 78: 449-453. 10.1016/j.ijrobp.2009.07.1708.
Jabbari S, Kim HM, Feng M, Lin A, Tsien C, Elshaikh M, Terrel JE, Murdoch-Kinch C, Eisbruch A: Matched case–control study of quality of life and xerostomia after intensity-modulated radiotherapy or standard radiotherapy for head-and-neck cancer: Initial report. Int J Radiat Oncol Biol Phys. 2005, 63: 725-731. 10.1016/j.ijrobp.2005.02.045.
Ramaekers BLT, Joore MA, Grutters JPC, van den Ende P, Jong J, Houben R, Lambin P, Christianen M, Beetz I, Pijls-Johannesma M, Langendijk JA: The impact of late treatment-toxicity on generic health-related quality of life in head and neck cancer patients after radiotherapy. Oral Oncol. 2011, 47: 768-774. 10.1016/j.oraloncology.2011.05.012.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/12/567/prepub
This study was supported financially, in part, by grant from the National Science Council (NSC) of the Executive Yuan of the Republic of China. (NSC-101-2221-E-151-007-MY3). The authors thank the reviewers for their helpful comments on the original manuscript.
Part of this study was presented on the International Conference on Intelligent Informatics in Biology and Medicine (ICIIBM 2012).
TFL: original idea, study design, and writing of manuscript. PJC, HYW and PSC: statistical analysis. HCH and WCC: data collection and technical supports. All authors read and approved the final manuscript.