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.

Voxel-based analysis: Roadmap for clinical translation.

Voxel-based analysis (VBA) allows the full, 3-dimensional, dose distribution to be considered in radiotherapy outcome analysis. This provides new insights into anatomical variability of pathophysiology and radiosensitivity by removing the need for a priori definition of organs assumed to drive the dose response associated with patient outcomes. This approach may offer powerful biological insights demonstrating the heterogeneity of the radiobiology across tissues and potential associations of the radiotherapy dose with further factors. As this methodological approach becomes established, consideration needs to be given to translating VBA results to clinical implementation for patient benefit. Here, we present a comprehensive roadmap for VBA clinical translation. Technical validation needs to demonstrate robustness to methodology, where clinical validation must show generalisability to external datasets and link to a plausible pathophysiological hypothesis. Finally, clinical utility requires demonstration of potential benefit for patients in order for successful translation to be feasible. For each step on the roadmap, key considerations are discussed and recommendations provided for best practice.

Stereotactic Radiation Therapy for Brain Metastases: Factors Affecting Outcomes and Radiation Necrosis.

Stereotactic radiation therapy (SRT) is a proven effective treatment for brain metastases (BM); however, symptomatic radiation necrosis (RN) is a late effect that may impact on patient's quality of life. The aim of our study was to retrospectively evaluate survival outcomes and characterize the occurrence of RN in a cohort of BM patients treated with ablative SRT at Federico II University Hospital. Clinical and dosimetric factors of 87 patients bearing a total of 220 BMs treated with SRT from 2016 to 2022 were analyzed. Among them, 46 patients with 127 BMs having clinical and MRI follow-up (FUP) ≥ 6 months were selected for RN evaluation. Dosimetric parameters of the uninvolved brain (brain without GTV) were extracted. The crude local control was 91% with neither clinical factors nor prescription dose correlating with local failure (LF). At a median FUP of 9 (1-68) months, the estimated median overall survival (OS), progression-free survival (PFS), and brain progression-free survival (bPFS) were 16, 6, and 9 months, respectively. The estimated OS rates at 1 and 3 years were 59.8% and 18.3%, respectively; bPFS at 1 and 3 years was 29.9% and 13.5%, respectively; PFS at 1 and 3 years was 15.7% and 0%, respectively; and local failure-free survival (LFFS) at 1 and 3 years was 87.2% and 83.8%, respectively. Extracranial disease status was an independent factor related to OS. Fourteen (30%) patients manifested RN. At multivariate analysis, adenocarcinoma histology, left location, and absence of chemotherapy were confirmed as independent risk factors for any-grade RN. Nine (20%) patients developed symptomatic (G2) RN, which improved or stabilized after 1-16 months of steroid therapy. With prompt recognition and, when necessary, medical therapy, RN radiological and clinical amelioration can be obtained.

Technical note: MAMBA-Multi-pAradigM voxel-Based Analysis: A computational cookbot.

BACKGROUND: Voxel-Based (VB) analysis embraces a multifaceted ensemble of sophisticated techniques, lying at the boundary between image processing and statistical modeling, that allow for a frequentist inference of pathophysiological properties anchored to an anatomical reference. VB methods has been widely adopted in neuroimaging studies and, more recently, they are gaining momentum in radiation oncology research. However, the price for the power of VB analysis is the complexity of the underlying mathematics and algorithms. PURPOSE: In this paper, we present the Multi-pAradigM voxel-Based Analysis (MAMBA) toolbox, which is intended for a flexible application of VB analysis in a wide variety of scenarios in medical imaging and radiation oncology. METHODS: The MAMBA toolbox is implemented in Matlab. It provides open-source functions to compute VB statistical models of the input data, according to a great variety of regression schemes, and to derive VB maps of the observed significance level, performing a non-parametric permutation inference. The toolbox allows for including VB and global outcomes, as well as an arbitrary amount of VB and global Explanatory Variables (EVs). In addition, the Matlab Parallel Computing Toolbox is exploited to take advantage of the perfect parallelizability of most workloads. RESULTS: The use of MAMBA was demonstrated by means of several realistic examples on a synthetic dataset mimicking a radiation oncology scenario. CONCLUSION: MAMBA is an open-source toolbox, freely available for academic and non-commercial purposes. It is designed to make state-of-the-art VB analysis accessible to research scientists without the programming resources needed to build from scratch their own software solutions. At the same time, the source code is handed out for more experienced users to complement their own tools, also customizing user-defined models. MAMBA guarantees high generality and flexibility in the design of the statistical models, significantly expanding on the features of available free tools for VB analysis. The presented toolbox aims at increasing the reach of VB studies as well as the sharing of research results.

Potential benefits of using radioactive ion beams for range margin reduction in carbon ion therapy.

Sharp dose gradients and high biological effectiveness make ions such as 12C an ideal tool to treat deep-seated tumors, however, at the same time, sensitive to errors in the range prediction. Tumor safety margins mitigate these uncertainties, but during the irradiation they lead to unavoidable damage to the surrounding healthy tissue. To fully exploit the Bragg peak benefits, a large effort is put into establishing precise range verification methods. Despite positron emission tomography being widely in use for this purpose in 12C therapy, the low count rates, biological washout, and broad activity distribution still limit its precision. Instead, radioactive beams used directly for treatment would yield an improved signal and a closer match with the dose fall-off, potentially enabling precise in vivo beam range monitoring. We have performed a treatment planning study to estimate the possible impact of the reduced range uncertainties, enabled by radioactive 11C ions treatments, on sparing critical organs in tumor proximity. Compared to 12C treatments, (i) annihilation maps for 11C ions can reflect sub- millimeter shifts in dose distributions in the patient, (ii) outcomes of treatment planning with 11C significantly improve and (iii) less severe toxicities for serial and parallel critical organs can be expected.

Radiation Therapy in Thoracic Tumors: Recent Trends and Current Issues.

Radiation therapy (RT) plays a fundamental role in the multidisciplinary treatment and management of thoracic cancers, and in particular, RT is the most used non-surgical treatment modality for lung cancer, which in turn is the most common type of thoracic malignancy [...].

Radiation-Induced Esophagitis in Non-Small-Cell Lung Cancer Patients: Voxel-Based Analysis and NTCP Modeling.

The aim of our study is to characterize the risk of radiation-induced esophagitis (RE) in a cohort of Non-Small-Cell Lung Cancer (NSCLC) patients treated with concurrent chemotherapy and photon/proton therapy. For each patient, the RE was graded according to the CTCAE v.3. The esophageal dose-volume histograms (DVHs) were extracted. Voxel-based analyses (VBAs) were performed to assess the spatial patterns of the dose differences between patients with and without RE of grade ≥ 2. Two hierarchical NTCP models were developed by multivariable stepwise logistic regression based on non-dosimetric factors and on the DVH metrics for the whole esophagus and its anatomical subsites identified by the VBA. In the 173 analyzed patients, 76 (44%) developed RE of grade ≥ 2 at a median follow-up time of 31 days. The VBA identified regions of significant association between dose and RE in a region encompassing the thoracic esophagus. We developed two NTCP models, including the RT modality and a dosimetric factor: V55Gy for the model related to the whole esophagus, and the mean dose for the model designed on the thoracic esophagus. The cross-validated performance showed good predictions for both models (ROC-AUC of 0.70 and 0.73, respectively). The only slight improvement provided by the analysis of the thoracic esophageal subsites might be due to the relevant sparing of cervical and lower thoracic esophagus in the analyzed cohort. Further studies on larger cohorts and a more heterogeneous set of dose distributions are needed to validate these preliminary findings and shed further light on the spatial patterns of RE development.

On the interplay between dosiomics and genomics in radiation-induced lymphopenia of lung cancer patients.

PURPOSE: To investigate the interplay between spatial dose patterns and single nucleotide polymorphisms in the development of radiation-induced lymphopenia (RIL) in 186 non-small-cell lung cancer (NSCLC) patients undergoing chemo-radiotherapy (RT). METHODS: This study included NSCLC patients enrolled in a randomized trial of protons vs. photons with available absolute lymphocyte counts at baseline and during RT and XRCC1-rs25487 genotyping data. After masking the GTV, planning CT scans and dose maps were spatially normalized to a common anatomical reference. A Voxel-Based Analysis (VBA) was performed to assess voxel-wise relationships of dosiomic and genomic explanatory variables with RIL. The underlying generalized linear model was designed to include both the explanatory variables (3D dose distributions and the XRCC1-rs25487 genotypes) and possible nuisance variables significantly correlated with RIL. The maps of model coefficients as well as their significance maps were generated. RESULTS: Measures for RIL definition during RT were characterized, including kinetic parameters for lymphocyte loss. The VBA generated three-dimensional maps of correlation between RIL and dose in lymphoid organs as well as organs with abundant blood pools. The identified voxel-wise relationships account for XRCC1-rs25487 polymorphism and demonstrate the variant AA genotype being detrimental to lymphocyte depletion (p = 0.03). CONCLUSION: The performed analyses blindly highlighted relevant anatomical regions that contributed most to lymphocyte depletion during RT and the interplay of the variant XRCC1-rs25487 AA genotype with the dose delivered to the primary lymphoid organs. These findings may help to guide the development of dosimetric RIL mitigation strategies for the application of effective individualized RT.

Local tumor control and treatment related toxicity after plaque brachytherapy for uveal melanoma: A systematic review and a data pooled analysis.

Uveal melanoma (UM) represents the most common primary intraocular tumor, and nowadays eye plaque brachytherapy (EPB) is the most frequently used visual acuity preservation treatment option for small to medium sized UMs. The excellent local tumor control (LTC) rate achieved by EPB may be associated with severe complications and adverse events. Several dosimetric and clinical risk factors for the development of EPB-related ocular morbidity can be identified. However, morbidity predictive models specifically developed for EPB are still scarce. PRISMA methodology was used for the present systematic review of articles indexed in PubMed in the last sixteen years on EPB treatment of UM which aims at determining the major factors affecting local tumor control and ocular morbidities. To our knowledge, for the first time in EPB field, local tumor control probability (TCP) and normal tissue complication probability (NTCP) modelling on pooled clinical outcomes were performed. The analyzed literature (103 studies including 21,263 UM patients) pointed out that Ru-106 EPB provided high local control outcomes while minimizing radiation induced complications. The use of treatment planning systems (TPS) was the most influencing factor for EPB outcomes such as metastasis occurrence, enucleation, and disease specific survival, irrespective of radioactive implant type. TCP and NTCP parameters were successfully extracted for 5-year LTC, cataract and optic neuropathy. In future studies, more consistent recordings of ocular morbidities along with accurate estimation of doses through routine use of TPS are needed to expand and improve the robustness of toxicity risk prediction in EPB.

Radiation-Induced Dyspnea in Lung Cancer Patients Treated with Stereotactic Body Radiation Therapy.

In this study, we investigated the prognostic factors for radiation-induced dyspnea after hypo-fractionated radiation therapy (RT) in 106 patients treated with Stereotactic Body RT for Non-Small-Cell Lung Cancer (NSCLC). The median prescription dose was 50 Gy (range: 40-54 Gy), delivered in a median of four fractions (range: 3-12). Dyspnea within six months after SBRT was scored according to CTCAE v.4.0. Biologically Effective Dose (α/ß = 3 Gy) volume histograms for lungs and heart were extracted. Dosimetric parameters along with patient-specific and treatment-related factors were analyzed, multivariable logistic regression method with Leave-One-Out (LOO) internal validation applied. Model performance was evaluated by the area under the receiver operating characteristic (ROC) curve (AUC) and calibration plot parameters. Fifty-seven patients (53.8%) out of 106 developed dyspnea of any grade after SBRT (25/57 grade ≥ 2 cases). A three-variable predictive model including patient comorbidity (COPD), heart volume and the relative lungs volume receiving more than 15 Gy was selected. The model displays an encouraging performance given by a training ROC-AUC = 0.71 [95%CI 0.61-0.80] and a LOO-ROC-AUC = 0.64 [95%CI 0.53-0.74]. Further modeling efforts are needed for dyspnea prediction in hypo-fractionated treatments in order to identify patients at high risk for developing lung toxicity more accurately.

Radiation Pneumonitis in Thoracic Cancer Patients: Multi-Center Voxel-Based Analysis.

This study investigates the dose-response patterns associated with radiation pneumonitis (RP) in patients treated for thoracic malignancies with different radiation modalities. To this end, voxel-based analysis (VBA) empowered by a novel strategy for the characterization of spatial properties of dose maps was applied. Data from 382 lung cancer and mediastinal lymphoma patients from three institutions treated with different radiation therapy (RT) techniques were analyzed. Each planning CT and biologically effective dose map (α/ß = 3 Gy) was spatially normalized on a common anatomical reference. The VBA of local dose differences between patients with and without RP was performed and the clusters of voxels with dose differences that significantly correlated with RP at a p-level of 0.05 were generated accordingly. The robustness of VBA inference was evaluated by a novel characterization for spatial properties of dose maps based on probabilistic independent component analysis (PICA) and connectograms. This lays robust foundations to the obtained findings that the lower parts of the lungs and the heart play a prominent role in the development of RP. Connectograms showed that the dataset can support a radiobiological differentiation between the main heart and lung substructures.

RESUMEN: A flexible class of multi-parameter qMRI protocols.

PURPOSE: To introduce a class of fast 3D quantitative MRI (qMRI) schemes (RESUMEN, for N=1,…,4) that allow for a thorough characterization of microstructural properties of brain tissues. METHODS: An arbitrary multi-echo GRE acquisition optimized for quantitative susceptibility mapping (QSM) is complemented with an appropriate low flip-angle GRE sequence drawn from four possible choices. The acquired signals are processed to analytically derive the longitudinal relaxation (R1) and free induction decay (R2∗) rates, as well as the proton density (PD) and QSM. A comprehensive modeling of the excitation and B1- profiles and of the RF-spoiling is included in the acquisition and processing pipeline. RESULTS: The RESUMEN maps appear homogeneous throughout the field-of-view and exhibit comparable values and high SNR across the considered range of N values. CONCLUSIONS: The introduced schemes represent a class of robust and flexible strategies to derive a thorough and fast qMRI study, suitable for a whole-brain acquisition with isotropic voxel resolution of 700 µm in less than 15 min.

Probing thoracic dose patterns associated to pericardial effusion and mortality in patients treated with photons and protons for locally advanced non-small-cell lung cancer.

PURPOSE: To investigate thoracic dose-response patterns for pericardial effusion (PCE) and mortality in patients treated for locally advanced Non-Small-Cell Lung Cancer (NSCLC) by Intensity Modulated RT (IMRT) or Passive-Scattering Proton Therapy (PSPT). METHODS: Among 178 patients, 43.5% developed grade ≥ 2 PCE. Clinical and dosimetric factors associated with PCE or overall survival (OS) were identified via multi-variable Cox proportional hazards modeling. The Voxel-Based Analyses (VBAs) of local dose differences between patients with and without PCE and mortality was performed. The robustness of VBA results was assessed by a novel characterization of spatial properties of dose distributions based on probabilistic independent component analysis (PICA) and connectograms. RESULTS: Several non-dosimetric variables were selected by the multivariable analysis for the considered outcomes, while the time-dependent PCE onset was uncorrelated with the OS (p = 0.34) at a multi-variable Cox analysis. Despite the significant PSPT dosimetric advantage, the RT technique did not affect the occurrence of PCE or OS. VBAs highlighted largely overlapping clusters significantly associated with PCE endpoints in heart and lungs. No significant dosimetric patterns related to mortality endpoints were found. PICA identified 43 components homogeneously scattered within thorax, while connectograms showed modest correlations between doses in main cardio-pulmonary substructures. CONCLUSIONS: Spatially resolved analysis highlighted dose patterns related to radiation-induced cardiac toxiciy and the observed organ-based dose-response mismatch in PSPT and IMRT. Indeed, the thoracic regions spared by PSPT poorly overlapped with the areas involved in PCE development, as highlited by VBA. PICA and connectograms proved valuable tools for assessing the robusteness of obtained VBA inferences.

NTCP Models for Severe Radiation Induced Dermatitis After IMRT or Proton Therapy for Thoracic Cancer Patients.

Radiation therapy (RT) of thoracic cancers may cause severe radiation dermatitis (RD), which impacts on the quality of a patient's life. Aim of this study was to analyze the incidence of acute RD and develop normal tissue complication probability (NTCP) models for severe RD in thoracic cancer patients treated with Intensity-Modulated RT (IMRT) or Passive Scattering Proton Therapy (PSPT). We analyzed 166 Non-Small-Cell Lung Cancer (NSCLC) patients prospectively treated at a single institution with IMRT (103 patients) or PSPT (63 patients). All patients were treated to a prescribed dose of 60 to 74 Gy in conventional daily fractionation with concurrent chemotherapy. RD was scored according to CTCAE v3 scoring system. For each patient, the epidermis structure (skin) was automatically defined by an in house developed segmentation algorithm. The absolute dose-surface histogram (DSH) of the skin were extracted and normalized using the Body Surface Area (BSA) index as scaling factor. Patient and treatment-related characteristics were analyzed. The Lyman-Kutcher-Burman (LKB) NTCP model recast for DSH and the multivariable logistic model were adopted. Models were internally validated by Leave-One-Out method. Model performance was evaluated by the area under the receiver operator characteristic curve, and calibration plot parameters. Fifteen of 166 (9%) patients developed severe dermatitis (grade 3). RT technique did not impact RD incidence. Total gross tumor volume (GTV) size was the only non dosimetric variable significantly correlated with severe RD (p = 0.027). Multivariable logistic modeling resulted in a single variable model including S 20Gy, the relative skin surface receiving more than 20 Gy (OR = 31.4). The cut off for S 20Gy was 1.1% of the BSA. LKB model parameters were TD50 = 9.5 Gy, m = 0.24, n = 0.62. Both NTCP models showed comparably high prediction and calibration performances. Despite skin toxicity has long been considered a potential limiting factor in the clinical use of PSPT, no significant differences in RD incidence was found between RT modalities. Once externally validated, the availability of NTCP models for prediction of severe RD may advance treatment planning optimization.

Dose-Response of TLD-100 in the Dose Range Useful for Hypofractionated Radiotherapy.

PURPOSE: The aim of the study was to exploit the feasibility of thermoluminescent dosimeters (TLDs) in radiation therapy techniques in which high dose per fraction is involved. METHODS: Dose-response of TLD-100 (LiF: Mg, Ti) was investigated in both 6-MV photon and 6-MeV electron beams. The element correction factor (ECF) generation method was applied to check the variability of the TLDs response. Two batches of 50 TLDs were divided into groups and exposed in the dose range 0 to 30 Gy. Regression analysis was performed with both linear and quadratic models. For each irradiation beam, the calibration curves were obtained in 3 dose range 0 to 8 Gy, 0 to 10 Gy, and 0 to 30 Gy. The best-fitting model was assessed by the Akaike Information Criterion test. RESULTS: The ECF process resulted a useful tool to reduce the coefficients of variation from original values higher than 5% to about 3.5%, for all the batches exposed. The results confirm the linearity of dose-response curve below the dose level of 10 Gy for photon and electron beam and the supralinear trend above. CONCLUSION: The TLDs are suitable dosimeters for dose monitoring and verification in radiation treatment involving dose up to 30 Gy in a single fraction.

Modelling the risk of radiation induced alopecia in brain tumor patients treated with scanned proton beams.

PURPOSE: To develop normal tissue complication probability (NTCP) models for radiation-induced alopecia (RIA) in brain tumor patients treated with proton therapy (PT). METHODS AND MATERIALS: We analyzed 116 brain tumor adult patients undergoing scanning beam PT (median dose 54 GyRBE; range 36-72) for CTCAE v.4 grade 2 (G2) acute (≤90 days), late (>90 days) and permanent (>12 months) RIA. The relative dose-surface histogram (DSH) of the scalp was extracted and used for Lyman-Kutcher-Burman (LKB) modelling. Moreover, DSH metrics (Sx: the surface receiving ≥ X Gy, D2%: near maximum dose, Dmean: mean dose) and non-dosimetric variables were included in a multivariable logistic regression NTCP model. Model performances were evaluated by the cross-validated area under the receiver operator curve (ROC-AUC). RESULTS: Acute, late and permanent G2-RIA was observed in 52%, 35% and 19% of the patients, respectively. The LKB models showed a weak dose-surface effect (0.09 ≤ n ≤ 0.19) with relative steepness 0.29 ≤ m ≤ 0.56, and increasing tolerance dose values when moving from acute and late (22 and 24 GyRBE) to permanent RIA (44 GyRBE). Multivariable modelling selected S21Gy for acute and S25Gy, for late G2-RIA as the most predictive DSH factors. Younger age was selected as risk factor for acute G2-RIA while surgery as risk factor for late G2-RIA. D2% was the only variable selected for permanent G2-RIA. Both LKB and logistic models exhibited high predictive performances (ROC-AUCs range 0.86-0.90). CONCLUSION: We derived NTCP models to predict G2-RIA after PT, providing a comprehensive modelling framework for acute, late and permanent occurrences that, once externally validated, could be exploited for individualized scalp sparing treatment planning strategies in brain tumor patients.

Normal tissue complication probability (NTCP) models for modern radiation therapy.

Mathematical models of normal tissue complication probability (NTCP) able to robustly predict radiation-induced morbidities (RIM) play an essential role in the identification of a personalized optimal plan, and represent the key to maximizing the benefits of technological advances in radiation therapy (RT). Most modern RT techniques pose, however, new challenges in estimating the risk of RIM. The aim of this report is to schematically review NTCP models in the framework of advanced radiation therapy techniques. Issues relevant to hypofractionated stereotactic body RT and ion beam therapy are critically reviewed. Reirradiation scenarios for new or recurrent malignances and NTCP are also illustrated. A new phenomenological approach to predict RIM is suggested.

Technological evolution of radiation treatment: Implications for clinical applications.

The contemporary approach to the management of a cancer patient requires an "ab initio" involvement of different medical domains in order to correctly design an individual patient's pathway toward cure. With new therapeutic tools in every medical field developing faster than ever before the patient care outcomes can be achieved if all surgical, drug, and radiation options are considered in the design of the appropriate therapeutic strategy for a given patient. Radiation therapy (RT) is a clinical discipline in which experts from different fields continuously interact in order to manage the multistep process of the radiation treatment. RT is found to be an appropriate intervention for diverse indications in about 50% of cancer patients during the course of their disease. Technologies are essential in dealing with the complexity of RT treatments and for driving the increasingly sophisticated RT approaches becoming available for the treatment of Cancer. High conformal techniques, namely intensity modulated or volumetric modulated arc techniques, ablative techniques (Stereotactic Radiotherapy and Stereotactic Radiosurgery), particle therapy (proton or carbon ion therapy) allow for success in treating irregularly shaped or critically located targets and for the sharpness of the dose fall-off outside the target. The advanced on-board imaging, including real-time position management systems, makes possible image-guided radiation treatment that results in substantial margin reduction and, in select cases, implementation of an adaptive approach. The therapeutic gains of modern RT are also due in part to the enhanced anticancer activity obtained by coadministering RT with chemotherapy, targeted molecules, and currently immune checkpoints inhibitors. These main clinically relevant steps forward in Radiation Oncology represent a change of gear in the field that may have a profound impact on the management of cancer patients.