Reporting and analyzing dose distributions: a concept of equivalent uniform dose

Prospects for online adaptive radiation therapy (ART) for head and neck cancer

BACKGROUND

The aim of the present study is to examine the impact of kV-CBCT-based online adaptive radiation therapy (ART) on dosimetric parameters in comparison to image-guided-radiotherapy (IGRT) in consecutive patients with tumors in the head and neck region from a prospective registry.

METHODS

The study comprises all consecutive patients with tumors in the head and neck area who were treated with kV-CBCT-based online ART or IGRT-modus at the linear-accelerator ETHOS™. As a measure of effectiveness, the equivalent-uniform-dose was calculated for the CTV (EUDCTV) and organs-at-risk (EUDOAR) and normalized to the prescribed dose. As an important determinant for the need of ART the interfractional shifts of anatomic landmarks related to the tongue were analyzed and compared to the intrafractional shifts. The latter determine the performance of the adapted dose distribution on the verification CBCT2 postadaptation.

RESULTS

Altogether 59 consecutive patients with tumors in the head-and-neck-area were treated from 01.12.2021 to 31.01.2023. Ten of all 59 patients (10/59; 16.9%) received at least one phase within a treatment course with ART. Of 46 fractions in the adaptive mode, irradiation was conducted in 65.2% of fractions with the adaptive-plan, the scheduled-plan in the remaining. The dispersion of the distributions of EUDCTV-values from the 46 dose fractions differed significantly between the scheduled and adaptive plans (Ansari-Bradley-Test, p = 0.0158). Thus, the 2.5th percentile of the EUDCTV-values by the adaptive plans amounted 97.1% (95% CI 96.6-99.5%) and by the scheduled plans 78.1% (95% CI 61.8-88.7%). While the EUDCTV for the accumulated dose distributions stayed above 95% at PTV-margins of ≥ 3 mm for all 8 analyzed treatment phases the scheduled plans did for margins ≥ 5 mm. The intrafractional anatomic shifts of all 8 measured anatomic landmarks were smaller than the interfractional with overall median values of 8.5 mm and 5.5 mm (p < 0.0001 for five and p < 0.05 for all parameters, pairwise comparisons, signed-rank-test). The EUDOAR-values for the larynx and the parotid gland were significantly lower for the adaptive compared with the scheduled plans (Wilcoxon-test, p < 0.001).

CONCLUSIONS

The mobile tongue and tongue base showed considerable interfractional variations. While PTV-margins of 5 mm were sufficient for IGRT, ART showed the potential of decreasing PTV-margins and spare dose to the organs-at-risk.

Development of an automated treatment planning approach for lattice radiation therapy

BACKGROUND

Lattice radiation therapy (LRT) alternates regions of high and low doses within the target. The heterogeneous dose distribution is delivered to a geometrical structure of vertices segmented inside the tumor. LRT is typically used to treat patients with large tumor volumes with cytoreduction intent. Due to the geometric complexity of the target volume and the required dose distribution, LRT treatment planning demands additional resources, which may limit clinical integration.

PURPOSE

We introduce a fully automated method to (1) generate an ordered lattice of vertices with various sizes and center-to-center distances and (2) perform dose optimization and calculation. We aim to report the dosimetry associated with these lattices to help clinical decision-making.

METHODS

Sarcoma cancer patients with tumor volume between 100 cm3 and 1500 cm3 who received radiotherapy treatment between 2010 and 2018 at our institution were considered for inclusion. Automated segmentation and dose optimization/calculation were performed by using the Eclipse Scripting Application Programming Interface (ESAPI, v16, Varian Medical Systems, Palo Alto, USA). Vertices were modeled by spheres segmented within the gross tumor volume (GTV) with 1 cm/1.5 cm/2 cm diameters (LRT-1 cm/1.5 cm/2 cm) and 2 to 5 cm center-to-center distance on square lattices alternating along the superior-inferior direction. Organs at risk were modeled by subtracting the GTV from the body structure (body-GTV). The prescription dose was that 50% of the vertice volume should receive at least 20 Gy in one fraction. The automated dose optimization included three stages. The vertices optimization objectives were refined during optimization according to their values at the end of the first and second stages. Lattices were classified according to a score based on the minimization of body-GTV max dose and the maximization of GTV dose uniformity (measured with the equivalent uniform dose [EUD]), GTV dose heterogeneity (measured with the GTV D90%/D10% ratio), and the number of patients with more than one vertex inserted in the GTV. Plan complexity was measured with the modulation complexity score (MCS). Correlations were assessed with the Spearman correlation coefficient (r) and its associated p-value.

RESULTS

Thirty-three patients with GTV volumes between 150 and 1350 cm3 (median GTV volume = 494 cm3 , IQR = 272-779 cm3 were included. The median time required for segmentation/planning was 1 min/21 min. The number of vertices was strongly correlated with GTV volume in each LRT lattice for each center-to-center distance (r > 0.85, p-values < 0.001 in each case). Lattices with center-to-center distance = 2.5 cm/3 cm/3.5 cm in LRT-1.5 cm and center-to-center distance = 4 cm in LRT-1 cm had the best scores. These lattices were characterized by high heterogeneity (median GTV D90%/D10% between 0.06 and 0.19). The generated plans were moderately complex (median MCS ranged between 0.19 and 0.40).

CONCLUSIONS

The automated LRT planning method allows for the efficacious generation of vertices arranged in an ordered lattice and the refinement of planning objectives during dose optimization, enabling the systematic evaluation of LRT dosimetry from various lattice geometries.

The contest between internal and external-beam dosimetry: The Zeno's paradox of Achilles and the tortoise

Radionuclide therapy, also called molecular radiotherapy (MRT), has come of age, with several novel radiopharmaceuticals being approved for clinical use or under development in the last decade. External beam radiotherapy (EBRT) is a well-established treatment modality, with about half of all oncologic patients expected to receive at least one external radiation treatment over their disease course. The efficacy and the toxicity of both types of treatment rely on the interaction of radiation with biological tissues. Dosimetry played a fundamental role in the scientific and technological evolution of EBRT, and absorbed doses to the target and to the organs at risk are calculated on a routine basis. In contrast, in MRT the usefulness of internal dosimetry has long been questioned, and a structured path to include absorbed dose calculation is missing. However, following a similar route of development as EBRT, MRT treatments could probably be optimized in a significant proportion of patients, likely based on dosimetry and radiobiology. In the present paper we describe the differences and the similarities between internal and external-beam dosimetry in the context of radiation treatments, and we retrace the main stages of their development over the last decades.

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.

Review and recommendations on deformable image registration uncertainties for radiotherapy applications

Deformable image registration (DIR) is a versatile tool used in many applications in radiotherapy (RT). DIR algorithms have been implemented in many commercial treatment planning systems providing accessible and easy-to-use solutions. However, the geometric uncertainty of DIR can be large and difficult to quantify, resulting in barriers to clinical practice. Currently, there is no agreement in the RT community on how to quantify these uncertainties and determine thresholds that distinguish a good DIR result from a poor one. This review summarises the current literature on sources of DIR uncertainties and their impact on RT applications. Recommendations are provided on how to handle these uncertainties for patient-specific use, commissioning, and research. Recommendations are also provided for developers and vendors to help users to understand DIR uncertainties and make the application of DIR in RT safer and more reliable.

Impact of Modern Low Dose Involved Site Radiation Therapy on Normal Tissue Toxicity in Cervicothoracic Non-Hodgkin Lymphomas: A Biophysical Study

This biophysical study aimed to determine fitting parameters for the Lyman-Kutcher-Burman (LKB) dose-response model for normal tissue complication probability (NTCP) calculations of acute side effects and to investigate the impact of reduced radiation doses on the probability of their occurrence in supradiaphragmatic non-Hodgkin lymphoma (NHL) irradiation. A cohort of 114 patients with NHL in the cervicothoracic region, treated between 2015 and 2021 at the University Hospitals of Münster, Hamburg, and Essen, with involved site radiation therapy (ISRT), were included. Among them, 68 patients with aggressive NHL (a-NHL) received consolidative radiation therapy with 24-54 Gy following (R-)CHOP chemotherapy. Additionally, 46 patients with indolent NHL (i-NHL) underwent radiotherapy with 22.5-45.0 Gy. Two treatment plans were prospectively created for each patient (a-NHL: 30.0/40.0 Gy; i-NHL: 24.0/30.0 Gy). NTCP were then calculated using the optimized LKB model. The adapted dose-response models properly predicted the patient's probability of developing acute side effects when receiving doses ≤ 50 Gy. In addition, it was shown that reduced radiation doses can influence the NTCP of acute side effects depending on the aggressiveness of NHL significantly. This study provided a foundation to prospectively assess the probability of adverse side effects among today's reduced radiation doses in the treatment of NHL.

Stereotactic Radiosurgery of Multiple Brain Metastases: A Review of Treatment Techniques

The advancement of systemic targeted treatments has led to improvements in the management of metastatic disease, particularly in terms of survival outcomes. However, brain metastases remain less responsive to systemic therapies, underscoring the significance of local interventions for comprehensive disease control. Over the past years, the threshold for treating brain metastases through stereotactic radiosurgery has risen. Yet, as the number of treated metastases increases, treatment complexity and duration also escalate. This trend has made multi-isocenter radiosurgery treatments, such as those with the Gamma Knife, challenging to plan and lengthy for patients. In contrast, single-isocenter approaches employing linear accelerators offer an efficient and expeditious treatment option. This review delves into the literature, comparing different linear-accelerator-based techniques with each other and in relation to dedicated systems, focusing on dosimetric considerations and feasibility.

Accounting for fractionation and heterogeneous dose distributions in the modelling of osteoradionecrosis in oropharyngeal carcinoma treatment

BACKGROUND AND PURPOSE

Osteoradionecrosis (ORN) of the mandible is a severe complication following radiotherapy (RT). With a renewed interest in hypofractionation for head and neck radiotherapy, more information concerning ORN development after high fraction doses is important. The aim of this explorative study was to develop a model for ORN risk prediction applicable across different fractionation schemes using Equivalent Uniform Doses (EUD).

MATERIAL AND METHODS

We performed a retrospective cohort study in 334 oropharyngeal squamous cell carcinoma (OPSCC) patients treated with either a hypofractionated Stereotactic Body Radiation Therapy (HF-SBRT) boost or conventional Intensity Modulated Radiation Therapy (IMRT). ORN was scored with the CTCAE v5.0. HF-SBRT and IMRT dose distributions were converted into equivalent dose in 2 Gy fractions (α/β = 0.85 Gy) and analyzed using EUD. The parameter a that led to an EUD that best discriminated patients with and without grade ≥ 2 ORN was selected. Patient and treatment-related risk factors of ORN were analyzed with uni- and multivariable regression analysis.

RESULTS

A total of 32 patients (9.6%) developed ORN grade ≥ 2. An EUD(a = 8) best discriminated between ORN and non-ORN (AUC = 0.71). In multivariable regression, pre-RT extractions (SHR = 2.34; p = 0.012), mandibular volume (SHR = 1.04; p = 0.003), and the EUD(a = 8) (SHR = 1.14; p < 0.001) were significantly associated with ORN.

CONCLUSION

Risk models for ORN based on conventional DVH parameters cannot be directly applied to HF-SBRT fractionation schemes and dose distributions. However, after correcting for fractionation and non-uniform dose distributions using EUD, a single model can distinguish between ORN and non-ORN after conventionally fractionated radiotherapy and hypofractionated boost treatments.

Model development of dose and volume predictors for esophagitis induced during chemoradiotherapy for lung cancer as a step towards radiobiological treatment planning

BACKGROUND

Currently, radiation therapy treatment planning system intends biological optimization that relies heavily upon plan metrics from tumor control probability (TCP) and normal tissue complication probability (NTCP) modeling. Implementation and expansion of TCP and NTCP models with alternative data is an important step towards reliable radiobiological treatment planning. In this retrospective single institution study, the treatment charts of 139 lung cancer patients treated with chemo-radiotherapy were reviewed and correlated dosimetric predictors with the incidence of esophagitis and established NTCP model of esophagitis grade 1 and 2 for lung cancer patients.

METHODS

Esophagus is an organ at risk (OAR) in lung cancer radiotherapy (RT). Esophagitis is a common toxicity induced by RT. In this study, dose volume parameters Vx (Vx: percentage esophageal volume receiving ≥ x Gy) and mean esophagus dose (MED) as quantitative dose-volume metrics, the esophagitis grade 1 and 2 as endpoints, were reviewed and derived from the treatment planning system and the electronic medical record system. Statistical analysis of binary logistic regression and probit were performed to have correlated the probability of grade 1 and 2 esophagitis to MED and Vx. IBM SPSS software version 24 at 5% significant level (α = 0.05) was used in the statistical analysis.

RESULTS

The probabilities of incidence of grade 1 and 2 esophagitis proportionally increased with increasing the values of Vx and MED. V20, V30, V40, V50 and MED are statistically significant good dosimetric predictors of esophagitis grade 1. 50% incidence probability (TD50) of MED for grade 1 and 2 esophagitis were determined. Lyman Kutcher Burman model parameters, such as, n, m and TD50, were fitted and compared with other published findings. Furthermore, the sigmoid shaped dose responding curve between probability of esophagitis grade 1 and MED were generated respecting to races, gender, age and smoking status.

CONCLUSIONS

V20, V30, V40 and V50 were added onto Quantitative Analysis of Normal Tissue Effects in the clinic, or QUANTEC group's dose constrains of V35, V50, V70 and MED. Our findings may be useful as both validation of 3-Dimensional planning era models and also additional clinical guidelines in treatment planning and plan evaluation using radiobiology optimization.

Selective sparing of bladder and rectum sub-regions in radiotherapy of prostate cancer combining knowledge-based automatic planning and multicriteria optimization

Background and Purpose

The association between dose to selected bladder and rectum symptom-related sub-regions (SRS) and late toxicity after prostate cancer radiotherapy has been evidenced by voxel-wise analyses. The aim of the current study was to explore the feasibility of combining knowledge-based (KB) and multi-criteria optimization (MCO) to spare SRSs without compromising planning target volume (PTV) dose delivery, including pelvic-node irradiation.

Materials and Methods

Forty-five previously treated patients (74.2 Gy/28fr) were selected and SRSs (in the bladder, associated with late dysuria/hematuria/retention; in the rectum, associated with bleeding) were generated using deformable registration. A KB model was used to obtain clinically suitable plans (KB-plan). KB-plans were further optimized using MCO, aiming to reduce dose to the SRSs while safeguarding target dose coverage, homogeneity and avoiding worsening dose volume histograms of the whole bladder, rectum and other organs at risk. The resulting MCO-generated plans were examined to identify the best-compromise plan (KB + MCO-plan).

Results

The mean SRS dose decreased in almost all patients for each SRS. D1% also decreased in the large majority, less frequently for dysuria/bleeding SRS. Mean differences were statistically significant (p < 0.05) and ranged between 1.3 and 2.2 Gy with maximum reduction of mean dose up to 3-5 Gy for the four SRSs. The better sparing of SRSs was obtained without compromising PTVs coverage.

Conclusions

Selectively sparing SRSs without compromising PTV coverage is feasible and has the potential to reduce toxicities in prostate cancer radiotherapy. Further investigation to better quantify the expected risk reduction of late toxicities is warranted.

Dosimetric and radiobiological comparison between conventional and hypofractionated breast treatment plans using the Halcyon system

Purpose

The objective of this research is to compare the efficacy of conventional and hypofractionated radiotherapy treatment plans for breast cancer patients, with a specific focus on the unique features of the Halcyon system.

Methods and materials

The study collected and analyzed dose volume histogram (DVH) data for two groups of treatment plans implemented using the Halcyon system. The first group consisted of 19 patients who received conventional fractionated (CF) treatment with a total dose of 50 Gy in 25 fractions, while the second group comprised 9 patients who received hypofractionated (HF) treatment with a total dose of 42.56 Gy in 16 fractions. The DVH data was used to calculate various parameters, including tumor control probability (TCP), normal tissue complication probability (NTCP), and equivalent uniform dose (EUD), using radiobiological models.

Results

The results indicated that the CF plan resulted in higher TCP but lower NTCP for the lungs compared to the HF plan. The EUD for the HF plan was approximately 49 Gy (114% of its total dose) while that for the CF plan was around 53 Gy (107% of its total dose).

Conclusions

The analysis suggests that while the CF plan is better at controlling tumors, it is not as effective as the HF plan in minimizing side effects. Additionally, it is suggested that there may be an optimal configuration for the HF plan that can provide the same or higher EUD than the CF plan.

Comparing Predicted Toxicities between Hypofractionated Proton and Photon Radiotherapy of Liver Cancer Patients with Different Adaptive Schemes

With the availability of MRI linacs, online adaptive intensity modulated radiotherapy (IMRT) has become a treatment option for liver cancer patients, often combined with hypofractionation. Intensity modulated proton therapy (IMPT) has the potential to reduce the dose to healthy tissue, but it is particularly sensitive to changes in the beam path and might therefore benefit from online adaptation. This study compares the normal tissue complication probabilities (NTCPs) for liver and duodenal toxicity for adaptive and non-adaptive IMRT and IMPT treatments of liver cancer patients. Adaptive and non-adaptive IMRT and IMPT plans were optimized to 50 Gy (RBE = 1.1 for IMPT) in five fractions for 10 liver cancer patients, using the original MRI linac images and physician-drawn structures. Three liver NTCP models were used to predict radiation-induced liver disease, an increase in albumin-bilirubin level, and a Child-Pugh score increase of more than 2. Additionally, three duodenal NTCP models were used to predict gastric bleeding, gastrointestinal (GI) toxicity with grades >3, and duodenal toxicity grades 2-4. NTCPs were calculated for adaptive and non-adaptive IMRT and IMPT treatments. In general, IMRT showed higher NTCP values than IMPT and the differences were often significant. However, the differences between adaptive and non-adaptive treatment schemes were not significant, indicating that the NTCP benefit of adaptive treatment regimens is expected to be smaller than the expected difference between IMRT and IMPT.

Comparing biological and physical cost functions in VMAT planning for pediatric nasopharyngeal cancer

BACKGROUND

Volumetric Modulated Arc Therapy (VMAT) is an option for the delivery of Radiotherapy treatment technique for pediatric nasopharyngeal cancer, VMAT is the most common treatment technique for pediatric nasopharyngeal cancer. The use of a combination of both biological and physical parameters in VMAT planning optimization may produce better target coverage and sparing of critical organs. This work was to compare Biological Cost Functions (BCFs) and Physical Cost Functions (PCFs) in the VMAT of pediatric nasopharyngeal cancer patients.

METHOD

VMAT plans for 20 nasopharyngeal pediatric cancer patients were created using Monaco 5.11® treatment planning system (TPS). Three VMAT plans were retrospectively generated for each patient using BCFs, PCFs and mixed plan with a total dose of 61.2 Gy in 34 fractions to planning target volume (PTV). All plans were adjusted to deliver 95% of the prescribed dose to 95% of the PTV. The calculated plans were qualitatively and quantitatively evaluated using the dose-volume histogram (DVH).

RESULTS

The coverage of the target and the maximum dose for the three plans were nearly the same, and better sparing was achieved in the serial organs (spinal cord and brain stem) with PCFs. On the contrary, more dose spring was observed using the BCFs in the organs at risk (OARs) that were not involved in the dose optimization, such as the optic nerve maximum dose, with a significant p-value (0.035 and 0.0001) respectively. Using the PCFs, both parotids received a lower mean dose, but not for the oral cavity, which had a lower mean dose using BCFs (p=<0.0001). The same values of tumor control probability (TCP) were found for both cost functions in PTVs and normal tissue complications probability (NTCP) (99%). The values reported were as follows: spinal cord = 0.5%, brain stem = 19.1%, and brain = 90.7% for BCFs, compared to spinal cord = 0.3%, brain stem = 14.9%, and brain = 90.7% for PCFs. The delivery time was found to be less in BCFs (p=0.005).

CONCLUSION

The BCFs are superior to the PCFs in conformity index and time of radiation delivery. However, PCFs were better at dose sparing for the serial organs and achieving a sharper falloff dose around the involved volumes. A patient-specific clinical compromise is recommended to gain the best plan that meets the clinical goals.

NTCP Calculations of Five Different Irradiation Techniques for the Treatment of Thymoma

This study provided normal tissue complication probability (NTCP) calculations from photon radiotherapy techniques in eleven patients with thymoma. Five plans were created for each participant using three-dimensional conformal radiotherapy (3D-CRT), five-field intensity modulated radiotherapy (5F-IMRT), seven-field IMRT (7F-IMRT), and volumetric modulated arc therapy with full arcs (FA-VMAT) and partial arcs (PA-VMAT). The target coverage, homogeneity index and conformation number for the planning target volume (PTV) and dosimetric parameters for the organs-at-risk (OARs) were taken from the fifty-five generated plans. The patient-specific NTCP of the lungs, heart and esophagus was calculated with an in-house software tool using differential dose-volume histograms and the equivalent uniform dose model. The PTV dose metrics from 3D-CRT were inferior to those from IMRT and VMAT plans. The dose constraints for the OARs were met in all treatment plans. The NTCP range of the lungs, heart and esophagus was 0.34-0.49%, 0.03-0.06% and 0.08-0.10%, respectively. The NTCPs of the heart for the incidence of peridarditis from IMRT and VMAT were significantly smaller than those from conformal treatment (p < 0.05). The 7F-IMRT was significantly superior to FA-VMAT in reducing the NTCP of the lungs and the risk of pneumonitis (p = 0.001). Similar superiority of 5F-IMRT over PA-VMAT for lung protection was found (p = 0.009). The presented results may be employed in the selection of the appropriate irradiation technique for restricting the complications in the adjacent OARs.

Influence on voxel-based dosimetry: noise effect on absorbed dose dosimetry at single time-point versus sequential single-photon emission computed tomography

OBJECTIVE

The purpose of this study was to evaluate how statistical fluctuation in single-photon emission computed tomography (SPECT) images propagate to absorbed dose maps.

METHODS

SPECT/computed tomography (CT) images of iodine-131 filled phantoms, using different acquisition and processing protocols, were evaluated using STRATOS software to assess the absorbed dose distribution at the voxel level. Absorbed dose values and coefficient of variation (COV) were analyzed for dosimetry based on single time-point SPECT images and time-integrated activities of SPECT sequences with low and high counts.

RESULTS

Considering dosimetry based on a single time-point, the mean absorbed dose was not significantly affected by total counts or reconstruction parameters, but the uniformity of the absorbed dose maps had an almost linear correlation with SPECT noise. When high- and low-count SPECT sequences were used to generate an absorbed dose map, the absorbed dose COV for each of the temporal sequences was slightly lower than the absorbed dose COV based on the single SPECT image with the highest count included in the sequence.

CONCLUSION

The impact of changes in SPECT counts and reconstruction parameters is almost linear when dosimetry is based on isolated SPECT images, but less pronounced when dosimetry is based on sequential SPECTs.

Built-in wavelet-induced smoothness to reduce plan complexity in intensity modulated radiation therapy (IMRT)

Objective.Reducing plan complexity in intensity modulated radiation therapy (IMRT) to ensure dosimetric accuracy and delivery efficiency of the radiation treatment plans. We propose a novel approach by representing the beamlet intensities using an incomplete wavelet basis that explicitly excludes fluctuating intensity maps from the decision space (explicit hard constraint). This technique provides a built-in wavelet-induced smoothness that improves both dosimetric plan quality and delivery efficiency.Approach.The beamlet intensity maps need to be especially smooth in the leaf travel direction (referred to as theX-direction). We treat the intensity map of each beam as a 2D image and represent it using the wavelets corresponding to low-frequency changes in theX-direction (i.e. approximation and horizontal). The absence of wavelets corresponding to high-frequency changes (i.e. vertical and diagonal) induces built-in smoothness. We still utilize a regularization term in the objective function to promote smoothness in theY-direction (perpendicular to theX-direction) and further possible smoothness in theX-direction. This technique has been tested on three patient cases of different disease sites (paraspinal, lung, prostate) and all final evaluations and comparisons have been performed on an FDA-approved commercial treatment planning system (Varian Eclipse^TM).Main results.Wavelet-induced smoothness reduced monitor units by about 10%, 45%, and 14% for paraspinal, lung, and prostate cases, respectively. It also improved organ at risk sparing, especially on the complex paraspinal case where it resulted in about 7%, 13%, and 14% less mean dose to esophagus, lung, and cord, respectively. Moreover, built-in wavelet-induced smoothness desensitizes the results to changing the weight associated to the regularization term, and thereby mitigates the weight fine-tuning difficulty.Significance.Fluence smoothness is often achieved by smoothing the beamlet intensity maps using a proper regularization term in the objective function aiming at disincentivizing fluctuation in the beamlet intensities (implicit soft constraint). This work reports a novel application of wavelets in imposing an explicit smoothness hard constraint in the search space using an incomplete wavelet basis. This idea has been successfully applied to exclude complex and clinically irrelevant radiation plans from the search space. The code and pertained models along with a sample dataset are released on our LowDimRT GitHub (https://github.com/PortPy-Project/LowDimRT).

Parameters of the Lyman Model for Calculation of Normal-Tissue Complication Probability: A Systematic Literature Review

PURPOSE

The Lyman model is one of the most used radiobiological models for calculation of normal-tissue complication probability (NTCP). Since its introduction in 1985, many authors have published parameter values for the model based on clinical data of different radiotherapeutic situations. This study attempted to collect the entirety of radiobiological parameter sets published to date and provide an overview of the data basis for different variations of the model. Furthermore, it sought to compare the parameter values and calculated NTCPs for selected endpoints with sufficient data available.

METHODS AND MATERIALS

A systematic literature analysis was performed, searching for publications that provided parameters for the different variations of the Lyman model in the Medline database using PubMed. Parameter sets were grouped into 13 toxicity-related endpoint groups. For 3 selected endpoint groups (≤25% reduction of saliva 12 months after irradiation of the parotid, symptomatic pneumonitis after irradiation of the lung, and bleeding of grade 2 or less after irradiation of the rectum), parameter values were compared and differences in calculated NTCP values were analyzed.

RESULTS

A total of 509 parameter sets from 130 publications were identified. Considerable heterogeneities were detected regarding the number of parameters available for different radio-oncological situations. Furthermore, for the 3 selected endpoints, large differences in published parameter values were found. These translated into great variations of calculated NTCPs, with maximum ranges of 35.2% to 93.4% for the saliva endpoint, of 39.4% to 90.4% for the pneumonitis endpoint, and of 5.4% to 99.3% for the rectal bleeding endpoint.

CONCLUSIONS

The detected heterogeneity of the data as well as the large variations of published radiobiological parameters underline the necessity for careful interpretation when using such parameters for NTCP calculations. Appropriate selection of parameters and validation of values are essential when using the Lyman model.

Effective optimization strategy for large optimization volume object, remaining volume at risk (RVR):α-value selection and usage from generalized equivalent uniform dose (gEUD) curve deviation perspective

Objective.A large optimization volume for intensity-modulated radiation therapy (IMRT), such as the remaining volume at risk (RVR), is traditionally unsuitable for dose-volume constraint control and requires planner-specific empirical considerations owing to the patient-specific shape. To enable less empirical optimization, the generalized equivalent uniform dose (gEUD) optimization is effective; however, the utilization of parametera-values remains elusive. Our study clarifies thea-value characteristics for optimization and to enable effectivea-value use.Approach.The gEUD can be obtained as a function of itsa-value, which is the weighted generalized mean; its curve has a continuous, differentiable, and sigmoid shape, deforming in its optimization state with retained curve characteristics. Using differential geometry, the gEUD curve changes in optimization is considered a geodesic deviation intervened by the forces between deforming and retaining the curve. The curvature and gradient of the curve are radically related to optimization. The vertex point (a=a_k) was set and thea-value roles were classified into the following three parts of the curve with respect to thea-value: (i) high gradient and middle curvature, (ii) middle gradient and high curvature, and (iii) low gradient and low curvature. Then, a strategy for multiplea-values was then identified using RVR optimization.Main results.Eleven head and neck patients who underwent static seven-field IMRT were used to verify thea-value characteristics and curvature effect for optimization. The lowera-value (i) (a= 1-3) optimization was effective for the whole dose-volume range; in contrast, the effect of highera-value (iii) (a= 12-20) optimization addressed strongly the high-dose range of the dose volume. The middlea-value (ii) (arounda=a_k) showed intermediate but effective high-to-low dose reduction. Thesea-value characteristics were observed as superimpositions in the optimization. Thus, multiple gEUD-based optimization was significantly superior to the exponential constraints normally applied to the RVR that surrounds the PTV, normal tissue objective (NTO), resulting in up to 25.9% and 8.1% improvement in dose-volume indices D2% and V10Gy, respectively.Significance.This study revealed an appropriatea-value for gEUD optimization, leading to favorable dose-volume optimization for the RVR region using fixed multiplea-value conditions, despite the very large and patient-specific shape of the region.

Automated Contouring and Planning in Radiation Therapy: What Is 'Clinically Acceptable'?

Developers and users of artificial-intelligence-based tools for automatic contouring and treatment planning in radiotherapy are expected to assess clinical acceptability of these tools. However, what is 'clinical acceptability'? Quantitative and qualitative approaches have been used to assess this ill-defined concept, all of which have advantages and disadvantages or limitations. The approach chosen may depend on the goal of the study as well as on available resources. In this paper, we discuss various aspects of 'clinical acceptability' and how they can move us toward a standard for defining clinical acceptability of new autocontouring and planning tools.

An integrated solution of deep reinforcement learning for automatic IMRT treatment planning in non-small-cell lung cancer

Purpose

To develop and evaluate an integrated solution for automatic intensity-modulated radiation therapy (IMRT) planning in non-small-cell lung cancer (NSCLC) cases.

Methods

A novel algorithm named as multi-objectives adjustment policy network (MOAPN) was proposed and trained to learn how to adjust multiple optimization objectives in commercial Eclipse treatment planning system (TPS), based on the multi-agent deep reinforcement learning (DRL) scheme. Furthermore, a three-dimensional (3D) dose prediction module was developed to generate the patient-specific initial optimization objectives to reduce the overall exploration space during MOAPN training. 114 previously treated NSCLC cases suitable for stereotactic body radiotherapy (SBRT) were selected from the clinical database. 87 cases were used for the model training, and the remaining 27 cases for evaluating the feasibility and effectiveness of MOAPN in automatic treatment planning.

Results

For all tested cases, the average number of adjustment steps was 21 ± 5.9 (mean ± 1 standard deviation). Compared with the MOAPN initial plans, the actual dose of chest wall, spinal cord, heart, lung (affected side), esophagus and bronchus in the MOAPN final plans reduced by 14.5%, 11.6%, 4.7%, 16.7%, 1.6% and 7.7%, respectively. The dose result of OARs in the MOAPN final plans was similar to those in the clinical plans. The complete automatic treatment plan for a new case was generated based on the integrated solution, with about 5-6 min.

Conclusion

We successfully developed an integrated solution for automatic treatment planning. Using the 3D dose prediction module to obtain the patient-specific optimization objectives, MOAPN formed action-value policy can simultaneously adjust multiple objectives to obtain a high-quality plan in a shorter time. This integrated solution contributes to improving the efficiency of the overall planning workflow and reducing the variation of plan quality in different regions and treatment centers. Although improvement is warranted, this proof-of-concept study has demonstrated the feasibility of this integrated solution in automatic treatment planning based on the Eclipse TPS.

Comparison of intensity-modulated radiation therapy (IMRT), 3D conformal proton therapy and intensity-modulated proton therapy (IMPT) for the treatment of metastatic brain cancer

The purpose of this study has been to compare photon intensity modulated radiation therapy (IMRT) against both conformal and intensity modulated proton therapy (IMPT) plans for metastatic brain cancer. Ten IMRT patients with brain cancer were chosen retrospectively, with prescription doses in the range of 20 to 40 Gy, delivered in 3 to 5 fractions using Varian TrueBeam STx machine. Three proton plans with proton double scattering, single collimation static-IMPT, and energy layer by layer collimation dynamic-IMPT were then generated for the same patients using the Mevion S250 system for conformal plans and the S250i system for IMPT plans. Each plan had respective treatment planning systems that include Brainlab iPlan for IMRT, Varian Eclipse for proton double scattering, and RaySearch Raystation for IMPT. Dosimetric and radiobiologic comparisons were made through dose-volume histogram (DVH) analysis of the target and the organs at risk (OAR); and with parameters of equivalent uniform dose (EUD), tumor control probability (TCP), and normal tissue complication probability (NTCP), respectively. A set of variables α/β ratio, survival fraction, and clonogenic cell density were selected and varied to observe their effect on the abovementioned parameters. Doses were observed to be more homogeneous for patients with brain malignancies with photon IMRT treatments, while dose conformity is improved with proton PBS treatments. Normal tissue is, on average, spared more through both proton treatment options. The minimum doses, closely approximated by dose to 98% of the target volume, are similar across treatment modalities with slight variations.

Machine-learning-based prediction of the effectiveness of the delivered dose by exhale-gated radiotherapy for locally advanced lung cancer: The additional value of geometric over dosimetric parameters alone

Purpose

This study aimed to assess interfraction stability of the delivered dose distribution by exhale-gated volumetric modulated arc therapy (VMAT) or intensity-modulated arc therapy (IMAT) for lung cancer and to determine dominant prognostic dosimetric and geometric factors.

Methods

Clinical target volume (CTV_Plan) from the planning CT was deformed to the exhale-gated daily CBCT scans to determine CTV_i, treated by the respective dose fraction. The equivalent uniform dose of the CTV_i was determined by the power law (gEUD_i) and cell survival model (EUD_iSF) as effectiveness measure for the delivered dose distribution. The following prognostic factors were analyzed: (I) minimum dose within the CTV_i (D_{min_i}), (II) Hausdorff distance (HDD_i) between CTV_i and CTV_Plan, (III) doses and deformations at the point in CTV_Plan at which the global minimum dose over all fractions per patient occurs (PD_{min_global_i}), and (IV) deformations at the point over all CTV_i margins per patient with the largest Hausdorff distance (HDPw_orst). Prognostic value and generalizability of the prognostic factors were examined using cross-validated random forest or multilayer perceptron neural network (MLP) classifiers. Dose accumulation was performed using back deformation of the dose distribution from CTV_i to CTV_Plan.

Results

Altogether, 218 dose fractions (10 patients) were evaluated. There was a significant interpatient heterogeneity between the distributions of the normalized gEUD_i values (p<0.0001, Kruskal-Wallis tests). Accumulated gEUD over all fractions per patient was 1.004-1.023 times of the prescribed dose. Accumulation led to tolerance of ~20% of fractions with gEUD_i <93% of the prescribed dose. Normalized D_min >60% was associated with predicted gEUD values above 95%. D_min had the highest importance for predicting the gEUD over all analyzed prognostic parameters by out-of-bag loss reduction using the random forest procedure. Cross-validated random forest classifier based on D_min as the sole input had the largest Pearson correlation coefficient (R=0.897) in comparison to classifiers using additional input variables. The neural network performed better than the random forest classifier, and the gEUD values predicted by the MLP classifier with D_min as the sole input were correlated with the gEUD values characterized by R=0.933 (95% CI, 0.913-0.948). The performance of the full MLP model with all geometric input parameters was slightly better (R=0.952) than that based on D_min (p=0.0034, Z-test).

Conclusion

Accumulated dose distributions over the treatment series were robust against interfraction CTV deformations using exhale gating and online image guidance. D_min was the most important parameter for gEUD prediction for a single fraction. All other parameters did not lead to a markedly improved generalizable prediction. Dosimetric information, especially location and value of D_min within the CTV _i , are vital information for image-guided radiation treatment.

Comparison of intensity-modulated proton therapy (IMPT) versus intensity-modulated radiation therapy (IMRT) for the treatment of head and neck cancer based on radiobiological modelling

Aim:

The aim of our study is to retrospectively report the radiobiological aspects for intensity-modulated proton therapy (IMPT) against intensity-modulated radiation therapy (IMRT) for patients with head and neck cancer treated at our institution. A secondary goal is to reinforce current model-based approaches to head and neck cancer patient selection for IMPT.

Materials and Methods:

Eighteen patients were evaluated with prescription doses ranging from 50 to 70 Gy delivered in 2 Gy per fraction. The dose volume histograms (DVH) were used to calculate equivalent uniform dose (EUD), tumour control probability (TCP) and normal tissue complication probability (NTCP) for biophysical comparison using mechanistic mathematical dose response models. Absolute values of TCP and NTCP were then compared between IMPT and IMRT.

Results:

The dose models demonstrate a minimal radiobiological advantage for IMPT compared to IMRT in treating head and neck cancers. Absolute values of TCP were slightly higher, while absolute values of NTCP were slightly lower for IMPT versus IMRT.

Conclusions:

Further studies are needed to determine if the radiobiological advantage indeed translates to a therapeutic advantage for patients.

Validation of esophageal cancer treatment methods from 3D-CRT, IMRT, and Rapid Arc plans using custom Python software to compare radiobiological plans to normal tissue integral dosage

Background

The aim was to develop in-house software that is able to calculate and generate the biological plan evaluation of the esophagus treatment plan using the Niemierko model for normal tissue complication probability and tumor control probability. The Niemierko model can be applied for esophagus cancer treatment plan to estimate the tumor control probability (TCP) and the normal tissue complication probability (NTCP) using different planning techniques. The equivalent uniform dose (EUD) and effective volume parameters were compared with organ at risk. Subsequently, EUD and TCP parameter were compared with tumor volume for all five different planning techniques.

Materials and methods

Ten cases for esophageal cancer were included in this study. For each patient, five treatment plans were generated. The Anisotropic analytical algorithms (AAA) were used for dose calculation for the three-dimensional conformal radiation therapy (3D-CRT), intensity modulated radiation therapy (IMRT) and volumetric modulated arc therapy (VMAT) techniques. The in-house developed radiobiological plan evaluation software using python programming is used for this study which takes a dose volume histogram (DVH) text file as an input file for biological plan evaluation.

Results and Conclusion

EUD, NTCP, TCP and effective volume were calculated from the Niemierko model using the in-house developed python based software and compared with treatment monitor units (MU) with all five different treatment plan. The best technique is quantified as benchmarked out of other different qualities of treatment. The four field 3D-CRT treatment plan is found to be the best suited from the perspective of biological plan index evaluation among the other planning techniques.

Escalation of radiotherapy dose in large locally advanced drug-resistant gastrointestinal stromal tumors by multi-shell simultaneous integrated boost intensity-modulated technique: a feasibility study

BACKGROUND

Resistance to conventional dose schemes and radiotoxicity of healthy tissue is a clinical challenge in the radiation therapy of large locally advanced drug-resistant gastrointestinal stromal tumor (LADR-GIST). This study aimed to assess the feasibility of using multi-shell Simultaneous Integrated Boost Intensity-Modulated modality (SIB-IMRT) strategy to provide a safe and effective escalation dose regimen for LADR-GIST.

METHODS

7 patients with LADR-GIST were selected in this study. The modified SIB-IMRT plans for all patients were generated by delivering different escalation-dose gradients to four ring shaped regions (shells) within the gross tumor volume (GTV). The doses of the central volume of the tumor (GTVcenter) were escalated up to 70-92.5 Gy (25 fractions), while the doses of planning target volume (PTV) and shell-1 were kept at 50.0 Gy. Based on different escalation-dose gradients, the modified SIB-IMRT plans were divided into four groups (SIB-IMRT groups). For comparison purposes, plans obtained by conventional IMRT technique (Con-IMRT) with 50 Gy (25 fractions) were also generated for all patients (Con-IMRT group). All plans were normalized to cover 95% of the PTV with the prescribed dose of 50.0 Gy. The equivalent uniform dose (EUD), relative equivalent uniform dose (rEUD), dose volume histogram (DVH), dose profile, conformity index (CI) and monitor unit (MU) were evaluated in five groups. The Friedman Test was performed to determine whether there were significant differences (P < 0.05).

RESULTS

Compared with the Con-IMRT group, the EUD of GTV (EUDGTV) and rEUD of SIB-IMRT groups were improved when escalation-dose gradient was increased, and the improvement became significant when the escalation-dose gradient reached 20% of the prescription dose. The rEUD tended to be stable as the escalation-dose gradient went up to 25% of the prescription dose. There were no significant differences in CIs and DVH metrics for OARs between the Con-IMRT group and any SIB-IMRT group, but the significant differences were observed between the SIB10-IMRT group and the SIB25-IMRT group. For the SIB-IMRT groups, as the dose gradient became steeper in the dose profiles, the higher dose was mainly accumulated in the inner part of GTV accompanied with a higher MU.

CONCLUSIONS

The proposed multi-shell SIB-IMRT strategy is feasible in dosimetry for LADR-GIST and can acquire higher therapeutic gain without sacrifice of healthy tissues. It appears that the scheme of delivering 20% of the prescribed escalation-dose gradient to the target volume can provide satisfactory dose irradiation for LADR-GIST, and it should be evaluated in future clinical study.

Incidental nodal irradiation in patients with esophageal cancer undergoing (chemo)radiation with 3D-CRT or VMAT

The extent of elective nodal irradiation (ENI) in patients undergoing definitive chemoradiotherapy (dCRT) for esophageal squamous cell carcinoma (ESCC) remains unclear. The aim of this dosimetric study was to evaluate the extent of incidental nodal irradiation using modern radiation techniques. A planning target volume (PTV) was generated for 30 patients with node-negative esophageal carcinoma (13 cervical/upper third, 7 middle third, 10 lower third/abdomen). Thereby, no elective nodal irradiation (ENI) was intended. Both three-dimensional conformal radiotherapy (3D-CRT) and volumetric-modulated arc therapy (VMAT) treatment plans (50 Gy in 25 fractions) were calculated for all patients. Fifteen nodal stations were contoured according to the definitions of the AJCC and investigated in regard to dosimetric parameters. Compared to 3D-CRT, VMAT was associated with lower dose distribution to the organs at risk (lower Dmean, V20 and V30 for the lungs and lower Dmean and V30 for the heart). For both techniques, the median Dmean surpassed 40 Gy in 12 of 15 (80%) nodal stations. However, VMAT resulted in significantly lower Dmeans and equivalent uniform doses (EUD) compared to 3D-CRT for eight nodal stations (1L, 2L, 2R, 4L, 7, 8L, 10L, 15), while differences did not reach significance for seven nodal station (1R, 4R, 8U, 8M, 10R, 16). For dCRT of ESCC, the use of VMAT was associated with significantly lower median (incidental) doses to eight of 15 regional lymph node areas compared to 3D-CRT. However, given the small absolute differences, these differences probably do not impair (regional) tumor control rates.

Radiobiological and dosimetric comparison of 60Co versus 192Ir high-dose-rate intracavitary-interstitial brachytherapy for cervical cancer

BACKGROUND

High-dose-rate (HDR) intracavitary-interstitial brachytherapy (IC-ISBT) is an effective treatment for bulky, middle, and advanced cervical cancer. In this study, we compared the differences between 60Co and 192Ir HDR IC-ISBT plans in terms of radiobiological and dosimetric parameters, providing a reference for clinical workers in brachytherapy.

METHODS

A total of 30 patients with cervical cancer receiving HDR IC-ISBT were included in this study, and IC-ISBT plans for each individual were designed with both 60Co and 192Ir at a prescribed dose of CTV D90 = 6 Gy while keeping the dose to OARs as low as possible. Physical dose and dose-volume parameters of CTV and OARs were extracted from TPS. The EQD2, EUBED, EUD, TCP, and NTCP were calculated using corresponding formulas. The differences between the 60Co and 192Ir IC-ISBT plans were compared using the paired t-test.

RESULTS

In each patient's 60Co and 192Ir IC-ISBT plan, the average physical dose and EQD2 of 60Co were lower than those of 192Ir, and there were statistically significant differences in D2cc and D1cc for the OARs (p < 0.05); there were statistically significant differences in D0.1 cc for the bladder (p < 0.05) and no significant differences in D0.1 cc for the rectum or intestines (p > 0.05). The EUBED ratio (60Co/192Ir) at the CTV was mostly close to 1 when neither 60Co or 192Ir passed their half-lives or when both passed two half-lives, and the difference between them was not significant; at the OARs, the mean value of 60Co was lower than that of 192Ir. There was no statistical difference between 60Co and 192Ir in the EUD (93.93 versus 93.92 Gy, p > 0.05) and TCP (97.07% versus 97.08%, p > 0.05) of the tumors. The mean NTCP value of 60Co was lower than that of 192Ir.

CONCLUSIONS

Considering the CTV and OARs, the dosimetric parameters of 60Co and 192Ir are comparable. Compared with 192Ir, the use of 60Co for HDR IC-ISBT can ensure a similar tumor control probability while providing better protection to the OARs. In addition, 60Co has obvious economic advantages and can be promoted as a good alternative to 192Ir.

Automatic Radiobiological Comparison of Radiation Therapy Plans: An Application to Gastric Cancer

(1) Aim: This study was conducted to radiobiologically compare radiotherapy plans for gastric cancer with a newly developed software tool. (2) Methods: Treatment planning was performed on two computational phantoms simulating adult male and female patients. Three-dimensional conformal radiotherapy (3D-CRT), intensity-modulated radiation therapy (IMRT) and volumetric modulated arc therapy (VMAT) plans for gastric cancer were generated with three-photon beam energies. The equivalent uniform dose (EUD), tumor control probability (TCP) of the target and normal tissue control probability (NTCP) of eight different critical organs were calculated. A new software was employed for these calculations using the EUD-based model and dose-volume-histogram data. (3) Results: The IMRT and VMAT plan led to TCPs of 51.3-51.5%, whereas 3D-CRT gave values up to 50.2%. The intensity-modulated techniques resulted in NTCPs of (5.3 × 10^-6-3.3 × 10^-1)%. The corresponding NTCPs from 3D-CRT were (3.4 × 10^-7-7.4 × 10^-1)%. The above biological indices were automatically calculated in less than 40 s with the software. (4) Conclusions: The direct and quick radiobiological evaluation of radiotherapy plans is feasible using the new software tool. The IMRT and VMAT reduced the probability of the appearance of late effects in most of the surrounding critical organs and slightly increased the TCP compared to 3D-CRT.

Machine learning-based normal tissue complication probability model for predicting albumin-bilirubin (ALBI) grade increase in hepatocellular carcinoma patients

PURPOSE

The aim of this study was to develop a normal tissue complication probability model using a machine learning approach (ML-based NTCP) to predict the risk of radiation-induced liver disease in hepatocellular carcinoma (HCC) patients.

MATERIALS AND METHODS

The study population included 201 HCC patients treated with radiotherapy. The patients' medical records were retrospectively reviewed to obtain the clinical and radiotherapy data. Toxicity was defined by albumin-bilirubin (ALBI) grade increase. The normal liver dose-volume histogram was reduced to mean liver dose (MLD) based on the fraction size-adjusted equivalent uniform dose (2 Gy/fraction and α/β = 2). Three types of ML-based classification models were used, a penalized logistic regression (PLR), random forest (RF), and gradient-boosted tree (GBT) model. Model performance was compared using the area under the receiver operating characteristic curve (AUROC). Internal validation was performed by 5-fold cross validation and external validation was done in 44 new patients.

RESULTS

Liver toxicity occurred in 87 patients (43.1%). The best individual model was the GBT model using baseline liver function, liver volume, and MLD as inputs and the best overall model was an ensemble of the PLR and GBT models. An AUROC of 0.82 with a standard deviation of 0.06 was achieved for the internal validation. An AUROC of 0.78 with a standard deviation of 0.03 was achieved for the external validation. The behaviors of the best GBT model were also in good agreement with the domain knowledge on NTCP.

CONCLUSION

We propose the methodology to develop an ML-based NTCP model to estimate the risk of ALBI grade increase.

Good Timing Matters: The Spatially Fractionated High Dose Rate Boost Should Come First

Monoplanar microbeam irradiation (MBI) and pencilbeam irradiation (PBI) are two new concepts of high dose rate radiotherapy, combined with spatial dose fractionation at the micrometre range. In a small animal model, we have explored the concept of integrating MBI or PBI as a simultaneously integrated boost (SIB), either at the beginning or at the end of a conventional, low-dose rate schedule of 5x4 Gy broad beam (BB) whole brain radiotherapy (WBRT). MBI was administered as array of 50 µm wide, quasi-parallel microbeams. For PBI, the target was covered with an array of 50 µm × 50 µm pencilbeams. In both techniques, the centre-to-centre distance was 400 µm. To assure that the entire brain received a dose of at least 4 Gy in all irradiated animals, the peak doses were calculated based on the daily BB fraction to approximate the valley dose. The results of our study have shown that the sequence of the BB irradiation fractions and the microbeam SIB is important to limit the risk of acute adverse effects, including epileptic seizures and death. The microbeam SIB should be integrated early rather than late in the irradiation schedule.

Comparison of dosimetric effects of MLC positional errors on VMAT and IMRT plans for SBRT radiotherapy in non-small cell lung cancer

The positional accuracy of multi-leaf collimators (MLC) is important in stereotactic body radiotherapy (SBRT). The aim of this study was to investigate the impact between MLC positional error and dosimetry of volume intensity modulated (VMAT) and general intensity modulated (IMRT) plans for non-small cell lung cancer (NSCLC). Fifteen patients with NSCLC were selected to design the 360 SBRT-VMAT plans and the 360 SBRT-IMRT error plans. The DICOM files for these treatment plans were imported into a proprietary computer program that introduced delivery errors. Random and systematic MLC position (0.1, 0.2, 0.5, 1.0, 1.5, and 2.0 mm) errors were introduced. The systematic errors were shift errors (caused by gravity), opening errors, and closing errors. The CI, GI, d2cm and generalized equivalent uniform dose (gEUD) were calculated for the original plan and all treatment plans, accounting for the errors. Dose sensitivity was calculated using linear regression for MLC position errors. The random MLC errors were relatively insignificant. MLC shift, opening, and closing errors had a significant effect on the dose distribution of the SBRT plan. VMAT was more significant than IMRT. To ensure that the gEUD variation of PTV is controlled within 2%, the shift error, opening error, and closing error of IMRT should be less than 2.4 mm, 1.15 mm, and 0.97 mm, respectively. For VMAT, the shift error, opening error, and closing error should be less than 0.95 mm, 0.32 mm, and 0.38 mm, respectively. The dose sensitivity results obtained in this study can be used as a guide for patient-based quality assurance efforts. The position error of the MLC system had a significant impact on the gEUD of the SBRT technology. The MLC systematic error has a greater dosimetric impact on the VMAT plan than on the IMRT plan for SBRT, which should be carefully monitored.

Dosimetric and biological comparisons of single planning and double plannings for bilateral lung cancer SBRT planning based on the Cyber-Knife system

Objective

The aim is to investigate the influence of single planning (Plan S) and double plannings (Plan D) on bilateral lung cancer stereotactic body radiation therapy planning from the perspective of dosimetry and biology respectively. Methods Cases with bilateral lung cancer patients who had undergone SBRT with the Cyber-Knife were enrolled, and a single planning and double plannings were designed in the Multiplan@4.2 treatment planning system equipped with the Cyber-Knife system. The single plan was to optimize the two target volumes in a separate plan, while the dual plan is to optimize two target volumes respectively in two separate plans, then perform dose superposition. Then based on the dosimetric results, the biological parameters were calculated. Thus the quality of SBRT plans for those bilateral lung cancer designed by the two methods were compared and evaluated according to the dosimetric and biological results.

Results

The dose distribution of both planning target volumes and surrounding organs at risk in Plan S and Plan D could meet the clinical prescription requirements. The target conformity index and the new conformity index of PTV were closer to 1 in the Double plannings, and the dose gradient GI in the Plan D was smaller than Plan S. For organs at risks, the doses received by the Plan D were relatively small. In terms of biological models, for the equivalent uniform dose of normal lung tissue, heart and esophagus, the Plan D was 6.51% (P=0.045), 19.8% (P=0.022), 27.08% (P&gt;0.05) lower than Plan S respectively. The results showed that the equivalent uniform dose of normal tissue in the Plan D was lower relative to Plan S.

Conclusions

Dosimetric and biological results show that both the use of Plan D have an advantage of protecting normal tissues, and it was suggested that to design double plannings for bilateral lung cancer stereotactic body radiation therapy planning based on Cyber-Knife in the clinical practice.

The impact of organ motion and the appliance of mitigation strategies on the effectiveness of hypoxia-guided proton therapy for non-small cell lung cancer

BACKGROUND AND PURPOSE

To investigate the impact of organ motion on hypoxia-guided proton therapy treatments for non-small cell lung cancer (NSCLC) patients.

MATERIALS AND METHODS

Hypoxia PET and 4D imaging data of six NSCLC patients were used to simulate hypoxia-guided proton therapy with different motion mitigation strategies including rescanning, breath-hold, respiratory gating and tumour tracking. Motion-induced dose degradation was estimated for treatment plans with dose painting of hypoxic tumour sub-volumes at escalated dose levels. Tumour control probability (TCP) and dosimetry indices were assessed to weigh the clinical benefit of dose escalation and motion mitigation. In addition, the difference in normal tissue complication probability (NTCP) between escalated proton and photon VMAT treatments has been assessed.

RESULTS

Motion-induced dose degradation was found for target coverage (CTV V_95% up to -4%) and quality of the dose-escalation-by-contour (Q_RMS up to 6%) as a function of motion amplitude and amount of dose escalation. The TCP benefit coming from dose escalation (+4-13%) outweighs the motion-induced losses (<2%). Significant average NTCP reductions of dose-escalated proton plans were found for lungs (-14%), oesophagus (-10%) and heart (-16%) compared to conventional VMAT plans. The best plan dosimetry was obtained with breath hold and respiratory gating with rescanning.

CONCLUSION

NSCLC affected by hypoxia appears to be a prime target for proton therapy which, by dose-escalation, allows to mitigate hypoxia-induced radio-resistance despite the sensitivity to organ motion. Furthermore, substantial reduction in normal tissue toxicity can be expected compared to conventional VMAT. Accessibility and standardization of hypoxia imaging and clinical trials are necessary to confirm these findings in a clinical setting.

Comparison of PRIMO Monte Carlo code and Eclipse treatment planning system in calculation of dosimetric parameters in brain cancer radiotherapy

Background

It is important to evaluate the dose calculated by treatment planning systems (TPSs) and dose distribution in tumor and organs at risk (OARs). The aim of this study is to compare dose calculated by the PRIMO Monte Carlo code and Eclipse TPS in radiotherapy of brain cancer patients.

Materials and methods

PRIMO simulation code was used to simulate a Varian Clinac 600C linac. The simulations were validated for the linac by comparison of the simulation and measured results. In the case of brain cancer patients, the dosimetric parameters obtained by the PRIMO code were compared with those calculated by Eclipse TPS. Gamma function analysis with 3%, 3 mm criteria was utilized to compare the dose distributions. The evaluations were based on the dosimetric parameters for the planning target volume (PTV) and OAR including D min, D mean, and D max, homogeneity index (HI), and conformity index (CI).

Results

The gamma function analysis showed a 98% agreement between the results obtained by the PRIMO code and measurement for the percent depth dose (PDD) and dose profiles. The corresponding value in comparing the dosimetric parameters from PRIMO code and Eclipse TPS for the brain patients was 94%, on average. The results of the PRIMO simulation were in good agreement with the measured data and Eclipse TPS calculations.

Conclusions

Based on the results of this study, the PRIMO code can be utilized to simulate a medical linac with good accuracy and to evaluate the accuracy of treatment plans for patients with brain cancer.

Stereotactic Body Radiotherapy Boost with the CyberKnife for Locally Advanced Cervical Cancer: Dosimetric Analysis and Potential Clinical Benefits

(1) Aim: To compare the treatment plans of stereotactic body radiotherapy (SBRT) with CyberKnife (CK) and high-dose-rate (HDR) intracavitary/interstitial brachytherapy (IC/ISBT) and examine the feasibility of CK-SBRT as a viable alternative to BT in patients with locally advanced cervical cancer (LACC). (2) Methods: A BT plan of 28 Gy in four fractions delivered previously to 20 patients with LACC was compared with a CK plan based on the same CT images with structures delineation for BT. The SBRT treatment plan was further divided according to two different approaches, with the high-risk planning target volume (HR-PTV) defined by the high-risk clinical target volume (HR-CTV) without and with a 5 mm margin, which were named CK-CTV plan and CK-PTV plan, respectively. The dose distributions and dosimetric parameters of the target volumes and organs at risk (OARs) were recorded and compared for the three boost plans. Radiobiological metrics were calculated based on the EUD for the hybrid plans. Additionally, the relationship between tumor volume and tolerance doses for the OARs in the BT plan and CK-PTV plan was investigated. (3) Results: Target coverage was better with the CK plan than with the BT plan, as the D95%, D98%, HI and CI of the CK-CTV plan and CK-PTV plan were higher than those of the BT plan; an exception was the D50%. Similarly, the TCP of the target was also significantly in favor of the CK hybrid plans (p < 0.01). For the OARs, the CK-CTV plan was superior to the BT plan as regards the rectum D2cc, bladder D2cc and bladder Dmax. The CK-PTV plan could achieve dosimetric parameters comparable to those of the BT plan for OARs concerning the small residual tumor volume. The NTCP of the rectum for the WPI+CK-CTV plans was significantly lower than that of the WPI+BT plans (p < 0.01). (4) Conclusions: CK-based SBRT can achieve better target coverage, dose sparing for the OARs and radiobiological effects compared with the BT plan for tumors that are not excessively large. CK-based SBRT could be an alternative option to administer a radiation boost for patients with LACC.

Longitudinally Heterogeneous Tumor Dose Optimizes Proton Broadbeam, Interlaced Minibeam, and FLASH Therapy

Simple Summary

The aim of any kind of external radiation therapy is to control a tumor with the highest possible probability of the lowest possible side effects. Here, we study further opportunities of reducing the side effects of proton therapy by applying longitudinally heterogeneous dose distributions in the tumor respecting the delivery of a minimum prescribed dose. In our simulations, the longitudinally heterogeneous dose distributions show a reduced dose in the healthy tissue already in the case of proton broadbeam irradiations, but a much higher (calculated) mean cell survival in the case of proton minibeam irradiation. This demonstrates its potential to substantially reduce side effects at a simultaneously higher tumor control probability, opening new opportunities of easier application when striving for high dose-rate applications of proton beams (>~10 Gy/s), in order to additionally profit from the so-called FLASH effects.

Abstract

The prerequisite of any radiation therapy modality (X-ray, electron, proton, and heavy ion) is meant to meet at least a minimum prescribed dose at any location in the tumor for the best tumor control. In addition, there is also an upper dose limit within the tumor according to the International Commission on Radiation Units (ICRU) recommendations in order to spare healthy tissue as well as possible. However, healthy tissue may profit from the lower side effects when waving this upper dose limit and allowing a larger heterogeneous dose deposition in the tumor, but maintaining the prescribed minimum dose level, particularly in proton minibeam therapy. Methods: Three different longitudinally heterogeneous proton irradiation modes and a standard spread-out Bragg peak (SOBP) irradiation mode are simulated for their depth-dose curves under the constraint of maintaining a minimum prescribed dose anywhere in the tumor region. Symmetric dose distributions of two opposing directions are overlaid in a 25 cm-thick water phantom containing a 5 cm-thick tumor region. Interlaced planar minibeam dose distributions are compared to those of a broadbeam using the same longitudinal dose profiles. Results and Conclusion: All longitudinally heterogeneous proton irradiation modes show a dose reduction in the healthy tissue compared to the common SOBP mode in the case of broad proton beams. The proton minibeam cases show eventually a much larger mean cell survival and thus a further reduced equivalent uniform dose (EUD) in the healthy tissue than any broadbeam case. In fact, the irradiation mode using only one proton energy from each side shows better sparing capabilities in the healthy tissue than the common spread-out Bragg peak irradiation mode with the option of a better dose fall-off at the tumor edges and an easier technical realization, particularly in view of proton minibeam irradiation at ultra-high dose rates larger than ~10 Gy/s (so-called FLASH irradiation modes).

Dosimetric analysis and biological evaluation between proton radiotherapy and photon radiotherapy for the long target of total esophageal squamous cell carcinoma

Objective

The purpose of this study is to compare the dosimetric and biological evaluation differences between photon and proton radiation therapy.

Methods

Thirty esophageal squamous cell carcinoma (ESCC) patients were generated for volumetric modulated arc therapy (VMAT) planning and intensity-modulated proton therapy (IMPT) planning to compare with intensity-modulated radiation therapy (IMRT) planning. According to dose–volume histogram (DVH), dose–volume parameters of the plan target volume (PTV) and homogeneity index (HI), conformity index (CI), and gradient index (GI) were used to analyze the differences between the various plans. For the organs at risk (OARS), dosimetric parameters were compared. Tumor control probability (TCP) and normal tissue complication probability (NTCP) was also used to evaluate the biological effectiveness of different plannings.

Results

CI, HI, and GI of IMPT planning were significantly superior in the three types of planning (p < 0.001, p < 0.001, and p < 0.001, respectively). Compared to IMRT and VMAT planning, IMPT planning improved the TCP (p<0.001, p<0.001, respectively). As for OARs, IMPT reduced the bilateral lung and heart accepted irradiation dose and volume. The dosimetric parameters, such as mean lung dose (MLD), mean heart dose (MHD), V₅, V₁₀, and V₂₀, were significantly lower than IMRT or VMAT. IMPT afforded a lower maximum dose (D_max) of the spinal cord than the other two-photon plans. What’s more, the radiation pneumonia of the left lung, which was caused by IMPT, was lower than IMRT and VMAT. IMPT achieved the pericarditis probability of heart is only 1.73% ± 0.24%. For spinal cord myelitis necrosis, there was no significant difference between the three different technologies.

Conclusion

Proton radiotherapy is an effective technology to relieve esophageal cancer, which could improve the TCP and spare the heart, lungs, and spinal cord. Our study provides a prediction of radiotherapy outcomes and further guides the individual treatment.

Dosimetry in Radiopharmaceutical Therapy

The application of radiopharmaceutical therapy for the treatment of certain diseases is well established, and the field is expanding. New therapeutic radiopharmaceuticals have been developed in recent years, and more are in the research pipeline. Concurrently, there is growing interest in the use of internal dosimetry as a means of personalizing, and potentially optimizing, such therapy for patients. Internal dosimetry is multifaceted, and the current state of the art is discussed in this continuing education article. Topics include the context of dosimetry, internal dosimetry methods, the advantages and disadvantages of incorporating dosimetry calculations in radiopharmaceutical therapy, a description of the workflow for implementing patient-specific dosimetry, and future prospects in the field.

A novel knowledge-based prediction model for estimating an initial equivalent uniform dose in semi-auto-planning for cervical cancer

Background

We developed a novel concept, equivalent uniform length (EUL), to describe the relationship between the generalized equivalent uniform dose (EUD) and the geometric anatomy around a tumor target. By correlating EUL with EUD, we established two EUD–EUL knowledge-based (EEKB) prediction models for the bladder and rectum that predict initial EUD values for generating quality treatment plans.

Methods

EUL metrics for the rectum and bladder were extracted and collected from the intensity-modulated radiotherapy therapy (IMRT) plans of 60 patients with cervical cancer. The two EEKB prediction models were built using linear regression to establish the relationships between EUL_r and EUD_r (EUL and EUD of rectum) and EUL_b, and EUD_b (EUL and EUD of bladder), respectively. The EE plans were optimized by incorporating the predicted initial EUD parameters for the rectum and bladder with the conventional pinnacle auto-planning (PAP) initial dose parameters for other organs. The efficiency of the predicted initial EUD values were then evaluated by comparing the consistency and quality of the EE plans, PAP plans (based on default PAP initial parameters), and manual plans (designed manually by different dosimetrists) for a sample of 20 patients.

Results

Linear regression analyses showed a significant correlation between EUL and EUD (R² = 0.79 and 0.69 for EUD_b and EUD_r, respectively). In a sample of 20 patients, the average bladder V40 and V50 derived from the EE plans were significantly lower (V40: 30.00 ± 5.76, V50: 14.36 ± 4.00) than the V40 and V50 values derived from manual plans (V40: 36.03 ± 8.02, V50: 19.02 ± 5.42). Compared with the PAP plans, the EE plans produced significantly lower average V30 and Dmean values for the bladder (V30: 50.55 ± 6.33, Dmean: 31.48 ± 1.97 Gy).

Conclusions

Our EEKB prediction models predicted reasonable initial EUD values for the rectum and bladder based on patient-specific geometric EUL values, thereby improving optimization and planning efficiency.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13014-022-02120-4.

Trade-off in healthy tissue sparing of FLASH and fractionation in stereotactic proton therapy of lung lesions with transmission beams

PURPOSE AND OBJECTIVE

Besides a dose-rate threshold of 40-100 Gy/s, the FLASH effect may require a dose >3.5-7 Gy. Even in hypofractioned treatments, with all beams delivered in each fraction (ABEF), most healthy tissue is irradiated to a lower fraction dose. This can be circumvented by single-beam-per-fraction (SBPF) delivery, with a loss of healthy tissue sparing by fractionation. We investigated the trade-off between FLASH and loss of fractionation in SBPF stereotactic proton therapy of lung cancer and determined break-even FLASH-enhancement ratios (FERs).

MATERIALS AND METHODS

Treatment plans for 12 patients were generated. GTV delineations were available and a 5 mm GTV-PTV margin was applied. Equiangular arrangements of 3, 5, 7, and 9 244 MeV proton transmission beams were used. To facilitate SBPF, the number of fractions was equal to the number of beams. Iso-effective fractionation schedules with a single field uniform dose prescription were used: D_95%,PTV = 100%D_pres per beam. All plans were evaluated in terms of dose to lung and conformity of dose to target of FLASH-enhanced biologically equivalent dose (EQD2).

RESULTS

Compared to ABEF, SBPF resulted in a median increase of EQD2_mean to healthy lung of 56%, 58%, 55% and 54% in plans with 3, 5, 7 and 9 fractions respectively and of 90%, 108%, 106% and 102% in V_{100% EQD2}, quantifying conformity. This can be compensated for by FERs of at least 1.28, 1.32, 1.30 and 1.23 respectively for EQD2_mean and 1.29, 1.18, 1.28 and 1.15 for V_100%,EQD2.

CONCLUSION

A FLASH effect outweighing the loss of fractionation in SBPF may be achieved in stereotactic lung treatments. The trade-off with fractionation depends on the conditions under which the FLASH effect occurs. Better understanding of the underlying biology and the impact of delivery conditions is needed.

Microbeam Irradiation as a Simultaneously Integrated Boost in a Conventional Whole-Brain Radiotherapy Protocol

Microbeam radiotherapy (MRT), an experimental high-dose rate concept with spatial fractionation at the micrometre range, has shown a high therapeutic potential as well as good preservation of normal tissue function in pre-clinical studies. We investigated the suitability of MRT as a simultaneously integrated boost (SIB) in conventional whole-brain irradiation (WBRT). A 174 Gy MRT SIB was administered with an array of quasi-parallel, 50 µm wide microbeams spaced at a centre-to-centre distance of 400 µm either on the first or last day of a 5 × 4 Gy radiotherapy schedule in healthy adult C57 BL/6J mice and in F98 glioma cell cultures. The animals were observed for signs of intracranial pressure and focal neurologic signs. Colony counts were conducted in F98 glioma cell cultures. No signs of acute adverse effects were observed in any of the irradiated animals within 3 days after the last irradiation fraction. The tumoricidal effect on F98 cell in vitro was higher when the MRT boost was delivered on the first day of the irradiation course, as opposed to the last day. Therefore, the MRT SIB should be integrated into a clinical radiotherapy schedule as early as possible.

Milestones in dosimetry for nuclear medicine therapy

Nuclear Medicine therapy has reached a critical juncture with an unprecedented number of patients being treated and an extensive list of new radiopharmaceuticals under development. Since the early applications of these treatments dosimetry has played a vital role in their development, in both aiding optimisation and enhancing safety and efficacy. To inform the future direction of this field, it is useful to reflect on the scientific and technological advances that have occurred since those early uses. In this review, we explore how dosimetry has evolved over the years and discuss why such initiatives were conceived and the importance of maintaining standards within our practise. Specific milestones and landmark publications are highlighted and a thematic review and significant outcomes during each decade are presented.

Radiobiological Comparison of Acuros External Beam and Anisotropic Analytical Algorithm on Esophageal Carcinoma Radiotherapy Treatment Plans

Objective

The present study aimed to investigate the dose differences and radiobiological assessment between Anisotropic Analytical Algorithm (AAA) and Acuros External Beam (AXB) with its 2 calculation models, namely, dose-to-water (AXB-Dw) and dose-to-medium (AXB-Dm), on esophageal carcinoma radiotherapy treatment plans.

Materials and methods

The AXB-Dw and AXB-Dm plans were generated by recalculating the initial 66 AAA plans using the AXB algorithm with the same monitor units and beam parameters as those in the original plan. The dosimetric and radiobiological assessment parameters were calculated for the planning target volume (PTV) and organs at risk (OARs). The gamma agreement for the PTV and the correlation between it and the volume of the air cavity and bone among the different algorithms were compared simultaneously. The dose discrepancy between the theoretical calculation and treatment planning system (TPS) when switching from AXB-Dm to AXB-Dw was analyzed according to the composition of the structures.

Results

The PTV dose of AXB-Dm plans was significantly smaller than that of the AAA and AXB-Dw plans (P < .05), except for D₂. The difference values for AAA vs AXB-Dm (∆D_{x,(AAA-AXB,Dm)}) and AXB-Dw vs AXB-Dm (∆D_{x,(AXB,Dw-AXB,Dm)}) were 1.94% [1.27%, 2.64%] and 1.95% [1.56%, 2.27%], respectively. For the spinal cord and heart, there were obvious differences between the AAA vs AXB-Dm (spinal cord: 1.15%, heart: 2.89%) and AXB-Dw vs AXB-Dm (spinal cord: 1.88%, heart: 3.25%) plans. For the lung, the differences between AAA vs AXB-Dm and AAA vs AXB-Dw were significantly larger than those of AXB-Dm vs AXB-Dw. Compared to the case of AAA and AXB-Dw, the decrease in biologically effective dose (BED¹⁰, αβ=10 ) of AXB-Dm due to dose non-uniformity exceeded 6.5%, even for a small σ. The average values of equivalent uniform dose in the AAA, AXB-Dw, and AXB-Dm plans were 52.03±.39 Gy, 52.24 ± .81 Gy, and 51.13 ± .47 Gy, respectively. The tumor control probability (TCP) results for PTV in the AAA, AXB-Dw, and AXB-Dm plans were 62.29 ± 1.57%, 62.82 ± 1.69%, and 58.68±1.88%, respectively. With the 2%/2 mm and 3%/3 mm acceptance criteria, the mean values of  ΔγAAAAXB−Dw, ΔγAAAAXB−Dm, and  ΔγAXB−DmAXB−Dw were 87.24, 63.3, and 64.81% vs 97.86, 91.77, and 89.25%, respectively. The dose discrepancy between the theoretical calculation and TPS when switching from AXB-Dm to AXB-Dw was approximately 1.63%.

Conclusions

The AAA and AXB-Dw algorithms overestimated the radiobiological parameters when the tumor particularly consisted of nonuniform tissues. A relatively small dose difference could cause a significant reduction in the corresponding TCP. Dose distribution algorithms should be carefully chosen by physicists and oncologists to improve tumor control, as well as to optimize OARs protection.

Evaluation of the Radiobiological Models predicting the Radiation-Induced Hypothyroidism in the Partially Irradiated Thyroid Gland of Patients with Breast Cancer

Background: Radiation-induced hypothyroidism (RHT) is one of the side effects that might have an impact on the quality of life of patients with breast cancer treated with radiotherapy. Objectives: The aim of the current study was to evaluate the performances of the Lyman-Kutcher-Burman (LKB) and Log-Logistic models in the prediction of hypothyroidism (HT) as well as the estimation of the model parameters for the incidence of RHT among patients with breast cancer. Methods: Fifty-two patients treated with radiation therapy (RT) for breast cancer were prospectively evaluated. Patients' serum samples [tri-iodothyronine, thyroxine, thyroid-stimulating hormone (TSH), free triiodothyronine, and free thyroxine] were measured before RT and also at a regular time interval until 1 year after the completion of RT. For each patient, dose-volume histograms (DVHs) of the thyroid gland were derived from their treatment planning dataset. Patients whose TSH levels were higher than normal with a decrease in FT4 levels were considered as cases with RHT. The LKB and Log-Logistic radiobiological models were evaluated by comparing them with the resultant follow-up data. The parameters for radiobiological models have been deduced by fitting the models to the follow-up data. The models were fitted in a Bayesian setting and compared according to the widely applicable information criterion (WAIC). Results: Twenty-one (40%) patients developed RHT at a follow-up of 1 year after the end of radiation treatment. The fitted values of D50 for the LKB and Log-Logistic models were 37.71 and 25.50 Gy, respectively for the partially irradiated thyroid of patients with breast cancer. The mean time to the incidence of RHT was obtained at 6.7 months in the studied group. Conclusions: A volumetric effect was found for the thyroid gland in the implemented normal tissue complication probability models. Compared to the follow-up data, the Log-Logistic model was ranked as the best model for predicting the rate of RHT in patients with breast cancer.

Can dose convolution modelling explain bath and shower effects in rat spinal cord?

Objective. ‘Bath and shower’ effects were first seen in proton irradiations of rat spinal cord, where a low dose ‘bath’ reduced the smaller field ‘shower’ dose needed for limb paralysis giving the appearance of sensitisation of the cord or disproportionate response. This was difficult to reconcile with existing tissue complication models. The purpose of this investigation is to explore a different approach using a dose convolution algorithm to model the 50% isoeffect endpoint. Approach. Bath and shower dose distributions were convolved with Gaussian functions with widths specified by the σ parameter. The hypothesis was that the maximum value from the convolved distributions was constant for isoeffect across the modelled scenarios. A simpler field length dependent relative biological effectiveness (FLRBE) approach was also used for a subset of the data which gave results independent of σ. Main results. The maximum values from the convolved distributions were constant within ±17% across the bath and shower experiments for σ = 3.5 mm, whereas the maximum dose varied by a factor of four. The FLRBE results were also within ±14% confirming the validity of the dose convolution approach. Significance. A simple approach using dose convolution modelling of the 50% isotoxicity gave compelling consistency with the full range of bath and shower results, while the FLRBE approach confirmed the results for the symmetric field data. Convolution modelling and the effect of time interval were consistent with a signalling factor diffusion mechanism such as the ‘bystander effect’. The results suggest biological effectiveness is reduced for very small field sizes, requiring a higher isoeffect dose. By implication, the bath dose does not sensitise the cord to the shower dose; when biological effectiveness is accounted for, a small increase in the bath dose requires a significantly larger reduction in shower dose for isoeffect.

Radial Data Mining to Identify Density-Dose Interactions That Predict Distant Failure Following SABR

Purpose

Lower dose outside the planned treatment area in lung stereotactic radiotherapy has been linked to increased risk of distant metastasis (DM) possibly due to underdosage of microscopic disease (MDE). Independently, tumour density on pretreatment computed tomography (CT) has been linked to risk of MDE. No studies have investigated the interaction between imaging biomarkers and incidental dose. The interaction would showcase whether the impact of dose on outcome is dependent on imaging and, hence, if imaging could inform which patients require dose escalation outside the gross tumour volume (GTV). We propose an image-based data mining methodology to investigate density-dose interactions radially from the GTV to predict DM with no a priori assumption on location.

Methods

Dose and density were quantified in 1-mm annuli around the GTV for 199 patients with early-stage lung cancer treated with 60 Gy in 5 fractions. Each annulus was summarised by three density and three dose parameters. For parameter combinations, Cox regressions were performed including a dose-density interaction in independent annuli. Heatmaps were created that described improvement in DM prediction due to the interaction. Regions of significant improvement were identified and studied in overall outcome models.

Results

Dose-density interactions were identified that significantly improved prediction for over 50% of bootstrap resamples. Dose and density parameters were not significant when the interaction was omitted. Tumour density variance and high peritumour density were associated with DM for patients with more cold spots (less than 30-Gy EQD2) and non-uniform dose about 3 cm outside of the GTV. Associations identified were independent of the mean GTV dose.

Conclusions

Patients with high tumour variance and peritumour density have increased risk of DM if there is a low and non-uniform dose outside the GTV. The dose regions are independent of tumour dose, suggesting that incidental dose may play an important role in controlling occult disease. Understanding such interactions is key to identifying patients who will benefit from dose-escalation. The methodology presented allowed spatial dose-density interactions to be studied at the exploratory stage for the first time. This could accelerate the clinical implementation of imaging biomarkers by demonstrating the impact of incidental dose for tumours of varying characteristics in routine data.

Radiobiological and Treatment-Related Aspects of Spatially Fractionated Radiotherapy

The continuously evolving field of radiotherapy aims to devise and implement techniques that allow for greater tumour control and better sparing of critical organs. Investigations into the complexity of tumour radiobiology confirmed the high heterogeneity of tumours as being responsible for the often poor treatment outcome. Hypoxic subvolumes, a subpopulation of cancer stem cells, as well as the inherent or acquired radioresistance define tumour aggressiveness and metastatic potential, which remain a therapeutic challenge. Non-conventional irradiation techniques, such as spatially fractionated radiotherapy, have been developed to tackle some of these challenges and to offer a high therapeutic index when treating radioresistant tumours. The goal of this article was to highlight the current knowledge on the molecular and radiobiological mechanisms behind spatially fractionated radiotherapy and to present the up-to-date preclinical and clinical evidence towards the therapeutic potential of this technique involving both photon and proton beams.