Warning: mkdir(): Permission denied in /home/virtual/lib/view_data.php on line 87 Warning: chmod() expects exactly 2 parameters, 3 given in /home/virtual/lib/view_data.php on line 88 Warning: fopen(/home/virtual/e-kjs/journal/upload/ip_log/ip_log_2024-10.txt): failed to open stream: No such file or directory in /home/virtual/lib/view_data.php on line 95 Warning: fwrite() expects parameter 1 to be resource, boolean given in /home/virtual/lib/view_data.php on line 96 Overall Survival of Non-Small Cell Lung Cancer With Spinal Metastasis: A Systematic Review and Meta-Analysis

Overall Survival of Non-Small Cell Lung Cancer With Spinal Metastasis: A Systematic Review and Meta-Analysis

Article information

Neurospine. 2023;20(2):567-576
Publication date (electronic) : 2023 June 30
doi : https://doi.org/10.14245/ns.2245026.513
1Division of Neurosurgery, Department of Surgery, National Taiwan University Hospital, Taipei, Taiwan
2Spine Tumor Center, National Taiwan University Hospital, Taipei, Taiwan
3Department of Neurosurgery, Hallym University College of Medicine, Chuncheon, Korea
4Department of Neurosurgery, Yee Zen General Hospital, Taoyuan, Taiwan
5School of Health Sciences, Faculty of Biology Medicine and Health, University of Manchester, Manchester, UK
Corresponding Author Chung Liang Chai Department of Neurosurgery, Yee Zen General Hospital, 30, Yangshin North Road. Lane 321, Yangmei, Taoyuan, Taiwan 32645 Email: xyforce@gmail.com
Received 2022 December 1; Revised 2023 February 22; Accepted 2023 March 23.

Abstract

Objective

The long-term survival data of lung cancer patients with spinal metastases are crucial for informed treatment decision-making. However, most studies in this field involve small sample sizes. Moreover, survival benchmarking and an analysis of changes in survival over time are required, but data are unavailable. To meet this need, we performed a meta-analysis of survival data from small studies to obtain a survival function based on large-scale data.

Methods

We performed a single-arm systematic review of survival function following a published protocol. Data of patients who received surgical, nonsurgical, and mixed modes of treatment were meta-analyzed separately. Survival data were extracted from published figures with a digitizer program and then processed in R. Median survival time was used as an effect size for moderator analysis to explain the heterogeneity.

Results

Sixty-two studies with 5,242 participants were included for pooling. The survival functions showed a median survival of 6.72 months for surgery (95% confidence interval [CI], 61.9–7.01; 2,367 participants; 36 studies), 5.99 months for nonsurgery (95% CI, 5.33–6.47; 891 participants; 12 studies), and 5.96 months for mixed (95% CI, 5.67–6.43; 1,984 participants; 18 studies). Patients enrolled since 2010 showed the highest survival rates.

Conclusion

This study provides the first large-scale data for lung cancer with spinal metastasis that allows survival benchmarking. Data from patients enrolled since 2010 had the best survival and thus may more accurately reflect current survival. Researchers should focus on this subset in future benchmarking and remain optimistic in the management of these patients.

INTRODUCTION

Survival benchmarking is fundamental in cancer management. First, the survival of spine metastases of lung cancer is highly diverse as it is affected by age, functional status, molecular profiles such as epidermal growth factor receptor (EGFR) mutation, and many difficult-to-measure variables. When the survival outcome of this highly diverse population is derived from a small number of cases, the likelihood of misleading conclusions can be high. Zairi et al. [1] studied 53 surgical patients with lung cancer with spine metastasis and concluded that surgery is not recommended on account of short life expectancy and frailty. Although employing a large sample size would potentially avoid this, small sample size studies such as this are common in the current literature. One potential solution is to use data from the Global Cancer Observatory [2], which provides large-scale survival data for lung cancer that inform cancer control, research, and survival benchmarking. Although lung cancer data are available in this format, data regarding patients with spine metastases are not. This defect should be filled in the current era of big data-oriented cancer statistics.

Second, surgeons typically consider fixation or fusion instrumentation appropriate for patients who survive for longer than 1 year after surgical treatment. Wide excision is also considered for patients who survive longer. However, decompression alone or vertebroplasty can be used as alternative treatments for patients with limited survival. Therefore, when treating patients with recurrent lesions after local treatment, surgeons and radio-oncologists need expected survival data for treatment decision-making.

Furthermore, the influence of time is significant and includes progressing aging demographics and the advancement of multimodal oncological treatments. Therefore, survival is not constant but could subtly improve over time. Following the introduction of targeted therapies (such as TKIs) in early 2000, the survival rates of patients with lung cancer have gradually improved [2]. Accordingly, improved survival is also expected among patients with spinal metastasis; however, this assumption has not been substantiated. The Global Spinal Tumour Study Group attempted—but failed—to statistically demonstrate improved survival among patients with spinal metastasis of lung cancer [3]. Although this multinational collaborative study included 263 patients, the confidence intervals of the survival curves were wide, and statistically significant differences could not be established (Fig. 1).

Fig. 1.

Survival function of lung cancer patients with spinal metastasis from the Global Spinal Tumour Study Group [3]. This Kaplan-Meier curve has been re-estimated based on graphical coordinates (see Methods). The original curve[3] comprised 4 groups. Here, we have removed 4 patients in the 1991–2000 period to avoid distraction in a log-rank test comparing the remaining curves. The p-value is the result of a log-rank test. Reprinted from Wright E, et al. World Neurosurg 2018;114:e809-17 [3], with permission from Elsevier.

One solution to the above problem is a meta-analysis, which can help obtain a survival function based on large-scale data. Two types of patient populations need to be addressed separately: those who underwent surgery and those who received nonsurgical treatments. In a randomized study conducted in 2005, Patchell et al. [4] reported that surgery provided a clear advantage in metastatic epidural spinal cord compression. Since then, surgery has become the first option for patients who can physically tolerate it. Moreover, the surgical community seems to fully accept this conclusion, as no other randomized studies have been reported in this field since the study was published [4]. However, patients who are unlikely to tolerate surgery are considered to have a poorer baseline condition, regardless of cord compression or spinal instability. Therefore, the surgical treatment itself may act as a selection criterion that differentiates the baseline condition. To separately accommodate the 2 types of patients, the current study meta-analyzed each group separately. This review was conducted in accordance with the PRISMA (Preferred Reporting Items for Systematic Reviews and Meta-Analyses) 2020 guidelines [5].

MATERIALS AND METHODS

All data are provided within the article and are freely available from the Open Science Framework [6]. The protocol was prospectively registered in the International Prospective Register of Systematic Reviews (PROSPERO; CRD42021256577) on July 8, 2021.

This single-arm meta-analysis analyzed the survival function of non-small cell lung cancer (NSCLC) patients with spinal metastasis. Pathologies other than NSCLC and patients without spinal metastasis were excluded. The patient population was divided into 3 groups: surgery, nonsurgery, and mixed. Surgery was defined as patients undergoing operative treatment such as palliative decompression, debulking, total vertebrectomy, separation surgery, or any form of minimally invasive surgery. Nonsurgery was defined as treatments such as radiotherapy, vertebroplasty, chemotherapy, or palliative care. The mixed category contained combined survival data of surgical and nonsurgical patients. Only studies with extractable survival data were included. Studies with unextractable data, those that included nonNSCLC pathologies or nonspinal metastases, and those that did not report the total number of events were excluded. The context was restricted to clinical settings such as tertiary institutions, medical centers, and oncological centers. There were no restrictions on the type of study (randomized controlled trials or nonrandomized studies) as long as the data were extractable.

1. Finding and Assessing Individual Studies

We searched MEDLINE, Embase, and Google Scholar, as illustrated in Supplementary material 1. We also included the survival data of 172 patients from our institution (the Spine Oncology Registry of National Taiwan University Hospital). The search strategy for each database is detailed in Supplementary material 2. The search was performed on June 26, 2021. There were no restrictions on language. Study titles, abstracts, and full text were independently screened for inclusion by 2 authors (CLC and FYT), and discrepancies were resolved based on a consensus with a third author (JPJ). To facilitate transparency and reproducibility, we used a prepiloted form (Supplementary material 3) in the process of study inclusion and data extraction. No automation tool was used.

For qualitative data, the bias in the body of evidence was assessed according to the Grading of Recommendations, Assessment, Development, and Evaluation (GRADE) version dedic ated to prognostic factors [7]. As the current GRADE framework has no dedicated version to accommodate single-arm reviews, we narrowed the definitions according to the fundamental intention of the assessment (Supplementary material 4). The risk of bias (study limitations) in individual studies was assessed according to the respective GRADE items. All 3 outcomes were determined a priori, as presented in a summary of findings table. The GRADE level was independently judged by CLC and FYT, and discrepancies were resolved based on a consensus with JPJ.

2. Synthesizing the Body of Evidence

Data extraction was straightforward for studies that provided the survival data of each patient within the publication. For studies that provided a Kaplan-Meier curve, we used WebPlotDigitizer [8] to obtain the approximate coordinates in the image. The coordinates were processed using the IPDfromKM package [9] in R statistical software (R Development Core Team 2020, www.rproject.org) to obtain the survival data of each patient. Survival time was measured in months.

The survival data of individuals were pooled to estimate the survival function separately for the surgery, nonsurgery, and mixed groups. The data were also grouped chronologically to facilitate the assessment of changes over time. However, we did not predetermine the exact date of the partition. Moderator analysis was applied to explore the sources of statistical heterogeneity in study-level covariates, with median survival as the effect size. The 5 predetermined study-level covariates were as follows: mean age; starting date of study/subgroup inclusion; surgery or alternative treatment; and proportions of patients who received tyrosine kinase inhibitor (TKI), had an EGFR mutation, or had synchronous/metachronous malignancies. Conventional definitions were used for the confidence interval (95% confidence interval [CI]), heterogeneity (I2< 40%: unimportant; 30%–60%: moderate; 50%–90%: substantial, 75%–100%: considerable), and p-value (<0.05: significant). The online GRADEpro GDT software was utilized for GRADE assessment. The R packages survival, survminer, metafor, and metaCART were used for analysis.

RESULTS

A total of 2,244 studies were identified and screened (Supplementary material 1), of which 196 studies were assessed for eligibility and 62 were included in the meta-analysis. We documented the reasons for excluding 134 eligible studies (Supplementary material 5). Three studies were exclusively obtained from Google Scholar, one from the National Taiwan University Hospital Spine Oncology Registry, and the remaining 58 from either MEDLINE or Embase. The characteristics of included studies, including individual citations, are detailed in Supplementary material 6. The interreviewer agreement using Cohen kappa was 0.927 (96.9% agreement; Supplementary material 7). Among the included studies, one [10] was in Chinese, one [11] in German, and one [12] in Japanese; the rest were in English.

1. Survival Between Treatment Types

The summary of findings table and the pooled survival function of the 3 outcomes are shown in Table 1 and Fig. 2A, respectively. The pooled survival function estimated from the data of multiple studies (Fig. 2A) showed a median survival time of 6.72 months in the surgery group (95% CI, 61.9–7.01; 2,367 participants; 36 studies), 5.99 months in the nonsurgery group (95% CI, 5.33–6.47; 891 participants; 12 studies), and 5.96 months in the mixed group (95% CI, 5.67–6.43; 1,984 participants; 18 studies).

Summary table of findings

Fig. 2.

(A) Pooled survival curves of surgery, nonsurgery, and mixed study/subgroup. The solid line represents the survival function, and the semi-translucent area represents the 95% confidence interval. The p-value is the result of a log-rank test. (B) Survival over time. Survival data were pooled according to the time of initial recruitment in each study. The means by which studies were grouped is clarified in the text (see Methods and Supplementary material 8).

The large-scale survival data narrowed the CI significantly compared to those shown in Fig. 1. Sections of the curves where the CIs did not overlap suggested a meaningful difference in the survival times of patients between the surgery and nonsurgery groups. The nonsurgery group showed a lower survival probability in the period between 6 months and 30 months. Although a log-rank test indicated significant differences (p = 0.012) between the surgery, nonsurgery, and mixed, the curves began to merge at 30 months and overlapped significantly thereafter (Fig. 2A).

2. Survival Changes Over Time

Although the study protocol was designed to allow the analysis of potential survival differences over time, the data classification was not predetermined. Therefore, we partitioned the data a posteriori so as to evenly distribute the number of studies into 5 groups (Fig. 2B). The partition in 1999 was determined based on moderator analysis, as discussed below. Because the studies had been conducted across varying time spans, a clean partition was impossible. As a compromise, the data were sorted in ascending order according to the time of initial recruitment in each study. We admit that this is an a posteriori classification. However, the partitions can be further investigated and refined, as the data are freely available from the Open Science Framework [6].

There were several findings from this analysis. First, patients enrolled since 2010 had the best survival, with a nonoverlapping 95% CI. Second, the partition in 1999 was an output of the moderator analysis, which showed a large survival gap when data were separated. Third, when patients enrolled in 1989 or earlier were omitted, survival times improved in the period before 30 months. Beginning at the 30-month timepoint, all survival curves merged, and the survival probability ranged from 8.5% (95% CI, 6.74%–10.79%) to 11.7% (95% CI, 9.63%–14.15%). The 30-month survival has remained seemingly constant for the recent 2 decades. Finally, starting in 1999, more studies reported TKI usage or EGFR mutation status (see Supplementary material 8).

3. Moderator Analysis

Along with the 5 study-level covariates mentioned in the Methods section, we also accounted for the exclusiveness of NSCLC. At the stage of study inclusion, 34 studies reported the pathology as NSCLC. However, the pathology was unavailable (i.e., labeled as “lung cancer”) in 36 studies. After discussion, we decided to include both types of studies and perform moderator analysis for this covariate. We divided studies into subgroups to accommodate the study-level covariates. For example, one [13] of the included studies with an ID of “Yang SZ 2017” reported 2 sets of survival data: surgery and nonsurgery. These were separated into 2 subgroups in the moderator analysis: “Yang SZ 2017-su,” and “Yang SZ 2017-no,” respectively. The subgroups and extracted covariates of all studies are listed in Supplementary material 6. The metaCART [14] package was used to identify the point of partition that best explained the heterogeneity (see Supplementary material 9). The covariate of synchronous metastasis was not influential according to the analysis in metaCART. The influential covariates were further visualized with survival data in Fig. 3. Overall, 5 moderators (excluding the percentage of patients receiving TKIs, Fig. 3F) were significantly related to survival differences among groups.

Fig. 3.

The survival data were separately visualized according to the result of the moderator analysis. This demonstrates the effect of covariates on the survival outcome. The p-value is the result of a log-rank test. Starting year of cohort inclusion (A), exclusiveness of non-small cell lung cancer (NSCLC) (B), mean age (C), treatment types (D), % epidermal growth factor receptor (EGFR) mutation (E), and % targeted therapy (F).

The R2 value indicates how well our separation model explains heterogeneity. As in regression analysis, R2 values closer to 1 indicate a good model. Although all 6 separation models indicated a significant difference, the R2 values were low (see Supplementary material 9). Thus, the analysis accounted for heterogeneity to only a limited extent, and most heterogeneity remained unexplained (residual heterogeneity). This suggests that there was serious heterogeneity in the survival data included in this review (in terms of median overall survival). Although median survival represents only one timepoint (that of 50% survival) in the survival function, this was a significant finding that downgraded several GRADE items, as discussed in the next section.

4. GRADE Assessment

The GRADE scores of all outcomes indicated very low certainty of evidence (Table 1 and Supplementary material 4). Three main issues are summarized here. First, there was considerable heterogeneity in the data in terms of median survival. Second, studies that reported a single survival curve for surgical and nonsurgical treatments (i.e., those in the mixed category) seldom indicated their inclusion criteria or clinical severity. This can cause sampling discrepancies or bias between studies and make it difficult to define the target population, thus limiting the generalizability of outcomes. However, this issue was less severe in studies in the surgery and nonsurgery categories. Third, we considered studies from a small number of countries or regions, representing a minority of lung cancer cases (according to data from The Cancer Atlas [15]). This may have introduced publication bias (Fig. 4).

Fig. 4.

World map showing the distribution of studies included in this review. (A) Number of patients whose survival data were pooled. (B) Age-standardized rates (ASRs) of lung cancer incidence in males in 2018. The data have been obtained from The Cancer Atlas [15]. One study (ID: Choi 2015) has been excluded because it analyzed patients from 10 countries.

DISCUSSION

There have been significant advancements in the treatment of lung cancer with spinal metastasis in the past decade. Biomarker-directed therapies, like Osimertinib treatment, are among the newer TKIs therapies replacing previous targeted therapies [16]. In addition, studies of PD-1 and PD-L1 inhibitors have led to the increased use of immunotherapy with pembrolizumab since 2015 [16]. Stereotactic radiosurgery is regarded as a game-changer for spinal metastasis in other pathologies, with ongoing trials investigating different sessions [17]. Radiosurgery has also led to the development of minimally invasive procedures like separation surgery [18]. Despite these advancements, the durability of these treatments is not ideal [16]. Hence, it is crucial to determine whether they translate to survival benefits, especially in patients with spinal metastasis.

1. Survival Benchmarking Is Now Possible

To date, this is the first and the largest dataset representing a collection of all currently published survival. The dataset is open source and is available in the Open Science Framework [6]. This allows researchers to subset the dataset to complement any individual comparison. To demonstrate this, we compared the survival rates of the surgical cohort of the National Taiwan University Hospital Spine Oncology Registry against those of the dataset. For a meaningful comparison, we subset the dataset that underwent surgery beginning in 2004 (Fig. 5). A log-rank test indicated a p-value of 0.16, suggesting a similarity to those of the dataset. Furthermore, as described in the Introduction, fusion is typically considered appropriate for patients who survive for longer than 1 year. The 1-year survival of our institute was 41.3%, slightly higher than those of the (subset) dataset of 32.9%. This suggests fusion should be highly considered. Additionally, the high heterogeneity of our dataset reasonably reflected the high diversity of this group of patients. Hence, researchers need to utilize the covariates highlighted in the moderator analysis (Fig. 3) to subset the appropriate dataset to achieve a fair and meaningful comparison.

Fig. 5.

Example demonstration of benchmarking against institutional data (data from the Spine Oncology Registry of National Taiwan University Hospital). NTUH, National Taiwan University Hospital.

For a highly diverse population, benchmarking with a large-scale dataset can avoid misleading conclusions, as described in the Introduction. For instance, Zairi et al. [1] studied 53 surgical patients with lung cancer with spine metastasis and concluded that surgery is not recommended, on account of short life expectancy and frailty. Although our data suggest that the outcome of these patients is poor, it was not as poor as that described by Zairi et al. [1]. Frequently, the survival outcome deriving from a small number of cases of highly diverse population results in a Kaplan-Meier curve with large ladder steps and a wide 95% CI. Moreover, we found that the survival of surgical patients was slightly better than those of nonsurgical patients and was statistically significant (p = 0.0046) (Fig. 3D), although each was meta-analyzed separately.

2. Survival Changes Over Time

Upon analyzing survival changes over time (Fig. 2B), data from patients enrolled since 2010 had the best survival. This finding could suggest that the overall improvement of survival in lung cancer, in general, also applies to those with spinal metastasis. Although this is expected, it has not previously been presented in statistical terms. Since the data from 2010 may more accurately reflect current survival, researchers should focus on this subset in future benchmarking and remain optimistic in the management of these patients.

However, we also found that improvements were evident only prior to 30 months. The data visualization suggests a time limit of 30 months for survival analysis using this dataset, although this interpretation is not based on any statistical calculations. Our analysis also suggests that future research on comparative effectiveness should aim to reach the minimal period of 30-month follow-up before concluding superiority. On the contrary, studies that use a point estimate, commonly in conjunction with a median survival time, would be unlikely to identify the 30-month impasse. Rothrock et al. [19] demonstrated a gradual improvement with regression of median survival time in 309 lung cancer patients from 1998 to 2017. Although these findings are consistent with this study (Fig. 2B) in terms of median survival time, the improvement in outcome at 30 months was not demonstrated. This is because the median survival only serves as a point estimate of the full survival curve and does not encapsulate trends in the first and third quartiles or at any other points of interest. Overall, due to the a posteriori classification of data, these findings can only be considered exploratory.

3. Limitations

There were significant limitations to this study. First, for survival comparisons between surgery and other groups, any causative conclusions (such as those commonly derived in comparative effectiveness research) would be inappropriate. This is because the survival data were derived from single-arm rather than double-arm trials. Although the grouped data can be illustrated side-by-side, we cannot conclude that either has a relative survival advantage. Second, there was serious heterogeneity in the data. As there is no standard technique for the assessment of heterogeneity in non-parametric survival data, we used median survival as an effect size for heterogeneity assessment; moderator analysis only accounted for a small amount of heterogeneity. Third, the findings related to survival changes over time can only be regarded as exploratory. Studies varied significantly in their follow-up durations, resulting in significant overlaps at each time span when survival data were obtained (Supplementary material 8). The sorting and partitioning of survival data (Fig. 2B) was an a posteriori attempt. Overall, this resulted in very low certainty of evidence, as assessed under the GRADE framework.

CONCLUSION

This is the first large-scale data regarding lung cancer with spinal metastasis that allows survival benchmarking. Data from patients enrolled since 2010 had the best survival and may more accurately reflect current survival. Researchers should focus on this subset in future benchmarking and remain optimistic in the management of these patients.

Supplementary Material

Supplementary materials 1-9 can be found via https://doi.org/10.14245/ns.2245026.513.

Supplementary material 1.

Flowchart in the PRISMA format

ns-2245026-513-Supplementary-1.pdf
Supplementary material 2.

Search strategy for studies included in this review

ns-2245026-513-Supplementary-2.pdf
Supplementary material 3.

Data extraction form

ns-2245026-513-Supplementary-3.pdf
Supplementary material 4.

GRADE definitions and judgment process

ns-2245026-513-Supplementary-4.pdf
Supplementary material 5.

Reasons for excluding 134 eligible studies

ns-2245026-513-Supplementary-5.pdf
Supplementary material 6.

Characteristics of included studies

ns-2245026-513-Supplementary-6.pdf
Supplementary material 7.

Study selection and inter-assessor agreement

ns-2245026-513-Supplementary-7.pdf
Supplementary material 8.

Time span of all included studies/subgroups

ns-2245026-513-Supplementary-8.pdf
Supplementary material 9.

Moderator analysis

ns-2245026-513-Supplementary-9.pdf

Notes

Conflict of Interest

The authors have nothing to disclose.

Funding/Support

This study was supported by Good Liver Foundation and Taiwan Brain Disease Foundation of Taiwan (MS476).

Author Contribution

Conceptualization: FYT, CLC; Data curation: JPJ, APH, CLC; Formal analysis: FYT, JPJ, CLC; Funding acquisition: FYT; Methodology: CLC; Project administration: FYT, JPJ, CLC; Visualization: CLC; Writing - original draft: FYT, CLC; Writing - review & editing: FYT, CLC.

References

1. Zairi F, Karnoub MA, Vieillard MH, et al. Evaluation of the relevance of surgery in a retrospective case series of patients who underwent the surgical treatment of a symptomatic spine metastasis from lung cancer. Eur Spine J 2016;25:4052–9.
2. Ferlay J, Ervik M, Lam F, et al. Global cancer observatory: cancer today [Internet]. Lyon (France): International Agency for Research on Cancer; 2020. [cited 2022 Apr 10]. Available from: https://gco.iarc.fr/today.
3. Wright E, Ricciardi F, Arts M, et al. Metastatic spine tumor epidemiology: comparison of trends in surgery across two decades and three continents. World Neurosurg 2018;114:e809–17.
4. Patchell RA, Tibbs PA, Regine WF, et al. Direct decompressive surgical resection in the treatment of spinal cord compression caused by metastatic cancer: a randomised trial. Lancet 2005;366:643–8.
5. Page MJ, McKenzie JE, Bossuyt PM, et al. The PRISMA 2020 statement: an updated guideline for reporting systematic reviews. BMJ 2021;372:n71.
6. Chai CL, Tsuang FY, Jeon JP. Overall survival of non-small cell lung cancer with spinal metastasis: a systematic review and meta-analysis [Internet]. OSF; 2022. [cited 2022 Apr 10]. Available from: https://osf.io/qyvne/.
7. Foroutan F, Guyatt G, Zuk V, et al. GRADE guidelines 28: use of GRADE for the assessment of evidence about prognostic factors: rating certainty in identification of groups of patients with different absolute risks. J Clin Epidemiol 2020;121:62–70.
8. Rohatgi A. WebPlotDigitizer [Internet]. 2021. [cited 2022 Apr 10]. Available from: https://automeris.io/WebPlotDigitizer.
9. Liu N, Zhou Y, Lee JJ. IPDfromKM: reconstruct individual patient data from published Kaplan-Meier survival curves. BMC Med Res Methodol 2021;21:111.
10. Li Y, Jiang L, Liu X, et al. Efficacy and survival rate analysis of lung cancer with spinal metastases. Beijing Da Xue Xue Bao 2014;46:138–43.
11. Ulmar B, Huch K, Kocak T, et al. The prognostic influence of primary tumour and region of the affected spinal segment in 217 surgical patients with spinal metastases of different entities. Z Orthop Ihre Grenzgeb 2007;145:31–8.
12. Fukuda Y, Yoneda S, Noguchi Y, et al. An analysis of 45 clinical cases of spinal cord compression by vertebral metastases from lung cancer. Jpn J Lung Cancer 1989;29:119–26.
13. Yang SZ, Tang Y, Zhang Y, et al. Prognostic factors and comparison of conservative treatment, percutaneous vertebroplasty, and open surgery in the treatment of spinal metastases from lung cancer. World Neurosurg 2017;108:163–75.
14. Li X, Dusseldorp E, Meulman JJ. Meta-CART: a tool to identify interactions between moderators in meta-analysis. Br J Math Stat Psychol 2017;70:118–36.
15. The Cancer Atlas. Lung cancer [Internet]. Atlanta (GA): American Cancer Society, Inc.; 2019. [cited 2022 Apr 10]. Available from: https://canceratlas.cancer.org/the-burden/lung-cancer/.
16. Arbour KC, Riely GJ. Systemic therapy for locally advanced and metastatic non-small cell lung cancer: a review. JAMA 2019;322:764–74.
17. ClinicalTrials.gov. Lung cancer spine metastasis [Internet]. Bethesda (MD): ClinicalTrials.gov; [cited 2022 Apr 10]. Available from: https://clinicaltrials.gov/ct2/results?cond=lung+cancer+spine+metastasis&term=&cntry=&state=&city=&dist=.
18. Barzilai O, Fisher CG, Bilsky MH. State of the art treatment of spinal metastatic disease. Neurosurgery 2018;82:757–69.
19. Rothrock RJ, Barzilai O, Reiner AS, et al. Survival trends after surgery for spinal metastatic tumors: 20-year cancer center experience. Neurosurgery 2021;88:402–12.

Article information Continued

Fig. 1.

Survival function of lung cancer patients with spinal metastasis from the Global Spinal Tumour Study Group [3]. This Kaplan-Meier curve has been re-estimated based on graphical coordinates (see Methods). The original curve[3] comprised 4 groups. Here, we have removed 4 patients in the 1991–2000 period to avoid distraction in a log-rank test comparing the remaining curves. The p-value is the result of a log-rank test. Reprinted from Wright E, et al. World Neurosurg 2018;114:e809-17 [3], with permission from Elsevier.

Fig. 2.

(A) Pooled survival curves of surgery, nonsurgery, and mixed study/subgroup. The solid line represents the survival function, and the semi-translucent area represents the 95% confidence interval. The p-value is the result of a log-rank test. (B) Survival over time. Survival data were pooled according to the time of initial recruitment in each study. The means by which studies were grouped is clarified in the text (see Methods and Supplementary material 8).

Fig. 3.

The survival data were separately visualized according to the result of the moderator analysis. This demonstrates the effect of covariates on the survival outcome. The p-value is the result of a log-rank test. Starting year of cohort inclusion (A), exclusiveness of non-small cell lung cancer (NSCLC) (B), mean age (C), treatment types (D), % epidermal growth factor receptor (EGFR) mutation (E), and % targeted therapy (F).

Fig. 4.

World map showing the distribution of studies included in this review. (A) Number of patients whose survival data were pooled. (B) Age-standardized rates (ASRs) of lung cancer incidence in males in 2018. The data have been obtained from The Cancer Atlas [15]. One study (ID: Choi 2015) has been excluded because it analyzed patients from 10 countries.

Fig. 5.

Example demonstration of benchmarking against institutional data (data from the Spine Oncology Registry of National Taiwan University Hospital). NTUH, National Taiwan University Hospital.

Table 1.

Summary table of findings

Outcomes Anticipated absolute number of patients survival* (95% CI)
Median overall survival (95% CI) No. of participants (studies) Certainty of the evidence (GRADE)
6 Months 1 Year 2 Years
Surgery 5.3 per 3.2 per 1.4 per 6.72 Months 2,367 Very low
10 (5.1–5.5) 10 (3.0–3.4) 10 (1.2–1.5) (6.19–7.01) (36 studies) ⦿◯◯◯
Nonsurgery 4.9 per 2.5 per 1.2 per 5.99 Months 891 Very low
10 (4.6–5.3) 10 (2.3–2.8) 10 (1.0–1.4) (5.33–6.47) (12 studies) ⦿◯◯◯
Mixed cohorts§ 5.0 per 3.0 per 1.4 per 5.96 Months 1,984 Very low
10 (4.7–5.2) 10 (2.8–3.2) 10 (1.2–1.5) (5.67–6.43) (18 studies) ⦿◯◯◯

CI, confidence interval.

*

The absolute data is estimated from the percentage of survival (and its 95% confidence interval).

Decompression, instrumentation, en-bloc, anterior/posterior approaches, or combinations.

Radiotherapy, stereotactic body ratiotherapy, chemotherapy, vertebroplasty, etc.

§

Varying proportion of surgical and nonsurgical data integrated into a single survival function.

Patient or population: Populations diagnosed with non-small cell lung cancer that presents with spinal metastasis.

Setting: Tertiary institutions, medical centers, oncological centers.

Exposure: There is no designated intervention or exposure in this study because the aim of this review is obtain the survival curve.

Comparison: There is no designated comparator.

GRADE Working Group grades of evidence risk associated with the prognostic factor.

High certainty: We are very confident that the variation in risk associated with the prognostic factor (probability of future events in those with/without the prognostic factor) lies close to that of the estimate.

Moderate certainty: We are moderately confident that the variation in risk associated with the prognostic factor (probability of future events in those with/without the prognostic factor) is likely to be close to the estimate, but there is a possibility that it is substantially different.

Low certainty: Our certainty in the estimate is limited: the variation in risk associated with the prognostic factor (probability of future events in those with/without the prognostic factor) may be substantially different from the estimate.

Certainty: We have very little certainty in the estimate: the variation in risk associated with the prognostic factor (probability of future events in those with/without the prognostic factor) is likely to be substantially different from the estimate.