Incorporating variable RBE in IMPT optimization for ependymoma.

PURPOSE: To study the dosimetric impact of incorporating variable relative biological effectiveness (RBE) of protons in optimizing intensity-modulated proton therapy (IMPT) treatment plans and to compare it with conventional constant RBE optimization and linear energy transfer (LET)-based optimization. METHODS: This study included 10 pediatric ependymoma patients with challenging anatomical features for treatment planning. Four plans were generated for each patient according to different optimization strategies: (1) constant RBE optimization (ConstRBEopt) considering standard-of-care dose requirements; (2) LET optimization (LETopt) using a composite cost function simultaneously optimizing dose-averaged LET (LETd ) and dose; (3) variable RBE optimization (VarRBEopt) using a recent phenomenological RBE model developed by McNamara et al.; and (4) hybrid RBE optimization (hRBEopt) assuming constant RBE for the target and variable RBE for organs at risk. By normalizing each plan to obtain the same target coverage in either constant or variable RBE, we compared dose, LETd , LET-weighted dose, and equivalent uniform dose between the different optimization approaches. RESULTS: We found that the LETopt plans consistently achieved increased LET in tumor targets and similar or decreased LET in critical organs compared to other plans. On average, the VarRBEopt plans achieved lower mean and maximum doses with both constant and variable RBE in the brainstem and spinal cord for all 10 patients. To compensate for the underdosing of targets with 1.1 RBE for the VarRBEopt plans, the hRBEopt plans achieved higher physical dose in targets and reduced mean and especially maximum variable RBE doses compared to the ConstRBEopt and LETopt plans. CONCLUSION: We demonstrated the feasibility of directly incorporating variable RBE models in IMPT optimization. A hybrid RBE optimization strategy showed potential for clinical implementation by maintaining all current dose limits and reducing the incidence of high RBE in critical normal tissues in ependymoma patients.

Impact of RBE variations on risk estimates of temporal lobe necrosis in patients treated with intensity-modulated proton therapy for head and neck cancer.

BACKGROUND: Temporal lobe necrosis (TLN) is a potential late effect after radiotherapy for skull base head and neck cancer (HNC). Several photon-derived dose constraints and normal tissue complication probability (NTCP) models have been proposed, however variation in relative biological effectiveness (RBE) may challenge the applicability of these dose constraints and models in proton therapy. The purpose of this study was therefore to investigate the influence of RBE variations on risk estimates of TLN after Intensity-Modulated Proton Therapy for HNC. MATERIAL AND METHODS: Seventy-five temporal lobes from 45 previously treated patients were included in the analysis. Sixteen temporal lobes had radiation associated Magnetic Resonance image changes (TLIC) suspected to be early signs of TLN. Fixed (RWDFix) and variable RBE-weighed doses (RWDVar) were calculated using RBE = 1.1 and two RBE models, respectively. RWDFix and RWDVar for temporal lobes were compared using Friedman's test. Based on RWDFix, six NTCP models were fitted and internally validated through bootstrapping. Estimated probabilities from RWDFix and RWDVar were compared using paired Wilcoxon test. Seven dose constraints were evaluated separately for RWDFix and RWDVar by calculating the observed proportion of TLIC in temporal lobes meeting the specific dose constraints. RESULTS: RWDVar were significantly higher than RWDFix (p < 0.01). NTCP model performance was good (AUC:0.79-0.84). The median difference in estimated probability between RWDFix and RWDVar ranged between 5.3% and 20.0% points (p < 0.01), with V60GyRBE and DMax at the smallest and largest differences, respectively. The proportion of TLIC was higher for RWDFix (4.0%-13.1%) versus RWDVar (1.3%-5.3%). For V65GyRBE ≤ 0.03 cc the proportion of TLIC was less than 5% for both RWDFix and RWDVar. CONCLUSION: NTCP estimates were significantly influenced by RBE variations. Dmax as model predictor resulted in the largest deviations in risk estimates between RWDFix and RWDVar. V65GyRBE ≤ 0.03 cc was the most consistent dose constraint for RWDFix and RWDVar.

Monte Carlo evaluation of target dose coverage in lung stereotactic body radiation therapy with flattening filter-free beams.

Aim: Previous studies showed that replacing conventional flattened beams (FF) with flattening filter-free (FFF) beams improves the therapeutic ratio in lung stereotactic body radiation therapy (SBRT), but these findings could have been impacted by dose calculation uncertainties caused by the heterogeneity of the thoracic anatomy and by respiratory motion, which were particularly high for target coverage. In this study, we minimized such uncertainties by calculating doses using high-spatial-resolution Monte Carlo and four-dimensional computed tomography (4DCT) images. We aimed to evaluate more reliably the benefits of using FFF beams for lung SBRT. Materials and methods: For a cohort of 15 patients with early stage lung cancer that we investigated in a previous treatment planning study, we recalculated dose distributions with Monte Carlo using 4DCT images. This included fifteen FF and fifteen FFF treatment plans. Results: Compared to Monte Carlo, the treatment planning system (TPS) over-predicted doses in low-dose regions of the planning target volume. For most patients, replacing FF beams with FFF beams improved target coverage, tumor control, and uncomplicated tumor control probabilities. Conclusions: Monte Carlo tends to reveal deficiencies in target coverage compared to coverage predicted by the TPS. Our data support previously reported benefits of using FFF beams for lung SBRT.

A comparison of two methodologies for radiotherapy treatment plan optimization and QA for clinical trials.

BACKGROUND AND PURPOSE: The efficacy of clinical trials and the outcome of patient treatment are dependent on the quality assurance (QA) of radiation therapy (RT) plans. There are two widely utilized approaches that include plan optimization guidance created based on patient-specific anatomy. This study examined these two techniques for dose-volume histogram predictions, RT plan optimizations, and prospective QA processes, namely the knowledge-based planning (KBP) technique and another first principle (FP) technique. METHODS: This analysis included 60, 44, and 10 RT plans from three Radiation Therapy Oncology Group (RTOG) multi-institutional trials: RTOG 0631 (Spine SRS), RTOG 1308 (NSCLC), and RTOG 0522 (H&N), respectively. Both approaches were compared in terms of dose prediction and plan optimization. The dose predictions were also compared to the original plan submitted to the trials for the QA procedure. RESULTS: For the RTOG 0631 (Spine SRS) and RTOG 0522 (H&N) plans, the dose predictions from both techniques have correlation coefficients of >0.9. The RT plans that were re-optimized based on the predictions from both techniques showed similar quality, with no statistically significant differences in target coverage or organ-at-risk sparing. The predictions of mean lung and heart doses from both methods for RTOG1308 patients, on the other hand, have a discrepancy of up to 14 Gy. CONCLUSIONS: Both methods are valuable tools for optimization guidance of RT plans for Spine SRS and Head and Neck cases, as well as for QA purposes. On the other hand, the findings suggest that KBP may be more feasible in the case of inoperable lung cancer patients who are treated with IMRT plans that have spatially unevenly distributed beam angles.

Geometric and dosimetric accuracy of deformable image registration between average-intensity images for 4DCT-based adaptive radiotherapy for non-small cell lung cancer.

PURPOSE: Re-planning for four-dimensional computed tomography (4DCT)-based lung adaptive radiotherapy commonly requires deformable dose mapping between the planning average-intensity image (AVG) and the newly acquired AVG. However, such AVG-AVG deformable image registration (DIR) lacks accuracy assessment. The current work quantified and compared geometric accuracies of AVG-AVG DIR and corresponding phase-phase DIRs, and subsequently investigated the clinical impact of such AVG-AVG DIR on deformable dose mapping. METHODS AND MATERIALS: Hybrid intensity-based AVG-AVG and phase-phase DIRs were performed between the planning and mid-treatment 4DCTs of 28 non-small cell lung cancer patients. An automated landmark identification algorithm detected vessel bifurcation pairs in both lungs. Target registration error (TRE) of these landmark pairs was calculated for both DIR types. The correlation between TRE and respiratory-induced landmark motion in the planning 4DCT was analyzed. Global and local dose metrics were used to assess the clinical implications of AVG-AVG deformable dose mapping with both DIR types. RESULTS: TRE of AVG-AVG and phase-phase DIRs averaged 3.2 ± 1.0 and 2.6 ± 0.8 mm respectively (p < 0.001). Using AVG-AVG DIR, TREs for landmarks with <10 mm motion averaged 2.9 ± 2.0 mm, compared to 3.1 ± 1.9 mm for the remaining landmarks (p < 0.01). Comparatively, no significant difference was demonstrated for phase-phase DIRs. Dosimetrically, no significant difference in global dose metrics was observed between doses mapped with AVG-AVG DIR and the phase-phase DIR, but a positive linear relationship existed (p = 0.04) between the TRE of AVG-AVG DIR and local dose difference. CONCLUSIONS: When the region of interest experiences <10 mm respiratory-induced motion, AVG-AVG DIR may provide sufficient geometric accuracy; conversely, extra attention is warranted, and phase-phase DIR is recommended. Dosimetrically, the differences in geometric accuracy between AVG-AVG and phase-phase DIRs did not impact global lung-based metrics. However, as more localized dose metrics are needed for toxicity assessment, phase-phase DIR may be required as its lower mean TRE improved voxel-based dosimetry.

Using FFF Beams to Improve the Therapeutic Ratio of Lung SBRT.

Aim: To investigate the extent to which lung stereotactic body radiotherapy (SBRT) treatment plans can be improved by replacing conventional flattening filter (FF) beams with flattening filter-free (FFF) beams. Material and Methods: We selected 15 patients who had received SBRT with conventional 6-MV photon beams for early-stage lung cancer. We imported the patients' treatment plans into the Eclipse 13.6 treatment planning system, in which we configured the AAA dose calculation model using representative beam data for a TrueBeam accelerator operated in 6-MV FFF mode. We then created new treatment plans by replacing the conventional FF beams in the original plans with FFF beams. Results: The FFF plans had better target coverage than the original FF plans did. For the planning target volume, FFF plans significantly improved the D98, D95, D90, homogeneity index, and uncomplicated tumor control probability. In most cases, the doses to organs at risk were lower in FFF plans. FFF plans significantly reduced the mean lung dose, V10, V20, V30, and normal tissue complication probability for the total lung and improved the dosimetric indices for the ipsilateral lung. For most patients, FFF beams achieved lower maximum doses to the esophagus, heart, and the spinal cord; and a lower chest wall V30. Findings: Compared with FF beams, FFF beams achieved lower doses to organs at risk, especially the lung, without compromising tumor coverage; in fact, FFF beams improved coverage in most cases. Thus, replacing FF beams with FFF beams can achieve a better therapeutic ratio.

Impact of Intra-Fractional Motion on Dose Distributions in Lung IMRT.

AIM: To investigate the impact of intra-fractional motion on dose distribution in patients treated with intensity-modulated radiation therapy (IMRT) for lung cancer. MATERIALS AND METHODS: Twenty patients who had undergone IMRT for non-small cell lung cancer were selected for this retrospective study. For each patient, a four-dimensional computed tomography (CT) image set was acquired and clinical treatment plans were developed using the average CT. Dose distributions were then re-calculated for each of the 10 phases of respiratory cycle and combined using deformable image registration to produce cumulative dose distributions that were compared with the clinical treatment plans. RESULTS: Intra-fractional motion reduced planning target volume (PTV) coverage in all patients. The median reduction of PTV volume covered by the prescription isodose was 3.4%; D98 was reduced by 3.1 Gy. Changes in the mean lung dose were within ±0.7 Gy. V20 for the lung increased in most patients; the median increase was 1.6%. The dose to the spinal cord was unaffected by intra-fractional motion. The dose to the heart was slightly reduced in most patients. The median reduction in the mean heart dose was 0.22 Gy, and V30 was reduced by 2.5%.The maximum dose to the esophagus was also reduced in most patients, by 0.74 Gy, whereas V50 did not change significantly. The median number of points in which dose differences exceeded 3%/3 mm was 6.2%. FINDINGS: Intra-fractional anatomical changes reduce PTV coverage compared to the coverage predicted by clinical treatment planning systems that use the average CT for dose calculation. Doses to organs at risk were mostly over-predicted.

Proton therapy for adults with mediastinal lymphomas: the International Lymphoma Radiation Oncology Group guidelines.

Among adult lymphoma survivors, radiation treatment techniques that increase the excess radiation dose to organs at risk (OARs) put patients at risk for increased side effects, especially late toxicities. Minimizing radiation to OARs in adults patients with Hodgkin and non-Hodgkin lymphomas involving the mediastinum is the deciding factor for the choice of treatment modality. Proton therapy may help to reduce the radiation dose to the OARs and reduce toxicities, especially the risks for cardiac morbidity and second cancers. Because proton therapy may have some disadvantages, identifying the patients and the circumstances that may benefit the most from proton therapy is important. We present modern guidelines to identify adult lymphoma patients who may derive the greatest benefit from proton therapy, along with an analysis of the advantages and disadvantages of proton treatment.

Impact of respiratory motion on variable relative biological effectiveness in 4D-dose distributions of proton therapy.

BACKGROUND: Organ motion during radiation therapy with scanned protons leads to deviations between the planned and the delivered physical dose. Using a constant relative biological effectiveness (RBE) of 1.1 linearly maps these deviations into RBE-weighted dose. However, a constant value cannot account for potential nonlinear variations in RBE suggested by variable RBE models. Here, we study the impact of motion on recalculations of RBE-weighted dose distributions using a phenomenological variable RBE model. MATERIAL AND METHODS: 4D-dose calculation including variable RBE was implemented in the open source treatment planning toolkit matRad. Four scenarios were compared for one field and two field proton treatments for a liver cancer patient assuming (α∕ß)x = 2 Gy and (α∕ß)x = 10 Gy: (A) the optimized static dose distribution with constant RBE, (B) a static recalculation with variable RBE, (C) a 4D-dose recalculation with constant RBE and (D) a 4D-dose recalculation with variable RBE. For (B) and (D), the variable RBE was calculated by the model proposed by McNamara. For (C), the physical dose was accumulated with direct dose mapping; for (D), dose-weighted radio-sensitivity parameters of the linear quadratic model were accumulated to model synergistic irradiation effects on RBE. RESULTS: Dose recalculation with variable RBE led to an elevated biological dose at the end of the proton field, while 4D-dose recalculation exhibited random deviations everywhere in the radiation field depending on the interplay of beam delivery and organ motion. For a single beam treatment assuming (α∕ß)x = 2 Gy, D95% was 1.98 Gy (RBE) (A), 2.15 Gy (RBE) (B), 1.81 Gy (RBE) (C) and 1.98 Gy (RBE) (D). The homogeneity index was 1.04 (A), 1.08 (B), 1.23 (C) and 1.25 (D). CONCLUSION: For the studied liver case, intrafractional motion did not reduce the modulation of the RBE-weighted dose postulated by variable RBE models for proton treatments.

Radiobiological issues in proton therapy.

BACKGROUND: The relative biological effectiveness (RBE) for particle therapy is a complex function of particle type, radiation dose, linear energy transfer (LET), cell type, endpoint, etc. In the clinical practice of proton therapy, the RBE is assumed to have a fixed value of 1.1. This assumption, along with the effects of physical uncertainties, may mean that the biologically effective dose distributions received by the patient may be significantly different from what is seen on treatment plans. This may contribute to unforeseen toxicities and/or failure to control the disease. Variability of Proton RBE: It has been shown experimentally that proton RBE varies significantly along the beam path, especially near the end of the particle range. While there is now an increasing acceptance that proton RBE is variable, there is an ongoing debate about whether to change the current clinical practice. Clinical Evidence: A rationale against the change is the uncertainty in the models of variable RBE. Secondly, so far there is no clear clinical evidence of the harm of assuming proton RBE to be 1.1. It is conceivable, however, that the evidence is masked partially by physical uncertainties. It is, therefore, plausible that reduction in uncertainties and their incorporation in the estimation of dose actually delivered may isolate and reveal the variability of RBE in clinical practice. Nevertheless, clinical evidence of RBE variability is slowly emerging as more patients are treated with protons and their response data are analyzed. Modelling and Incorporation of RBE in the Optimization of Proton Therapy: The improvement in the knowledge of RBE could lead to better understanding of outcomes of proton therapy and in the improvement of models to predict RBE. Prospectively, the incorporation of such models in the optimization of intensity-modulated proton therapy could lead to improvements in the therapeutic ratio of proton therapy.

Comparison of linear and nonlinear programming approaches for "worst case dose" and "minmax" robust optimization of intensity-modulated proton therapy dose distributions.

Robust optimization of intensity-modulated proton therapy (IMPT) takes uncertainties into account during spot weight optimization and leads to dose distributions that are resilient to uncertainties. Previous studies demonstrated benefits of linear programming (LP) for IMPT in terms of delivery efficiency by considerably reducing the number of spots required for the same quality of plans. However, a reduction in the number of spots may lead to loss of robustness. The purpose of this study was to evaluate and compare the performance in terms of plan quality and robustness of two robust optimization approaches using LP and nonlinear programming (NLP) models. The so-called "worst case dose" and "minmax" robust optimization approaches and conventional planning target volume (PTV)-based optimization approach were applied to designing IMPT plans for five patients: two with prostate cancer, one with skull-based cancer, and two with head and neck cancer. For each approach, both LP and NLP models were used. Thus, for each case, six sets of IMPT plans were generated and assessed: LP-PTV-based, NLP-PTV-based, LP-worst case dose, NLP-worst case dose, LP-minmax, and NLP-minmax. The four robust optimization methods behaved differently from patient to patient, and no method emerged as superior to the others in terms of nominal plan quality and robustness against uncertainties. The plans generated using LP-based robust optimization were more robust regarding patient setup and range uncertainties than were those generated using NLP-based robust optimization for the prostate cancer patients. However, the robustness of plans generated using NLP-based methods was superior for the skull-based and head and neck cancer patients. Overall, LP-based methods were suitable for the less challenging cancer cases in which all uncertainty scenarios were able to satisfy tight dose constraints, while NLP performed better in more difficult cases in which most uncertainty scenarios were hard to meet tight dose limits. For robust optimization, the worst case dose approach was less sensitive to uncertainties than was the minmax approach for the prostate and skull-based cancer patients, whereas the minmax approach was superior for the head and neck cancer patients. The robustness of the IMPT plans was remarkably better after robust optimization than after PTV-based optimization, and the NLP-PTV-based optimization outperformed the LP-PTV-based optimization regarding robustness of clinical target volume coverage. In addition, plans generated using LP-based methods had notably fewer scanning spots than did those generated using NLP-based methods.

Perturbation of water-equivalent thickness as a surrogate for respiratory motion in proton therapy.

Respiratory motion is traditionally assessed using tumor motion magnitude. In proton therapy, respiratory motion causes density variations along the beam path that result in uncertainties of proton range. This work has investigated the use of water-equivalent thickness (WET) to quantitatively assess the effects of respiratory motion on calculated dose in passively scattered proton therapy (PSPT). A cohort of 29 locally advanced non-small cell lung cancer patients treated with 87 PSPT treatment fields were selected for analysis. The variation in WET (ΔWET) along each field was calculated between exhale and inhale phases of the simulation four-dimensional computed tomography. The change in calculated dose (ΔDose) between full-inhale and full-exhale phase was quantified for each field using dose differences, 3D gamma analysis, and differential area under the curve (ΔAUC) analysis. Pearson correlation coefficients were calculated between ΔDose and ΔWET. Three PSPT plans were redesigned using field angles to minimize variations in ΔWET and compared to the original plans. The median ΔWET over 87 treatment fields ranged from 1-9 mm, while the ΔWET 95th percentile value ranged up to 42 mm. The ΔWET was significantly correlated (p < 0.001) to the ΔDose for all metrics analyzed. The patient plans that were redesigned using ΔWET analysis to select field angles were more robust to the effects of respiratory motion, as ΔAUC values were reduced by more than 60% in all three cases. The tumor motion magnitude alone does not capture the potential dosimetric error due to respiratory motion because the proton range is sensitive to the motion of all patient anatomy. The use of ΔWET has been demonstrated to identify situations where respiratory motion can impact the calculated dose. Angular analysis of ΔWET may be capable of designing radiotherapy plans that are more robust to the effects of respiratory motion.

Is there an impact of heart exposure on the incidence of radiation pneumonitis? Analysis of data from a large clinical cohort.

BACKGROUND: The goal of the present study was to determine, in a large clinical cohort, whether incidental radiation exposure to the heart during definitive radiotherapy of inoperable non-small cell lung cancer (NSCLC) detectably increased the risk of radiation pneumonitis (RP) beyond that resulting from radiation exposure to lung. MATERIAL AND METHODS: Data were analyzed from all patients who received definitive three-dimensional (3D) concurrent radiotherapy or intensity-modulated radiotherapy for the treatment of NSCLC over a 10-year period at our institution, except those who had previous lung cancer or for whom radiation treatment plans were unavailable for calculation of heart and lung dose-volume histograms (DVHs). Parameters computed from heart and lung DVHs included mean lung dose (MLD), effective lung dose computed using volume parameter n = 0.5 (Deff), mean heart dose (MHD), percentage of heart receiving > 65 Gy (V65), and minimum dose to the hottest 10% of heart (D10). Univariate and multivariate normal-tissue complication probability (NTCP) models were used to analyze incidence of Grade ≥ 2 or Grade ≥ 3 RP as a function of these and other parameters. RESULTS: The study cohort included 629 patients, with crude rates of Grade ≥ 2 RP and Grade ≥ 3 RP of N = 263 (42%) and N = 124 (20%), respectively. Univariate NTCP models based on dosimetric lung parameters (MLD and Deff) fit the data better than models based on univariate heart parameters (heart D10, heart V65 or MHD). In multivariate modeling, incorporation of heart parameters did not significantly improve the fit of RP risk models based on lung parameters alone (p > 0.38 in each case). CONCLUSIONS: In this large clinical cohort, there was no evidence that incidental heart exposure during radiotherapy of NSCLC had a detectable impact on the occurrence of moderate or severe RP.

Radiation-Induced Lymphopenia is a Causal Mediator of Survival After Chemoradiation Therapy for Esophagus Cancer.

Purpose: Radiation-induced lymphopenia (RIL) is common during chemoradiation therapy. Severe lymphopenia is associated with reduced survival. Proton beam therapy (PBT), with its substantially more compact dose distributions, spares circulating lymphocytes and immune organs at risk to a greater extent than photon therapy. Recent studies comparing PBT to photon radiation therapy, specifically intensity-modulated radiation therapy (IMRT) for esophageal cancer (EC), showed that the incidence of grade 4 RIL (G4RIL) is significantly reduced among patients receiving PBT for EC. However, whether the extent of this reduction has a direct causative link with improved survival is unknown. This study applies causal mediation analysis to answer this question. Methods and Materials: We retrospectively assessed 734 patients treated with concurrent chemoradiation therapy for biopsy-proven EC from 2004 to 2017. To address the potential for bias in the choice of radiation modality, propensity score analysis was used to evaluate and reduce imbalances between the PBT and IMRT cohorts. Causal mediation analysis was applied to decompose the total effect of radiation modality on overall survival (OS) into indirect (mediated through G4RIL) and direct effects. Results: We found that PBT was associated with a significantly lower incidence of G4RIL and prolonged OS compared with IMRT (odds ratio, 0.41; 95% CI, 0.28-0.60; P < .001). In the propensity-matched cohort of 506 patients (253 PBT, 253 IMRT), G4RIL risk reduction with PBT versus IMRT translated into a 5% reduction in the relative rate of death (P = .032). Mediation of G4RIL explained ∼14.5% of the difference in OS. Conclusions: G4RIL was found to mediate survival; however, a statistically significant direct effect of PBT on survival was not observed. In other words, the statistical significance of survival benefit from protons over photons in this EC cohort was lost in the absence of G4RIL risk reduction.

Personalized Composite Dosimetric Score-Based Machine Learning Model of Severe Radiation-Induced Lymphopenia Among Patients With Esophageal Cancer.

PURPOSE: Radiation-induced lymphopenia (RIL) is common among patients undergoing radiation therapy (RT)' Severe RIL has been linked to adverse outcomes. The severity and risk of RIL can be predicted from baseline clinical characteristics and dosimetric parameters. However, dosimetric parameters, e.g. dose-volume (DV) indices, are highly correlated with one another and are only weakly associated with RIL. Here we introduce the novel concept of "composite dosimetric score" (CDS) as the index that condenses the dose distribution in immune tissues of interest to study the dosimetric dependence of RIL. We derived an improved multivariate classification scheme for risk of grade 4 RIL (G4RIL), based on this novel RT dosimetric feature, for patients receiving chemo RT for esophageal cancer. METHODS AND MATERIALS: DV indices were extracted for 734 patients who received chemo RT for biopsy-proven esophageal cancer. Nonnegative matrix factorization was used to project the DV indices of lung, heart, and spleen into a single CDS; XGBoost was employed to explore significant interactions among predictors; and logistic regression was applied to combine the resultant CDS with baseline clinical factors and interaction terms to facilitate individualized prediction of immunotoxicity. Five-fold cross-validation was applied to evaluate the model performance. RESULTS: The CDS for selected immune organs at risk (ie, heart, lung, and spleen) (OR 1.791; 95 CI [1.350, 2.377]) was a statistically significant risk determinant for G4RIL. Pearson correlation coefficients for CDS versus G4RIL risk for individual immune organs at risk were greater than any single DV indicx. Personalized prediction of G4RIL based on CDS and 4 clinical risk factors yielded an area under the curve value of 0.78. Interaction between age and CDS revealed that G4RIL risk increased more sharply with increasing CDS for patients aged ≥65 years. CONCLUSIONS: Risk of immunotoxicity for patients undergoing chemo RT for esophageal cancer can be predicted by CDS. The CDS concept can be extended to immunotoxicity in other cancer types and in dose-response models currently based on DV indices. Personalized treatment planning should leverage composite dosimetric scoring methods rather than using individual or subsets of DV indices.

Severe Lymphopenia During Chemoradiation Therapy for Esophageal Cancer: Comprehensive Analysis of Randomized Phase 2B Trial of Proton Beam Therapy Versus Intensity Modulated Radiation Therapy.

PURPOSE: Lymphocytes play an important role in antitumor immunity; however, they are also especially vulnerable to depletion during chemoradiation therapy (CRT). The purpose of this study was to compare the incidence of grade 4 lymphopenia (G4L) between proton beam therapy (PBT) and intensity modulated photon radiation therapy (IMRT) in patients with esophageal cancer treated with CRT in a completed randomized trial and to ascertain patient heterogeneity to G4L risk based on treatment and established prognostic factors. METHODS AND MATERIALS: Between April 2012 and March 2019, a single-institution, open-label, nonblinded, phase 2 randomized trial (NCT01512589) was conducted at the University of Texas MD Anderson Cancer Center. Patients were randomly assigned to IMRT or PBT, either definitively or preoperatively. This secondary analysis of the randomized trial was G4L during concurrent CRT according to Common Terminology Criteria for Adverse Events version 5.0. RESULTS: Among 105 patients evaluable for analysis, 44 patients (42%) experienced G4L at a median of 28 days after the start date of concurrent CRT. Induction chemotherapy (P = .003), baseline absolute lymphocyte count (P < .001), radiation therapy modality (P = .002), and planning treatment volume (P = .033) were found to be significantly associated with G4L. Multivariate classification analysis partitioned patients into 5 subgroups for whom the incidence of G4L was observed in 0%, 14%, 35%, 70%, and 100% of patients. The benefit of PBT over IMRT was most pronounced in patients with an intermediate baseline absolute lymphocyte count and large planning treatment volume (P = .011). CONCLUSIONS: This is the first prospective evidence that limiting dose scatter by PBT significantly reduced the incidence of G4L, especially in the intermediate-risk patients. The implication of this immune-sparing effect of PBT, especially in the context of standard adjuvant immunotherapy, needs further examination in the current phase 3 randomized trials.

Comparing 90-Day Postoperative Mortality After Neoadjuvant Proton-Based Versus Photon-Based Chemoradiotherapy for Esophageal Cancer.

Purpose: Evidence suggests that proton-beam therapy (PBT) results in less toxicity and postoperative complications compared to photon-based radiotherapy in patients who receive chemoradiotherapy followed by esophagectomy for cancer. Ninety-day mortality (90DM) is an important measure of the postoperative (nononcologic) outcome as proxy of quality-of-care. We hypothesize that PBT could reduce 90DM compared to photon-based radiotherapy. Materials and Methods: From a single-center retrospective database patients treated with chemoradiotherapy before esophagectomy for cancer were selected (1998-2022). Univariable logistic regression was used to study the association of radiotherapy modality with 90DM. Three separate methods were applied to adjust for confounding bias, including multivariable logistic regression, propensity score matching, and inverse probability of treatment weighting. Stratified analysis for the age threshold that maximized the difference in 90DM (ie, ≥67 vs <67 years) was performed. Results: A total of 894 eligible patients were included and 90DM was 5/202 (2.5%) in the PBT versus 29/692 (4.2%) in the photon-based radiotherapy group (P = .262). After adjustment for age and tumor location, PBT versus photon-based radiotherapy was not significantly associated with 90DM (P = .491). The 90DM was not significantly different for PBT versus photon-based radiotherapy in the propensity score matching (P = .379) and inverse probability of treatment weighting cohort (P = .426). The stratified analysis revealed that in patients aged ≥67 years, PBT was associated with decreased 90DM (1.3% vs 8.8%; P = .026). Higher age significantly increased 90DM risk within the photon-based radiotherapy (8.8% vs 2.7%; P = .001), but not within the PBT group (1.3% vs 3.2%; P = .651). Conclusion: No statistically significant difference was observed in postoperative 90DM after esophagectomy for cancer between PBT and photon-based neoadjuvant chemoradiotherapy. However, among older patients a signal was observed that PBT may reduce 90DM risk.

Clinical Translation of a Deep Learning Model of Radiation-Induced Lymphopenia for Esophageal Cancer.

Purpose: Radiation-induced lymphopenia is a common immune toxicity that adversely impacts treatment outcomes. We report here our approach to translate a deep-learning (DL) model developed to predict severe lymphopenia risk among esophageal cancer into a strategy for incorporating the immune system as an organ-at-risk (iOAR) to mitigate the risk. Materials and Methods: We conducted "virtual clinical trials" utilizing retrospective data for 10 intensity-modulated radiation therapy (IMRT) and 10 passively-scattered proton therapy (PSPT) esophageal cancer patients. For each patient, additional treatment plans of the modality other than the original were created employing standard-of-care (SOC) dose constraints. Predicted values of absolute lymphocyte count (ALC) nadir for all plans were estimated using a previously-developed DL model. The model also yielded the relative magnitudes of contributions of iOARs dosimetric factors to ALC nadir, which were used to compute iOARs dose-volume constraints, which were incorporated into optimization criteria to produce "IMRT-enhanced" and "intensity-modulated proton therapy (IMPT)-enhanced" plans. Results: Model-predicted ALC nadir for the original IMRT (IMRT-SOC) and PSPT plans agreed well with actual values. IMPT-SOC showed greater immune sparing vs IMRT and PSPT. The average mean body doses were 13.10 Gy vs 7.62 Gy for IMRT-SOC vs IMPT-SOC for patients treated with IMRT-SOC; and 8.08 Gy vs 6.68 Gy for PSPT vs IMPT-SOC for patients treated with PSPT. For IMRT patients, the average predicted ALC nadir of IMRT-SOC, IMRT-enhanced, IMPT-SOC, and IMPT-enhanced was 281, 327, 351, and 392 cells/µL, respectively. For PSPT patients, the average predicted ALC nadir of PSPT, IMPT-SOC, and IMPT-enhanced was 258, 316, and 350 cells/µL, respectively. Enhanced plans achieved higher predicted ALC nadir, with an average improvement of 40.8 cells/µL (20.6%). Conclusion: The proposed DL model-guided strategy to incorporate the immune system as iOAR in IMRT and IMPT optimization has the potential for radiation-induced lymphopenia mitigation. A prospective clinical trial is planned.

Quantifying the Effect of 4-Dimensional Computed Tomography-Based Deformable Dose Accumulation on Representing Radiation Damage for Patients with Locally Advanced Non-Small Cell Lung Cancer Treated with Standard-Fractionated Intensity-Modulated Radiation Therapy.

PURPOSE: The aim of this study was to investigate the dosimetric and clinical effects of 4-dimensional computed tomography (4DCT)-based longitudinal dose accumulation in patients with locally advanced non-small cell lung cancer treated with standard-fractionated intensity-modulated radiation therapy (IMRT). METHODS AND MATERIALS: Sixty-seven patients were retrospectively selected from a randomized clinical trial. Their original IMRT plan, planning and verification 4DCTs, and ∼4-month posttreatment follow-up CTs were imported into a commercial treatment planning system. Two deformable image registration algorithms were implemented for dose accumulation, and their accuracies were assessed. The planned and accumulated doses computed using average-intensity images or phase images were compared. At the organ level, mean lung dose and normal-tissue complication probability (NTCP) for grade ≥2 radiation pneumonitis were compared. At the region level, mean dose in lung subsections and the volumetric overlap between isodose intervals were compared. At the voxel level, the accuracy in estimating the delivered dose was compared by evaluating the fit of a dose versus radiographic image density change (IDC) model. The dose-IDC model fit was also compared for subcohorts based on the magnitude of NTCP difference (|ΔNTCP|) between planned and accumulated doses. RESULTS: Deformable image registration accuracy was quantified, and the uncertainty was considered for the voxel-level analysis. Compared with planned doses, accumulated doses on average resulted in <1-Gy lung dose increase and <2% NTCP increase (up to 8.2 Gy and 18.8% for a patient, respectively). Volumetric overlap of isodose intervals between the planned and accumulated dose distributions ranged from 0.01 to 0.93. Voxel-level dose-IDC models demonstrated a fit improvement from planned dose to accumulated dose (pseudo-R2 increased 0.0023) and a further improvement for patients with ≥2% |ΔNTCP| versus for patients with <2% |ΔNTCP|. CONCLUSIONS: With a relatively large cohort, robust image registrations, multilevel metric comparisons, and radiographic image-based evidence, we demonstrated that dose accumulation more accurately represents the delivered dose and can be especially beneficial for patients with greater longitudinal response.

Discordance in acute gastrointestinal toxicity between synchrotron-based proton and linac-based electron ultra-high dose rate irradiation.

Purpose: Proton FLASH has been investigated using cyclotron and synchrocyclotron beamlines but not synchrotron beamlines. We evaluated the impact of dose rate (ultra-high [UHDR] vs. conventional [CONV]) and beam configuration (shoot-through [ST] vs. spread-out-Bragg-peak [SOBP]) on acute radiation-induced gastrointestinal toxicity (RIGIT) in mice. We also compared RIGIT between synchrotron-based protons and linac-based electrons with matched mean dose rates. Methods and Materials: We administered abdominal irradiation (12-14 Gy single fraction) to female C57BL/6J mice with an 87 MeV synchrotron-based proton beamline (2 cm diameter field size as a lateral beam). Dose rates were 0.2 Gy/s (S-T pCONV), 0.3 Gy/s (SOBP pCONV), 150 Gy/s (S-T pFLASH), and 230 Gy/s (SOBP pFLASH). RIGIT was assessed by the jejunal regenerating crypt assay and survival. We also compared responses to proton [pFLASH and pCONV] with responses to electron CONV (eCONV, 0.4 Gy/s) and electron FLASH (eFLASH, 188-205 Gy/s). Results: The number of regenerating jejunal crypts at each matched dose was lowest for pFLASH (similar between S-T and SOBP), greater and similar between pCONV (S-T and SOBP) and eCONV, and greatest for eFLASH. Correspondingly, mice that received pFLASH SOBP had the lowest survival rates (50% at 50 days), followed by pFLASH S-T (80%), and pCONV SOBP (90%), but 100% of mice receiving pCONV S-T survived (log-rank P = 0.047 for the four groups). Conclusions: Our findings are consistent with an increase in RIGIT after synchrotron-based pFLASH versus pCONV. This negative proton-specific FLASH effect versus linac-based electron irradiation underscores the importance of understanding the physical and biological factors that will allow safe and effective clinical translation.