Evaluation of the clinical impact of the differences between planned and delivered dose in prostate cancer radiotherapy based on CT-on-rails IGRT and patient-reported outcome scores.

PURPOSE: To estimate the clinical impact of differences between delivered and planned dose using dose metrics and normal tissue complication probability (NTCP) modeling. METHODS: Forty-six consecutive patients with prostate adenocarcinoma between 2010 and 2015 treated with intensity-modulated radiation therapy (IMRT) and who had undergone computed tomography on rails imaging were included. Delivered doses to bladder and rectum were estimated using a contour-based deformable image registration method. The bladder and rectum NTCP were calculated using dose-response parameters applied to planned and delivered dose distributions. Seven urinary and gastrointestinal symptoms were prospectively collected using the validated prostate cancer symptom indices patient reported outcome (PRO) at pre-treatment, weekly treatment, and post-treatment follow-up visits. Correlations between planned and delivered doses against PRO were evaluated in this study. RESULTS: Planned mean doses to bladder and rectum were 44.9 ± 13.6 Gy and 42.8 ± 7.3 Gy, while delivered doses were 46.1 ± 13.4 Gy and 41.3 ± 8.7 Gy, respectively. D10cc for rectum was 64.1 ± 7.6 Gy for planned and 60.1 ± 9.3 Gy for delivered doses. NTCP values of treatment plan were 22.3% ± 8.4% and 12.6% ± 5.9%, while those for delivered doses were 23.2% ± 8.4% and 9.9% ± 8.3% for bladder and rectum, respectively. Seven of 25 patients with follow-up data showed urinary complications (28%) and three had rectal complications (12%). Correlations of NTCP values of planned and delivered doses with PRO follow-up data were random for bladder and moderate for rectum (0.68 and 0.67, respectively). CONCLUSION: Sensitivity of bladder to clinical variations of dose accumulation indicates that an automated solution based on a DIR that considers inter-fractional organ deformation could recommend intervention. This is intended to achieve additional rectum sparing in cases that indicate higher than expected dose accumulation early during patient treatment in order to prevent acute severity of bowel symptoms.

Evaluation of a commercial DIR platform for contour propagation in prostate cancer patients treated with IMRT/VMAT.

PURPOSE: To assess the performance and limitations of contour propagation with three commercial deformable image registration (DIR) algorithms using fractional scans of CT-on-rails (CTOR) and Cone Beam CT (CBCT) in image guided prostate therapy patients treated with IMRT/VMAT. METHODS: Twenty prostate cancer patients treated with IMRT/VMAT were selected for analysis. A total of 453 fractions across those patients were analyzed. Image data were imported into MIM (MIM Software, Inc., Cleveland, OH) and three DIR algorithms (DIR Profile, normalized intensity-based (NIB) and shadowed NIB DIR algorithms) were applied to deformably register each fraction with the planning CT. Manually drawn contours of bladder and rectum were utilized for comparison against the DIR propagated contours in each fraction. Four metrics were utilized in the evaluation of contour similarity, the Hausdorff Distance (HD), Mean Distance to Agreement (MDA), Dice Similarity Coefficient (DSC), and Jaccard indices. A subfactor analysis was performed per modality (CTOR vs. CBCT) and time (fraction). Point estimates and 95% confidence intervals were assessed via a Linear Mixed Effect model for the contour similarity metrics. RESULTS: No statistically significant differences were observed between the DIR Profile and NIB algorithms. However, statistically significant differences were observed between the shadowed NIB and NIB algorithms for some of the DIR evaluation metrics. The Hausdorff Distance calculation showed the NIB propagated contours vs. shadowed NIB propagated contours against the manual contours were 14.82 mm vs. 8.34 mm for bladder and 15.87 mm vs. 11 mm for rectum, respectively. Similarly, the Mean Distance to Agreement calculation comparing the NIB propagated contours vs. shadowed NIB propagated contours against the manual contours were 2.43 mm vs. 0.98 mm for bladder and 2.57 mm vs. 1.00 mm for rectum, respectively. The Dice Similarity Coefficients comparing the NIB propagated contours and shadowed NIB propagated contours against the manual contours were 0.844 against 0.936 for bladder and 0.772 against 0.907 for rectum, respectively. The Jaccard indices comparing the NIB propagated contours and shadowed NIB propagated contours against the manual contours were 0.749 against 0.884 for bladder and 0.637 against 0.831 for rectum, respectively. The shadowed NIB DIR, which showed the closest agreement with the manual contours performed significantly better than the DIR Profile in all the comparisons. The OAR with the greatest agreement varied substantially across patients and image guided radiation therapy (IGRT) modality. Intra-patient variability of contour metric evaluation was insignificant across all the DIR algorithms. Statistical significance at α = 0.05 was observed for manual vs. deformably propagated contours for bladder for all the metrics except Hausdorff Distance (P = 0.01 for MDA, P = 0.02 for DSC, P = 0.01 for Jaccard), whereas the corresponding values for rectum were: P = 0.03 for HD, P = 0.01 for MDA, P < 0.01 for DSC, P < 0.01 for Jaccard. The performance of the different metrics varied slightly across the fractions of each patient, which indicates that weekly contour propagation models provide a reasonable approximation of the daily contour propagation models. CONCLUSION: The high variance of Hausdorff Distance across all automated methods for bladder indicates widely variable agreement across fractions for all patients. Lower variance across all modalities, methods, and metrics were observed for rectum. The shadowed NIB propagated contours were substantially more similar to the manual contours than the DIR Profile or NIB contours for both the CTOR and CBCT imaging modalities. The relationship of each algorithm to similarity with manual contours is consistent across all observed metrics and organs. Screening of image guidance for substantial differences in bladder and rectal filling compared with the planning CT reference could aid in identifying fractions for which automated DIR would prove insufficient.

Assessment of PlanIQ Feasibility DVH for head and neck treatment planning.

INTRODUCTION: Designing a radiation plan that optimally delivers both target coverage and normal tissue sparing is challenging. There are limited tools to determine what is dosimetrically achievable and frequently the experience of the planner/physician is relied upon to make these determinations. PlanIQ software provides a tool that uses target and organ at risk (OAR) geometry to indicate the difficulty of achieving different points for organ dose-volume histograms (DVH). We hypothesized that PlanIQ Feasibility DVH may aid planners in reducing dose to OARs. METHODS AND MATERIALS: Clinically delivered head and neck treatments (clinical plan) were re-planned (re-plan) putting high emphasis on maximally sparing the contralateral parotid gland, contralateral submandibular gland, and larynx while maintaining routine clinical dosimetric objectives. The planner was blinded to the results of the clinically delivered plan as well as the Feasibility DVHs from PlanIQ. The re-plan treatments were designed using 3-arc VMAT in Raystation (RaySearch Laboratories, Sweden). The planner was then given the results from the PlanIQ Feasibility DVH analysis and developed an additional plan incorporating this information using 4-arc VMAT (IQ plan). The DVHs across the three treatment plans were compared with what was deemed "impossible" by PlanIQ's Feasibility DVH (Impossible DVH). The impossible DVH (red) is defined as the DVH generated using the minimal dose that any voxel outside the targets must receive given 100% target coverage. RESULTS: The re-plans performed blinded to PlanIQ Feasibilty DVH achieved superior sparing of aforementioned OARs compared to the clinically delivered plans and resulted in discrepancies from the impossible DVHs by an average of 200-700 cGy. Using the PlanIQ Feasibility DVH led to additionalOAR sparing compared to both the re-plans and clinical plans and reduced the discrepancies from the impossible DVHs to an average of approximately 100 cGy. The dose reduction from clinical to re-plan and re-plan to IQ plan were significantly different even when taking into account multiple hypothesis testing for both the contralateral parotid and the larynx (P < 0.004 for all comparisons). No significant differences were observed between the three plans for the contralateral parotid when considering multiple hypothesis testing. CONCLUSIONS: Clinical treatment plans and blinded re-plans were found to suboptimally spare OARs. PlanIQ could aid planners in generating treatment plans that push the limits of OAR sparing while maintaining routine clinical target coverage goals.

Thermal dosimetry characteristics of deep regional heating of non-muscle invasive bladder cancer.

PURPOSE: The aim of this paper is to report thermal dosimetry characteristics of external deep regional pelvic hyperthermia combined with intravesical mitomycin C (MMC) for treating bladder cancer following transurethral resection of bladder tumour, and to use thermal data to evaluate reliability of delivering the prescribed hyperthermia dose to bladder tissue. MATERIALS AND METHODS: A total of 14 patients were treated with MMC and deep regional hyperthermia (BSD-2000, Sigma Ellipse or Sigma 60). The hyperthermia objective was 42° ± 2 °C to bladder tissue for ≥40 min per treatment. Temperatures were monitored with thermistor probes and recorded values were used to calculate thermal dose and evaluate treatment. Anatomical characteristics were examined for possible correlations with heating. RESULTS: Combined with BSD-2000 standard treatment planning and patient feedback, real-time temperature monitoring allowed thermal steering of heat sufficient to attain the prescribed thermal dose to bladder tissue within patient tolerance in 91.6% of treatments. Mean treatment time for bladder tissue >40 °C was 61.9 ± 11.4 min and mean thermal dose was 21.3 ± 16.5 CEM43. Average thermal doses obtained in normal tissues were 1.6 ± 1.2 CEM43 for the rectum and 0.8 ± 1.3 CEM43 in superficial normal tissues. No significant correlation was seen between patient anatomical characteristics and thermal dose achieved in bladder tissue. CONCLUSIONS: This study demonstrates that a hyperthermia prescription of 42° ± 2 °C for 40-60 min can be delivered safely to bladder tissue with external radiofrequency phased array applicators for a typical range of patient sizes. Using the available thermometry and treatment planning, the BSD-2000 hyperthermia system was shown to be an effective method of focusing heat regionally around the bladder with good patient tolerance.

Quantifying Regional Radiation-Induced Lung Injury in Patients Using Hyperpolarized 129Xe Gas Exchange Magnetic Resonance Imaging.

PURPOSE: Radiation-induced lung injury has been shown to alter regional ventilation and perfusion in the lung. However, changes in regional pulmonary gas exchange have not previously been measured. METHODS AND MATERIALS: Ten patients receiving conventional radiation therapy (RT) for lung cancer underwent pre-RT and 3-month post-RT magnetic resonance imaging (MRI) using an established hyperpolarized 129Xe gas exchange technique to map lung function. Four patients underwent an additional 8-month post-RT MRI. The MR signal from inhaled xenon was measured in the following 3 pulmonary compartments: the lung airspaces, the alveolar membrane tissue, and the pulmonary capillaries (interacting with red blood cells [RBCs]). Thoracic 1H MRI scans were acquired, and deformable registration was used to transfer 129Xe functional maps to the RT planning computed tomography scan. The RT-associated changes in ventilation, membrane uptake, and RBC transfer were computed as a function of regional lung dose (equivalent dose in 2-Gy fractions). Pearson correlations and t tests were used to determine statistical significance, and weighted sum of squares linear regression subsequently characterized the dose dependence of each functional component. The pulmonary function testing metrics of forced vital capacity and diffusing capacity for carbon monoxide were also acquired at each time point. RESULTS: Compared with pre-RT baseline, 3-month post-RT ventilation decreased by an average of -0.24 ± 0.05%/Gy (ρ = -0.88; P < .001), membrane uptake increased by 0.69 ± 0.14%/Gy (ρ = 0.94; P < .001), and RBC transfer decreased by -0.41 ± 0.06%/Gy (ρ = -0.92; P < .001). Membrane uptake maintained a strong positive correlation with regional dose at 8 months post-RT, demonstrating an increase of 0.73 ± 0.11%/Gy (ρ = 0.92; P = .006). Changes in membrane uptake and RBC transfer appeared greater in magnitude (%/Gy) for individuals with low heterogeneity in their baseline lung function. An increase in whole-lung membrane uptake showed moderate correlation with decreases in forced vital capacity (ρ = -0.50; P = .17) and diffusing capacity for carbon monoxide (ρ = -0.44; P = .23), with neither correlation reaching statistical significance. CONCLUSIONS: Hyperpolarized 129Xe MRI measured and quantified regional, RT-associated, dose-dependent changes in pulmonary gas exchange. This tool could enable future work to improve our understanding and management of radiation-induced lung injury.

PET with 62Cu-ATSM and 62Cu-PTSM is a useful imaging tool for hypoxia and perfusion in pulmonary lesions.

OBJECTIVE: Hypoxia is a characteristic of many tumors and portends a worse prognosis in lung, cervical, prostate, and rectal cancers. Unlike the others, lung cancers present a unique challenge in measuring hypoxia, with invasive biopsies and higher rates of complications. Noninvasive imaging studies detecting hypoxia using isotopes of copper-diacetyl-bis(N4-methylthiosemicarbazone) ((62)Cu-ATSM) have predicted prognosis and treatment outcomes in some small feasibility trials. These images, however, may not identify all areas of hypoxia. Hence, we hypothesize that the addition of another PET imaging agent, copper-pyruvaldehyde-bis(N4-methylthiosemicarbazone) ((62)Cu-PTSM), which can detect areas of perfusion, can augment the information obtained in (62)Cu-ATSM PET scans. SUBJECTS AND METHODS: To characterize tumors on the basis of both perfusion and hypoxia, 10 patients were studied using both (62)Cu-ATSM and (62)Cu-PTSM PET scans. In addition, proteomic arrays looking at specific proangiogenic, survival, and proinflammatory targets were assessed. RESULTS: Six of 10 patients had evaluable PET scans. Our initial experience of characterizing lung tumor hypoxia using (62)Cu-ATSM and (62)Cu-PTSM PET scans showed that visualization of areas with hypoxia normalized for perfusion is feasible. All studied tumors exhibited some hypoxia. Despite the small sample size, a positive relationship was noted between epidermal growth factor levels and (62)Cu-ATSM-detected hypoxia. CONCLUSION: This initial series of (62)Cu-ATSM and (62)Cu-PTSM PET scans shows that evaluating lung masses by visualizing hypoxia and perfusion is a feasible and novel technique to provide more information. Further investigation is warranted to assess the potential role of (62)Cu-ATSM and (62)Cu-PTSM PET techniques combined with proteomics as alternatives to invasive biopsy techniques in clinical care.

Augmenting Quality Assurance Measures in Treatment Review with Machine Learning in Radiation Oncology.

Purpose: Pretreatment quality assurance (QA) of treatment plans often requires a high cognitive workload and considerable time expenditure. This study explores the use of machine learning to classify pretreatment chart check QA for a given radiation plan as difficult or less difficult, thereby alerting the physicists to increase scrutiny on difficult plans. Methods and Materials: Pretreatment QA data were collected for 973 cases between July 2018 and October 2020. The outcome variable, a degree of difficulty, was collected as a subjective rating by physicists who performed the pretreatment chart checks. Potential features were identified based on clinical relevance, contribution to plan complexity, and QA metrics. Five machine learning models were developed: support vector machine, random forest classifier, adaboost classifier, decision tree classifier, and neural network. These were incorporated into a voting classifier, where at least 2 algorithms needed to predict a case as difficult for it to be classified as such. Sensitivity analyses were conducted to evaluate feature importance. Results: The voting classifier achieved an overall accuracy of 77.4% on the test set, with 76.5% accuracy on difficult cases and 78.4% accuracy on less difficult cases. Sensitivity analysis showed features associated with plan complexity (number of fractions, dose per monitor unit, number of planning structures, and number of image sets) and clinical relevance (patient age) were sensitive across at least 3 algorithms. Conclusions: This approach can be used to equitably allocate plans to physicists rather than randomly allocate them, potentially improving pretreatment chart check effectiveness by reducing errors propagating downstream.

Idiopathic Pneumonitis Syndrome After Total Body Irradiation in Pediatric Patients Undergoing Myeloablative Hematopoietic Stem Cell Transplantation: A PENTEC Comprehensive Review.

PURPOSE: Pulmonary complications, especially idiopathic pneumonitis syndrome (IPS), are potentially life altering or fatal sequelae of hematopoietic cell transplantation (HCT). Total body irradiation (TBI) as part of the conditioning regimen has been implicated in IPS. A comprehensive PENTEC (Pediatric Normal Tissues in the Clinic) review was performed to increase our understanding of the role of TBI in the development of acute, noninfectious IPS. METHODS AND MATERIALS: A systematic literature search was conducted using the MEDLINE, PubMed, and Cochrane library databases for articles describing pulmonary toxicity in children treated with HCT. Data pertaining to TBI and pulmonary endpoints were extracted. Risk of IPS was analyzed in relation to patient age, TBI dose, fractionation, dose rate, lung shielding, timing, and type of transplant, with the goal to better understand factors associated with this complication in children undergoing HCT. A logistic regression model was developed using a subset of studies with comparable transplant regimens and sufficient TBI data. RESULTS: Six studies met criteria for modeling of the correlation of TBI parameters with IPS; all consisted of pediatric patients undergoing allogeneic HCT with a cyclophosphamide-based chemotherapy regimen. IPS was variably defined, but all studies that reported IPS were included in this analysis. The mean incidence of post-HCT IPS was 16% (range, 4%-41%). Mortality from IPS, when it occurred, was high (median, 50%; range, 45%-100%). Fractionated TBI prescription doses encompassed a narrow range of 9 to 14 Gy. Many differing TBI methods were reported, and there was an absence of 3-dimensional dose analysis of lung blocking techniques. Thus, a univariate correlation between IPS and total TBI dose, dose fractionation, dose rate, or TBI technique could not be made. However, a model, built from these studies based on prescribed dose using a normalized dose parameter of equivalent dose in 2-Gy fractions (EQD2), adjusted for dose rate, suggested correlation with the development of IPS (P = .0004). The model-predicted odds ratio for IPS was 24.3 Gy-1 (95% confidence interval, 7.0-84.3). Use of TBI lung dose metrics (eg, midlung point dose) could not be successfully modeled, potentially because of dosimetric uncertainties in the actual delivered volumetric lung dose and imperfections in our modeling process. CONCLUSIONS: This PENTEC report is a comprehensive review of IPS in pediatric patients receiving fractionated TBI regimens for allogenic HCT. IPS was not clearly associated with 1 single TBI factor. Modeling using dose-rate adjusted EQD2 showed a response with IPS for allogeneic HCT using a cyclophosphamide-based chemotherapy regimen. Therefore, this model suggests IPS mitigation strategies can focus on not just the dose and dose per fraction but also the dose rate used in TBI. More data are needed to confirm this model and to determine the influence of chemotherapy regimens and contribution from graft-versus-host disease. The presence of confounding variables (eg, systemic chemotherapies) that affect risk, the narrow range of fractionated TBI doses found in the literature, and limitations of other reported data (eg, lung point dose) may have prevented a more straightforward link between IPS and total dose from being observed.

Utility of treatment planning for thermochemotherapy treatment of nonmuscle invasive bladder carcinoma.

PURPOSE: A recently completed Phase I clinical trial combined concurrent Mitomycin-C chemotherapy with deep regional heating using BSD-2000 Sigma-Ellipse applicator (BSD Corporation, Salt Lake City, UT, U.S.A.) for the treatment of nonmuscle invasive bladder cancer. This work presents a new treatment planning approach, and demonstrates potential impact of this approach on improvement of treatment quality. METHODS: This study retrospectively analyzes a subset of five patients on the trial. For each treatment, expert operators selected "clinical-optimal" settings based on simple model calculation on the BSD-2000 control console. Computed tomography (CT) scans acquired prior to treatment were segmented to create finite element patient models for retrospective simulations with Sigma-HyperPlan (Dr. Sennewald Medizintechnik GmbH, Munchen, Germany). Since Sigma-HyperPlan does not account for the convective nature of heat transfer within a fluid filled bladder, an effective thermal conductivity for bladder was introduced. This effective thermal conductivity value was determined by comparing simulation results with clinical measurements of bladder and rectum temperatures. Regions of predicted high temperature in normal tissues were compared with patient complaints during treatment. Treatment results using "computed-optimal" settings from the planning system were compared with clinical results using clinical-optimal settings to evaluate potential of treatment improvement by reducing hot spot volume. RESULTS: For all five patients, retrospective treatment planning indicated improved matches between simulated and measured bladder temperatures with increasing effective thermal conductivity. The differences were mostly within 1.3 °C when using an effective thermal conductivity value above 10 W/K/m. Changes in effective bladder thermal conductivity affected surrounding normal tissues within a distance of ∼1.5 cm from the bladder wall. Rectal temperature differences between simulation and measurement were large due to sensitivity to the sampling locations in rectum. The predicted bladder T90 correlated well with single-point bladder temperature measurement. Hot spot locations predicted by the simulation agreed qualitatively with patient complaints during treatment. Furthermore, comparison between the temperature distributions with clinical and computed-optimal settings demonstrated that the computed-optimal settings resulted in substantially reduced hot spot volumes. CONCLUSIONS: Determination of an effective thermal conductivity value for fluid filled bladder was essential for matching simulation and treatment temperatures. Prospectively planning patients using the effective thermal conductivity determined in this work can potentially improve treatment efficacy (compared to manual operator adjustments) by potentially lower discomfort from reduced hot spots in normal tissue.

Assessing the impact of radiation-induced changes in soft tissue density ∕ thickness on the study of radiation-induced perfusion changes in the lung and heart.

PURPOSE: Abnormalities in single photon emission computed tomography (SPECT) perfusion within the lung and heart are often detected following radiation for tumors in∕around the thorax (e.g., lung cancer or left-sided breast cancer). The presence of SPECT perfusion defects is determined by comparing pre- and post-RT SPECT images. However, RT may increase the density of the soft tissue surrounding the lung∕heart (e.g., chest wall∕breast) that could possibly lead to an "apparent" SPECT perfusion defect due to increased attenuation of emitted photons. Further, increases in tissue effective depth will also increase SPECT photon attenuation and may lead to "apparent" SPECT perfusion defects. The authors herein quantitatively assess the degree of density changes and effective depth in soft tissues following radiation in a series of patients on a prospective clinical study. METHODS: Patients receiving thoracic RT were enrolled on a prospective clinical study including pre- and post-RT thoracic computed tomography (CT) scans. Using image registration, changes in tissue density and effective depth within the soft tissues were quantified (as absolute change in average CT Hounsfield units, HU, or tissue thickness, cm). Changes in HU and tissue effective depth were considered as a continuous variable. The potential impact of these tissue changes on SPECT images was estimated using simulation data from a female SPECT thorax phantom with varying tissue densities. RESULTS: Pre- and serial post-RT CT images were quantitatively studied in 23 patients (4 breast cancer, 19 lung cancer). Data were generated from soft tissue regions receiving doses of 20-50 Gy. The average increase in density of the chest was 5 HU (range 46 to -69). The average change in breast density was a decrease of -1 HU (range 13 to -13). There was no apparent dose response in neither the dichotomous nor the continuous analysis. Seventy seven soft tissue contours were created for 19 lung cancer patients. The average change in tissue effective depth was +0.2 cm (range -1.9 to 2.2 cm). The changes in HU represent a <2% average change in tissue density. Based on simulation, the small degree of density and tissue effective depth change is unlikely to yield meaningful changes in either SPECT lung or heart perfusion. CONCLUSIONS: RT doses of 20-50 Gy can cause up to a 46 HU increase in soft tissue density 6 months post-RT. Post-RT soft tissue effective depth may increase by 2.0 cm. These modest increases in soft tissue density and effective depth are unlikely to be responsible for the perfusion changes seen on post-RT SPECT lung or heart scans. Further, there was no clear dose response of the soft tissue density changes. Ultimately, the authors findings suggest that prior perfusion reports do reflect changes in the physiology of the lungs and heart.

Impact of Oral Cavity Dosimetry on Patient Reported Xerostomia and Dysgeusia in the Setting of Deintensified Chemoradiotherapy.

Purpose: To determine the relationship between mean oral cavity (OC) dose (treated as a singular organ at risk) to patient reported xerostomia and dysgeusia. In addition, we will examine the relationship between oral cavity substructure doses to patient reported xerostomia and dysgeusia. All patients were treated in the setting of deintensification (60 Gy). Methods and Materials: In the study, 184 and 177 prospectively enrolled patients for de-escalated chemoradiotherapy (CRT) for human papillomavirus (HPV)-positive oropharyngeal cancer submitted PROs at 6 and 12 months, respectively using Patient-Reported Outcomes version of the Common Terminology Criteria for Adverse Events questionnaire. Patient's OC consisting of the following substructures were segmented: oral tongue, base of tongue, floor of mouth, hard and soft palate, cheek mucosa, and upper and lower lip mucosa. Ordinal logistic regression (no/mild vs moderate vs severe/very severe symptoms) was used to compare organs at risk dosimetry to patient reported xerostomia and dysgeusia at 6 and 12 months. Multivariate ordinal logistic regression models were generated. Results: Mean dose to the contralateral parotid (P = .04), OC (P = .04), and baseline patient reported xerostomia (P = .009) were significantly associated with xerostomia severity at 6 months. Only baseline xerostomia (P = .02) and mean dose to the contralateral submandibular gland (P = .0001) were significantly associated with xerostomia severity at 12 months. The only significant factor related to dysgeusia at either time point was mean dose to the OC at 12 months (P = .009). On examining substructures, the mean dose to the floor of mouth was implicated for the dose relationship to 6-month xerostomia (P = .04), and the oral tongue was found to be implicated for the relationship for 12-month dysgeusia (P = .04). Conclusions: The mean dose to the OC was found to relate to xerostomia symptoms at 6 months post-CRT and dysgeusia symptoms at 12 months post-CRT. The mean dose to the floor of mouth and oral tongue appeared to drive this relationship for xerostomia and dysgeusia symptoms, respectively. This work suggests the floor of mouth and oral tongue should be prioritized during planning over the rest of the OC. The effect of OC dose relative to other salivary structures for xerostomia appeared to depend on time post-CRT.

Clinical Use of A Priori Knowledge of Organ-At-Risk Sparing During Radiation Therapy Treatment for Oropharyngeal Cancer: Dosimetric and Patient Reported Outcome Improvements.

PURPOSE: This study aimed to prospectively assess dosimetric and clinical effects of treatment planners having a priori knowledge of the maximum achievable dose sparing for organs at risk (OARs) for patients with oropharynx cancer receiving intensity modulated radiation therapy (RT). METHODS AND MATERIALS: We examined patients with oropharynx cancer who were treated in prospective clinical trials between February 2012 and April 2019 at our institution. A tool generating estimates of maximum achievable dose sparing for OARs (feasibility dose-volume histogram [FDVH]) was used clinically starting July 2016. Patients were divided into 2 cohorts: Before (ie, baseline) and after (ie, FDVH-guided) FDVH. Doses received by various OARs were compared with those estimated to be achievable per FDVH, and that difference was defined as the excess of feasible dose (EFD). Patient-reported outcome (PRO) questionnaires were completed at 3, 6, and 12 months after treatment. The baseline and FDVH-guided cohorts were compared in terms of EFD, plan quality metrics, and post-RT PRO assessments. RESULTS: A total of 139 patients were included in the analysis (60 in the baseline cohort, 79 in the FDVH-guided cohort). The FDVH-guided cohort had lower EFD to the contralateral parotid by 4.1 Gy, the ipsilateral parotid by 10.6 Gy, the larynx by 4.3 Gy, the oral cavity by 1.5 Gy, and the contralateral submandibular gland by 0.4 Gy. Plan quality metrics were similar between the cohorts. Less variation of EFD was seen in the FDVH-guided cohort for the parotid glands and contralateral submandibular gland (P < .05). The average post-RT PROs were better in the FVHD cohort versus baseline (particularly at the 6-month timepoint for dry mouth frequency, sticky saliva, meal enjoyment, severity of pain, and Eating Assessment Tool 10 composite [swallowing]; P < .05). CONCLUSIONS: Use of FDVH was associated with improved and less variable OAR sparing for clinically delivered plans. FDVH-guided patients had improved PROs compared with baseline with a variety of outcomes significantly improved at 6 months after treatment.

Prospective assessment of sparing the parotid ducts via MRI sialography for reducing patient reported xerostomia.

PURPOSE: To assess the impact of prospectively sparing the parotid ducts via MRI sialography on patient reported xerostomia for those receiving definitive radiotherapy (RT) for oropharyngeal squamous cell carcinoma. METHODS AND MATERIALS: Thirty-eight patients with oropharynx cancer to be treated with definitive RT underwent pre-treatment MRI sialograms to localize their parotid ducts. The parotid ducts were maximally spared during treatment planning. Patients reported symptoms (PRO-CTCAE and QLQ-H&N35) were collected at 6 and 12 months post-RT and compared to a historical cohort who underwent conventional parotid gland mean dose sparing. Regression models were generated using parotid and submandibular gland doses with and without incorporating the dose to the parotid ducts to determine the impact of parotid duct dose on patient reported xerostomia. RESULTS: At 6 months post-RT, 12/26 (46%) patients reported ≥moderate xerostomia when undergoing parotid ductal sparing compared to 43/61 (70%) in the historical cohort (p = 0.03). At 12 months post-RT, 8/22 (36%) patients reported ≥moderate xerostomia when undergoing parotid ductal sparing compared to 34/68(50%) in the historical cohort (p = 0.08). Using nested logistic regression models, the mean parotid duct dose was found to significantly relate to patient reported xerostomia severity at 6 months post-RT (p = 0.04) and trended towards statistical significance at 12 months post-RT (p = 0.09). At both 6 and 12 months post-RT, the addition of mean parotid duct dose significantly improved model fit (p < 0.05). CONCLUSIONS: MRI sialography guided parotid duct sparing appears to reduce the rates of patient-reported xerostomia. Further, logistic regression analysis found parotid duct dose to be significantly associated with patient reported xerostomia. A significant improvement in model fit was observed when adding mean parotid duct dose compared to models that only contain mean parotid gland dose and mean contralateral submandibular gland dose.

Feature Engineering for Interpretable Machine Learning for Quality Assurance in Radiation Oncology.

Chart checking is a time intensive process with high cognitive workload for physicists. Previous studies have partially automated and standardized chart checking, but limited studies implement data-driven approaches to reduce cognitive workload for quality assurance processes. This study aims to evaluate feature selection methods to improve the interpretability and transparency of machine learning models in predicting the degree of difficulty for a pretreatment physics chart check. We compare chi-square, mutual information, feature importance thresholding, and greedy feature selection for four different classifiers. Random forest has the highest performance with SMOTE oversampling using mutual information for feature selection (accuracy 84.0%, AUC 87.0%, precision 80.0%, recall 80.0%). This study demonstrates that feature selection methods can improve model interpretability and transparency.

Knowledge-based IMRT treatment planning for prostate cancer.

PURPOSE: To demonstrate the feasibility of using a knowledge base of prior treatment plans to generate new prostate intensity modulated radiation therapy (IMRT) plans. Each new case would be matched against others in the knowledge base. Once the best match is identified, that clinically approved plan is used to generate the new plan. METHODS: A database of 100 prostate IMRT treatment plans was assembled into an information-theoretic system. An algorithm based on mutual information was implemented to identify similar patient cases by matching 2D beam's eye view projections of contours. Ten randomly selected query cases were each matched with the most similar case from the database of prior clinically approved plans. Treatment parameters from the matched case were used to develop new treatment plans. A comparison of the differences in the dose-volume histograms between the new and the original treatment plans were analyzed. RESULTS: On average, the new knowledge-based plan is capable of achieving very comparable planning target volume coverage as the original plan, to within 2% as evaluated for D98, D95, and D1. Similarly, the dose to the rectum and dose to the bladder are also comparable to the original plan. For the rectum, the mean and standard deviation of the dose percentage differences for D20, D30, and D50 are 1.8% +/- 8.5%, -2.5% +/- 13.9%, and -13.9% +/- 23.6%, respectively. For the bladder, the mean and standard deviation of the dose percentage differences for D20, D30, and D50 are -5.9% +/- 10.8%, -12.2% +/- 14.6%, and -24.9% +/- 21.2%, respectively. A negative percentage difference indicates that the new plan has greater dose sparing as compared to the original plan. CONCLUSIONS: The authors demonstrate a knowledge-based approach of using prior clinically approved treatment plans to generate clinically acceptable treatment plans of high quality. This semiautomated approach has the potential to improve the efficiency of the treatment planning process while ensuring that high quality plans are developed.

Hyperpolarized 129Xe Magnetic Resonance Imaging for Functional Avoidance Treatment Planning in Thoracic Radiation Therapy: A Comparison of Ventilation- and Gas Exchange-Guided Treatment Plans.

PURPOSE: To present a methodology to use pulmonary gas exchange maps to guide functional avoidance treatment planning in radiation therapy (RT) and evaluate its efficacy compared with ventilation-guided treatment planning. METHODS AND MATERIALS: Before receiving conventional RT for non-small cell lung cancer, 11 patients underwent hyperpolarized 129Xe gas exchange magnetic resonance imaging to map the distribution of xenon in its gas phase (ventilation) and transiently bound to red blood cells in the alveolar capillaries (gas exchange). Both ventilation and gas exchange maps were independently used to guide development of new functional avoidance treatment plans for every patient, while adhering to institutional dose-volume constraints for normal tissues and target coverage. Furthermore, dose-volume histogram (DVH)-based reoptimizations of the clinical plan, with reductions in mean lung dose (MLD) equal to the functional avoidance plans, were created to serve as the control group. To evaluate each plan (regardless of type), gas exchange maps, representing end-to-end lung function, were used to calculate gas exchange-weighted MLD (fMLD), gas exchange-weighted volume receiving ≥20 Gy (fV20), and mean dose in the highest gas exchanging 33% and 50% volumes of lung (MLD-f33% and MLD-f50%). Using each clinically approved plan as a baseline, the reductions in functional metrics were compared for ventilation-optimization, gas exchange optimization, and DVH-based reoptimization. Statistical significance was determined using the Freidman test, with subsequent subdivision when indicated by P values less than .10 and post hoc testing with Wilcoxon signed rank tests to determine significant differences (P < .05). Toxicity modeling was performed using an established function-based model to estimate clinical significance of the results. RESULTS: Compared with DVH-based reoptimization of the clinically approved plans, gas exchange-guided functional avoidance planning more effectively reduced the gas exchange-weighted metrics fMLD (average ± SD, -78 ± 79 cGy, compared with -45 ± 34 cGy; P = .03), MLD-f33% (-135 ± 136 cGy, compared with -52 ± 47 cGy; P = .004), and MLD-f50% (-96 ± 95 cGy, compared with -47 ± 40 cGy; P = .01). Comparing the 2 functional planning types, Gas Exchange-Guided planning more effectively reduced MLD-f33% compared with ventilation-guided planning (-64 ± 95; P = .009). For some patients, Gas Exchange-Guided functional avoidance plans demonstrated clinically significant reductions in model-predicted toxicity, more so than the accompanying ventilation-guided plans and DVH-based reoptimizations. CONCLUSION: Gas Exchange-Guided planning effectively reduced dose to high gas exchanging regions of lung while maintaining clinically acceptable plan quality. In many patients, ventilation-guided planning incidentally reduced dose to higher gas exchange regions, to a lesser extent. This methodology enables future prospective trials to examine patient outcomes.

NTCP modeling and dose-volume correlations for acute xerostomia and dry eye after whole brain radiation.

BACKGROUND: Whole brain radiation (WBRT) may lead to acute xerostomia and dry eye from incidental parotid and lacrimal exposure, respectively. We performed a prospective observational study to assess the incidence/severity of this toxicity. We herein perform a secondary analysis relating parotid and lacrimal dosimetric parameters to normal tissue complication probability (NTCP) rates and associated models. METHODS: Patients received WBRT to 25-40 Gy in 10-20 fractions using 3D-conformal radiation therapy without prospective delineation of the parotids or lacrimals. Patients completed questionnaires at baseline and 1 month post-WBRT. Xerostomia was assessed using the University of Michigan xerostomia score (scored 0-100, toxicity defined as ≥ 20 pt increase) and xerostomia bother score (scored from 0 to 3, toxicity defined as ≥ 2 pt increase). Dry eye was assessed using the Subjective Evaluation of Symptom of Dryness (SESoD, scored from 0 to 4, toxicity defined as ≥ 2 pt increase). The clinical data were fitted by the Lyman-Kutcher-Burman (LKB) and Relative Seriality (RS) NTCP models. RESULTS: Of 55 evaluable patients, 19 (35%) had ≥ 20 point increase in xerostomia score, 11 (20%) had ≥ 2 point increase in xerostomia bother score, and 13 (24%) had ≥ 2 point increase in SESoD score. For xerostomia, parotid V10Gy-V20Gy correlated best with toxicity, with AUC 0.68 for xerostomia score and 0.69-0.71 for bother score. The values for the D50, m and n parameters of the LKB model were 22.3 Gy, 0.84 and 1.0 for xerostomia score and 28.4 Gy, 0.55 and 1.0 for bother score, respectively. The corresponding values for the D50, γ and s parameters of the RS model were 23.5 Gy, 0.28 and 0.0001 for xerostomia score and 32.0 Gy, 0.45 and 0.0001 for bother score, respectively. For dry eye, lacrimal V10Gy-V15Gy were found to correlate best with toxicity, with AUC values from 0.67 to 0.68. The parameter values of the LKB model were 53.5 Gy, 0.74 and 1.0, whereas of the RS model were 54.0 Gy, 0.37 and 0.0001, respectively. CONCLUSIONS: Xerostomia was most associated with parotid V10Gy-V20Gy, and dry eye with lacrimal V10Gy-V15Gy. NTCP models were successfully created for both toxicities and may help clinicians refine dosimetric goals and assess levels of risk in patients receiving palliative WBRT.

The utility of non-axial treatment beam orientations for lower lobe lung cancers.

Traditional treatment beams for non-small-cell lung cancer are limited to the axial plane. For many tumor geometries, non-axial orientations appear to reduce the dose to normal tissues (e.g. heart, liver). We hypothesize that non-axial beams provide a significant reduction in incidental irradiation of the heart and liver, while maintaining adequate target coverage. CT scans of twenty-four patients with lower lobe lung cancers were studied. For each patient, an opposed oblique axial beam pair and a competing non-axial opposed oblique pair were generated, both off-cord. The competing plans delivered comparable doses/margins to the GTV. DVHs and integral doses were computed for all structures of interest for the two competing plans. The integral dose was compared for axial and non-axial beams for each contoured organ using a paired t-test. Dose to the heart was significantly lower for the non-axial plans ( p = .0001). For 20/24 patients, the integral heart dose was reduced by using non-axial beams. In those patients with tumors located in the inferior right lower lobe, a lower dose to the liver was achieved when non-axial beams were used. There were no meaningful differences in dose to the GTV, lungs, or skin between axial and non-axial beams. Non-axial beams can reduce the dose to the heart and liver in patients with lower lobe lung cancers. Non-axial beams may be clinically beneficial in these patients and should be considered as an option during planning.

A method to dynamically balance intensity modulated radiotherapy dose between organs-at-risk.

The IMRT treatment planning process typically follows a path that is based on the manner in which the planner interactively adjusts the target and organ-at-risk (OAR) constraints and priorities. The time-intensive nature of this process restricts the planner from fully understanding the dose tradeoff between structures, making it unlikely that the resulting plan fully exploits the extent to which dose can be redistributed between anatomical structures. Multiobjective Pareto optimization has been used in the past to enable the planner to more thoroughly explore alternatives in dose trade-off by combining pre-generated Pareto optimal solutions in real time, thereby potentially tailoring a plan more exactly to requirements. However, generating the Pareto optimal solutions can be nonintuitive and computationally time intensive. The author presents an intuitive and fast non-Pareto approach for generating optimization sequences (prior to planning), which can then be rapidly combined by the planner in real time to yield a satisfactory plan. Each optimization sequence incrementally reduces dose to one OAR at a time, starting from the optimization solution where dose to all OARs are reduced with equal priority, until user-specified target coverage limits are violated. The sequences are computationally efficient to generate, since the optimization at each position along a sequence is initiated from the end result of the previous position in the sequence. The pre-generated optimization sequences require no user interaction. In real time, a planner can more or less instantaneously visualize a treatment plan by combining the dose distributions corresponding to user-selected positions along each of the optimization sequences (target coverage is intrinsically maintained in the combination). Interactively varying the selected positions along each of the sequences enables the planner to rapidly understand the nature of dose trade-off between structures and, thereby, arrive at a suitable plan in a short time. This methodology is demonstrated on a prostate cancer case and olfactory neuroblastoma case.

Integral dose conservation in radiotherapy.

Treatment planners frequently modify beam arrangements and use IMRT to improve target dose coverage while satisfying dose constraints on normal tissues. The authors herein analyze the limitations of these strategies and quantitatively assess the extent to which dose can be redistributed within the patient volume. Specifically, the authors hypothesize that (1) the normalized integral dose is constant across concentric shells of normal tissue surrounding the target (normalized to the average integral shell dose), (2) the normalized integral shell dose is constant across plans with different numbers and orientations of beams, and (3) the normalized integral shell dose is constant across plans when reducing the dose to a critical structure. Using the images of seven patients previously irradiated for cancer of brain or prostate cancer and one idealized scenario, competing three-dimensional conformal and IMRT plans were generated using different beam configurations. Within a given plan and for competing plans with a constant mean target dose, the normalized integral doses within concentric "shells" of surrounding normal tissue were quantitatively compared. Within each patient, the normalized integral dose to shells of normal tissue surrounding the target was relatively constant (1). Similarly, for each clinical scenario, the normalized integral dose for a given shell was also relatively constant regardless of the number and orientation of beams (2) or degree of sparing of a critical structure (3). 3D and IMRT planning tools can redistribute, rather than eliminate dose to the surrounding normal tissues (intuitively known by planners). More specifically, dose cannot be moved between shells surrounding the target but only within a shell. This implies that there are limitations in the extent to which a critical structure can be spared based on the location and geometry of the critical structure relative to the target.