Dosimetric effects of anatomical changes during fractionated photon radiation therapy in pancreatic cancer patients

Comparing different boost concepts and beam configurations for proton therapy of pancreatic cancer

Background and Purpose

Interfractional geometrical and anatomical variations impact the accuracy of proton therapy for pancreatic cancer. This study investigated field-in-field (FIF) and simultaneous integrated boost (SIB) concepts for scanned proton therapy treatment with different beam configurations.

Materials and Methods

Robustly optimized treatment plans for fifteen patients were generated using FIF and SIB techniques with two, three, and four beams. The prescribed dose in 20 fractions was 60 Gy(RBE) for the internal gross tumor volume (IGTV) and 46 Gy(RBE) for the internal clinical target volume. Verification computed tomography (vCT) scans was performed on treatment days 1, 7, and 16. Initial treatment plans were recalculated on the rigidly registered vCTs. V100% and D95% for targets and D2cm3 for the stomach and duodenum were evaluated. Robustness evaluations (range uncertainty of 3.5 %) were performed to evaluate the stomach and duodenum dose-volume parameters.

Results

For all techniques, IGTV V100% and D95% decreased significantly when recalculating the dose on vCTs (p < 0.001). The median IGTV V100% and D95% over all vCTs ranged from 74.2 % to 90.2 % and 58.8 Gy(RBE) to 59.4 Gy(RBE), respectively. The FIF with two and three beams, and SIB with two beams maintained the highest IGTV V100% and D95%. In robustness evaluations, the ΔD2cm3 of stomach was highest in two beams plans, while the ΔD2cm3 of duodenum was highest in four beams plans, for both concepts.

Conclusion

Target coverage decreased when recalculating on CTs at different time for both concepts. The FIF with three beams maintained the highest IGTV coverage while sparing normal organs the most.

Deformable image registration for composite planned doses during adaptive radiation therapy

INTRODUCTION

Some patients have significant anatomic changes during radiotherapy, necessitating an adaptive repeat CT-simulation and re-planning. This yields two unique planning datasets that introduce uncertainty into total dose records. This study explored the impact of using deformable image registration (DIR) to spatially align repeat CT-simulation images and calculate total planned dose distributions.

MATERIALS & METHODS

Data from 5 head-and-neck, 5 lung, and 5 sarcoma patients who had unanticipated re-planning during radiotherapy were analyzed in a treatment planning system (RayStation v6.1 RaySearch Laboratories). Total planned doses to normal tissues were calculated using two methods and the previously generated manual contours defined on each CT. The first method, termed 'parameter addition', simply sums the relevant DVH metrics from the initial and re-planned distributions without spatially registering the CTs. The second, termed 'dose accumulation', uses a validated hybrid contour/intensity-based DIR algorithm to deform initial CT and dose distribution onto the repeat CT and re-planning dose distribution. DVH metrics from the summed distribution on the repeat CT are then calculated. Dose differences for organs-at-risk between parameter addition and dose accumulation ≥100 cGy were assumed to be clinically relevant. To elucidate whether relevant differences were due to registration accuracy or contouring variability between CTs, the analysis was repeated using contours on the first CT and the same contours deformed to the repeat CT with DIR.

RESULTS

For all patients, high overall DIR accuracy was verified visually (qualitatively) and numerically (quantitatively) using image similarity and contour-based metrics. All head-and-neck and lung patients, and one sarcoma patient (11 of 15 total) had dose differences between parameter addition and dose accumulation ≥100 cGy, with absolute mean differences of 160 cGy (range 101-436 cGy) seen in 41 of 205 total DVH criteria. In 22 of these 41 criteria, these differences were attributed to contouring variability between CTs. After correcting for contouring variations using DIR, the mean absolute differences in 7 of these 22 criteria with a relevant result (across 6 patients) was 146 cGy (range 100-502 cGy). In only 4 DVH criteria, the DIR mapped contours had higher variations than the original contours. One lung patient had a DVH criteria exceeding the clinical dose constraint by 125 cGy with parameter addition, and with accurate DIR and dose accumulation, the criteria was actually 97 cGy lower than the constraint.

CONCLUSIONS

The use of DIR to generate total planned dose records revealed substantial dose differences in most cases compared to commonly used clinical methods (i.e. parameter addition), and altered the planned acceptance criteria in a minority. DIR is recommended to be used for future adaptive re-plans to generate total planned dose records and facilitate accurate re-contouring. More accurate dose records may also improve our understanding of clinical outcomes.

Assessment of interfraction dose variation in pancreas SBRT using daily simulation MR images

Pancreatic Cancer is associated with poor treatment outcomes compared to other cancers. High local control rates have been achieved by using hypofractionated stereotactic body radiotherapy (SBRT) to treat pancreatic cancer. Challenges in delivering SBRT include close proximity of several organs at risk (OARs) and target volume inter and intra fraction positional variations. Magnetic resonance image (MRI) guided radiotherapy has shown potential for online adaptive radiotherapy for pancreatic cancer, with superior soft tissue contrast compared to CT. The aim of this study was to investigate the variability of target and OAR volumes for different treatment approaches for pancreatic cancer, and to assess the suitability of utilizing a treatment-day MRI for treatment planning purposes. Ten healthy volunteers were scanned on a Siemens Skyra 3 T MRI scanner over two sessions (approximately 3 h apart), per day over 5 days to simulate an SBRT daily simulation scan for treatment planning. A pretreatment scan was also done to simulate patient setup and treatment. A 4D MRI scan was taken at each session for internal target volume (ITV) generation and assessment. For each volunteer a treatment plan was generated in the Raystation treatment planning system (TPS) following departmental protocols on the day one, first session dataset (D1S1), with bulk density overrides applied to enable dose calculation. This treatment plan was propagated through other imaging sessions, and the dose calculated. An additional treatment plan was generated on each first session of each day (S1) to simulate a daily replan process, with this plan propagated to the second session of the day. These accumulated mock treatment doses were assessed against the original treatment plan through DVH comparison of the PTV and OAR volumes. The generated ITV showed large variations when compared to both the first session ITV and daily ITV, with an average magnitude of 22.44% ± 13.28% and 25.83% ± 37.48% respectively. The PTV D95 was reduced by approximately 23.3% for both plan comparisons considered. Surrounding OARs had large variations in dose, with the small bowel V30 increasing by 128.87% when compared to the D1S1 plan, and 43.11% when compared to each daily S1 plan. Daily online adaptive radiotherapy is required for accurate dose delivery for pancreas cancer in the absence of additional motion management and tumour tracking techniques.

Quantification of tissue property and perfusion uncertainties in hyperthermia treatment planning: Multianalysis using polynomial chaos expansion

INTRODUCTION

Hyperthermia treatment planning (HTP) tools can guide treatment delivery, particularly with locoregional radiative phased array systems. Uncertainties in tissue and perfusion property values presently lead to quantitative inaccuracy of HTP, leading to sub-optimal treatment. Assessment of these uncertainties would allow for better judgement of the reliability of treatment plans and improve their value for treatment guidance. However, systematically investigating the impact of all uncertainties on treatment plans is a complex, high-dimensional problem and too computationally expensive for traditional Monte Carlo approaches. This study aims to systematically quantify the treatment-plan impact of tissue property uncertainties by investigating their individual contribution to, and combined impact on predicted temperature distributions.

METHODS

A novel Polynomial Chaos Expansion (PCE)-based HTP uncertainty quantification was developed and applied for locoregional hyperthermia of modelled tumours in the pancreatic head, prostate, rectum, and cervix. Patient models were based on the Duke and Ella digital human models. Using Plan2Heat, treatment plans were created to optimise tumour temperature (represented by T90) for treatment using the Alba4D system. For all 25-34 modelled tissues, the impact of tissue property uncertainties was analysed individually i.e., electrical and thermal conductivity, permittivity, density, specific heat capacity and perfusion. Next, combined analyses were performed on the top 30 uncertainties with the largest impact.

RESULTS

Uncertainties in thermal conductivity and heat capacity were found to have negligible impact on the predicted temperature ( < 1 × 10-10 °C), density and permittivity uncertainties had a small impact (< 0.3 °C). Uncertainties in electrical conductivity and perfusion can lead to large variations in predicted temperature. However, variations in muscle properties result in the largest impact at locations that could limit treatment quality, with a standard deviation up to almost 6 °C (pancreas) and 3.5 °C (prostate) for perfusion and electrical conductivity, respectively. The combined influence of all significant uncertainties leads to large variations with a standard deviation up to 9.0, 3.6, 3.7 and 4.1 °C for the pancreatic, prostate, rectal and cervical cases, respectively.

CONCLUSION

Uncertainties in tissue and perfusion property values can have a large impact on predicted temperatures from hyperthermia treatment planning. PCE-based analysis helps to identify all major uncertainties, their impact and judge the reliability of treatment plans.

The role of multiple anatomical scenarios in plan optimization for carbon ion radiotherapy of pancreatic cancer

PURPOSE /OBJECTIVE

To quantify benefits of robust optimization on multiple 4DCT acquisitions combined with off-line treatment adaptation for neoadjuvant carbon ion therapy (CIRT) of pancreatic cancer.

MATERIAL/METHODS

For 10 previously treated patients, 4DCTs were acquired around -15 (CT_Plan), -5 (RE₁), -1 (RE₂) and +6 (RE₃) days from RT start. Treatment plans were newly optimized to a dose prescription of 38.4 Gy(RBE) (8 fractions) with a constraint of 38 Gy(RBE) to 1% of the gastrointestinal organs at risk volume (D_1%). Three strategies were tested: (A) robust optimization on CT_Plan maximum exhale (0Ex) with 3 mm set-up, 3% range uncertainty, including 30%-inhale; (B) addition of the RE₁-0Ex scenario; (C) plan recalculation at each RE_i and adaptation (RP_i) according to deviation thresholds from clinical goals. The cumulative variation of target coverage and GI-OARs doses was evaluated. Duodenum contours of all 4DCTs of each patient were registered on CT_Plan-0Ex. The capacity of pre-RT acquisitions to predict duodenum position was investigated by computing the intersection of contours at CT_plan, RE₁, or their union, with respect to subsequent 4DCTs and the CTV, coupled with increasing margin.

RESULTS

(A) No recalculation exceeded the D_1% constraint. (B) The inclusion of RE₁-0Ex in the optimization problem improved inter-fraction robustness on a patient-specific basis, but was non-significant on average. (C) Half of the plans would be re-optimized to recover target coverage and/or minimize duodenum dose, at least once. A significant difference was observed between pre-RT duodenum contours when intersecting subsequent contours, either with a margin expansion.

CONCLUSION

Anatomical variations highlighted at multiple RE_i proved that a fast and efficient online adaptation is essential to optimize treatment quality of CIRT for pancreatic cancer.

Dosimetric effect of variable rectum and sigmoid colon filling during carbon ion radiotherapy to sacral chordoma

PURPOSE

Carbon ion radiotherapy (CIRT) is sensitive to anatomical density variations. We examined the dosimetric effect of variable intestinal filling condition during CIRT to ten sacral chordoma patients.

METHODS

For each patient, eight virtual computed tomography scans (vCTs) were generated by varying the density distribution within the rectum and the sigmoid in the planning computed tomography (pCT) with a density override approach mimicking a heterogeneous combination of gas and feces. Totally full and empty intestinal preparations were modelled. In addition, five different intestinal filling conditions were modelled by a mixed density pattern derived from two combined and weighted Gaussian distributions simulating gas and feces respectively. Finally, a patient-specific mixing proportion was estimated by evaluating the daily amount of gas detected in the cone beam computed tomography (CBCT). Dose distribution was recalculated on each vCT and dose volume histograms (DVHs) were examined.

RESULTS

No target coverage degradation was observed at different vCTs. Rectum and sigma dose degradation ranged respectively between: [-6.7; 21.6]GyE and [-0.7; 15.4]GyE for D_50%; [-377.4; 1197.9] and [-95.2; 1027.5] for AUC; [-1.2; 10.7]GyE and [-2.6; 21.5]GyE for D_1%.

CONCLUSIONS

Variation of intestinal density can greatly influence the penetration depth of charged particle and might compromise dose distribution. In particular cases, with large clinical target volume in very close proximity to rectum and sigmoid colon, it is appropriate to evaluate the amount of gas present in the daily CBCT images even if it is totally included in the reference planning structures.

Effect of gastrointestinal gas on the temperature distribution in pancreatic cancer hyperthermia treatment planning

PURPOSE

In pancreatic cancer treatment, hyperthermia can be added to increase efficacy of chemo- and/or radiotherapy. Gas in stomach, intestines and colon is often in close proximity to the target volume. We investigated the impact of variations in gastrointestinal gas (GG) on temperature distributions during simulated hyperthermia treatment (HT).

METHODS

We used sets of one CT and eight cone-beam CT (CBCT) scans obtained prior to/during fractionated image-guided radiotherapy in four pancreatic cancer patients. In Plan2Heat, we simulated locoregional heating by an ALBA-4D phased array radiofrequency system and calculated temperature distributions for (i) the segmented CT (sCT), (ii) sCT with GG replaced by muscle (sCT₀), (iii) sCT₀ with eight different GG distributions as visible on CBCT inserted (sCT_CBCT). We calculated cumulative temperature-volume histograms for the clinical target volume (CTV) for all ten temperature distributions for each patient and investigated the relationship between GG volume and change in ΔT₅₀ (temperature increase at 50% of CTV volume). We determined location and volume of normal tissue receiving a high thermal dose.

RESULTS

GG volume on CBCT varied greatly (9-991 cm³). ΔT₅₀ increased for increasing GG volume; maximum ΔT₅₀ difference per patient was 0.4-0.6 °C. The risk for GG-associated treatment-limiting hot spots appeared low. Normal tissue high-temperature regions mostly occurred anteriorly; their volume and maximum temperature showed moderate positive correlations with GG volume, while fat-muscle interfaces were associated with higher risks for hot spots.

CONCLUSIONS

Considerable changes in volume and position of gastrointestinal gas can occur and are associated with clinically relevant tumor temperature differences.

Quantifying and Assessing the Dosimetric Impact of Changing Gas Volumes Throughout the Course of VMAT Radiation Therapy of Upper Gastrointestinal Tumors

Purpose

This retrospective patient study assessed the consistency of abdominal gas presence throughout radiation therapy for patients with upper gastrointestinal cancer and determined the impact of variations in gas volume on the calculated dose distribution of volumetric modulated arc therapy.

Methods and Materials

Eight patients with pancreatic cancer were included for analysis. A plan library consisting of 3 reference plans per patient (Ref_0.0, Ref_0.5, and Ref_1.0) was created based on planning computed tomography (CT) with density overrides of 0.0, 0.5, and 1.0 applied to gas volumes, respectively. Corresponding cone beam CT (CBCT) data sets were obtained and density overrides were applied to enable fractional dose calculation. Variation in gas volume relative to initial volume determined from CT was assessed. Dose metrics for targets and organs at risk were compared between the accumulated CBCT dose and the planned dose of the 3 reference plans for each patient.

Results

There was a significant decrease in gas present from CT to treatment CBCT, with a mean decrease in volume of 48.6% for the entire cohort. Dosimetrically, all accumulated target and organ-at-risk parameters, aside from the kidneys, exhibited the smallest mean deviation from the Ref_0.0 plan and largest mean deviation from the Ref_1.0 plan. A statistically significant difference in mean accumulated dose to Ref_0.0 and Ref_1.0 was observed for the dose delivered to 95% of the planning target volume.

Conclusions

Significant variation in gas volumes from CT to treatment can occur throughout volumetric modulated arc therapy for pancreatic cancer. Through the use of a plan library, it was determined that initial assessment of a patient's treatment plan with an assigned gas density of 0.0 provided the most accurate prediction of the accumulated dose.

Impact of interfractional anatomical variation and setup correction methods on interfractional dose variation in IMPT and VMAT plans for pancreatic cancer patients: A planning study

Purpose

To investigate the impact of interfractional anatomical changes and setup correction methods on dose distributions in pancreatic cancer patients under breath‐hold conditions.

Methods

Three intensity‐modulated proton therapy (IMPT) plans with different beam arrangements and one volumetric‐modulated arc therapy (VMAT) plan prescribing 54 Gy in 30 fractions were created for 10 patients who underwent three additional CT scans performed at an interval of 1–2 weeks. The additional CT sets were rigidly registered to the simulation CT set using both bone‐matching (BM) and organ‐matching (OM) methods in each patient. Recalculated dose distributions and dose–volume indices on the additional CT sets using either the BM or the OM method were compared with the simulation values.

Results

Differences in the gross tumor volume D_98% value from the simulation sets ranged from −0.8 to −5.9% on average. In addition, the variations were larger with OM compared with BM for two IMPT plans. Meanwhile, differences in the D_98% value in the region isotropically enlarged by 5 mm from the gross tumor volume were significantly improved with OM on two IMPT plans and the VMAT plan. Among the organs at risk, the dose–volume indices were significantly improved with OM only in the duodenum on all plans.

Conclusion

Organ‐matching may be a better setup correction technique than BM for both photon therapy and IMPT plans. However, in some beam arrangements of IMPT, the dose distribution may be somewhat worse using OM, due to interfractional anatomical variation. Therefore, it is important to choose beam angles that are less likely to be influenced by changes in the gastrointestinal gas volume, especially in IMPT plans.

Method for a motion model based automated 4D dose calculation

The Vero system can treat intra-fractionally moving tumors with gimbaled dynamic tumor tracking (DTT) by rotating the treatment beam so that it follows the motion of the tumor. However, the changes in the beam geometry and the constant breathing motion of the patient influence the dose applied to the patient. This study aims to perform a full 4D dose reconstruction for thirteen patients treated with DTT at the Vero system at the Universitätsklinikum Erlangen and investigates the temporal resolution required to perform an accurate 4D dose reconstruction. For all patients, a 4DCT was used to train a 4D motion model, which is able to calculate pseudo-CT images for arbitrary breathing phases. A new CT image was calculated for every 100 ms of treatment and a dose calculation was performed according to the current beam geometry (i.e. the rotation of the treatment beam at this moment in time) by rotating according to the momentary beam rotation, which is extracted from log-files. The resulting dose distributions were accumulated on the planning CT and characteristic parameters were extracted and compared. [Formula: see text]-evaluations of dose accumulations with different spatial-temporal resolutions were performed to determine the minimal required resolution. In total 173 700 dose calculations were performed. The accumulated 4D dose distributions show a reduced mean GTV dose of 0.77% compared to the static treatment plan. For some patients larger deviations were observed, especially in the presence of a poor 4DCT quality. The [Formula: see text]-evaluation showed that a temporal resolution of 500 ms is sufficient for an accurate dose reconstruction. If the tumor motion is regarded as well, a spatial-temporal sampling of 1400 ms and 2 mm yields accurate results, which reduces the workload by 84%.

Effect of Changing Phantom Thickness on Helical Radiotherapy Plan: Dosimetric Analysis

Purpose: The aim of this study is to investigate the effect of changing phantom thickness on high dose region of interest (HD_ROI) and low dose ROI’s (LW_ROI’s) doses during helical radiotherapy (RT) by utilizing Adaptive RT (ART) technique.

Materials and Methods: The cylindrical phantom (CP) is wrapped with different thickness boluses and scanned in the kilovoltage computed tomography (KVCT). HD_ROI and LW_ROI’s were created in contouring system and nine same plans (1.8 Gy/Fr) were made with images of different thicknesses CP. The point dose measurements were performed using ionization chamber in Helical Tomotherapy (HT) treatment machine. For detecting thickness reduction effect, CP was irradiated using bolus-designed plans and it was irradiated using without bolus plan. The opposite of this scenario was applied to determine the thickness increase. KVCT and megavoltage CT (MVCT) images were used for dose comparison. The HT Planned Adaptive Software was used to see the differences in the planning and verification doses at dose volume histograms (DVH).

Results: Point dose measurements showed a 4.480% dose increase in 0.5 cm depth reduction for HD_ROI. These differences reached 8.508% in 2 cm depth and 15,279% in 5 cm depth. At the same time, a dose reduction of 0.665% was determined for a 0.5cm depth increase, a dose reduction of 1.771% was determined for a 2 cm depth increase, a dose reduction of 5.202% was determined for a 5 cm depth increase for the HD_ROI. The ART plan results show that the dose changes in the HD_ROI was greater than the LW_ROI’s.

Conclusion: Phantom thicknesses change can lead to a serious dose increase or decrease in the HD_ROI and LW_ROI’s.