Early radiologic and metabolic tumour response assessment during combined chemo-radiotherapy for locally advanced NSCLC.

Background: The role of early treatment response for patients with locally advanced non-small cell lung cancer (LA-NSCLC) treated with concurrent chemo-radiotherapy (cCRT) is unclear. The study aims to investigate the predictive value of response to induction chemotherapy (iCX) and the correlation with pattern of failure (PoF). Materials and methods: Patients with LA-NSCLC treated with cCRT were included for analyses (n = 276). Target delineations were registered from radiotherapy planning PET/CT to diagnostic PET/CT, in between which patients received iCX. Volume, sphericity, and SUVpeak were extracted from each scan. First site of failure was categorised as loco-regional (LR), distant (DM), or simultaneous LR+M (LR+M). Fine and Gray models for PoF were performed: a baseline model (including performance status (PS), stage, and histology), an image model for squamous cell carcinoma (SCC), and an image model for non-SCC. Parameters included PS, volume (VOL) of tumour, VOL of lymph nodes, ΔVOL, sphericity, SUVpeak, ΔSUVpeak, and oligometastatic disease. Results: Median follow-up was 7.6 years. SCC had higher sub-distribution hazard ratio (sHR) for LRF (sHR = 2.771 [1.577:4.87], p < 0.01) and decreased sHR for DM (sHR = 0.247 [0.125:0.485], p  <  0.01). For both image models, high diagnostic SUVpeak increased risk of LRF (sHR = 1.059 [1.05:1.106], p < 0.01 for SCC, sHR = 1.12 [1.03:1.21], p < 0.01 for non-SCC). Patients with SCC and less decrease in VOL had higher sHR for DM (sHR = 1.025[1.001:1.048] pr. % increase, p = 0.038). Conclusion: Poor response in disease volume was correlated with higher sHR of DM for SCC, no other clear correlation of response and PoF was observed. Histology significantly correlated with PoF with SCC prone to LRF and non-SCC prone to DM as first site of failure. High SUVpeak at diagnosis increased the risk of LRF for both histologies.

Independent Monte Carlo dose calculation identifies single isocenter multi-target radiosurgery targets most likely to fail pre-treatment measurement.

PURPOSE: For individual targets of single isocenter multi-target (SIMT) Stereotactic radiosurgery (SRS), we assess dose difference between the treatment planning system (TPS) and independent Monte Carlo (MC), and demonstrate persistence into the pre-treatment Quality Assurance (QA) measurement. METHODS: Treatment plans from 31 SIMT SRS patients were recalculated in a series of scenarios designed to investigate sources of discrepancy between TPS and independent MC. Targets with > 5% discrepancy in DMean[Gy] after progressing through all scenarios were measured with SRS MapCHECK. A matched pair analysis was performed comparing SRS MapCHECK results for these targets with matched targets having similar characteristics (volume & distance from isocenter) but no such MC dose discrepancy. RESULTS: Of 217 targets analyzed, individual target mean dose (DMean[Gy]) fell outside a 5% threshold for 28 and 24 targets before and after removing tissue heterogeneity effects, respectively, while only 5 exceeded the threshold after removing effect of patient geometry (via calculation on StereoPHAN geometry). Significant factors affecting agreement between the TPS and MC included target distance from isocenter (0.83% decrease in DMean[Gy] per 2 cm), volume (0.15% increase per cc), and degree of plan modulation (0.37% increase per 0.01 increase in modulation complexity score). SRS MapCHECK measurement had better agreement with MC than with TPS (2%/1 mm / 10% threshold gamma pass rate (GPR) = 99.4 ± 1.9% vs. 93.1 ± 13.9%, respectively). In the matched pair analysis, targets exceeding 5% for MC versus TPS also had larger discrepancies between TPS and measurement with no GPR (2%/1 mm / 10% threshold) exceeding 90% (71.5% ± 16.1%); whereas GPR was high for matched targets with no such MC versus TPS difference (96.5% ± 3.3%, p = 0.01). CONCLUSIONS: Independent MC complements pre-treatment QA measurement for SIMT SRS by identifying problematic individual targets prior to pre-treatment measurement, thus enabling plan modifications earlier in the planning process and guiding selection of targets for pre-treatment QA measurement.

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.

Different benefits of adaptive radiotherapy for different histologies of NSCLC.

BACKGROUND: Adenocarcinoma (AC) and squamous cell carcinoma (SCC) are the most frequent histological subtypes of non-small cell lung cancer (NSCLC). The aim of this study was to investigate how patients with AC and SCC benefit from image-guided adaptive radiotherapy (ART) with tumour match. MATERIAL AND METHODS: Consecutive patients diagnosed with AC or SCC of the lung treated with definitive chemo-radiotherapy before and after the implementation of ART and tumour match were retrospectively included for analyses. Data collection included baseline patient and treatment characteristics in addition to clinical data on radiation pneumonitis (RP), failure, and survival. Patients were divided into four categories based on their histology and treatment before (n = 173 [89 AC and 84 SCC]) and after implementation of ART (n = 240 [141 AC and 99 SCC]). RESULTS: Median follow-up was 5.7 years for AC and 6.3 years for SCC. Mean lung dose decreased for both histologies with ART, whereas mean heart dose only decreased for patients with AC. Incidences of grade 3 and 5 RP decreased for both histologies with ART. Loco-regional failure (LRF) rates decreased significantly for patients with SCC after ART (p = .04), no significant difference was observed for AC. Overall survival (OS) increased significantly for SCC after ART (p < .01): the 2-year OS increased from 31.0% (95% confidence interval [CI] [22.5-42.6]) to 54.5% (95% CI [45.6-65.3]). No significant effect on OS was observed for patients with AC. CONCLUSION: ART and tumour match in the radiotherapeutic treatment of patients with locally advanced NSCLC primarily led to decreased LRF and improved OS for patients with SCC.

Generation, validation, and benchmarking of a commercial independent Monte Carlo calculation beam model for multi-target SRS.

BACKGROUND: Dosimetric validation of single isocenter multi-target radiosurgery plans is difficult due to conditions of electronic disequilibrium and the simultaneous irradiation of multiple off-axis lesions dispersed throughout the volume. Here we report the benchmarking of a customizable Monte Carlo secondary dose calculation algorithm specific for multi-target radiosurgery which future users may use to guide their commissioning and clinical implementation. PURPOSE: To report the generation, validation, and clinical benchmarking of a volumetric Monte Carlo (MC) dose calculation beam model for single isocenter radiosurgery of intracranial multi-focal disease. METHODS: The beam model was prepared within SciMoCa (ScientificRT, Munich Germany), a commercial independent dose calculation software, with the aim of broad availability via the commercial software for use with single isocenter radiosurgery. The process included (1) definition & acquisition of measurement data required for beam modeling, (2) tuning model parameters to match measurements, (3) validation of the beam model via independent measurements and end-to-end testing, and finally, (4) clinical benchmarking and validation of beam model utility in a patient specific QA setting. We utilized a 6X Flattening-Filter-Free photon beam from a TrueBeam STX linear accelerator (Siemens Healthineers, Munich Germany). RESULTS: In addition to the measured data required for standard IMRT/VMAT (depth dose, central axis profiles & output factors, leaf gap), beam modeling and validation for single-isocenter SRS required central axis and off axis (5 cm & 9 cm) small field output factors and comparison between measurement and simulation of backscatter with aperture for jaw much greater than MLCs. Validation end-to-end measurements included SRS MapCHECK in StereoPHAN geometry (2%/1 mm Gamma = 99.2% ±â€¯2.2%), and OSL & scintillator measurements in anthropomorphic STEEV phantom (6 targets, volume = 0.1-4.1cc, distance from isocenter = 1.2-7.9 cm) for which mean difference was -1.9% ±â€¯2.2%. For 10 patient cases, MC for individual PTVs was -0.8% ±â€¯1.5%, -1.3% ±â€¯1.7%, and -0.5% ±â€¯1.8% for mean dose, D95%, and D1%, respectively. This corresponded to custom passing rates action limits per AAPM TG-218 guidelines of ±5.2%, ±6.4%, and ±6.3%, respectively. CONCLUSIONS: The beam modeling, validation, and clinical action criteria outlined here serves as a benchmark for future users of the customized beam model within SciMoCa for single isocenter radiosurgery of multi-focal disease.

A time-resolved clonogenic assay for improved cell survival and RBE measurements.

Purpose: The in vitro clonogenic assay (IVCA) is the mainstay of quantitative radiobiology. Here, we investigate the benefit of a time-resolved IVCA version (trIVCA) to improve the quantification of clonogenic survival and relative biological effectiveness (RBE) by analyzing cell colony growth behavior. Materials & Methods: In the IVCA, clonogenicity classification of cell colonies is performed based on a fixed colony size threshold after incubation. In contrast, using trIVCA, we acquire time-lapse microscopy images during incubation and track the growth of each colony using neural-net-based image segmentation. Attributes of the resulting growth curves are then used as predictors for a decision tree classifier to determine clonogenicity of each colony. The method was applied to three cell lines, each irradiated with 250 kV X-rays in the range 0-8 Gy and carbon ions of high LET (100 keV/µm, dose-averaged) in the range 0-2 Gy. We compared the cell survival curves determined by trIVCA to those from the classical IVCA across different size thresholds and incubation times. Further, we investigated the impact of the assaying method on RBE determination. Results: Size distributions of abortive and clonogenic colonies overlap consistently, rendering perfect separation via size threshold unfeasible at any readout time. This effect is dose-dependent, systematically inflating the steepness and curvature of cell survival curves. Consequently, resulting cell survival estimates show variability between 3% and 105%. This uncertainty propagates into RBE calculation with variability between 8% and 25% at 2 Gy.Determining clonogenicity based on growth curves has an accuracy of 95% on average. Conclusion: The IVCA suffers from substantial uncertainty caused by the overlap of size distributions of delayed abortive and clonogenic colonies. This impairs precise quantification of cell survival and RBE. By considering colony growth over time, our method improves assaying clonogenicity.

Leaf-individual calibration for a double stack multileaf collimator in photon radiotherapy.

Background and Purpose: In online adaptive stereotactic body radiotherapy treatments, linear accelerator delivery accuracy is essential. Recently introduced double stack multileaf collimators (MLCs) have new facets in their calibration. We established a radiation-based leaf-individual calibration (LIMCA) method for double stack MLCs. Materials and Methods: MLC leaf positions were evaluated from four cardinal angles with test patterns at measurement positions throughout the radiation field on EBT3 radiochromic film for each single stack. The accuracy of the method and repeatability of the results were assessed. The effect of MLC positioning errors was characterized for a measured output factor curve and a clinical patient plan. Results: All positions in the motor step - position calibration file were optimized in the established LIMCA method. The resulting double stack mean accuracy for all angles was 0.2 ± 0.1 mm for X1 (left bank) and 0.2 ± 0.2 mm for X2 (right bank). The accuracy of the leaf position evaluation was 0.2 mm (95% confidence level). The MLC calibration remained stable over four months. Small MLC leaf position errors (e.g. 1.2 mm field size reduction) resulted in important dose errors (-5.8 %) for small quadratic fields of 0.83 × 0.83 cm2. Single stack position accuracy was essential for highly modulated treatment plans. Conclusions: LIMCA is a new double stack MLC calibration method that increases treatment accuracy from four angles and for all moving leaves.

A body mass index-based method for "MR-only" abdominal MR-guided adaptive radiotherapy.

PURPOSE: Dose calculation for MR-guided radiotherapy (MRgRT) at the 0.35 T MR-Linac is currently based on deformation of planning CTs (defCT) acquired for each patient. We present a simple and robust bulk density overwrite synthetic CT (sCT) method for abdominal treatments in order to streamline clinical workflows. METHOD: Fifty-six abdominal patient treatment plans were retrospectively evaluated. All patients had been treated at the MR-Linac using MR datasets for treatment planning and plan adaption and defCT for dose calculation. Bulk density CTs (4M-sCT) were generated from MR images with four material compartments (bone, lung, air, soft tissue). The relative electron densities (RED) for bone and lung were extracted from contoured CT structure average REDs. For soft tissue, a correlation between BMI and RED was evaluated. Dose was recalculated on 4M-sCT and compared to dose distributions on defCTs assessing dose differences in the PTV and organs at risk (OAR). RESULTS: Mean RED of bone was 1.17 ±â€¯0.02, mean RED of lung 0.17 ±â€¯0.05. The correlation between BMI and RED for soft tissue was statistically significant (p < 0.01). PTV dose differences between 4M-sCT and defCT were Dmean: -0.4 ±â€¯1.0%, D1%: -0.3 ±â€¯1.1% and D95%: -0.5 ±â€¯1.0%. OARs showed D2%: -0.3 ±â€¯1.9% and Dmean: -0.1 ±â€¯1.4% differences. Local 3D gamma index pass rates (2%/2mm) between dose calculated using 4M-sCT and defCT were 96.8 ±â€¯2.6% (range 89.9-99.6%). CONCLUSION: The presented method for sCT generation enables precise dose calculation for MR-only abdominal MRgRT.

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.

Radiation-induced contrast enhancement following proton radiotherapy for low-grade glioma depends on tumor characteristics and is rarer in children than adults.

BACKGROUND AND PURPOSE: Proton beam radiotherapy (PRT) is used in the treatment of low-grade glioma (LGG) to mitigate long-term sequelae. Following PRT, increased rates of radiation-induced contrast enhancements (RICE) are suspected but poorly understood. MATERIALS AND METHODS: We analyzed consecutive 227 patients (42 children and 185 adults) treated with PRT (54 Gy RBE) for LGG from 2010 to 2020 and followed with serial clinical exams and magnetic resonance imaging for in median 5.6 years. RESULTS: Tumors were graded WHO 1 in a minority (n = 22, 12%) of adults, but a majority of children (n = 29, 69%). In contrast, tumors were graded WHO 2 in the majority (n = 160, 87%) of adults and a minority of children (n = 10, 24%). Five-year overall survival following PRT was 81% in adults and 91% in children. The risk of RICE was 5-fold more frequent in adults (25%) vs. children (5%; p = 0.0043). In children and adults, RICE were symptomatic in 50% and 55% (n = 1 and 26) of cases with CTCAE grade 0 in 47% (n = 23), grade 1 in 25% (n = 12), 0% grade 2 (n = 0) and 29% grade 3 (n = 14), respectively. In adults, RICE risk was associated to WHO grading (8% in WHO grade 1 vs. 24% in WHO grade 2, p = 0.026), independent of age (p = 0.44) and irradiation dose (p = 0.005), but not independent of IDH mutational status. CONCLUSIONS: These data demonstrate effectiveness of PRT for LGG in both children and adults. The RICE risk is lower in children which are a main target group for PRT and differs with WHO grading.

Dose escalation to hypoxic subvolumes in head and neck cancer: A randomized phase II study using dynamic [18F]FMISO PET/CT.

BACKGROUND AND PURPOSE: Tumor hypoxia is a major cause of resistance to radiochemotherapy in locally advanced head-and-neck cancer (LASCCHN). We present results of a randomized phase II trial on hypoxia dose escalation (DE) in LASCCHN based on dynamic [18F]FMISO (dynFMISO) positron emission tomography (PET). The purpose was to confirm the prognostic value of hypoxia PET and assess feasibility, toxicity and efficacy of hypoxia-DE. MATERIALS AND METHODS: Patients with LASCCHN underwent baseline dynFMISO PET/CT. Hypoxic volumes (HV) were derived from dynFMISO data. Patients with hypoxic tumors (HV > 0) were randomized into standard radiotherapy (ST: 70Gy/35fx) or dose escalation (DE: 77Gy/35fx) to the HV. Patients with non-hypoxic tumors were treated with ST. After a minimum follow-up of 2 years feasibility, acute/late toxicity and local control (LC) were analyzed. RESULTS: The study was closed prematurely due to slow accrual. Between 2009 and 2017, 53 patients were enrolled, 39 (74%) had hypoxic tumors and were randomized into ST or DE. For non-hypoxic patients, 100% 5-year LC was observed compared to 74% in patients with hypoxic tumors (p = 0.039). The difference in 5-year LC between DE (16/19) and ST (10/17) was 25%, p = 0.150. No relevant differences related to acute and late toxicities between the groups were observed. CONCLUSION: This study confirmed the prognostic value of hypoxia PET in LASCCHN for LC. Outcome after hypoxia DE appears promising and may support the concept of DE. Slow accrual and premature closure may partly be due to a high complexity of the study setup which needs to be considered for future multicenter trials.

Survival benefits for non-small cell lung cancer patients treated with adaptive radiotherapy.

INTRODUCTION: Tumor match and adaptive radiotherapy based on on-treatment imaging increases the precision of RT. This allows a reduction of treatment volume and, consequently, of the dose to organs at risk. We investigate the clinical benefits of tumor match and adaptive radiotherapy for a cohort of non-small cell lung cancer patients (NSCLC). METHODS: In 2013, tumor match and adaptive radiotherapy based on daily cone-beam CT scans was introduced to ensure adaption of the radiotherapy treatment plan for all patients with significant anatomical changes during radiotherapy. Before 2013, the daily cone-beam CT scans were matched on the vertebra and anatomical changes were not evaluated systematically. To estimate the effect of tumor match and adaptive radiotherapy, 439 consecutive NSCLC patients treated with definitive chemo-radiotherapy (50-66 Gy/25-33 fractions, 2010-2018) were investigated retrospectively. They were split in two groups, pre-ART (before tumor match and adaptive radiotherapy, 184 patients), and ART (after tumor match and adaptive radiotherapy, 255 patients) and compared with respect to clinical, treatment-specific and dosimetric variables (χ2 tests, Mann Whitney U tests), progression, survival and radiation pneumonits (CTCAEv3). Progression-free and overall survival as well as radiation pneumonitis were compared with log-rank tests. Hazard ratios were estimated from Cox proportional hazard regression. RESULTS: No significant differences in stage (p = 0.36), histology (p = 0.35), PS (p = 0.12) and GTV volumes (p = 0.24) were observed. Concomitant chemotherapy was administered more frequently in the ART group (78%) compared to preART (64%), p < 0.001. Median[range] PTV volumes decreased from 456 [71;1262] cm3 (preART) to 270 [31;1166] cm3 (ART), p < 0.001, thereby significantly reducing mean doses to lungs (median, preART 16.4 [1.9;24.7] Gy, ART 12.1 [1.7;19.4] Gy, p < 0.001) and heart (median, preART 8.0 [0.1;32.1] Gy, ART 4.4 [0.1;33.9] Gy, p < 0.001). The incidence of RP at nine months decreased significantly with ART (50% to 20% for symptomatic RP (≥G2), 21% to 7% for severe RP (≥G3), 6% to 0.4% for lethal RP (G5), all p < 0.001). The two-year progression free survival increased from 22% (preART) to 30% (ART), while the overall survival increased from 43% (preART) to 56% (ART). The median overall survival time increased from 20 (preART) to 28 months (ART). CONCLUSION: Tumor match and adaptive radiotherapy significantly decreased radiation pneumonitis, while maintaining loco-regional control. Further, we observed a significantly improved progression-free and overall survival.

Quality assurance for on-table adaptive magnetic resonance guided radiation therapy: A software tool to complement secondary dose calculation and failure modes discovered in clinical routine.

Online adaption of treatment plans on a magnetic resonance (MR)-Linac enables the daily creation of new (adapted) treatment plans using current anatomical information of the patient as seen on MR images. Plan quality assurance (QA) relies on a secondary dose calculation (SDC) that is required because a pretreatment measurement is impossible during the adaptive workflow. However, failure mode and effect analysis of the adaptive planning process shows a large number of error sources, and not all of them are covered by SDC. As the complex multidisciplinary adaption process takes place under time pressure, additional software solutions for pretreatment per-fraction QA need to be used. It is essential to double-check SDC input to ensure a safe treatment delivery. Here, we present an automated treatment plan check tool for adaptive radiotherapy (APART) at a 0.35 T MR-Linac. It is designed to complement the manufacturer-provided adaptive QA tool comprising SDC. Checks performed by APART include contour analysis, electron density map examinations, and fluence modulation complexity controls. For nine of 362 adapted fractions (2.5%), irregularities regarding missing slices in target volumes and organs at risks as well as in margin expansion of target volumes have been found. This demonstrates that mistakes occur and can be detected by additional QA measures, especially contour analysis. Therefore, it is recommended to implement further QA tools additional to what the manufacturer provides to facilitate an informed decision about the quality of the treatment plan.

Radiation induced contrast enhancement after proton beam therapy in patients with low grade glioma - How safe are protons?

PURPOSE: The optimal treatment strategy for low-grade glioma (LGG) is still a matter of controversy. Considering that the prognosis is typically favorable, the prevention of late sequelae is of particular importance. Proton beam therapy (PRT) has the potential to further reduce the burden of treatment related side effects. We set out to evaluate the clinical outcome of proton irradiation with a particular focus on morphologic features on magnetic resonance imaging (MRI). METHODS: We assessed prospectively 110 patients who received radiotherapy with protons for histologically proven LGG. Clinical and radiological information were analyzed resulting in more than 1200 available MRI examinations with a median follow-up of 39 months. Newly diagnosed contrast-enhancing lesions on MRI were delineated and correlated with parameters of the corresponding treatment plan. A voxel-based dose-matched paired analysis of the linear energy transfer (LET) inside vs outside lesions was performed. RESULTS: Proton beam irradiation of patients with low-grade glioma results in overall survival (OS) of 90% after seven years. Median progression free survival had not yet been reached with surviving fraction of 54% after seven years. The incidence of temporary or clinically silent radiation induced contrast enhancement was significantly higher than previously assumed, however, symptomatic radiation necrosis was only detected in one patient. These radiation-induced contrast-enhancing lesions were almost exclusively seen at the distal beam end of the proton beam. In 22 out of 23 patients, the average LET of voxels inside contrast-enhancing lesions was significantly increased, compared to dose-matched voxels outside the lesions. CONCLUSION: Symptomatic radiation necrosis following PRT was as rare as conventional photon-based treatment series suggest. However, the increased incidence of asymptomatic radiation-induced brain injuries with an increased average LET observed in this cohort provides strong clinical evidence to support the hypothesis that the relative biological effectiveness of protons is variable and different to the fixed factor of 1.1 currently used worldwide.

Biologically consistent dose accumulation using daily patient imaging.

BACKGROUND: This work addresses a basic inconsistency in the way dose is accumulated in radiotherapy when predicting the biological effect based on the linear quadratic model (LQM). To overcome this inconsistency, we introduce and evaluate the concept of the total biological dose, bEQDd. METHODS: Daily computed tomography imaging of nine patients treated for prostate carcinoma with intensity-modulated radiotherapy was used to compute the delivered deformed dose on the basis of deformable image registration (DIR). We compared conventional dose accumulation (DA) with the newly introduced bEQDd, a new method of accumulating biological dose that considers each fraction dose and tissue radiobiology. We investigated the impact of the applied fractionation scheme (conventional/hypofractionated), uncertainties induced by the DIR and by the assigned α/ß-value. RESULTS: bEQDd was systematically higher than the conventionally accumulated dose with difference hot spots of 3.3-4.9 Gy detected in six out of nine patients in regions of high dose gradient in the bladder and rectum. For hypofractionation, differences are up to 8.4 Gy. The difference amplitude was found to be in a similar range to worst-case uncertainties induced by DIR and was higher than that induced by α/ß. CONCLUSION: Using bEQDd for dose accumulation overcomes a potential systematic inaccuracy in biological effect prediction based on accumulated dose. Highest impact is found for serial-type late responding organs at risk in dose gradient regions and for hypofractionation. Although hot spot differences are in the order of several Gray, in dose-volume parameters there is little difference compared with using conventional or biological DA. However, when local dose information is used, e.g. dose surface maps, difference hot spots can potentially change outcomes of dose-response modelling and adaptive treatment strategies.

Prospectively scored pulmonary toxicities in non-small cell lung cancer: Results from a randomized phase II dose escalation trial.

PURPOSE: Prospectively scored radiation pneumonitis (RP) observed in a national, randomized phase II dose-escalation trial for patients with locally advanced non-small cell lung cancer (NSCLC) was investigated. METHODS: Patients with stage IIB-IIIB histologically proven NSCLC were treated with concomitant chemo-radiotherapy (oral Vinorelbine 3times/week) at 60 Gy/30fx (A-59pts) and 66 Gy/33fx (B-58pts) from 2009 to 2013 at five Danish RT centers. Grade 2 RP (CTCAEv3.0) was investigated with univariate analysis for association with clinical and dosimetric parameters, including dyspnea and cough at baseline and during RT. Multivariable logistic regression and Cox regression with regularization were used to find a multivariable model for RP ≥ G2. RESULTS: Despite a tendency of higher mean lung dose in the high-dose arm (median[range] A = 14.9 Gy[5.8,23.1], B = 17.5 Gy[8.6,24.8], p = 0.075), pulmonary toxicities were not significantly different (RP ≥ G2 41%(A) and 52%(B), p = 0.231). A Kaplan Meier analysis of the time to RP ≥ G2 between the two arms did not reach statistical significance (p = 0.180). Statistically significant risk factors for RP ≥ G2 were GTV size (OR = 2.091/100 cm3, p = 0.002), infection at baseline or during RT (OR = 8.087, p = 0.026), dyspnea at baseline (OR = 2.184, p = 0.044) and increase of cough during RT (OR = 2.787, p = 0.008). In the multivariable logistic regression and the Cox regression analysis, the deviances of the most predictive models were within one standard deviation of the null model. CONCLUSION: No statistical difference between the high- and low dose arm was found in the risk of developing RP. The univariate analysis identified target volume, infection, dyspnea at baseline, and increase of cough during RT as risk factors for RP. The number of patients was too small to establish a statistically sound multivariable model.

How can scanned proton beam treatment planning for low-grade glioma cope with increased distal RBE and locally increased radiosensitivity for late MR-detected brain lesions?

A novel risk model has recently been proposed for the occurrence of late contrast-enhancing brain lesions (CEBLs) after proton irradiation of low-grade glioma (LGG) patients. It predicts a strong dependence on dose-weighted linear-energy transfer (LETd effect) and an increased radiosensitivity of the ventricular proximity, a 4-mm fringe surrounding the ventricular system (VP4mm effect). On this basis, we investigated (A) how these two risk factors and patient-specific anatomical and treatment plan (TP) features contribute to normal tissue complication probability (NTCP) and (B) if conventional LETd -reduction techniques like multiple-field TP are able to reduce NTCP. (A) The LGG model cohort (N = 110) was stratified with respect to prescribed dose, tumor grade, and treatment field configuration. NTCP predictions and CEBL occurrence rates per strata were analyzed. (B) The effect of multiple-field TP was investigated in two patient groups: (i) nine high-risk subjects with extended lateral target volumes who had developed CEBLs after single-beam treatments were retrospectively replanned with a clinical standard two-field setting using almost orthogonal fields and strictly opposing fields, (ii) single-field treatments were simulated for seven low-risk patients with small central target volumes clinically treated with two strictly opposing fields. (A) In the model cohort, we identified the exposure of the radiosensitive VP4mm fringe with proton field components of increased biological effectiveness as dominant NTCP driving factor. We observed that larger target volumes and location lateral to the main ventricles, both being characteristic for WHO°II tumors, presented with the highest complication risks. Among subjects of an equal dose prescription of 54 Gy(RBE), the highest median NTCP was obtained for the WHO°II group treated with two fields using sharp angles. (B) Regarding the effect of multiple-field plans, we found that an NTCP reduction was only achievable in the low-risk group where the LETd effect dominates and the VP4mm effect is small. NTCP of the single-field plans was 23% higher compared to the clinical opposing field plan. In the high-risk group, where the VP4mm effect dominates the risk, both two-field scenarios yielded 44% higher NTCP predictions compared to the clinical single-field plans. The interplay of an increased radiosensitivity in the VP4mm fringe with proton field components of increased biological effectiveness creates a geometric complexity that can hardly be managed by current clinical TP. Our results underline that advanced biologically guided TP approaches become crucial for an effective risk minimization in proton therapy of LGG.

A framework for automated time-resolved analysis of cell colony growth after irradiation.

Understanding dose-dependent survival of irradiated cells is a pivotal goal in radiotherapy and radiobiology. To this end, the clonogenic assay is the standard in vitro method, classifying colonies into either clonogenic or non-clonogenic based on a size threshold at a fixed time. Here we developed a methodological framework for the automated analysis of time course live-cell image data to examine in detail the growth dynamics of large numbers of colonies that occur during such an experiment. We developed a segmentation procedure that exploits the characteristic composition of phase-contrast images to identify individual colonies. Colony tracking allowed us to characterize colony growth dynamics as a function of dose by extracting essential information: (a) colony size distributions across time; (b) fractions of differential growth behavior; and (c) distributions of colony growth rates across all tested doses. We analyzed three data sets from two cell lines (H3122 and RENCA) and made consistent observations in line with already published results: (i) colony growth rates are normally distributed with a large variance; (ii) with increasing dose, the fraction of exponentially growing colonies decreases, whereas the fraction of delayed abortive colonies increases; as a novel finding, we observed that (iii) mean exponential growth rates decrease linearly with increasing dose across the tested range (0-10 Gy). The presented method is a powerful tool to examine live colony growth on a large scale and will help to deepen our understanding of the dynamic, stochastic processes underlying the radiation response in vitro.

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.

Volume-dependent dose-response of the intestinal stem cell niche and lymphoid tissue.

BACKGROUND AND PURPOSE: Plasticity of the intestinal stem cell compartment in response to radiation injury is regulated by a stem cell niche. We present here the first experimental observations of a dose-volume effect of the intestinal stem cell niche and of the solitary intestinal lymphoid tissues (SILT). MATERIALS AND METHODS: Regeneration of intestinal crypts in mice was studied following irradiation of millimetre-size jejunal sections with single doses of 6 to 24 Gy and compared to total body irradiation (TBI). The statistical distribution of cells per crypt was scored and regressed to a biomathematical model. The number of SILTs was scored for different doses and field sizes and crypt regeneration was correlated with SILT proximity. RESULTS: We observed a differential dose-response of the intestinal stem cell niche at the centres of the irradiated sections, but only for field sizes below 10 mm. Irradiation of 5 mm jejunum results in an increase in crypt survival by up to an order of magnitude, compared to TBI. Distributions of cell-per-crypt numbers and comparison to biomathematical modelling suggest that these observations stem from a field size-dependent regeneration rate. The density of SILTs also exhibits a volume-dependent dose-response and increased crypt survival correlates with a proximity to SILTs. CONCLUSION: Our findings present the first observation of a field-size dependent dose-response of the intestinal stem cell niche. Its regeneration process does apparently not rely on distant radiation-sensitive resources of the organism, such as the bone marrow. Yet, our observations suggest that the niche interacts with intact tissue in millimetres distance, leading to faster crypt regeneration. The field-size dependent dose-response of SILTs posits a role of the immune system on the dose-volume effect.