A multi-plane source model for out-of-field head scatter dose calculations in external beam photon therapy

Deep-Learning for Rapid Estimation of the Out-of-Field Dose in External Beam Photon Radiation Therapy - A Proof of Concept

PURPOSE

The dose deposited outside of the treatment field during external photon beam radiation therapy treatment, also known as out-of-field dose, is the subject of extensive study as it may be associated with a higher risk of developing a second cancer and could have deleterious effects on the immune system that compromise the efficiency of combined radio-immunotherapy treatments. Out-of-field dose estimation tools developed today in research, including Monte Carlo simulations and analytical methods, are not suited to the requirements of clinical implementation because of their lack of versatility and their cumbersome application. We propose a proof of concept based on deep learning for out-of-field dose map estimation that addresses these limitations.

METHODS AND MATERIALS

For this purpose, a 3D U-Net, considering as inputs the in-field dose, as computed by the treatment planning system, and the patient's anatomy, was trained to predict out-of-field dose maps. The cohort used for learning and performance evaluation included 3151 pediatric patients from the FCCSS database, treated in 5 clinical centers, whose whole-body dose maps were previously estimated with an empirical analytical method. The test set, composed of 433 patients, was split into 5 subdata sets, each containing patients treated with devices unseen during the training phase. Root mean square deviation evaluated only on nonzero voxels located in the out-of-field areas was computed as performance metric.

RESULTS

Root mean square deviations of 0.28 and 0.41 cGy/Gy were obtained for the training and validation data sets, respectively. Values of 0.27, 0.26, 0.28, 0.30, and 0.45 cGy/Gy were achieved for the 6 MV linear accelerator, 16 MV linear accelerator, Alcyon cobalt irradiator, Mobiletron cobalt irradiator, and betatron device test sets, respectively.

CONCLUSIONS

This proof-of-concept approach using a convolutional neural network has demonstrated unprecedented generalizability for this task, although it remains limited, and brings us closer to an implementation compatible with clinical routine.

Analytical models for external photon beam radiotherapy out-of-field dose calculation: a scoping review

A growing body of scientific evidence indicates that exposure to low dose ionizing radiation (< 2 Gy) is associated with a higher risk of developing radio-induced cancer. Additionally, it has been shown to have significant impacts on both innate and adaptive immune responses. As a result, the evaluation of the low doses inevitably delivered outside the treatment fields (out-of-field dose) in photon radiotherapy is a topic that is regaining interest at a pivotal moment in radiotherapy. In this work, we proposed a scoping review in order to identify evidence of strengths and limitations of available analytical models for out-of-field dose calculation in external photon beam radiotherapy for the purpose of implementation in clinical routine. Papers published between 1988 and 2022 proposing a novel analytical model that estimated at least one component of the out-of-field dose for photon external radiotherapy were included. Models focusing on electrons, protons and Monte-Carlo methods were excluded. The methodological quality and potential limitations of each model were analyzed to assess their generalizability. Twenty-one published papers were selected for analysis, of which 14 proposed multi-compartment models, demonstrating that research efforts are directed towards an increasingly detailed description of the underlying physical phenomena. Our synthesis revealed great inhomogeneities in practices, in particular in the acquisition of experimental data and the standardization of measurements, in the choice of metrics used for the evaluation of model performance and even in the definition of regions considered out-of-the-field, which makes quantitative comparisons impossible. We therefore propose to clarify some key concepts. The analytical methods do not seem to be easily suitable for massive use in clinical routine, due to the inevitable cumbersome nature of their implementation. Currently, there is no consensus on a mathematical formalism that comprehensively describes the out-of-field dose in external photon radiotherapy, partly due to the complex interactions between a large number of influencing factors. Out-of-field dose calculation models based on neural networks could be promising tools to overcome these limitations and thus favor a transfer to the clinic, but the lack of sufficiently large and heterogeneous data sets is the main obstacle.

A publicly available dataset of out-of-field dose profiles of a 6 MV linear accelerator

An increase in radiotherapy-induced secondary malignancies has led to recent developments in analytical modelling of out-of-field dose. These models must be validated against measurements, but currently available datasets are outdated or limited in scope. This study aimed to address these shortcomings by producing a large dataset of out-of-field dose profiles measured with modern equipment. A novel method was developed with the intention of allowing physicists in all clinics to perform these measurements themselves using commonly available dosimetry equipment. A standard 3D scanning water tank was used to collect 36 extended profiles. Each profile was measured in two sections, with the inner section measured with the beam directly incident on the tank, and the outer section with the beam incident on a water-equivalent phantom abutted next to the tank. The two sections were then stitched using a novel feature-matching approach. The profiles were compared against linac commissioning data and manually inspected for discontinuities in the overlap region. The dataset is presented as a publicly accessible comma separated variable file containing off-axis ratios at a range of off-axis distances. This dataset may be applied to the development and validation of analytical models of out-of-field dose. Additionally, it may be used to inform dose estimates to radiosensitive implants and anatomy. Physicists are encouraged to perform these out-of-field measurements in their own clinics and share their results with the community.

Systematic Screening of COVID-19 Disease Based on Chest CT and RT-PCR for Cancer Patients Undergoing Radiation Therapy in a Coronavirus French Hotspot

PURPOSE

Patients with cancer are presumed to be more vulnerable to COVID-19. We evaluated a screening strategy combining chest computed tomography (CT) and reverse-transcription polymerase chain reaction (RT-PCR) for patients treated with radiation therapy at our cancer center located in a COVID-19 French hotspot during the first wave of the pandemic.

METHODS AND MATERIALS

Chest CT images were proposed during radiation therapy CT simulation. Images were reviewed by an expert radiologist according to the COVID-19 Reporting and Data System classification. Nasal swabs with RT-PCR assay were initially proposed in cases of suspicious imaging or clinical context and were eventually integrated into the systematic screening. A dedicated radiation therapy workflow was proposed for COVID-19 patients to limit the risk of contamination.

RESULTS

From March 18, 2020 to May 1, 2020, 480 patients were screened by chest CT, and 313 patients had both chest CT and RT-PCR (65%). The cumulative incidence of COVID-19 was 5.4% (95% confidence interval [CI], 3.6-7.8; 26 of 480 patients). Diagnosis of COVID-19 was made before radiation therapy for 22 patients (84.6%) and during RT for 4 patients (15.3%). Chest CT directly aided the diagnosis of 7 cases in which the initial RT-PCR was negative or not feasible, out of a total of 480 patients (1.5%) and 517 chest CT acquisitions. Four patients with COVID-19 at the time of the chest CT screening had a false negative CT. Sensitivity and specificity of chest CT screening in patients with both RT-PCR and chest CT testing were estimated at 0.82 (95% CI, 0.60-0.95) and 0.98 (95% CI, 0.96-0.99), respectively. Adaptation of the radiation therapy treatment was made for all patients, with 7 postponed treatments (median: 5 days; interquartile range, 1.5-14.8).

CONCLUSIONS

The benefit of systematic use of chest CT screening during CT simulation for patients undergoing radiation therapy during the COVID-19 pandemic seemed limited.

Development of clinical application program for radiotherapy induced cancer risk calculation using Monte Carlo engine in volumetric-modulated arc therapy

Background

The purpose of this study is to develop a clinical application program that automatically calculates the effect for secondary cancer risk (SCR) of individual patient. The program was designed based on accurate dose calculations using patient computed tomography (CT) data and Monte Carlo engine. Automated patient-specific evaluation program was configured to calculate SCR.

Methods

The application program is designed to re-calculate the beam sequence of treatment plan using the Monte Carlo engine and patient CT data, so it is possible to accurately calculate and evaluate scatter and leakage radiation, difficult to calculate in TPS. The Monte Carlo dose calculation system was performed through stoichiometric calibration using patient CT data. The automatic SCR evaluation program in application program created with a MATLAB was set to analyze the results to calculate SCR. The SCR for organ of patient was calculated based on Biological Effects of Ionizing Radiation (BEIR) VII models. The program is designed to sequentially calculate organ equivalent dose (OED), excess absolute risk (EAR), excess relative risk (ERR), and the lifetime attributable risk (LAR) in consideration of 3D dose distribution analysis. In order to confirm the usefulness of the developed clinical application program, the result values from clinical application program were compared with the manual calculation method used in the previous study.

Results

The OED values calculated in program were calculated to be at most approximately 13.3% higher than results in TPS. The SCR result calculated by the developed clinical application program showed a maximum difference of 1.24% compared to the result of the conventional manual calculation method. And it was confirmed that EAR, ERR and LAR values can be easily calculated by changing the biological parameters.

Conclusions

We have developed a patient-specific SCR evaluation program that can be used conveniently in the clinic. The program consists of a Monte Carlo dose calculation system for accurate calculation of scatter and leakage radiation and a patient-specific automatic SCR evaluation program using 3D dose distribution. The clinical application program that improved the disadvantages of the existing process can be used as an index for evaluating a patient treatment plan.

Out-of-field scattering from the Integral Quality Monitor® in megavolt photon beams

PURPOSE

To investigate out-of-field scattered doses from the Integral Quality Monitor (IQM) transmission detector in megavoltage photon beams.

MATERIALS AND METHODS

We measured out-of-field point doses for 20 × 20 cm² 6-15 MV photon beams using 10 × 2 cm² Gafchromic EBT3 film strips placed across the surfaces of 5-cm thick water-equivalent RW3 slabs. The films were placed at 10 cm intervals from the central axis (CAX) of each beam, up to 1.0 m away on opposite sides of the CAX. The measurements were conducted at 80 cm and 100 cm source-to-surface distances (SSD) without the IQM, and were repeated with the IQM in the paths of the beams. Measurements were also performed at 90 cm SSD for 20 × 20 and 30 × 30 cm² 15 MV fields. Surface dose profiles were then constructed from the measurements for each beam setup with and without the IQM to examine the differences in scattered dose off-axis. The dose profile for each beam setup was normalised to dose maximum measured on the CAX.

RESULTS

Overall, surface dose profiles acquired with the IQM in the paths of the beams were higher than the corresponding profiles without the IQM. The out-of-field dose increased with increase in photon energy, field size, and shorter SSDs, and decreased with off-axis distance. At 80 cm SSD for the 20 × 20 cm² field, the IQM-induced surface dose ranged from -0.6% ÷ 1.9%, -0.3% ÷ 3.0%, and 0.3% ÷ 6.8% for 6, 10, and 15 MV beams, respectively.

CONCLUSION

The higher surface dose profiles measured with the IQM attached to the linac compared to the profiles without the IQM indicates that the device is acting as an additional source of scattered radiation.

An objective method to evaluate radiation dose distributions varying by three orders of magnitude

PURPOSE

Modern radiotherapy practices typically report the absorbed dose (D) within the 5% relative isodose volume (i.e., the therapeutic dose region) to an accuracy of 3%-5%. Gamma-index analysis, the most commonly used method to evaluate dosimetric accuracy, has low sensitivity to discrepancies that occur outside of this region. The objective of this study was to develop an evaluation method with high sensitivity across dose distributions spanning three orders of magnitude.

METHODS

We generalized the gamma index to include an additional criterion, the absolute absorbed dose difference, specifically for the low-dose region (i.e., D ≤ 5%). We also proposed a method to objectively select the appropriate magnitudes for relative-dose-difference, absolute-dose-difference, and distance-to-agreement criteria. We demonstrated the generalized gamma-index method by first finding the appropriate generalized gamma-index agreement criteria at an interval of specified passing rates. Next, we used the generalized gamma index to evaluate one-, two-, and three-dimensional absorbed dose distributions in a water-box phantom and voxelized patient geometry.

RESULTS

Generalized gamma-index passing rates for one-, two-, and three-dimensional dose distributions were 55.4%, 44.5%, and 8.9%, respectively. Traditional gamma-index passing rates were 100%, 97.8%, and 96.4%, respectively. These results reveal that the generalized method has adequate sensitivity in all regions (i.e., therapeutic and low dose). Additionally, the algorithmic determination of triplets of agreement criteria revealed that they are strong functions of the specified passing rate.

CONCLUSIONS

The major finding of this work is that the proposed method provides an objective evaluation of the agreement of dose distributions spanning three orders of magnitude. In particular, this generalized method correctly characterized dosimetric agreement in the low-dose region, which was not possible by traditional methods. The proposed algorithmic selection of agreement criteria decreased subjectivity and requirements of user judgment and skill. This method could find utility in a variety of applications including dose-algorithm development and translation.

Supplemental computational phantoms to estimate out-of-field absorbed dose in photon radiotherapy

The purpose of this study was to develop a straightforward method of supplementing patient anatomy and estimating out-of-field absorbed dose for a cohort of pediatric radiotherapy patients with limited recorded anatomy. A cohort of nine children, aged 2-14 years, who received 3D conformal radiotherapy for low-grade localized brain tumors (LBTs), were randomly selected for this study. The extent of these patients' computed tomography simulation image sets were cranial only. To approximate their missing anatomy, we supplemented the LBT patients' image sets with computed tomography images of patients in a previous study with larger extents of matched sex, height, and mass and for whom contours of organs at risk for radiogenic cancer had already been delineated. Rigid fusion was performed between the LBT patients' data and that of the supplemental computational phantoms using commercial software and in-house codes. In-field dose was calculated with a clinically commissioned treatment planning system, and out-of-field dose was estimated with a previously developed analytical model that was re-fit with parameters based on new measurements for intracranial radiotherapy. Mean doses greater than 1 Gy were found in the red bone marrow, remainder, thyroid, and skin of the patients in this study. Mean organ doses between 150 mGy and 1 Gy were observed in the breast tissue of the girls and lungs of all patients. Distant organs, i.e. prostate, bladder, uterus, and colon, received mean organ doses less than 150 mGy. The mean organ doses of the younger, smaller LBT patients (0-4 years old) were a factor of 2.4 greater than those of the older, larger patients (8-12 years old). Our findings demonstrated the feasibility of a straightforward method of applying supplemental computational phantoms and dose-calculation models to estimate absorbed dose for a set of children of various ages who received radiotherapy and for whom anatomies were largely missing in their original computed tomography simulations.

A descriptive and broadly applicable model of therapeutic and stray absorbed dose from 6 to 25 MV photon beams

PURPOSE

To develop a simple model of therapeutic and stray absorbed dose for a variety of treatment machines and techniques without relying on proprietary machine-specific parameters.

METHODS

Dosimetry measurements conducted in this study and from the literature were used to develop an analytical model of absorbed dose from a variety of treatment machines and techniques in the 6 to 25 MV interval. A modified one-dimensional gamma-index analysis was performed to evaluate dosimetric accuracy of the model on an independent dataset consisting of measured dose profiles from seven treatment units spanning four manufacturers.

RESULTS

The average difference between the calculated and measured absorbed dose values was 9.9% for those datasets on which the model was trained. Additionally, these results indicate that the model can provide accurate calculations of both therapeutic and stray radiation dose from a wide variety of radiotherapy units and techniques.

CONCLUSIONS

We have developed a simple analytical model of absorbed dose from external beam radiotherapy treatments in the 6 to 25 MV beam energy range. The model has been tested on measured data from multiple treatment machines and techniques, and is broadly applicable to contemporary external beam radiation therapy.

A review of uncertainties in radiotherapy dose reconstruction and their impacts on dose-response relationships

Proper understanding of the risk of radiation-induced late effects for patients receiving external photon beam radiotherapy requires the determination of reliable dose-response relationships. Although significant efforts have been devoted to improving dose estimates for the study of late effects, the most often questioned explanatory variable is still the dose. In this work, based on a literature review, we provide an in-depth description of the radiotherapy dose reconstruction process for the study of late effects. In particular, we focus on the identification of the main sources of dose uncertainty involved in this process and summarise their impacts on the dose-response relationship for radiotherapy late effects. We provide a number of recommendations for making progress in estimating the uncertainties in current studies of radiotherapy late effects and reducing these uncertainties in future studies.

A simple and fast physics-based analytical method to calculate therapeutic and stray doses from external beam, megavoltage x-ray therapy

State-of-the-art radiotherapy treatment planning systems provide reliable estimates of the therapeutic radiation but are known to underestimate or neglect the stray radiation exposures. Most commonly, stray radiation exposures are reconstructed using empirical formulas or lookup tables. The purpose of this study was to develop the basic physics of a model capable of calculating the total absorbed dose both inside and outside of the therapeutic radiation beam for external beam photon therapy. The model was developed using measurements of total absorbed dose in a water-box phantom from a 6 MV medical linear accelerator to calculate dose profiles in both the in-plane and cross-plane direction for a variety of square field sizes and depths in water. The water-box phantom facilitated development of the basic physical aspects of the model. RMS discrepancies between measured and calculated total absorbed dose values in water were less than 9.3% for all fields studied. Computation times for 10 million dose points within a homogeneous phantom were approximately 4 min. These results suggest that the basic physics of the model are sufficiently simple, fast, and accurate to serve as a foundation for a variety of clinical and research applications, some of which may require that the model be extended or simplified based on the needs of the user. A potentially important advantage of a physics-based approach is that the model is more readily adaptable to a wide variety of treatment units and treatment techniques than with empirical models.

Childhood cancer survivor cohorts in Europe

With the advent of multimodality therapy, the overall five-year survival rate from childhood cancer has improved considerably now exceeding 80% in developed European countries. This growing cohort of survivors, with many years of life ahead of them, has raised the necessity for knowledge concerning the risks of adverse long-term sequelae of the life-saving treatments in order to provide optimal screening and care and to identify and provide adequate interventions. Childhood cancer survivor cohorts in Europe. Considerable advantages exist to study late effects in individuals treated for childhood cancer in a European context, including the complementary advantages of large population-based cancer registries and the unrivalled opportunities to study lifetime risks, together with rich and detailed hospital-based cohorts which fill many of the gaps left by the large-scale population-based studies, such as sparse treatment information. Several large national cohorts have been established within Europe to study late effects in individuals treated for childhood cancer including the Nordic Adult Life after Childhood Cancer in Scandinavia study (ALiCCS), the British Childhood Cancer Survivor Study (BCCSS), the Dutch Childhood Oncology Group (DCOG) LATER study, and the Swiss Childhood Cancer Survivor Study (SCCSS). Furthermore, there are other large cohorts, which may eventually become national in scope including the French Childhood Cancer Survivor Study (FCCSS), the French Childhood Cancer Survivor Study for Leukaemia (LEA), and the Italian Study on off-therapy Childhood Cancer Survivors (OTR). In recent years significant steps have been taken to extend these national studies into a larger pan-European context through the establishment of two large consortia - PanCareSurFup and PanCareLIFE. The purpose of this paper is to present an overview of the current large, national and pan-European studies of late effects after childhood cancer. This overview will highlight the strong cooperation across Europe, in particular the EU-funded collaborative research projects PanCareSurFup and PanCareLIFE. Overall goal. The overall goal of these large cohort studies is to provide every European childhood cancer survivor with better care and better long-term health so that they reach their full potential, and to the degree possible, enjoy the same quality of life and opportunities as their peers.

Use of synchrotron medical microbeam irradiation to investigate radiation-induced bystander and abscopal effects in vivo

The question of whether bystander and abscopal effects are the same is unclear. Our experimental system enables us to address this question by allowing irradiated organisms to partner with unexposed individuals. Organs from both animals and appropriate sham and scatter dose controls are tested for expression of several endpoints such as calcium flux, role of 5HT, reporter assay cell death and proteomic profile. The results show that membrane related functions of calcium and 5HT are critical for true bystander effect expression. Our original inter-animal experiments used fish species whole body irradiated with low doses of X-rays, which prevented us from addressing the abscopal effect question. Data which are much more relevant in radiotherapy are now available for rats which received high dose local irradiation to the implanted right brain glioma. The data were generated using quasi-parallel microbeams at the biomedical beamline at the European Synchrotron Radiation Facility in Grenoble France. This means we can directly compare abscopal and "true" bystander effects in a rodent tumour model. Analysis of right brain hemisphere, left brain and urinary bladder in the directly irradiated animals and their unirradiated partners strongly suggests that bystander effects (in partner animals) are not the same as abscopal effects (in the irradiated animal). Furthermore, the presence of a tumour in the right brain alters the magnitude of both abscopal and bystander effects in the tissues from the directly irradiated animal and in the unirradiated partners which did not contain tumours, meaning the type of signal was different.

Joint research towards a better radiation protection-highlights of the Fifth MELODI Workshop

MELODI is the European platform dedicated to low-dose radiation risk research. From 7 October through 10 October 2013 the Fifth MELODI Workshop took place in Brussels, Belgium. The workshop offered the opportunity to 221 unique participants originating from 22 countries worldwide to update their knowledge and discuss radiation research issues through 118 oral and 44 poster presentations. In addition, the MELODI 2013 workshop was reaching out to the broader radiation protection community, rather than only the low-dose community, with contributions from the fields of radioecology, emergency and recovery preparedness, and dosimetry. In this review, we summarise the major scientific conclusions of the workshop, which are important to keep the MELODI strategic research agenda up-to-date and which will serve to establish a joint radiation protection research roadmap for the future.

A collimated detection system for assessing leakage dose from medical linear accelerators at the patient plane

Leakage radiation from linear accelerators can make a significant contribution to healthy tissue dose in patients undergoing radiotherapy. In this work thermoluminescent dosimeters (LiF:Mg,Cu,P TLD chips) were used in a focused lead cone loaded with TLD chips for the purpose of evaluating leakage dose at the patient plane. By placing the TLDs at one end of a stereotactic cone, a focused measurement device is created; this was tested both in and out of the primary beam of a Varian 21-iX linac using 6 MV photons. Acrylic build up material of 1.2 cm thickness was used inside the cone and measurements made with either one or three TLD chips at a given distance from the target. Comparing the readings of three dosimeters in one plane inside the cone offered information regarding the orientation of the cone relative to a radiation source. Measurements in the patient plane with the linac gantry at various angles demonstrated that leakage dose was approximately 0.01% of the primary beam out of field when the cone was pointed directly towards the target and 0.0025% elsewhere (due to scatter within the gantry). No specific 'hot spots' (e.g., insufficient shielding or gaps at abutments) were observed. Focused cone measurements facilitate leakage dose measurements from the linac head directly at the patient plane and allow one to infer the fraction of leakage due to 'direct' photons (along the ray-path from the bremsstrahlung target) and that due to scattered photons.