Abstract
Background/Aim: The aim of this study was to develop a normal tissue complication probability (NTCP) model for trismus in head and neck cancer (HNC) patients treated with radiotherapy (RT). Patients and Methods: Prospective measurements of maximum inter-incisal opening (MIO) were performed at baseline and 6 months after definitive RT in 132 HNC patients. The primary endpoint of this study was defined when a patient fulfilled both of the following criteria: 1) MIO at 6 months after RT ≤35 mm and 2) MIO at 6 months after RT ≤80% of baseline MIO. Eleven clinical factors and a wide range of dosimetric factors (mean dose, maximum dose, V5, V10, V20, and V40) in twelve organs at risk (OARs) were chosen as candidate prognostic variables. Results: Thirty out of 132 patients (23%) developed the primary endpoint. Multivariate logistic regression analysis revealed that the mean dose to the contralateral mandible joint (p=0.001) and baseline MIO (p=0.027) were independent prognostic factors. Conclusion: A multivariable NTCP model for trismus in HNC patients treated with RT was established including the mean dose to contralateral mandible joint and baseline MIO.
Trismus after radical radiotherapy has a major impact on quality of life (QoL) of head and neck cancer (HNC) patients (1-4). Previous studies have shown that rehabilitation of radiation-induced trismus is difficult to treat (5-8). Trismus is considered to result from slowly progressive fibrosis of the mandible joints and muscles involved in mandibular movements (9-11). Some data exist on the relationship between dosimetric factors of organs at risk (OARs) and trismus in univariate analysis and multivariate analysis (12-19). However, limited data exist on the relation among dosimetric factors of OARs, clinical factors, and trismus in multivariate analysis (20). Therefore, it remains unclear which risk OARs, dose-volume parameters and patient factors contribute most to the development of trismus. Therefore, identification of the most important prognostic factors and establishment of a multivariable prediction model for trismus is of major importance. Therefore, the aim of this prospective study was to develop a multivariable normal tissue complication probability (NTCP) model for trismus in HNC patients treated with RT using dosimetric and clinical factors by multivariate logistic regression analysis.
Patients and Methods
Patients and eligibility criteria. Between July 2007 and August 2010, 132 HNC patients who received definitive RT with or without systemic treatment at the department of Radiation Oncology of the UMCG were included in this retrospective analysis of prospectively collected data. All patients were subjected to the standardised follow-up program (SFP), as previously described (21), including a prospective evaluation of the maximum inter-incisal opening (MIO) routinely before and 6 months after curative (chemo-) radiotherapy (CH)RT. To be included in the analysis, patients had to fulfill the following eligibility criteria: 1) HNC originating in the oral cavity, nasopharynx, oropharynx, hypopharynx, or larynx; 2) treated with definitive RT either alone or in combination with chemotherapy or cetuximab; 3) no previous surgery, RT and/or chemotherapy; 4) no previous malignancies; 5) no distant metastases; 6) planning-computed tomography (CT) and three-dimensional dose distributions available prior to and 6 months after completion of (CH)RT; 7) no efforts of reducing dose in potential organs at risks (OARs) for trismus; 8) alive with no recurrences 6 months after RT, and: 9) no rehabilitation or medication for trismus. The study population used for this analysis was composed of 132 patients who fulfilled all these eligibility criteria. The majority of the patients were male (77%) and the mean age of the study population was 62 years, ranging from 33 to 89 years. The number of the patients treated with concurrent systemic therapy was 40 (30%). The demographic and tumor characteristics of this study population are listed in Table I.
Endpoints. MIO, was assessed using a commercially available device named Terabyte® (Atos Medical, Hörby, Sweden). In patients with an edentulous maxilla and not wearing dentures, the distance between the incisal edge of the alveolar ridge in the mandibular central incisor of the right (location 41) and that in the maxillary central incisor of the right (location 11) can be measured. In edentulous patients wearing dentures, the distance between the upper and lower dentures (location 41 and 11) can be measured.
To evaluate trismus due to RT, MIO was measured at baseline and at 6 months after RT. The 6 months interval was chosen because radiation-induced trismus usually becomes apparent between 3 and 6 months after completion of RT (22).
Trismus, the primary endpoint of this study, was defined when a patient fulfilled both following criteria:
MIO at six months after RT ≤35 mm,
MIO at six months after RT ≤80% of baseline MIO.
The first criterion was chosen because it demonstrated the most significant cut-off point of trismus in a previous study and is considered clinically relevant (23). The second criterion was chosen because the mean MIO at 6 months after RT was 80% of that of baseline by the report of Wang et al. (24). In the first diagnosis, some patient's MIO was less than or equal to 35 mm. Therefore, the second criterion was needed based on the reference.
Treatment. RT was delivered with linear accelerators (LINIAC) using megavoltage equipment. In all patients, a planning CT scan was made in supine position. All patients were treated with three dimensional conformal RT (3D-CRT) or intensity modulated radiation therapy (IMRT). Patient position was fixed with a five-point individual thermoplastic mask (Posicast® thermoplastics, CIVICO, Orange, IA, USA) in combination with a standard head support (Posifix® supine headrest, CIVICO). Position verification was carried out by using a shrinking action level correction protocol (SAL-protocol), using an electronic portal imaging device (EPID). Patients with early glottic cancer were treated with a fractionation dose of 2.5 Gy (5 times/week) up to a total dose of 60 Gy in 5 weeks or with a fraction dose of 2.0 Gy (5 or 6 times/week) up to a total dose of 66 Gy. These patients were only irradiated at the primary site. The 26 laryngeal cancer patients received local therapy. Their data are an important part of this study as a control, although their OARs for trismus did not receive a significant dose. Patients treated with concomitant CHRT were treated with conventional fractionation (2.0 Gy per fraction, 5 times per week up to 70 Gy in 7 weeks). In case of primary RT of the more advanced cases, which were considered not eligible for CHRT, an accelerated schedule with concomitant boost technique was used, either or not combined with cetuximab. These patients were generally treated with 6 fractions per week with a second fraction on Friday afternoon with minimum interval of 6 h, up to a total dose of 70 Gy in 6 weeks. Most of the patients received bilateral elective irradiation of the neck nodes to a total dose of 46 Gy and a boost on the primary tumor and swelling lymph nodes to a total dose of 70 Gy.
Patients characteristics and treatment parameters.
Contouring of organs at risk. OARs potentially related to trismus were delineated according to an institutional atlas and according to the report by Teguh et al. (12) including the ipsilateral and contralateral medial pterygoid muscles, lateral pterygoid muscles, masseter muscles, mandible joints, mandible joints + 5 mm, and temporal muscles. OARs were delineated by radiation oncologists (M. M and HP. B). The dose-volume parameters (mean dose, maximum dose, V5, V10, V20, and V40) in each OARs were extracted.
Statistics. Pearson correlation matrixes were produced twice (Table II and Figure 1) in order to check for high correlation between potential prognostic determinants. In case of high Pearson correlation coefficients, these variables were reduced to avoid the problem of multicollinearity, which may negatively affect the generalizability of the model. Table II was used to decide the representative dose-volume parameter of each OARs. Figure 1 was used to decide the clinical factors and the dosimetric factors for univariate regression analysis.
Next, the variables after the reference of Pearson correlation matrixes (Table II and Figure 1) initially included in the univariate logistic regression model are shown in Table III. Then, a multivariate logistic regression analysis was carried out including only variables that were significant in the univariate logistic regression model. For each patient, predictive values (i.e., NTCP values) were calculated for each set of prognostic variables based on the regression coefficients according to the formula:
Results
Prevalence of trismus. At 6 months after completion of treatment, 30 out of 132 patients (23%) showed trismus six months after RT.
Variable selection procedure and NTCP models for trismus. We produced the first Pearson correlation matrix to identify the dose-volume parameters of all OARs that were strongly (Pearson correlation coefficient ≥0.8) or moderately (Pearson correlation coefficient ≥0.7) correlated (Table II). Strong or moderate correlation was observed among many dose-volume parameters within each OAR. Therefore, we decided to make the mean dose of each OARs represent the dose-volume parameters of each OARs, while the V5, V10, V20, V40 and maximum dose of each OARs were excluded from further analyses. The mean dose, in general, showed good correlations with other dose-volume parameters. We decided to not enter the mean dose to temporal muscles in the second Pearson correlation matrix, because the average of the mean doses to the ipsilateral and contralateral temporal muscles was relatively low, i.e., only 4.4 Gy and 3.6 Gy, respectively, and not expected to be important.
Next, the variables (which were reduced by the above procedure) with the second Pearson correlation matrix are shown in Figure 1. We selected the variables that showed less than four very strong correlations (Pearson correlation coefficient ≥0.85) for the univariate regression analysis. Eventually, clinical factors and three dosimetric factors, including the mean doses to contralateral mandible joint, contralateral masseter muscle and the ipsilateral medial pterygoid muscle were selected for entering the univariate regression analysis.
In the univariate logistic regression analysis (Table III), all three dosimetric factors showed a significant association with trismus. In addition, conventional RT [versus (vs.) accelerated RT], oral cavity, nasopharyngeal, oropharyngeal, or hypopharyngeal cancer (vs. laryngeal cancer), N2b-3 (vs. N0-2a), female (vs. male), Stage III-IV (vs. Stage I-II), and systemic therapy + RT (vs. RT alone) and baseline MIO were also significant factors associated with trismus in the univariate analysis.
A multivariate logistic regression analysis was carried out including only the significant variables shown above in the univariate logistic regression analysis. The multivariate regression analysis eventually revealed two independent prognostic factors, including the mean dose to the contralateral mandible joint and baseline MIO. The average Dmean to the contralateral mandible joint in patients with trismus was 18 Gy [95% confidence interval (CI)=12.4-23.6 Gy], which was significantly higher than that observed among those without trismus, which was 7.7 Gy (95%CI=5.8-9.7 Gy) (t-test: p=0.001). The average baseline MIO among patients with trismus was 38.3 mm (95% CI=35.4-41.3 mm), which was significantly lower compared to that observed among those without trismus, which was 43.2 mm (95%CI=41.4-44.9 mm) (t-test: p=0.008). The distribution of all 132 mean doses to the contralateral mandible joint and baseline MIO are shown in Table IV. The NTCP-value for each variable is shown in Table V. The bar graph of the NTCP model based on the two variables and the quick reference are shown in Figure 2 and Figure 3, respectively. The goodness of fit of this model was confirmed in Hosmer and Lemeshow test (p=0.182). The calibration of this model is shown in Figure 4 and confirmed the goodness of fit of this model. The NTCP-value for each individual patient can be calculated by the following regression formula:
Discussion
The innovation and intensification of RT techniques and schedules together with the addition of systemic treatment has led to improved prognosis of HNC (25-30). Consequently, the prevalence of patients suffering from late toxicity is rising and becomes increasingly relevant. Xerostomia and dysphagia have been considered as the representative toxicities affecting QoL (31-33). However, recently it has been shown that trismus may affect QoL in HNC patients as well (1-4). Morimoto et al. recorded that the prevalence of trismus after RT was 13% in Japanese HNC patients. The prevalence ranges from 5% to 50% (34-43). In the current study, the prevalence was within this range.
The first Pearson correlation matrix of the dosimetric factors (mean dose, maximum dose, V5, V10, V20 and V40) of the twelve organs at risk (OARs). Strong (≥0.8) and moderate (≥0.7) correlations are shown in bold and italics, respectively. The mean dose in each OARs was chosen for the representative of each OARs, while the maximum dose, V5, V10, V20 and V40 of each OARs were excluded from further analyses.
Fibrosis of the mandible joints and/or muscles involved in mouth opening is considered as the main cause of radiation-induced trismus. The mandible joint is covered with a capsule and divided into the upper and lower joint cavity by the articular disk. The articular disk consists of fibrocartilage. While opening the mouth, depression (hinge motion in lower joint cavity) and protrusion (sliding motion in upper joint cavity) in the mandible joint take place simultaneously (44). Movement of the mandible joint is extremely sensitive to relative minor changes of the anatomy and therefore this complicated movement could be disturbed by fibrosis of the irradiated mandible joint. This hypothesis is further confirmed by the findings of magnetic resonance imaging (MRI) in patients with trismus, which demonstrated muscle fibrosis, mandible joint, mandibular condyle sclerosis, mandibular ramus signal change, and joint capsular thickening (45). Fibroblasts, which may normally serve as the primary collagen-producing cells, are responsible for excess extracellular matrix collagen production, which may eventually lead to loss of smooth movement of the mandible joint. In chronic inflammation by irradiation, tumor necrosis factor alpha (TNFα) expression has been shown to be involved in the activation of macrophages in injured tissue, which leads to the release of downstream fibrogenic cytokines. These cytokines include fibroblast growth factor 2 (FGF2) and transforming growth factor beta 1 (TGFβ1). FGF2 is chemotactic and mitotic for fibroblasts, whereas TGFβ1 stimulates fibroblast proliferation and premature end-differentiation (9).
There are some studies on the relationship between trismus and dosimetric factors of OARs. Teguh et al. showed a significant correlation between the mean dose to the masseter muscles and the pterygoid muscles with trismus (12). Goldstein et al. found a linear relationship between dose to mandible joint and pterygoid muscle and MIO. In that study, functional impairment of the mandible joint already appeared at a relatively low dose level (from 14.9 Gy) (13). These investigators reported that in the univariate analysis the dose to pterygoid muscle, masseter muscle or mandibular joint was significantly associated with trismus.
The five reports demonstrated a significant dose-effect relationship between ipsilateral masseter or pterygoid muscle and trismus in the multivariate analysis (14-18). The three reports demonstrated that the clinical factors; tumor location, gender, baseline MIO, baseline age, dentition, free soft tissue transfer after surgery, reirradiation, chemotherapy, natural logarithm of time post-RT, and overall treatment time of RT correlated with trismus in multivariate analysis (43, 46, 47).
Using dosimetric and clinical factors, Kraaijenga et al. recently reported that the NTCP model was established by the mean dose to ipsilateral masseter muscle, the mean dose to ipsilateral medial pterygoid muscle, and the baseline MIO (20). Thus, limited data exist on the relation among dosimetric factors of OARs, clinical factors, and trismus in a multivariate NTCP model.
The second Pearson correlation matrix of the clinical and dosimetric variables for trismus. Very strong correlations (≥0.85) are colored in red. In the dosimetric factors, MDCMJ, MDIMPM, and MDCMM showed less than 4 red colors. The three dosimetric factors were chosen for univariate logistic regression analysis. In addition, all clinical factors were chosen for univariate analysis. Lo: Location; T: T classification; N: N classification; His: histology; Base: baseline maximum inter-incisal opening; Accel: accelerated radiotherapy; Bil: bilateral neck radiotherapy; Sys: systemic therapy; MDIMJ: mean dose to ipsilateral mandible joint; MDCMJ: mean dose to contralateral mandible joint; MDIMJ5: mean dose to ipsilateral mandible joint + 5 mm; MDCMJ5: mean dose to ipsilateral mandible joint + 5 mm; MDIMM: mean dose to ipsilateral masseter muscle; MDCMM: mean dose to contralateral masseter muscle; MDILPM: mean dose to ipsilateral lateral pterygoid muscle; MDCLPM: mean dose to contralateral lateral pterygoid muscle; MDIMPM: mean dose to ipsilateral medial pterygoid muscle; MDCMPM: mean dose to contralateral medial pterygoid muscle.
In the current study, the mean dose to contralateral mandible joint and baseline MIO were independent prognostic factors for trismus in the multivariate analysis. In addition, the NTCP model for trismus after RT was established with these two variables. Baseline MIO should be noted that it may be related to our definition of trismus. It might not be intuitive that mean dose to contralateral mandible joint was more significant than those of the ipsilateral masseter or the pterygoid muscle mentioned above. One possible reason for why the mean dose to the mandible joint was more important than that to other OARs is that the joint tissues are more sensitive to radiation induced-damage and/or that movement of the mandible joint is more complicated than muscle function, and are thus more prone to the effects of a lower dose (13). In addition, it could be hypothesized that the contralateral mandible joint played the most important role as this joint is the main compensator for trismus in the long-term process of fibrosis. We would like to emphasize that, to the best of our knowledge, this study is a unique prospective research that detected the significant clinical and dosimetric factors related to trismus using an actual measurement of MIO with multiple logistic regression analysis.
In order to avoid or improve trismus in HNC survivors, some strategies may be considered. First, based on the results of the current studies, constraint for the contralateral mandible joint could be employed, which may be different among individual patients based on their baseline MIO. For example, in a patient with a baseline MIO of 40 mm, the probability of trismus would be 33% if the dose to the contralateral mandible joint is 20 Gy, but can be decreased to 21% and 17% if the dose is reduced to 10 Gy and 5 Gy, respectively. Although there are a few reports on reducing trismus by IMRT (38, 48, 49), no optimal dose constraints for the mandible joints have been established. Our study adds information to such RT practice.
Distributions of the mean dose to the contralateral mandible joint and the baseline maximum inter-incisal opening.
Second, our prediction model may identify patients at high risk for developing trismus and who may benefit from a preventive exercise program early during or directly after treatment. Such rehabilitation programs are currently available for HNC patients (5, 6, 50-52).
Normal tissue complication probability (NTCP) model for trismus in head and neck cancer patients treated with radiotherapy.
Quick reference matrix of normal tissue complication probability (NTCP) for trismus in head and neck cancer patients treated with radiotherapy.
Third, pentoxifylline, which is a suppressor of TNFα (9, 10), or microcurrent therapy may be useful to prevent or relieve trismus in the same high-risk group for developing trismus (53). To verify efficacy of the above-mentioned approach, future prospective trials are warranted.
Conclusion
We developed a multivariable prediction model for trismus, after definitive RT with or without systemic treatment for HNC, consisting of the mean dose to contralateral mandible joint and baseline MIO. This NTCP model can be used to define dose constraint in clinical practice.
Multivariable logistic regression model for trismus.
Calibration of a predictive model for patient-related trismus based on two prognostic factors, which were the baseline of maximum inter-incisal opening and the mean dose to contralateral mandible joint. Dashed line indicates perfect agreement. Solid line represents the fit of the observed versus the expected probability.
Acknowledgements
This work was supported by a grant of the European Union (ALLEGRO-project), the Japan Society for Promotion Science (JSPS) Core-to-Core Program (number 23003) and JSPS KAKENHI Grant Number 18K15616.
Footnotes
Authors' Contributions
Conception and design: Henk P. Bijl, Arjen van der Schaaf, and Johannes A. Langendijk. Acquisition of data: Masahiro Morimoto, Henk P. Bijl, Roel J.H.M. Steenbakkers, Olga Chouvalova, and Johannes A. Langendijk. Analysis and interpretation of data: Masahiro Morimoto, Henk P. Bijl, Arjen van der Schaaf, Cheng-Jian Xu, Yasuo Yoshioka, Teruki Teshima, and Johannes A. Langendijk. Writing, review, and/or revision of the manuscript: Masahiro Morimoto, Yasuo Yoshioka, and Johannes A. Langendijk.
This article is freely accessible online.
Conflicts of Interest
The are no conflicts of interest to declare regarding this study.
- Received October 18, 2019.
- Revision received November 6, 2019.
- Accepted November 7, 2019.
- Copyright© 2019, International Institute of Anticancer Research (Dr. George J. Delinasios), All rights reserved