Skip to main content

ORIGINAL RESEARCH article

Front. Oncol., 24 August 2021
Sec. Radiation Oncology
This article is part of the Research Topic Automation and Artificial Intelligence in Radiation Oncology View all 22 articles

Replacing Manual Planning of Whole Breast Irradiation With Knowledge-Based Automatic Optimization by Virtual Tangential-Fields Arc Therapy

Roberta Castriconi&#x;Roberta Castriconi1†Pier Giorgio Esposito&#x;Pier Giorgio Esposito1†Alessia TuddaAlessia Tudda1Paola MangiliPaola Mangili1Sara BroggiSara Broggi1Andrei FodorAndrei Fodor2Chiara L. DeantoniChiara L. Deantoni2Barbara LongobardiBarbara Longobardi1Marcella PasettiMarcella Pasetti2Lucia PernaLucia Perna1Antonella del VecchioAntonella del Vecchio1Nadia Gisella Di MuzioNadia Gisella Di Muzio2Claudio Fiorino*Claudio Fiorino1*
  • 1Medical Physics, San Raffaele Scientific Institute, Milano, Italy
  • 2Radiotherapy, San Raffaele Scientific Institute, Milano, Italy

Purpose: To implement Knowledge Based (KB) automatic planning for right and left-sided whole breast treatment through a new volumetric technique (ViTAT, Virtual Tangential-fields Arc Therapy) mimicking conventional tangential fields (TF) irradiation.

Materials and Method: A total of 193 clinical plans delivering TF with wedged or field-in-field beams were selected to train two KB-models for right(R) and left(L) sided breast cancer patients using the RapidPlan (RP) tool implemented in the Varian Eclipse system. Then, a template for ViTAT optimization, incorporating individual KB-optimized constraints, was interactively fine-tuned. ViTAT plans consisted of four arcs (6 MV) with start/stop angles consistent with the TF geometry variability within our population; the delivery was completely blocked along the arcs, apart from the first and last 20° of rotation for each arc. Optimized fine-tuned KB templates for automatic plan optimization were generated. Validation tests were performed on 60 new patients equally divided in R and L breast treatment: KB automatic ViTAT-plans (KB-ViTAT) were compared against the original TF plans in terms of OARs/PTVs dose-volume parameters. Wilcoxon-tests were used to assess the statistically significant differences.

Results: KB models were successfully generated for both L and R sides. Overall, 1(3%) and 7(23%) out of 30 automatic KB-ViTAT plans were unacceptable compared to TF for R and L side, respectively. After the manual refinement of the start/stop angles, KB-ViTAT plans well fitted TF-performances for these patients as well. PTV coverage was comparable, while PTV D1% was improved with KB-ViTAT by R:0.4/L:0.2 Gy (p < 0.05); ipsilateral OARs Dmean were similar with a slight (i.e., few % volume) improvement/worsening in the 15–35 Gy/2–15 Gy range, respectively. KB-ViTAT better spared contralateral OARs: Dmean of contralateral OARs was 0.1 Gy lower (p < 0.05); integral dose was R:5%/L:8% lower (p < 0.05) than TF. The overall time for the automatic plan optimization and final dose calculation was 12 ± 2 minutes.

Conclusions: Fully automatic KB-optimization of ViTAT can efficiently replace manually optimized TF planning for whole breast irradiation. This approach was clinically implemented in our institute and may be suggested as a large-scale strategy for efficiently replacing manual planning with large sparing of time, elimination of inter-planner variability and of, seldomly occurring, sub-optimal manual plans.

Introduction

Post-operative irradiation of the whole breast is a well assessed and effective therapeutic option in the treatment of localized breast cancer (1). Typically, more than 70% of women submitted to breast-conserving surgery receives post-operative radiotherapy and a large fraction of them is treated to sterilize the whole breast, typically delivering 38–40 Gy in 15–16 fractions or 50 Gy in 25 fractions. Due to the large and increasing incidence of breast cancer in the female population (2), this treatment represents a quite relevant fraction of the patients daily treated in the radiation oncology departments worldwide. Despite the evolution toward more personalized approaches including the possibility to reduce the treated volumes (i.e., partial breast irradiation) or to deliver higher dose to the tumor bed or to include selected nodal regions at risk, whole breast irradiation still maintains a central role in the treatment of breast cancer. Nowadays, different techniques are used; most Institutes still prefer the conventional tangential fields (TF) arrangement, either using 3-dimensional conformal radiation therapy (3DCRT) with wedges (physical or dynamic) or intensity modulated radiotherapy (IMRT), often delivered with a few segments per beam, mostly manually optimized (36). Rotational techniques (7, 8) generally showed better performances in better tailoring the dose distribution to the Planning Target Volume (PTV) shape with a consequent improved sparing of organs at risk (OARs) at high-intermediate dose levels, especially in the case of concave-shaped PTVs. Nevertheless, the issues related to the potential clinical meaning of the higher low-dose spread to the heart, lungs, and contralateral breast with rotational techniques is still open (9, 10). Due to this, the TF arrangement, limiting any relevant low-dose spread, is expected to stably remain among the most used techniques to treat breast cancer also in the next decade. On the other hand, manual (and also inverse planned) optimization is time consuming and highly dependent on the planner skill (11). Standardization in radiotherapy treatment planning is an important goal aimed to guarantee to all patients a high quality treatment, independent of the planner time and skills; this seems still more urgent in countries with a rapidly growing incidence of cancer and low/middle income (2). Automated planning solutions were recently introduced (1225) with the aim of reducing planning time and inter-operator variability while conserving (or improving) a high quality plan (23, 2629). Many systems were largely investigated for a different clinical situation including the breast site (11, 17, 3035). Regarding the TF approach for whole breast irradiation, a relatively weak point of auto-planning is the difficulty to take into account the inter-patient variations in assessing the best position of the fields to limit the dose to the adjacent organs, concomitantly assuring PTV coverage and highly homogenous dose distribution within PTV. Consequently, automatic solutions for this application were rarely reported using in-house systems, intrinsically difficult to extend on a larger scale (11, 30, 32).

Automating the optimization of TF by knowledge-based (KB) approaches would be suitable as the (largely available) past information could be modelled to be optimally applied on new patients. The RapidPlan® system (Varian Inc.) is available commercially and widely tested in many clinical scenarios, including breast VMAT (17, 31, 33). However, no clinical examples for the TF irradiation of whole breast are reported likely due to the configuration of the system implemented for IMRT/VMAT inverse-planning optimization. In order to obtain this objective, within a project for the large-scale implementation of automatic KB plan optimization (MIKAPOCo, Multi-Institutional Knowledge-based Approach to Plan Optimization for the Community), we previously demonstrated the possibility to mimic (and slightly improve) the performances of TF irradiation through a partially blocked multi-arcs approach using RapidArc®, named ViTAT [Virtual Tangential-fields Arc Therapy (36)]. Aims of the present work were:

1. to develop KB-models based on TF plans aimed to drive inverse-planned ViTAT plans; and

2. to demonstrate the possibility of replacing TF manual optimization with an automatic KB-ViTAT plan optimization.

Material and Methods

Clinical and Planning Procedures for Tangential Field Planning

At our institute, breast cancer patients treated with whole breast irradiation receives 40 Gy in 15 fractions (2.67 Gy/fr), prescribed as a mean dose to PTV. CTV and PTV are defined according to the AIRO national guidelines (36): in short, CTV is contoured excluding the skin, with a 5-mm margin from the surface while PTV is obtained by a 5-mm isotropic expansion. PTV is finally cropped with a 5-mm margin from the body surface. During the last 10 years, patients were treated mostly with opposed or quasi-opposed, BEV-based (Beam’s Eye View) optimized, fields, using physical wedges. Manually optimized segments were generally added to improve the PTV coverage and dose distribution homogeneity. Starting from 2018, this technique was gradually replaced by the manually optimized field-in-field technique (3), avoiding wedges. The number of segments for each field ranged between one (only open fields) and four, with most patients optimized using 2–3 segments per field. Moreover, for both techniques, the gantry position of few segments could be slightly changed (between 5° and 10°) with respect to the “prevalent” field, aiming to limit hot spots and improve homogeneity within PTV and OARs sparing. Based on the internally conducted plan comparison performed in 2017, the differences between FIF and our TF approach using physical wedges were clinically negligible, in slight favor of the FIF. Patients were treated with 6 MV beams, with few exceptions when 18 MV and 6 MV fields/segments were combined. During the last 5 years, all patients were treated with a Varian CLINAC-IX 2300 equipped with a 5-mm Millenium-MLC system; daily image-guidance with CBCT was performed for all patients; plans were optimized using the Varian Eclipse TPS system (v. 13.6) using the AcurosXB® algorithm for dose calculation.

The planning goals for PTV were: V95% > 95% (the fraction of PTV volume receiving more than 95% of the prescribed dose higher that 95%) and Dmax < 108%. OARs were defined according to the AIRO national guidelines and always included the contralateral breast, the contralateral lung, the ipsilateral lung, and the heart: constraints were V20Gy < 20% for the ipsilateral lung and V5Gy < 40%, Dmean < 5 Gy for the heart in the left breast case. Independently from these constraints, planners always tried to reduce the dose to OARs as much as possible while respecting the goals for PTV. Contralateral breast was always avoided by the medial entry beam, while the lateral beam could include small portions at the exit. In the case the constraints could not be respected, PTV coverage received a higher priority against OARs sparing, with few exceptions. Alternatively, a VMAT plan could also be optimized and chosen by the physician.

The ViTAT Technique

As largely explained in Esposito et al. (36), we previously demonstrated the feasibility to mimic the TF irradiation performances through a multi-arcs approach where the delivery was partly blocked using RapidArc® (RA), named Virtual Tangential-fields Arc Therapy (ViTAT). Figure 1 shows the planning ViTAT set-up.

FIGURE 1
www.frontiersin.org

Figure 1 ViTAT setup: 3D view of a right side ViTAT plan showing the geometry of the four arcs used (on the left)—the segments show the beam delivery while the rest of the arcs are blocked for delivery; beam eye view of the medial angle (on the right).

In short, ViTAT consisted of four arcs (6 MV) with a collimator rotation angle of ±10° and optimized with a RA technique with start/stop angles of 60°/220° for the right side and 300°/135° for the left one, consistent with the TF geometry in our population. Irradiation through the arcs was completely blocked, apart in the first and the last 20° of rotation. As for the breast VMAT-optimization (17), Target structure for plan optimization is obtained expanding the PTV outside the body (1.5 cm expansion in the external-lateral and anterior direction) to account for any breathing and inter-fraction deformation effect. A virtual bolus (1.5 cm thickness, -500 HU density) (17, 37) is linked to the fields during optimization to fill the Target structure during the optimization in order to avoid effects due to the electronic non-equilibrium and to assure a proper safety set-up margin (38, 39). Of course, the bolus is removed before the final dose calculation. Field dimensions are BEV-based adapted to the individual patient anatomy. Plans are optimized with the Rapid Arc (RA) optimization module by inverse planning. We implemented the ViTAT technique in this study with the goal of replacing the manual planning optimization in TF irradiation through a KB-optimization approach.

KB-Model Generation

In order to guarantee a better plan homogeneity, TF clinical plans of the period 2016–2019 were considered for building the model, resulting in 90 plans (70 with wedged fields and 20 with FIF) for the right-side breast. In the same way, 103 (88 with wedged fields and 15 FIF) TF plans in the range 2016–2020 were selected for the left-side breast. These plans were used to train two KB-models, one for the right-side and one for the left-side breast using the RapidPlan (RP) tool implemented in the Eclipse TPS (v 13.6, Varian Inc.). This choice also balanced the need to train a sufficient number of plans and to keep a sufficiently large validation group, considering both the recent FIF (2019–2020) and older TF (before 2016) plans, as explained later. RP is configured to model inverse plans delivered with the IMRT/VMAT technique. Hence, each TF dose distribution was linked to a virtual RA-plan (40) consisting of two partially reverse arcs in the range 60°/220° for the right-side and 300°/135° for the left-side with a collimator rotation angle of ±10°. The choice of the start/stop angles was optimized according to the distribution of the beam angles of the original plans, as explained in (36). During KB modelling, arcs were not considered to be blocked (no avoidance sector were used), so the entire geometry of the structure is seen from the Beams Eye View (BEV) of the arc. Before accomplishing model configuration, the impact of using multiple arcs for building the models on DVH-estimates was investigated. As the partitions of the OAR during the training phase are the same for equal arcs, no differences were observed in the DVH-estimates for model configuration using two arcs instead of four ones (as used for the ViTAT techniques). The OARs considered for training the models were the ipsilateral lung and contralateral breast for the right side plus the heart for the left side.

Right-side heart and contralateral lung were not modelled as they were considered to be more easily managed in the final template (with the same fixed constraints for all patients) without any support of the prediction of the model, due to their anatomical position treatment and to the ballistic of the arc orientations. The tuning of the models was performed by using statistical tools available in the RP system and the Model Analytic platform: the methods followed to limit the impact of outliers were reported elsewhere (25, 41, 42). In short, for each OARs, the features that exceeded by >2 SD the principal components fitted with the KB-models were identified and the original clinical plan were individually re-evaluated; they were excluded as “dosimetric outliers” only in case of recognized sub-optimal planning. All the other potential outliers were found to be representative of an “uncommon” but clinically suitable geometry/anatomy condition and were kept in the models.

Template for Automatic ViTAT Plan Optimization

Based on the resulting prediction models, the system generates the confidence intervals of the expected DVH for each OAR and suggests the lower confidence values to be used for plan optimization. In order to obtain a robust and efficient template, the choice of the position and penalty of the generated DVH constraints needs to be optimized. The optimization of the KB-based template for each side was carried out by several fine-tuning tests based on repetitive automatic plans for five sample patients for each of the two situations (right and left), testing the impact of the position of DVH constraints and their penalties, as similarly reported in previous studies (2325). Based on our experience, once the PTV priorities are fixed, the OARs priorities were gradually increased through three KB-test templates. Selected relevant dose-volume parameters were analyzed for the five sample patients and compared against the original clinical plans to assess the performances of the three KB-based templates and to assess the one showing the best performances.

Validation Tests

Validation tests were performed on 60 new patients equally divided between right-side and left-side breast. The external validation of the models was performed on patients selected in the range 2019–2020 and 2013–2016.In this way, we also intended to validate the model separately for wedged fields plans (the oldest group) and for the more recent FIF plans (delivered in 2019–2020). The number of FiF plans is 15 for right-sided breast patients and 11 for the left-sided ones, while the number of wedged fields is 15 for right-sided breast patients and 19 for left-sided ones. Importantly, the validation assessed the performances of the model to properly adapt the anatomical/morphological features of each individual patient and the feasibility of using TF DVH-prediction to fully automate the ViTAT optimization. According to others (17, 19, 23, 40), the validation was performed by re-optimizing a number of clinical plans and by comparing them against the original ones. All KB fully-automatic ViTAT-plans (KB-ViTAT) were compared against the original plans (TF) in terms of OARs/PTVs dose-volume parameters. The comparison was based on the analysis of the mean dose, maximum dose (D1% for PTV and D2% for OARs), and selected dose-volume parameters extracted from DVH. All selected parameters and DVHs were semi-automatically exported via ESAPI scripts and saved in spreadsheets for analysis. Wilcoxon-tests were performed to assess the statistically significant differences (p < 0.05).

Results

In total, 10 patients (out of 90) for the right side and 18 patients (out of 103) for the left side were excluded by the model resulting in sub-optimal plans for at least one of the considered OARs or for PTV coverage. The final model returned χ2 and R2 parameters for the predicted OARs is shown in Table 1. The resulting DVH-estimates were used to generate an individually optimized KB-template for the ViTAT optimization. For the not trained OARs, the position of DVH constraints and their penalties were fixed and tuned as previously explained. Given the overall results, the KB-based template for automatic planning optimization was finally generated and is shown in Table 2 for the right-side case and in Table 3 for the left-side one.

TABLE 1
www.frontiersin.org

Table 1 Final model goodness parameter for models for both right and left sided breast cancer patients.

TABLE 2
www.frontiersin.org

Table 2 The KB-based template for automatic planning optimization for the ViTAT technique for right-sided breast cancer treatment.

TABLE 3
www.frontiersin.org

Table 3 The KB-based template for automatic planning optimization for the ViTAT technique for left-sided breast cancer treatment.

Mean DVHs comparison between KB-ViTAT plans and clinical ones are shown in Figures 2 and 3 for right and left side, respectively. On 30 KB automatic re-optimized cases per each side, only one for the right case and seven for the left case resulted in plans unacceptability in terms of the PTV coverage and/or ipsilateral lung constraints. The reason for this has to be found in the shape and position of the PTV with respect to OARs, causing an uncovering in the medial part of the PTV. However, after the manual refinement of the start/stop angles of 5°/10° by an expert planner, all automatic KB-ViTAT plans well fitted TF-performances. A manual refinement of 5° for the medial angle was necessary for the right-breast case and four of the left-side patients getting the following angles respectively: 65°/220° and 295°/135°. In the remaining three cases for left-side patients, it was necessary to also change the distal angle by 5° and 10° for two and one patients, respectively. Making the modification in the start/stop angles and re-starting the automated optimization, the PTV coverage and the ipsilateral lung constraints were satisfied and comparable to original clinical TF plan. As an example, Supplementary Figures 1 and 2 (Supplementary Material) show the comparison of the originally automatic KB-ViTAT plan with the original start/stop angle and the refined KB-ViTAT with start angle modified by 5° (lower) for the right-side patient and for one of the left-side patients with manual modification of gantry angles.

FIGURE 2
www.frontiersin.org

Figure 2 Mean-DVH comparison of 30 right-sided breast patient tests between the original clinical TF (solid lines) and fully automatic KB-ViTAT (dashed lines) plans.

FIGURE 3
www.frontiersin.org

Figure 3 Mean-DVH comparison of 30 left-sided breast patient tests between the original clinical TF (solid lines) and fully automatic KB-ViTAT (dashed lines) plans.

Quantitative analyses are shown in Tables 4 and 5. Overall, differences between TF and KB-ViTAT were small and in slight favor of ViTAT. PTV coverage was similar, while PTV D1% was improved with automatic optimization (p-value < 0.05). KB-ViTAT better spared contralateral OARs with respect to TF: mean dose was lowered by 0.1 Gy for all contralateral OARs, resulting in a decrease of 33% of the mean dose to the contralateral lung, 20% for the right and 20% for the left side contralateral breast, and 14% for the right-side heart.

TABLE 4
www.frontiersin.org

Table 4 Dose-volume parameters comparison (TF vs ViTAT) for the validation cohort test of 30 plans for the right-sided breast case.

TABLE 5
www.frontiersin.org

Table 5 Dose-volume parameters comparison (TF vs ViTAT) for the validation cohort test of 30 plans for the left-sided breast case.

The KB-ViTAT integral (i.e., body) dose was 5% lower than TF for the right case and 8% for the left case. Ipsilateral lung mean dose was identical (Right: 6.8 Gy, p > 0.05; Left: 5.7 Gy, p > 0.05): there was a modest worsening (i.e., few % volume) for the ipsilateral lung in the 2–15 Gy range and a slight improvement in the 20-35 Gy range when considering KB-ViTAT vs TF. The same behaviour was found for the left-sided case for the heart with a modest worsening in the range 2–10 Gy and a slight improvement in the range 15–35 Gy, while on average, delivering the same mean heart dose (2.7 Gy).

Furthermore, in Figures 4 and 5, the histograms representing the distribution of the differences between TF and KB-ViTAT were reported for the selected PTV/OARs dose-volume parameters. Overall, the time for automatic plan optimization and final dose calculation was registered and found to be 12 ± 2 minutes.

FIGURE 4
www.frontiersin.org

Figure 4 Population histograms of the differences between the clinical tangential field TF and automated re-optimized KB-ViTAT plans for the investigated dosimetric parameters for right-sided breast case, showing the V95% parameter for PTV and mean doses for OARs.

FIGURE 5
www.frontiersin.org

Figure 5 Population histograms of the differences between the clinical tangential field TF and automated re-optimized KB-ViTAT plans for the investigated dosimetric parameters for left-sided breast case, showing the V95% parameter for PTV and mean doses for OARs.

Importantly, the model performances were also evaluated separately for wedged fields plans (the oldest group) and the more recent FIF plans (delivered in 2019–2020). No statistically significant differences resulted in terms of the PTV/OARs dose-volume parameters when comparing the differences between KB-ViTAT and TF in the two cohorts of plans.

The KB-based automatic approach for ViTAT was clinically implemented, first for the right-sided breast and, more recently, for the left-sided one. Five plans for the right-side and five ones for left-side breast were preliminarily verified in terms of dose distribution in a planar phantom using a detector matrix. All clinical KB-VITAT plans underwent the same dosimetric verification before treatment. The gamma passing rate in comparing the calculated vs delivered dose maps was larger than 98% for all plans (in total n = 30), considering 3% - 3 mm as criteria, in agreement with our experience (43).

Discussion

Scope of this work was to fully automate the planning of the whole breast treatment through the KB-optimization approach to mimic the performances of the tangential beam technique. The study aimed at assessing the capability of the RapidPlan tool to handle the dose distributions treated without using the VMAT technique. Moreover, we also explored the feasibility of translating the DVH-prediction model into a fully automatic optimized template for ViTAT planning. As previously shown (36), the ViTAT technique was able to generate dose distributions comparable to TF performances, with some slight improvement in the PTV coverage and homogeneity, in the sparing of contralateral OARs and with a mild reduction of the integral body dose. Then, as demonstrated here, once configured and validated, a KB model trained with TF dose distribution can be efficiently implemented for an automatic ViTAT optimization, completely and efficiently replacing manual plan optimization. The entire workflow was optimized aiming to automatize the selection of the beam angles through the VMAT approach but at the same time to avoid the dose-bath at intermediate–low doses typical of the rotational techniques (7). As a matter of fact, the concerns related to the increase of the mean dose to OARs dealt with the risk of increasing radiation-induced secondary malignancies and late cardiac events (1, 7, 44, 45). Few authors suggested a compromise between better conformity/OARs sparing at high doses against OARs sparing at intermediate–low doses to limit the dose bath associated to VMAT by the partial blocking of arcs (35, 46). Although all these approaches obtained a significant reduction of the low-dose bath, none of them followed the goal of mimicking tangential field irradiation, differently from our ViTAT approach. Other authors demonstrated the possibility of automatizing the angle selection for tangential fields (11) but none investigated the possibility of using the KB-approach. With the approach presented here, we demonstrated the feasibility of the complete replacement of manual TF plans with automatic plans with huge improvements in efficiency, in reducing/eliminating inter-planner variability, and in avoiding sub-optimal plans. Concerning this last point, despite a careful quantitative evaluation was not accomplished, it is relevant that about 15% of the clinical TF plans were a priori excluded when building the KB models as due to sub-optimal planning.

As a matter of fact, the automatic workflow that involved the fixed selection of the start/stop angles had proven to be an efficient way to reproduce the TF performances well. Only in about 13% of cases this approach had failed: on 60 KB automatic re-optimized cases, only eight resulted in plans that were unacceptable in terms of the PTV coverage and ipsilateral lung constraints. Seven out of eight referred to the left side were found to be more challenging from the point of view of the choice of the start/stop angles. However, with the manual refinement of the start/stop angles by an expert planner, the resulting automatic KB-ViTAT plans well fitted TF-performances also for these patients. Furthermore, plans were deliverable showing excellent dosimetric verification performances in phantoms. Importantly, a large sparing of planning time was obtained, with an overall time for automatic plan optimization and final dose calculation of 12 ± 2 minutes.

A highly relevant point, worthy to be underlined here, concerns the demonstration of the automation of a largely used technique by using a commercially available tool, making its potential adoption easy for Varian users. This issue is of primary importance, opening the possibility to a large-scale implementation with a consequent large reduction of the repetitive, manual procedures usually followed during a whole breast plan optimization. The possible sharing of KB-models, permitted by this system, should also be considered as an additional opportunity to rapidly spread this approach primarily to Varian users but in principle extendible/adaptable to other delivery systems (24, 40).

From a wider point of view, the large availability of TF plans and the “limited” inter-institute variability of PTV/OARs can make the development of robust KB DVH prediction models easier and, more importantly, the possibility of extending this approach on a multi-institutional scale. This issue is worthy of investigation and is currently under study within the MIKAPOCo (Multi-Institutional Knowledge-based approach to plan optimization for the community) consortium, joining several Italian Institutes.

Another relevant point concerns the possibility to exploit the ViTAT approach more to obtain further improved plan performances by the possibility of training models using the planning data obtained from clinical plans optimized and delivered at our Institute by applying more stressed intensity-modulation techniques: a project to train KB-models using the planning data of patients treated with Tomotherapy (using the Tomo-Direct module, resulting in IMRT-like TF plans) is currently ongoing.

Conclusions

The approach followed here demonstrated the possibility of the complete replacement of manual tangential breast planning with automatic planning, including beam angle choice. Automatic fixed selection of the start/stop angles and KB-driven optimization were found to be an efficient way to well fit TF planning, with evident advantages in terms of time sparing, elimination of inter-planner variability, and of sub-optimal planning. Manual refinement of the start/stop angles was necessary in 13% of patients, with a large unbalance between the right (3%) and left (23%) sides, resulting in an additional 15–20 minutes more compared to the 12 ± 2 minutes spent for planning optimization and dose calculation in the remaining 87% of the patients. Due to its versatility and the use of a commercial system, this approach shows promising applications for a large-scale implementation.

Data Availability Statement

The datasets presented in this article are not readily available because they are data of patients treated in a single center and cannot be shared without transfer agreement. Requests to access the datasets should be directed to Zmlvcmluby5jbGF1ZGlvQGhzci5pdA==.

Ethics Statement

Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author Contributions

RC, PE, and CF contributed with the conception of the study. RC, PE, and AT worked on data base and modeling. PM, SB, BL, and LP contributed in data integrity evaluation. CF, AF, SB, PM, NM, RC, and PE discussed the interpretation of the results. AF, CDA, NM, and MP were clinical responsible for the treatment of the patient. RC, PE, and CF wrote the manuscript. AV, NM, and CF reviewed and edited the manuscript. CF was responsible for the scientific coordination, managing, and project funding. All authors contributed to the article and approved the submitted version.

Funding

The work was supported by an AIRC (Associazione Italiana Ricerca sul Cancro) grant (IG 23150).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2021.712423/full#supplementary-material

Supplementary Figure 1 | Comparison of plans between the originally automatic KB-ViTAT plan with start/stop angle of 60°/220° (upper) and the refined KB-ViTAT with start angle of 65° (lower) for the right-sided breast case. Starting angles are modified in order to obtain an acceptable coverage to PTV.

Supplementary Figure 2 | Comparison of plans between the originally automatic KB-ViTAT with start/stop angle of 300°/135° (upper) and the refined KB-ViTAT with start angle of 295° (lower)for the left-sided breast case. Starting angles are modified in order to obtain an acceptable coverage to PTV and better hearth and lung dose distribution similarity between TF and KB-ViTAT case.

References

1. Darby S, McGale P, Correa C, Taylor C, Arriagada R, Clarke M, et al. Effect of Radiotherapy After Breast-Conserving Surgery on 10-Year Recurrence and 15-Year Breast Cancer Death: Meta-Analysis of Individual Patient Data for 10 801 Women in 17 Randomised Trials. Lancet (2011) 378:1707–16. doi: 10.1016/S0140-6736(11)61629-2

PubMed Abstract | CrossRef Full Text | Google Scholar

2. WHO. Who Report on Cancer: Setting Priorities, Investing Wisely and Providing Care for All. Geneva: World Health Organization (2020) Licence: CC BY-NC-SA 3.0 IGO.

Google Scholar

3. Chui C-S, Hong L, Hunt M, McCormick B. A Simplified Intensity Modulated Radiation Therapy Technique for the Breast. Med Phys (2002) 29:522–9. doi: 10.1118/1.1460875

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Pignol J-P, Olivotto I, Rakovitch E, Gardner S, Sixel K, Beckham W, et al. A Multicenter Randomized Trial of Breast Intensity-Modulated Radiation Therapy to Reduce Acute Radiation Dermatitis. J Clin Oncol (2008) 26:2085–92. doi: 10.1200/JCO.2007.15.2488

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Donovan EM, Yarnold JR, Adams EJ, Morgan A, Warrington APJ, Evans PM. An Investigation Into Methods of IMRT Planning Applied to Breast Radiotherapy. Br J Radiol (2008) 81:311–22. doi: 10.1259/bjr/28583675

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Smith W, Menon G, Wolfe N, Ploquin N, Trotter T, Pudney D. IMRT for the Breast: A Comparison of Tangential Planning Techniques. Phys Med Biol (2010) 55:1231–41. doi: 10.1088/0031-9155/55/4/022

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Nichols GP, Fontenot JD, Gibbons JP, Sanders ME. Evaluation of Volumetric Modulated Arc Therapy for Postmastectomy Treatment. Radiat Oncol (2014) 9:66. doi: 10.1186/1748-717X-9-66

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Johansen S, Cozzi L, Olsen DR. A Planning Comparison of Dose Patterns in Organs at Risk and Predicted Risk for Radiation Induced Malignancy in the Contralateral Breast Following Radiation Therapy of Primary Breast Using Conventional, IMRT and Volumetric Modulated Arc Treatment Techniqu. Acta Oncol (Madr) (2009) 48:495–503. doi: 10.1080/02841860802657227

CrossRef Full Text | Google Scholar

9. Darby SC, Ewertz M, McGale P, Bennet AM, Blom-Goldman U, Brønnum D, et al. Risk of Ischemic Heart Disease in Women After Radiotherapy for Breast Cancer. N Engl J Med (2013) 368:987–98. doi: 10.1056/nejmoa1209825

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Johansen H, Kaae S, Johansen H, Kaae S, Jensen M-B, Mouridsen HT. Extended Radical Mastectomy Versus Simple Mastectomy Followed by Radiotherapy in Primary Breast Cancer. A Fifty-Year Follow-Up to the Copenhagen Breast Cancer Randomised Study. Acta Oncol (Madr) (2008) 47:633–8. doi: 10.1080/02841860801989753

CrossRef Full Text | Google Scholar

11. Penninkhof J, Spadola S, Breedveld S, Baaijens M, Lanconelli N, Heijmen B. Individualized Selection of Beam Angles and Treatment Isocenter in Tangential Breast Intensity Modulated Radiation Therapy. Int J Radiat Oncol Biol Phys (2017) 98:447–53. doi: 10.1016/j.ijrobp.2017.02.008

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Yan H, Yin F-F, Willett C. Evaluation of an Artificial Intelligence Guided Inverse Planning System: Clinical Case Study. Radiother Oncol (2007) 83:76–85. doi: 10.1016/j.radonc.2007.02.013

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Monz M, Küfer KH, Bortfeld TR, Thieke C. Pareto Navigation - Algorithmic Foundation of Interactive Multi-Criteria IMRT Planning. Phys Med Biol (2008) 53:985–98. doi: 10.1088/0031-9155/53/4/011

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Teichert K, Süss P, Serna JI, Monz M, Küfer KH, Thieke C. Comparative Analysis of Pareto Surfaces in Multi-Criteria IMRT Planning. Phys Med Biol (2011) 56:3669–84. doi: 10.1088/0031-9155/56/12/014

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Breedveld S, Storchi PRM, Voet PWJ, Heijmen BJM. ICycle: Integrated, Multicriterial Beam Angle, and Profile Optimization for Generation of Coplanar and Noncoplanar IMRT Plans. Med Phys (2012) 39:951–63. doi: 10.1118/1.3676689

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Xhaferllari I, Wong E, Bzdusek K, Lock M, Chen JZ. Automated IMRT Planning With Regional Optimization Using Planning Scripts. J Appl Clin Med Phys (2013) 14:176–91. doi: 10.1120/jacmp.v14i1.4052

CrossRef Full Text | Google Scholar

17. Fogliata A, Nicolini G, Bourgier C, Clivio A, De Rose F, Fenoglietto P, et al. Performance of a Knowledge-Based Model for Optimization of Volumetric Modulated Arc Therapy Plans for Single and Bilateral Breast Irradiation. PloS One (2015) 10:e0145137. doi: 10.1371/journal.pone.0145137

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Wu H, Jiang F, Yue H, Zhang H, Wang K, Zhang Y. Applying a RapidPlan Model Trained on a Technique and Orientation to Another: A Feasibility and Dosimetric Evaluation. Radiat Oncol (2016) 11:1–7. doi: 10.1186/s13014-016-0684-9

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Hussein M, South CP, Barry MA, Adams EJ, Jordan TJ, Stewart AJ, et al. Clinical Validation and Benchmarking of Knowledge-Based IMRT and VMAT Treatment Planning in Pelvic Anatomy. Radiother Oncol (2016) 120:473–9. doi: 10.1016/j.radonc.2016.06.022

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Delaney AR, Tol JP, Dahele M, Cuijpers J, Slotman BJ, Verbakel WFAR. Effect of Dosimetric Outliers on the Performance of a Commercial Knowledge-Based Planning Solution. Int J Radiat Oncol Biol Phys (2016) 94:469–77. doi: 10.1016/j.ijrobp.2015.11.011

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Moore KL. Automated Radiotherapy Treatment Planning. Semin Radiat Oncol (2019) 29:209–18. doi: 10.1016/j.semradonc.2019.02.003

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Ge Y, Wu QJ. Knowledge-Based Planning for Intensity-Modulated Radiation Therapy: A Review of Data-Driven Approaches. Med Phys (2019) 46:2760–75. doi: 10.1002/mp.13526

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Castriconi R, Fiorino C, Broggi S, Cozzarini C, Di Muzio N, Calandrino R, et al. Comprehensive Intra-Institution Stepping Validation of Knowledge-Based Models for Automatic Plan Optimization. Phys Med (2019) 57:231–7. doi: 10.1016/j.ejmp.2018.12.002

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Castriconi R, Cattaneo GM, Mangili P, Esposito PG, Broggi S, Cozzarini C, et al. Clinical Implementation of Knowledge-Based (KB) Automatic Plan Optimization for Helical Tomotherapy. Pract Radiat Oncol (2021) 11(2):E236–44. doi: 10.1016/j.prro.2020.09.012

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Castriconi R, Fiorino C, Passoni P, Broggi S, Di Muzio NG, Cattaneo GM, et al. Knowledge-Based Automatic Optimization of Adaptive Early-Regression-Guided VMAT for Rectal Cancer. Phys Med (2020) 70:58–64. doi: 10.1016/j.ejmp.2020.01.016

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Scaggion A, Fusella M, Roggio A, Bacco S, Pivato N, Rossato MA, et al. Reducing Inter- and Intra-Planner Variability in Radiotherapy Plan Output With a Commercial Knowledge-Based Planning Solution. Phys Med (2018) 53:86–93. doi: 10.1016/j.ejmp.2018.08.016

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Nawa K, Haga A, Nomoto A, Sarmiento RA, Shiraishi K, Yamashita H, et al. Evaluation of a Commercial Automatic Treatment Planning System for Prostate Cancers. Med Dosim (2017) 42:203–9. doi: 10.1016/j.meddos.2017.03.004

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Schubert C, Waletzko O, Weiss C, Voelzke D, Toperim S, Roeser A, et al. Intercenter Validation of a Knowledge Based Model for Automated Planning of Volumetric Modulated Arc Therapy for Prostate Cancer. The Experience of the German RapidPlan Consortium. PloS One (2017) 12(5):e0178034. doi: 10.1371/journal.pone.0178034

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Good D, Lo J, Lee WR, Wu QJ, Yin F-F, Das SK, et al. A Knowledge-Based Approach to Improving and Homogenizing Intensity Modulated Radiation Therapy Planning Quality Among Treatment Centers: An Example Application to Prostate Cancer Planning. Int J Radiat Oncol Biol Phys (2013) 87:176–81. doi: 10.1016/j.ijrobp.2013.03.015

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Zhao X, Kong D, Jozsef G, Chang J, Wong EK, Formenti SC, et al. Automated Beam Placement for Breast Radiotherapy Using a Support Vector Machine Based Algorithm. Med Phys (2012) 39:2536–43. doi: 10.1118/1.3700736

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Wang W, Purdie TG, Rahman M, Marshall A, Liu F-F, Fyles A. Rapid Automated Treatment Planning Process to Select Breast Cancer Patients for Active Breathing Control to Achieve Cardiac Dose Reduction. Int J Radiat Oncol Biol Phys (2012) 82:386–93. doi: 10.1016/j.ijrobp.2010.09.026

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Purdie TG, Dinniwell RE, Fyles A, Sharpe MB. Automation and Intensity Modulated Radiation Therapy for Individualized High-Quality Tangent Breast Treatment Plans. Int J Radiat Oncol Biol Phys (2014) 90:688–95. doi: 10.1016/j.ijrobp.2014.06.056

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Nicolini G, Clivio A, Fogliata A, Vanetti E, Cozzi L. Simultaneous Integrated Boost Radiotherapy for Bilateral Breast: A Treatment Planning and Dosimetric Comparison for Volumetric Modulated Arc and Fixed Field Intensity Modulated Therapy. Radiat Oncol (2009) 4:27. doi: 10.1186/1748-717X-4-27

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Fogliata A, Seppälä J, Reggiori G, Lobefalo F, Palumbo V, De Rose F, et al. Dosimetric Trade-Offs in Breast Treatment With VMAT Technique. Br J Radiol (2017) 90:20160701. doi: 10.1259/bjr.20160701

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Tsai PF, Lin S-M, Lee S-H, Yeh C-Y, Huang Y-T, Lee C-C, et al. The Feasibility Study of Using Multiple Partial Volumetric-Modulated Arcs Therapy in Early Stage Left-Sided Breast Cancer Patients. J Appl Clin Med Phys (2012) 13:62–73. doi: 10.1120/jacmp.v13i5.3806

CrossRef Full Text | Google Scholar

36. Esposito PG, Castriconi R, Mangili P, Fodor A, Pasetti M, Di Muzio NG, et al. Virtual Tangential-Fields Arc Therapy (ViTAT) for Whole Breast Irradiation: Technique Optimization and Validation. Phys Med (2020) 77:160–8. doi: 10.1016/j.ejmp.2020.08.011

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Lizondo M, Latorre-Musoll A, Ribas M, Carrasco P, Espinosa N, Coral A, et al. Pseudo Skin Flash on VMAT in Breast Radiotherapy: Optimization of Virtual Bolus Thickness and HU Values. Phys Med (2019) 63:56–62. doi: 10.1016/j.ejmp.2019.05.010

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Ricotti R, Ciardo D, Fattori G, Leonardi MR, Morra A, Dicuonzo S, et al. Intra-Fraction Respiratory Motion and Baseline Drift During Breast Helical Tomotherapy. Radiother Oncol (2017) 122:79–86. doi: 10.1016/j.radonc.2016.07.019

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Goddu SM, Yaddanapudi S, Pechenaya OL, Chaudhari SR, Klein EE, Khullar D, et al. Dosimetric Consequences of Uncorrected Setup Errors in Helical Tomotherapy Treatments of Breast-Cancer Patients. Radiother Oncol (2009) 93:64–70. doi: 10.1016/j.radonc.2009.07.013

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Cagni E, Botti A, Micera R, Galeandro M, Sghedoni R, Orlandi M, et al. Knowledge-Based Treatment Planning: An Inter-Technique and Inter-System Feasibility Study for Prostate Cancer. Phys Med (2017) 36:38–45. doi: 10.1016/j.ejmp.2017.03.002

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Alpuche Aviles JE, Cordero Marcos MI, Sasaki D, Sutherland K, Kane B, Kuusela E. Creation of Knowledge-Based Planning Models Intended for Large Scale Distribution: Minimizing the Effect of Outlier Plans. J Appl Clin Med Phys (2018) 19:215–26. doi: 10.1002/acm2.12322

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Cagni E, Botti A, Wang Y, Iori M, Petit SF, Heijmen BJM. Pareto-Optimal Plans as Ground Truth for Validation of a Commercial System for Knowledge-Based DVH-Prediction. Phys Med (2018) 55:98–106. doi: 10.1016/j.ejmp.2018.11.002

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Rinaldin G, Perna L, Agnello G, Pallazzi G, Cattaneo GM, Fiorino C, et al. Quality Assurance of Rapid Arc Treatments: Performances and Pre-Clinical Verifications of a Planar Detector (MapCHECK2). Phys Med (2014) 30:184–90. doi: 10.1016/j.ejmp.2013.05.004

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Grantzau T, Thomsen MS, Væth M, Overgaard J. Risk of Second Primary Lung Cancer in Women After Radiotherapy for Breast Cancer. Radiother Oncol (2014) 111:366–73. doi: 10.1016/j.radonc.2014.05.004

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Virén T, Heikkilä J, Myllyoja K, Koskela K, Lahtinen T, Seppälä J. Tangential Volumetric Modulated Arc Therapy Technique for Left-Sided Breast Cancer Radiotherapy. Radiat Oncol (2015) 10:79. doi: 10.1186/s13014-015-0392-x

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Yu AF, Jones LW. Breast Cancer Treatment-Associated Cardiovascular Toxicity and Effects of Exercise Countermeasures. Cardiooncology (2016) 2:1. doi: 10.1186/s40959-016-0011-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: breast cancer, radiation oncology, automation, plan optimization, tangential field, knowledge-based

Citation: Castriconi R, Esposito PG, Tudda A, Mangili P, Broggi S, Fodor A, Deantoni CL, Longobardi B, Pasetti M, Perna L, del Vecchio A, Di Muzio NG and Fiorino C (2021) Replacing Manual Planning of Whole Breast Irradiation With Knowledge-Based Automatic Optimization by Virtual Tangential-Fields Arc Therapy. Front. Oncol. 11:712423. doi: 10.3389/fonc.2021.712423

Received: 20 May 2021; Accepted: 02 August 2021;
Published: 24 August 2021.

Edited by:

Savino Cilla, Gemalli Molise Hospital, Italy

Reviewed by:

Raquel Bar-Deroma, Rambam Health Care Campus, Israel
Tsair-Fwu Lee, National Kaohsiung University of Science and Technology, Taiwan

Copyright © 2021 Castriconi, Esposito, Tudda, Mangili, Broggi, Fodor, Deantoni, Longobardi, Pasetti, Perna, del Vecchio, Di Muzio and Fiorino. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Claudio Fiorino, Zmlvcmluby5jbGF1ZGlvQGhzci5pdA==

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.