Dosimetric variations due to interfraction organ deformation in cervical cancer brachytherapy

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.

The impact of an educational tool in cervix image registration across three imaging modalities

OBJECTIVES

Accurate image registration is vital in cervical cancer where changes in both planning target volume (PTV) and organs at risk (OARs) can make decisions regarding image registration complicated. This work aims to determine the impact of a dedicated educational tool compared with experience gained in MR-guided radiotherapy (MRgRT).

METHODS

10 therapeutic radiographers acted as observers and were split into two groups based on previous experience with MRgRT and Monaco treatment planning system. Three CBCT-CT, three MR-CT and two MR-MR registrations were completed per patient by each observer. Observers recorded translations, time to complete image registration and confidence. Data were collected in two phases; prior to and following the introduction of a cervix registration guide.

RESULTS

No statistically significant differences were noted between imaging modalities. Each group was assessed independently pre- and post-education, no statistically significant differences were noted in either CBCT-CT or MR-CT imaging. Group 1 MR-MR imaging showed a statistically significant reduction in interobserver variability (p=0.04), in Group 2, the result was not statistically significant (p=0.06). Statistically significant increases in confidence were seen in all three modalities (p≤0.05).

CONCLUSIONS

At The Christie NHS Foundation Trust, radiographers consistently registered images across three different imaging modalities regardless of their previous experience. The implementation of an image registration guide had limited impact on inter- and intraobserver variability. Radiographers' confidence showed statistically significant improvements following the use of the registration manual.

ADVANCES IN KNOWLEDGE

This work helps evaluate training methods for novel roles that are developing in MRgRT.

Estimating the accumulative dose uncertainty for intracavitary and interstitial brachytherapy

BACKGROUND

Image-guided adaptive brachytherapy shows the ability to deliver high doses to tumors while sparing normal tissues. However, interfraction dose delivery introduces uncertainties to high dose estimation, which relates to normal tissue toxicity. The purpose of this study was to investigate the high-dose regions of two applicator approaches in brachytherapy.

METHOD

For 32 cervical cancer patients, the CT images from each fraction were wrapped to a reference image, and the displacement vector field (DVF) was calculated with a hybrid intensity-based deformable registration algorithm. The fractional dose was then accumulated to calculate the position and the overlap of high dose (D2cc) during multiple fractions.

RESULT

The overall Dice similarity coefficient (DSC) of the deformation algorithm for the bladder and the rectum was (0.97 and 0.91). No significant difference was observed between the two applicators. However, the location of the intracavitary brachytherapy (ICBT) high-dose region was relatively concentrated. The overlap volume of bladder and rectum D2cc was 0.42 and 0.71, respectively, which was higher than that of interstitial brachytherapy (ISBT) (0.26 and 0.31). The cumulative dose was overestimated in ISBT cases when using the GEC-recommended method. The ratio of bladder and rectum D2cc to the GEC method was 0.99 and 1, respectively, which was higher than that of the ISBT method (0.96 and 0.94).

CONCLUSION

High-dose regions for brachytherapy based on different applicator types were different. The 3D-printed ICBT has better high-dose region consistency than freehand ISBT and hence is more predictable.

Tensor Regression-based Model to Investigate Heterogeneous Spatial Radiosensitivity After I-125 Seed Implantation for Prostate Cancer

BACKGROUND/AIM

We established a data-driven method for extracting spatial patterns of dose distribution associated with radiation injuries, based on patients with prostate cancer who underwent iodine-125 (I-125) seed implantation.

PATIENTS AND METHODS

Seventy-five patients underwent I-125 seed implantation for prostate cancer. We modeled the severity of lower urinary tract symptoms (LUTS) to be estimated using a linear model, which is formulated as an inner product between the dose distribution D and voxel-wise radiosensitivity B inside the prostate. For the estimation, tensor regression based on a low-rank decomposition with generalized fused lasso penalty was applied.

RESULTS

The spatial distribution of B was visually assessed. Positive parameters appeared dominantly in the region close to the urethra and the prostate base.

CONCLUSION

Our tensor regression-based model can predict intra-organ radiosensitivity in a data-driven manner, providing a compelling parameter distribution associated with the development of LUTS after I-125 seed implantation for prostate cancer.

Single Application Multifractionated Image Guided Adaptive High-Dose-Rate Brachytherapy for Cervical Cancer: Dosimetric and Clinical Outcomes

PURPOSE

A prospective phase 2 study was conducted to evaluate the feasibility and safety of single-application multifractionated (SA-MF), high-dose-rate (HDR), image guided adaptive brachytherapy (IGABT) for cervical cancer.

METHODS AND MATERIALS

Patients (N = 41) with International Federation of Gynaecology and Obstetrics 2009 stage IIB-IVA disease recruited between 2017 and 2019 underwent SA-MF. After completion of external beam radiation therapy of 50 Gy in 25 fractions, patients received magnetic resonance IGABT. The IGABT protocol consisted of a single brachytherapy (BT) application and treatment with 3 fractions of HDR (9 Gy on day 1; 2 fractions of 7 Gy with a minimum 6-hour gap on day 2) after achieving planning aims of the high-risk clinical target volume (HRCTV) receiving >84 Gy EQD2 and 2 cm³ of the bladder and rectum/sigmoid receiving ≤85 Gy and <71 Gy, respectively. Interfraction variation was addressed by performing computed tomography planning and coregistration using a mutual information-based coordinate system on day 2 before the second fraction. Organ at risk contouring was done on computed tomography, and doses were re-evaluated and reoptimized if required.

RESULTS

Thirty-eight patients were treated as per the protocol. All patients underwent Intracavitary + Interstitial BT with needles (median, 4; range, 3-11). The mean ± standard deviation HRCTV volume was 41 ± 21 cm³ and HRCTV D90 dose was 87.2 ± 3.6Gy. The 0.1 cm³ and 2 cm³ to bladder, rectum, and sigmoid were -103.2 ± 10.6 Gy and -84.6 ± 6.8 Gy, 82.2 ± 9.5 Gy and -68.3 ± 5.7 Gy, and 83.5 ± 9.8 Gy and -69.5 ± 5.9 Gy, respectively. Six patients required reoptimization before the second fraction to meet planning aims. Mean overall treatment time was 47 ± 6 days. With a median follow up of 22 months (range, 2-37), 2-year local control and disease-free and overall survival were 90.1%, 85%, and 94.5%, respectively. So far 1 patient with grade II and 2 patients with grade III rectal toxicities have been reported.

CONCLUSION

Magnetic resonance IGABT with SA-MF BT was feasible in 95% of patients. The dosimetric parameters and clinical results achieved so far look promising.

Comparison of predictive performance for toxicity by accumulative dose of DVH parameter addition and DIR addition for cervical cancer patients

We compared predictive performance between dose volume histogram (DVH) parameter addition and deformable image registration (DIR) addition for gastrointestinal (GI) toxicity in cervical cancer patients. A total of 59 patients receiving brachytherapy and external beam radiotherapy were analyzed retrospectively. The accumulative dose was calculated by three methods: conventional DVH parameter addition, full DIR addition and partial DIR addition. ${D}_{2{cm}^3}$, ${D}_{1{cm}^3}$ and ${D}_{0.1{cm}^3}$ (minimum doses to the most exposed 2 cm3, 1cm3 and 0.1 cm3 of tissue, respectively) of the rectum and sigmoid were calculated by each method. V50, V60 and V70 Gy (volume irradiated over 50, 60 and 70 Gy, respectively) were calculated in full DIR addition. The DVH parameters were compared between toxicity (≥grade1) and non-toxicity groups. The area under the curve (AUC) of the receiver operating characteristic (ROC) curves were compared to evaluate the predictive performance of each method. The differences between toxicity and non-toxicity groups in ${D}_{2{cm}^3}$ were 0.2, 5.7 and 3.1 Gy for the DVH parameter addition, full DIR addition and partial DIR addition, respectively. The AUCs of ${D}_{2{cm}^3}$ were 0.51, 0.67 and 0.57 for DVH parameter addition, full DIR addition and partial DIR addition, respectively. In full DIR addition, the difference in dose between toxicity and non-toxicity was the largest and AUC was the highest. AUCs of V50, V60 and V70 Gy were 0.51, 0.63 and 0.62, respectively, and V60 and V70 were high values close to the value of ${D}_{2{cm}^3}$ of the full DIR addition. Our results suggested that the full DIR addition may have the potential to predict toxicity more accurately than the conventional DVH parameter addition, and that it could be more effective to accumulate to all pelvic irradiation by DIR.

Predictive value of Excel forms based on an automatic calculation of dose equivalent in 2 Gy per fraction in adaptive brachytherapy for cervical cancer

Purpose

External beam radiotherapy (EBRT) combined with brachytherapy (BT) is the standard mode of radical radiotherapy for locally advanced cervical cancer. The cumulative equivalent doses in 2 Gy per fraction (EQD₂) is an important basis for estimating the probability of local control of tumors and monitoring the occurrence of side effects in normal tissues. The purpose of this study was to explore the predictive value of Excel forms based on an automatic calculation in radical adaptive BT for cervical cancer.

Material and methods

A retrospective analysis of 119 patients suffering from cervical cancer, treated with radical radiotherapy. All patients were treated with EBRT and adaptive BT. EBRT prescribed dose was 42.0-50.4 Gy in 21-28 fractions. BT nominal prescribed dose was 28 Gy in 4 fractions, separated by one week. Total EQD₂ prediction at nth (n = 1-3) BT (TEPB_n) or actual cumulative EQD₂ (ACEQD₂) can be calculated automatically by inputting the physical dose based on an in-house designed application. The relationship between TEPB_n and ACEQD₂ was evaluated, and the predictive value of Excel forms based on the automatic calculation was analyzed.

Results

For the volume of high-risk clinical target, there was a significant decrease between BT1 and BT2. Similarly, for the volume of intermediate-risk clinical target, there was a significant decrease between BT2 and BT3. The sensitivity ranges of TEPB₁, TEPB₂, and TEPB₃ prediction were 74.5-91.3%, 83.7-95.7%, and 92.9-99.1%, respectively, and the specificity ranges were 46.7-80.0%, 53.3-90.5%, and 66.7-90.5%, respectively.

Conclusions

The in-house designed application has the function of quickly reading dose-volume histogram (DVH) parameters from the treatment planning system, which allows for balance between the total dose to target volumes and organs at risk (OARs). Excel forms based on EQD₂ automatic calculation presents high predictive accuracy.

A Population-based Statistical Model for Investigating Heterogeneous Intraprostatic Sensitivity to Radiation Toxicity After 125I Seed Implantation

AIM

To develop a population-based statistical model in order to find a spatial pattern of dose distribution which is related to lower urinary tract symptoms (LUTS) after iodine-125 (¹²⁵I) seed implantation for prostate cancer.

PATIENTS AND METHODS

A total of 75 patients underwent ¹²⁵I seed implantation for prostate cancer. Principal component analysis was applied to the standardized dose array and for each patient dose distribution was uniquely characterized by a combination of weighted eigenvectors. The correlation between eigenvectors and the severity of LUTS was investigated with linear regression analysis.

RESULTS

Eight eigenvectors were identified as being significantly associated with the severity of LUTS (p<0.05). Multivariate regression model identified that intraprostatic parameters, which were positively associated with the severity of LUTS, were distributed around a portion of the urethral base and a peripheral region of the prostate.

CONCLUSION

We established a population-based statistical model that may indicate a significant dose pattern associated with the severity of radiation toxicity.

Evaluation of deformable image registration algorithm for determination of accumulated dose for brachytherapy of cervical cancer patients

Purpose

This study was designed to assess the dose accumulation (DA) of bladder and rectum between brachytherapy fractions using hybrid-based deformable image registration (DIR) and compare it with the simple summation (SS) approach of GEC-ESTRO in cervical cancer patients.

Material and methods

Patients (n = 137) with cervical cancer treated with 3D conformal radiotherapy and three fractions of high-dose-rate brachytherapy were selected. CT images were acquired to delineate organs at risk and targets according to GEC-ESTRO recommendations. In order to determine the DA for the bladder and rectum, hybrid-based DIR was done for three different fractions of brachytherapy and the results were compared with the standard GEC-ESTRO method. Also, we performed a phantom study to calculate the uncertainty of the hybrid-based DIR algorithm for contour matching and dose mapping.

Results

The mean ± standard deviation (SD) of the Dice similarity coefficient (DICE), Jaccard, Hausdorff distance (HD) and mean distance to agreement (MDA) in the DIR process were 0.94 ±0.02, 0.89 ±0.03, 8.44 ±3.56 and 0.72 ±0.22 for bladder and 0.89 ±0.05, 0.80 ±0.07, 15.46 ±10.14 and 1.19 ±0.59 for rectum, respectively. The median (Q1, Q3; maximum) Gy_EQD2 differences of total D_2cc between DIR-based and SS methods for the bladder and rectum were reduced by –1.53 (–0.86, –2.98; –9.17) and –1.38 (–0.80, –2.14; –7.11), respectively. The mean ± SD of DICE, Jaccard, HD, and MDA for contour matching were 0.98 ±0.008, 0.97 ±0.01, 2.00 ±0.70 and 0.20 ±0.04, respectively for large deformation. Maximum uncertainty of dose mapping was about 3.58%.

Conclusions

The hybrid-based DIR algorithm demonstrated low registration uncertainty for both contour matching and dose mapping. The DA difference between DIR-based and SS approaches was statistically significant for both bladder and rectum and hybrid-based DIR showed potential to assess DA between brachytherapy fractions.

Deformable image registration for radiation therapy: principle, methods, applications and evaluation

Background: Deformable image registration (DIR) is increasingly used in the field of radiation therapy (RT) to account for anatomical deformations. The aims of this paper are to describe the main applications of DIR in RT and discuss current DIR evaluation methods. Methods: Articles on DIR published from January 2000 to October 2018 were extracted from PubMed and Science Direct. Our search was restricted to articles that report data obtained from humans, were written in English, and address DIR methods for RT. A total of 207 articles were selected from among 2506 identified in the search process. Results: At planning, DIR is used for organ delineation using atlas-based segmentation, deformation-based planning target volume definition, functional planning and magnetic resonance imaging-based dose calculation. In image-guided RT, DIR is used for contour propagation and dose calculation on per-treatment imaging. DIR is also used to determine the accumulated dose from fraction to fraction in external beam RT and brachytherapy, both for dose reporting and adaptive RT. In the case of re-irradiation, DIR can be used to estimate the cumulated dose of the two irradiations. Finally, DIR can be used to predict toxicity in voxel-wise population analysis. However, the evaluation of DIR remains an open issue, especially when dealing with complex cases such as the disappearance of matter. To quantify DIR uncertainties, most evaluation methods are limited to geometry-based metrics. Software companies have now integrated DIR tools into treatment planning systems for clinical use, such as contour propagation and fraction dose accumulation. Conclusions: DIR is increasingly important in RT applications, from planning to toxicity prediction. DIR is routinely used to reduce the workload of contour propagation. However, its use for complex dosimetric applications must be carefully evaluated by combining quantitative and qualitative analyses.

Deformable image registration for dose mapping between external beam radiotherapy and brachytherapy images of cervical cancer

For locally advanced cervical cancer (LACC), anatomy correspondence with and without BT applicator needs to be quantified to merge the delivered doses of external beam radiation therapy (EBRT) and brachytherapy (BT). This study proposed and evaluated different deformable image registration (DIR) methods for this application. Twenty patients who underwent EBRT and BT for LACC were retrospectively analyzed. Each patient had a pre-BT CT at EBRT boost (without applicator) and a CT and MRI at BT (with applicator). The evaluated DIR methods were the diffeomorphic Demons, commercial intensity and hybrid methods, and three different biomechanical models. The biomechanical models considered different boundary conditions (BCs). The impact of the BT devices insertion on the anatomy was quantified. DIR method performances were quantified using geometric criteria between the original and deformed contours. The BT dose was deformed toward the pre-CT BT by each DIR method. The impact of boundary conditions to drive the biomechanical model was evaluated based on the deformation vector field and dose differences. The GEC-ESTRO guideline dose indices were reported. Large organ displacements, deformations, and volume variations were observed between the pre-BT and BT anatomies. Rigid registration and intensity-based DIR resulted in poor geometric accuracy with mean Dice similarity coefficient (DSC) inferior to 0.57, 0.63, 0.42, 0.32, and 0.43 for the rectum, bladder, vagina, cervix and uterus, respectively. Biomechanical models provided a mean DSC of 0.96 for all the organs. By considering the cervix-uterus as one single structure, biomechanical models provided a mean DSC of 0.88 and 0.94 for the cervix and uterus, respectively. The deformed doses were represented for each DIR method. Caution should be used when performing DIR for this application as standard techniques may have unacceptable results. The biomechanical model with the cervix-uterus as one structure provided the most realistic deformations to propagate the BT dose toward the EBRT boost anatomy.

Dose comparison between prescription methods according to anatomical variations in intracavitary brachytherapy for cervical cancer

PURPOSE

We compared how doses delivered via two-dimensional (2D) intracavitary brachytherapy (ICBT) and three-dimensional (3D) ICBT varied anatomically.

Materials and METHODS

A total of 50 patients who received 30 Gy of 3D ICBT after external radiotherapy (RT) were enrolled. We compared the doses of the actual 3D and 2D ICBT plans among patients grouped according to six anatomical variations: differences in a small-bowel V2Gy, small bowel circumference, the direction of bladder distension, bladder volume, sigmoid V3.5Gy, and sigmoid circumference. Seven dose parameters were measured in line with the EMBRACE recommendations.

RESULTS

In terms of bladder volume, the bladder and small-bowel D2cc values were lower in the 150-250 mL bladder volume subgroup; and the rectum, sigmoid, and bladder D2mL values were all lower in the >250 mL subgroup, for 3D vs. 2D ICBT. In the sigmoid V3.5Gy >2 mL subgroup, the sigmoid and bladder D2mL values were significantly lower for 3D than 2D ICBT. The bladder D2mL value was also significantly lower for 3D ICBT, as reflected by the sigmoid circumference. In patients with a small bowel V2.0Gy >10 mL or small bowel circumference >15%, most dose parameters were significantly lower for 3D than 2D ICBT. The bladder distension direction did not significantly affect the doses.

CONCLUSION

Compared to 2D ICBT, a greater bladder volume can reduce the internal 3D ICBT organ dose without affecting the target dose.

Radiobiological considerations in combining doses from external beam radiotherapy and brachytherapy for cervical cancer

The recommended radio-therapeutic treatment for cervix cancer consists of a first phase of external beam radiotherapy (EBRT) plus a second phase of brachytherapy (BT), the combined treatment being delivered within 8 weeks. In order to assess a comprehensive dosimetry of the whole treatment, it is necessary to take into account that these two phases are characterized by different spatial and temporal dosimetric distributions, which complicates the task of the summation of the two contributions, EBRT and BT. Radiobiology allows to tackle this issue pragmatically by means of the LQ model and, in fact, this is the usual tool currently in use for this matter. In this work, we describe the rationale behind the summation of the dosimetric contributions of the two phases of the treatment, EBRT and BT, for cervix cancer, as carried out with the LQ model. Besides, we address, from a radiobiological point of view, several important considerations regarding the use of the LQ model for this task. One of them is the analysis of the effect of the overall treatment time in the result of the global treatment. Another important question considered is related to the fact that the capacity of LQ to predict the treatment outcomes is deteriorated when the dose per fraction of the radiotherapic scheme exceeds 6-10 Gy, which is a typical brachytherapy fractionation. Finally, we analyze the influence of the uncertainty and the variability of the main parameters utilized in the LQ model formulation in the assessment of the global dosimetry.

Accumulation of rectum dose-volume metrics for prostate external beam radiotherapy combined with brachytherapy: Evaluating deformably registered dose distribution addition using parameter-based addition

INTRODUCTION

To investigate the accuracy of deriving dose-volume histogram (DVH) parameters from deformably registered data by comparing values with the simple addition of DVHs from each phase of a combined external beam radiotherapy (EBRT)/high-dose-rate (HDR-BT) brachytherapy prostate treatment.

METHODS

Eighty-two patients received EBRT in 23 fractions of 2 Gy and HDR-BT TG43 in three fractions of 6.5 Gy. The HDR-BT CT was deformably registered to the EBRT CT. The rectum D_0.1cc , D_1cc , D_2cc and D_10cc were calculated in two ways. (i) Parameter-adding: the EBRT DVH parameters (or the EBRT prescription dose) were added to the unregistered HDR-BT DVH parameters. (ii) Distribution-adding: the parameters were extracted after the EBRT doses were 3D-summed with the registered HDR-BT doses. Resulting differences between the parameters were investigated.

RESULTS

The D_0.1cc , D_1cc and D_2cc from parameter-adding were 21.3% (P < 0.001), 6.3% (P < 0.001) and 3.5% (P < 0.001) smaller than those from distribution-adding. The D_10cc was 2.2% (P = 0.015) larger for distribution-adding.

CONCLUSION

Distribution-adding was confounded by unsystematic inter/intra-observer rectum-contouring errors and registration accuracy near the anterior rectal wall. Consequently, clinical use of distribution-adding to assess rectal doses requires careful contour and registration evaluation.

A comprehensive evaluation of adaptive daily planning for cervical cancer HDR brachytherapy

The purpose of this study was to evaluate adaptive daily planning for cervical cancer patients who underwent high‐dose‐rate intracavitary brachytherapy (HDR‐BT) using comprehensive interfractional organ motion measurements. This study included 22 cervical cancer patients who underwent 5 fractions of HDR‐BT. Regions of interest (ROIs) including high‐risk clinical tumor volume (HR‐CTV) and organs at risk (OARs) were manually contoured on daily CT images. All patients were clinically treated with adaptive daily plans (ADP), which involved ROI delineation and dose optimization at each treatment fraction. Single treatment plans (SP) were retrospectively generated by applying the first treatment fraction's dwell times adjusted for decay and dwell positions of the applicator to subsequent treatment fractions. Various existing similarity metrics were calculated for the ROIs to quantify interfractional organ variations. A novel similarity (JRARM) score was established, which combined both volumetric overlap metrics (DSC, JSC, and RVD) and distance metrics (ASD, MSD, and RMSD). Linear regression was performed to determine a relationship between interfractional organ variations of various similarity metrics and D2cc variations from both plans. Wilcoxon signed‐rank tests were used to assess ADP and SP by comparing EQD2D2cc(α/β=3) for OARs. For interfractional organ variations, the sigmoid demonstrated the greatest variations based on the JRARM, DSC, and RMSD metrics. Comparisons between paired ROIs showed differences in metrics at each treatment fraction. RVD, MSD, and RMSD were found to be significantly correlated to D2cc variations for bladder and sigmoid. The comparison between plans found ADP provided lower EQD2 D2cc of OARs than SP. Specifically, the sigmoid demonstrated statistically significant dose variations (p=0.015). Substantial interfractional organ motion occurs during HDR‐BT based on comprehensive measurements and may significantly affect D2cc of OARs. Adaptive daily planning provides improved dose sparing for OARs compared to single planning with the extent of sparing being different among OARs.

PACS number(s): 87.55.D, 87.55.de, 87.55.kh, 87.57.nj

In vivo spatial correlation between (18)F-BPA and (18)F-FDG uptakes in head and neck cancer

BACKGROUND AND PURPOSE

Borono-2-(18)F-fluoro-phenylalanine ((18)F-BPA) has been used to estimate the therapeutic effects of boron neutron capture therapy (BNCT), while (18)F-fluorodeoxyglucose ((18)F-FDG) is the most commonly used positron emission tomography (PET) radiopharmaceutical in a routine clinical use. The aim of the present study was to evaluate spatial correlation between (18)F-BPA and (18)F-FDG uptakes using a deformable image registration-based technique.

MATERIAL AND METHODS

Ten patients with head and neck cancer were recruited from January 2014 to December 2014. All patients underwent whole-body (18)F-BPA PET/computed tomography (CT) and (18)F-FDG PET/CT within a 2-week period. For each patient, (18)F-BPA PET/CT and (18)F-FDG PET/CT images were aligned based on a deformable image registration framework. The voxel-by-voxel spatial correlation of standardized uptake value (SUV) within the tumor was analyzed.

RESULTS

Our image processing framework achieved accurate and validated registration results for each PET/CT image. In 9/10 patients, the spatial distribution of SUVs between (18)F-BPA and (18)F-FDG showed a significant, positive correlation in the tumor volume.

CONCLUSIONS

Deformable image registration-based voxel-wise analysis demonstrated a spatial correlation between (18)F-BPA and (18)F-FDG uptakes in the head and neck cancer. A tumor sub-volume with a high (18)F-FDG uptake may predict high accumulation of (18)F-BPA.

Dose reconstruction technique using non-rigid registration to evaluate spatial correspondence between high-dose region and late radiation toxicity: a case of tracheobronchial stenosis after external beam radiotherapy combined with endotracheal brachytherapy for tracheal cancer

Purpose

Small organ subvolume irradiated by a high-dose has been emphasized to be associated with late complication after radiotherapy. Here, we demonstrate a potential use of surface-based, non-rigid registration to investigate how high-dose volume topographically correlates with the location of late radiation morbidity in a case of tracheobronchial radiation stenosis.

Material and methods

An algorithm of point set registration was implemented to handle non-rigid registration between contour points on the organ surfaces. The framework estimated the global correspondence between the dose distribution and the varying anatomical structure. We applied it to an 80-year-old man who developed tracheobronchial stenosis 2 years after high-dose-rate endobronchial brachytherapy (HDR-EBT) (24 Gy in 6 Gy fractions) and external beam radiotherapy (EBRT) (40 Gy in 2 Gy fractions) for early-stage tracheal cancer.

Results and conclusions

Based on the transformation function computed by the non-rigid registration, irradiated dose distribution was reconstructed on the surface of post-treatment tracheobronchial stenosis. For expressing the equivalent dose in a fractional dose of 2 Gy in HDR-EBT, α/β of linear quadratic model was assumed as 3 Gy for the tracheal bronchus. The tracheobronchial surface irradiated by more than 100 Gy_αβ3 tended to develop severe stenosis, which attributed to a more than 50% decrease in the luminal area. The proposed dose reconstruction technique can be a powerful tool to predict late radiation toxicity with spatial consideration.