Impact of radiation dose distribution on nutritional supplementation needs in head and neck cancer radiotherapy: a voxel-based machine learning approach

Objectives: To investigate the relationship between nutritional supplementation and radiation dose to the pharyngeal constrictor muscles and larynx for head and neck (HN) cancer patients undergoing radiotherapy. Methods: We retrospectively analyzed radiotherapy (RT) dose for 231 HN cancer patients, focusing on the pharyngeal constrictors and larynx. We defined nutritional supplementation as feeding tube utilization or >10% weight loss from baseline within 90 days after radiotherapy completion. Using deformable image registration (DIR), we mapped each patient's anatomical structures to a reference coordinate system, and corresponding deformations were applied to dose matrices. Voxel doses were utilized as features for ridge logistic regression models, optimized through 5-fold cross-validation. Model performance was assessed with area under the curve of a receiver operating curve (AUC) and F1 score. We built and compared models using 1) pharyngeal constrictor voxels, 2) larynx voxels, 3) clinical factors and mean regional dose metrics, and 4) clinical factors and dose-volume histogram metrics. Test set AUCs were compared among the models, and feature importance was evaluated. Results: DIR of the pharyngeal constrictors and larynx yielded mean Dice coefficients of 0.80 and 0.84, respectively. Pharyngeal constrictors voxels and larynx voxel models had AUC of 0.88 and 0.82, respectively. Voxel-based dose modeling identified the superior to middle regions of the pharyngeal constrictors and the superior region of larynx as most predictive of feeding tube use/weight loss. Univariate analysis found treatment setting, treatment laterality, chemotherapy, baseline dysphagia, weight, and socioeconomic status predictive of outcome. An aggregated model using mean doses of pharyngeal constrictors and larynx sub regions had an AUC of 0.87 and the model using conventional DVH metrics had an AUC of 0.85 with p-value of 0.04. Feature importance calculations from the regional dose model indicated that mean doses to the superior-middle pharyngeal constrictor muscles followed by mean dose to the superior larynx were most predictive of nutritional supplementation. Conclusions: Machine learning modeling of voxel-level doses enables identification of sub regions within organs that correlate with toxicity. For HN radiotherapy, doses to the superior-middle pharyngeal constrictors are most predictive of feeding tube use/weight loss followed by the doses to superior portion of the larynx.