Sensitivity and specificity of Monte Carlo based independent secondary dose computation for detecting modulation-related dose errors in intensity modulated radiotherapy.

BACKGROUND: The recent availability of Monte Carlo based independent secondary dose calculation (ISDC) for patient-specific quality assurance (QA) of modulated radiotherapy requires the definition of appropriate, more sensitive action levels, since contemporary recommendations were defined for less accurate ISDC dose algorithms. PURPOSE: The objective is to establish an optimum action level and measure the efficacy of a Monte Carlo ISDC software for pre-treatment QA of intensity modulated radiotherapy treatments. METHODS: The treatment planning system and the ISDC were commissioned by their vendors from independent base data sets, replicating a typical real-world scenario. In order to apply Receiver-Operator-Characteristics (ROC), a set of treatment plans for various case classes was created that consisted of 190 clinical treatment plans and 190 manipulated treatment plans with dose errors in the range of 1.5-2.5%. All 380 treatment plans were evaluated with ISDC in the patient geometry. ROC analysis was performed for a number of Gamma (dose-difference/distance-to-agreement) criteria. QA methods were ranked according to Area under the ROC curve (AUC) and optimum action levels were derived via Youden's J statistics. RESULTS: Overall, for original treatment plans, the mean Gamma pass rate (GPR) for Gamma(1%, 1 mm) was close to 90%, although with some variation across case classes. The best QA criterion was Gamma(2%, 1 mm) with GPR > 90% and an AUC of 0.928. Gamma criteria with small distance-to-agreement had consistently higher AUC. GPR of original treatment plans depended on their modulation degree. An action level in terms of Gamma(1%, 1 mm) GPR that decreases with modulation degree was the most efficient criterion with sensitivity = 0.91 and specificity = 0.95, compared with Gamma(3%, 3 mm) GPR > 99%, sensitivity = 0.73 and specificity = 0.91 as a commonly used action level. CONCLUSIONS: ISDC with Monte Carlo proves highly efficient to catch errors in the treatment planning process. For a Monte Carlo based TPS, dose-difference criteria of 2% or less, and distance-to-agreement criteria of 1 mm, achieve the largest AUC in ROC analysis.

Statistical breathing curve sampling to quantify interplay effects of moving lung tumors in a 4D Monte Carlo dose calculation framework.

PURPOSE: The interplay between respiratory tumor motion and dose application by intensity modulated radiotherapy (IMRT) techniques can potentially lead to undesirable and non-intuitive deviations from the planned dose distribution. We developed a 4D Monte Carlo (MC) dose recalculation framework featuring statistical breathing curve sampling, to precisely simulate the dose distribution for moving target volumes aiming at a comprehensive assessment of interplay effects. METHODS: We implemented a dose accumulation tool that enables dose recalculations of arbitrary breathing curves including the actual breathing curve of the patient. This MC dose recalculation framework is based on linac log-files, facilitating a high temporal resolution up to 0.1 s. By statistical analysis of 128 different breathing curves, interplay susceptibility of different treatment parameters was evaluated for an exemplary patient case. To facilitate prospective clinical application in the treatment planning stage, in which patient breathing curves or linac log-files are not available, we derived a log-file free version with breathing curves generated by a random walk approach. Interplay was quantified by standard deviations σ in D5%, D50% and D95%. RESULTS: Interplay induced dose deviations for single fractions were observed and evaluated for IMRT and volumetric arc therapy (σD95% up to 1.3 %) showing a decrease with higher fraction doses and an increase with higher MU rates. Interplay effects for conformal treatment techniques were negligible (σ<0.1%). The log-file free version and the random walk generated breathing curves yielded similar results (deviations in σ< 0.1 %) and can be used as substitutes for interplay assessment. CONCLUSION: It is feasible to combine statistically sampled breathing curves with MC dose calculations. The universality of the presented framework allows comprehensive assessment of interplay effects in retrospective and prospective clinically relevant scenarios.

Commissioning and clinical implementation of the first commercial independent Monte Carlo 3D dose calculation to replace CyberKnife M6™ patient-specific QA measurements.

PURPOSE: To report on the commissioning and clinical validation of the first commercially available independent Monte Carlo (MC) three-dimensional (3D) dose calculation for CyberKnife robotic radiosurgery system® (Accuray, Sunnyvale, CA). METHODS: The independent dose calculation (IDC) by SciMoCa® (Scientific RT, Munich, Germany) was validated based on water measurements of output factors and dose profiles (unshielded diode, field-size dependent corrections). A set of 84 patient-specific quality assurance (QA) measurements for multi-leaf collimator (MLC) plans, using an Octavius two-dimensional SRS1000 array (PTW, Freiburg, Germany), was compared to results of respective calculations. Statistical process control (SPC) was used to detect plans outside action levels. RESULTS: Of all output factors for the three collimator systems of the CyberKnife, 99% agreed within 2% and 81% within 1%, with a maximum deviation of 3.2% for a 5-mm fixed cone. The profiles were compared using a one-dimensional gamma evaluation with 2% dose difference and 0.5 mm distance-to-agreement (Γ(2,0.5)). The off-centre ratios showed an average pass rate >99% (92-100%). The agreement of the depth dose profiles depended on field size, with lowest pass rates for the smallest MLC field sizes. The average depth dose pass rate was 88% (35-99%). The IDCs showed a Γ(2,1) pass rate of 98%. Statistical process control detected six plans outside tolerance levels in the measurements, all of which could be attributed the measurement setup. Independent dose calculations showed problems in five plans, all due to differences in the algorithm between TPS and IDC. Based on these results changes were made in the class solution for treatment plans. CONCLUSION: The first commercially available MC 3D dose IDC was successfully commissioned and validated for the CyberKnife and replaced all routine patient-specific QA measurements in our clinic.

Volumetric and actuarial analysis of brain necrosis in proton therapy using a novel mixture cure model.

BACKGROUND AND PURPOSE: High-dose fractionated radiotherapy is often necessary to achieve long-term tumor control in several types of tumors involving or within close proximity to the brain. There is limited data to guide on optimal constraints to the adjacent nontarget brain. This investigation explored the significance of the three-dimensional (3D) dose distribution of passive scattering proton therapy to the brain with other clinicopathological factors on the development of symptomatic radiation necrosis. MATERIALS AND METHODS: All patients with head and neck, skull base, or intracranial tumors who underwent proton therapy (minimum prescription dose of 59.4 Gy(RBE)) with collateral moderate to high dose radiation exposure to the nontarget brain were retrospectively reviewed. A mixture cure model with respect to necrosis-free survival was used to derive estimates for the normal tissue complication probability (NTCP) model while adjusting for potential confounding factors. RESULTS: Of 179 identified patients, 83 patients had intracranial tumors and 96 patients had primary extracranial tumors. The optimal dose measure obtained to describe the occurrence of radiation necrosis was the equivalent uniform dose (EUD) with parameter a = 9. The best-fit parameters of logistic NTCP models revealed D50 = 57.7 Gy for intracranial tumors, D50 = 39.5 Gy for extracranial tumors, and γ50 = 2.5 for both tumor locations. Multivariable analysis revealed EUD and primary tumor location to be the strongest predictors of brain radiation necrosis. CONCLUSION: In the current clinical volumetric data analyses with multivariable modelling, EUD was identified as an independent and strong predictor for brain radiation necrosis from proton therapy.

Comparison of planned dose on different CT image sets to four-dimensional Monte Carlo dose recalculation using the patient's actual breathing trace for lung stereotactic body radiation therapy.

PURPOSE: The need for four-dimensional (4D) treatment planning becomes indispensable when it comes to radiation therapy for moving tumors in the thoracic and abdominal regions. The primary purpose of this study is to combine the actual breathing trace during each individual treatment fraction with the Linac's log file information and Monte Carlo 4D dose calculations. We investigated this workflow on multiple computed tomography (CT) datasets in a clinical environment for stereotactic body radiation therapy (SBRT) treatment planning. METHODS: We have developed a workflow, which allows us to recalculate absorbed dose to a 4DCT dataset using Monte Carlo calculation methods and accumulate all 4D doses in order to compare them to the planned dose using the Linac's log file, a 4DCT dataset, and the patient's actual breathing curve for each individual fraction. For five lung patients, three-dimensional-conformal radiation therapy (3D-CRT) and volumetric modulated arc treatment (VMAT) treatment plans were generated on four different CT image datasets: a native free-breathing 3DCT, an average intensity projection (AIP) and a maximum intensity projection (MIP) CT both obtained from a 4DCT, and a 3DCT with density overrides based on the 3DCT (DO). The Monte Carlo 4D dose has been calculated on each 4DCT phase using the Linac's log file and the patient's breathing trace as a surrogate for tumor motion and dose was accumulated to the gross tumor volume (GTV) at the 50% breathing phase (end of exhale) using deformable image registration. RESULTS: Δ D 98 % and Δ D 2 % between 4D dose and planned dose differed largely for 3DCT-based planning and also for DO in three patients. Least dose differences between planned and recalculated dose have been found for AIP and MIP treatment planning which both tend to be superior to DO, but the results indicate a dependency on the breathing variability, tumor motion, and size. An interplay effect has not been observed in the small patient cohort. CONCLUSIONS: We have developed a workflow which, to our best knowledge, is the first incorporation of the patient breathing trace over the course of all individual treatment fractions with the Linac's log file information and 4D Monte Carlo recalculations of the actual treated dose. Due to the small patient cohort, no clear recommendation on which CT can be used for SBRT treatment planning can be given, but the developed workflow, after adaption for clinical use, could be used to enhance a priori 4D Monte Carlo treatment planning in the future and help with the decision on which CT dataset treatment planning should be carried out.

Statistical motion modelling for robust evaluation of clinically delivered accumulated dose distributions after curative radiotherapy of locally advanced prostate cancer.

BACKGROUND AND PURPOSE: Planned doses are used as surrogate for the actually delivered dose in radiotherapy. We have estimated the delivered dose in a dose-escalation trial of locally advanced prostate cancer by statistical dose-accumulation and by DVH-summation, and compared to planned dose. MATERIALS AND METHOD: Prescribed dose-escalation to the prostate was 67.5 Gy/25fr., corresponding to 81GyEQD2 assuming α/ß = 1.5. The 21 patients had three targets (i.e. CTV67.5 + 2 mm, CTV60 + 5 mm, CTV50 + 10 mm) irradiated by a simultaneous-integrated-boost technique. Analysis was based on 213 CT scans and 5-years of follow-up. For statistical dose-accumulation, we modelled 10000 possible treatment courses based on planned dose and deformation-vector-fields from contour-based registration. For DVH-summation we recalculated dose on repeat-CTs and estimated median D98%/EUD. Groups with/without disease recurrence were compared. RESULTS: Discrepancies between planned and accumulated dose were mostly seen for CTV67.5, where under-dosage was found at different locations in the prostate in 12/21 patients. Delivered dose-escalation (D98%) was on average 73.9GyEQD2 (range: 68.3-78.7GyEQD2). No significant difference in accumulated-D98% was found in patients with (n = 8) and without (n = 13) recurrence (p > 0.05). Average D98%/EUD with statistical dose-accumulation vs DVH-summation was significantly different in CTV60, CTV50, rectum and bladder but not in CTV67.5. CONCLUSION: The planned dose escalation was not received by more than half-of-the patients. Robustness of the prostate target (CTV67.5) should therefore be better prioritized in these patients given the low toxicity profile. Estimates of delivered dose were less conservative for dose-accumulation due to interaction of random organ motion with the dose matrix.

Evaluation of machine log files/MC-based treatment planning and delivery QA as compared to ArcCHECK QA.

PURPOSE: A treatment planning/delivery QA tool using linac log files (LF) and Monte Carlo (MC) dose calculation is investigated as a standalone alternative to phantom-based patient-specific QA (ArcCHECK (AC)). METHODS: Delivering a variety of fields onto MapCHECK2 and ArcCHECK, diode sensitivity dependence on dose rate (in-field) and energy (primarily out-of-field) was quantified. AC and LF QAs were analyzed with respect to delivery complexity by delivering 12 × 12 cm static fields/arcs comprised of varying numbers of abutting sub-fields onto ArcCHECK. About 11 clinical dual-arc VMAT patients planned using Pinnacle's convolution-superposition (CS) were delivered on ArcCHECK and log file dose (LF-CS and LF-MC) calculated. To minimize calculation time, reduced LF-CS sampling (1/2/3/4° control point spacing) was investigated. Planned ("Plan") and LF-reconstructed CS and MC doses were compared with each other and AC measurement via statistical [mean ± StdDev(σ)] and gamma analyses to isolate dosimetric uncertainties and quantify the relative accuracies of AC QA and MC-based LF QA. RESULTS: Calculation and ArcCHECK measurement differed by up to 1.5% in-field due to variation in dose rate and up to 5% out-of-field. For the experimental segment-varying plans, despite CS calculation deviating by as much as 13% from measurement, Plan-MC and LF-MC doses generally matched AC measurement within 3%. Utilizing 1° control point spacing, 2%/2 mm LF-CS vs AC pass rates (97%) were slightly lower than Plan-CS vs AC pass rates (97.5%). Utilizing all log file samples, 2%/2 mm LF-MC vs AC pass rates (97.3%) were higher than Plan-MC vs AC (96.5%). Phantom-dependent, calculation algorithm-dependent (MC vs CS), and delivery error-dependent dose uncertainties were 0.8 ± 1.2%, 0.2 ± 1.1%, and 0.1 ± 0.9% respectively. CONCLUSION: Reconstructing every log file sample with no increase in computational cost, MC-based LF QA is faster and more accurate than CS-based LF QA. Offering similar dosimetric accuracy compared to AC measurement, MC-based log files can be used for treatment planning QA.

Clinical workflow optimization to improve 4DCT reconstruction for Toshiba Aquilion CT scanners.

Respiratory motion remains a source of major uncertainties in radiotherapy. Respiratory correlated computed tomography (referred to as 4DCT) serves as one way of reducing breathing artifacts in 3D-CTs and allows the investigation of tumor motion over time. The quality of the 4DCT images depends on the data acquisition scheme, which in turn is dependent on the vendor. Specifically, the only way Toshiba Aquilion LB CT scanners can reconstruct 4DCTs is a cycle-based reconstruction using triggers provided by an external surrogate signal. The accuracy is strongly dependent on the method of trigger generation. Two consecutive triggers are used to define a breathing cycle which is divided into respiratory phases of equal duration. The goal of this study is to identify if there are advantages in the usage of local-amplitude based sorting (LAS) of the respiration motion states, in order to reduce image artifacts and improve 4DCT quality. Furthermore, this study addresses the generation and optimization of a clinical workflow using as surrogate motion monitoring system the Sentinel™ (C-RAD AB, Sweden) optical surface scanner in combination with a Toshiba Aquilion LB CT scanner. For that purpose, a phantom study using 10 different breathing waveforms and a retrospective patient study using the 4DCT reconstructions of 10 different patients has been conducted. The error in tumor volume has been reduced from 2.9±3.7% to 2.7±2.6% using optimal cycle-based triggers (manipulated CBS) and to 2.7±2.2% using LAS in the phantom study. Moreover, it was possible to decrease the tumor volume variability from 5.0±3.6% using the original cycle-based triggers (original CBS) to 3.5±2.5% using the optimal triggers and to 3.7±2.7% using LAS in the patient data analysis. We therefore propose the usage of the manipulated CBS, also with regard to an accurate and safe clinical workflow.

Multi-criterial patient positioning based on dose recalculation on scatter-corrected CBCT images.

BACKGROUND AND PURPOSE: Our aim was to evaluate the feasibility and potential advantages of dose guided patient positioning based on dose recalculation on scatter corrected cone beam computed tomography (CBCT) image data. MATERIAL AND METHODS: A scatter correction approach has been employed to enable dose calculations on CBCT images. A recently proposed tool for interactive multicriterial dose-guided patient positioning which uses interpolation between pre-calculated sample doses has been utilized. The workflow was retrospectively evaluated for two head and neck patients with a total of 39 CBCTs. Dose-volume histogram (DVH) parameters were compared to rigid image registration based isocenter corrections (clinical scenario). RESULTS: The accuracy of the dose interpolation was found sufficient, facilitating the implementation of dose guided patient positioning. Compared to the clinical scenario, the mean dose to the parotid glands could be improved for 2 out of 5 fractions for the first patient while other parameters were preserved. For the second patient, the mean coverage over all fractions of the high dose PTV could be improved by 4%. For this patient, coverage improvements had to be traded against organ at risk (OAR) doses within their clinical tolerance limits. CONCLUSIONS: Dose guided patient positioning using in-room CBCT data is feasible and offers increased control over target coverage and doses to OARs.

Practical implications for the quality assurance of modulated radiation therapy techniques using point detector arrays.

PURPOSE: Linac parameters potentially influencing the delivery quality of IMRT and VMAT plans are investigated with respect to threshold ranges, consequently to be considered in a linac based quality assurance procedure. Three commercially available 2D arrays are used to further investigate the influence of the measurement device. METHODS: Using three commercially available 2D arrays (Mx: MatriXXevolution , Oc: Octavius1500 , Mc: MapCHECK2), simple static measurements, measurements for MLC characterization and dynamic interplay of gantry movement, MLC movement and variable dose rate were performed. The results were evaluated with respect to each single array as well as among each other. RESULTS: Simple static measurements showed different array responses to dose, dose rate and profile homogeneity and revealed instabilities in dose delivery and profile shape during linac ramp up. Using the sweeping gap test, all arrays were able to detect small leaf misalignments down to ±0.1 mm, but this test also demonstrated up to 15% dose deviation due to profile instabilities and fast accelerating leaves during linac ramp up. Tests including gantry rotation showed different stability of gantry mounts for each array. Including gantry movement and dose rate variability, differences compared to static delivery were smaller compared to dose differences when simultaneously controling interplay of gantry movement, leaf movement and dose rate variability. CONCLUSION: Linac based QA is feasible with the tested commercially available 2D arrays. Limitations of each array and the linac ramp up characteristics should be carefully considered during individual plan generation and regularly checked in linac QA. Especially the dose and dose profile during linac ramp up should be checked regularly, as well as MLC positioning accuracy using a sweeping gap test. Additionally, dynamic interplay tests including various gantry rotation speeds and angles, various leaf speeds and various dose rates should be included.

Hippocampal sparing radiotherapy for glioblastoma patients: a planning study using volumetric modulated arc therapy.

BACKGROUND: The purpose of this study is to investigate the potential to reduce exposure of the contralateral hippocampus in radiotherapy for glioblastoma using volumetric modulated arc therapy (VMAT). METHODS: Datasets of 27 patients who had received 3D conformal radiotherapy (3D-CRT) for glioblastoma with a prescribed dose of 60Gy in fractions of 2Gy were included in this planning study. VMAT plans were optimized with the aim to reduce the dose to the contralateral hippocampus as much as possible without compromising other parameters. Hippocampal dose and treatment parameters were compared to the 3D-CRT plans using the Wilcoxon signed-rank test. The influence of tumour location and PTV size on the hippocampal dose was investigated with the Mann-Whitney-U-test and Spearman's rank correlation coefficient. RESULTS: The median reduction of the contralateral hippocampus generalized equivalent uniform dose (gEUD) with VMAT was 36 % compared to the original 3D-CRT plans (p < 0.05). Other dose parameters were maintained or improved. The median V30Gy brain could be reduced by 17.9 % (p < 0.05). For VMAT, a parietal and a non-temporal tumour localisation as well as a larger PTV size were predictors for a higher hippocampal dose (p < 0.05). CONCLUSIONS: Using VMAT, a substantial reduction of the radiotherapy dose to the contralateral hippocampus for patients with glioblastoma is feasible without compromising other treatment parameters. For larger PTV sizes, less sparing can be achieved. Whether this approach is able to preserve the neurocognitive status without compromising the oncological outcome needs to be investigated in the setting of prospective clinical trials.

Impact of MLC properties and IMRT technique in meningioma and head-and-neck treatments.

PURPOSE: The impact of multileaf collimator (MLC) design and IMRT technique on plan quality and delivery improvements for head-and-neck and meningioma patients is compared in a planning study. MATERIAL AND METHODS: Ten previously treated patients (5 head-and-neck, 5 meningioma) were re-planned for step-and-shoot IMRT (ssIMRT), sliding window IMRT (dMLC) and VMAT using the MLCi2 without (-) and with (+) interdigitation and the Agility-MLC attached to an Elekta 6MV linac. This results in nine plans per patient. Consistent patient individual optimization parameters are used. Plans are generated using the research tool Hyperion V2.4 (equivalent to Elekta Monaco 3.2) with hard constraints for critical structures and objectives for target structures. For VMAT plans, the improved segment shape optimization is used. Critical structures are evaluated based on QUANTEC criteria. PTV coverage is compared by EUD, Dmean, homogeneity and conformity. Additionally, MU/plan, treatment times and number of segments are evaluated. RESULTS: As constrained optimization is used, all plans fulfill the hard constraints. Doses to critical structures do not differ more than 1 Gy between the nine generated plans for each patient. Only larynx, parotids and eyes differ up to 1.5 Gy (Dmean or Dmax) or 7% (volume-constraint) due to (1) increased scatter, (2) not avoiding structures when using the full range of gantry rotation and (3) improved leaf sequencing with advanced segment shape optimization for VMAT plans. EUD, Dmean, homogeneity and conformity are improved using the Agility-MLC. However, PTV coverage is more affected by technique. MU increase with the use of dMLC and VMAT, while the MU are reduced by using the Agility-MLC. Fastest treatments are always achieved using Agility-MLC, especially in combination with VMAT. CONCLUSION: Fastest treatments with the best PTV coverage are found for VMAT plans with Agility-MLC, achieving the same sparing of healthy tissue compared to the other combinations of ssIMRT, dMLC and VMAT with either MLCi2(-/+) or Agility.

Intra-fractional bladder motion and margins in adaptive radiotherapy for urinary bladder cancer.

BACKGROUND: The bladder is a tumour site well suited for adaptive radiotherapy (ART) due to large inter-fractional changes, but it also displays considerable intra-fractional motion. The aim of this study was to assess target coverage with a clinically applied method for plan selection ART and to estimate population-based and patient-specific intra-fractional margins, also relevant for a future re-optimisation strategy. MATERIAL AND METHODS: Nine patients treated in a clinical phase II ART trial of daily plan selection for bladder cancer were included. In the library plans, 5 mm isotropic margins were added to account for intra-fractional changes. Pre-treatment and weekly repeat magnetic resonance imaging (MRI) series were acquired in which a full three-dimensional (3D) volume was scanned every second min for 10 min (a total of 366 scans in 61 series). Initially, the bladder clinical target volume (CTV) was delineated in all scans. The t = 0 min scan was then rigidly registered to the planning computed tomography (CT) and plan selections were simulated using the CTV_0 (at t = 0 min). To assess intra-fractional motion, coverage of the CTV_10 (at t = 10 min) was quantified using the applied PTV. Population-based margins were calculated using the van Herk margin recipe while patient-specific margins were calculated using a linear model. RESULTS: For 49% of the cases, the CTV_10 extended more than 5 mm outside the CTV_0. However, in 58 of the 61 cases (97%) CTV_10 was covered by the selected PTV. Population-based margins of 14 mm Sup/Ant, 9 mm Post and 5 mm Inf/Lat were sufficient to cover the bladder. Using patient-specific margins, the overlap between PTV and bowel-cavity was reduced from 137 cm(3) with the plan selection strategy to 24 cm(3). CONCLUSION: In this phase II ART trial, 5 mm isotropic margin for intra-fractional motion was sufficient even though considerable intra-fractional motion was observed. In online re-optimised ART, population-based margin can be applied although patient-specific margins are preferable.

Technical evaluation of different respiratory monitoring systems used for 4D CT acquisition under free breathing.

Respiratory monitoring systems are required to supply CT scanners with information on the patient's breathing during the acquisition of a respiration-correlated computer tomography (RCCT), also referred to as 4D CT. The information a respiratory monitoring system has to provide to the CT scanner depends on the specific scanner. The purpose of this study is to compare two different respiratory monitoring systems (Anzai Respiratory Gating System; C-RAD Sentinel) with respect to their applicability in combination with an Aquilion Large Bore CT scanner from Toshiba. The scanner used in our clinic does not make use of the full time dependent breathing signal, but only single trigger pulses indicating the beginning of a new breathing cycle. Hence the attached respiratory monitoring system is expected to deliver accurate online trigger pulse for each breathing cycle. The accuracy of the trigger pulses sent to the CT scanner has to be ensured by the selected respiratory monitoring system. Since a trigger pulse (output signal) of a respiratory monitoring system is a function of the measured breathing signal (input signal), the typical clinical range of the input signal is estimated for both examined respiratory monitoring systems. Both systems are analyzed based on the following parameters: time resolution, signal amplitude, noise, signal-to-noise ratio (SNR), signal linearity, trigger compatibility, and clinical examples. The Anzai system shows a better SNR (≥ 28 dB) than the Sentinel system (≥ 14.6 dB). In terms of compatibility with the cycle-based image sorting algorithm of the Toshiba CT scanner, the Anzai system benefits from the possibility to generate cycle-based triggers, whereas the Sentinel system is only able to generate amplitude-based triggers. In clinical practice, the combination of a Toshiba CT scanner and the Anzai system will provide better results due to the compatibility of the image sorting and trigger release methods.

Normal tissue complication models for clinically relevant acute esophagitis (≥ grade 2) in patients treated with dose differentiated accelerated radiotherapy (DART-bid).

BACKGROUND: One of the primary dose-limiting toxicities during thoracic irradiation is acute esophagitis (AE). The aim of this study is to investigate dosimetric and clinical predictors for AE grade ≥ 2 in patients treated with accelerated radiotherapy for locally advanced non-small cell lung cancer (NSCLC). PATIENTS AND METHODS: 66 NSCLC patients were included in the present analysis: 4 stage II, 44 stage IIIA and 18 stage IIIB. All patients received induction chemotherapy followed by dose differentiated accelerated radiotherapy (DART-bid). Depending on size (mean of three perpendicular diameters) tumors were binned in four dose groups: <2.5 cm 73.8 Gy, 2.5-4.5 cm 79.2 Gy, 4.5-6 cm 84.6 Gy, >6 cm 90 Gy. Patients were treated in 3D target splitting technique. In order to estimate the normal tissue complication probability (NTCP), two Lyman models and the cutoff-logistic regression model were fitted to the data with AE ≥ grade 2 as statistical endpoint. Inter-model comparison was performed with the corrected Akaike information criterion (AICc), which calculates the model's quality of fit (likelihood value) in relation to its complexity (i.e. number of variables in the model) corrected by the number of patients in the dataset. Toxicity was documented prospectively according to RTOG. RESULTS: The median follow up was 686 days (range 84-2921 days), 23/66 patients (35 %) experienced AE ≥ grade 2. The actuarial local control rates were 72.6 % and 59.4 % at 2 and 3 years, regional control was 91 % at both time points. The Lyman-MED model (D50 = 32.8 Gy, m = 0.48) and the cutoff dose model (Dc = 38 Gy) provide the most efficient fit to the current dataset. On multivariate analysis V38 (volume of the esophagus that receives 38 Gy or above, 95 %-CI 28.2-57.3) was the most significant predictor of AE ≥ grade 2 (HR = 1.05, CI 1.01-1.09, p = 0.007). CONCLUSION: Following high-dose accelerated radiotherapy the rate of AE ≥ grade 2 is slightly lower than reported for concomitant radio-chemotherapy with the additional benefit of markedly increased loco-regional tumor control. In the current patient cohort the most significant predictor of AE was found to be V38. A second clinically useful parameter in treatment planning may be MED (mean esophageal dose).

Stereotactic radiotherapy of intrapulmonary lesions: comparison of different dose calculation algorithms for Oncentra MasterPlan®.

BACKGROUND: The use of high accuracy dose calculation algorithms, such as Monte Carlo (MC) and Collapsed Cone (CC) determine dose in inhomogeneous tissue more accurately than pencil beam (PB) algorithms. However, prescription protocols based on clinical experience with PB are often used for treatment plans calculated with CC. This may lead to treatment plans with changes in field size (FS) and changes in dose to organs at risk (OAR), especially for small tumor volumes in lung tissue treated with SABR. METHODS: We re-evaluated 17 3D-conformal treatment plans for small intrapulmonary lesions with a prescription of 60 Gy in fractions of 7.5 Gy to the 80% isodose. All treatment plans were initially calculated in Oncentra MasterPlan® using a PB algorithm and recalculated with CC (CCre-calc). Furthermore, a CC-based plan with coverage similar to the PB plan (CCcov) and a CC plan with relaxed coverage criteria (CCclin), were created. The plans were analyzed in terms of Dmean, Dmin, Dmax and coverage for GTV, PTV and ITV. Changes in mean lung dose (MLD), V10Gy and V20Gy were evaluated for the lungs. The re-planned CC plans were compared to the original PB plans regarding changes in total monitor units (MU) and average FS. RESULTS: When PB plans were recalculated with CC, the average V60Gy of GTV, ITV and PTV decreased by 13.2%, 19.9% and 41.4%, respectively. Average Dmean decreased by 9% (GTV), 11.6% (ITV) and 14.2% (PTV). Dmin decreased by 18.5% (GTV), 21.3% (ITV) and 17.5% (PTV). Dmax declined by 7.5%. PTV coverage correlated with PTV volume (p < 0.001). MLD, V10Gy, and V20Gy were significantly reduced in the CC plans. Both, CCcov and CCclin had significantly increased MUs and FS compared to PB. CONCLUSIONS: Recalculation of PB plans for small lung lesions with CC showed a strong decline in dose and coverage in GTV, ITV and PTV, and declined dose in the lung. Thus, switching from a PB algorithm to CC, while aiming to obtain similar target coverage, can be associated with application of more MU and extension of radiotherapy fields, causing greater OAR exposition.

Comparison and limitations of DVH-based NTCP models derived from 3D-CRT and IMRT data for prediction of gastrointestinal toxicities in prostate cancer patients by using propensity score matched pair analysis.

PURPOSE: This study compared normal tissue complication probability (NTCP) modeling of chronic gastrointestinal toxicities following prostate cancer treatment for 2 treatment modalities. Possible factors causing discrepancies in optimal NTCP model parameters between 3-dimensional conformal radiation therapy (3D-CRT) and intensity modulated RT (IMRT) were analyzed and discussed, including the impact of patient characteristics, image guidance, toxicity scoring bias, and NTCP model limitations. METHODS AND MATERIALS: Rectal wall dose-volume histograms of 1115 patients treated for prostate cancer under an adaptive radiation therapy protocol were used to model gastrointestinal toxicity grade ≥2 (according to Common Terminology Criteria for Adverse Events). A total of 457 patients were treated with 3D-CRT and 658 with IMRT. 3D-CRT patients were matched to IMRT patients based on various patient characteristics, using a propensity score-based algorithm. Parameters of the Lyman equivalent uniform dose and cut-off dose logistic regression NTCP models were estimated for the 2 matched treatment modalities and the combined group. RESULTS: After they were matched, the 3D-CRT and IMRT groups contained 275 and 550 patients with a large discrepancy of 28.7% versus 7.8% toxicities, respectively (P<.001). For both NTCP models, optimal parameters found for the 3D-CRT groups did not fit the IMRT patients well and vice versa. Models developed for the combined data overestimated NTCP for the IMRT patients and underestimated NTCP for the 3D-CRT group. CONCLUSIONS: Our analysis did not reveal a single definitive cause for discrepancies of model parameters between 3D-CRT and IMRT. Patient characteristics and bias in toxicity scoring, as well as image guidance alone, are unlikely causes of the large discrepancy of toxicities. Whether the cause was inherent to the specific NTCP models used in this study needs to be verified by future investigations. Because IMRT is increasingly used clinically, it is important that appropriate NTCP model parameters are determined for this treatment modality.

Phantom based evaluation of CT to CBCT image registration for proton therapy dose recalculation.

The ability to perform dose recalculation on the anatomy of the day is important in the context of adaptive proton therapy. The objective of this study was to investigate the use of deformable image registration (DIR) and cone beam CT (CBCT) imaging to generate the daily stopping power distribution of the patient. We investigated the deformation of the planning CT scan (pCT) onto daily CBCT images to generate a virtual CT (vCT) using a deformable phantom designed for the head and neck (H & N) region. The phantom was imaged at a planning CT scanner in planning configuration, yielding a pCT and in deformed, treatment day configuration, yielding a reference CT (refCT). The treatment day configuration was additionally scanned at a CBCT scanner. A Morphons DIR algorithm was used to generate a vCT. The accuracy of the vCT was evaluated by comparison to the refCT in terms of corresponding features as identified by an adaptive scale invariant feature transform (aSIFT) algorithm. Additionally, the vCT CT numbers were compared to those of the refCT using both profiles and regions of interest and the volumes and overlap (DICE coefficients) of various phantom structures were compared. The water equivalent thickness (WET) of the vCT, refCT and pCT were also compared to evaluate proton range differences. Proton dose distributions from the same initial fluence were calculated on the refCT, vCT and pCT and compared in terms of proton range. The method was tested on a clinical dataset using a replanning CT scan acquired close in time to a CBCT scan as reference using the WET evaluation. Results from the aSIFT investigation suggest a deformation accuracy of 2-3 mm. The use of the Morphon algorithm did not distort CT number intensity in uniform regions and WET differences between vCT and refCT were of the order of 2% of the proton range. This result was confirmed by proton dose calculations. The patient results were consistent with phantom observations. In conclusion, our phantom study suggests the vCT approach is adequate for proton dose recalculation on the basis of CBCT imaging.

Hippocampal EUD in primarily irradiated glioblastoma patients.

BACKGROUND: Radiation delivery for malignant brain tumors is gradually becoming more precise. Particularly the possibilities of sparing adjacent normal structures such as the hippocampus are increasing. To determine its radiation exposure more exactly, the equivalent uniform dose (EUD) of the hippocampus was compared with further treatment parameters. This way sparing options could be found. METHODS: From the database of the University hospital of Munich 61 glioblastoma patients were selected who received primary radiotherapy in 2011. General data about the etiology, treatment course, survival of the patients and dose parameters were retrieved. RESULTS: In a linear regression analysis the side of the tumor (left hippocampus: p < 0.001/right hippocampus: p = 0.009) and its temporal location (left hippocampus: p = 0.015/right hippocampus: p = 0.033) were identified as factors with a significant influence on the EUD of the respective hippocampus. Besides this, the size of the planning target volume (PTV) and the EUD of the hippocampus correlated significantly (p = 0.027; Pearson correlation = 0.291). The median PTV size of the tumor in the right hemisphere was 386.1 ml (range 131.2-910.7 ml), and in the left hemisphere 291.3 ml (range 146.0-588.9 ml) (Kruskal-Wallis test: p = 0.048). A dose quartile analysis showed that 31 patients had a high dose exposure of the hippocampus on one side while having a moderate dose exposure in the other side. CONCLUSIONS: The radiation exposure of the respective hippocampus is dependent on the side where the tumor is located as well as on whether it is temporally located. The exposure of the contralateral hippocampus is further dependent on multiple additional factors - nevertheless a reasonable protection seems to be possible in about half of all cases.