Dosimetric benefits of daily treatment plan adaptation for prostate cancer stereotactic body radiotherapy
Radiation Oncology volume 16, Article number: 145 (2021)
Hypofractionation is increasingly being applied in radiotherapy for prostate cancer, requiring higher accuracy of daily treatment deliveries than in conventional image-guided radiotherapy (IGRT). Different adaptive radiotherapy (ART) strategies were evaluated with regard to dosimetric benefits.
Treatments plans for 32 patients were retrospectively generated and analyzed according to the PACE-C trial treatment scheme (40 Gy in 5 fractions). Using a previously trained cycle-generative adversarial network algorithm, synthetic CT (sCT) were generated out of five daily cone-beam CT. Dose calculation on sCT was performed for four different adaptation approaches: IGRT without adaptation, adaptation via segment aperture morphing (SAM) and segment weight optimization (ART1) or additional shape optimization (ART2) as well as a full re-optimization (ART3). Dose distributions were evaluated regarding dose-volume parameters and a penalty score.
Compared to the IGRT approach, the ART1, ART2 and ART3 approaches substantially reduced the V37Gy(bladder) and V36Gy(rectum) from a mean of 7.4cm3 and 2.0cm3 to (5.9cm3, 6.1cm3, 5.2cm3) as well as to (1.4cm3, 1.4cm3, 1.0cm3), respectively. Plan adaptation required on average 2.6 min for the ART1 approach and yielded doses to the rectum being insignificantly different from the ART2 approach. Based on an accumulation over the total patient collective, a penalty score revealed dosimetric violations reduced by 79.2%, 75.7% and 93.2% through adaptation.
Treatment plan adaptation was demonstrated to adequately restore relevant dose criteria on a daily basis. While for SAM adaptation approaches dosimetric benefits were realized through ensuring sufficient target coverage, a full re-optimization mainly improved OAR sparing which helps to guide the decision of when to apply which adaptation strategy.
Through the widespread introduction of image-guidance in radiotherapy (IGRT) its delivery accuracy has continuously improved over the last decade . Nevertheless, the acquired image information is still rarely being used to adapt the treatment sequence to the given patient morphology of the day. The standard image-guided treatment workflow for fractionated radiotherapy at a conventional medical linear accelerator (linac) may not adequately account for interfractional organ variations. Methods of adaptive radiotherapy (ART) can address this issue , for example by using daily cone-beam computed tomography (CBCT) images to adapt the initial treatment plan according to daily organ deformations, tissue motion or weight loss [3,4,5]. In this context, inferior CBCT image quality poses a major challenge for the clinical implementation of ART procedures at a conventional linac [6, 7]. Nowadays, this can mostly be overcome by deep learning algorithms which are trained a priori with specific multimodal image datasets (treatment planning CT (pCT) and CBCT) and are afterwards used to generate corrected synthetic CT (sCT) images. These sCT were proven to successfully allow for reasonably fast workflow steps like automated image segmentation and accurate dose calculation for different tumor sites [8,9,10,11].
Prostate cancer patients could benefit from ART-techniques since the dosimetric accuracy is compromised when applying one treatment plan over the entire treatment course [12, 13]. While it would be desirable to consider also intrafractional organ motion, ART can at least compensate for interfractional changes such as varying organ fillings and relative distances between the clinical target volume (CTV) and the organs at risk (OAR) which are predominant causes for dosimetric deviations [14, 15]. Not only because of the Covid-19 pandemic, ultra-hypofractionated stereotactic body radiotherapy (SBRT) for prostate cancer has been increasingly used recently . Exhibiting promising outcomes in terms of toxicity and survival , clinical trials such as PACE-C , pHART3 , HYPO-RT-PC , hypo-FLAME  or HYPOSTAT  were found to be alternatives to the conventional normo-fractionated or moderately hypo-fractionated treatment schemes [23,24,25].
Generating a fully optimized treatment plan for every treatment fraction is typically too time-consuming after IGRT and while the patient is positioned on the treatment couch. Alternatively, treatment plan libraries have been proposed [26,27,28] which yet need an accurate estimate of possible organ deformations. A fast modification of the initial treatment sequence  or applying segment aperture morphing (SAM) based on the “morphology-of-the-day”, eventually followed by an adaptation of the weights (monitor units) and shapes of the initial segments [30,31,32,33], may already adequately solve the problem.
Aims of this work are, first, to determine the maximum overall benefit of daily treatment plan adaptation for hypofractionated prostate SBRT. Second, dose distributions created by different treatment plan adaptation procedures ranging from SAM to a full optimization are compared to the standard IGRT approach. Lastly, a recommendation under which circumstances a certain approach may be most appropriate based on image segmentation information is expressed.
32 patients with primary prostate cancer of low or intermediate risk were included in this retrospective treatment planning study and analyzed after IRB approval (2018-836R-MA) obtained by the ethics committee II of the University of Heidelberg. The guidelines of the PACE-C trial for prostate SBRT  were chosen exemplary as a framework for this analysis. No further consent to participation in the SBRT study was necessary as all patients were treated with the established normo-fractionated regimen. Additional procedures compared to the standard radiation therapy treatment were not performed. Relevant patient demographics are displayed in Table 1. Patients in prone position or carrying endoprosthesis were excluded from the study. All methods and treatment procedures of the present study were carried out in accordance with regulations and relevant guidelines for imaging, structure delineation, treatment planning and treatment delivery in radiation therapy.
All patients were initially treated with volumetric modulated arc therapy (VMAT) at a linac (VersaHD, Elekta AB, Stockholm, Sweden) with a cumulative dose of 75 Gy delivered in a normo-fractionated treatment scheme . Daily kV-CBCT (XVI 5.0, Elekta AB) were acquired prior to each fraction with 120 kV tube voltage, 132mAs exposure time product, an axial length of 27.7 cm, a gantry velocity of 1 rpm and a slice thickness of 2 mm. After rigid image registration of the CBCT and the pCT in 3 positional and 3 rotational degrees of freedom (in XVI), the CBCT were retrospectively reconstructed in the reference system of the pCT with a slice thickness of 1 mm and imported into the treatment planning system (Monaco 5.17, Elekta AB).
Although being “strongly recommended” by the PACE-C guidelines , no fiducial markers were available for image guidance due to our clinical routine treatment setup. One scanner specific CT-number-to-electron-density calibration curve was applied to all pCT and sCT. According to the PACE-C treatment scheme (40 Gy in 5 fractions) , CBCT of every second treatment fraction were selected up to a total of 5 CBCT, starting with the first treatment fraction. The average time between the pCT and the first fraction CBCT was 4.2 ± 1.7 days. The validity of this approach is based on the assumption that anatomical changes in normo- and hypofractionated prostate RT are indistinguishable. With the help of an individually trained cycle-generative adversarial network (cycle-GAN) based pelvis sCT model in the research software ADMIRE (Elekta AB) all CBCT were converted into sCT, taking on average 30 s per dataset. Specifications of the GAN framework and the pelvis sCT model can be taken from previously published literature .
SBRT treatment planning
An expert physician delineated prostate and seminal vesicles (SV) as well as the OAR bladder, rectum, bowel, bilateral femoral heads and penile bulb on pCT and sCT. According to PACE-C guidelines , the clinical target volume (CTV) was defined as the prostate plus the proximal 1 cm of the SV. The planning target volume (PTV) was generated by adding an isotropic margin of 4 mm to the CTV. Based on the dosimetric criteria of the PACE trial summarized in Table 2, a treatment planning template was generated and applied to the pCT of each patient with minor adjustments. A dual 360° arc VMAT delivery with a nominal beam acceleration potential of 10MV with flattening filter was used for all treatments. Dose calculation was Monte Carlo-based with a statistical uncertainty of 1% and a grid size of 2 mm.
Adaptive treatment planning approaches
Warm start optimization (WSO) is a common approach to efficiently adapt a previously generated treatment plan. By using information of the reference treatment plan such as fluence, beam angles, MLC shapes and/or weights the re-optimization includes a priori information in the optimization process to reduce the required calculation time. One possible realization of WSO is segment aperture morphing (SAM) which was previously established by Ahunbay et al. . SAM performs a morphing of the MLC positions in beam’s eye view of every segment or control point according to the respective change between the reference and actual shape of the target volume. Afterwards, an additional modification of the segment weights (monitor units) of the original treatment plan segments or control points can be performed in order to reach the optimum conditions of a pre-defined goal dose distribution (e.g. reference on the pCT). Further improvements can be gained through a re-optimization of the MLC shapes of the SAM-adjusted treatment plan.
In total, four different treatment planning approaches were retrospectively generated on each sCT:
IGRT approach: In analogy to regular IGRT treatments, the reference treatment plan was copied unaltered to the registered sCT. The isocenter position (ISO) on the pCT was the center-of-mass of the CTV while the ISO position on the sCT was shifted based on the translational and rotational corrections for all dose calculations in all scenarios.
ART1 approach: SAM was performed on the original treatment plan, followed by a segment weights optimization. A limit of maximum iterations of 50, 20 and 5 was used for the segment weight, segment shape and the combined optimization, respectively, with an activated gradient control. According to PACE-C criteria the dose distribution was rescaled to a prescription dose of 40 Gy to cover 95% of the CTV in approaches (II)–(IV).
ART2 approach: Same as ART1 but a WSO of the segment weights as well as the segment shapes was performed.
ART3 approach: No prior treatment sequence information was used. A full optimization with dose constraints identical to the ones on the pCT was performed.
All treatment planning steps were performed on a workstation with a dual core Intel Xeon E5-2687 W v4 3.0 GHz processor, 64 GB RAM and an Nvidia Quadro P6000 GPU. To further evaluate clinical feasibility, the calculation time of every treatment planning procedure was recorded.
Dosimetric and statistical treatment plan analysis
Four different treatment plans on each of the five sCT for each of the 32 patients were generated (in total 4 × 5 × 32 = 640), adding up to a total of 160 treatment plans per IGRT, ART1, ART2 and ART3 approach. Mean dose-volume-histograms (DVH) with point-wise standard deviations (SD) over all patients and treatment fractions for the four treatment planning approaches were generated. Furthermore, the most relevant dose-volume constraints of the PACE guidelines  were compared. A treatment plan quality scoring system, as proposed by the ESTRO QUASIMODO group , was applied to characterize the overall benefit of plan adaptation per dose criterion and per treatment planning approach. The penalty score S was based on the percentage difference between an actual value of a dose-volume parameter M for a given dose distribution and the corresponding PACE-C based dose-volume constraint C (see Table 2). Only violations of dose criteria contributed to S which resulted in an optimal plan having a count of S = 0:
The summation index n refers in the following to one of two different summations: First, a summation over all 160 treatment plans per dose criterion (inter-modality comparison) and second, a summation over all four dose criteria per patient (inter-patient comparison). To test whether the results of a certain treatment adaptation approach were statistically different from the ones of a different treatment adaptation approach paired t-tests between the resulting dose-volume parameters were performed. A p-value of p < 0.05 was considered statistically significant.
With regard to possible consequences of daily organ deformations on the dose distribution, the overlap-volumes between the PTV and the OAR bladder VPTV ∩ bladder as well as rectum VPTV ∩ rectum were determined in (cm3). Pearson’s correlation coefficient was calculated between the overlap-volumes of the bladder and rectum, respectively, and the total penalty score S per treatment plan. Subsequently, in order to identify a threshold for potential OAR overdose by means of overlap-volumes, the following ratio was calculated for the IGRT, ART1, ART2 and ART3 approach: The amount of treatment plans with over-average overlap-volumes and coincident dose criteria V37Gy(bladder) > 9 cc or V36Gy(rectum) > 1.5 cc was divided by the amount of total treatment plans with V37Gy(bladder) > 9 cc or V36Gy(rectum) > 1.5 cc.
Dose distributions of the reference plan (panel 1a) and of the ART3 and IGRT plans (panel 1a and 1b) as well as dose difference maps between the four treatment planning approaches in the sagittal isocenter plane (panels 1d,1e,1f) of one exemplary patient are shown in Fig. 1. In contrast to the pCT, a slightly deformed target structure and additional rectal flatus can be identified on the sCT. For the IGRT approach, normal tissue and OAR received doses higher than the PTV prescription dose of 36.25 Gy. These overdosages were reduced by up to 8 Gy (indicated by dark blue areas) by applying the ART1 or ART2 approach. The ART3 approach showed the highest dose sparing, especially along the PTV-bladder interface and in the cranial region of the rectum.
Figure 2 displays the averaged DVH over 160 treatment plans per adaptation approach including point-wise SD (light colored ribbons equate to one standard deviation) of the bladder, rectum, CTV and PTV. The IGRT approach showed the largest SD-ribbons for all four structures. The SD decreased with an increasing degree of adaptation, especially for the OAR. Compared to the ART1 and ART2 approach, the ART3 approach (panel 2d) had the narrowest SD-ribbons above 30 Gy for both OAR and near 45 Gy for the target volumes.
During the initial treatment plan optimization process, the most relevant (i.e. dose-shaping) dose-volume criteria were identified to be V40Gy(CTV), V36.25 Gy(PTV), D98%(PTV), V37Gy(bladder), V36Gy(rectum). Boxplots for these five dose-volume criteria are shown in Fig. 3 with respect to the reference dose distribution and the four adaptation approaches. With regard to CTV and PTV coverage (panel 3a and 3b), the IGRT approach yielded the largest variations with inter-quartile ranges (IQR) of 5.9% and 5.2% but stayed within the acceptable volume region of > 90%. Only the three re-optimization approaches (ART1-ART3) accomplished dose coverages in the optimum volume range close to 95%. Regarding the criterion D98%(PTV) (panel 3c), the mean value and large parts of the IQR of the IGRT approach were considerably below the goal value of 34.4 Gy while the ART3 approach achieved the highest value with 35.1 Gy with the smallest IQR of all three re-optimization approaches. Doses to OAR were highest for the IGRT approach with a mean of 7.4 cc and 2.0 cc for the V37Gy(bladder) and V36Gy(rectum) (panel 3d and 3e), respectively. Several outliers over 15 cc (bladder) and 3 cc (rectum) were noticeable for all three re-optimization approaches. However, for the bladder, all IQR of the re-optimization approaches were within the mandatory volume region below 8 cc and had mean values of < 6.1 cc. In general, the most protocol violations were found for V36Gy(rectum) > 2 cc where only the ART3 approach achieved a mean V36Gy of < 1.3 cc.
The paired t-tests between the four adaptation approaches for the dose-volume criteria presented in Fig. 3 all yielded significant differences except for the correlation of V36Gy(rectum) between ART1 and ART2 approaches with p = 0.32. Due to the rescaling of D95%(CTV) to the prescription dose of 40 Gy for the adaptation approaches, V40Gy(CTV) was neither included in the statistical tests nor in the calculation of the penalty score S.
The results of the penalty score S summed over all 160 treatment plans per dose criterion as well as summed over all four dose criteria per patient/treatment plan is shown in Table 3. Regarding the first summation per dose criterion, the total penalty scores of the three re-optimization approaches were considerably lower than the related score of the IGRT approach (reduced by 79.2%, 75.7% and 93.2%). For the OAR dose criteria, the ART3 approach achieved the largest improvements compared to the IGRT approach with penalty scores reduced by 94.9% and 95.7% for bladder and rectum constraints, respectively. Averaged over the two PTV/OAR criteria, S was reduced by 86.3%/79.9%, 89.9%/76.0% and 75.0%/95.3% for the ART1, ART2 and ART3 approaches. For all four adaptation approaches, the violation of the criterion V36Gy(rectum) < 2 cc constituted the largest proportion of the total penalty score.
With respect to the summation per patient/treatment plan, the largest mean penalty scores averaged over all five sCT were obtained for the IGRT approach with patient P14, P24 and P29 having Smean = 157.6, Smean = 288.7 and Smean = 210.2. Together with patient P6, they showed the largest SD exceeding 100 and also led to comparably high values of Smean and SD for the ART1 and ART2 approaches. For patient P27 and P29 only a full re-planning (ART3) was able to achieve scores of Smean < 4 and SD < 4 whereas P4, 14 and P24 remained as the patients with most dosimetric violations having Smean of at least 19.9 for the ART3 approach. Averaged over the entire patient collective, the mean penalty scores for the IGRT, ART1, ART2 and ART3 approaches were 57.1, 11.9, 13.9 and 3.9.
The overlap-volumes VPTV ∩ bladder and VPTV ∩ rectum were 7.2 ± 3.4 cc and 2.6 ± 1.3 cc averaged over all sCT. Pearson’s correlation coefficients between each of the two overlap-volumes and the total penalty score S per treatment plan for the planning approaches IGRT, ART1, ART2 and ART3 were 0.53, 0.49, 0.50, 0.46 (bladder) and 0.47, 0.52, 0.55, 0.43 (rectum), respectively.
For the IGRT approach, 68.5% and 69.1% of all treatment plans exhibiting high OAR doses of V37Gy(bladder) > 9 cc and V36Gy(rectum) > 1.5 cc coincided with an over-average overlap-volume for bladder and rectum, respectively. This proportion of correctly anticipated overdosage to bladder/rectum amounted to 100%/86.9%, 100%/90.1% and 100%/96.3% for the ART1, ART2 and ART3 approach.
Treatment plan generation times (dose calculation and segment adaptation, where applicable) below three minutes were obtained for the IGRT and ART1 approach with 2.5 ± 0.3 min and 2.6 ± 0.3 min. Due to the additional segment adaptation, the ART2 and ART3 scenarios required 12.1 ± 3.2 min and 19.4 ± 4.0 min to generate a treatment plan on the daily sCT. Since time is of crucial importance for daily ART, plan calculation was also performed with a grid size of 3 mm which led to total plan generation times of 0.8 ± 0.1 min, 0.9 ± 0.1 min, 4.8 ± 1.9 min and 6.7 ± 1.7 min for the IGRT, ART1, ART2 and ART3 approaches, respectively.
Previously published comparative studies on radiotherapy for prostate cancer with and without image guidance concluded the existence of dosimetric advantages of IGRT over non-IGRT [1, 37]. In this sense, the goal of this work was to evaluate whether further dosimetric benefits of different treatment plan adaptation approaches exist in comparison to the conventional IGRT approach. Although this study focused on the specific guidelines of the PACE-C trial, the presented outcomes emphasize the necessity of daily adapting the initial treatment plan and could be applied to other prostate cancer treatment schemes at conventional (non-MRI-based) linacs in an analogous manner.
For the IGRT approach large dose outliers as well as large overall penalty scores and comparably high SD in the DVH were obtained. This implies that the large variation of interfractional organ fillings and deformations cannot be compensated entirely by simply rotating and translating the patient as it has been previously reported [38, 39]. However, the presented IGRT approach did not achieve the precision as obtained by a workflow including fiducial markers. Furthermore, due to the more robust matching process, the time for image registration and manual delineation could have been reduced. In general, two reasons for treatment plan adaptation exist : Adaptation of the treatment sequence to positional changes of the anatomy and/or to a deformed shape of the anatomy. In the first scenario, ART is feasible but not required and can be replaced by a marker-based IGRT which corrects for any translational and rotational degrees of freedom . If, on the other hand, the anatomy at the time of treatment is noteworthy deformed, for example due to rectal or bladder filling, no translation or rotation will be sufficient to adequately compensate for these deficits and ART is a possible treatment approach.
But even in the presence of anatomical deformations at the time of image guidance, alternative methods to compensate for potential dosimetric deviations in relevant areas such as rectal spacers , rectal balloons  or strict bladder filling  exist which present an alternative to ART. While ART is typically more comfortable for the patient and, once the relevant technology is robustly available, is potentially easier to realize it might only be efficient for treatments with a low number of fractions or adaptation procedures. While combining motion mitigation techniques with SAM-based adaptation approaches like ART1 or ART2 could lead to dosimetric improvements and faster optimization times, a combination with a full re-optimization (ART3) may not generate further clinically relevant dosimetric benefits.
The evaluation of the proposed adaptation approaches revealed that an additional benefit over the IGRT dose distribution can already be achieved by modifying the weights of the original plan segments (ART1). This leads to a short treatment plan adaptation time of below three minutes and, due to the insignificant difference of the rectum dose-volume criterion compared to the ART2 approach, makes the ART1 approach the most favorable method of choice for an efficient use in daily clinical routine. An additional optimization of the initial segment shapes (ART2) could lead to improved target coverage and reduced hotspots, on the expenses of increasing the calculation time by a factor of 4.7 (5.3 for the 3 mm grid). These values can only be compared with some tolerance range to related research due to different plan delivery modes, hardware components, calculation settings and image resolutions. Plan generation times of other prostate and lymph node studies were reported as approximately 3 min (SAM plus SWO and dose calculation) , < 2 min (adaptive sequencer plus SWO including different shift methods)  and 11-119 s (different SAM methods at the MR linac with a 3 mm grid)  which compares well to our findings of 2.6 min for a grid of 2 mm for the ART1 approach (SAM plus WSO).
The ART3 approach yielded the highest treatment plan quality and lowest total penalty score since a full re-optimization can adapt best to the daily anatomy of the patient. A recent study also concluded that the most promising adaptive approaches for prostate treatments are re-planning procedures rather than MLC adjustments or pre-optimized plan libraries .
An explanation for the lowest obtained OAR doses of the ART3 approach can be found in the way the SBRT treatment planning template was built. It prioritized a sparing of the OARs over maximum target coverage, while ART1 and ART2 approaches aimed at a reconstruction of the initial target coverage . Thus, a full re-optimization is only recommended for daily situations in which target structure deformations require large re-shaping of the MLC or OAR unexpectedly move into high dose areas. This outcome is in-line with findings of previous studies [45, 46] and, of course, dependent on the applied PTV margin .
The distribution of the mean penalty score S and the respective SD among the total population revealed The ART1 approach was already able to reduce the total Smean of the IGRT approach from 57.1 to 11.9, including a reduction of the SD from 87.4 to 25.7. Most substantial improvements were concentrated rather on treatment plans of single patients (P14, P24 and P29) than scattered throughout the total treatment plan collective. Even the ART3 approach exhibited a few extraordinary criteria violations with a large SD > 25 for patients P4, P14 and P24 which is related to their comparably exceptional anatomy. All three patients showed over-average to very large mean VPTV ∩ bladder and VPTV ∩ rectum of (8.7 cc, 13.9 cc, 15.5 cc) and (4.9 cc, 3.0 cc, 4.1 cc), respectively. Thus, only a re-optimized treatment plan with altered dose constraints or an exclusion from the SBRT treatment would have solved this issue. The largest benefit of ART3 over the adaptive approaches of ART1 and ART2 for all five sCT was noticeable for patient P24. In this case, a high mean VPTV ∩ rectum of 4.1 cc and comparable large daily CTV deformations in anterior direction of up to 0.60 cm due to rectal filling could not be accounted for by both SAM-based approaches ART1 and ART2. Despite its advantage over simpler pass/fail scoring systems, the presented penalty scoring has limitations by being more time consuming to calculate and prioritizing certain dose criteria like OAR violations over PTV violations. Although the penalty score S cannot replace the visual inspection of the spatial dose distribution, it offers a quick treatment plan evaluation at one glance. Noteworthy deviations from the reference score on the pCT can serve as an action trigger to re-evaluate the adapted plan and immediately detect exceptional dosimetric violations.
The obtained Pearson coefficients indicated an intermediate correlation between the overlap-volumes of the OAR and the treatment plan quality measured by the penalty score S. Within this patient collective, accumulated violations of relevant dosimetric criteria per single treatment plan showed a comparable correlation of large overlaps of both OAR with the PTV (maximum difference of 0.06 Pearson’s correlation between bladder and rectum).
In addition, more than two thirds (IGRT) and at least 86% (re-optimization approaches) of treatment plans with OAR overdoses could be predicted by VPTV ∩ bladder and VPTV ∩ rectum. Consequently, overlap-volumes could serve as a fast and reliable indicator for triggering treatment plan adaptation. As recently reported, ad-hoc offline re-planning still constitutes the major part of all performed ART techniques for prostate cancer . Pre-defined action levels based on simple structure metrics and different WSO approaches thus could make a difference in the transition process from offline to online ART treatments.
While deformable image registration (DIR)-based sCT generation includes a potential deformation of the CBCT anatomy based on the matching to the reference pCT , the characteristic feature of the presented neural network-based sCT generation is that the geometry remains identical to the CBCT and solely Hounsfield units or, respectively, electron densities are changed by the sCT algorithm. Different outcomes of image registration between sCT or CBCT and the pCT are negligibly small since a preliminary comparison of all treatment plans between the IGRT approach and the approach having an ISO defined as the center-of-mass of the target volume on the sCT merely changed the mean D98%(PTV) parameter from 32.9 to 33.3 Gy. For DIR-based sCT generation the sCT should be generated prior to image registration to minimize errors.
With regard to clinical feasibility in an online treatment workflow, the lack of tracking devices such as markers/transponder beacons , ultrasound  or EPID megavoltage imaging  constitute a shortcoming of the presented study since they are crucial for the monitoring of time-dependent intrafractional organ motion during the daily treatment plan delivery. All steps of an online adaptive workflow including image processing, segmentation and plan adaptation have to be kept as short as possible in order to avoid losing the benefits generated through this adaptation. Thus, the presented adaptive re-planning methods need to be combined with intrafractional monitoring devices or, for example, an acquisition of a second CBCT prior to the delivery of the adapted treatment plan for an eventual clinical implementation [35, 53, 54].
Despite the study’s retrospective setup, the presented adaptive planning approaches can ensure accurate daily dose delivery and contribute to the general motivation of performing daily online ART for prostate cancer radiotherapy . Several offline and hybrid approaches already proved to be efficient for compensating organ movement but focused on a normo-fractionated treatment scheme with larger PTV margins, more than five fractions or solely on maintaining target coverage [3, 15, 56]. This study revealed the deficiencies of an IGRT approach and, more importantly, offered several adaptation techniques focusing on available planning time (ART1) and increasing target coverage (ART2) or reducing OAR dosage (ART3). Although there has been a long-time debate about the actual benefit of adaptive strategies with respect to their costs [57,58,59], promising results for daily treatment modifications have been obtained by fast and automated tools for image correction [10, 60], segmentation [61, 62] and treatment planning . As previously reported, converting CBCT into sCT and subsequent image segmentation including manual correction of the generated structures takes up to 30 s, 30 s and 5.2 ± 1.6 min . Together with the presented re-optimization approaches an end-to-end adaptive workflow can eventually become feasible within a reasonable timeframe of minimum 0.5 min + 0.5 min + 5.2 min + 2.6 min (ART1 approach) = 8.8 min, being in line with reported adaptation times of 10 min including online plan QA .
This study demonstrated that the three adaptation approaches were able to restore and achieve the dosimetric goals of a prostate SBRT protocol and thus substantially improved treatment plan penalty in comparison to the conventional IGRT approach. Besides standard dose-volume metrics, a penalty score and overlap-volumes could identify the differences of dosimetric benefits among three different adaptation approaches, facilitating the decision when to apply which adaptation strategy.
Availability of data and materials
The datasets analyzed in the presented study are available from the corresponding author on reasonable request.
Clinical target volume
Cycle-generative adversarial network
Deformable image registration
Dose volume histogram
Electronic portal imaging device
Organs at risk
Prostate advances in comparative evidence
Planning target volume
Segment aperture morphing
Stereotactic body radiotherapy
Volumetric modulated arc therapy
Warm start optimization
X-ray volumetric imaging
Bujold A, Craig T, Jaffray D, Dawson LA. Image-guided radiotherapy: has it influenced patient outcomes? Semin Radiat Oncol. 2012;22(1):50–61.
Thornqvist S, Hysing LB, Tuomikoski L, Vestergaard A, Tanderup K, Muren LP, et al. Adaptive radiotherapy strategies for pelvic tumors - a systematic review of clinical implementations. Acta Oncol. 2016;55(8):943–58.
Qin A, Sun Y, Liang J, Yan D. Evaluation of online/offline image guidance/adaptation approaches for prostate cancer radiation therapy. Int J Radiat Oncol Biol Phys. 2015;91(5):1026–33.
Stankiewicz M, Li W, Rosewall T, Tadic T, Dickie C, Velec M. Patterns of practice of adaptive re-planning for anatomic variances during cone-beam CT guided radiotherapy. Tech Innov Patient Support Radiat Oncol. 2019;12:50–5.
Elsayad K, Kriz J, Reinartz G, Scobioala S, Ernst I, Haverkamp U, et al. Cone-beam CT-guided radiotherapy in the management of lung cancer: Diagnostic and therapeutic value. Strahlenther Onkol. 2016;192(2):83–91.
Schulze R, Heil U, Gross D, Bruellmann DD, Dranischnikow E, Schwanecke U, et al. Artefacts in CBCT: a review. Dentomaxillofac Radiol. 2011;40(5):265–73.
Poludniowski GG, Evans PM, Webb S. Cone beam computed tomography number errors and consequences for radiotherapy planning: an investigation of correction methods. Int J Radiat Oncol Biol Phys. 2012;84(1):e109–14.
Maspero M, Houweling AC, Savenije MHF, van Heijst TCF, Verhoeff JJC, Kotte ANTJ, et al. A single neural network for cone-beam computed tomography-based radiotherapy of head-and-neck, lung and breast cancer. Phys Imaging Radiat Oncol. 2020;14:24–31.
Barateau A, De Crevoisier R, Largent A, Mylona E, Perichon N, Castelli J, et al. Comparison of CBCT-based dose calculation methods in head and neck cancer radiotherapy: from Hounsfield unit to density calibration curve to deep learning. Med Phys. 2020.
Eckl M, Hoppen L, Sarria GR, Boda-Heggemann J, Simeonova-Chergou A, Steil V, et al. Evaluation of a cycle-generative adversarial network-based cone-beam CT to synthetic CT conversion algorithm for adaptive radiation therapy. Physica Med. 2020;80:308–16.
Spadea MF, Maspero M, Zaffino P, Seco J. Deep learning-based synthetic-CT generation in radiotherapy and PET: a review. arXiv preprint arXiv:210202734. 2021.
Scaife JE, Thomas SJ, Harrison K, Romanchikova M, Sutcliffe MP, Forman JR, et al. Accumulated dose to the rectum, measured using dose-volume histograms and dose-surface maps, is different from planned dose in all patients treated with radiotherapy for prostate cancer. Br J Radiol. 2015;88(1054):20150243.
Bostel T, Sachpazidis I, Splinter M, Bougatf N, Fechter T, Zamboglou C, et al. Dosimetric impact of interfractional variations in prostate cancer radiotherapy-implications for imaging frequency and treatment adaptation. Front Oncol. 2019;9:940.
Jmour O, Benna M, Champagnol P, Ben Mrad M, Hamrouni A, Obeid L, et al. CBCT evaluation of inter- and intra-fraction motions during prostate stereotactic body radiotherapy: a technical note. Radiat Oncol. 2020;15(1):85.
Bockelmann F, Putz F, Kallis K, Lettmaier S, Fietkau R, Bert C. Adaptive radiotherapy and the dosimetric impact of inter- and intrafractional motion on the planning target volume for prostate cancer patients. Strahlenther Onkol. 2020.
Zaorsky NG, Yu JB, McBride SM, Dess RT, Jackson WC, Mahal BA, et al. Prostate cancer radiation therapy recommendations in response to COVID-19. Adv Radiat Oncol. 2020;5(4):659–65.
Jackson WC, Silva J, Hartman HE, Dess RT, Kishan AU, Beeler WH, et al. Stereotactic body radiation therapy for localized prostate cancer: a systematic review and meta-analysis of over 6000 patients treated on prospective studies. Int J Radiat Oncol Biol Phys. 2019;104(4):778–89.
Brand DH, Tree AC, Ostler P, van der Voet H, Loblaw A, Chu W, et al. Intensity-modulated fractionated radiotherapy versus stereotactic body radiotherapy for prostate cancer (PACE-B): acute toxicity findings from an international, randomised, open-label, phase 3, non-inferiority trial. Lancet Oncol. 2019;20(11):1531–43.
Loblaw A, Cheung P, D’Alimonte L, Deabreu A, Mamedov A, Zhang L, et al. Prostate stereotactic ablative body radiotherapy using a standard linear accelerator: toxicity, biochemical, and pathological outcomes. Radiother Oncol. 2013;107(2):153–8.
Widmark A, Gunnlaugsson A, Beckman L, Thellenberg-Karlsson C, Hoyer M, Lagerlund M, et al. Ultra-hypofractionated versus conventionally fractionated radiotherapy for prostate cancer: 5-year outcomes of the HYPO-RT-PC randomised, non-inferiority, phase 3 trial. Lancet. 2019;394(10196):385–95.
Draulans C, van der Heide UA, Haustermans K, Pos FJ, van der Voort van Zyp J, De Boer H, et al. Primary endpoint analysis of the multicentre phase II hypo-FLAME trial for intermediate and high risk prostate cancer. Radiother Oncol. 2020;147:92–8.
Jiang P, Krockenberger K, Vonthein R, Tereszczuk J, Schreiber A, Liebau S, et al. Hypo-fractionated SBRT for localized prostate cancer: a German bi-center single treatment group feasibility trial. Radiat Oncol. 2017;12(1):138.
Lehrer EJ, Kishan AU, Yu JB, Trifiletti DM, Showalter TN, Ellis R, et al. Ultrahypofractionated versus hypofractionated and conventionally fractionated radiation therapy for localized prostate cancer: a systematic review and meta-analysis of phase III randomized trials. Radiother Oncol. 2020;148:235–42.
Franzese C, Badalamenti M, Di Brina L, D’Agostino G, Franceschini D, Comito T, et al. Linac-based stereotactic body radiation therapy for low and intermediate-risk prostate cancer. Strahlenther Onkol. 2020;196(7):608–16.
Hocht S, Aebersold DM, Albrecht C, Bohmer D, Flentje M, Ganswindt U, et al. Hypofractionated radiotherapy for localized prostate cancer. Strahlenther Onkol. 2017;193(1):1–12.
Heijkoop ST, Langerak TR, Quint S, Bondar L, Mens JW, Heijmen BJ, et al. Clinical implementation of an online adaptive plan-of-the-day protocol for nonrigid motion management in locally advanced cervical cancer IMRT. Int J Radiat Oncol Biol Phys. 2014;90(3):673–9.
Antico M, Prinsen P, Cellini F, Fracassi A, Isola AA, Cobben D, et al. Real-time adaptive planning method for radiotherapy treatment delivery for prostate cancer patients, based on a library of plans accounting for possible anatomy configuration changes. PLoS One. 2019;14(2):e0213002.
de Jong R, Crama KF, Visser J, van Wieringen N, Wiersma J, Geijsen ED, et al. Online adaptive radiotherapy compared to plan selection for rectal cancer: quantifying the benefit. Radiat Oncol. 2020;15(1):162.
Crijns W, Van Herck H, Defraene G, Van den Bergh L, Slagmolen P, Haustermans K, et al. Dosimetric adaptive IMRT driven by fiducial points. Med Phys. 2014;41(6):061716.
Winkel D, Bol GH, Werensteijn-Honingh AM, Kiekebosch IH, van Asselen B, Intven MPW, et al. Evaluation of plan adaptation strategies for stereotactic radiotherapy of lymph node oligometastases using online magnetic resonance image guidance. Phys Imaging Radiat Oncol. 2019;9:58–64.
Ahunbay EE, Peng C, Holmes S, Godley A, Lawton C, Li XA. Online adaptive replanning method for prostate radiotherapy. Int J Radiat Oncol Biol Phys. 2010;77(5):1561–72.
Ahunbay EE, Peng C, Chen GP, Narayanan S, Yu C, Lawton C, et al. An on-line replanning scheme for interfractional variations. Med Phys. 2008;35(8):3607–15.
van Timmeren JE, Chamberlain M, Krayenbuehl J, Wilke L, Ehrbar S, Bogowicz M, et al. Treatment plan quality during online adaptive re-planning. Radiat Oncol. 2020;15(1):203.
PACE: International randomised study of prostatectomy vs stereotactic body radiotherapy (SBRT) and conventional radiotherapy vs SBRT for organ-confined prostate cancer. Accessed August 30 2020. Available from: https://www.icr.ac.uk/our-research/centres-and-collaborations/centres-at-the-icr/clinical-trials-and-statistics-unit/clinical-trials/pace
Sihono DSK, Ehmann M, Heitmann S, von Swietochowski S, Grimm M, Boda-Heggemann J, et al. Determination of intrafraction prostate motion during external beam radiation therapy with a transperineal 4-dimensional ultrasound real-time tracking system. Int J Radiat Oncol Biol Phys. 2018;101(1):136–43.
Bohsung J, Gillis S, Arrans R, Bakai A, De Wagter C, Knoos T, et al. IMRT treatment planning: a comparative inter-system and inter-centre planning exercise of the ESTRO QUASIMODO group. Radiother Oncol. 2005;76(3):354–61.
Wortel RC, Incrocci L, Pos FJ, Lebesque JV, Witte MG, van der Heide UA, et al. Acute toxicity after image-guided intensity modulated radiation therapy compared to 3D conformal radiation therapy in prostate cancer patients. Int J Radiat Oncol Biol Phys. 2015;91(4):737–44.
Godley A, Ahunbay E, Peng C, Li XA. Accumulating daily-varied dose distributions of prostate radiation therapy with soft-tissue-based kV CT guidance. J Appl Clin Med Phys. 2012;13(3):3859.
Peng C, Ahunbay E, Chen G, Anderson S, Lawton C, Li XA. Characterizing interfraction variations and their dosimetric effects in prostate cancer radiotherapy. Int J Radiat Oncol Biol Phys. 2011;79(3):909–14.
Winkel D, Bol GH, Kroon PS, van Asselen B, Hackett SS, Werensteijn-Honingh AM, et al. Adaptive radiotherapy: the Elekta unity MR-linac concept. Clin Transl Radiat Oncol. 2019;18:54–9.
Deutschmann H, Kametriser G, Steininger P, Scherer P, Scholler H, Gaisberger C, et al. First clinical release of an online, adaptive, aperture-based image-guided radiotherapy strategy in intensity-modulated radiotherapy to correct for inter- and intrafractional rotations of the prostate. Int J Radiat Oncol Biol Phys. 2012;83(5):1624–32.
Greco C, Stroom J, Vieira S, Mateus D, Cardoso MJ, Soares A, et al. Reproducibility and accuracy of a target motion mitigation technique for dose-escalated prostate stereotactic body radiotherapy. Radiother Oncol. 2021;160:240–9.
Yoon S, Cao M, Aghdam N, Shabsovich D, Kahlon S, Ballas L, et al. Prostate bed and organ-at-risk deformation: prospective volumetric and dosimetric data from a phase II trial of stereotactic body radiotherapy after radical prostatectomy. Radiother Oncol. 2020;148:44–50.
Kontaxis C, Bol GH, Kerkmeijer LGW, Lagendijk JJW, Raaymakers BW. Fast online replanning for interfraction rotation correction in prostate radiotherapy. Med Phys. 2017;44(10):5034–42.
Thongphiew D, Wu QJ, Lee WR, Chankong V, Yoo S, McMahon R, et al. Comparison of online IGRT techniques for prostate IMRT treatment: adaptive vs repositioning correction. Med Phys. 2009;36(5):1651–62.
Wahl M, Descovich M, Shugard E, Pinnaduwage D, Sudhyadhom A, Chang A, et al. Interfraction anatomical variability can lead to significantly increased rectal dose for patients undergoing stereotactic body radiotherapy for prostate cancer. Technol Cancer Res Treat. 2017;16(2):178–87.
Battista JJ, Johnson C, Turnbull D, Kempe J, Bzdusek K, Van Dyk J, et al. Dosimetric and radiobiological consequences of computed tomography-guided adaptive strategies for intensity modulated radiation therapy of the prostate. Int J Radiat Oncol Biol Phys. 2013;87(5):874–80.
Bertholet J, Anastasi G, Noble D, Bel A, van Leeuwen R, Roggen T, et al. Patterns of practice for adaptive and real-time radiation therapy (POP-ART RT) part II: Offline and online plan adaption for interfractional changes. Radiother Oncol. 2020.
Kida S, Nakamoto T, Nakano M, Nawa K, Haga A, Kotoku J, et al. Cone Beam Computed Tomography Image Quality Improvement Using a Deep Convolutional Neural Network. Cureus. 2018;10(4):e2548.
Lovelock DM, Messineo AP, Cox BW, Kollmeier MA, Zelefsky MJ. Continuous monitoring and intrafraction target position correction during treatment improves target coverage for patients undergoing SBRT prostate therapy. Int J Radiat Oncol Biol Phys. 2015;91(3):588–94.
Li M, Ballhausen H, Hegemann NS, Reiner M, Tritschler S, Gratzke C, et al. Comparison of prostate positioning guided by three-dimensional transperineal ultrasound and cone beam CT. Strahlenther Onkol. 2017;193(3):221–8.
Azcona JD, Li R, Mok E, Hancock S, Xing L. Automatic prostate tracking and motion assessment in volumetric modulated arc therapy with an electronic portal imaging device. Int J Radiat Oncol Biol Phys. 2013;86(4):762–8.
Vanhanen A, Poulsen P, Kapanen M. Dosimetric effect of intrafraction motion and different localization strategies in prostate SBRT. Phys Med. 2020;75:58–68.
Hewson EA, Nguyen DT, O'Brien R, Kim JH, Montanaro T, Moodie T, et al. The accuracy and precision of the KIM motion monitoring system used in the multi-institutional TROG 15.01 Stereotactic Prostate Ablative Radiotherapy with KIM (SPARK) trial. Med Phys. 2019;46(11):4725–37.
McVicar N, Popescu IA, Heath E. Techniques for adaptive prostate radiotherapy. Phys Med. 2016;32(3):492–8.
Liu H, Wu Q. Dosimetric and geometric evaluation of a hybrid strategy of offline adaptive planning and online image guidance for prostate cancer radiotherapy. Phys Med Biol. 2011;56(15):5045–62.
Li XA, Wu Q, Orton CG. Point/Counterpoint. Online adaptive planning for prostate cancer radiotherapy is necessary and ready now. Med Phys. 2014;41(8):080601.
Fiorino C, Cozzarini C, Passoni P. The promise of adaptive radiotherapy for pelvic tumors: “too high cost for too little result” or “a low cost for a significant result”? Acta Oncol. 2016;55(8):939–42.
Ghilezan M, Yan D, Martinez A. Adaptive radiation therapy for prostate cancer. Semin Radiat Oncol. 2010;20(2):130–7.
Chen L, Liang X, Shen C, Jiang S, Wang J. Synthetic CT generation from CBCT images via deep learning. Med Phys. 2020;47(3):1115–25.
Cardenas CE, Yang J, Anderson BM, Court LE, Brock KB. Advances in auto-segmentation. Semin Radiat Oncol. 2019;29(3):185–97.
Schreier J, Genghi A, Laaksonen H, Morgas T, Haas B. Clinical evaluation of a full-image deep segmentation algorithm for the male pelvis on cone-beam CT and CT. Radiother Oncol. 2020;145:1–6.
Lee H, Kim H, Kwak J, Kim YS, Lee SW, Cho S, et al. Fluence-map generation for prostate intensity-modulated radiotherapy planning using a deep-neural-network. Sci Rep. 2019;9(1):15671.
Zhang J, Ahunbay E, Li XA. Technical Note: Acceleration of online adaptive replanning with automation and parallel operations. Med Phys. 2018;45(10):4370–6.
The authors would like to thank Gustav Meedt, Peter Voet and Virgil Willcut from Elekta AB for their support on the adaptation approaches. The research software ADMIRE and an implementation of the treatment adaptation approaches were provided by Elekta in a research version of Monaco.
Open Access funding enabled and organized by Projekt DEAL.
Ethics approval and consent to participate
IRB approval (2018-836R-MA) was obtained by the ethics committee II of the University of Heidelberg. All involved patients gave their informed consent to the analysis of their medical datasets and personal data.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Eckl, M., Sarria, G.R., Springer, S. et al. Dosimetric benefits of daily treatment plan adaptation for prostate cancer stereotactic body radiotherapy. Radiat Oncol 16, 145 (2021). https://doi.org/10.1186/s13014-021-01872-9