Evaluation of offline adaptive planning techniques in image-guided brachytherapy of cervical cancer

Method to assess the need for re-planning HDR brachytherapy tandem and ring treatments

High dose rate (HDR) brachytherapy procedures for cervical cancer require multiple applicator insertions for multiple (typically 5) fractions of a single plan, which carries a risk for variability in applicator position between fractions. Due to applicator displacement relative to patient anatomy, the dose to nearby organs-at-risk (OARs) may vary significantly from one fraction to the next. The purpose of this study was to evaluate the effect of changes in HDR tandem and ring (T&R) applicator position on doses to nearby OARs and to present a quick and simple method to estimate doses to OARs inter-fractionally without having to perform a re-plan. Ninety CT image sets for 20 patients, ages 44 to 86, undergoing T&R-based HDR for cervical cancer were used retrospectively for this study. Measures of applicator positional and angular changes relative to the bony anatomy were obtained using image fusion in MIM software, between the planning CT (plan CT) and the CT on the treatment day (CT-TX). Dosimetric data were determined, also using MIM software, using the original (first fraction) dose distribution applied to organs at risk (rectum and bladder), transferred via rigid registration from the plan CT to each CT-TX. Bladder and rectum contours were also transferred from each plan CT to each CT-TX and were tweaked manually to match anatomy on each CT-TX and examined visually for appropriateness. Differences in translation and rotation of the T&R applicator between the planning CT and subsequent individual fractions were recorded and plotted against dose differences between each fraction of treatment and the original (first) fraction. Absolute dose (D2cc) and volume (V50) differences vs positional shifts were calculated and plotted, and the Pearson Product-Moment correlation coefficient between dose parameters and measured positional shifts was determined. Average dosimetric differences between planned dose and subsequent fractional doses obtained through rigid registration were 1.48 ± 1.92 Gy, 14.91 ± 11.92 cm3, 0.56 ± 0.93 Gy, and 1.77 ± 2.18 cm3 for Bladder D2cc, Bladder V50, Rectum D2cc, and Rectum V50, respectively. Correlation between Bladder V50 and sagittal plane rotation gave an r2 of 0.4, showing the most correlation of all parameters studied. Bladder dose and volume increased by a maximum of about 2.7 Gy and 50 cm3 overall for Bladder D2cc and Bladder V50, respectively. Bladder V50 was most sensitive to T&R applicator displacements. We have quantified the effects of applicator positional changes on dose changes for the bladder and rectum. Even large changes in applicator position between fractions did not result in significant changes in dose to these normal tissues, indicating that adaptive re-planning is not necessary.

Effects of Interfraction Dose Variations of Target and Organs at Risk on Clinical Outcomes in High Dose Rate Brachytherapy for Cervical Cancer

Meeting dose prescription is critical to control tumors in radiation therapy. Interfraction dose variations (IDVs) from the prescribed dose in high dose rate brachytherapy (HDR) would cause the target dose to deviate from the prescription but their clinical effect has not been widely discussed in the literature. Our previous study found that IDVs followed a left-skewed distribution. The clinical effect of the IDVs in 100 cervical cancer HDR patients will be addressed in this paper. An in-house Monte Carlo (MC) program was used to simulate clinical outcomes by convolving published tumor dose response curves with IDV distributions. The optimal dose and probability of risk-free local control (RFLC) were calculated using the utility model. The IDVs were well-fitted by the left-skewed Beta distribution, which caused a 3.99% decrease in local control probability and a 1.80% increase in treatment failure. Utility with respect to IDV uncertainty increased the RFLC probability by 6.70% and predicted an optimal dose range of 83 Gy-91 Gy EQD2. It was also found that a 10 Gy dose escalation would not affect toxicity. In conclusion, HRCTV IDV uncertainty reduced LC probabilities and increased treatment failure rates. A dose escalation may help mitigate such effects.

Statistical Analysis of Interfraction Dose Variations of High-Risk Clinical Target Volume and Organs at Risk for Cervical Cancer High-Dose-Rate Brachytherapy

Purpose

High-dose-rate (HDR) brachytherapy for cervical cancer treatment includes significant uncertainties. The aim of this study was to quantify the interfraction dosimetric variation (IDV) of the high-risk clinical target volume (HRCTV) from the prescribed dose and the corresponding effect on organ-at-risk (OAR) dose based on a comprehensive statistical analysis.

Methods and Materials

Fifty patients with cervical cancer treated with high-dose-rate intracavity brachytherapy from October 2019 to December 2020 were retrospectively analyzed. The OARs of interest were the rectum, bladder, sigmoid, and bowel. The dosimetric parameters evaluated for all patients was the dose absorbed by 90% of the HRCTV ( D 90 ) and the dose absorbed by 0.1 ( D 0.1 c c ) and 2 cm3 ( D 2 c c ) of each respective OAR. The HRCTV variations were from the prescribed dose and the OAR variations were from the corresponding tolerance dose. Distribution fitting of the HRCTV variations was determined to quantify the IDV. Comparative statistics of the HRCTV variations with the OAR variations were conducted to determine correlations.

Results

The mean HRCTV variation from the prescribed dose was -2.53% ± 8.74%. The HRCTV variations and OAR variations showed moderate to weak linear correlations despite the variations being relative to each other, with the bladder D 2 c c having the strongest correlation. There was a 30.0% (±2.62%, 95% confidence interval) probability of underdosing the HRCTV (-5% variation from prescription) and a 23.3% (±2.62%, 95% confidence interval) probability of overdosing the HRCTV (+5% variation from prescription). This tendency to underdose the HRCTV was a consequence of HRCTV IDV not being normally distributed.

Conclusions

HRCTV dosimetric variations and OAR variations were complexly correlated with the bladder D 2 c c having the strongest correlation. HRCTV IDV was best described as a left-skewed distribution that indicates a tendency of underdosing the HRCTV. The clinical significance of such dose variations is expected and will be further investigated.