Skip to main content
Download PDF

Dosimetric benefits of adaptive radiation therapy for patients with stage III non-small cell lung cancer

Radiation Oncology volume 18, Article number: 34 (2023) Cite this article

Abstract

Background

Daily adaptive radiation therapy (ART) of patients with non-small cell lung cancer (NSCLC) lowers organs at risk exposure while maintaining the planning target volume (PTV) coverage. Thus, ART allows an isotoxic approach with increased doses to the PTV that could improve local tumor control. Herein we evaluate daily online ART strategies regarding their impact on relevant dose-volume metrics.

Methods

Daily cone-beam CTs (1 × n = 28, 1 × n = 29, 11 × n = 30) of 13 stage III NSCLC patients were converted into synthetic CTs (sCTs). Treatment plans (TPs) were created retrospectively on the first-fraction sCTs (sCT1) and subsequently transferred unaltered to the sCTs of the remaining fractions of each patient (sCT2−n) (IGRT scenario). Two additional TPs were generated on sCT2−n: one minimizing the lung-dose while preserving the D95%(PTV) (isoeffective scenario), the other escalating the D95%(PTV) with a constant V20Gy(lungipsilateral) (isotoxic scenario).

Results

Compared to the original TPs predicted dose, the median D95%(PTV) in the IGRT scenario decreased by 1.6 Gy ± 4.2 Gy while the V20Gy(lungipsilateral) increased in median by 1.1% ± 4.4%. The isoeffective scenario preserved the PTV coverage and reduced the median V20Gy(lungipsilateral) by 3.1% ± 3.6%. Furthermore, the median V5%(heart) decreased by 2.9% ± 6.4%. With an isotoxic prescription, a median dose-escalation to the gross target volume of 10.0 Gy ± 8.1 Gy without increasing the V20Gy(lungipsilateral) and V5%(heart) was feasible.

Conclusions

We demonstrated that even without reducing safety margins, ART can reduce lung-doses, while still reaching adequate target coverage or escalate target doses without increasing ipsilateral lung exposure. Clinical benefits by means of toxicity and local control of both strategies should be evaluated in prospective clinical trials.

Background

The aim of radiation therapy (RT) is to achieve a high tumor control probability (TCP) with the lowest possible normal tissue complication probability (NTCP) [1]. For patients with locally advanced non-small cell lung cancer (NSCLC) the dose limiting organs at risk (OAR) are usually lungs, heart and esophagus, which tend to develop pneumonitis [2, 3], cardiovascular toxicity [4,5,6,7], or acute esophagitis [8, 9] due to high radiation doses, respectively. Most NSCLC radiation therapies are performed with 60–70 Gy total dose in 1.8–2 Gy-fractions [5, 10, 11]. Survival rates of patients with locally advanced (stage III) NSCLC have improved considerably during the last years. Modern concurrent chemoradiotherapy (cCRT) results in 5-year overall survival (OS) rates of 32–33% [12, 13]; results can be further improved in patients who are eligible for immunotherapy (absence of toxic events ≥ grade 2) with durvalumab as a consolidation therapy resulting in 5-year OS rates of approximately 43% [13]. Despite these improvements, only one third of patients treated with cCRT or immunotherapy survives or remains free from progression five years after treatment [13].

Currently, treatment plans are typically generated on a pre-treatment computed tomography (CT) scan and delivered in all treatment fractions by means of daily image guided radiation therapy (IGRT), sometimes in combination with an active motion management concept. Intrafractional motion management is mostly realized by either an internal target volume (ITV) margin, beam tracking or gating concepts [14]. The latter two reduce the planning target volume (PTV) size and aim at a quasi-static situation. However, all these workflows ensure adequate patient positioning but neglect interfractional morphologic changes in the patient’s anatomy which may lead to discrepancies between the delivered and planned dose [1, 15, 16]. This is of particular relevance for NSCLC since these tumors are generally early responders to radiation [16, 17]. On average the gross target volume (GTV) of NSCLC patients shrinks 0.6–2.4% per day during RT [3].

In online adaptive radiation therapy (ART), treatment plans are adjusted to a possibly changed patient anatomy [1, 15, 16] on daily pre-treatment imaging. This is mostly performed with cone-beam computed tomography (CBCT)-based approaches. Compared to fan-beam CTs, CBTCs present limitations such as, inferior image quality and lack of a unique CT-number-to-electron density calibration [18,19,20,21]. Furthermore, they have only a limited field of view (FOV) capacity, which may not enclose the structure set in its entirety. For image segmentation and dose calculation, the insufficient image quality has recently been overcome by generating synthetic CTs (sCTs) with an artificial intelligence (AI) approach [22].

The aim of ART is to minimize the discrepancies between high dose volumes and the actual clinical target volumes (CTV) and therefore, to further broaden the therapeutic window. In ART for patients with NSCLC, this can be realized with one of the following adaptation strategies: sparing the OARs while leaving the dose delivered to the PTV unaltered (isoeffective scenario) or escalating the dose to the PTV without increasing the ipsilateral lung-dose (isotoxic scenario) [16]. By minimizing OAR-doses, the risk of toxic events can be reduced [3, 9, 23]. This increases the probability that patients will be eligible for immunotherapy [7, 24]. Several analyses showed that increased OAR-doses are associated with worse OS [4, 5, 25]. Moreover, no correlation between effect and different tumor types or locations has been yet determined. Several studies concluded that in lung cancer higher fractional doses lead to higher TCP and OS [26,27,28]. This was, however, questioned by the results of the RTOG-0617 study which showed that patients treated with 60 Gy had better progression-free survival (PFS) and OS rates than patients treated with 74 Gy [5]. In the 74 Gy arm higher doses to the heart, lung, and esophagus were delivered compared to the 60 Gy arm [5, 25]. Higher lung- and heart-doses were associated with worsened OS [4, 5, 25, 29], which may explain the better outcome of the 60-Gy arm. Further analysis of the patient population indicated that dose-escalation may improve OS rates for patients with radioresistant genotypes [30].

This study aims to determine the differences between planned and delivered doses. Furthermore, the maximum dosimetric benefits of isotoxic and isoeffective daily online ART approaches with deep inspiration breath-hold technique (DIBH) for target immobilization shall be determined.

Methods

Patient population and treatment workflow

In this retrospective treatment planning study 13 patients with stages III/IV NSCLC, previously treated with cCRT were identified (Table 1). The primary tumors of two patients with distant metastases were included in the study, while the secondary lesions were not.

Table 1 Patient characteristics
Full size table

The mean treatment period was 47.5 days ± 4.0 days (range: 40–54 days). All patients were treated on a linear accelerator (VersaHD, Elekta AB, Stockholm, Sweden) with 10MV volumetric modulated arc therapy (VMAT). Treatment planning CTs (pCTs) and all dose deliveries were performed using a computer-controlled DIBH for target immobilization (ABC (Active Breathing Coordinator, Elekta AB, Sweden): 6, Catalyst (Catalyst, C-RAD, Sweden): 7). Therefore we did not consider an ITV but define GTV-PTV margins (axial: 10 mm, inferior-superior: 15 mm) [31]. For each patient, daily kV-based CBCT (XVI 5.0, Elekta AB, Sweden) scans were acquired for patient positioning. These scans were also obtained in multiple breath-hold-phases (“stop-and-go” breath-hold-only approach), immediately prior every treatment.

Study image data preparation

The CBCT scans were reconstructed, rigidly registered to the pCTs using the clinically used registrations, in which a tumor match was favored and exported with a slice thickness of 3 mm to the treatment planning system (Monaco 5.11, Elekta). These CBCTs were converted into sCTs in a dedicated research software (ADMIRE, Elekta). For sCT generation this software utilizes a pre-trained artificial neural network [22]. Retrospectively, an expert physician contoured the GTV, heart, lungs, esophagus, spinal cord (sc), and the patient outline on the sCTs of the first fractions (sCT1) for all patients. These structures were deformably registered in an unsorted manner from the sCT1 of a patient to the remaining sCTs (sCT2−n) of the patient using deformable image registration (DIR). The deformed structures were retrospectively reviewed and corrected for all sCTs of every patient, if necessary, by the same expert physician. The same window and level settings were used for correcting all structures in all sCTs of the same patient, respectively. An additional structure ipsilateral lung minus GTV (lungipsilateral) was created for all sCTs. According to the clinical protocol all PTV2−n were generated out of the GTV2−n with GTV-PTV margins of 10 mm in axial and 15 mm in inferior-superior directions. These margins were used unaltered throughout the study and are in agreement with the margins which were used in the RTOG-0617 study protocol [5].

Treatment planning and dosimetric analysis

Initial VMAT treatment plans were generated in a Monte Carlo based treatment planning system (Monaco 5.11, Elekta) with a grid spacing of 3 mm and a statistical uncertainty of 1% per dose calculation. Dosimetric constraints for inverse treatment planning were: V20Gy(lungipsilateral) ≤ 37%, V20Gy(lungtotal) ≤ 30%, ipsilateral mean lung dose (MLDipsilateral) ≤ 20 Gy, V5Gy(heart) as low as possible, V35Gy(heart) < 10%, mean heart dose (MHD) < 10 Gy, D0.1%(sc) < 50.5 Gy and the mean esophagus dose (MED) < 34 Gy, while 95% of the PTV should be covered with the prescription dose (Dpresc). Except for the heart the dosimetric constraints were adopted from RTOG-0617 [5]. The tolerance doses of the heart were further reduced, because V5Gy(heart), V35Gy(heart), and MHD are potentially associated with cardiac events [23, 32, 33]. All treatment plans were clinically accepted for patient treatment and optimized for maximum dose conformality.

The initial treatment plan of each patient was recalculated with identical control point settings, grid size, statistical uncertainty, and an isocenter which resulted from the image registration process onto sCT2−n. Assuming only interfractional and neglecting residual intrafractional motion in DIBH [34] this scenario represents the daily delivered dose to the patient with IGRT (IGRT scenario). For the two adaptation approaches, new treatment plans were re-optimized on sCT2−n. Therefore, a new isocenter was set in the center-of-mass of the daily PTVs. For the isoeffective scenario the ipsilateral lung constraints of the initial prescription template were iteratively reduced as long as the target coverage remained adequate. Each treatment plan was normalized to cover 95% of the actual PTV2−n with Dpresc. For the treatment plans of the isotoxic approach an equivalent procedure was applied with increasing the target dose. Subsequently, the V20Gy,2−n(lungipsilateral) was normalized to the initial treatment plan V20Gy,1(lungipsilateral) … with a maximum normalization between 90% and 110%. For three patients, the lungs were cropped due to a reduced FOV in the CBCTs. For these patients, the absolute initial ipsilateral lung volume was taken for normalization. In addition, a maximum dose for the D95% (PTV) of 3.3 Gy per fraction was implemented in accordance with the 2.2–3.8 Gy of RTOG-1106 study protocol. Our lower maximal fractional dose of 3.3 Gy was chosen to preferably stay below the maximum total dose of 80.4 Gy of the RTOG-1106. Their study showed no adverse effects for this dose-escalation [35]. To avoid uncertainties due to a DIR-based dose accumulation, the resulting dose-volume histogram (DVH)-parameters of each fraction were evaluated independently. To enable using the standard DVH-parameters for cumulative dose distributions (like V20Gy) each fraction was evaluated with the total Dpresc. Thus, the fractional dose limit of 3.3 Gy (D95%(PTV)) corresponded to a limit of 100 Gy (D95%(PTV)) in the total Dpresc. In isotoxic cases where D95%(PTV) > 100 Gy was feasible, the treatment plan was normalized, so that 100 Gy covered 95% of the PTV. The mean value of all fractions of each patient was considered as an estimation of the dose accumulation.

The dose distributions of all fractions were evaluated regarding the GTV and PTV coverage, V20Gy(lungipsilateral), MLDipsilateral, MLDcontralateral, MHD, V5Gy(heart), MED, and the D0.1%(sc). Furthermore, the equivalent dose in 2 Gy-fractions (EQD2) was determined using the linear quadratic model [36] with a specific tissue characterization ratio α/β = 8.2 Gy for the GTV and PTV [37]. The dosimetric parameters of the initial treatment plans were compared to the resulting plans of the IGRT, isoeffective, and isotoxic scenarios.

In Fig. 1 the dose distributions in the axial isocenter plane of the initial treatment plan (Fig. 1a) on sCT1 of a representative patient is compared to the IGRT-, isoeffective- and isotoxic-dose distributions on sCT30. One can see the regressed GTV in sCT30 (Fig. 1b–d).

Fig. 1
figure 1

Exemplary dose distributions of a representative patient. Shown are dose distributions (a) of the original treatment plan, (b) of the last treatment fraction without adaptive radiotherapy (ART), (c) of an isoeffective ART, and (d) of an isotoxic ART treatment plan

Full size image

Statistical analysis

The aforementioned dose-volume relationships for IGRT, isoeffective and isotoxic ART approaches were compared to the initial treatment plan results and tested for statistically significant (p < 0.05) differences using a non-parametric Wilcoxon signed-rank test in R (RStudio 1.4.1717, PBC, Boston, MA). Pearson’s correlation coefficients (PCC) were determined to analyze potential correlations between GTV-regression and analyzed DVH-parameters.

Ethics statement

This investigation was performed according to the principles of the Declaration of Helsinki and after Institutional Review Board (IRB) approval (2018-836R-MA). All data was anonymized prior to inclusion.

Results

GTV-regression

The mean decrease in GTV on the final sCTs was 59.9% ± 15.8% (range: 24.0–81.8%) compared with the initial GTV, corresponding to a mean reduction of 1.4% ± 0.6% per treatment fraction and 0.9% ± 0.3% per day (Fig. 2).

Fig. 2
figure 2

Daily volume changes of the gross target volume (GTV). Displayed are GTV dynamics for 13 patients (P1–P13) with locally advanced non-small cell lung cancer

Full size image

Dosimetric analysis

The median DVH-parameters of the initial treatment plans and the three scenarios are presented in Table 2.

Table 2 DVH-parameters
Full size table

While the GTV-coverage remained intact in the IGRT scenario, the median D95%(PTV) was 1.6 Gy ± 4.2 Gy and the median EQD2(D95%(PTV)) 1.7 Gy ± 5.0 Gy lower than initially planned. The D95%(GTV) and D95%(PTV) were below 95% volume coverage in 21 (5.4%) and 108 (27.9%) fractions. Averaged over all treatment fractions of each patient the GTV and PTV coverage of 1 (7.7%) and 5 patients (38.5%) were too low, respectively. In the IGRT scenario all analyzed OAR-parameters were slightly higher than initially predicted. Compared to the initial treatment plan the median V20Gy(lungipsilateral), MLDipsilateral, V5Gy(heart) and MHD were 1.1% ± 4.4%, 0.5 Gy ± 1.8 Gy, 0.2% ± 6.7%, and 0.2 Gy ± 1.2 Gy higher. The V20Gy(lungipsilateral), MLDipsilateral, and the MHD were violating the previously mentioned constraints of 37%, 20 Gy, and 10 Gy in 69 (17.8%), 6 (1.6%), and 41 (10.6%) fractions and on average for 3 (23.1%), 0, and 1 patients (7.7%), while these tolerances were never exceeded in the initial treatment plans.

The PTV-coverage was restored in the isoeffective ART plans. The above-mentioned constraints for the D95%(GTV) and D95%(PTV) were not violated in any fraction in the isoeffective treatment plans. Simultaneously, the ipsilateral lung and heart were spared by 3.1% ± 3.6% (median V20Gy(lungipsilateral)), 1.4 Gy ± 1.3 Gy (median MLDipsilateral), 2.9% ± 6.4% (median V5Gy(heart)), and 0.5 Gy ± 1.4 Gy (median MHD) compared to the initial plan. The OAR-sparing of the isoeffective approach is even more pronounced compared to the IGRT scenario with 3.7% ± 6.9% (median V20Gy(lungipsilateral)), 1.9 Gy ± 2.7 Gy (median MLDipsilateral), 3.0% ± 5.7% (median V5Gy(heart)), and 0.6 Gy ± 1.3 Gy (median MHD). The constraints for the V20Gy(lungipsilateral), MLDipsilateral, and the MHD were exceeded for 7 (1.8%), 0, and 13 (3.4%) fractions and on average for no patient, respectively.

In the dose-escalation scenario median dose-escalations of 10.0 Gy ± 8.1 Gy (D95%(GTV)), 12.4 Gy ± 10.3 Gy (EQD2(D95%(GTV))), 6.6 Gy ± 8.9 Gy (D95%(PTV)), and 8.1 Gy ± 13.7 Gy (EQD2(D95%(PTV))) were achieved. The aforementioned constraints for the D95%(GTV) and D95%(PTV) were violated in 0 and 12 (3.1%) cases, respectively. The constraints for the V20Gy(lungipsilateral), MLDipsilateral, and the MHD were violated in 0, 7 (1.8%), and 27 (7.0%) fractions, respectively. In the accumulated dose no constraints were violated for any patients.

The IGRT scenario increased the median D0.1%(sc), MED, and the MLDcontralateral compared with the initial TPs. In the two ART approaches these median DVH-parameters were lower compared with the IGRT scenario. The previously defined constraints (D0.1%(sc) < 50.5 Gy, MED < 34 Gy, V20Gy(lungtotal) ≤ 30%) were not violated in any fraction of the three scenarios. The ipsilateral lung volume changes resulted in higher MLDipsilateral-doses in some fractions. That explains the higher maximal MLDipsilateral-doses in the IGRT, isoeffective, and isotoxic scenarios. Higher maximum doses for the V5Gy(heart) and MHD in the three scenarios were caused by GTV movements towards the heart in 2 patients.

DVH-parameter correlation with GTV-regression

A PCC between the GTV-volume and a DVH-parameter ≥ ± 0.5, which indicates an at least moderately strong correlation, was found for the V20Gy(lungipsilateral) = −0.6 in the IGRT scenario. This indicates that a decreasing GTV-volume leads to an increase in ipsilateral lung-dose with IGRT-only. For the isoeffective scenario, on the other hand, the PCC of V20Gy(lungipsilateral) was + 0.6 which implies a decrease of ipsilateral lung-dose with a decrease in GTV-volume. In the isotoxic scenario with a PCC of the D95%(GTV) = −0.5, an increasing GTV-dose with decreasing GTV-volume was observed.

Fraction-wise analysis of a representative patient

Fraction-wise dosimetric analysis of the D95%(GTV), D95%(PTV), V20Gy(lungipsilateral), and V5Gy(heart) over all treatment fractions of one representative patient (with median GTV-reduction) with Dpresc = 60 Gy is presented in Fig. 3. In the last sCT of this patient, the GTV-volume was reduced to 59.3% of the initial GTV. The mean D95%(PTV) was 69.8 Gy in the isotoxic scenario, whereas in the isoeffective scenario the mean V20Gy(lungipsilateral) was reduced by 3.0%. One can see the aforementioned dose-escalation (D95%(GTV) and D95%(PTV)) and the lung-dose sparing (V20Gy) in the isotoxic and isoeffective scenario, respectively. For the V5%(heart) a strong fluctuation is noticeable in every scenario. In the IGRT scenario an increase in V20Gy(lungipsilateral) by 1.6% was found while the target coverage remained intact over all fractions.

Fig. 3
figure 3

Fraction-wise analysis of dose-volume histogram (DVH)-parameters of a representative patient. The DVH-parameters of the gross target volume (GTV) (D95%(GTV)), planning target volume (PTV) (D95%(PTV)), ipsilateral lung (V20Gy(lungipsilateral)) and the heart (V5Gy (heart)) are shown as a function of treatment fractions for the three scenarios - without adaptive radiotherapy (ART) (blue), isoeffective ART (green), and isotoxic ART (red). The dashed grey line represents the respective DVH-parameter of the initial treatment plan

Full size image

Discussion

To our knowledge, this is the first study evaluating the dosimetric benefits of daily isotoxic and isoeffective ART approaches for patients with stage III NSCLC on daily CBCT-based sCTs. The generated sCTs have the advantage of preserving the daily anatomy and providing accurate electron density allocations, comparable to those of the pCTs. While dosimetric errors were reported to be typically less than 1% with this method [22], the errors were further minimized by comparing dosimetric results to treatment plans optimized on the sCT1 and not pCTs.

Dose accumulation is mostly based on deformation vector fields (DVFs) generated with a DIR. The accuracy of this method depends strongly on the image resolution, artifacts, and image distortion [38, 39]. Furthermore, the generation of the DVFs is a problem with no unique solution due to too many degrees of freedom [40]. This eventually results in dose inaccuracy in DIR-based dose accumulation [39]. This inaccuracy becomes more pronounced for regressing structures because their mass is not preserved [41, 42]. That makes dose accumulations challenging when shrinking tumors are included [41, 42]. Zhong et al. reported a reduced MLD by around 2 Gy for isoeffective ART of patients with NSCLC from 17.3 Gy to either 15.2 Gy, 14.5 Gy or 14.8 Gy when using three different DIR algorithms for the same data. That yields an algorithm-based MLD fluctuation of 0.7 Gy (4%) between the three DIR algorithms [42].

Without adaptation, the GTV-coverage remained adequate, but the median D95%(PTV) decreased by 1.6 Gy ± 4.2 Gy. In a study by Luo et al. in which the treatment plans of 24 NSCLC patients were retrospectively recalculated on weekly CBCTs, only 81.5% of the PTVs received the Dpresc instead of 95.6% of the initial plans [43]. The evaluated OAR-doses were exceeded for all analyzed DVH-parameters. Since the known OAR tolerances are today solely based on the initial treatment plans [44], OAR tolerances are likely to change for ART. The analysis of Hoegen et al. yielded a lung dose increase (MLD, V20Gy) for 10 patients with NSCLC by daily recalculation on CBCTs without plan adaptations comparable to the results of this work [45].

Benefits of isoeffective ART

In isoeffective ART, treatment plan adaptations to the regressed tumor volume allowed a restoration of the PTV coverage while minimizing the OAR-doses. The median MLDipsilateral and V20Gy (lungipsilateral) were 1.4 Gy ± 1.3 Gy and 3.1% ± 3.6% lower than initially planned. In the study of Hoegen et al. with weekly isoeffective plan adaptations for 10 patients with NSCLC, this effect is less evident, which can be explained by less frequent adaptations [45]. Guckenberger et al. obtained weekly CT for 13 patients with NSCLC to simulate ART. At week 3 or 5, or both, the treatment plans were adjusted. The results of their isoeffective ART at weeks 3 and 5 were comparable to our results for MLD. The effect that they achieved similar results with only 2 adaptations could be due to larger mean GTV-volumes and higher GTV-regression rates in their patient collective [46].

In a recent retrospective study, an isoeffective ART group had a 4.3 Gy lower median MLD and 3.6 Gy lower median MHD than a non-ART group. It is of notice that this group in contrast to our study used smaller margin sizes for the ART group, which explains the large MLD and MHD sparing [24]. Despite their retrospective nature, these results demonstrate the possibilities of isoeffective ART. Only 20%, 7%, and 0.4% of the patients in the ART group experienced pneumonitis with a grade ≥ 2, ≥ 3, or lethal instead of 50%, 21%, and 6% of the patient population without ART. The OS and PFS after two years increased by 13% and 8% in the ART group compared to the non-ART control group, respectively [24]. In this study, new 4D-CTs for each adaptation were required. If CBCTs, which are already used for patient positioning, or CBCT-based sCTs could be used for the treatment plan adaptation, the extra dose of the 4D-CTs could be avoided.

In some treatment fractions the OAR-constraints were violated due to tumor movement towards the heart and lung in two patients. Furthermore, lung volume changes caused higher V20Gy of the ipsilateral lung in some fractions. To gain the initial coverage, the OAR doses were exceeded in the isoeffective ART scenario.

Benefits of isotoxic ART

It is still a matter of discussion which patients could benefit from an isotoxic ART approach with a target dose-escalation and which upper dose limits are appropriate. Improved outcomes in dose-escalation studies for patients who received RT alone or after inducting chemotherapy up to 84 Gy [47] and 103 Gy [26] were found. Ramroth et al. showed a prolonged median survival with dose-escalation for patients without chemotherapy in a meta-analysis involving 3795 patients from 25 trials. The time-corrected EQD2 ranged from 36.4 Gy to 80.8 Gy. With cCRT the opposite occurred which might be due to higher levels of toxicity [48]. In the study trial RTOG-0617 patients with cCRT who were treated with a Dpresc of 60 Gy had better OS rates than patients who were treated with 74 Gy [5]. Schild et al. showed the more frequent occurrence of adverse events grade III or greater with dose-escalation in a meta-analysis involving 3600 patients from 16 trials [49]. Different to the RTOG-0617 and the study of Schild et al. for the isotoxic ART scenario, a fractional dose-escalation was chosen in this work to avoid tumor repopulation due to treatment prolongation [11, 50, 51] and the OAR-doses should not exceed the initially planned values to avoid the increased occurrence of adverse events.

The median D95%(GTV) was 10.0 Gy ± 8.1 Gy higher compared to the initial plans. That is a larger dose increase than the result of the isotoxic ART in the prior mentioned study by Guckenberger et al. with a mean escalation of the D95%(GTV) by 5.7 Gy. Remarkably, this dose-escalation was achieved by adapting only twice in the treatment period [46]. The median EQD2(D95%(PTV)) was 8.1 Gy ± 13.7 Gy higher than the initial treatment plans. That is lower than the mean EQD2(D95%(PTV)) dose-escalation of 13.4 Gy which was obtained by Weiss et al. by adapting treatment plans of 10 patients with NSCLC in weeks 2 and 4. That they achieved a higher dose-escalation with only 2 adaptations might be explainable by our dosimetric limit of the D95%(PTV), the larger GTV-reduction they reported, the chosen α/β-ratio, and furthermore, they escalated the dose until the MLD exceeded the initial dose by 1 Gy [52]. The ongoing RTOG-1106 study compares normofractionated patients with patients for whom a mid-treatment PET-CT was used for individualized dose-escalation. The median target dose with one adaptation was 71 Gy. The 2-year in-field local–regional control and primary tumor control of the patients in the adaptation arm was 10.7% and 17.1% higher than for the normofractionated patients [35]. An adaptation based on a PET-CT to the residual active volume might lead to more promising results, but is time-consuming and leads to additional doses which cannot be neglected and therefore, are unsuitable for daily adaptations [45, 53].

The normalization of the isotoxic ART treatment plans to the V20Gy of the ipsilateral lung could lead to an excess of the remaining DVH parameters in cases where the V20Gy was not the limiting OAR-constraint. Furthermore, we allowed dose normalization ratios between 90% and 110%, which could also cause a violation of OAR constraints in the isotoxic ART scenario.

Limitations and benefits of the proposed approach

Having one physician performing all image segmentation reduces possible bias introduced by inter-individual assessments. Furthermore, contouring retrospectively in an unsorted sCT-order minimized a potential bias in favour of any method. Nevertheless, the volume changes of the GTV seen in Fig. 2 cannot be solely attributed to daily changes in tumor size, but could also be influenced by the resolution of the sCT and intraobserver variability of the structure delineation, which, however, is already covered by the GTV-PTV margin. Furthermore, accurate image segmentation of the heart was expected to be challenging since it is mostly located at the periphery of CBCT scans, where soft tissue discrimination is weak. Therefore, dosimetric consequences in this structure will be more likely prone to errors. Daily fluctuations of certain DVH-parameters as displayed in Fig. 3 demonstrate the difficulty of a reproducible treatment scenario and subsequent image segmentation. Especially daily fluctuations in the structure which is used for dose normalization (PTV or ipsilateral lung) will lead to dosimetric deviations. This is of particular interest when adaptations are performed on a regular but not daily basis, e.g. once weekly [28].

The relatively small FOV of conventional CBCT scanners resulted in only partially visible lungs for three patients. A larger FOV would be advantageous but might imply higher imaging doses.

Although daily ART becomes feasible with modern linear accelerators it is typically still time-consuming. The benefit of ART decreases from midterm to weekly and daily ART [28]. Woodford et al. stated that a tumor regression of at least 30% after 20 treatment fractions is worth replanning. A predictor of the time points for every patient, where the replanning is worth the effort or a predictor for OAR toxicities would be beneficial. This however, implies that pre-treatment volume images are available with a sufficiently high image quality.

Changes in the fractional dose of the isotoxic treatment plans could affect the patient’s anatomy differently than the normofractionated ones. Higher fractional doses could further increase the effect of the GTV-regression and therefore, as well of the dose-escalation [42]; nonetheless, fractionation-related effects remain yet uncertain. In addition, the tested ART strategies are pre-treatment adaptations, which address only interfractional and not intrafractional anatomic changes, thus the effect is larger than that of offline ART strategies, but could underestimate the effect of real-time online ART Strategies.

Our approach for two different treatment strategies namely, isoeffective ART and isotoxic ART was feasible, confirmed by an accurate dosimetric analysis of daily ART. All treatment plans were calculated for the same patient collective which makes them reproducible. In the isoeffective scenario the OAR-doses were spared with adequate target-coverage in all treatment fractions. In the isotoxic treatment plans the dose to the target volume could be escalated without relevantly increasing the OAR exposure.

Potential future implementations

sCT-guided ART can be implemented on conventional linear accelerators that are equipped with a CBCT scanner. However, a separate CBCT to sCT conversion algorithm might need to be trained for different CBCT scanners [22]. Alternatively, some of the new linear accelerators, such as the Ethos (Varian Medical Systems, Palo Alto, CA), are already equipped for online ART [54].

Other possibilities are adaptations based on PET-CT or CT, which, however, are associated with additional radiation exposure [45, 53], or MR-guided adaptations. The advantages of adapted MR-guided radiotherapy are that MR imaging is possible without ionizing radiation [55, 56], MR images have a better soft tissue contrast than in CT, which facilitates delineation of the OAR and the target [57,58,59,60], and MR-guided ART allows for real-time monitoring [58, 61] and can be combined with functional imaging to assess tumor response to therapy. This offers the possibility to adjust the therapy based on the biological information [58]. The disadvantages of MR-guided ART are the low availability of MR linacs [62], the required manpower that has to be present during the entire treatment [60], and the treatment duration [61,62,63].

Conclusion

Lung-IGRT results in appropriate interfractional target-coverage, at expenses of higher OAR exposure. Both online issoeffective and isotoxic ART achieve adequate OAR-sparing dosimetry, whereas the latter allows dose-escalated PTV coverage. Ongoing research will elucidate the clinical applicability of this approach.

Availability of data and materials

The datasets used and analysed during the current study are available from the corresponding author on reasonable request.

References

  1. Sonke JJ, Aznar M, Rasch C. Adaptive radiotherapy for anatomical changes. Semin Radiat Oncol. 2019;29(3):245–57.

    Article  PubMed  Google Scholar 

  2. Yamamoto T, et al. Changes in regional ventilation during treatment and dosimetric advantages of CT ventilation image guided radiation therapy for locally advanced lung cancer. Int J Radiat Oncol Biol Phys. 2018;102(4):1366–73.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Sonke JJ, Belderbos J. Adaptive radiotherapy for lung cancer. Semin Radiat Oncol. 2010;20(2):94–106.

    Article  PubMed  Google Scholar 

  4. Speirs CK, et al. Heart dose is an independent dosimetric predictor of overall survival in locally advanced non-small cell lung cancer. J Thorac Oncol. 2017;12(2):293–301.

    Article  PubMed  Google Scholar 

  5. Bradley JD, et al. Standard-dose versus high-dose conformal radiotherapy with concurrent and consolidation carboplatin plus paclitaxel with or without cetuximab for patients with stage IIIA or IIIB non-small-cell lung cancer (RTOG 0617): a randomised, two-by-two factorial phase 3 study. Lancet Oncol. 2015;16(2):187–99.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Wang K, et al. Heart dosimetric analysis of three types of cardiac toxicity in patients treated on dose-escalation trials for stage III non-small-cell lung cancer. Radiother Oncol. 2017;125(2):293–300.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Appel S, et al. Image-guidance triggered adaptive replanning of radiation therapy for locally advanced lung cancer: an evaluation of cases requiring plan adaptation. Br J Radiol. 2020;93(1105):20190743.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Caglar HB, Othus M, Allen AM. Esophagus in-field: a new predictor for esophagitis. Radiother Oncol. 2010;97(1):48–53.

    Article  PubMed  Google Scholar 

  9. Kaymak-Cerkesli ZA, Ozkan EE, Ozseven A. The esophageal dose-volume parameters for predicting Grade I–II acute esophagitis correlated with weight loss and serum albumin decrease in lung cancer radiotherapy. J Cancer Res Ther. 2021;17(1):94–8.

    Article  PubMed  Google Scholar 

  10. Moller DS, et al. Heterogeneous FDG-guided dose-escalation for locally advanced NSCLC (the NARLAL2 trial): design and early dosimetric results of a randomized, multi-centre phase-III study. Radiother Oncol. 2017;124(2):311–7.

    Article  PubMed  Google Scholar 

  11. Ma L, et al. A current review of dose-escalated radiotherapy in locally advanced non-small cell lung cancer. Radiol Oncol. 2019;53(1):6–14.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Bradley JD, et al. Long-term results of NRG oncology RTOG 0617: standard- versus high-dose chemoradiotherapy with or without cetuximab for unresectable stage III non-small-cell lung cancer. J Clin Oncol. 2020;38(7):706–14.

    Article  CAS  PubMed  Google Scholar 

  13. Spigel DR, et al. Five-year survival outcomes from the PACIFIC trial: durvalumab after chemoradiotherapy in stage III non-small-cell lung cancer. J Clin Oncol. 2022;40(12):1301–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Benedict SH, et al. Stereotactic body radiation therapy: the report of AAPM Task Group 101. Med Phys. 2010;37(8):4078–101.

    Article  PubMed  Google Scholar 

  15. Lim-Reinders S, et al. Online adaptive radiation therapy. Int J Radiat Oncol Biol Phys. 2017;99(4):994–1003.

    Article  PubMed  Google Scholar 

  16. Moller DS, et al. Adaptive radiotherapy for advanced lung cancer ensures target coverage and decreases lung dose. Radiother Oncol. 2016;121(1):32–8.

    Article  PubMed  Google Scholar 

  17. Kaplan LP, et al. Cone beam CT based dose calculation in the thorax region. Phys Imaging Radiat Oncol. 2018;7:45–50.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Giacometti V, Hounsell AR, McGarry CK. A review of dose calculation approaches with cone beam CT in photon and proton therapy. Phys Med. 2020;76:243–76.

    Article  PubMed  Google Scholar 

  19. Schröder L, et al. Evaluating the impact of cone-beam computed tomography scatter mitigation strategies on radiotherapy dose calculation accuracy. Phys Imaging Radiat Oncol. 2019;10:35–40.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Hansen DC, et al. ScatterNet: a convolutional neural network for cone-beam CT intensity correction. Med Phys. 2018;45(11):4916–26.

    Article  CAS  PubMed  Google Scholar 

  21. Jiang Y, et al. Scatter correction of cone-beam CT using a deep residual convolution neural network (DRCNN). Phys Med Biol. 2019;64(14):145003.

    Article  PubMed  Google Scholar 

  22. Eckl M, et al. Evaluation of a cycle-generative adversarial network-based cone-beam CT to synthetic CT conversion algorithm for adaptive radiation therapy. Phys Med. 2020;80:308–16.

    Article  PubMed  Google Scholar 

  23. Wang K, et al. Cardiac toxicity after radiotherapy for stage III non-small-cell lung cancer: pooled analysis of dose-escalation trials delivering 70 to 90 Gy. J Clin Oncol. 2017;35(13):1387–94.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Moller DS., et al. Survival benefits for non-small cell lung cancer patients treated with adaptive radiotherapy. Radiother Oncol. 2022.

  25. Vivekanandan S, et al. The impact of cardiac radiation dosimetry on survival after radiation therapy for non-small cell lung cancer. Int J Radiat Oncol Biol Phys. 2017;99(1):51–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Kong FM, et al. High-dose radiation improved local tumor control and overall survival in patients with inoperable/unresectable non-small-cell lung cancer: long-term results of a radiation dose escalation study. Int J Radiat Oncol Biol Phys. 2005;63(2):324–33.

    Article  PubMed  Google Scholar 

  27. Machtay M, et al. Higher biologically effective dose of radiotherapy is associated with improved outcomes for locally advanced non-small cell lung carcinoma treated with chemoradiation: an analysis of the Radiation Therapy Oncology Group. Int J Radiat Oncol Biol Phys. 2012;82(1):425–34.

    Article  PubMed  Google Scholar 

  28. Dial C, et al. Benefits of adaptive radiation therapy in lung cancer as a function of replanning frequency. Med Phys. 2016;43(4):1787.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Haslett K, et al. Isotoxic intensity modulated radiation therapy in stage III non-small cell lung cancer: a feasibility study. Int J Radiat Oncol Biol Phys. 2021;109(5):1341–8.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Kong FM, et al. RTOG0617 to externally validate blood cell ERCC1/2 genotypic signature as a radiosensitivity biomarker for both tumor and normal tissue for individualized dose prescription. Int J Radiat Oncol Biol Phys. 2020;108(3):S2.

    Article  Google Scholar 

  31. Boda-Heggemann J, et al. Deep inspiration breath hold-based radiation therapy: a clinical review. Int J Radiat Oncol Biol Phys. 2016;94(3):478–92.

    Article  PubMed  Google Scholar 

  32. Ning MS, et al. Incidence and predictors of pericardial effusion after chemoradiation therapy for locally advanced non-small cell lung cancer. Int J Radiat Oncol Biol Phys. 2017;99(1):70–9.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Atkins KM, et al. Cardiac radiation dose, cardiac disease, and mortality in patients with lung cancer. J Am Coll Cardiol. 2019;73(23):2976–87.

    Article  PubMed  Google Scholar 

  34. Vogel L, et al. Intra-breath-hold residual motion of image-guided DIBH liver-SBRT: an estimation by ultrasound-based monitoring correlated with diaphragm position in CBCT. Radiother Oncol. 2018;129(3):441–8.

    Article  PubMed  Google Scholar 

  35. Kong F-MS., et al. NRG-RTOG 1106/ACRIN 6697: a phase IIR trial of standard versus adaptive (mid-treatment PET-based) chemoradiotherapy for stage III NSCLC: results and comparison to NRG-RTOG 0617 (non-personalized RT dose escalation). J. Clin. Oncol. 2021; 39(15_suppl): 8548.

  36. Loap P, Fourquet A, Kirova Y. The limits of the linear quadratic (LQ) model for late cardiotoxicity prediction: example of hypofractionated rotational intensity modulated radiation therapy (IMRT) for breast cancer. Int J Radiat Oncol Biol Phys. 2020;106(5):1106–8.

    Article  PubMed  Google Scholar 

  37. van Leeuwen CM, et al. The alfa and beta of tumours: a review of parameters of the linear-quadratic model, derived from clinical radiotherapy studies. Radiat Oncol. 2018;13(1):96.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Schultheiss TE, Tome WA, Orton CG. Point/counterpoint: it is not appropriate to “deform” dose along with deformable image registration in adaptive radiotherapy. Med Phys. 2012;39(11):6531–3.

    Article  PubMed  Google Scholar 

  39. Chetty IJ, Rosu-Bubulac M. Deformable registration for dose accumulation. Semin Radiat Oncol. 2019;29(3):198–208.

    Article  PubMed  Google Scholar 

  40. Amstutz F., et al. An approach for estimating dosimetric uncertainties in deformable dose accumulation in pencil beam scanning proton therapy for lung cancer. Phys. Med. Biol. 2021;66(10).

  41. Zhong H, Chetty IJ. Caution must be exercised when performing deformable dose accumulation for tumors undergoing mass changes during fractionated radiation therapy. Int J Radiat Oncol Biol Phys. 2017;97(1):182–3.

    Article  PubMed  Google Scholar 

  42. Zhong H, et al. Evaluation of adaptive treatment planning for patients with non-small cell lung cancer. Phys Med Biol. 2017;62(11):4346–60.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Luo J. et al. Study of the cumulative dose between fractions of lung cancer radiotherapy based on CT and CBCT image deformable registration technology. Front. Phys. 2020; 8.

  44. Marks LB, et al. Use of normal tissue complication probability models in the clinic. Int J Radiat Oncol Biol Phys. 2010;76(3 Suppl):S10–9.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Hoegen P, et al. Cone-beam-CT guided adaptive radiotherapy for locally advanced non-small cell lung cancer enables quality assurance and superior sparing of healthy lung. Front Oncol. 2020;10:564857.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Guckenberger M, et al. Potential of adaptive radiotherapy to escalate the radiation dose in combined radiochemotherapy for locally advanced non-small cell lung cancer. Int J Radiat Oncol Biol Phys. 2011;79(3):901–8.

    Article  PubMed  Google Scholar 

  47. Rosenzweig KE, et al. Results of a phase I dose-escalation study using three-dimensional conformal radiotherapy in the treatment of inoperable nonsmall cell lung carcinoma. Cancer. 2005;103(10):2118–27.

    Article  PubMed  Google Scholar 

  48. Ramroth J, et al. Dose and fractionation in radiation therapy of curative intent for non-small cell lung cancer: meta-analysis of randomized trials. Int J Radiat Oncol Biol Phys. 2016;96(4):736–47.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Schild SE, et al. Toxicity related to radiotherapy dose and targeting strategy: a pooled analysis of cooperative group trials of combined modality therapy for locally advanced non-small cell lung cancer. J Thorac Oncol. 2019;14(2):298–303.

    Article  CAS  PubMed  Google Scholar 

  50. Zehentmayr F, et al. Radiation dose escalation with modified fractionation schedules for locally advanced NSCLC: a systematic review. Thorac Cancer. 2020;11(6):1375–85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. McMillan MT, et al. Radiation treatment time and overall survival in locally advanced non-small cell lung cancer. Int J Radiat Oncol Biol Phys. 2017;98(5):1142–52.

    Article  PubMed  Google Scholar 

  52. Weiss E, et al. Dose escalation for locally advanced lung cancer using adaptive radiation therapy with simultaneous integrated volume-adapted boost. Int J Radiat Oncol Biol Phys. 2013;86(3):414–9.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Smith-Bindman R, et al. Radiation dose associated with common computed tomography examinations and the associated lifetime attributable risk of cancer. Arch Intern Med. 2009;169(22):2078–86.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Mao W, et al. Evaluation of auto-contouring and dose distributions for online adaptive radiation therapy of patients with locally advanced lung cancers. Pract Radiat Oncol. 2022;12(4):e329–38.

    Article  PubMed  Google Scholar 

  55. Ahmad SB, et al. Evaluation of a commercial MRI Linac based Monte Carlo dose calculation algorithm with GEANT4. Med Phys. 2016;43(2):894–907.

    Article  PubMed  Google Scholar 

  56. Randall JW., et al. Towards accurate and precise image-guided radiotherapy: clinical applications of the MR-Linac. J Clin Med, 2022; 11(14).

  57. Kong FM, et al. Effect of Midtreatment PET/CT-adapted radiation therapy with concurrent chemotherapy in patients with locally advanced non-small-cell lung cancer: a phase 2 clinical trial. JAMA Oncol. 2017;3(10):1358–65.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Crockett CB, et al. Initial clinical experience of MR-guided radiotherapy for non-small cell lung cancer. Front Oncol. 2021;11:617681.

    Article  PubMed  PubMed Central  Google Scholar 

  59. Schiff JP, et al. A pilot study of same-day MRI-only simulation and treatment with MR-guided adaptive palliative radiotherapy (MAP-RT). Clin Transl Radiat Oncol. 2023;39:100561.

    Article  CAS  PubMed  Google Scholar 

  60. Murray J, Tree AC. Prostate cancer: advantages and disadvantages of MR-guided RT. Clin Transl Radiat Oncol. 2019;18:68–73.

    Article  PubMed  PubMed Central  Google Scholar 

  61. Botman R, et al. The clinical introduction of MR-guided radiation therapy from a RTT perspective. Clin Transl Radiat Oncol. 2019;18:140–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Garcia Schuler HI, et al. Operating procedures, risk management and challenges during implementation of adaptive and non-adaptive MR-guided radiotherapy: 1-year single-center experience. Radiat Oncol. 2021;16(1):217.

    Article  PubMed  PubMed Central  Google Scholar 

  63. Henke LE, et al. Stereotactic MR-guided online adaptive radiation therapy (SMART) for ultracentral thorax malignancies: results of a phase 1 trial. Adv Radiat Oncol. 2019;4(1):201–9.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

The authors would like to thank Gustav Meedt, Nicolette O’Connell, Peter Voet and Jiaofeng Xu from Elekta AB for their continuous support of this project.

Funding

Open Access funding enabled and organized by Projekt DEAL. The project was partially funded by an Elekta AB research grant. The research software ADMIRE was provided by Elekta in a research version of Monaco.

Author information

Authors and Affiliations

  1. Department of Radiation Oncology, University Medical Center Mannheim, University of Heidelberg, Theodor-Kutzer-Ufer 1-3, 68167, Mannheim, Germany

    Lea Hoppen, Chung S. Kwok, Judit Boda-Heggemann, Daniel Buergy, Michael Ehmann, Frank A. Giordano & Jens Fleckenstein

  2. Department of Radiation Oncology, University Hospital Bonn, University of Bonn, Bonn, Germany

    Gustavo R. Sarria

Authors
  1. Lea Hoppen
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Gustavo R. Sarria
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Chung S. Kwok
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Judit Boda-Heggemann
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Daniel Buergy
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Michael Ehmann
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Frank A. Giordano
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Jens Fleckenstein
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

L.H. and J.F. conceived the study and prepared the manuscript. L.H. and C.S.K. exported datasets, generated sCTs and treatment plans and evaluated the results. The structures were delineated by G.R.S. J.B., V.S., M.E., F.A.G. and J.F. contributed intellectually to the manuscript and critically reviewed it. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Lea Hoppen.

Ethics declarations

Ethics approval

This investigation was performed according to the principles of the Declaration of Helsinki and after IRB approval (2018-836R-MA). All data was anonymized prior to inclusion.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hoppen, L., Sarria, G.R., Kwok, C.S. et al. Dosimetric benefits of adaptive radiation therapy for patients with stage III non-small cell lung cancer. Radiat Oncol 18, 34 (2023). https://doi.org/10.1186/s13014-023-02222-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13014-023-02222-7

Keywords

Radiation Oncology

ISSN: 1748-717X

Contact us