Dose/Volume histogram patterns in Salivary Gland subvolumes influence xerostomia injury and recovery

Xerostomia is a common consequence of radiotherapy in head and neck cancer. The objective was to compare the regional radiation dose distribution in patients that developed xerostomia within 6 months of radiotherapy and those recovered from xerostomia within 18 months post-radiotherapy. We developed a feature generation pipeline to extract dose volume histogram features from geometrically defined ipsilateral/contralateral parotid glands, submandibular glands, and oral cavity surrogates for each patient. Permutation tests with multiple comparisons were performed to assess the dose difference between injury vs. non-injury and recovery vs. non-recovery. Ridge logistic regression models were applied to predict injury and recovery using clinical features along with dose features (D10-D90) of the subvolumes extracted from oral cavity and salivary gland contours +3 mm peripheral shell. Model performances were assessed by the area under the receiver operating characteristic curve (AUC) using nested cross-validation. We found that different regional dose/volume metrics patterns exist for injury vs. recovery. Compared to injury, recovery has increased importance to the subvolumes receiving lower dose. Within the subvolumes, injury tends to have increased importance towards D10 from D90. This suggests that different threshold for xerostomia injury and recovery. Injury is induced by the subvolumes receiving higher dose, and the ability to recover can be preserved by further reducing the dose to subvolumes receiving lower dose.