Roadmap for precision preclinical x-ray radiation studies.

This Roadmap paper covers the field of precision preclinical x-ray radiation studies in animal models. It is mostly focused on models for cancer and normal tissue response to radiation, but also discusses other disease models. The recent technological evolution in imaging, irradiation, dosimetry and monitoring that have empowered these kinds of studies is discussed, and many developments in the near future are outlined. Finally, clinical translation and reverse translation are discussed.

A novel data management platform to improve image-guided precision preclinical biological research.

OBJECTIVE:: Preclinical biological research is mandatory for developing new drugs to investigate the toxicity and efficacy of the drug. In this paper, the focus is on radiobiological research as an example of advanced preclinical biological research. In radiobiology, recent technological advances have produced novel research platforms which can precisely irradiate targets in animals and use advanced onboard image-guidance, mimicking the clinical radiotherapy environment. These platforms greatly facilitate complex research combining several agents simultaneously (in our example, radiation and non-radiation agents). Since these modern platform can produce a large amount of wide-ranging data, one of the main impediments in preclinical research platforms is a proper data management system for preclinical studies. METHODS:: A preclinical data management system, inspired by current radiotherapy clinical data management systems was designed. The system was designed with InterSystems technology, i.e. a programmable Enterprise Service Bus solution. New DICOM animal imaging standards are used such as DICOM suppl. 187 for storing small animal acquisition context and the DICOM second generation course model. RESULTS:: A small animal big data warehouse environment for research is designed to work with modern image-guided precision research platforms. Its modular design includes (1) a study workflow manager, (2) a data manager, and (3) a storage manager. The system provides interfaces to, e.g. preclinical treatment planning systems and data analysis plug-ins, and guides the user efficiently through the many steps involved in preclinical research. The system manages various data source locations, and arranges access to the data centrally. CONCLUSION:: A novel preclinical data management system can be designed to improve preclinical workflow, facilitate data exchange between researchers, and support translation to clinical trials. ADVANCES IN KNOWLEDGE:: A preclinical data management system such as the one proposed here would greatly benefit preparation, execution and analysis of biological experiments, and will eventually facilitate translation to clinical trials.

Time-resolved versus time-integrated portal dosimetry: the role of an object's position with respect to the isocenter in volumetric modulated arc therapy.

The aim of this work is to compare time-resolved (TR) and time-integrated (TI) portal dosimetry, focussing on the role of an object's position with respect to the isocenter in volumetric modulated arc therapy (VMAT). Portal dose images (PDIs) are simulated and measured for different cases: a sphere (1), a bovine bone (2) and a patient geometry (3). For the simulated case (1) and the experimental case (2), several transformations are applied at different off-axis positions. In the patient case (3), three simple plans with different isocenters are created and pleural effusion is simulated in the patient. The PDIs before and after the sphere transformations, as well as the PDIs with and without simulated pleural effusion, are compared using a TI and TR gamma analysis. In addition, the performance of the TI and TR gamma analyses for the detection of real geometric changes in patients treated with clinical plans is investigated and a correlation analysis is performed between gamma fail rates and differences in dose volume histogram (DVH) metrics. The TI gamma analysis can show large differences in gamma fail rates for the same transformation at different off-axis positions (or for different plan isocenters). The TR gamma analysis, however, shows consistent gamma fail rates. For the detection of real geometric changes in patients treated with clinical plans, the TR gamma analysis has a higher sensitivity than the TI gamma analysis. However, the specificity for the TR gamma analysis is lower than for the TI gamma analysis. Both the TI and TR gamma fail rates show no correlation with changes in DVH metrics. This work shows that TR portal dosimetry is fundamentally superior to TI portal dosimetry, because it removes the strong dependence of the gamma fail rate on the off-axis position/plan isocenter. However, for 2D TR portal dosimetry, it is still difficult to interpret gamma fail rates in terms of changes in DVH metrics for patients treated with VMAT.

The influence of gastric filling instructions on dose delivery in patients with oesophageal cancer: A prospective study.

PURPOSE: To evaluate whether adaptive radiotherapy for unaccounted stomach changes in patients with adenocarcinoma of the gastroesophageal junction (GEJ) is necessary and whether dose differences could be prevented by giving patients food and fluid instructions before treatment simulation and radiotherapy. MATERIAL AND METHODS: Twenty patients were randomly assigned into two groups: patients with and without instructions about restricting food and fluid intake prior to radiotherapy simulation and treatment. Redelineation and offline recalculation of dose distributions based on cone-beam computed tomography (n=100) were performed. Dose-volume parameters were analysed for the clinical target volume extending into the stomach. RESULTS: Four patients who did not receive instructions had a geometric miss (0.7-12 cm(3)) in only one fraction. With instructions, 3 out of 10 patients had a geometric miss (0.1-1.9 cm(3)) in one (n=2) or two (n=1) fractions. The V95% was reduced by more than 5% for one patient, but this underdosage was in an in-air region without further clinical importance. CONCLUSIONS: Giving patients food and fluid instructions for the treatment of GEJ cancer offers no clinical benefit. Using a planning target volume margin of 1cm implies that there is no need for adaptive radiotherapy for GEJ tumours.

Weekly kilovoltage cone-beam computed tomography for detection of dose discrepancies during (chemo)radiotherapy for head and neck cancer.

BACKGROUND: Use of highly conformal radiotherapy in patients with head and neck carcinoma may lead to under-/overdosage of gross target volume (GTV) and organs at risk (OAR) due to changes in patients' anatomy. A method to achieve more effective radiation treatment combined with less toxicity is dose-guided radiotherapy (DGRT). The aim of this study was to evaluate discrepancies between planned and actually delivered radiation dose in head and neck patients and to identify predictive factors. METHODS: In this retrospective analysis, 20 patients with cT2-4 N0-3 M0 carcinoma originating from oropharynx, oral cavity, larynx and hypopharynx (Cohort 1), and seven patients with cT1-4 N0-3 M0 nasopharyngeal carcinoma (Cohort 2) treated with primary (chemo)radiotherapy and undergoing weekly kV-CBCT scans were included. Radiation dose was recalculated on 184 kV-CBCT images, which was quantified by D95% (GTV), Dmean (parotid and submandibular glands) and D2% (spinal cord). Predictive factors investigated for changes in these dose metrics were: gender, age, cT/N-stage, tumor grade, HPV-status, systemic therapy, body mass index at start of treatment, weight loss and volume change over the duration of the radiotherapy. RESULTS: There was no significant difference between the planned and delivered dose for GTV and OARs of Week 1 to subsequent weeks for Cohort 1. In Cohort 2, actually delivered Dmean to parotid glands was significant higher than planned dose (1.1 Gy, p = 0.002). No clinically relevant correlations between dose changes and predictive factors were found. CONCLUSION: Weekly dose calculations do not seem to improve dose delivery for patients with tumors of the oral cavity, oropharynx, larynx and hypopharynx. In patients with nasopharyngeal carcinoma, however, mid-treatment imaging may facilitate DGRT.

Is integrated transit planar portal dosimetry able to detect geometric changes in lung cancer patients treated with volumetric modulated arc therapy?

BACKGROUND: Geometric changes are frequent during the course of treatment of lung cancer patients. This may potentially result in deviations between the planned and actual delivered dose. Electronic portal imaging device (EPID)-based integrated transit planar portal dosimetry (ITPD) is a fast method for absolute in-treatment dose verification. The aim of this study was to investigate if ITPD could detect geometric changes in lung cancer patients. MATERIALS AND METHODS: A total of 460 patients treated with volumetric modulated arc therapy (VMAT) following daily cone beam computed tomography (CT)-based setup were visually inspected for geometrical changes on a daily basis. Forty-six patients were subject to changes and had a re-CT and an adaptive treatment plan. The reasons for adaptation were: change in atelectasis (n = 18), tumor regression (n = 9), change in pleural effusion (n = 8) or other causes (n = 11). The ITPDs were calculated on both the initial planning CT and the re-CT and compared with a global gamma (γ) evaluation (criteria: 3%\3mm). A treatment fraction failed when the percentage of pixels failing in the radiation fields exceeded 10%. Dose-volume histograms (DVHs) were compared between the initial plan versus the plan re-calculated on the re-CT. RESULTS: The ITPD threshold method detected 76% of the changes in atelectasis, while only 50% of the tumor regression cases and 42% of the pleural effusion cases were detected. Only 10% of the cases adapted for other reasons were detected with ITPD. The method has a 17% false-positive rate. No significant correlations were found between changes in DVH metrics and γ fail-rates. CONCLUSIONS: This study showed that most cases with geometric changes caused by atelectasis could be captured by ITPD, however for other causes ITPD is not sensitive enough to detect the clinically relevant changes and no predictive power of ITPD was found.

What level of accuracy is achievable for preclinical dose painting studies on a clinical irradiation platform?

Advancements made over the past decades in both molecular imaging and radiotherapy planning and delivery have enabled studies that explore the efficacy of heterogeneous radiation treatment ("dose painting") of solid cancers based on biological information provided by different imaging modalities. In addition to clinical trials, preclinical studies may help contribute to identifying promising dose painting strategies. The goal of this current study was twofold: to develop a reproducible positioning and set-up verification protocol for a rat tumor model to be imaged and treated on a clinical platform, and to assess the dosimetric accuracy of dose planning and delivery for both uniform and positron emission tomography-computed tomography (PET-CT) based heterogeneous dose distributions. We employed a syngeneic rat rhabdomyosarcoma model, which was irradiated by volumetric modulated arc therapy (VMAT) with uniform or heterogeneous 6 MV photon dose distributions. Mean dose to the gross tumor volume (GTV) as a whole was kept at 12 Gy for all treatment arms. For the nonuniform plans, the dose was redistributed to treat the 30% of the GTV representing the biological target volume (BTV) with a dose 40% higher than the rest of the GTV (GTV - BTV) (~15 Gy was delivered to the BTV vs. ~10.7 Gy was delivered to the GTV - BTV). Cone beam computed tomography (CBCT) images acquired for each rat prior to irradiation were used to correctly reposition the tumor and calculate the delivered 3D dose. Film quality assurance was performed using a water-equivalent rat phantom. A comparison between CT or CBCT doses and film measurements resulted in passing rates >98% with a gamma criterion of 3%/2 mm using 2D dose images. Moreover, between the CT and CBCT calculated doses for both uniform and heterogeneous plans, we observed maximum differences of <2% for mean dose to the tumor and mean dose to the biological target volumes. In conclusion, we have developed a robust method for dose painting in a rat tumor model on a clinical platform, with a high accuracy achieved in the delivery of complex dose distributions. Our work demonstrates the technical feasibility of this approach and enables future investigations on the therapeutic effect of preclinical dose painting strategies using a state-of-the-art clinical platform.

A novel time dependent gamma evaluation function for dynamic 2D and 3D dose distributions.

Modern external beam radiotherapy requires detailed verification and quality assurance so that confidence can be placed on both the delivery of a single treatment fraction and on the consistency of delivery throughout the treatment course. To verify dose distributions, a comparison between prediction and measurement must be made. Comparisons between two dose distributions are commonly performed using a Gamma evaluation which is a calculation of two quantities on a pixel by pixel basis; the dose difference, and the distance to agreement. By providing acceptance criteria (e.g. 3%, 3 mm), the function will find the most appropriate match within its two degrees of freedom. For complex dynamic treatments such as IMRT or VMAT it is important to verify the dose delivery in a time dependent manner and so a gamma evaluation that includes a degree of freedom in the time domain via a third parameter, time to agreement, is presented here. A C++ (mex) based gamma function was created that could be run on either CPU and GPU computing platforms that would allow a degree of freedom in the time domain. Simple test cases were created in both 2D and 3D comprising of simple geometrical shapes with well-defined boundaries varying over time. Changes of varying magnitude in either space or time were introduced and repeated gamma analyses were performed varying the criteria. A clinical VMAT case was also included, artificial air bubbles of varying size were introduced to a patient geometry, along with shifts of varying magnitude in treatment time. For all test cases where errors in distance, dose or time were introduced, the time dependent gamma evaluation could accurately highlight the errors.The time dependent gamma function presented here allows time to be included as a degree of freedom in gamma evaluations. The function allows for 2D and 3D data sets which are varying over time to be compared using appropriate criteria without penalising minor offsets of subsequent radiation fields in time.

Time dependent pre-treatment EPID dosimetry for standard and FFF VMAT.

Methods to calibrate Megavoltage electronic portal imaging devices (EPIDs) for dosimetry have been previously documented for dynamic treatments such as intensity modulated radiotherapy (IMRT) using flattened beams and typically using integrated fields. While these methods verify the accumulated field shape and dose, the dose rate and differential fields remain unverified. The aim of this work is to provide an accurate calibration model for time dependent pre-treatment dose verification using amorphous silicon (a-Si) EPIDs in volumetric modulated arc therapy (VMAT) for both flattened and flattening filter free (FFF) beams. A general calibration model was created using a Varian TrueBeam accelerator, equipped with an aS1000 EPID, for each photon spectrum 6 MV, 10 MV, 6 MV-FFF, 10 MV-FFF. As planned VMAT treatments use control points (CPs) for optimization, measured images are separated into corresponding time intervals for direct comparison with predictions. The accuracy of the calibration model was determined for a range of treatment conditions. Measured and predicted CP dose images were compared using a time dependent gamma evaluation using criteria (3%, 3 mm, 0.5 sec). Time dependent pre-treatment dose verification is possible without an additional measurement device or phantom, using the on-board EPID. Sufficient data is present in trajectory log files and EPID frame headers to reliably synchronize and resample portal images. For the VMAT plans tested, significantly more deviation is observed when analysed in a time dependent manner for FFF and non-FFF plans than when analysed using only the integrated field. We show EPID-based pre-treatment dose verification can be performed on a CP basis for VMAT plans. This model can measure pre-treatment doses for both flattened and unflattened beams in a time dependent manner which highlights deviations that are missed in integrated field verifications.

Towards a semantic PACS: Using Semantic Web technology to represent imaging data.

The DICOM standard is ubiquitous within medicine. However, improved DICOM semantics would significantly enhance search operations. Furthermore, databases of current PACS systems are not flexible enough for the demands within image analysis research. In this paper, we investigated if we can use Semantic Web technology, to store and represent metadata of DICOM image files, as well as linking additional computational results to image metadata. Therefore, we developed a proof of concept containing two applications: one to store commonly used DICOM metadata in an RDF repository, and one to calculate imaging biomarkers based on DICOM images, and store the biomarker values in an RDF repository. This enabled us to search for all patients with a gross tumor volume calculated to be larger than 50 cc. We have shown that we can successfully store the DICOM metadata in an RDF repository and are refining our proof of concept with regards to volume naming, value representation, and the applications themselves.

First clinical results of adaptive radiotherapy based on 3D portal dosimetry for lung cancer patients with atelectasis treated with volumetric-modulated arc therapy (VMAT).

UNLABELLED: Atelectasis in lung cancer patients can change rapidly during a treatment course, which may displace the tumor/healthy tissues, or change tissue densities locally. This may result in differences between the planned and the actually delivered dose. With complex delivery techniques treatment verification is essential and inter-fractional adaptation may be necessary. We present the first clinical results of treatment adaptation based on an in-house developed three-dimensional (3D) portal dose measurement (PDM) system. MATERIAL AND METHODS: A method was developed for 3D PDM combined with cone beam computed tomography (kV-CBCT) imaging. Lung cancer patients are monitored routinely with this imaging technique. During treatment, the first three fractions are analyzed with 3D PDM and weekly thereafter. The reconstructed measured dose is compared to the planned dose using dose-volume histograms and a γ evaluation. Patients having |γ|> 1 in more than 5% of the (primary tumor or organ at risk) volume were subjected to further analysis. In this study we show the PDM dose changes for five patients. RESULTS: We detected relevant dose changes induced by changes in atelectasis in the presented cases. Two patients received two treatment adaptations after being detected with PDM confirmed by visual inspection of the kV-CBCTs, and in two other patients the radiation treatment plan was adapted once. In one case no dose delivery change was detected with PDM. CONCLUSION: The first clinical patients show that 3D PDM combined with kV-CBCT is a valuable quality assurance tool for detecting anatomical alterations and their dosimetric consequences during the course of radiotherapy. In our clinic, 3D PDM is fully automated for ease and speed of the procedure, and for minimization of human error. The technique is able to flag patients with suspected dose discrepancies for potential adaptation of the treatment plan.

Benefits of a clinical data warehouse with data mining tools to collect data for a radiotherapy trial.

INTRODUCTION: Collecting trial data in a medical environment is at present mostly performed manually and therefore time-consuming, prone to errors and often incomplete with the complex data considered. Faster and more accurate methods are needed to improve the data quality and to shorten data collection times where information is often scattered over multiple data sources. The purpose of this study is to investigate the possible benefit of modern data warehouse technology in the radiation oncology field. MATERIAL AND METHODS: In this study, a Computer Aided Theragnostics (CAT) data warehouse combined with automated tools for feature extraction was benchmarked against the regular manual data-collection processes. Two sets of clinical parameters were compiled for non-small cell lung cancer (NSCLC) and rectal cancer, using 27 patients per disease. Data collection times and inconsistencies were compared between the manual and the automated extraction method. RESULTS: The average time per case to collect the NSCLC data manually was 10.4 ± 2.1 min and 4.3 ± 1.1 min when using the automated method (p<0.001). For rectal cancer, these times were 13.5 ± 4.1 and 6.8 ± 2.4 min, respectively (p<0.001). In 3.2% of the data collected for NSCLC and 5.3% for rectal cancer, there was a discrepancy between the manual and automated method. CONCLUSIONS: Aggregating multiple data sources in a data warehouse combined with tools for extraction of relevant parameters is beneficial for data collection times and offers the ability to improve data quality. The initial investments in digitizing the data are expected to be compensated due to the flexibility of the data analysis. Furthermore, successive investigations can easily select trial candidates and extract new parameters from the existing databases.

Results of a multicentric in silico clinical trial (ROCOCO): comparing radiotherapy with photons and protons for non-small cell lung cancer.

INTRODUCTION: This multicentric in silico trial compares photon and proton radiotherapy for non-small cell lung cancer patients. The hypothesis is that proton radiotherapy decreases the dose and the volume of irradiated normal tissues even when escalating to the maximum tolerable dose of one or more of the organs at risk (OAR). METHODS: Twenty-five patients, stage IA-IIIB, were prospectively included. On 4D F18-labeled fluorodeoxyglucose-positron emission tomography-computed tomography scans, the gross tumor, clinical and planning target volumes, and OAR were delineated. Three-dimensional conformal radiotherapy (3DCRT) and intensity-modulated radiotherapy (IMRT) photon and passive scattered conformal proton therapy (PSPT) plans were created to give 70 Gy to the tumor in 35 fractions. Dose (de-)escalation was performed by rescaling to the maximum tolerable dose. RESULTS: Protons resulted in the lowest dose to the OAR, while keeping the dose to the target at 70 Gy. The integral dose (ID) was higher for 3DCRT (59%) and IMRT (43%) than for PSPT. The mean lung dose reduced from 18.9 Gy for 3DCRT and 16.4 Gy for IMRT to 13.5 Gy for PSPT. For 10 patients, escalation to 87 Gy was possible for all 3 modalities. The mean lung dose and ID were 40 and 65% higher for photons than for protons, respectively. CONCLUSIONS: The treatment planning results of the Radiation Oncology Collaborative Comparison trial show a reduction of ID and the dose to the OAR when treating with protons instead of photons, even with dose escalation. This shows that PSPT is able to give a high tumor dose, while keeping the OAR dose lower than with the photon modalities.

A fast three-dimensional gamma evaluation using a GPU utilizing texture memory for on-the-fly interpolations.

PURPOSE: A widely accepted method to quantify differences in dose distributions is the gamma (gamma) evaluation. Currently, almost all gamma implementations utilize the central processing unit (CPU). Recently, the graphics processing unit (GPU) has become a powerful platform for specific computing tasks. In this study, we describe the implementation of a 3D gamma evaluation using a GPU to improve calculation time. METHODS: The gamma evaluation algorithm was implemented on an NVIDIA Tesla C2050 GPU using the compute unified device architecture (CUDA). First, several cubic virtual phantoms were simulated. These phantoms were tested with varying dose cube sizes and set-ups, introducing artificial dose differences. Second, to show applicability in clinical practice, five patient cases have been evaluated using the 3D dose distribution from a treatment planning system as the reference and the delivered dose determined during treatment as the comparison. A calculation time comparison between the CPU and GPU was made with varying thread-block sizes including the option of using texture or global memory. RESULTS: A GPU over CPU speed-up of 66 +/- 12 was achieved for the virtual phantoms. For the patient cases, a speed-up of 57 +/- 15 using the GPU was obtained. A thread-block size of 16 x 16 performed best in all cases. The use of texture memory improved the total calculation time, especially when interpolation was applied. Differences between the CPU and GPU gammas were negligible. CONCLUSIONS: The GPU and its features, such as texture memory, decreased the calculation time for gamma evaluations considerably without loss of accuracy.

Design of and technical challenges involved in a framework for multicentric radiotherapy treatment planning studies.

This report introduces a framework for comparing radiotherapy treatment planning in multicentric in silico clinical trials. Quality assurance, data incompatibility, transfer and storage issues, and uniform analysis of results are discussed. The solutions that are given provide a useful guide for the set-up of future multicentric planning studies or public repositories of high quality data.