Skip to main content

ORIGINAL RESEARCH article

Front. Oncol., 27 March 2023
Sec. Radiation Oncology
This article is part of the Research Topic Prospective Utilization and Clinical Applications of Artificial Intelligence and Data-driven Automation for Radiotherapy View all 8 articles

Development and validation of a unifying pre-treatment decision tool for intracranial and extracranial metastasis-directed radiotherapy

Roman KowalchukRoman Kowalchuk1Trey C. MullikinTrey C. Mullikin2William BreenWilliam Breen2Hunter C. GitsHunter C. Gits2Marcus FlorezMarcus Florez3Brian DeBrian De3William S. HarmsenWilliam S. Harmsen4Peter Sean RosePeter Sean Rose5Brittany L. SiontisBrittany L. Siontis6Brian A. CostelloBrian A. Costello6Jonathan M. MorrisJonathan M. Morris7John J. LucidoJohn J. Lucido8Kenneth R. OlivierKenneth R. Olivier1Brad StishBrad Stish1Nadia N. LaackNadia N. Laack1Sean ParkSean Park1Dawn OwenDawn Owen1Amol J. GhiaAmol J. Ghia3Paul D. BrownPaul D. Brown1Kenneth Wing Merrell*Kenneth Wing Merrell1*
  • 1Department of Radiation Oncology, Mayo Clinic, Rochester, MN, United States
  • 2Department of Radiation Oncology, Duke University, Durham, NC, United States
  • 3Department of Radiation Oncology, Houston, MD Anderson Cancer Center, Houston, TX, United States
  • 4Mayo Clinic, Department of Statistics, Rochester, MN, United States
  • 5Mayo Clinic, Department of Orthopedic Surgery, Rochester, MN, United States
  • 6Mayo Clinic, Department of Radiology, Rochester, MN, United States
  • 7Mayo Clinic, Department of Medical Oncology, Rochester, MN, United States
  • 8Mayo Clinic, Department of Medical Physics, Rochester, MN, United States

Background: Though metastasis-directed therapy (MDT) has the potential to improve overall survival (OS), appropriate patient selection remains challenging. We aimed to develop a model predictive of OS to refine patient selection for clinical trials and MDT.

Patients and methods: We assembled a multi-institutional cohort of patients treated with MDT (stereotactic body radiation therapy, radiosurgery, and whole brain radiation therapy). Candidate variables for recursive partitioning analysis were selected per prior studies: ECOG performance status, time from primary diagnosis, number of additional non-target organ systems involved (NOS), and intracranial metastases.

Results: A database of 1,362 patients was assembled with 424 intracranial, 352 lung, and 607 spinal treatments (n=1,383). Treatments were split into training (TC) (70%, n=968) and internal validation (IVC) (30%, n=415) cohorts. The TC had median ECOG of 0 (interquartile range [IQR]: 0-1), NOS of 1 (IQR: 0-1), and OS of 18 months (IQR: 7-35). The resulting model components and weights were: ECOG = 0, 1, and > 1 (0, 1, and 2); 0, 1, and > 1 NOS (0, 1, and 2); and intracranial target (2), with lower scores indicating more favorable OS. The model demonstrated high concordance in the TC (0.72) and IVC (0.72). The score also demonstrated high concordance for each target site (spine, brain, and lung).

Conclusion: This pre-treatment decision tool represents a unifying model for both intracranial and extracranial disease and identifies patients with the longest survival after MDT who may benefit most from aggressive local therapy. Carefully selected patients may benefit from MDT even in the presence of intracranial disease, and this model may help guide patient selection for MDT.

Introduction

Recent studies have explored oligometastatic disease, the idea that select patients with limited metastatic disease burden may benefit from metastasis-directed therapy (MDT) (1). MDT typically includes interventional radiology techniques, surgical resection, and radiotherapy, often using stereotactic planning techniques involving a sharp dose fall-off to a focal target (2). Clinical trials are ongoing to optimize patient selection for MDT to maximize cost-effectiveness, quality of life, and survival while minimizing side effects of treatment (36). The SABR-COMET trial randomized 99 patients with one to five metastatic sites to standard of care with or without stereotactic body radiation therapy (SBRT) to each of the metastatic lesions, demonstrating an improvement in overall survival (OS) with SBRT (2, 79). The ORIOLE and STOMP trials specifically assessed patients with primary prostate cancer and up to 3 sites of metastatic disease, and MDT improved 6 month progression-free and androgen deprivation therapy-free survival (10, 11). These encouraging results have led to the design of numerous other clinical trials studying MDT, including NRG-BR002 and SABR-COMET-10 (12, 13).

Despite differences in inclusion of primary histologies and type of MDT delivered, most clinical trials share the same fundamental patient inclusion criteria. First, the majority of clinical trials evaluating outcomes in oligometastatic cancer include patients with 5 or fewer metastatic lesions, with the assumption that the number of lesions, not the number of organs involved, is the ideal stratification. Second, clinical trials commonly exclude patients with intracranial disease based on the assumption that patients with brain metastases have a poor prognosis (14). However, recent analyses have questioned these assumptions and stratification factors, with emerging evidence that the number of organ systems involved by metastatic disease may better predict OS than the number of individual metastatic lesions (1518). Furthermore, prognostication by the number of metastatic lesions is challenging in the setting of new imaging modalities that identify previously occult sites of disease, whose influence on a patient’s prognosis remains to be clearly defined (19). Some patients treated for brain metastasis may also have durable OS similar to metastatic disease to other organs, and applicability of intracranial-specific models to other sites of MDT (e.g. spine SBRT) is yet untested (20, 21). Substantial work involving the generation of high-fidelity prognostic models has been pioneered by Sperduto et al., and our group aims to extend this work regarding intracranial disease by focusing on extracranial disease (20).

Given the clinical uncertainty of the ideal patient selection for aggressive MDT, we analyzed a large, multi-institutional database of patients treated with MDT involving the brain, lungs, or spine to develop and validate a unifying prognostic model applicable for both intracranial and extracranial disease to guide treatment selection and future clinical trial design.

Methods

Patient cohort

The dataset was comprised of 4 separate databases from 2 large academic centers. The first institution provided 3 cohorts of patients who received MDT (1): brain, (2) lung, and (3) spine metastases. The second institution contributed a single cohort of patients who received MDT for spine metastases. More detailed information concerning each underlying dataset is described in the corresponding publications, as well as in Supplemental Table 1 (17, 22, 23). MDT for brain metastases included whole brain radiation therapy (WBRT) and/or stereotactic radiosurgery (SRS) (24). This decision was made to allow for comparison to the prognostic models generated by Sperduto et al. (20). Furthermore, given the lack of consensus concerning the threshold (either volumetric or number of metastatic lesions) for recommendation of WBRT and SRS, both were considered within the analysis. MDT for spine and lung metastasis included only SBRT. Treatments with conventional palliative dose-fractionation schemes for spine and lung targets were excluded. The decision to treat a patient as having “oligometastatic disease” was defined per the clinical discretion of the treating multi-disciplinary team, as there is no validated consensus definition. Decisions regarding treating all sites of disease with MDT or alternative approaches (e.g. treatment of only a dominant lesion) were also left to the discretion of the treating physicians, secondary to a lack of consensus clinical standard of practice. Institutional review board (IRB) approval was obtained at both institutions. Benign primary tumors were excluded. If a patient had multiple treatments targeting the same organ (e.g. two spine SBRT treatments), only the first such treatment was included.

Treatment

Treatment planning and delivery were conducted per institutional protocol. Rigid immobilization was utilized for stereotactic radiotherapy treatments, with narrow PTV margins of 0-5 mm. Motion management was considered for all cases of lung SBRT, including use of breath hold and phase gating, per the discretion of the treating physician. Stereotactic treatments were generally delivered using 1-5 daily treatments. Treatments of intracranial targets involved standard use of a 3-point mask for immobilization. Treatments with Gamma Knife were also conducted in accordance with institutional protocol. For WBRT, hippocampal avoidance and memantine were used whenever possible, with 10 daily fractions offered for most patients.

Model generation

Candidate variable selection was determined per prior work involved the development of a prognostic model after SBRT for spinal metastases, with variable selection targeted for the primary endpoint of OS (17). This included thorough analysis of other potential prognostic factors, including but not limited to patient characteristics (e.g. age), local tumor characteristics (e.g. gross tumor volume), systemic therapy, and histology. After assessment of a wide range of candidate variables, the final model included: Eastern Cooperative Oncology Group (ECOG) performance status, number of non-target organ systems (NOS), and time from primary diagnosis (TPD). NOS was defined as the number of organ systems involved with disease outside of the organ system of the treatment target (i.e. bone for spine SBRT, lung for lung SBRT, and brain for intracranial treatments). Further, intracranial target was added as a fourth potential variable. The overall cohort was split into training (70%, n=968) and internal validation (30%, n=415) cohorts via a random number generator, and comparable patient characteristics were confirmed in both sets. Feature importance testing was conducted and a correlation heatmap was generated. The overarching training set was used for model generation and RPA analysis while the internal validation set was only used afterwards to validate the prognostic value of the resulting risk score.

Recursive partitioning analysis (RPA) was conducted with open-source packages in Python (version 3.8.0), with similar methodology to that reported in the literature (2530). The overarching training set was used for all modeling efforts. RPA was conducted through decision-tree analysis using the overarching training set (70% of the overall cohort), split into training (60%), validation (20%), and test sets (20%). Potential RPA models were generated with the training and validation sets, and model performance was assessed via the independent test set (4). A minimum of 15 patients in each group was required for node splitting, but this criterion was analyzed as a range from 15-20 minimum patients. The log-rank test was used as the criterion for node splitting, as the primary endpoint of the analysis was OS. RPA trees were only allowed to contain a maximum of 4 groups to minimize overfitting. The highest-fidelity models were then considered, and the prognostic factors used in these models were considered for inclusion in the final risk score.

Model validation

Additional statistical analyses were conducted with SAS version 9.4 (SAS Institute, Cary, NC). Cox proportional hazards analyses were conducted, and the resulting concordance value was reported as a measure of the fidelity of the resulting model and underlying variables. Such analyses were conducted for both the overarching training and internal validation sets to assess the prognostic value of the final risk score and individual component variables. Cox proportional hazards analyses were also performed for each distinct treatment target (brain, lung, and spine). Hazard ratios (HR) were reported, along with 95% confidence intervals (CI). For all analyses, p<0.05 was the threshold for statistical significance. Finally, the Kaplan-Meier method was used to plot OS as stratified by the resulting risk score.

Results

Patient cohort

A total of 1,260 individual patients were included with 424 brain, 352 lung, and 607 spinal targets (n=1,383). Further details concerning the patient demographics and primary histologies of each treatment target site are shown in Supplemental Table 1. For spine SBRT, 326 (54%) treatments involved single-fraction radiotherapy, and the median dose delivered for all fractionation schemes was 24 Gy (IQR: 20-30). For lung SBRT, a median of 5 fractions (IQR: 3-5) were used to deliver a median dose of 50 Gy (IQR: 50-54). Concerning intracranial treatments, 145 (34%) patients received initial WBRT, and 349 (82%) had initial SRS with or without WBRT. 78 patients (18%) had more than 10 brain metastases at the time of brain metastasis diagnosis, and of the remaining patients, a median 2 (IQR: 1-3) brain metastases were identified. Many patients with intracranial disease also had extracranial disease (75%), most commonly bone (43%) and lung (38%) lesions (Supplemental Table 2). The training set had median ECOG of 0 (IQR: 0-1), NOS of 1 (IQR: 0-1), and OS of 18 months (IQR: 7-35). The internal validation set demonstrated comparable characteristics to the training set (Table 1).

TABLE 1
www.frontiersin.org

Table 1 Patient characteristics utilized for recursive partitioning analysis are demonstrated, as divided into training and validation sets.

Risk score generation

Feature importance testing failed to exclude any of the four candidate variables (ECOG, NOS, TPD, and intracranial) for consideration in the resulting risk score. A correlation heatmap also did not reveal substantial correlation between the variables, so all 4 were included as inputs into the RPA. Though ECOG, NOS, and intracranial were frequently included in the highest-fidelity models, TPD was not. It was, therefore, excluded from further consideration in the final model. Two separate thresholds for ECOG and NOS were identified via RPA: ECOG = 1 and > 1, and NOS = 1 and > 1. Cox proportional hazards analysis was performed and each of these thresholds demonstrated a statistically significant association with OS (Table 2). Corresponding weights were then approximated for development of the corresponding risk score: ECOG = 0, 1, and > 1 (0, 1, and 2, respectively); NOS = 0, 1, and > 1 (0, 1, and 2); and intracranial target (2) (Table 3). While ECOG > 1 had a higher hazard ratio (HR=4.22, 95% CI: 3.18-5.61, p<0.0001) than NOS > 1 (HR=2.12, 95% CI: 1.73-2.60, p<0.0001) and intracranial (HR=2.60, 95% CI: 2.16-3.13, p<0.0001), there were only 39 patients with ECOG > 1 in the internal validation set. For this reason, a higher relative weight for ECOG > 1 than NOS > 1 or intracranial was not assigned. The resulting risk score was plotted using the Kaplan-Meier method, revealing an association between the risk score and overall survival in the training and internal validation sets (Figure 1, p<0.0001). 2-year OS was 76%, 64%, 41%, 28%, 21%, 11%, and 12% for scores of 0, 1, 2, 3, 4, 5, and 6, respectively.

TABLE 2
www.frontiersin.org

Table 2 The predictive power of the resulting model is demonstrated for the training and validation sets, with each individual model component also demonstrating statistical significance.

TABLE 3
www.frontiersin.org

Table 3 The model components and their corresponding weights are tabulated.

FIGURE 1
www.frontiersin.org

Figure 1 Overall survival is demonstrated via the prognostic score is demonstrated via the training set (A) and validation set (B).

Risk score validation

Cox proportional hazards analysis revealed a concordance value of 0.73 (95% CI: 0.71-0.75) for the risk score in predicting OS in the training set. Similarly, the risk score had a concordance value of 0.72 (95% CI: 0.68-0.74) for predicting OS in the internal validation set. Each individual variable within the risk score also demonstrated statistical significance in both the training and internal validation sets (p<0.0001) (Table 2). Each stratification of these variables also revealed a statistically significant association with OS (p ≤ 0.01). Next, the applicability of the internal validation set for each individual score of the resulting risk score was assessed (Figure 2). The observed OS for each risk group was within the 95% confidence interval of the predicted OS, except for the small group of score 6. The score of 6 only occurred 11 times in the training set and was, therefore, only a very small fraction of the final prognostic score.

FIGURE 2
www.frontiersin.org

Figure 2 A demonstration of the adequacy of the validation set for each score of the final prognostic model. For all scores but one, the number of events in the validation set fit within the 95% confidence interval of the predicted event frequencies. The only exception was the score of 6, which only occurred 11 times in the training set and was, therefore, only a very small fraction of the final prognostic score.

Treatment target analyses

The resulting risk score was next individually applied to each of the 3 treatment targets involved in the overall cohort. Concordance values were 0.72 (95% CI: 0.68-0.74), 0.71 (95% CI: 0.68-0.75), and 0.65 (95% CI: 0.62-0.68) for lung, spine, and brain, respectively. Considering patients with intracranial targets, statistical significance was maintained when considering patients with only a single brain metastasis (n=150, p=0.0003) or for patients with multiple brain metastases (n=274, p<0.0001). Both ECOG and NOS demonstrated statistical significance in each subset (p<0.0001) (Supplemental Table 3). Comparable hazard ratios were identified in the 3 subsets for each stratification of ECOG and NOS. Each stratification of these variables also revealed a statistically significant association with OS, except for NOS = 1 for the lung cohort (HR=1.31, 95% CI: 0.89-1.94, p=0.17). Similarly, the model demonstrated prognostic value for each of the most common primary histologies (breast, lung, and GI), and each individual variable within the risk score demonstrated a statistically significant association with OS within each primary histology (Supplemental Table 4).

Discussion

Using a large, multi-institutional cohort, we developed and validated a pre-treatment prognostic risk score for OS to aid in selection of patients most likely to benefit from metastasis-directed radiotherapy. Recent advances in systemic therapy have substantially improved OS for patients with metastatic disease, allowing additional opportunities to incorporate MDT for patients (31, 32). Therefore, our model provides a tool to identify patients who may derive the most benefit from MDT, and it also identifies those who may be spared MDT and the associated risks of toxicity (4, 8). Unique to our study is the inclusion of both intracranial and extracranial targets, and our resulting concordance index was higher in both the training and validation sets than prior models assessing only extracranial disease (15, 18). The inclusion of patients with brain metastasis adds critical data to the topic of patient selection for MDT reflective of true clinical practice, as approximately a quarter of patients who die of cancer develop brain metastasis (33, 34). Therefore, our model may be useful to guide future MDT clinical trial inclusion criteria design, as well as to provide a unifying model for patients with intracranial and extracranial disease.

Quantification of a patient’s overall metastatic disease burden is a useful tool to guide clinical decision making for MDT and characterize a patient’s prognosis. One clear example is in the setting of prostate cancer, where clinical trials suggest that the role of radiotherapy differs in patients with low volume vs. high volume metastatic disease (35). The optimal manner to stratify other metastatic patients, however, remains unclear. Past approaches include reporting the number of total metastatic lesions, total volume of metastatic lesions, or number of organ systems involved by metastatic disease (36). Other data suggest that the number of metastatic lesions alone is an incomplete prognostic variable, and the number of non-target organ systems involved by disease or the presence of only a solitary metastasis may be superior prognostic variables (15, 16). For example, the classification from the CHAARTED trial included two criteria for high volume metastatic burden: 4 or more bone metastases with one or more outside the vertebral bodies or pelvis, or visceral metastases (or both) (37). That is, five metastases within vertebral bodies would only classify as low metastatic burden, just as our model would have accounted for this scenario with an NOS of 0. Alternatively, the presence of a single vertebral body lesion and a single visceral metastasis would have qualified as high metastatic burden, and our model would have incorporated an NOS of 1. Stratification by only the number of metastatic lesions, however, would have allowed the latter patient to enroll in many oligometastatic clinical trials due to the presence of only 2 sites of metastatic disease. This difference may allow for more accurate prognostication using our model and NOS, in line with the CHAARTED definition.

An additional concern with stratification by only the number of metastatic lesions is that as imaging sensitivity and sophistication increases, the number of previously occult metastatic lesions identified are likely to increase, making comparison and applicability of prior studies using this classification challenging (19). Therefore, prior analyses would have underestimated the true number of metastatic sites, and it remains unknown whether the presence of additional, small, clinically occult lesions within the same organ system substantially impacts patient prognosis.

Perhaps the most unique aspect of our study was the inclusion of intracranial targets. A substantial fraction (34%) of these patients were treated with WBRT, which involves treating the entire brain to cover subclinical disease. Furthermore, 75% of patients treated for intracranial disease also had extracranial disease, and the model still effectively stratified such patients by OS. This result suggests that patients with intracranial metastatic disease and otherwise favorable characteristics (e.g. ECOG of 0-1 and NOS of 0-1) could still be candidates for MDT, despite a larger number of individual metastatic lesions than included in most clinical trials. Therefore, we propose that carefully selected patients with intracranial metastatic disease be further studied to consider potential inclusion in future clinical trials (20).

ECOG performance status was identified as an important prognostic variable in our analysis, consistent with prior studies (17). In the setting of brain metastases, rigorous graded prognostic assessment has consistently identified performance status as a crucial prognostic variable across multiple primary tumor histologies (20). Multiple prognostic systems have been constructed for spine SBRT incorporating performance status as a key variable (15, 16, 38, 39). Systematic reviews have noted encouraging inter-rater reliability for performance status scoring, with no clear superiority of one system over the other (40, 41). Overall, we strongly encourage utilization of performance status to aid in prognostication of OS and ultimately patient selection for MDT.

A limitation of our analysis is that the model was developed using retrospective data. However, the large size of our patient cohort (1,260 patients) and multi-institutional involvement, are substantial strengths of our analysis. Patients were also treated with modern imaging and radiotherapy treatment techniques, allowing us to adequately capture current clinical practice. Next, while we did not directly test the number of metastatic lesions as a potential prognostic variable, the number of metastatic lesions would have demonstrated strong correlation with NOS. We acknowledge that clinical judgement will be required to incorporate this model in a given clinical context. Finally, this model is most generalizable to patients with GI, breast, and NSCLC primary tumors, which constituted the histologies of patients with intracranial metastases included in this study. Even so, patients with spine SBRT and lung SBRT had a range of primary histologies, and primary tumor histology was not found to be a key prognostic variable in that setting (42). Therefore, the authors feel that it would be appropriate to generalize this model outside of only primary GI, breast, and NSCLC.

Conclusions

Our pre-treatment decision tool identifies patients most likely to benefit from metastasis-directed radiotherapy secondary to favorable OS and applies for both intracranial and extracranial disease. We encourage future clinical trial design to consider ECOG and NOS as prognostic variables for patients with metastatic disease. Carefully selected patients may benefit from MDT even in the presence of intracranial disease.

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.

Ethics statement

The studies involving human participants were reviewed and approved by Mayo Clinic Institutional Review Board. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author contributions

All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.

Conflict of interest

The wife of RK is a senior technical product manager for GE Healthcare.

The remaining 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.2023.1095170/full#supplementary-material

References

1. Weichselbaum RR, Hellman S. Oligometastases revisited. Nat Rev Clin Oncol (2011) 8(6):378–82. doi: 10.1038/nrclinonc.2011.44

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Gomez DR, Tang C, Zhang J, Blumenschein GR, Hernandez M, Lee JJ, et al. Local consolidative therapy vs. maintenance therapy or observation for patients with oligometastatic non–Small-Cell lung cancer: Long-term results of a multi-institutional, phase II, randomized study. J Clin Oncol (2019) 37(18):1558–65. doi: 10.1200/JCO.19.00201

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Kowalchuk RO, Mullikin TC, Kim DK, Jonathan MM, Ebner DK, Harmsen WS, et al. Cost-effectiveness of treatment strategies for spinal metastases. Pract Radiat Oncol (2022) 12(3):236–44. doi: 10.1016/j.prro.2021.12.012

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Kowalchuk RO, Johnson-Tesch BA, Marion JT, Mullikin T, Harmsen W, Rose P, et al. Development and assessment of a predictive score for vertebral compression fracture after stereotactic body radiation therapy for spinal metastases. JAMA Oncol (2022) 8(3):412–419. doi: 10.1001/jamaoncol.2021.7008

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Kowalchuk RO, Niranjan A, Hess J, Antonios JP, Zhang MY, Braunstein S, et al. Stereotactic radiosurgery and local control of brain metastases from triple-negative breast cancer. J Neurosurg (2022), 1–7. doi: 10.3171/2022.10.JNS221900

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Lehrer EJ, Kowalchuk RO, Gurewitz J, Bernstein K, Kondziolka D, Niranjan A, et al. Concurrent administration of immune checkpoint inhibitors and single fraction stereotactic radiosurgery in patients with non-small cell lung cancer, melanoma, and renal cell carcinoma brain metastases is not associated with an increased risk of radiation necrosis over non-concurrent treatment: An international multicenter study of 657 patients. Int J Radiat Oncol Biol Phys (2023) S0360-3016(23):00057-3. doi: 10.1016/j.ijrobp.2023.01.017

CrossRef Full Text | Google Scholar

7. Palma DA, Olson R, Harrow S, Gaede S, Louie AV, Haasbeek C, et al. Stereotactic ablative radiotherapy for the comprehensive treatment of oligometastatic cancers: Long-term results of the SABR-COMET randomized trial. Int J Radiat Oncol Biol Phys (2020) 108(3):S88–S9. doi: 10.1016/j.ijrobp.2020.07.2251

CrossRef Full Text | Google Scholar

8. Palma DA, Olson R, Harrow S, Gaede S, Louie AV, Haasbeek C, et al. Stereotactic ablative radiotherapy versus standard of care palliative treatment in patients with oligometastatic cancers (SABR-COMET): a randomised, phase 2, open-label trial. Lancet (2019) 393(10185):2051–8. doi: 10.1016/S0140-6736(18)32487-5

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Iyengar P, Wardak Z, Gerber DE, Tumati V, Ahn C, Hughes RS, et al. Consolidative radiotherapy for limited metastatic non–Small-Cell lung cancer. JAMA Oncol (2018) 4(1):e173501. doi: 10.1001/jamaoncol.2017.3501

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Ost P, Reynders D, Decaestecker K, Fonteyne V, Lumen N, De Bruycker A, et al. Surveillance or metastasis-directed therapy for oligometastatic prostate cancer recurrence: A prospective, randomized, multicenter phase II trial. J Clin Oncol (2018) 36(5):446–53. doi: 10.1200/JCO.2017.75.4853

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Phillips R, Shi WY, Deek M, Radwan N, Lim SJ, Antonarakis ES, et al. Outcomes of observation vs stereotactic ablative radiation for oligometastatic prostate cancer: The ORIOLE phase 2 randomized clinical trial. JAMA Oncol (2020) 6(5):650–9. doi: 10.1001/jamaoncol.2020.0147

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Chmura SJ, Winter KA, Al-Hallaq HA, Borges VF, Jaskowiak NT, Matuszak M, et al. NRG-BR002: A phase IIR/III trial of standard of care therapy with or without stereotactic body radiotherapy (SBRT) and/or surgical ablation for newly oligometastatic breast cancer (NCT02364557) American Society of clinical oncology. (2019). doi: 10.1200/JCO.2019.37.15_suppl.TPS1117

CrossRef Full Text | Google Scholar

13. Palma DA, Olson R, Harrow S, Correa RJM, Schneiders F, Haasbeek CJA, et al. Stereotactic ablative radiotherapy for the comprehensive treatment of 4–10 oligometastatic tumors (SABR-COMET-10): study protocol for a randomized phase III trial. BMC Cancer (2019) 19(1). doi: 10.1186/s12885-019-5977-6

CrossRef Full Text | Google Scholar

14. Andrews DW, Scott CB, Sperduto PW, Flanders AE, Gaspar LE, Schell MC, et al. Whole brain radiation therapy with or without stereotactic radiosurgery boost for patients with one to three brain metastases: Phase III results of the RTOG 9508 randomised trial. Lancet (2004) 363(9422):1665–72. doi: 10.1016/S0140-6736(04)16250-8

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Jensen G, Tang C, Hess KR, Bishop AJ, Pan HY, Li J, et al. Internal validation of the prognostic index for spine metastasis (PRISM) for stratifying survival in patients treated with spinal stereotactic radiosurgery. J Radiosurg SBRT (2017) 5(1):25. doi: 10.1016/j.ijrobp.2017.02.085

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Tang C, Hess K, Bishop AJ, Pan HY, Christensen EN, Yang JN, et al. Creation of a prognostic index for spine metastasis to stratify survival in patients treated with spinal stereotactic radiosurgery: Secondary analysis of mature prospective trials. Int J Radiat OncologyBiologyPhysics (2015) 93(1):118–25. doi: 10.1016/j.ijrobp.2015.04.050

CrossRef Full Text | Google Scholar

17. Kowalchuk RO, Mullikin TC, Harmsen WS, Rose PS, Siontis BL. Development and internal validation of an RPA-based pre-treatment decision tool for spinal SBRT. ESTRO (2021). doi: 10.1016/S0167-8140(21)06892-4

CrossRef Full Text | Google Scholar

18. Chen H, Poon I, Atenafu EG, Badellino S, Biswas T, Dagan R, et al. Development of a prognostic model for overall survival in patients with extracranial oligometastatic disease treated with stereotactic body radiation therapy. Int J Radiat Oncol Biol Phys (2021) 114(5):892–901. doi: 10.1016/j.ijrobp.2021.12.001

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Hofman MS, Hicks RJ, Maurer T, Eiber M. Prostate-specific membrane antigen PET: clinical utility in prostate cancer, normal patterns, pearls, and pitfalls. Radiographics (2018) 38(1):200–17. doi: 10.1148/rg.2018170108

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Sperduto PW, Mesko S, Li J, Cagney D, Aizer A, Lin NU, et al. Survival in patients with brain metastases: Summary report on the updated diagnosis-specific graded prognostic assessment and definition of the eligibility quotient. J Clin Oncol (2020) 38(32):3773–84. doi: 10.1200/JCO.20.01255

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Gaspar LE, Scott C, Murray K, Curran W. Validation of the RTOG recursive partitioning analysis (RPA) classification for brain metastases. Int J Radiat Oncol Biol Physics (2000) 47(4):1001–6. doi: 10.1016/S0360-3016(00)00547-2

CrossRef Full Text | Google Scholar

22. Sperduto PW, Mesko S, Li J, Cagney D, Aizer A, Lin NU, et al. Beyond an updated graded prognostic assessment (breast GPA): a prognostic index and trends in treatment and survival in breast cancer brain metastases from 1985 to today. Int J Radiat Oncol Biol Physics (2020) 107(2):334–43. doi: 10.1016/j.ijrobp.2020.01.051

CrossRef Full Text | Google Scholar

23. Sperduto PW, Fang P, Li J, Breen W, Brown PD, Cagney D, et al. Estimating survival in patients with gastrointestinal cancers and brain metastases: An update of the graded prognostic assessment for gastrointestinal cancers (GI-GPA). Clin Trans Radiat Oncol (2019) 18:39–45. doi: 10.1016/j.ctro.2019.06.007

CrossRef Full Text | Google Scholar

24. Oncology N. Stereotactic radiosurgery compared with hippocampal-avoidant whole brain radiotherapy (HA-WBRT) plus memantine for 5-15 brain metastases.

Google Scholar

25. Davidson-Pilon C. Lifelines: Survival analysis in Python. J Open Source Software (2019) 4(40):1317. doi: 10.21105/joss.01317

CrossRef Full Text | Google Scholar

26. Kowalchuk RO, Waters MR, Dutta SW, Mack ML, Richardson KM, Spencer K, et al. Recursive partitioning analysis for local control achieved with stereotactic body radiation therapy for the liver, spine, or lymph nodes. Adv Radiat Oncol (2021) 6(3):100612. doi: 10.1016/j.adro.2020.10.025

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Kowalchuk RO, Shepard MJ, Sheehan K, Sheehan D, Faramand A, Niranjan A, et al. Treatment of WHO grade 2 meningiomas with stereotactic radiosurgery. Identification an Optimal Group SRS Using RPA (2021) 110(3):804–814. doi: 10.1016/j.ijrobp.2021.01.048

CrossRef Full Text | Google Scholar

28. Pedregosa F, Varoquaux G, Gramfort A, Michel V, Thirion B, Grisel O, et al. Scikit-learn: Machine learning in Python. J Mach Learn Res (2011) 12:2825–30.

Google Scholar

29. Tosi S. Matplotlib for Python developers: Packt publishing Ltd. (2009).

Google Scholar

30. Kowalchuk R, Mullikin T, Harmsen W, Rose P, Siontis B, Kim D, et al. Development and internal validation of an RPA-based model predictive of pain flare incidence after spine SBRT. Int J Radiat Oncol Biol Phys (2021) 111(3):S60–S1. doi: 10.1016/j.ijrobp.2021.07.154

CrossRef Full Text | Google Scholar

31. Gandhi L, Rodríguez-Abreu D, Gadgeel S, Esteban E, Felip E, De Angelis F, et al. Pembrolizumab plus chemotherapy in metastatic non–Small-Cell lung cancer. New Engl J Med (2018) 378(22):2078–92. doi: 10.1056/NEJMoa1801005

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Mok TS, Wu Y-L, Kudaba I, Kowalski DM, Cho BC, Turna HZ, et al. Pembrolizumab versus chemotherapy for previously untreated, PD-L1-expressing, locally advanced or metastatic non-small-cell lung cancer (KEYNOTE-042): As randomised, open-label, controlled, phase 3 trial. Lancet (2019) 393(10183):1819–30. doi: 10.1016/S0140-6736(18)32409-7

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Amsbaugh MJ, Kim CS. Brain metastasis (2017).

Google Scholar

34. Nayak L, Lee EQ, Wen PY. Epidemiology of brain metastases. Curr Oncol Rep (2012) 14(1):48–54. doi: 10.1007/s11912-011-0203-y

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Parker CC, James ND, Brawley CD, Clarke NW, Hoyle AP, Ali A, et al. Radiotherapy to the primary tumour for newly diagnosed, metastatic prostate cancer (STAMPEDE): A randomised controlled phase 3 trial. Lancet (2018) 392(10162):2353–66. doi: 10.1016/S0140-6736(18)32486-3

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Oh Y, Taylor S, Bekele BN, Debnam JM, Allen PK, Suki D, et al. Number of metastatic sites is a strong predictor of survival in patients with nonsmall cell lung cancer with or without brain metastases. Cancer (2009) 115(13):2930–8. doi: 10.1002/cncr.24333

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Sweeney CJ, Chen Y-H, Carducci M, Liu G, Jarrard DF, Eisenberger M, et al. Chemohormonal therapy in metastatic hormone-sensitive prostate cancer. New Engl J Med (2015) 373(8):737–46. doi: 10.1056/NEJMoa1503747

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Chao ST, Koyfman SA, Woody N, Angelov L, Soeder SL, Reddy CA, et al. Recursive partitioning analysis index is predictive for overall survival in patients undergoing spine stereotactic body radiation therapy for spinal metastases. Int J Radiat Oncol Biol Physics (2012) 82(5):1738–43. doi: 10.1016/j.ijrobp.2011.02.019

CrossRef Full Text | Google Scholar

39. Balagamwala EH, Miller JA, Reddy CA, Angelov L, Suh JH, Tariq MB, et al. Recursive partitioning analysis is predictive of overall survival for patients undergoing spine stereotactic radiosurgery. J neuro-oncology (2018) 137(2):289–93. doi: 10.1007/s11060-017-2716-1

CrossRef Full Text | Google Scholar

40. Chow R, Bruera E, Temel JS, Krishnan M, Im J, Lock M. Inter-rater reliability in performance status assessment among healthcare professionals: An updated systematic review and meta-analysis. Supportive Care Cancer (2020) 28(5):2071–8. doi: 10.1007/s00520-019-05261-7

CrossRef Full Text | Google Scholar

41. Chow R, Chiu N, Bruera E, Krishnan M, Chiu L, Lam H, et al. Inter-rater reliability in performance status assessment among health care professionals: A systematic review. Ann Palliat Med (2016) 5(2):83–92. doi: 10.21037/apm.2016.03.02

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Kowalchuk RO, Mullikin TC, Florez M, De BS, Spears GM, Rose PS, et al. Development and validation of a recursive partitioning analysis–based pretreatment decision-making tool identifying ideal candidates for spine stereotactic body radiation therapy. Cancer (2022) 129(6):956–965. doi: 10.1002/cncr.34626

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: metastatic disease, metastasis-directed radiotherapy, oligometastasis, modeling, outcomes

Citation: Kowalchuk R, Mullikin TC, Breen W, Gits HC, Florez M, De B, Harmsen WS, Rose PS, Siontis BL, Costello BA, Morris JM, Lucido JJ, Olivier KR, Stish B, Laack NN, Park S, Owen D, Ghia AJ, Brown PD and Merrell KW (2023) Development and validation of a unifying pre-treatment decision tool for intracranial and extracranial metastasis-directed radiotherapy. Front. Oncol. 13:1095170. doi: 10.3389/fonc.2023.1095170

Received: 10 November 2022; Accepted: 15 March 2023;
Published: 27 March 2023.

Edited by:

Xinglei Shen, University of Kansas Medical Center, United States

Reviewed by:

Yoseful Haque, University of Kansas, United States
Jay Shiao, University of Kansas, United States

Copyright © 2023 Kowalchuk, Mullikin, Breen, Gits, Florez, De, Harmsen, Rose, Siontis, Costello, Morris, Lucido, Olivier, Stish, Laack, Park, Owen, Ghia, Brown and Merrell. 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: Kenneth Wing Merrell, Merrell.Kenneth@mayo.edu

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.