A framework for modeling radiation induced lymphopenia in radiotherapy

Lymphocyte radiosensitivity: An extension to the linear-quadratic model?

BACKGROUND AND PURPOSE

The linear-quadratic (LQ) model has been pivotal for evaluating the effects of radiation on cells, but it is primarily characterized by linear responses, which has exhibited limitations when applied to lymphocyte data. The present research aims to address these limitations and to explore an alternative model extended from the conventional LQ model.

MATERIALS AND METHODS

Literature providing lymphocyte counts from assays investigating apoptosis and survival after in vitro irradiation was selected. To address the nonlinearity in lymphocyte responses to radiation, we developed a saturation model characterized by a negative exponential relationship between radiation dose and cellular response. We compared the performance of this saturation model against that of conventional models, including the LQ model and its variants (linear model LM and linear-quadratic-cubic model LQC), as well as the repair-misrepair (RMR) model. The models were evaluated based on prediction-residual plots, residual standard errors, and the Akaike information criterion (AIC). We applied the saturation model to two additional datasets: (1) a dataset from the existing literature that assessed stimulated and unstimulated human lymphocytes exposed to gamma irradiation in vitro and (2) a novel dataset involving T lymphocytes from rodent spleens after exposure to various radiation types (X-rays and protons).

RESULTS

The literature (n = 15 out of 2342) showed that lymphocyte apoptosis varies with dose, time and experimental conditions. The saturation model had a lower AIC of 718 compared to the LM, LQ, LQC and RMR models (AIC of 728, 720, 720 and 734, respectively). The saturation model had a lower residual error and more consistent error distribution. Integrating time as a covariate, the saturation model also had a better AIC for demonstrating time-dependent variations in lymphocyte responses after irradiation. For datasets involving unstimulated lymphocytes before irradiation, the saturation model provided a more accurate fit than did the LM, LQ, and RMR models. In these cases, the fit of the saturation model was comparable to that of the LQC model but offered an advantage when extrapolating to higher doses, where the LQC model might underestimate survival. For stimulated lymphocytes, which are radioresistant, all the models approximated the LM. Both the LQ and saturation models indicated greater radiosensitivity to protons in vitro.

CONCLUSION

The new "saturation model" performed better than the LQ model in quantifying lymphocyte apoptosis and survival, estimating time dependency and assessing the role of radiation modalities or lymphocyte stimulation. Further experiments are warranted to experimentally explore the validity of the saturation model as a promising alternative in the clinical setting.

Prognostic value of lymphocytes in patients with breast cancer receiving radiotherapy after breast-conserving surgery: A post hoc analysis of a phase III randomized trial

PURPOSE

To evaluate the prognostic value of peripheral lymphocyte count (PLC) in the breast cancer patients after breast-conserving surgery (BCS) with radiotherapy (RT).

METHODS AND MATERIALS

This post hoc analysis was performed using data of 628 patients from a phase III, randomized controlled trial comparing hypofractionated RT (HFRT) with conventional fractionated RT (CFRT) after BCS. PLCs were obtained before, during, and after RT until the 1-year follow-up. The optimal cut-off PLCs were determined using the maxstat package in R. Survival rates were estimated using the Kaplan-Meier method and compared with the log-rank test.

RESULTS

A total of 275 (46.1 %) patients developed lymphopenia during RT, among them, 17 (2.8 %) had grade 3 lymphopenia and no one developed grade 4 lymphopenia. With a median follow-up of 110.8 months, patients with pre-RT PLCs of < 1.77 × 109/L had a significantly lower 10-year breast cancer-specific survival (BCSS) rate (P = 0.013) and overall survival (OS) rate (P = 0.026). Patients with a nadir PLC of < 1.35 × 109/L had a significantly poorer 10-year OS rate (P = 0.048). Multivariate analysis showed that a pre-RT PLC of < 1.77 × 109/L was an independent factor influencing BCSS and OS, while the effect of the nadir PLC did not remain significant. Neither PLC nor lymphopenia recovery at post-RT 1, 3, and 6 months and 1 year was associated with survival.

CONCLUSIONS

Radiation-induced lymphopenia in patients with breast cancer after BCS tends to be mild. The lower pre-RT PLC predicted poorer survival.

Treatment-Related Lymphopenia is Possibly a Marker of Good Prognosis in Nasopharyngeal Carcinoma: a Propensity-Score Matching Analysis

Purpose

The aims of the study were to monitor circulating lymphocyte subset counts before and after therapy for nasopharyngeal carcinoma (NPC), and investigate their relationships with patient outcomes.

Patients and Methods

Subjects comprised patients with TNM stage I-IVA NPC who underwent radiotherapy. Peripheral venous blood samples were collected before and after treatment. Lymphocyte subset counts were analyzed by flow cytometry. Differences between post-treatment and baseline counts were calculated to determine Δ values. Patients were divided into high and low groups, based on median lymphocyte subset counts; propensity score matching was applied to balance groups. Progression-free survival (PFS) and overall survival (OS) were plotted using Kaplan-Meier curves and compared using a Log rank test. Relationships between lymphocyte subset counts and patient survival were subjected to Cox regression analysis.

Results

Patients with NPC (n=746) were enrolled from 2012-2022. Higher CD8+ and total T cell baseline counts were associated with better 5-year PFS (73.7% vs 63.1%, P=0.002 and 73.8% vs 64.1%, P=0.005, respectively). Similarly, higher Δ values of CD4+ and total T cells were associated with higher 5-year PFS (76.2% vs 63.5%, P=0.001; 74.3% vs 65.4%, P=0.010) and OS (89.8% vs 81.6%, P=0.005; 88.6% vs 82.5%, P=0.009). Multivariate Cox regression revealed that CD8+ (hazard ratio (HR) 0.651, P=0.002) and total T (HR 0.600, P<0.001) cells were significantly associated with PFS. CD4+ (HR 0.708, P=0.038) and total T (HR 0.639, P=0.031) cells were independent prognostic factors for OS.

Conclusion

NPC patients with low total or CD8+ T cell counts before treatment had worse prognosis; however, those with more significant decreases in total or CD4+ T cells possibly had better outcomes. T cell counts can be reliable indicators to predict prognosis.

Predicting radiation-induced immune suppression in lung cancer patients treated with stereotactic body radiation therapy

BACKGROUND

Stereotactic body radiation therapy (SBRT) is known to modulate the immune system and contribute to the generation of anti-tumor T cells and stimulate T cell infiltration into tumors. Radiation-induced immune suppression (RIIS) is a side effect of radiation therapy that can decrease immunological function by killing naive T cells as well as SBRT-induced newly created effector T cells, suppressing the immune response to tumors and increasing susceptibility to infections.

PURPOSE

RIIS varies substantially among patients and it is currently unclear what drives this variability. Models that can accurately predict RIIS in near real time based on treatment plan characteristics would allow treatment planners to maintain current protocol specific dosimetric criteria while minimizing immune suppression. In this paper, we present an algorithm to predict RIIS based on a model of circulating blood using early stage lung cancer patients treated with SBRT.

METHODS

This Python-based algorithm uses DICOM data for radiation therapy treatment plans, dose maps, patient CT data sets, and organ delineations to stochastically simulate blood flow and predict the doses absorbed by circulating lymphocytes. These absorbed doses are used to predict the fraction of lymphocytes killed by a given treatment plan. Finally, the time dependence of absolute lymphocyte count (ALC) following SBRT is modeled using longitudinal blood data up to a year after treatment. This model was developed and evaluated on a cohort of 64 patients with 10-fold cross validation.

RESULTS

Our algorithm predicted post-treatment ALC with an average error of0.24 ± 0.21 × 10 9 $0.24 \pm 0.21 \times {10}^9$ cells/L with 89% of the patients having a prediction error below 0.5 × 109 cells/L. The accuracy was consistent across a wide range of clinical and treatment variables. Our model is able to predict post-treatment ALC < 0.8 (grade 2 lymphopenia), with a sensitivity of 81% and a specificity of 98%. This model has a ∼38-s end-to-end prediction time of post treatment ALC.

CONCLUSION

Our model performed well in predicting RIIS in patients treated using lung SBRT. With near-real time model prediction time, it has the capability to be interfaced with treatment planning systems to prospectively reduce immune cell toxicity while maintaining national SBRT conformity and plan quality criteria.

An in-silico planning study of stereotactic body radiation therapy for polymetastatic patients with more than ten extra-cranial lesions

Background and purpose

Limited data is available about the feasibility of stereotactic body radiation therapy (SBRT) for treating more than five extra-cranial metastases, and almost no data for treating more than ten. The aim of this study was to investigate the feasibility of SBRT in this polymetatstatic setting.

Materials and methods

Consecutive metastatic melanoma patients with more than ten extra-cranial metastases and a maximum lesion diameter below 11 cm were selected from a single-center prospective registry for this in-silico planning study. For each patient, SBRT plans were generated to treat all metastases with a prescribed dose of 5x7Gy, and dose-limiting organs (OARs) were analyzed. A cell-kill based inverse planning approach was used to automatically determine the maximum deliverable dose to each lesion individually, while respecting all OARs constraints.

Results

A total of 23 polymetastatic patients with a medium of 17 metastases (range, 11-51) per patient were selected. SBRT plans with sufficient target coverage and respected OARs dose constraints were achieved in 16 out of 23 patients. In the remaining seven patients, the lungs V5Gy < 80 % and the liver D700 cm3 < 15Gy were most frequently the dose-limiting constraints. The cell-kill based planning approach allowed optimizing the dose administration depending on metastases total volume and location.

Conclusion

This retrospective planning study shows the feasibility of definitive SBRT for 70% of polymetastatic patients with more than ten extra-cranial lesions and proposes the cell-killing planning approach as an approach to individualize treatment planning in polymetastatic patients'.

Dose-volume metric-based prediction of radiotherapy-induced lymphocyte loss in patients with non-small-cell lung cancer treated with modern radiotherapy techniques

Background and Purpose

Radiation-induced lymphopenia (RIL) is a common side effect of radiotherapy (RT) that may negatively impact survival. We aimed to identify RIL predictors in patients with non-small-cell lung cancer (NSCLC) treated intensity modulated radiotherapy (IMRT) and volumetric modulated arc therapy (VMAT).

Materials and Methods

We retrospectively analysed data of 306 patients who underwent radical RT for NSCLC. Absolute lymphocyte count (ALC) loss was evaluated for each patient by fitting an exponential decay curve to data from first 45 days since treatment start, and percentage ALC loss relative to baseline was calculated based on area under the decay curve and baseline ALC. We compared IMRT and VMAT treatment plans and used linear regression to predict ALC loss.

Results

ALC decreased during RT in the whole patient group, while neutrophil counts remained stable and decreased only in those treated with concurrent chemoradiotherapy (CRT). Percentage ALC loss ranged between 11 and 78 % and was more strongly than lymphocyte nadir correlated with dose-volume metrics for relevant normal structures. We found evidence for the association of high radiation dose to the lungs, heart and body with percentage ALC loss, with lung volume exposed to 20-30 Gy being most important predictors in patients treated with IMRT. A multivariable model based on CRT use, baseline ALC and first principal component (PC1) of the dose-volume predictors showed good predictive performance (bias-corrected R2 of 0.40).

Conclusion

Percentage lymphocyte loss is a robust measure of RIL that is predicted by baseline ALC, CRT use and dose-volume parameters to the lungs, heart and body.

Modeling frameworks for radiation induced lymphopenia: A critical review

Radiation-induced lymphopenia (RIL) is a frequent, and often considered unavoidable, side effect of radiation therapy (RT), whether or not chemotherapy is included. However, in the last few years several studies have demonstrated the detrimental effect of RIL on therapeutic outcomes, with conflicting findings concerning possible inferior patient survival. In addition, since immunotherapeutic treatment has become an integral part of cancer therapy, preserving the immune system is recognized as crucial. Given this background, various research groups have reported on different frameworks for modelling RIL, frequently based on different definitions of RIL itself, and discordant results have been reported. Our aim is to critically review the current literature on RIL modelling and summarize the different approaches recently proposed to improve the prediction of RIL after RT and aimed at immunity-sparing RT. A detailed description of these approaches will be outlined and illustrated through their applications as found in the literature from the last five years. Such a critical analysis represents the necessary starting step to develop an effective strategy that ultimately could harmonize the diverse modelling methods.

A stochastic model of blood flow to calculate blood dose during radiotherapy

Purpose. Lymphopenia is a common side effect in patients treated with radiotherapy, potentially caused by direct cell killing of circulating lymphocytes in the blood. To investigate this hypothesis, a method to assess dose to circulating lymphocytes is needed.Methods. A stochastic model to simulate systemic blood flow in the human body was developed based on a previously designed compartment model. Blood dose was obtained by superimposing the spatiotemporal distribution of blood particles with a time-varying dose rate field, and used as a surrogate for dose to circulating lymphocytes. We discuss relevant theory on compartmental modeling and how to combine it with models of explicit organ vasculature.Results. A general workflow was established which can be used for any anatomical site. Stochastic compartments can be replaced by explicit models of organ vasculatures for improved spatial resolution, and tumor compartments can be dynamically assigned. Generating a patient-specific blood flow distribution takes about one minute, fast enough to investigate the effect of varying treatment parameters such as the dose rate. Furthermore, the anatomical structures contributing most to the overall blood dose can be identified, which could potentially be used for lymphocyte-sparing treatment planning.Conclusion. The ability to report the blood dose distribution during radiotherapy is imperative to test and act upon the current paradigm that radiation-induced lymphopenia is caused by direct cell killing of lymphocytes in the blood. We have built a general model that can do so for various treatment sites. The presented framework is publicly available athttp://github.com/mghro/hedos.

A review on lymphocyte radiosensitivity and its impact on radiotherapy

It is well known that radiation therapy causes lymphopenia in patients and that this is correlated with a negative outcome. The mechanism is not well understood because radiation can have both immunostimulatory and immunosuppressive effects. How tumor dose conformation, dose fractionation, and selective lymph node irradiation in radiation therapy does affect lymphopenia and immune response is an active area of research. In addition, understanding the impact of radiation on the immune system is important for the design and interpretation of clinical trials combining radiation with immune checkpoint inhibitors, both in terms of radiation dose and treatment schedules. Although only a few percent of the total lymphocyte population are circulating, it has been speculated that their increased radiosensitivity may contribute to, or even be the primary cause of, lymphopenia. This review summarizes published data on lymphocyte radiosensitivity based on human, small animal, and in vitro studies. The data indicate differences in radiosensitivity among lymphocyte subpopulations that affect their relative contribution and thus the dynamics of the immune response. In general, B cells appear to be more radiosensitive than T cells and NK cells appear to be the most resistant. However, the reported dose-response data suggest that in the context of lymphopenia in patients, aspects other than cell death must also be considered. Not only absolute lymphocyte counts, but also lymphocyte diversity and activity are likely to be affected by radiation. Taken together, the reviewed data suggest that it is unlikely that radiation-induced cell death in lymphocytes is the sole factor in radiation-induced lymphopenia.

Impact of Radiation Dose to Circulating Immune Cells on Tumor Control and Survival in Esophageal Cancer

Background: The immune system is well known to exert tumor immunosurveillance effects, and that immune cells circulating in the peripheral blood affect tumor prognosis. The study investigated the effect of estimated dose of radiation on circulating immune cells (EDRIC) and tumor control for esophageal cancer patients treated with concurrent chemo-radiotherapy. Materials and Methods: A total of 146 esophageal cancer patients treated with radiotherapy between January 2016 and June 2020 were retrospectively identified. We determined EDRIC, known prognostic factors, and the association of these factors with progression-free survival (PFS) and overall survival (OS). Results: The median follow-up was 17.9 months (2.7-60.4 months). The 3-year OS was 39.2%. Severe post-treatment lymphopenia was observed in 84.2% of patients. A negative correlation between EDRIC and absolute lymphocyte count was found (r = -0.679; p < 0.001). Patients with EDRIC ≥10.3 Gy were more likely to demonstrate grade 4 lymphopenia (55.2% vs. 4.5%; p < 0.001). Patients with grade 4 lymphopenia had a worse OS and PFS. On multivariate analysis, EDRIC was independently associated with OS (hazard ratio [HR], 1.142; p = 0.016) and PFS (HR, 1.121; p = 0.019). Conclusions: EDRIC can predict lymphocyte reduction and poor prognosis for esophageal cancer patients treated with radiotherapy.

SDC1-dependent TGM2 determines radiosensitivity in glioblastoma by coordinating EPG5-mediated fusion of autophagosomes with lysosomes

Glioblastoma multiforme (GBM) is the most common brain malignancy insensitive to radiotherapy (RT). Although macroautophagy/autophagy was reported to be a fundamental factor prolonging the survival of tumors under radiotherapeutic stress, the autophagic biomarkers coordinated to radioresistance of GBM are still lacking in clinical practice. Here we established radioresistant GBM cells and identified their protein profiles using tandem mass tag (TMT) quantitative proteomic analysis. It was found that SDC1 and TGM2 proteins were overexpressed in radioresistant GBM cells and tissues and they contributed to the poor prognosis of RT. Knocking down SDC1 and TGM2 inhibited the fusion of autophagosomes with lysosomes and thus enhanced the radiosensitivity of GBM cells. After irradiation, TGM2 bound with SDC1 and transported it from the cell membrane to lysosomes, and then bound to LC3 through its two LC3-interacting regions (LIRs), coordinating the encounter between autophagosomes and lysosomes, which should be a prerequisite for lysosomal EPG5 to recognize LC3 and subsequently stabilize the STX17-SNAP29-VAMP8 QabcR SNARE complex assembly. Moreover, when combined with RT, cystamine dihydrochloride (a TGM2 inhibitor) extended the lifespan of GBM-bearing mice. Overall, our findings demonstrated the EPG5 tethering mode with SDC1 and TGM2 during the fusion of autophagosomes with lysosomes, providing new insights into the molecular mechanism and therapeutic target underlying radioresistant GBM.Abbreviations: BafA₁: bafilomycin A₁; CQ: chloroquine; Cys-D: cystamine dihydrochloride; EPG5: ectopic P-granules 5 autophagy tethering factor; GBM: glioblastoma multiforme; GFP: green fluorescent protein; LAMP2: lysosomal associated membrane protein 2; LIRs: LC3-interacting regions; MAP1LC3/LC3: microtubule associated protein 1 light chain 3; NC: negative control; RFP: red fluorescent protein; RT: radiotherapy; SDC1: syndecan 1; SNAP29: synaptosome associated protein 29; SQSTM1/p62: sequestosome 1; STX17: syntaxin 17; TGM2: transglutaminase 2; TMT: tandem mass tag; VAMP8: vesicle associated membrane protein 8; WT: wild type.

Modeling of radiation effects to immune system: a review

Cancer metastasis is the major cause of cancer mortality and accounts for about 90% of cancer death. Although radiation therapy has been considered to reduce the localized cancer burden, emerging evidence that radiation can potentially turn tumors into an in situ vaccine has raised significant interest in combining radiation with immunotherapy. However, the combination approach might be limited by the radiation-induced immunosuppression. Assessment of radiation effects on the immune system at the patient level is critical to maximize the systemic antitumor response of radiation. In this review, we summarize the developed solutions in three different categories for systemic radiation therapy: blood dose, radiation-induced lymphopenia, and tumor control. Furthermore, we address how they could be combined to optimize radiotherapy regimens and maximize their synergy with immunotherapy.

Potential Determinants for Radiation-Induced Lymphopenia in Patients With Breast Cancer Using Interpretable Machine Learning Approach

Radiation-induced lymphopenia is known for its survival significance in patients with breast cancer treated with radiation therapy. This study aimed to evaluate the impact of radiotherapy on lymphocytes by applying machine learning strategies. We used Extreme Gradient Boosting (XGboost) to predict the event of lymphopenia (grade≥1) and conduced an independent validation. Then, we induced feature attribution analysis (Shapley additive explanation, SHAP) in explaining the XGboost models to explore the directional contribution of each feature to lymphopenia. Finally, we implemented the proof-of-concept clinical validation. The results showed that the XGboost models had rigorous generalization performances (accuracies 0.764 and ROC-AUC 0.841, respectively) in the independent cohort. The baseline lymphocyte counts are the most protective feature (SHAP = 5.226, direction of SHAP = -0.964). Baseline platelets and monocytes also played important protective roles. The usage of taxane only chemotherapy was less risk on lymphopenia than the combination of anthracycline and taxane. By the contribution analysis of dose, we identified that firstly lymphocytes were sensitive to a radiation dose less than 4Gy; secondly the irradiation volume was more important in promoting lymphopenia than the irradiation dose; thirdly the irradiation dose promoted the event of lymphopenia when the irradiation volume was fixed. Overall, our findings paved the way to clarifying the radiation dose volume effect. To avoid radiation-induced lymphopenia, irradiation volume should be kept to a minimum during the planning process, as long as the target coverage is not compromised.

Radiation Induced Lymphopenia Is Associated With the Effective Dose to the Circulating Immune Cells in Breast Cancer

Background

Lymphopenia is a known significant factor for treatment outcome in cancer patients, with underlying risk factor poorly understood in breast cancer. We hypothesize that the effective dose to the circulating immune cells (EDIC) which was related with lymphopenia in lung cancer will also have significant effect for radiation induced lymphopenia (RIL) in patients with breast cancer.

Material and Methods

Patients treated with adjuvant radiotherapy (RT) and with complete blood tests within one week from RT end/start (post/preRT) were eligible in this study. Radiation dosimetric factors were collected retrospectively, and EDIC for each patient was calculated based on the doses to lung, heart and total body according to the model description, as previously reported. RIL was defined by the CTCAE5.0 based on postRT peripheral lymphocyte count (PLC). Linear regression was first used to test the correlation between EDIC with post/preRT PLC ratio and postRT PLC, using all these as continuous variables. Normal tissue complication probability (NTCP) was used to develop models that predict the CTCAE graded RIL from EDIC.

Results

A total of 735 patients were eligible. The mean post/preRT PLC ratio was 0.66 (95% CI: 0.64-0.68) and mean EDIC of breast cancer was 1.70Gy (95% CI: 1.64-1.75). Both post/preRT PLC ratio and postRT PLC were significantly correlated with EDIC (P<0.001), with R² of 0.246. For patients with normal preRT PLC, the post/preRT PLC ratio was better associated with EDIC, and postRT PLC was expressed as PLC_preRT × (0.89 – 0.16 × EDIC). For patients with preRT lymphopenia, postRT PLC was better associated with EDIC and it was 1.1 – 0.17 × EDIC. Using binned EDIC as the dose variable, the bootstrap validated NTCPs fit the data nicely with R² of 0.93, 0.96, and 0.94 for grade-1, grade-2, and grade-3 RIL, respectively. The corresponding EDIC to induce 50% of grade-1, grade-2 and grade-3 RIL was 1.2, 2.1 and 3.7 Gy, respectively.

Conclusion

EDIC is a significant factor for RIL in patients with breast cancer, and may be used to compute the risk of lymphopenia in each individual patient with the use of the conventional NTCP modeling. External validation is needed before the EDIC can be used to guide RT plan.

Mathematical Modeling to Simulate the Effect of Adding Radiation Therapy to Immunotherapy and Application to Hepatocellular Carcinoma

PURPOSE

To develop a comprehensive framework to simulate the response to immune checkpoint inhibitors (ICIs) in combination with radiation therapy (RT) and to apply the framework for investigating ICI-RT combination regimen in patients with hepatocellular carcinoma (HCC).

METHODS AND MATERIALS

The mechanistic mathematical model is based on dynamic biological interactions between the immune system and the tumor using input data from patient blood samples and outcomes of clinical trials. The cell compartments are described by ordinary differential equations and represent irradiated and nonirradiated tumor cells and lymphocytes. The effect of ICI is modeled using an immune activation term that is based on tumor size changes observed in a phase 1/2 clinical trial for HCC. Simulated combination regimen are based on ongoing ICI-RT trials.

RESULTS

The proposed framework successfully describes tumor volume trajectories observed in early-stage clinical trials of durvalumab monotherapy in patients with HCC. For ICI-RT treatment regimen the irradiated tumor fraction is the most important parameter for the efficacy. For 90% of the tumor cells being irradiated, adding RT to ICI yields an increase in clinical benefit from 33% to 71% in nonirradiated tumor sites. The model agrees with clinical data showing an association of outcome with initial tumor volume and lymphocyte counts. We demonstrate model application in clinical trial design to predict progression-free survival curves, showing that the cohort size to show significant improvement heavily depends on the irradiated tumor fraction.

CONCLUSIONS

We present a framework extending radiation cell kill models to include circulating lymphocytes and the effect of ICIs and enable simulation of combination strategies. The simulations predict that a significant amount of the benefit from RT in combination with ICI stems from the reduction in irradiated tumor burden and associated immune suppression. This aspect needs to be included in the interpretation of outcomes and the design of novel combination trials.

A dynamic blood flow model to compute absorbed dose to circulating blood and lymphocytes in liver external beam radiotherapy

We have developed a novel 4D dynamic liver blood flow model, capable of accurate dose estimation to circulating blood cells during liver-directed external beam radiotherapy, accounting for blood recirculation and radiation delivery time structure. Adult male and adult female liver computational phantoms with detailed vascular trees were developed to include the hepatic arterial, hepatic portal venous, and hepatic venous trees. A discrete time Markov Chain approach was applied to determine the spatiotemporal distribution of 10⁵blood particles (BP) in the human body based on reference values for cardiac output and organ blood volumes. For BPs entering the liver, an explicit Monte Carlo simulation was implemented to track their propagation along ∼2000 distinct vascular pathways through the liver. The model tracks accumulated absorbed dose from time-dependent radiation fields with a 0.1 s time resolution. The computational model was then evaluated for 3 male and 3 female patients receiving photon (VMAT and IMRT) and proton (passive SOBP and active PBS) treatments. The dosimetric impact of treatment modality, delivery time, and fractionation on circulating blood cells was investigated and quantified using the mean dose (μdose,b),V_>0Gy,V_>0.125Gy,andD_2%. Average reductions inμdose,b,V_>0Gy,V_>0.125GyandD_2%of 45%, 6%, 53%, 19% respectively, were observed for proton treatments as compared to photon treatments. Our simulation also showed thatV_>0Gy,V_>0.125Gy, andD_2%were highly sensitive to the beam-on time. BothV_>0GyandV_>0.125Gyincreased with beam-on time, whereasD_2%decreased with increasing beam-on time, demonstrating the tradeoff between low dose to a large fraction of blood cells and high dose to a small fraction of blood cells. Consequently, proton treatments are not necessarily advantageous in terms of dose to the blood simply based on integral dose considerations. Instead, both integral dose and beam-on time can substantially impact relevant dosimetric indices.

Radiation-Induced Lymphopenia Risks of Photon Versus Proton Therapy for Esophageal Cancer Patients

Purpose

To assess possible differences in radiation-induced lymphocyte depletion for esophageal cancer patients being treated with the following 3 treatment modalities: intensity-modulated radiation therapy (IMRT), passive scattering proton therapy (PSPT), and intensity-modulated proton therapy (IMPT).

Methods and Materials

We used 2 prediction models to estimate lymphocyte depletion based on dose distributions. Model I used a piecewise linear relationship between lymphocyte survival and voxel-by-voxel dose. Model II assumes that lymphocytes deplete exponentially as a function of total delivered dose. The models can be fitted using the weekly absolute lymphocyte counts measurements collected throughout treatment. We randomly selected 45 esophageal cancer patients treated with IMRT, PSPT, or IMPT at our institution (15 per modality) to demonstrate the fitness of the 2 models. A different group of 10 esophageal cancer patients who had received PSPT were included in this study of in silico simulations of multiple modalities. One IMRT and one IMPT plan were created, using our standards of practice for each modality, as competing plans to the existing PSPT plan for each patient. We fitted the models by PSPT plans used in treatment and predicted absolute lymphocyte counts for IMRT and IMPT plans.

Results

Model validation on each modality group of patients showed good agreement between measured and predicted absolute lymphocyte counts nadirs with mean squared errors from 0.003 to 0.023 among the modalities and models. In the simulation study of IMRT and IMPT on the 10 PSPT patients, the average predicted absolute lymphocyte count (ALC) nadirs were 0.27, 0.35, and 0.37 K/μL after IMRT, PSPT, and IMPT treatments using Model I, respectively, and 0.14, 0.22, and 0.33 K/μL using Model II.

Conclusions

Proton plans carried a lower predicted risk of lymphopenia after the treatment course than did photon plans. Moreover, IMPT plans outperformed PSPT in terms of predicted lymphocyte preservation.

A Rational Approach to Unilateral Neck RT for Head and Neck Cancers in the Era of Immunotherapy

Radiotherapy plays an important role in the definitive and adjuvant treatment of head and neck squamous cell carcinoma (HNSCC). However, standard courses of radiation therapy may contribute to the depletion of circulating lymphocytes and potentially attenuate optimal tumor antigen presentation that may be detrimental to the efficacy of novel immunotherapeutic agents. This review explores the advantages of restricting radiation to the primary tumor/tumor bed and ipsilateral elective neck as it pertains to the evolving field of immunotherapy.

Chronological Analysis of Acute Hematological Outcomes after Proton and Photon Beam Craniospinal Irradiation in Pediatric Brain Tumors

Purpose

This study aimed to compare the early hematological dynamics and acute toxicities between proton beam craniospinal irradiation (PrCSI) and photon beam craniospinal irradiation (PhCSI) for pediatric brain tumors.

Materials and Methods

We retrospectively reviewed patients with pediatric brain tumors who received craniospinal irradiation (CSI). The average change in hemoglobin levels (ΔHb_avg), absolute lymphocyte counts (ΔALC_avg), and platelet counts (ΔPLT_avg) from baseline values was evaluated and compared between the PrCSI and PhCSI groups at 1 and 2 weeks after the initiation of CSI, 1 week before and at the end of radiotherapy, and 3–4 weeks after the completion of radiotherapy using t-test and mixed-model analysis.

Results

The PrCSI and PhCSI groups consisted of 36 and 30 patients, respectively. There were no significant differences in ΔHb_avg between the two groups at any timepoint. However, ΔALC_avg and ΔPLT_avg were significantly lower in the PhCSI group than in PrCSI group at every timepoint, demonstrating that PrCSI resulted in a significantly lower rate of decline and better recovery of absolute lymphocyte and platelet counts. The rate of grade 3 acute anemia was significantly lower in the PrCSI group than in in the PhCSI group.

Conclusion

PrCSI showed a lower rate of decline and better recovery of absolute lymphocyte and platelet counts than PhCSI in the CSI for pediatric brain tumors. Grade 3 acute anemia was significantly less frequent in the PrCSI group than in the PhCSI group. Further large-scale studies are warranted to confirm these results.

Risk factors for radiation induced lymphopenia in patients with breast cancer receiving adjuvant radiotherapy

Background

This study aimed to investigate radiation-induced lymphopenia and its potential risk factors in patients with breast cancer receiving adjuvant radiotherapy.

Methods

Breast cancer patients received adjuvant radiotherapy (RT) at our hospital with peripheral lymphocyte counts (PLC) at pre-and immediately after RT (post-RT) were eligible. The primary endpoints were any grade of lymphopenia post-RT and nadir-PLC/pre-PLC <0.8. Patient characteristics, tumor factors, and treatment factors were collected for risk assessment. Data are presented as mean and 95% confidence interval (CI) unless otherwise specified. Matched analysis was used to compare the statistical significance between different RT techniques.

Results

A total of 735 consecutive patients met the study criteria. The mean PLC was 1.58×10⁹/L before and 0.99×10⁹/L post-RT (P<0.001). At the end of RT, 60.5% of patients had lymphopenia. Univariate and multivariable logistic analyses showed that RT technique involving RapidArc, mean lung dose, and chemotherapy were significant risk factors (P<0.05) for lymphopenia. RT technique was the only significant risk factor (P<0.05) for nadir-PLC/pre-PLC <0.8. Patients treated with RapidArc had a significantly greater reduction of PLC along with greater V5 of the lungs, even after matching mean lung dose and radiated volume.

Conclusions

Lymphopenia is common in patients with breast cancer after adjuvant RT. RT technique is the only significant factor for lymphopenia and nadir-PLC/pre-PLC <0.8, suggesting the significance of RT technique choice to minimize lymphopenia and improve treatment outcomes.

HEDOS-a computational tool to assess radiation dose to circulating blood cells during external beam radiotherapy based on whole-body blood flow simulations

We have developed a time-dependent computational framework, hematological dose (HEDOS), to estimate dose to circulating blood cells from radiation therapy treatment fields for any treatment site. Two independent dynamic models were implemented in HEDOS: one describing the spatiotemporal distribution of blood particles (BPs) in organs and the second describing the time-dependent radiation field delivery. A whole-body blood flow network based on blood volumes and flow rates from ICRP Publication 89 was simulated to produce the spatiotemporal distribution of BPs in organs across the entire body using a discrete-time Markov process. Constant or time-varying transition probabilities were applied and their impact on transition time was investigated. The impact of treatment time and anatomical site were investigated using imaging data and dose distributions from a liver cancer and a brain cancer patient. The simulations revealed different dose levels to the circulating blood for brain irradiation compared to liver irradiation even for similar field sizes due to the different blood flow properties of the two organs. The volume of blood receiving any dose (V>0 Gy) after a single radiation fraction increases from 1.2% for a 1 s delivery time to 20.9% for 120 s delivery time for the brain cancer treatment, and from 10% (1 s) to 48.7% (120 s) for a liver cancer treatment. Other measures of the low-dose bath to the circulating blood such as the dose to small volumes of blood (D2%) decreases with longer delivery time. Furthermore, we demonstrate that the blood dose-volume histogram is highly sensitive to changes in the treatment time, indicating that dynamic modeling of blood flow and radiation fields is necessary to evaluate dose to circulating blood cells for the assessment of radiation-induced lymphopenia. HEDOS is publicly available and allows for the estimation of patient-specific dose to circulating blood cells based on organ DVHs, thus enabling the study of the impact of different treatment plans, dose rates, and fractionation schemes.

Radiotherapy planning parameters correlate with changes in the peripheral immune status of patients undergoing curative radiotherapy for localized prostate cancer

Purpose

The influence of radiotherapy on patient immune cell subsets has been established by several groups. Following a previously published analysis of immune changes during and after curative radiotherapy for prostate cancer, this analysis focused on describing correlations of changes of immune cell subsets with radiation treatment parameters.

Patients and methods

For 13 patients treated in a prospective trial with radiotherapy to the prostate region (primary analysis) and five patients treated with radiotherapy to prostate and pelvic nodal regions (exploratory analysis), already published immune monitoring data were correlated with clinical data as well as radiation planning parameters such as clinical target volume (CTV) and volumes receiving 20 Gy (V20) for newly contoured volumes of pelvic blood vessels and bone marrow.

Results

Most significant changes among immune cell subsets were observed at the end of radiotherapy. In contrast, correlations of age and CD8⁺ subsets (effector and memory cells) were observed early during and 3 months after radiotherapy. Ratios of T cells and T cell proliferation compared to baseline correlated with CTV. Early changes in regulatory T cells (Treg cells) and CD8⁺ effector T cells correlated with V20 of blood vessels and bone volumes.

Conclusions

Patient age as well as radiotherapy planning parameters correlated with immune changes during radiotherapy. Larger irradiated volumes seem to correlate with early suppression of anti-cancer immunity. For immune cell analysis during normofractionated radiotherapy and correlations with treatment planning parameters, different time points should be looked at in future projects.

Trial registration number:

NCT01376674, 20.06.2011

Supplementary Information

The online version contains supplementary material available at 10.1007/s00262-021-03002-6.

Modeling Radioimmune Response-Current Status and Perspectives

The combination of immune therapy with radiation offers an exciting and promising treatment modality in cancer therapy. It has been hypothesized that radiation induces damage signals within the tumor, making it more detectable for the immune system. In combination with inhibiting immune checkpoints an effective anti-tumor immune response may be established. This inversion from tumor immune evasion raises numerous questions to be solved to support an effective clinical implementation: These include the optimum immune drug and radiation dose time courses, the amount of damage and associated doses required to stimulate an immune response, and the impact of lymphocyte status and dynamics. Biophysical modeling can offer unique insights, providing quantitative information addressing these factors and highlighting mechanisms of action. In this work we review the existing modeling approaches of combined ‘radioimmune’ response, as well as associated fields of study. We propose modeling attempts that appear relevant for an effective and predictive model. We emphasize the importance of the time course of drug and dose delivery in view to the time course of the triggered biological processes. Special attention is also paid to the dose distribution to circulating blood lymphocytes and the effect this has on immune competence.

Radiotherapy-Related Lymphopenia Affects Overall Survival in Patients With Lung Cancer

INTRODUCTION

Lymphopenia after radiotherapy has an adverse effect on the patient's outcome. However, the relationship between radiotherapy dose delivery and lymphopenia is not fully understood. This work used image-based data mining to identify anatomical regions where the received dose is correlated with severe lymphopenia.

METHODS

A total of 901 patients with lung cancer were analyzed. A Cox model was used to assess prognostic factors of overall survival (OS). Two matched groups were defined-patients with lymphopenia of grade 3 or higher and patients without lymphopenia of grade 3-based on tumor volume, baseline lymphocytes, and prescribed dose. Then, data mining was used to identify regions where dose correlates significantly with lymphopenia of grade 3 or higher. For this, dose matrices were aligned using registration of the computed tomography images to one reference patient. Mean dose distributions were obtained for the two groups, and organs of significance were detected. Dosimetric parameters from the identified organs that had the highest correlation with lymphocytes at nadir were selected. Multivariable analysis was conducted for lymphopenia of grade 3 or higher on the full lung cohort, and the model was tested on 305 patients with esophageal cancer.

RESULTS

Adjusted Cox regression revealed that lymphopenia of grade 3 or higher is an independent factor of OS. The anatomical regions identified were the heart, lung, and thoracic vertebrae. Dosimetric parameters for lymphopenia included thoracic vertebrae V20, mean lung dose, and mean heart dose, which were further validated in the esophageal cancer cohort.

CONCLUSIONS

We report that severe lymphopenia during radiotherapy is a poor prognostic factor for OS in patients with lung cancer and could be mitigated by minimizing thoracic vertebrae V20, mean lung dose, and mean heart dose to limit the irradiation of stem cells and blood pool.

Novel Methodology to Investigate the Effect of Radiation Dose to Heart Substructures on Overall Survival

PURPOSE

For patients with lung cancer treated with radiation therapy, a dose to the heart is associated with excess mortality; however, it is often not feasible to spare the whole heart. Our aim is to define cardiac substructures and dose thresholds that optimally reduce early mortality.

METHODS AND MATERIALS

Fourteen cardiac substructures were delineated on 5 template patients with representative anatomies. One thousand one hundred sixty-one patients with non-small cell lung cancer were registered nonrigidly to these 5 template anatomies, and their radiation therapy doses were mapped. Mean and maximum dose to each substructure were extracted, and the means were evaluated as input to prediction models. The cohort was bootstrapped into 2 variable reduction techniques: elastic net least absolute shrinkage and selection operator and the random survival forest model. Each method was optimized to extract variables contributing most to overall survival, and model coefficients were evaluated to select these substructures. The most important variables common to both models were selected and evaluated in multivariable Cox-proportional hazard models. A threshold dose was defined, and Kaplan-Meier survival curves plotted.

RESULTS

Nine hundred seventy-eight patients remained after visual quality assurance of the registration. Ranking the model coefficients across the bootstraps selected the maximum dose to the right atrium, right coronary artery, and ascending aorta as the most important factors associated with survival. The maximum dose to the combined cardiac region showed significance in the multivariable model, a hazard ratio of 1.01/Gy, and P = .03 after accounting for tumor volume (P < .001), N stage (P < .01), and performance status (P = .01). The optimal threshold for the maximum dose, equivalent dose in 2-Gy fractions, was 23 Gy. Kaplan-Meier survival curves showed a significant split (log-rank P = .008).

CONCLUSIONS

The maximum dose to the combined cardiac region encompassing the right atrium, right coronary artery, and ascending aorta was found to have the greatest effect on patient survival. A maximum equivalent dose in 2-Gy fractions of 23 Gy was identified for consideration as a dose limit in future studies.

Ultra-high dose rate effect on circulating immune cells: A potential mechanism for FLASH effect?

PURPOSE

"FLASH" radiotherapy (RT) is a potential paradigm-changing RT technology with marked tumor killing and normal tissue sparing. However, the mechanism of the FLASH effect is not well understood. We hypothesize that the ultra-high dose rate FLASH-RT significantly reduces the killing of circulating immune cells which may partially contribute to the reported FLASH effect.

METHODS

This computation study directly models the effect of radiation dose rate on the killing of circulating immune cells. The model considers an irradiated volume that takes up A% of cardiac output and contains B% of total blood. The irradiated blood volume and dose were calculated for various A%, B%, blood circulation time, and irradiation time (which depends on the dose rate). The linear-quadratic model was used to calculate the extent of killing of circulating immune cells at ultra-high vs. conventional dose rates.

RESULTS

A strong sparing effect on circulating blood cells by FLASH-RT was noticed; i.e., killing of circulating immune cells reduced from 90% to 100% at conventional dose rates to 5-10% at ultra-high dose rates. The threshold FLASH dose rate was determined to be ~40 Gy/s for mice in an average situation (A% = 50%), consistent with the reported FLASH dose rate in animal studies, and it was approximately one order of magnitude lower for humans than for mice. The magnitude of this sparing effect increased with the dose/fraction, reached a plateau at 30-50 Gy/fraction, and almost completely vanished at 2 Gy/fraction.

CONCLUSION

We have calculated a strong sparing effect on circulating immune cells by FLASH-RT, which may contribute to the reported FLASH effects in animal studies.

Radiation-induced lymphopenia: the main aspects to consider in immunotherapy trials for endometrial and cervical cancer patients

BACKGROUND

Although the chemotherapy-induced depletion of circulating white blood cells (WBC) is well recognized, the impact of exclusive radiotherapy (RT) on the different subpopulations of WBC remains unexplored. This may be important for immunotherapy administrated in combination with radiation, especially in malignant tumors usually treated with RT or chemoradiotherapy (CRT), and characterized by a high mutational burden, such as endometrial (EC) or cervical cancer (CC). We aimed to evaluate the impact of RT and CRT on circulating WBC in uterine cancers and its correlation with survival.

MATERIAL AND METHODS

A total of 202 consecutive patients with uterine cancers treated with RT or CRT between 2009 and 2016 in a large European center and with available basal and post-treatment blood tests were retrospectively evaluated. EC and CC patients were analyzed separately. The differences between pre- and post- treatment WBC mean values were evaluated independently in patients treated with CRT and exclusive RT. Two-sided T test for paired samples and Kaplan-Meier curves were applied for analysis (p value < 0.05, SPSS v.23).

RESULTS

Among EC patients, 29 received CRT and 34 exclusive postoperative RT, while in CC cohort, 105 were treated with CRT and 34 with RT. In both cohorts, CRT affected significantly all WBC subtypes, whereas exclusive RT decreased only lymphocytes population (p = 0.000). Radiation-induced lymphopenia (RIL) had no impact on survival outcomes.

CONCLUSIONS

The selective depletion of lymphocytes after RT was significant in both EC and CC. Our results are of interest for further research on RIL and for design of immunotherapy-based clinical trials.

Radiation-Induced Lymphopenia Predicts Poorer Prognosis in Patients With Breast Cancer: A Post Hoc Analysis of a Randomized Controlled Trial of Postmastectomy Hypofractionated Radiation Therapy

PURPOSE

The aim of this study was to determine whether radiation-induced lymphopenia affects the survival of patients with breast cancer.

METHODS AND MATERIALS

Post hoc analysis was conducted on data from 598 patients with breast cancer from a randomized controlled trial comparing postmastectomy hypofractionated radiation therapy (HFRT; 43.5 Gy in 15 fractions over 3 weeks) with conventional fractionated radiation therapy (CFRT; 50 Gy in 25 fractions over 5 weeks). Mean peripheral lymphocyte count (PLC) at different time points in the 2 groups was compared by the t test. Disease-free survival and overall survival were analyzed by the Kaplan-Meier method and compared between groups by the log-rank test.

RESULTS

Baseline PLC (pre-PLC) was comparable between HFRT and CFRT patients (1.60 ± 0.57 × 10⁹/L vs 1.56 ± 0.52 × 10⁹/L; P = .33). In both groups, the PLC declined steadily during the course of radiation therapy but started to recover at 1 month after radiation therapy. Incidence of lymphopenia was significantly lower in HFRT patients (45.4% vs 55.7%; P = .01). Nadir-PLC was significantly higher in HFRT patients (1.08 ± 0.37 × 10⁹/L vs 0.97 ± 0.31× 10⁹/L; P < .001), as was the nadir-PLC/pre-PLC ratio (0.72 ± 0.28 vs 0.67 ± 0.28; P = .02). Median follow-up was 57.6 months (interquartile range, 38.5-81.4). The 5-year disease-free survival was significantly lower in patients with a nadir-PLC/pre-PLC ratio <0.8 than in those with a ratio ≥0.8 (71.8% vs 82.6%; P = .01); however, overall survival was comparable between the groups (85.8% vs 90.6%; P = .24).

CONCLUSIONS

The risk of radiation-induced lymphopenia in patients with breast cancer is lower with HFRT than with CFRT. A low nadir-PLC/pre-PLC ratio may predict poor prognosis.