- Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA, United States
Radiotherapy (RT) doses to cardiac substructures from the definitive treatment of locally advanced non-small cell lung cancers (NSCLC) have been linked to post-RT cardiac toxicities. With modern treatment delivery techniques, it is possible to focus radiation doses to the planning target volume while reducing cardiac substructure doses. However, it is often challenging to design such treatment plans due to complex tradeoffs involving numerous cardiac substructures. Here, we built a cardiac-substructure-based knowledge-based planning (CS-KBP) model and retrospectively evaluated its performance against a cardiac-based KBP (C-KBP) model and manually optimized patient treatment plans. CS-KBP/C-KBP models were built with 27 previously-treated plans that preferentially spare the heart. While the C-KBP training plans were created with whole heart structures, the CS-KBP model training plans each have 15 cardiac substructures (coronary arteries, valves, great vessels, and chambers of the heart). CS-KBP training plans reflect cardiac-substructure sparing preferences. We evaluated both models on 28 additional patients. Three sets of treatment plans were compared: (1) manually optimized, (2) C-KBP model-generated, and (3) CS-KBP model-generated. Plans were normalized to receive the prescribed dose to at least 95% of the PTV. A two-tailed paired-sample t-test was performed for clinically relevant dose-volume metrics to evaluate the performance of the CS-KBP model against the C-KBP model and clinical plans, respectively. Overall results show significantly improved cardiac substructure sparing by CS-KBP in comparison to C-KBP and the clinical plans. For instance, the average left anterior descending artery volume receiving 15 Gy (V15 Gy) was significantly lower (p < 0.01) for CS-KBP (0.69 ± 1.57 cc) compared to the clinical plans (1.23 ± 1.76 cc) and C-KBP plans (1.05 ± 1.68 cc). In conclusion, the CS-KBP model significantly improved cardiac-substructure sparing without exceeding the tolerances of other OARs or compromising PTV coverage.
Introduction
Radiation therapy (RT) for patients receiving definitive treatment for locally advanced non-small cell lung cancers (NSCLC) has been linked to post-RT cardiac toxicities (1). A number of promising phase II trials reported increased tumor control through radiation dose escalation via stereotactic fractionations (2–4). The Radiation Therapy Oncology Group (RTOG) 0617 phase III dose escalation trial compared 74 Gy with the standard 60 Gy dose with concurrent chemotherapy (5). In this study, the higher prescription dose resulted in worse overall survival. Furthermore, cardiac volumes receiving 5 Gy and 30 Gy were found to be linked with higher death rates. In a multivariate analysis, an increased radiation dose to the heart was independently linked with worse overall survival (6). More recent studies have reported the link between radiation doses to cardiac substructures and major adverse cardiac events, including mortality (1, 7). Atkins et al. reported coronary artery disease progression and the risk of ischemic events due to radiation exposure to arterial vessels, especially the left anterior descending artery (LADA) (7). Skytta et al. reported an increase of serum troponin T>30% as an endpoint to greater than 10% of LADA receiving 15 Gy (8). McWilliam et al. identified the base of the heart as a dose-sensitive region linked with poorer survival in lung cancer patients. In this study, patients receiving a mean dose of more than 8.5 Gy to this region were found to have significantly poorer survival (9). In another study, the maximum dose to the LADA of greater than 53 Gy was linked to increased risks of cardiac morbidity (10). The maximum dose to the right atrium, right coronary artery, and ascending aorta were found to be the most important factors associated with survival (11). Gore et al. reported that the pericardium mean dose and atria/ventricles volume covered by 45 Gy have detrimental effects on overall survival (12).
Clinical treatment planning is a trial-and-error process: the planner creates a treatment plan by iteratively updating optimization constraints in order to meet the target coverage while minimizing the dose to surrounding OARs. Despite best efforts from the planners, many clinical plans can be further improved in terms of OAR sparing. Such sub-optimal dose sparing can lead to increased doses to cardiac substructures and surrounding critical organs. In the past decade, knowledge-based planning (KBP) methods have been proposed to leverage knowledge embedded in prior high-quality plans to generate treatment plans for new patients (13, 14). Briefly, RapidPlan™ KBP estimates achievable DVHs for various structures included into the model, which are then translated into optimization objectives. The algorithm is divided into two parts, the model configuration and DVH estimation components. In model configuration, data extraction phase prepares the dataset for model training by extracting different geometric features of each structure. The model training phase produces a DVH estimation model for each structure. This is followed by generating estimated DVHs and optimization objectives through DVH estimation components. Previously, Harms et al. utilized RapidPlan™ KBP to improve cardiac sparing by incorporating the whole heart structure into the model training phase (15). While including the heart structure in the model would result in useful DVH predictions for the whole heart, it does not provide spatial information regarding dose distribution within the heart. We hypothesized that including cardiac substructures in the training phase of KBP would provide further dose reduction to critical cardiac substructures due to the additional spatial and geometrical features encoded into to each cardiac structure models. Therefore, the focus of this study was to build an improved KBP containing the features of each substructure (CS-KBP), in addition to other relevant OARs, to further enhance cardiac substructure sparing without compromising the target coverage. To assess the performance of a trained CS-KBP model, we used a validation dataset consisting of 28 additional patients that were not used in training the CS-KBP model.
Materials and methods
Patient dataset
Two separate datasets, training and validation, consisting of patients previously treated with 60 Gy in 30 fractions for NSCLC, were retrospectively collected for this study. The training cohort consisted of 27 patients who were treated with 60 Gy in 30 fractions for non-small cell lung cancer; 26 of 27 patients were included in a previous study from our institution (16), whereas the validation cohort consisted of 28 additional patients that were not included in the training phase. Patient characteristics of training and validation datasets have been listed in Table S1. For each patient, 15 cardiac substructures were contoured on free-breathing CT scans by resident physicians based on guidelines provided by Feng et al. (17) and reviewed by an attending physician. These substructures include the coronary arteries [LADA, left circumflex (LCFLX), left main coronary artery (LMCA), right coronary artery (RCA)], the great vessels [ascending aorta (AA), pulmonary artery (PA), superior vena cava (SVC)], the valves [atrial (AV), mitral (MV), pulmonary (PV), and tricuspid (TV)], and the chambers of the heart [left atrium (LA), right atrium (RA), right ventricle (RV), and left ventricle (LV)] The target volumes consisted of the clinical target volume (CTV) and planning target volume (PTV). These cardiac substructure contours are derived from a previous study in our group (15). However, datasets are randomly re-selected for training and validation cohorts in this study. The model constraints have also been re-tuned to improve model robustness.
RapidPlan knowledge-based planning model
RapidPlan is a knowledge-based planning (KBP) module integrated within the Eclipse treatment planning systems (Varian medical Systems, Palo Alto, CA). RapidPlan has been widely studied and validated for various treatment sites (13, 14, 18–23). Briefly, RapidPlan utilizes anatomical and geometrical features from previously optimized plans to build a dose prediction model. The trained model is then used to predict dose volume histograms (DVHs) of various OARs and target structures for a new case of the same disease site. Detailed information on knowledge-based planning can be found in previous publications (24–26). To build a cardiac-substructure-specific KBP model (CS-KBP), we first optimized the 27 patients’ plans with their cardiac substructure preferences. To optimize plans with cardiac substructure preference, we used a line constraint for each cardiac substructure in the optimization template. Each line constraint is a collection of dose-volume point constraints placed based on RapidPlan DVH predictions.
To train the CS-KBP model, these CS-optimized plans were then used as an input along with structures including CTV, PTV, cardiac substructures, and additional organs at risk (OAR) structures (ipsilateral and contralateral lung, whole lungs, lungs cropped out of CTV, esophagus, and spinal cord). For validation, the trained CS-KBP model was then used to retrospectively generate a separate cohort of the 28 NSCLC patients’ treatment plans. In the treatment planning system (TPS), the DVH estimation tool was used to initialize the optimization to generate plans without any intervention.
Dosimetric results of this validation cohort from CS-KBP were compared against two other set of plans: clinically optimized plans and the plans generated from cardiac-based KBP model (C-KBP). Clinically optimized plans were optimized by our dosimetry team to meet desired clinical dosimetric goals for relevant OARs including the whole heart as a single structure. The C-KBP model was trained with cardiac-optimized treatment plans with the whole heart and other clinically relevant OARs as the structural input. The plan parameters including number of arcs and control points were held constant for each plan in three datasets. Table 1 shows the summary of the three different set of plans with corresponding preferences and model inputs.
Plan evaluations
The treatment plans generated by the CS-KBP model were compared against clinically optimized plans and plans generated by the C-KBP model, respectively. All dose calculations were performed in the Eclipse treatment planning software by using the Anisotropic Analytical Algorithm, version 15.6.05. To perform the dose calculation, contours were transferred through rigid registration to the average 4D CT. For dosimetric evaluations of OARs structures, clinically relevant dose volume constraints were compared across the three optimization methods. In addition to the maximum and mean doses to each cardiac substructure, we also investigated the changes in cardiac substructures’ specific dose volume metrics that have been shown to be associated with post-RT cardiac toxicities (8, 11, 12, 27, 28). We performed a paired two-tailed t-test to determine whether there is statistical evidence that the mean difference between paired observations is significantly different from zero. The statistical analysis was performed using Matlab’s statistical toolkit (Matlab, Mathworks Inc.). We compared the CS-KBP vs. C-KBP and CS-KBP vs. clinical plans with 95% and 99% confidence intervals for each dosimetric parameter. This test was performed with statistical significance set at p < 0.05. Finally, it is important to note that all plans were normalized to cover 100% PTV by ≥95% of the prescribed dose (57 Gy).
Results
We performed qualitative (Figure 1) and quantitative (Figures 2, 3, 4; S1, 2) comparisons among the treatment plans generated by three methods. Figure 1 shows the dose distribution comparison among three different methods for two sample cases, each with the same CT scan slice, from the validation cohort. Clinically optimized plans are more conformal than the plans produced by C-KBP and CS-KBP for both cases, but at the expense of higher doses to coronary arteries such as LADA, LCFLX, and LMCA. While the C-KBP plan leads to better substructure sparing than clinical plan, the left circumflex coronary artery receives up to 33 Gy dose. CS-KBP plans spare cardiac substructures to a large extent with a tolerable tradeoff in plan conformity. Table 2 shows the dosimetric comparisons among the three techniques for the two cases shown in Figure 1.
Figure 1 Dose distribution comparison between CS-KBP, C-KBP, and manually optimized plans for two sample cases from the validation cohort.
Figure 2 Box plots for clinically relevant dose volume metrics for LADA (A, B), LCFLX (C), LMCA (D), and RCA (E) over 28 treatment plans generated by clinically optimized (orange), CKBP (purple), and CS-KBP (yellow). The central mark indicates the median, and the bottom and top edges of the box indicate 25th and 75th percentile, respectively. The whiskers extend to the most extreme data points not considered outliers, and the outliers are plotted individually using the ‘+’ marker symbol. p < 0.05 is denoted as ‘*’; p < 0.01 is denoted as ‘**’; NS = Not Significant.
Figure 3 Box plots for clinically relevant dose volume metrics for LA (A, B), LV (C, D), RA (E), and RV (F) over 28 treatment plans generated by clinically optimized (orange), CKBP (purple), and CS-KBP (yellow). The central mark indicates the median, and the bottom and top edges of the box indicate 25th and 75th percentile, respectively. The whiskers extend to the most extreme data points not considered outliers, and the outliers are plotted individually using the ‘+’ marker symbol. p < 0.05 is denoted as ‘*’; p < 0.01 is denoted as ‘**’; NS = Not Significant.
Figure 4 Box plots for volume of PA covered by 40 Gy isodose line (A) , maximum and mean dose to PA (B) , AA (C) , and SVC (D) over 28 treatment plans generated by clinically optimized (orange), C-KBP (purple), and CS-KBP (yellow). The central mark indicates the median, and the bottom and top edges of the box indicate 25th and 75th percentile, respectively. The whiskers extend to the most extreme data points not considered outliers, and the outliers are plotted individually using the ‘+’ marker symbol. p < 0.05 is denoted as ‘*’; p < 0.01 is denoted as ‘**’; NS = Not Significant.
Table 2 Comparison of dosimetric parameters among three techniques for two cases shown in Figure 1.
Figure 2A shows that the CS-KBP model significantly reduced the volume covered by 15, 20, 30, 40, and 45 Gy isodose volumes for LADA compared to two competing methods (Figure 2B). The CS-KBP model also produced plans with significantly lower mean and maximum doses (9.02 Gy and 20.04 Gy) than the plans generated by the C-KBP model (12.88 Gy and 25.94 Gy) and manual optimization (16.05 Gy and 29.65 Gy). A similar trend is observed for mean doses to LCFLX (12.74 Gy vs. 17.68 Gy and 20.68 Gy), LMCA (11.36 Gy vs. 20.44 Gy and 26.74 Gy), and RCA (11.36 Gy vs. 20.44 Gy and 26.74 Gy) (Figures 2C, D). In general, CS-KBP outperformed competing methods for reducing doses to coronary arteries.
Figure 3 outlines the doses to the chambers of the heart from the three techniques. While there was no significant differences in maximum dose and V63 Gy for LA among the three techniques, CS-KBP produced plans with the significant sparing of LA measured in terms of mean dose (13.34 Gy vs. 16.22 Gy and 17.79 Gy) and V45 Gy (6.74 cc vs. 7.88 cc and 11.13 cc) in comparison to treatment plans generated by the C-KBP model and manual optimizations (Figures 3A, B). For LV (3.55 Gy vs. 4.72 Gy and 6.52 Gy), RV (2.80 Gy vs. 4.11 Gy and 5.95 Gy) and RA (3.95 Gy vs. 5.59 Gy and 6.55 Gy), CS-KBP produced significant sparing across all dose volume metrics (Figures 3C-F).
Figure 4 summarizes the results for the great vessels across all three methods. Figure 4 shows a significant difference between CS-KBP and clinical plans for AA. However, Figure 4 shows that all three methods performed comparably (p > 0.05) for maximum doses to PA and SVC. However, CS-KBP plans resulted in significantly lower mean doses to all three great vessels compared to two competing methods. As shown in Figure S2, CS-KBP also produced plans with significant sparing of all four cardiac valves compared to C-KBP and manually optimized plans. For whole heart structures, mean doses were 14.87 ± 5.01 Gy, 10.97 ± 3.41 Gy, and 8.95 ± 3.89 Gy for clinical, C-KBP, and CS-KBP plans, respectively. In terms of doses to the lungs, spinal cord, and esophagus, manually optimized plans largely performed better than both KBP models (Figures S1, S3). However, it is important to note that there was no violation of the clinically important dose volume constraints for these structures by treatment plans generated by the C-KBP and CS-KBP model.
Discussion
A large number of studies in the literature have shown a correlation between post-radiotherapy cardiac toxicities and cardiac substructure irradiation (1, 3, 6, 10, 11, 29, 30). As a result, recent treatment planning efforts for thoracic cancers have been placed on sparing cardiac substructures surrounding the PTV compared to the traditional approach of sparing the heart as a single organ. Given a large number of substructures, it can be time-consuming to reach a desired treatment plan while effectively sparing substructures and other OARs. In general, it is challenging to build a KBP model for lung cancer due to variations in the location, shape, size, and orientation of PTV with respect to cardiac substructures and other critical OARs. This study developed a CS-KBP model and demonstrated its utility in the significant sparing of cardiac substructures for the VMAT of patients treated with definitive RT for locally advanced NSCLC.
For thoracic cancers, original QUANTEC recommendations are to limit the volume of the heart receiving ≥25 Gy less than 10% to keep the risk of cardiac mortality under 1% (31). Most current protocols now, however, aim to keep the mean heart dose below 20 Gy due to its correlation with a significantly higher rate of cardiac events (32). Recent works recommended substructure-specific dose volume constraints for various endpoints. For instance, Dmax of 19.5 Gy to RCA has been associated with worse overall survival (11). The CS-KBP model (Dmax = 8.1 ± 10.2 Gy) was able to reduce doses to RCA by 5.1 Gy and 10.1 Gy, respectively, compared to plans generated by C-KBP and manual optimization. For LADA, Dmean ≥23.8 Gy and V15 Gy ≥10.0% has been associated with an increased percentage of serum troponin-T, which is an indicator of myocardial infarction. The CS-KBP model was able to generate plans with significantly lower mean doses to LADA compared to two competing methods (9.02 Gy vs. 12.88 Gy and 16.05 Gy), and also a lower volume of LADA covered with 15 Gy isodose than two competing methods (0.69 cc vs. 1.04 cc and 1.23 cc). In terms of great vessels, Dmean ≥8.5 Gy to AA, PA, and SVC has been linked to worse overall survival (9). While all three methods met these constraints, the CS-KBP model significantly reduced radiation doses to great vessels (Figure 4). However, it is important to note that great vessels of the heart were not entirely contoured within the whole heart structure for clinical and C-KBP plans, which is one limitation of this study. For LV, the volume covered by the 5 Gy (V5) isodose has been linked to overall survival and acute coronary events (12, 33). The CS-KBP model-produced plans significantly improved this constraint compared to the C-KBP and optimized plans (Figure 3). For LV, Dmean ≥ 6.7 Gy has been shown to increase the percentage of serum troponin (8). While all three techniques satisfied this constraint, the CS-KBP model resulted in lower Dmean to LV than C-KBP and manual optimization (3.55 Gy vs. 6.52 Gy and 4.72 Gy). For LV, V30 Gy is also linked with symptomatic cardiac events (34). Plans generated by the CS-KBP model resulted in lower values of V30 Gy compared to two competing methods (1.54 cc vs. 3.46 cc and 9.44 cc). In terms of atrial valve (AV), Dmax ≥19.5 Gy has been associated with worse overall survival (11). While all three methods met this Dmax constraint for AV, the CS-KBP model produced plans with the significantly lower Dmax and Dmean to the valves of the heart (Figure S2). Pericardial effusion is commonly reported as one of the cardiac toxicities post thoracic radiotherapy (35). The occurrence of pericardial effusion of any grade has been reported to be higher in patients with a mean heart dose >23.45 Gy (2). While the mean heart dose constraint is met with all three techniques in this study, the proposed CS-KBP outperforms the other two techniques (Clinical: 10.68 Gy; C-KBP: 8.57 Gy; CS-KBP: 6.36 Gy).
Sparing of cardiac substructures comes with some necessary tradeoffs. Radiation doses to 65% of the lungs were significantly lower in manually optimized plans than both KBP plans. In terms of mean doses to lungs, plans generated by C-KBP (14.23 Gy) resulted in slightly lower doses than CS-KBP (14.85 Gy) and manual optimization (14.31 Gy). In terms of the dose received by 0.03 cc of the spinal cord, manually optimized plans resulted in the lowest dose (28.06 Gy) compared to C-KBP (30.77 Gy) and CS-KBP (34.26 Gy). Nonetheless, all three techniques satisfied the clinical constraints of keeping the maximum cord dose below 40 Gy. Plans generated with CS-KBP had higher plan complexity (630MU) than the ones generated via the C-KBP model (565MU) and manual optimization (486 MU). This tradeoff among different structures and in increased plan complexity is consistent with that of Harms et al., who found that their clinical RapidPlan™ (RP) model significantly outperformed the cardiac sparing RP model with increased plan complexity by the cardiac sparing RP model in terms of the number of monitor units (15). While CS-KBP led to the lowest overall doses to the cardiac substructures, the two competing methods also achieved the dose-volume constraints. The clinical impact of recently introduced cardiac substructure-specific dose constraints is actively being studied based on randomized trials and from real-world datasets on lung cancer (36). In contrast, there are standardized dose volume constraints for other OARs (i.e., lung and spinal cord). This study demonstrates the ability of CS-KBP to meet these standard dose volume constraints of other OARs while significantly reducing the doses to cardiac substructures. Though the direct clinical benefits from reducing the dose to cardiac substructures is yet to be validated with patient outcome data, lower doses to these cardiac substructures are expected to translate to reduced probability of radiation-induced cardiac toxicities post radiotherapy.
A limitation of the present study is that all the manually optimized plans were generated by a few dosimetrists and all thoracic OARs were contoured by a single observer. Another limitation of this study is the use of conventional CT for contouring cardiac substructures, which may introduce uncertainty due to the poor discernibility of substructures on CT images. This limitation can potentially be addressed by previous studies that leveraged the enhanced soft tissue contrast of MRI with conventional non-contrast CT to develop a deep learning-based framework (37) and a hybrid MR/CT segmentation atlas-based framework (38) for cardiac substructure segmentation. Though training the sample size was relatively small, it was within a range that had been shown to be adequate for a KBP model for thoracic radiation therapy (39). Another minor limitation of this study is a slight mismatch between the contours and the plans. Cardiac substructures were contoured on free-breathing scans for higher spatial resolution and then transferred via rigid registration to the average 4D CT for treatment planning. Although all patients in the cohort were treated with a free breathing technique, slight differences between the free breathing scans and the 4D average scans can be expected. Nonetheless, since all three sets of plans were generated on the same datasets, the results reflect the true differences among the three techniques. Following preliminary cardiac substructure segmentations by a physician resident, all cardiac substructures were verified by an attending physician. While we did not perform a robust evaluation for the quality of segmentation due to the nature of this study; different strategies such as percent overlap of observer volumes, a robust atlas-based method, etc. can be utilized clinically to ensure high quality segmentation of cardiac substructures.
All treatment plans were generated within one institution, lacking variations in terms of contouring and dosimetric protocols from multi-institutions. A validation cohort from multi-institutional data is necessary to further ensure the robustness of this model and also consider differences between planners and dosimetry protocols. Nonetheless, there was a large variation in the size of the PTV volumes in the validation cohort (454.05 ± 249.42 cc) of this study, demonstrating the feasibility of the proposed model for variable tumor sizes. In this study, all the structures were manually contoured, which may be time consuming. Various deep learning frameworks have been introduced to increase treatment planning efficiency (40). The impact of auto-contouring versus manual contouring on the prediction model will be investigated in the future studies. While KBP models have been built for various treatment sites, it has not been widely adopted for the purpose of cardiac substructure sparing (13). Overall results demonstrated the ability of proposed CS-KBP model to achieve significant cardiac substructure sparing compared to the traditional method of manual optimization as well as the KBP model based on cardiac structure alone (C-KBP).
Conclusion
The proposed approach can aid future researchers to implement the knowledge-based planning module to effectively spare cardiac substructures while satisfying the clinical constraints of other OARs and target coverage and also further enhance treatment planning efficiency. This model may offer reduced planning time, improve plan quality, and should be studied in the prospective setting as a way to improve outcomes in patients with locally advanced lung cancer.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Author contributions
SM: First authorship. JW, JR, YL, TL, JB, KH, XY: Co-authorship. JZ: Senior authorship. All authors contributed to the article and approved the submitted version.
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.2022.1055428/full#supplementary-material
References
1. Venkatesulu BP, Mahadevan LS, Aliru ML, Yang X, Bodd MH, Singh PK, et al. Radiation-induced endothelial vascular injury: a review of possible mechanisms. JACC: Basic to Trans Sci (2018) 3(4):563–72. doi: 10.1016/j.jacbts.2018.01.014
2. Bradley JD, Moughan J, Graham MV, Byhardt R, Govindan R, Fowler J, et al. A phase I/II radiation dose escalation study with concurrent chemotherapy for patients with inoperable stages I to III non-small-cell lung cancer: Phase I results of RTOG 0117. Int J Radiat Oncol Biol Phys (2010) 77(2):367–72. doi: 10.1016/j.ijrobp.2009.04.029
3. Stinchcombe TE, Lee CB, Moore DT, Rivera MP, Halle J, Limentani S, et al. Long-term follow-up of a phase I/II trial of dose escalating three-dimensional conformal thoracic radiation therapy with induction and concurrent carboplatin and paclitaxel in unresectable stage IIIA/B non-small cell lung cancer. J Thorac Oncol (2008) 3(11):1279–85. doi: 10.1097/JTO.0b013e31818b1971
4. Socinski MA, Blackstock AW, Bogart JA, Wang X, Munley M, Rosenman J, et al. Randomized phase II trial of induction chemotherapy followed by concurrent chemotherapy and dose-escalated thoracic conformal radiotherapy (74 gy) in stage III non–small-cell lung cancer: CALGB 30105. J Clin Oncol (2008) 26(15):2457–63. doi: 10.1200/JCO.2007.14.7371
5. Bradley JD, Paulus R, Komaki R, Masters G, Blumenschein G, Schild S, 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. doi: 10.1016/S1470-2045(14)71207-0
6. Faivre-Finn C. Dose escalation in lung cancer: Have we gone full circle? Lancet Oncol (2015) 16(2):125–7. doi: 10.1016/S1470-2045(15)70001-X
7. Atkins KM, Chaunzwa TL, Lamba N, Bitterman DS, Rawal B, Bredfeldt J, et al. Association of left anterior descending coronary artery radiation dose with major adverse cardiac events and mortality in patients with non–small cell lung cancer. JAMA Oncol (2021) 7(2):206–19. doi: 10.1001/jamaoncol.2020.6332
8. Skyttä T, Tuohinen S, Boman E, Virtanen V, Raatikainen P, Kellokumpu-Lehtinen P-L, et al. Troponin T-release associates with cardiac radiation doses during adjuvant left-sided breast cancer radiotherapy. Radiat Oncol (2015) 10(1):1–8. doi: 10.1186/s13014-015-0436-2
9. McWilliam A, Kennedy J, Hodgson C, Osorio EV, Faivre-Finn C, Van Herk M, et al. Radiation dose to heart base linked with poorer survival in lung cancer patients. Eur J Cancer (2017) 85:106–13. doi: 10.1016/j.ejca.2017.07.053
10. Poitevin-Chacón A, Chávez-Nogueda J, Prudencio RR, Fernández AC, Laguna AR, Linares J, et al. Dosimetry of the left anterior descending coronary artery in left breast cancer patients treated with postoperative external radiotherapy. Rep Pract Oncol Radiother (2018) 23(2):91–6. doi: 10.1016/j.rpor.2018.01.003
11. McWilliam A, Khalifa J, Osorio EV, Banfill K, Abravan A, Faivre-Finn C, et al. Novel methodology to investigate the effect of radiation dose to heart substructures on overall survival. Int J Radiat Oncol Biol Phys (2020) 108(4):1073–81. doi: 10.1016/j.ijrobp.2020.06.031
12. Gore E, Hu C, Ad VB, Robinson C, Wheatley M, Bogart J, et al. Impact of incidental cardiac radiation on cardiopulmonary toxicity and survival for locally advanced non-small cell lung cancer: Reanalysis of NRG oncology/RTOG 0617 with centrally contoured cardiac structures. Int J Radiat Oncol Biol Phys (2016) 96(2):S129–30. doi: 10.1016/j.ijrobp.2016.06.316
13. Momin S, Fu Y, Lei Y, Ropley J, Bradley JD, Curran WJ, et al. Knowledge-based radiation treatment planning: A data-driven method survey. J App Clin Med Phys (2021), 22(8):16–44. doi: 10.1002/acm2.13337
14. Ge Y, Wu QJ. Knowledge-based planning for intensity-modulated radiation therapy: a review of data-driven approaches. Med Phys (2019) 46(6):2760–75. doi: 10.1002/mp.13526
15. Harms J, Zhang J, Kayode O, Wolf J, Tian S, McCall N, et al. Implementation of a knowledge-based treatment planning model for cardiac-sparing lung radiation therapy. Adv Radiat Oncol (2021) 6(6):100745. doi: 10.1016/j.adro.2021.100745
16. Ferris MJ, Martin KS, Switchenko JM, Kayode OA, Wolf J, Dang Q, et al. Sparing cardiac substructures with optimized volumetric modulated arc therapy and intensity modulated proton therapy in thoracic radiation for locally advanced non-small cell lung cancer. Pract Radiat Oncol (2019) 9(5):e473–81. doi: 10.1016/j.prro.2019.04.013
17. Feng M, Moran JM, Koelling T, Chughtai A, Chan JL, Freedman L, et al. Development and validation of a heart atlas to study cardiac exposure to radiation following treatment for breast cancer. Int J Rad Onc Biol Phys (2011) 79(1):10–8. doi: 10.1016/j.ijrobp.2009.10.058
18. Tol JP, Delaney AR, Dahele M, Slotman BJ, Verbakel WFAR. Evaluation of a knowledge-based planning solution for head and neck cancer. Int J Radiat Oncol Biol Phys (2015) 91(3):612–20. doi: 10.1016/j.ijrobp.2014.11.014
19. Cornell M, Kaderka R, Hild SJ, Ray XJ, Murphy JD, Atwood TF, et al. Noninferiority study of automated knowledge-based planning versus human-driven optimization across multiple disease sites. Int J Radiat Oncol Biol Phys (2020) 106(2):430–9. doi: 10.1016/j.ijrobp.2019.10.036
20. Chang ATY, Hung AWM, Cheung FWK, Lee MCH, Chan OSH, Philips H, et al. Comparison of planning quality and efficiency between conventional and knowledge-based algorithms in nasopharyngeal cancer patients using intensity modulated radiation therapy. Int J Radiat Oncol Biol Phys (2016) 95(3):981–90. doi: 10.1016/j.ijrobp.2016.02.017
21. Wu H, Jiang F, Yue H, Li S, Zhang Y. A dosimetric evaluation of knowledge-based VMAT planning with simultaneous integrated boosting for rectal cancer patients. J Appl Clin Med Phys (2016) 17(6):78–85. doi: 10.1120/jacmp.v17i6.6410
22. Yu G, Li Y, Feng Z, Tao C, Yu Z, Li B, et al. Knowledge-based IMRT planning for individual liver cancer patients using a novel specific model. Radiat Oncol (2018) 13(1):52. doi: 10.1186/s13014-018-0996-z
23. Kubo K, Monzen H, Ishii K, Tamura M, Kawamorita R, Sumida I, et al. Dosimetric comparison of RapidPlan and manually optimized plans in volumetric modulated arc therapy for prostate cancer. Phys Med (2017) 44:199–204. doi: 10.1016/j.ejmp.2017.06.026
24. Zhang J, Ge Y, Wu QJ. Knowledge-based treatment planning. In: Machine and deep learning in oncology, medical physics and radiology. Springer: Cham (2022). p. 307–34.
25. Fogliata A, Cozzi L, Reggiori G, Stravato A, Lobefalo F, Franzese C, et al. RapidPlan knowledge based planning: Iterative learning process and model ability to steer planning strategies. Radiat Oncol (2019) 14(1):1–12. doi: 10.1186/s13014-019-1403-0
26. DVH estimation algorithm for RapidPlan. In: Eclipse photon and electron algorithms reference guide. Palo Alto, CA: Varian Medical Systems.
27. Piroth MD, Baumann R, Budach W, Dunst J, Feyer P, Fietkau R, et al. Heart toxicity from breast cancer radiotherapy. Strahlentherapie und Onkol (2019) 195(1):1–12. doi: 10.1007/s00066-018-1378-z
28. Ma J-T, Sun L, Sun X, Xiong Z-C, Liu Y, Zhang S-L, et al. Is pulmonary artery a dose-limiting organ at risk in non-small cell lung cancer patients treated with definitive radiotherapy? Radiat Oncol (2017) 12(1):1–10. doi: 10.1186/s13014-017-0772-5
29. Dess RT, Sun Y, Matuszak MM, Sun G, Soni PD, Bazzi L, et al. Cardiac events after radiation therapy: Combined analysis of prospective multicenter trials for locally advanced non–small-cell lung cancer. J Clin Oncol (2017) 35(13):1395. doi: 10.1200/JCO.2016.71.6142
30. Wang K, Pearlstein KA, Patchett ND, Deal AM, Mavroidis P, Jensen BC, 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. doi: 10.1016/j.radonc.2017.10.001
31. Gagliardi G, Constine LS, Moiseenko V, Correa C, Pierce LJ, Allen AM, et al. Radiation dose–volume effects in the heart. Int J Radiat Oncol Biol Phys (2010) 76(3):S77–85. doi: 10.1016/j.ijrobp.2009.04.093
32. Vojtíšek R. Cardiac toxicity of lung cancer radiotherapy. Rep Pract Oncol Radiother (2020) 25(1):13–9. doi: 10.1016/j.rpor.2019.10.007
33. van den Bogaard VA, Khalifa J, Osorio EV, Banfill K, Abravan A, Faivre-Finn C, et al. Validation and modification of a prediction model for acute cardiac events in patients with breast cancer treated with radiotherapy based on three-dimensional dose distributions to cardiac substructures. J Clin Oncol (2017) 35(11):1171. doi: 10.1200/JCO.2016.69.8480
34. Wang K, Ta BD, van der Schaaf A, Bouma AB, Middag AM, Bantema-Joppe EJ, 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. doi: 10.1200/JCO.2016.70.0229
35. Pao T-H, Belosi F, Clivio A, Navarria P, Nicolini G, Scorsetti M, et al. Pericardial effusion after definitive concurrent chemotherapy and intensity modulated radiotherapy for esophageal cancer. Radiat Oncol (2020) 15(1):1–10. doi: 10.1186/s13014-020-01498-3
36. Abravan A, Eblan MJ, Deal AM, Lipner M, Zagar TM, Wang Y, et al. Role of real-world data in assessing cardiac toxicity after lung cancer radiotherapy. Front Oncol (2022) 12. doi: 10.3389/fonc.2022.934369
37. Morris ED, Chang W-L, Chiang N-J, Lin C-Y, Lai W-W, Tseng Y-L, et al. Cardiac substructure segmentation with deep learning for improved cardiac sparing. Med Phys (2020) 47(2):576–86. doi: 10.1002/mp.13940
38. Morris ED, Price G, Banfill K, Marchant T, Craddock M, Wood J, et al. Cardiac substructure segmentation and dosimetry using a novel hybrid magnetic resonance and computed tomography cardiac atlas. Int J Radiat Oncol Biol Phys (2019) 103(4):985–93. doi: 10.1016/j.ijrobp.2018.11.025
39. Fogliata A, Ghanem AI, Dong M, Pantelic MV, Walker EM, Glide‐Hurst CK, et al. On the pre-clinical validation of a commercial model-based optimisation engine: Application to volumetric modulated arc therapy for patients with lung or prostate cancer. Radiother Oncol (2014) 113(3):385–91. doi: 10.1016/j.radonc.2014.11.009
Keywords: knowledge based treatment planning, cardiac substructures, non-small cell lung cancer radiotherapy, machine learning, cardiac toxicities
Citation: Momin S, Wolf J, Roper J, Lei Y, Liu T, Bradley JD, Higgins K, Yang X and Zhang J (2022) Enhanced cardiac substructure sparing through knowledge-based treatment planning for non-small cell lung cancer radiotherapy. Front. Oncol. 12:1055428. doi: 10.3389/fonc.2022.1055428
Received: 27 September 2022; Accepted: 10 November 2022;
Published: 02 December 2022.
Edited by:
James Chow, University of Toronto, CanadaReviewed by:
Joshua Niedzielski, University of Texas MD Anderson Cancer Center, United StatesAhmed I. Ghanem, Henry Ford Health System, United States
Copyright © 2022 Momin, Wolf, Roper, Lei, Liu, Bradley, Higgins, Yang and Zhang. 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: Shadab Momin, c2hhZGFiLm1vbWluQGVtb3J5LmVkdQ==