Abstract
Purpose/Objective(s):
The additional personnel and imaging procedures required for Adaptive Radiation Therapy (ART) pose a challenge for a broad implementation. We hypothesize that a change in transit fluence during the treatment course is correlated with the change of quality of life and thus can be used as a replanning trigger.
Materials/Methods:
Twenty-one head and neck cancer (HNC) patients filled out an MD Anderson Dysphagia Inventory (MDADI) questionnaire, before-and-after the radiotherapy treatment course. The transit fluence was measured by the Watchdog (WD) in-vivo portal dosimetry system. The patients were monitored with daily WD and weekly CBCTs. The region of interest (ROI) of each patient was defined as the outer contour of the patient between approximate spine levels C1 to C4, essentially the neck and mandible inside the beam’s eye view. The nth day integrated transit fluence change, Δϕn, and the volume change, ΔVROI, of the ROI of each patient was calculated from the corresponding WD and CBCT measurements. The correlation between MDADI scores and age, gender, planning mean dose to salivary glands <Dsg>, weight change ΔW, ΔVROI, and Δϕn, were analyzed using the ranked-Pearson correlation.
Results:
No statistically significant correlation was found for age, gender and ΔW. <Dsg> was found to have clinically important correlation with functional MDADI (ρ = −0.39, P = 0.081). ΔVROI was found to have statistically significant correlation of 0.44, 0.47 and 0.44 with global, physical and functional MDADI (P-value < 0.05). Δϕn was found to have statistically significant ranked-correlation (−0.46, −0.46 and −0.45) with physical, functional and total MDADI (P-value < 0.05).
Conclusion:
A transit fluence based decision support metric (DSM) is statistically correlated with the dysphagia risk. It can not only be used as an early signal in assisting clinicians in the ART patient selection for replanning, but also lowers the resource barrier of ART implementation.
Keywords
Introduction
Adaptive radiotherapy (ART), the adaption of the treatment plan to account for significant anatomical changes 1 -5 occurring during the course of radiation therapy, is a viable clinical strategy to both maintain the PTV coverage and to control the degradation risk of quality of life for head and neck cancer (HNC) patients. The degradation is a result of radiation induced conditions such as dysphagia and xerostomia 6 -8 due to unplanned dose increases to the salivary glands. To capture longitudinal anatomical changes, regular 3-dimensional imaging procedures, 1 -3,5,9,10 such as cone beam computed tomography (CBCT) and magnetic resonance imaging (MRI), are performed and reviewed during the course of treatment. Although this technique has demonstrated the effectiveness in reducing dose to the parotid glands (PG), 1,3,5 the high operational cost, such as the additional time required for imaging procedures and physician assessment, is preventing ART from being more broadly implemented. 5 As the scans are usually taken on only weekly, the best time window to adapt potentially may be missed. Furthermore, significant variability in the anatomic longitudinal change during a treatment course is observed among patients. 11,12 Both real-time and offline electronic portal image detector (EPID) based in-vivo dosimetry systems that utilize exit radiation, or transit fluence, from patients during treatments have been investigated and implemented. 13 -16 Our recent study has established the feasibility of using transit fluence to monitor volume change in HNC patients 17 without additional patient exposure or significant workload increase to the therapists. However, no prior study has investigated the relationship between the patients’ self-reported outcome and corresponding transit fluence change as of the time of this study. This study focuses on the investigation of the effectiveness of using transit fluence change as a replanning trigger in ART. To evaluate the effectiveness, we hypothesize that the change of transit fluence is associated with the changes in dose deposition in critical structures, such as the salivary glands and pharyngeal constrictor muscles (PCM), and the dose change is associated with the risk of dysphagia.
Method
Transit Fluence Monitoring
Transit fluence is the sum of the attenuated primary fluence and the scatter fluence of any treatment field leaving a patient. In this study, the transit fluence was measured with the Watchdog (WD) transit dosimetry system 16 utilizing the on board electronic portal imager (EPID) in portal dosimetry mode. About 3000 to 4000 cine images per treatment were captured by WD. Every pixel of all the images for a given treatment session were summed across all treatment arcs to give an integrated transit fluence. Twenty-four patients, who were on a retrospective protocol, were randomly selected and monitored with WD daily and CBCT weekly. All patients were treated with normal fractionated 70 Gy treatments. About 76% of the patients also received chemotherapy. All treatments were delivered with a 6 MV photon beam using a Volumetric Modulated Arc Therapy (VMAT) technique with millennium 120 (M120) multi-leaf collimators (MLC). Patients were immobilized with custom thermoplastic facemasks. A pair of kV-kV images was used to register the bony landmark to the digitally reconstructed radiographs from the planning computed tomography (CT) on daily treatment setup. The residual daily setup variation is typically small but contributes to the uncertainty of the measurements. The region of interest (ROI) of each patient was defined by the outer contour of the patient between approximate spine levels C1 and C4, essentially the neck and the mandible inside the beam’s eye view. The average integrated transit fluence intersecting the projected ROI of each patient at each fraction, ϕ, was calculated from the corresponding WD measurement. The CBCT, performed on the same day, was used to determine the volume of the ROI. The changes on the nth day relative to the day 0 baseline in the transit fluence, Δϕn, and the volume of the ROI, ΔVROI, were defined as the ϕn − ϕo and VROI, n − VROI,0 respectively. It has been reported that a change of 10% of neck separation, 18 approximately equivalent to 5% volume change, is associated with the increase risk of grade 2 xerostomia which is a major contributing factor to dysphagia. 19 In this study, a 5% volume change, as determined by the CBCT, was thus used as the ART replanning trigger point, Trig(ΔVROI, n). When this Trig(ΔVROI, n) was reached, it would initiate the replanning process which typically include planning review by the clinicians to determine the necessity of a new adapted plan. The logistic regression was performed between the ART replanning trigger point and Δϕn. − min(Δϕn) where min(Δϕn) is the minimum fluence change in the patient cohort. The fitted parameter, β1 and the intercept, β0, were used to determine the Δϕn replanning trigger point, Trig(Δϕn).
The correlation and the area under the curve (AUC) of the receiver operating characteristic (ROC) curve of the logistic regression with a 5% volume change trigger between Δϕn and ΔVROI were computed and compared with our previous study 17 to assess the reproducibility. Additional volumetric triggers, 0.0%, −2.5% and −7.5%, were also used to calculate assess the stability of the AUC.
Quality of Life (QoL)
Twenty-one of the 24 patients filled out an MD Anderson Dysphagia Inventory (MDADI) questionnaire, before and after the radiation treatment course. Each set of questions was sub-divided into 4 subsets 20 : global, emotional, physical and functional. In brief, these are respectively the patient’s self-reported overall assessment of their daily routine, their affective response, their self-perception, and the impact on daily activities of their swallowing disorder. Each was graded according to MDADI protocol. The score of each subset is defined as the sum of the scores of the questions in the corresponding subset. The total MDADI score is the sum of all the subset question scores. To reduce the impact of outliers, the association between each of the 4 subsets of MDADI and each relevant factor were analyzed using ranked Pearson correlation. These factors are age, gender, planned mean dose to a structure composed of parotids and sub-mandibular glands <Dsg>, weight change ΔW, ΔVROI, and Δϕn . The P-value of the ranked Pearson correlation factors, ρ, were used to assess the significance of the association. P-values of less 0.05 and 0.10 were used for the thresholds for the classification of the statistically significant and clinically important association 21 respectively. A lower post-treatment MDADI score means a decrease of QoL resulting from dysphagia 20,22 and a change of 10 points or greater means a clinically significant deterioration in QoL. 23,24
Results
Transit Fluence Monitoring
The change in the ROI transit fluence was found to be strongly correlated with ROI volume. Nine of the 24 patients exhibited larger than 5.0% ROI volumetric reduction at some point in their treatment courses (Range 6% to 11%). For these patients a corresponding increase in transit fluence was observed (range 2% to 7%). The time observed for the volume to decrease by 5.0% or more was variable, ranging from 10 to 29 days from the start of treatment. For patients with a ROI volume change of less than 5%, the increase in transit fluence was also smaller (less than 3.5%). Excluding the baseline points, a total of 110 pairs of ΔVROI and Δϕn were obtained. The correlation was found to be −0.84 with the P-value less than 0.001. Figure 1 shows the scatter plot of ΔVROI and Δϕn of all the patients The R2 of the linear regression between ΔVROI and Δϕn is 0.70 with a slope of the regression of −1.39. The χ2 statistic of the model is 17.4 (P-value < 0.001). The β0 and β1 of the logistic regression are −10.04 and 99.38 respectively. The corresponding P-values of β0 and β1 are 0.001 and 0.003 respectively indicating that both factors are statistically significant. The logarithmic odds ratio (OR) of Δϕn was found to be 87.0 with a 95% confident interval of [27.4, 146.6]. With 5% volume change threshold, the Trig(Δϕn) was found to be 4.6%. Using Trig (ΔVROI, n) as the ground truth and Trig(Δϕn), the receiver operating curve (ROC) of Δϕn was calculated. The AUC of the ROC is 0.91 (Figure 2). Volumetric threshold of 0.0%, −2.5% and −7.5% were also plotted in Figure 2 for comparison. The AUC were found to be 0.88, 0.89 and 0.95 for 0.0%, −2.5% and −7.5% volume threshold. This indicates a statistically significant association between the transit fluence change and the corresponding volumetric change based replanning trigger. The observed AUC did also not vary significantly within the threshold level of this study.

Scatter plot of the transit fluence variation and CBCT volumetric changes of all 24 the patients to assess the correlation between fluence and volume change during the course of treatment.

The solid line represents the ROC curve (AUC = 0.91) of Δϕn using 24 patients and −5.0% volumetric change. The green, blue and black dashed lines represent the volumetric change threshold of 0.0%, −2.5% and −7.5%, respectively; the dotted line represents the reference ROC curve with AUC of 0.50.
Quality of Life
Twenty-one of the 24 patients in this study answered the pre- and post-treatment MDADI questionnaire. The average and standard deviation of the total score were found to be 62.8 and 6.8 respectively. The gender mix was found to be 29% female and 71% male. The median age was 63 years old. The maximum and minimum were 86 and 32 years old respectively. The median weight and ROI volume were found to be 80.3 kg (range 57.3 to 118.3 kg) and 1530.4 cm3 (range 658.7 to 1942.8 cm3). The average weight change was found to be −6.0 kg (range +5.7 kg to −12.8 kg).
Table 1 shows the summary of the ranked correlation coefficients of age for global, emotion, physical, functional and total ranked correlation factors with the parameters and the corresponding P-values. Age, gender and weight change were found to have no statistically significant correlation with MDADI. The relationship between <Dsg> and total MDADI were also found to be insignificant, although a clinically important association was found between the functional subset and <Dsg>. The most important correlates with change in MDADI were found to be ΔVROI and Δϕn. Three of the 4 MDADI subsets, global, physical and functional, were found to be statistically significant correlated with ΔVROI. The MDADI total was also found to have a relational trend with ΔVROI with possible clinically importance (P-value = 0.073). Global, emotion, physical and functional subsets ranked correlation with Δϕn were found to be −0.35, 0.23, −0.46 and −0.46, respectively. The total MDADI change ranked correlation with Δϕn was found to be −0.45 (Figure 3A). Figure 3B shows the relationship between 2 of the MDADI subset changes, physical and functional, and the Δϕn. The significant ranked correlations were found to be between 0.036 and 0.038 with P-values less than 0.05 indicating a statistically significant correlation between transit fluence change and HNC patient self-reported outcome. Three out of the 21 patients were found to have a decrease of more than 10 points in their total scores, indicating a clinically significant deterioration of QoL. The associated Δϕn were found to be in range between 2.9% to 6.1%.

The relationship between transit fluence change, Δϕn, and (A) total MDAD score change; (B) functional and physical MDADI subset scores change. The negative score change indicates the increase risk of quality of life degradation.
Ranked Correlation Between MDADI and Factors.
Discussion
With the increased patient number since the previous study, 17 Δϕn continues to show statistically significant correlation with ΔVROI. Here the correlation coefficient was found to be slightly higher indicating a stronger relationship between Δϕn and ΔVROI. In the cases with ΔVROI less than 5%, amounting to about 63% of all cases, both CBCT and transit fluence change show only minor volumetric change. This is an indication that ART is likely not necessary for these patients. The observed wide variation in the replanning trigger point timing indicates that a pre-determined replanning time is not an optimal strategy, as also suggested by van Dijk et al 12 and Marzi et al. 11 Together with the AUC of the ROC, this data set indicates that Δϕn can give useful information to physicians to support their ART decisions in terms of adaptation timing and patient selection without incurring a large increase in clinical workload.
With the increased number of patients in the current study, about 71% of the transit fluence change can be attributed to the ROI volume change indicating a closer relationship between transit fluence and volume change from the previous study. 17 However, the day-to-day variability in patient setup and machine performance will contribute some noise to the transit fluence signal. As OR > 0 is associated with higher odds of outcome, 21 the high OR of Δϕn indicates that a higher fluence transmission is associated with higher probability 21 of replanning. The study also demonstrates that the Δϕn signal has good accuracy with an AUC of 0.91. The AUC of a typical 25 ultrasound and a CT are 0.81 and 0.93 respectively. As volumetric change is a good predictor for the degradation of QoL, change in transit fluence has the potential to predict the volumetric changes during RT and act as a DSM for ART.
Similar to other studies, 26,27 age and gender were found to have no significant correlation with the change of QoL. Surprisingly, weight change during treatment was also found to have no significant association in this study. <Dsg> and normal tissue complication probability 28,29 (NTCP) are very effective in predicting high grade xerostomia prior to treatment. 8,26,30,31 Early dysphagia was also found to be significantly related to early xerostomia. 32,33 Interestingly, our study found <Dsg> is not an effective predictor of QoL. Van Dijk et al 34 and You et al 18 have also found that <Dsg> is not an effective predictor for xerostomia during treatment. As the planned dose of these patients were already within accepted planning guidelines 8 to minimize xerostomia, <Dsg> was found to be a weak predictor of xerostomia risk in line with Galbry’s study. 35 These factors may contribute to the weak association between <Dsg> and the outcome in this study. As this study does not measure xerostomia directly, a more specific patient self-reporting measure, such as European Organization for Research and Treatment of Cancer Quality of Life Questionnaire—H&N35 (EORTC QLQ-H&N35), 36 and more sophisticated dosimetric metrics, such as dose gradient or parotid dose change, may be more useful. 37 Such investigation falls outside of the scope of the current study.
For ART, volumetric changes determined from parotid gland contours 12,22,27,38 and neck separation 18 are found to be stronger and earlier predictors of QoL degradation than dosimetric metrics. 12 Volumetric changes are also believed to be more strongly associated with stem cell sterilization. 39 The association between ΔVROI, i and the MDADI found in this study is consistent with these published findings. Interestingly, significant correlation was also found between Δϕn and the MDADI total score change. Looking at the subset decomposition, both DSM’s were also found to be associated with the physical and functional subsets of MDADI. This implies that using either DSM as an early ART replanning signal can potentially improve outcome in HNC treatment. Given the relatively less computational demanding implementation of Δϕn, it may hold an advantage over the more computationally intensive salivary specific version in the case of clinical and automation implementation. In terms of trigger level, the replanning trigger of 5%, based on the xerostomia study, may have been too relaxed as 2 of the 3 patients, who exhibited more than 10 points drop in MDADI, had less than 5.0% increase in the Δϕn. A lower Δϕn trigger, such as 2.5%, may be a more appropriate replanning trigger point without the significant decrease in the AUC of the ROC.
Although the current study is based on relatively small sample of 24 patients to demonstrate the feasibility of the transit fluence based ART replanning trigger DSM, the correlation with patient reported outcomes, associated with QoL degradation, was very encouraging. A larger data set will be beneficial to further verify the relationship between clinical outcome and the DSM. Automation of this analysis process can expedite the current labor-intensive ART workflow and broaden its clinical implementation. Dysphagia has also been reported to correlate significantly with the dose to PCM. 40,41 The relationship between the transit fluence change and the increase in dose deposition to the PCM will be further studied.
Conclusion
To the best of our knowledge, this is the first study to demonstrate the statistically significant correlation between the change in transit fluence and the change in quality of life which is associated with dysphagia. It can not only be used as an early signal in assisting clinicians in the ART patient selection for replanning, but also lowers the resource barrier of ART implementation.
Footnotes
Authors’ Note
Statistical analysis was performed by Lim SB, PhD. Our study was approved by Memorial Sloan Kettering Cancer Center Institution Review Board (approval no. 18-257). This is a retrospective study and does not require informed consent. Patient data will not be shared with outside parties.
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was partially supported by the MSK Cancer Center Support Grant/Core Grant (P30 CA008748).
