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.

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.

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.

Renormalization of radiobiological response functions by energy loss fluctuations and complexities in chromosome aberration induction: deactivation theory for proton therapy from cells to tumor control.

We employ a multi-scale mechanistic approach built upon our recent phenomenological/computational methodologies [R. Abolfath et al., Sci. Rep. 7, 8340 (2017)] to investigate radiation induced cell toxicities and deactivation mechanisms as a function of linear energy transfer in hadron therapy. Our theoretical model consists of a system of Markov chains in microscopic and macroscopic spatio-temporal landscapes, i.e., stochastic birth-death processes of cells in millimeter-scale colonies that incorporates a coarse-grained driving force to account for microscopic radiation induced damage. The coupling, hence the driving force in this process, stems from a nano-meter scale radiation induced DNA damage that incorporates the enzymatic end-joining repair and mis-repair mechanisms. We use this model for global fitting of the high-throughput and high accuracy clonogenic cell-survival data acquired under exposure of the therapeutic scanned proton beams, the experimental design that considers γ-H2AX as the biological endpoint and exhibits maximum observed achievable dose and LET, beyond which the majority of the cells undergo collective biological deactivation processes. An estimate to optimal dose and LET calculated from tumor control probability by extension to ~106 cells per mm-size voxels is presented. We attribute the increase in degree of complexity in chromosome aberration to variabilities in the observed biological responses as the beam linear energy transfer (LET) increases, and verify consistency of the predicted cell death probability with the in vitro cell survival assay of approximately 100 non-small cell lung cancer (NSCLC) cells. The present model provides an interesting interpretation to variabilities in α and ß indices via perturbative expansion of the cell survival fraction (SF) in terms of specific and lineal energies, z and y, corresponding to continuous transitions in pair-wise to ternary, quaternary and more complex recombination of broken chromosomes from the entrance to the end of the range of proton beam.