Stewart Gaede

Stereotactic ablative radiotherapy for comprehensive treatment of oligometastatic tumors (SABR-COMET): Study protocol for a randomized phase II trial

BackgroundStereotactic ablative radiotherapy (SABR) has emerged as a new treatment option for patients with oligometastatic disease. SABR delivers precise, high-dose, hypofractionated radiotherapy, and achieves excellent rates of local control. Survival outcomes for patients with oligometastatic disease treated with SABR appear promising, but conclusions are limited by patient selection, and the lack of adequate controls in most studies. The goal of this multicenter randomized phase II trial is to assess the impact of a comprehensive oligometastatic SABR treatment program on overall survival and quality of life in patients with up to 5 metastatic cancer lesions, compared to patients who receive standard of care treatment alone.MethodsAfter stratification by the number of metastases (1-3 vs. 4-5), patients will be randomized between Arm 1: current standard of care treatment, and Arm 2: standard of care treatment + SABR to all sites of known disease. Patients will be randomized in a 1:2 ratio to Arm 1:Arm 2, respectively. For patients receiving SABR, radiotherapy dose and fractionation depends on the site of metastasis and the proximity to critical normal structures. This study aims to accrue a total of 99 patients within four years. The primary endpoint is overall survival, and secondary endpoints include quality of life, toxicity, progression-free survival, lesion control rate, and number of cycles of further chemotherapy/systemic therapy.DiscussionThis study will provide an assessment of the impact of SABR on clinical outcomes and quality of life, to determine if long-term survival can be achieved for selected patients with oligometastatic disease, and will inform the design of a possible phase III study.Trial registrationClinicaltrials.gov identifier: NCT01446744

Inter-observer and intra-observer reliability for lung cancer target volume delineation in the 4D-CT era

BACKGROUND AND PURPOSE To investigate inter-observer and intra-observer target volume delineation (TVD) error in 4D-CT imaging of thoracic tumours. MATERIALS AND METHODS Primary and nodal gross tumour volumes (GTV) of 10 lung tumours on the 10 respiratory phases of a 4D-CT scan were contoured by six radiation oncologist observers. Inter-observer and intra-observer variability were assessed by the coefficient of variation (COV) and the volume overlap index (VOI). ANOVA was performed to assess differences in inter-observer and intra-observer variability based on patient case difficulty, respiratory phase, physician seniority, and physician observer. RESULTS VOI analysis determined that inter-observer was a more significant source of error than intra-observer variability. VOI improved with the use of 4D-CT as compared to conventional CT. ANOVA analysis for COVs found case difficulty (easy versus difficult) to be significant for inter-observer primary tumour and intra-observer nodal disease delineation. Physician seniority, respiratory phase, and individual physician were not found to be significant for TVD error. CONCLUSION Variability in TVD is a major source of error in 4D-CT treatment planning. Development of measures to reduce inter-observer and intra-observer TVD variability are necessary in order to deliver high quality radiotherapy.

Detection of longitudinal lung structural and functional changes after diagnosis of radiation-induced lung injury using hyperpolarized 3He magnetic resonance imaging.

PURPOSE Therapeutic radiation doses for thoracic tumors are significantly restricted to decrease the risk of nontumor tissue damage, yet radiation-induced lung injury (RILI) still occurs in over 1/3 of thoracic radiation treatment cases. Although RILI can be clinically monitored using pulmonary function measurements, the regional functional effects of the injury are not well understood. Hyperpolarized 3He magnetic resonance imaging provides measurements of regional lung function and structure with high spatial and temporal resolution; the authors use this tool longitudinally for the first time in seven subjects after clinical diagnosis of RILI in order to better understand regional changes in lung function and structure post-RILI. METHODS All subjects underwent spirometry, plethysmography, and MRI at 3.0 T 35.1 +/- 12.2 weeks after radiation therapy commenced. Thoracic 1H, static 3He ventilation, and 3He diffusion-weighted images were acquired to generate the 3He apparent diffusion coefficient (ADC) and 3He percent ventilated volume (PVV). Four subjects returned 22.0 +/- 0.8 weeks after baseline imaging for follow-up spirometry and 3He MRI measurements of ADC and PVV. RESULTS At baseline, PVV was significantly different (p = 0.025) and lower in the ipsilateral diseased lung (55 +/- 29%) compared to the contralateral nondiseased lung (88 +/- 5%). Longitudinally, significant increases were observed for 3He MRI PVV (16% +/- 6%, p = 0.012) and 3He MRI ADC (0.02 +/- 0.01 cm2/s, p = 0.003) in the contralateral lung only, in the four subjects who returned for follow-up, while no changes in the ipsilateral lung were reported. CONCLUSIONS Hyperpolarized 3He MRI was well tolerated in all subjects with moderate to severe RILI. Functional improvements and microstructural changes were observed in the contralateral lung, while the ipsilateral lung remained stable, suggesting that functional compensatory changes may have occurred in the contralateral lung due to ipsilateral lung radiation-induced injury.

Inter and Intrafraction Uncertainty in Prostate Bed Image-Guided Radiotherapy

PURPOSE The goals of this study were to measure inter- and intrafraction setup error and prostate bed motion (PBM) in patients undergoing post-prostatectomy image-guided radiotherapy (IGRT) and to propose appropriate population-based three-dimensional clinical target volume to planning target volume (CTV-PTV) margins in both non-IGRT and IGRT scenarios. METHODS AND MATERIALS In this prospective study, 14 patients underwent adjuvant or salvage radiotherapy to the prostate bed under image guidance using linac-based kilovoltage cone-beam CT (kV-CBCT). Inter- and intrafraction uncertainty/motion was assessed by offline analysis of three consecutive daily kV-CBCT images of each patient: (1) after initial setup to skin marks, (2) after correction for positional error/immediately before radiation treatment, and (3) immediately after treatment. RESULTS The magnitude of interfraction PBM was 2.1 mm, and intrafraction PBM was 0.4 mm. The maximum inter- and intrafraction prostate bed motion was primarily in the anterior-posterior direction. Margins of at least 3-5 mm with IGRT and 4-7 mm without IGRT (aligning to skin marks) will ensure 95% of the prescribed dose to the clinical target volume in 90% of patients. CONCLUSIONS PBM is a predominant source of intrafraction error compared with setup error and has implications for appropriate PTV margins. Based on inter- and estimated intrafraction motion of the prostate bed using pre- and post-kV-CBCT images, CBCT IGRT to correct for day-to-day variances can potentially reduce CTV-PTV margins by 1-2 mm. CTV-PTV margins for prostate bed treatment in the IGRT and non-IGRT scenarios are proposed; however, in cases with more uncertainty of target delineation and image guidance accuracy, larger margins are recommended.

An evaluation of an automated 4D-CT contour propagation tool to define an internal gross tumour volume for lung cancer radiotherapy

BACKGROUND AND PURPOSE To evaluate an automated 4D-CT contouring propagation tool by its impact on the inter- and intra-physician variability in lung tumour delineation. MATERIALS AND METHODS In a previous study, six radiation oncologists contoured the gross tumour volume (GTV) and nodes on 10 phases of the 4D-CT dataset of 10 lung cancer patients to examine the intra- and inter-physician variability. In this study, a model-based deformable image registration algorithm was used to propagate the GTV and nodes on each phase of the same 4D-CT datasets. A blind review of the contours was performed by each physician and edited. Inter- and intra-physician variability for both the manual and automated methods was assessed by calculating the centroid motion of the GTV using the Pearson correlation coefficient and the variability in the internal gross tumour volume (IGTV) overlap using the Dice similarity coefficient (DSC). RESULTS The time for manual delineation was (42.7±18.6)min versus (17.7±5.4)min when the propagation tool was used. A significant improvement in the mean Pearson correlation coefficient was also observed. There was a significant decrease in mean DSC in only 1 out of 10 primary IGTVs and 2 out of 10 nodal IGTVs. Intra-physician variability was not significantly impacted (DSC>0.742). CONCLUSIONS Automated 4D-CT propagation tools can significantly decrease the IGTV delineation time without significantly decreasing the inter- and intra-physician variability.

Technology Assessment of Automated Atlas Based Segmentation in Prostate Bed Contouring

BackgroundProstate bed (PB) contouring is time consuming and associated with inter-observer variability. We evaluated an automated atlas-based segmentation (AABS) engine in its potential to reduce contouring time and inter-observer variability.MethodsAn atlas builder (AB) manually contoured the prostate bed, rectum, left femoral head (LFH), right femoral head (RFH), bladder, and penile bulb of 75 post-prostatectomy cases to create an atlas according to the recent RTOG guidelines. 5 other Radiation Oncologists (RO) and the AABS contoured 5 new cases. A STAPLE contour for each of the 5 patients was generated. All contours were anonymized and sent back to the 5 RO to be edited as clinically necessary. All contouring times were recorded. The dice similarity coefficient (DSC) was used to evaluate the unedited- and edited- AABS and inter-observer variability among the RO. Descriptive statistics, paired t-tests and a Pearson correlation were performed. ANOVA analysis using logit transformations of DSC values was calculated to assess inter-observer variability.ResultsThe mean time for manual contours and AABS was 17.5- and 14.1 minutes respectively (p = 0.003). The DSC results (mean, SD) for the comparison of the unedited-AABS versus STAPLE contours for the PB (0.48, 0.17), bladder (0.67, 0.19), LFH (0.92, 0.01), RFH (0.92, 0.01), penile bulb (0.33, 0.25) and rectum (0.59, 0.11). The DSC results (mean, SD) for the comparison of the edited-AABS versus STAPLE contours for the PB (0.67, 0.19), bladder (0.88, 0.13), LFH (0.93, 0.01), RFH (0.92, 0.01), penile bulb (0.54, 0.21) and rectum (0.78, 0.12). The DSC results (mean, SD) for the comparison of the edited-AABS versus the expert panel for the PB (0.47, 0.16), bladder (0.67, 0.18), LFH (0.83, 0.18), RFH (0.83, 0.17), penile bulb (0.31, 0.23) and rectum (0.58, 0.09). The DSC results (mean, SD) for the comparison of the STAPLE contours and the 5 RO are PB (0.78, 0.15), bladder (0.96, 0.02), left femoral head (0.87, 0.19), right femoral head (0.87, 0.19), penile bulb (0.70, 0.17) and the rectum (0.89, 0.06). The ANOVA analysis suggests inter-observer variability among at least one of the 5 RO (p value = 0.002).ConclusionThe AABS tool results in a time savings, and when used to generate auto-contours for the femoral heads, bladder and rectum had superior to good spatial overlap. However, the generated auto-contours for the prostate bed and penile bulb need improvement.

An in-depth Monte Carlo study of lateral electron disequilibrium for small fields in ultra-low density lung: implications for modern radiation therapy

Modern radiation therapy techniques such as intensity-modulated radiation therapy (IMRT) and stereotactic body radiation therapy (SBRT) use tightly conformed megavoltage x-ray fields to irradiate a tumour within lung tissue. For these conditions, lateral electron disequilibrium (LED) may occur, which systematically perturbs the dose distribution within tumour and nearby lung tissues. The goal of this work is to determine the combination of beam and lung density parameters that cause significant LED within and near the tumour. The Monte Carlo code DOSXYZnrc (National Research Council of Canada, Ottawa, ON) was used to simulate four 20 × 20 × 25 cm(3) water-lung-water slab phantoms, which contained lung tissue only, or one of three different centrally located small tumours (sizes: 1 × 1 × 1, 3 × 3 × 3, 5 × 5 × 5 cm(3)). Dose calculations were performed using combinations of six beam energies (Co-60 up to 18 MV), five field sizes (1 × 1 cm(2) up to 15 × 15 cm(2)), and 12 lung densities (0.001 g cm(-3) up to 1 g cm(-3)) for a total of 1440 simulations. We developed the relative depth-dose factor (RDDF), which can be used to characterize the extent of LED (RDDF <1.0). For RDDF <0.7 severe LED occurred, and both lung and tumour dose were drastically reduced. For example, a 6 MV (3 × 3 cm(2)) field was used to irradiate a 1 cm(3) tumour embedded in lung with ultra-low density of 0.001 g cm(-3) (RDDF = 0.2). Dose in up-stream lung and tumour centre were reduced by as much as 80% with respect to the water density calculation. These reductions were worse for smaller tumours irradiated with high energy beams, small field sizes, and low lung density. In conclusion, SBRT trials based on dose calculations in homogeneous tissue are misleading as they do not reflect the actual dosimetric effects due to LED. Future clinical trials should only use dose calculation engines that can account for electron scatter, with special attention given to patients with low lung density (i.e. emphysema). In cases where tissue inhomogeneity corrections are applied, the nature of the correction used may be inadequate in predicting the correct level of LED. In either case, the dose to the tumour is not the prescribed dose and clinical response data are uncertain. The new information from this study can be used by radiation oncologists who wish to perform advanced radiation therapy techniques while avoiding the deleterious predictable dosimetric effects of LED.

The Use of CT Density Changes at Internal Tissue Interfaces to Correlate Internal Organ Motion with an External Surrogate

The purpose of this paper is to describe a non-invasive method to monitor the motion of internal organs affected by respiration without using external markers or spirometry, to test the correlation with external markers, and to calculate any time shift between the datasets. Ten lung cancer patients were CT scanned with a GE LightSpeed Plus 4-Slice CT scanner operating in a ciné mode. We retrospectively reconstructed the raw CT data to obtain consecutive 0.5 s reconstructions at 0.1 s intervals to increase image sampling. We defined regions of interest containing tissue interfaces, including tumour/lung interfaces that move due to breathing on multiple axial slices and measured the mean CT number versus respiratory phase. Tumour motion was directly correlated with external marker motion, acquired simultaneously, using the sample coefficient of determination, r(2). Only three of the ten patients showed correlation higher than r(2) = 0.80 between tumour motion and external marker position. However, after taking into account time shifts (ranging between 0 s and 0.4 s) between the two data sets, all ten patients showed correlation better than r(2) = 0.8. This non-invasive method for monitoring the motion of internal organs is an effective tool that can assess the use of external markers for 4D-CT imaging and respiratory-gated radiotherapy on a patient-specific basis.

Experimental measurements and Monte Carlo simulations for dosimetric evaluations of intrafraction motion for gated and ungated intensity modulated arc therapy deliveries.

Respiratory gated radiation therapy allows for a smaller margin expansion for the planning target volume (PTV) to account for respiratory induced motion and is emerging as a common method to treat lung and liver tumors. We investigated the dosimetric effect of free motion and gated delivery for intensity modulated arc therapy (IMAT) with experimental measurements and Monte Carlo simulations. The impact of PTV margin and duty cycle for gated delivery is studied with Monte Carlo simulations. A motion phantom is used for this study. Two sets of contours were drawn on the mid-inspiration CT scan of this motion phantom. For each set of contours, an IMAT plan to be delivered with constant dose rate was created. The plans were generated on a CT scan of the phantom in the static condition with 3 mm PTV margin and applied to the motion phantom under four conditions: static, full superior-inferior (SI) motion (A = 1 cm, T = 4 s) and gating conditions (25% and 50% duty cycles) with full SI motion. A 6 by 15 cm piece of radiographic film was placed in the sagittal plane of the phantom and then irradiated under all measurement conditions. Film calibration was performed with a step-wedge method to convert optical density to dose. Gated IMAT delivery was first validated in 2D by comparing static film with that from gating and full motion. A previously verified simulation tool for IMRT that takes the log files from the multileaf collimator (MLC) controller and the gating system were adapted to simulate the delivered IMAT treatment for full 3D dosimetric analysis. The IMAT simulations were validated against the 2D film measurements. The resultant IMAT simulations were evaluated with dose criteria, dose-volume histograms and 3D gamma analysis. We validated gated IMAT deliveries when we compared the static film with the one from gating using 25% duty cycle using 2D gamma analysis. Within experimental and setup uncertainties, film measurements agreed with their corresponding simulated plans using 2D gamma analysis. Finally, when planning with margins designed for gating with 25% duty cycle and applying 50% or no gating during treatment, the dose differences in D(min,) D(99%) and D(95%) of the clinical target volume can be up to 27 cGy, 20 cGy and 18 cGy, respectively, for a plan with 200 cGy prescription dose. We have experimentally delivered gated IMAT with constant dose rate to a motion phantom and assessed their accuracies with film dosimetry and Monte Carlo simulations. Film dosimetry demonstrated that 25% gating and static plans are within 2%, 2 mm. The Monte Carlo simulation method was employed to generate dose delivered in 3D to a motion phantom, and the dosimetric results were reported. Since our film measurements agreed well with Monte Carlo simulations, we can reliably use this simulation tool to further study the dosimetric effects of target motion and effectiveness of gating for IMAT deliveries.

The effect of an inconsistent breathing amplitude on the relationship between an external marker and internal lung deformation in a porcine model

PURPOSE Investigate the relationship between the motion of the Varian Real-time Position Management (RPM) device and the internal motion of a pig during induced inconsistencies in the amplitude of breathing. METHODS Twelve studies were performed on four ventilated female Landrace cross pigs using a GE Healthcare, Discovery CT 750 HD scanner. In each study, a 4.0 cm section (64 slices) of the pigs lungs was repeatedly scanned 20 times using cine mode, each time lasting more than one breathing cycle. During these cine scans, a Varian RPM device was used to collect respiratory amplitudes and the ventilator air return tube was periodically crimped to induce inconsistent breathing amplitudes. Each breathing cycle and its associated cine scan were categorized as either consistent or inconsistent, based on thresholds of the minimum expiration and maximum inspiration amplitudes. From the group of consistent amplitude cine scans in a study, a reference scan was chosen. The effect of inconsistent breathing amplitudes on the relationship between the motion of the RPM marker and the motion within three regions of interest (in each lung and the chest wall) was investigated with two methods: (1) A 4D-CT sorting algorithm based on RPM amplitude was used to sort volumes into 4D-CT phase bins. Within each phase bin, the nonlinear deformation of volumes collected during consistent and inconsistent breathing amplitude was calculated with respect to the reference volume. The magnitude of the deformations (in mm) were compared to determine if inconsistent breathing amplitude caused greater deformations. (2) Nonlinear deformations between each CT volume from a cine scan and the maximum expiration volume of the reference scan were calculated. Regression analyses between the nonlinear deformations within three regions of interest (in each lung and the chest wall) and the RPM amplitudes were performed to test the effect of inconsistent breathing amplitudes on the linearity of the relationship between the 3D motion of internal anatomy and the 1D motion of the RPM external marker. RESULTS (1) Inconsistent versus consistent breathing amplitudes caused a significant increase in deformation relative to the reference scan within the left lung (1.40 +/- 0.42 versus 1.29 +/- 0.36 mm, p < 0.05). (2) One-to-one correspondences between motions of internal anatomies and motion of the RPM external marker did not exist. The regression lines between the two types of motions did not yield an identity relationship (unity slope and zero intercept). Inconsistent breathing produced significantly different regression lines than consistent breathing in ten of the 12 studies within a left lung region of interest. CONCLUSIONS The results of these two studies indicate that inconsistency in the amplitude of breathing disrupted the correspondence between the motion of the external marker and internal anatomies. As a consequence, radiation therapy of tumors embedded in lung tissue may be prone to significant errors if inconsistent breathing amplitudes occur during treatment.