Next Article in Journal
Modern Immunotherapy in the Treatment of Triple-Negative Breast Cancer
Next Article in Special Issue
Stereotactic Radiosurgery for Lung Cancer with a Risk-Adapted Strategy Using the Volumetric Modulated Arc Therapy Technique: A Single Arm Phase II Study
Previous Article in Journal
Circulating SOD2 Is a Candidate Response Biomarker for Neoadjuvant Therapy in Breast Cancer
Previous Article in Special Issue
Relationship between Dose Prescription Methods and Local Control Rate in Stereotactic Body Radiotherapy for Early Stage Non-Small-Cell Lung Cancer: Systematic Review and Meta-Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Establishment of a Prediction Model for Overall Survival after Stereotactic Body Radiation Therapy for Primary Non-Small Cell Lung Cancer Using Radiomics Analysis

Department of Radiology, University of Tokyo Hospital, 7-3-1, Hongo, Bunkyo-ku, Tokyo 113-8655, Japan
*
Author to whom correspondence should be addressed.
Cancers 2022, 14(16), 3859; https://doi.org/10.3390/cancers14163859
Submission received: 30 June 2022 / Revised: 4 August 2022 / Accepted: 9 August 2022 / Published: 10 August 2022

Abstract

:

Simple Summary

Lung cancer remains the leading cause of cancer-related mortality worldwide. Although early-stage non-small cell lung cancer (NSCLC) is likely to be controlled with stereotactic body radiation therapy (SBRT), approximately 18% of patients lead to recurrence. The aim of this study was to identify prognostic factors and establish a predictive model for survival outcomes of patients with non-metastatic NSCLC treated with SBRT. Several radiomic features were selected as predictive factors and two prediction models were established from the pre-treatment computed tomography images of 250 patients in the training cohort. One radiomic factor remained a significant prognostic factor of overall survival (OS) (p = 0.044), and one predicting model could estimate OS time (mean: 37.8 months) similar to the real OS time (33.7 months). In this study, we identified one radiomic factor and one prediction model that can be widely used.

Abstract

Stereotactic body radiation therapy (SBRT) for early-stage non-small cell lung cancer (NSCLC) leads to recurrence in approximately 18% of patients. We aimed to extract the radiomic features, with which we predicted clinical outcomes and to establish predictive models. Patients with primary non-metastatic NSCLC who were treated with SBRT between 2002 and 2022 were retrospectively reviewed. The 358 primary tumors were randomly divided into a training cohort of 250 tumors and a validation cohort of 108 tumors. Clinical features and 744 radiomic features derived from primary tumor delineation on pre-treatment computed tomography were examined as prognostic factors of survival outcomes by univariate and multivariate analyses in the training cohort. Predictive models of survival outcomes were established from the results of the multivariate analysis in the training cohort. The selected radiomic features and prediction models were tested in a validation cohort. We found that one radiomic feature showed a significant difference in overall survival (OS) in the validation cohort (p = 0.044) and one predicting model could estimate OS time (mean: 37.8 months) similar to the real OS time (33.7 months). In this study, we identified one radiomic factor and one prediction model that can be widely used.

1. Introduction

Lung cancer is one of the most common causes of cancer-related mortality worldwide [1,2,3]. Non-small cell lung cancer (NSCLC) is the major subtype of lung cancer, accounting for approximately 85% of all lung cancer cases [4]. Stereotactic body radiation therapy (SBRT) is currently one of the standard curative treatment options for patients with early-stage NSCLC who are medically inoperable or refuse surgery [5,6,7,8,9]. While a very high local tumor control rate of 92–98% can be expected [10,11,12,13,14], SBRT for early-stage NSCLC leads to local, regional, or distant metastatic recurrence in 18–20% of cases [13,14,15,16].
Appropriate prediction of recurrence risk before treatment is essential for personalized treatment. If we can select patients with a high possibility of recurrence in advance, we may increase the treatment intensity by, for example, changing the total radiation dose and fractions. It has been reported that mediastinal pleural contact, maximum tumor diameter, and the maximum standard uptake value (SUV-max) in positron emission tomography (PET) are associated with distant metastasis as clinical features in patients with stage I-II NSCLC treated with SBRT [17]. Other groups have reported that the biologically effective dose (BED) is also related to the prognosis of patients treated with SBRT for early-stage NSCLC [18]. Clinical T stage, tumor size, and total radiation dose also had a statistically significant impact on survival outcomes of patients with early-stage NSCLC treated with SBRT [19,20]. Nowadays, many quantitative features can be extracted from radiographic images using radiomics. Radiomic features potentially include information about tumor heterogeneity, tumor microenvironment, and underlying gene-expression patterns, which may be associated with prognosis [21,22]. Radiomics approaches are widely used to detect prognostic factors in several types of cancers, such as hepatic cancer, brain tumor, pancreatic cancer, and head and neck cancer [21,23,24,25,26]. Several articles have reported prognostic factors in patients with early-stage NSCLC who underwent SBRT using radiomics [27,28,29,30,31,32,33]. They concluded that some radiomic features were prognostic factors and showed predictive models of prognosis; however, there is currently no widely used prediction model. It is possible that the number of patients, which was approximately 100, was not large enough to establish a versatile prediction model.

2. Materials and Methods

2.1. Patients

This study was reviewed and approved by the institutional review board and ethics committee (examination number 3372). We retrospectively reviewed an institutional database to identify patients with primary non-metastatic NSCLC who were treated with SBRT with curative intent between April 2002 and March 2022 with or without salvage treatment. Patients who had plain chest CT scans within two weeks prior to SBRT were included. Patients who had never been examined to confirm the clinical outcomes after SBRT were excluded. We detected 358 primary tumors in 338 patients. A total of 250 primary tumors were randomly assigned to the training cohort, and 108 tumors were assigned to the validation cohort.

2.2. Radiotherapy

All patients with no pathological diagnosis by biopsy required a continuous increase in the primary tumor size or abnormal accumulation of fluorodeoxyglucose in the tumor before SBRT. Tumor markers (neuron-specific enolase, NSE, and pro-gastrin-releasing peptide, ProGRP) needed to be checked to exclude small cell carcinoma. The patient’s forced expiratory volume in one second needed to be more than 750 cc to receive SBRT. In addition, indications for SBRT must be decided by an institutional cancer board consisting of pulmonologists, thoracic surgeons, radiologists, and radiation oncologists.
Six or 10 MV photon linear accelerators were used for SBRT. The radiation technique used was three-dimensional conformal radiotherapy (3D-CRT) or volumetric modulated arc therapy (VMAT). The fixed multi-port technique consisting of 7–12 beams was used as 3D-CRT. A total dose of 42–64 Gy in 4–10 fractions was administered. When the α/β ratio is assumed to be 10, BED was 75–166 Gy. The total dose and fractions were determined by radiation oncologists by considering the location of the tumor and the expected dose of organs at risk. Gross tumor volume (GTV) was contoured on the planning CT images displayed at a lung window level. The clinical target volume (CTV) was defined to be the same as the GTV. An internal margin was added to the CTV to cover the respiratory motion of the tumor and generate the internal target volume (ITV) using four-dimensional CT images. A tracking system was not used. Then, a 5 mm margin in every direction was added to the ITV to generate the planning target volume (PTV). Radiation dose was prescribed to 95% of the PTV, and radiation was administered every weekday.

2.3. Clinical Endpoints

As background features, age, Karnofsky performance status (KPS), sex, total radiation dose, radiotherapy (RT) technique (fixed multi-ports, VMAT, or both), primary tumor location, maximum tumor diameter, and SUV-max in PET were investigated.
Overall survival (OS), local relapse-free survival (LRFS), and progression-free survival (PFS) were used as the endpoints. These were counted from the date of SBRT initiation. The patients were usually followed-up and underwent a CT scan every 2–3 months for the first year and repeated 3–6 months thereafter. Local recurrence, which meant a recurrence in the radiation field, was usually identified on PET when SUV-max of the tumor was higher than 2.5 or by biopsy in case of obvious enlargement of the local tumor on the follow-up CT scan. However, some relapses were judged only by continuous enlargement on follow-up CT.

2.4. CT Image and Tumor Contouring

All patients underwent plain chest CT for treatment planning. We used Aquilion LB (Canon Medical Systems Corporation, Otawara, Japan) as the CT scanner. The CT tube voltage was 120 kV, the tube current was 350 mA, and the slice thickness was 2 mm. Every primary tumor was manually delineated as the GTV for this study on the expiratory phase of the planning CT displayed with a lung window level using Monaco (Elekta AB, Stockholm, Sweden) by a single expert radiation oncologist. Tumor outline was contoured, and pleura, vessels, heart, bones, and chest wall were excluded from the GTV.

2.5. Radiomic Feature Extraction

We used PyRadiomics v3.0.1 (Boston, MA, USA) to extract radiomic features from contoured regions of interest (ROIs) of GTVs [34]. A fixed bin count of 64 was used for discretization of the image gray level. All CT voxels were resampled to 1 × 1 × 1 mm3 using a B-Spline interpolation function. A total of 107 quantitative features were automatically extracted, including 19 first-order statistics features (intensity histogram, IH), 26 shape-based histogram features, and texture features (gray-level co-occurrence matrix, GLCM, 24 features; gray-level run-length matrix, GLRLM, 16 features; gray-level size-zone matrix, GLSZM, 16 features; neighboring gray-tone difference matrix, NGTDM, 5 features; and gray-level dependence matrix, GLDM, 14 features). All quantitative features are listed in Table S1. From these original radiomic features, 744 quantitative features were obtained using wavelet transformation. Coiflets 1, which is one of the methods of wavelet transformation, was used in x, y and z axes in this study. As a result, 8 wavelet-transformed images (HHH, HHL, HLH, HLL, LHH, LHL, LLH, and LLL) were obtained, where H represents high frequency and L represents low frequency. Mathematical definitions of these radiomic features have previously been described and are available at https://pyradiomics.readthedocs.io/en/latest/features.html (accessed on 11 June 2022).

2.6. Statistical Analysis

The Kaplan–Meier method was used for survival time analysis. Clinical features of continuous variables consisted of maximum tumor diameter, age, KPS, SUV-max in PET, and total radiation dose, and radiomic features were divided into two groups according to the median values in the training cohort. Univariate analysis was performed using the log-rank test to assess the difference in clinical outcomes between the high- and low-value groups. Multivariate analysis was conducted using the Cox proportional hazards model among the factors with a significant difference in univariate analysis. The stepwise method with the Bayesian information criterion was used for variable selection in the training cohort (direction = “backward/forward”). We created a prediction score from the results of the Cox proportional hazards model and multiple linear regression analysis of the training cohort. The selected radiomic factors and prediction scores were validated in the validation cohort. All statistical analyses were performed using R v4.1.3. The results of the univariate analysis of radiomic features were considered statistically significant at p < 0.01, whereas the other results were considered statistically significant at p < 0.05.

3. Results

3.1. Patient Characteristics

Patients were randomly divided into training (n = 250) and validation (n = 108) groups. As shown in Table 1, there was no bias in the proportions or means of both the patient and tumor background factors in either group. In the training group, 68% were male, mean age was 77 (range, 42–93) years, 84% had peripheral lesions, mean maximum tumor diameter was 20.5 mm (4.0–90.0), mean total dose was 52.5 Gy (42–64), 70% received VMAT, and 12.8% received salvage treatment after recurrence was recognized. No patient underwent adjuvant treatment after SBRT. The median follow-up time for patients who survived at the last follow-up of the training cohort and the validation cohort was 33.0 (0.9–167.9) and 33.9 (0.5–149.4) months, respectively.
The median crude survival of the training cohort and the validation cohort was 30.4 (interquartile range (IQR), 15.9–50.9) and 28.2 (14.2–44.1) months, respectively. Similarly, the median crude LRFS was 24.6 (10.4–45.8) vs. 23.2 (11.4–38.1) months and the median crude PFS was 23.3 (9.7–45.0) vs. 21.3 (11.0–38.4).

3.2. Univariate Analysis

As shown in Table 2, as background factors affecting OS, KPS, maximum tumor diameter, sex, peripheral lesions, and salvage treatment showed a significant difference of p < 0.05 in the univariate analysis using the log-rank test. Regarding LRFS, KPS, maximum tumor diameter, sex, peripheral lesions, consolidation/tumor (C/T) ratio, salvage treatment, and SUV-max showed a significant difference of p < 0.05 in the univariate analysis. For PFS, KPS, peripheral lesions, and salvage treatment showed a significant difference of p < 0.05 in the univariate analysis.
Among the 744 radiomic features, 70 had a significant difference of p < 0.01 in OS between the high- and low-value groups in the training cohort in the log-rank test. For LRFS, 223 factors, and for PFS, 210 factors showed a significant difference of p < 0.01.

3.3. Multivariate Analysis

As shown in Table 3, multiple factors with significant differences in univariate analysis were subjected to the Cox proportional hazards regression model. As a result, four factors, (1) “90 Percentile_HHH,” (2) “LargeAreaEmphasis_LHH,” (3) “Mean_HHH,” and (4) “Median_HLL” remained as factors related to OS. Similarly, regarding LRFS, two factors, (1) “InverseVariance_HLL” and (2) “SmallDependenceHighGrayLevelEmphasis_HHH” remained. Regarding PFS, three factors, (1) “SmallDependenceHighGrayLevelEmphasis_HHH,” (2) “TotalEnergy_HHL,” and (3) “JointEntropy_HLL” remained.

3.4. Prediction Score

From the results of the Cox proportional hazards regression model, we created a prediction score, which was named the Cox-score, calculated by the following formula using the coef value as a coefficient.
Cox-scoreOS = −0.6545 × [90 Percentile_HHH] + 0.5871 × [LargeAreaEmphasis_LHH] + 0.6795 × [Mean_HHH] + 0.6852 × [Median_HLL]
Cox-scoreLRFS = 1.2655 × [InverseVariance_HLL] − 1.5453 × [SmallDependenceHighGrayLevelEmphasis_HHH]
Cox-scorePFS = −0.9901 × [InverseVariance_HLL] + 0.9006 × [InverseVariance_HLL] − 0.6994 × [InverseVariance_HLL]
For the [x] value, “+1” is substituted if the value of x is higher than the median value, and “−1” is substituted if the value is lower than the median value.
We estimated survival time using multiple linear regression analysis in the training cohort. The estimated survival time was calculated using the following formula.
Estimated OS time (months) = 49.128 + 9.528 × [90 Percentile_HHH] − 11.909 × [LargeAreaEmphasis_LHH] − 6.781 × [Mean_HHH] − 11.658 × [Median_HLL]
Estimated LRFS time (months) = 61.431 − 22.649 × [InverseVariance_HLL] + 31.459 × [SmallDependenceHighGrayLevelEmphasis_HHH]
Estimated PFS time (months) = 38.400 + 8.804 × [JointEntropy_HLL] + 21.647 × [SmallDependenceHighGrayLevelEmphasis_HHH] − 14.024 × [TotalEnergy_HHL]
For the [x] value, “+1” is substituted if the value of x is higher than the median value, and “0” is substituted if the value is lower than the median value.

3.5. Validation

As a result of the log-rank test of 108 cases of the validation cohort for each of the radiomic factors remaining in the multivariate analysis, only one factor of “LargeAreaEmphasis_LHH” showed a significant difference in OS (p = 0.044, 5-year OS of 70.7% vs. 50.3%) (Figure 1). “LargeAreaEmphasis_LHH” was a parameter that was categorized in GLSZM features. Regarding LRFS and PFS, none of the factors that showed a significant difference in the multivariate analysis of the training cohort showed a significant difference in the univariate analysis of the validation cohort.
Regarding the Cox-score, the value for each case was calculated using the distribution of the radiomic features of 108 cases in the validation group. As a result, 38 cases with a score of less than −0.1 were divided into a low value group, 36 cases with a score of +1.0 or more were divided into a high value group, and the other 34 cases were divided into a median value group for OS. In the log-rank test, the difference in OS between the three groups was not significant (p = 0.63, 5-year OS of 62.6% in the high score group, 59.6% in the median, and 65.4% in the low score group). The difference in LRFS between the two groups, if +1.0 was the cutoff value, was not significant (p = 0.086, 5-year OS of 60.8% in the high, and 82.8% in the low). The difference in PFS between the two groups, if −1.0 was the cutoff value, was not significant (p = 0.070, 5-year OS of 52.0% in the high, and 94.7% in the low).
We compared the estimated survival time calculated using the formula from multiple linear regression analysis with the real survival time in the validation cohort. As a result of the estimated OS time, the difference in the average values was 4.1 months (33.7 months in the real survival time vs. 37.8 months in the estimated survival time) and the t value and p value were 1.619 and 0.11 by the paired t-test, respectively (Figure 2). Regarding LRFS and PFS, the p value was less than 0.0001 using the paired t-test.

4. Discussion

In this study, we extracted radiomic features from the GTV on SBRT planning CT of patients with non-metastatic primary NSCLC and selected features correlated with survival outcomes. We established a prediction model using these selected features. The selected features and prediction models were validated in an independent external cohort. “LargeAreaEmphasis_LHH” categorized in GLSZM features remained as a factor significantly correlated with OS in the validation cohort. The estimated OS time calculated using the formula from the multiple linear regression analysis was similar to the actual OS time in the validation cohort. There was no correlation between the radiomic features and LRFS or PFS.
Prediction of recurrence risk before treatment for early-stage NSCLC is so important. The treatment intensity such as the total radiation dose, fractions, and radiation field can be increased for high-risk patients. Another measure to increase the treatment intensity is combined therapy. Ernani et al. [35] showed the promising OS data of SBRT and adjuvant chemotherapy for patients with tumors less than 4 cm and node-negative NSCLC. Use of immune checkpoint inhibitors as an adjuvant therapy after curative chemoradiotherapy for advanced stage NSCLC have rapidly spread worldwide [36,37]. Several studies have also reported the effectiveness of immune checkpoint inhibitors together with SBRT for high-risk patients with early-stage NSCLC [38]. The selected radiomic features and OS prediction model derive in the present study may contribute to selecting patients who benefit from these treatment options.
Several studies have investigated the prognostic factors of patients who received SBRT for early-stage NSCLC using radiomics analysis. Oikonomou et al. [32] examined CT and PET-derived radiomic features in 150 patients with non-metastatic NSCLC treated with SBRT. They examined several principal components, consisting of six to eight radiomic features. Four principal components were found to be significantly associated with the survival outcomes. In the combined analysis with SUV-max, radiomic features remained the only predictors of OS, disease specific survival, and regional control. Huynh et al. [30] analyzed 113 patients with stage I–II NSCLC who were treated with SBRT. Among the 1605 radiomic features extracted from the tumors in the free breathing CT images, one radiomic feature was a significant prognostic factor for distant metastasis and four radiomic features were prognostic for OS. Starkov et al. [33] retrospectively analyzed 116 patients treated with SBRT for biopsy-confirmed primary NSCLC at a single institution. They showed that quantitative image features from NSCLC nodules on pre-treatment CT images were correlated with OS after SBRT. Similar to the present study, no significant relationship was observed between radiomic features and survival outcomes except for OS. No external validation cohort was established in any of these studies. In contrast, Dissaux et al. [28] built a prediction model of local control derived from PET among a training set of 64 early-stage NSCLC patients treated with SBRT. The model maintained a significant correlation with local control in the validation set of 23 patients. The fact that patients in the training set were collected from two institutions and those in the validation set were collected from another single institution may have affected the results of the study. Yu et al. [39] retrospectively extracted radiomic features from contrast-enhanced CT images of a training cohort (n = 147) and an independent validation cohort (n = 295) of patients with stage I NSCLC. Although the cohort size was large, there was an important difference between the two cohorts. The prediction model consisting of two radiomic features was significantly associated with both OS and distant metastasis in the validation cohort. Patients in the training cohort were treated with surgery, and those in the validation cohort were treated with SBRT. In addition, the patients in the validation cohort were significantly older than those in the training cohort, which seemed to be due to the difference in the indication for treatment. These differences may have affected the accuracy of the validation. In contrast, the training and validation cohorts were randomly divided in the present study. Thus, the clinical background features were not significantly different between the two cohorts. Franceschini et al. [29] also analyzed radiomic features in the treatment planning CT images of 70 patients treated with SBRT for early-stage NSCLC for training and 32 patients for validation. The two cohorts were randomly divided as in the present study. A set of 45 textural features was extracted from the tumor volumes on the treatment planning CT images. They built a predictive model that was prognostic of PFS and disease-specific survival. Although these studies included a validation process, the sample size may not have been large enough. In our study, we examined much more NSCLC nodules, which were randomly divided into training and validation cohorts to minimize the risk of bias. In addition to background factors, crude OS, LRFS, and PFS were similar between the two cohorts. Several radiomic features remained as factors correlated with PFS and LRFS in the multivariate analysis of the training cohort; however, none of them were found to be prognostic factors in the external validation process. Unless validation is conducted, these radiomic features could be identified as prognostic factors.
The present study has some limitations. First, although PET or biopsy was recommended to identify local recurrence, there were some cases in which local recurrence was judged using only CT. Thus, there may have been an overdiagnosis of local recurrence. Second, all patients were recruited from a single institution and all tumors were contoured by a single radiation oncologist. External validation at various institutions is expected to be conducted to develop the established prediction model into a widely usable model. Third, the fact that all GTVs were manually contoured may have reduced reproducibility of this study. Homogeneity of this study was kept since all contouring tasks were conducted by a single expert radiation oncologist.

5. Conclusions

We found that the radiomic features categorized as GLSZM features derived from the GTV on the pre-treatment CT image were an independent prognostic factor of OS in non-metastatic NSCLC patients treated with curative SBRT. Our findings suggest that the estimated OS time derived from the multiple linear regression analysis was similar to the real OS time. These findings may contribute to selecting patients who benefit from increasing the treatment intensity by, for example, increasing radiation dose and using combined therapy.
The selected radiomic factor and OS prediction model derived in this study can be widely used in the future. External validation at various institutions is expected.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers14163859/s1, Table S1: List of all 107 radiomic features.

Author Contributions

Conceptualization, S.S. and H.Y.; Methodology, S.S. and H.Y.; Software, Y.N.; Validation, S.S., H.Y. and Y.N.; Formal Analysis, H.Y.; Investigation, S.S., H.Y., R.T., Y.M., K.M., H.U. and A.K.; Resources, S.S.; Data Curation, S.S. and T.O.; Writing—Original Draft Preparation, S.S. and H.Y.; Writing—Review and Editing, S.S., H.Y., R.T. and A.K.; Visualization, S.S.; Supervision, H.Y.; Project Administration, H.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Institutional Ethics Committee of the University of Tokyo Hospital (3372, 26 April 2011).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

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

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bray, F.; Laversanne, M.; Weiderpass, E.; Soerjomataram, I. The ever-increasing importance of cancer as a leading cause of premature death worldwide. Cancer 2021, 127, 3029–3030. [Google Scholar] [CrossRef] [PubMed]
  2. Fitzmaurice, C.; Akinyemiju, T.F.; Al Lami, F.H.; Alam, T.; Alizadeh-Navaei, R.; Allen, C.; Alsharif, U.; Alvis-Guzman, N.; Amini, E.; Anderson, B.O.; et al. Global, Regional, and National Cancer Incidence, Mortality, Years of Life Lost, Years Lived with Disability, and Disability-Adjusted Life-Years for 29 Cancer Groups, 1990 to 2016. JAMA Oncol. 2018, 4, 1553. [Google Scholar] [CrossRef] [PubMed]
  3. Sung, H.; Ferlay, J.; Siegel, R.L.; Laversanne, M.; Soerjomataram, I.; Jemal, A.; Bray, F. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J. Clin. 2021, 71, 209–249. [Google Scholar] [CrossRef] [PubMed]
  4. Lemjabbar-Alaoui, H.; Hassan, O.U.; Yang, Y.-W.; Buchanan, P. Lung cancer: Biology and treatment options. Biochim. Biophys. Acta -Rev. Cancer 2015, 1856, 189–210. [Google Scholar] [CrossRef] [Green Version]
  5. Chang, J.Y.; Senan, S.; Paul, M.A.; Mehran, R.J.; Louie, A.V.; Balter, P.; Groen, H.J.M.; McRae, S.E.; Widder, J.; Feng, L.; et al. Stereotactic ablative radiotherapy versus lobectomy for operable stage I non-small-cell lung cancer: A pooled analysis of two randomised trials. Lancet Oncol. 2015, 16, 630–637. [Google Scholar] [CrossRef] [Green Version]
  6. Chi, A.; Fang, W.; Sun, Y.; Wen, S. Comparison of Long-term Survival of Patients with Early-Stage Non–Small Cell Lung Cancer after Surgery vs Stereotactic Body Radiotherapy. JAMA Netw. Open 2019, 2, e1915724. [Google Scholar] [CrossRef] [PubMed]
  7. Razi, S.S.; Kodia, K.; Alnajar, A.; Block, M.I.; Tarrazzi, F.; Nguyen, D.; Villamizar, N. Lobectomy Versus Stereotactic Body Radiotherapy in Healthy Octogenarians with Stage I Lung Cancer. Ann. Thorac. Surg. 2021, 111, 1659–1665. [Google Scholar] [CrossRef] [PubMed]
  8. Schneider, B.J.; Daly, M.E.; Kennedy, E.B.; Antonoff, M.B.; Broderick, S.; Feldman, J.; Jolly, S.; Meyers, B.; Rocco, G.; Rusthoven, C.; et al. Stereotactic Body Radiotherapy for Early-Stage Non–Small-Cell Lung Cancer: American Society of Clinical Oncology Endorsement of the American Society for Radiation Oncology Evidence-Based Guideline. J. Clin. Oncol. 2018, 36, 710–719. [Google Scholar] [CrossRef] [Green Version]
  9. Timmerman, R.D.; Paulus, R.; Pass, H.I.; Gore, E.M.; Edelman, M.J.; Galvin, J.; Straube, W.L.; Nedzi, L.A.; McGarry, R.C.; Robinson, C.G.; et al. Stereotactic Body Radiation Therapy for Operable Early-Stage Lung Cancer. JAMA Oncol. 2018, 4, 1263. [Google Scholar] [CrossRef] [PubMed]
  10. Chang, J.Y.; Liu, H.; Balter, P.; Komaki, R.; Liao, Z.; Welsh, J.; Mehran, R.J.; Roth, J.A.; Swisher, S.G. Clinical outcome and predictors of survival and pneumonitis after stereotactic ablative radiotherapy for stage I non-small cell lung cancer. Radiat. Oncol. 2012, 7, 152. [Google Scholar] [CrossRef] [Green Version]
  11. Onishi, H.; Shirato, H.; Nagata, Y.; Hiraoka, M.; Fujino, M.; Gomi, K.; Karasawa, K.; Hayakawa, K.; Niibe, Y.; Takai, Y.; et al. Stereotactic body radiotherapy (SBRT) for operable stage I non-small-cell lung cancer: Can SBRT be comparable to surgery? Int. J. Radiat. Oncol. Biol. Phys. 2011, 81, 1352–1358. [Google Scholar] [CrossRef]
  12. Timmerman, R. Stereotactic Body Radiation Therapy for Inoperable Early Stage Lung Cancer. JAMA 2010, 303, 1070. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Dong, B.; Zhu, X.; Shu, Z.; Ji, Y.; Lu, F.; Wang, J.; Chen, M. Video-Assisted Thoracoscopic Lobectomy Versus Stereotactic Body Radiotherapy Treatment for Early-Stage Non-Small Cell Lung Cancer: A Propensity Score-Matching Analysis. Front. Oncol. 2020, 10, 585709. [Google Scholar] [CrossRef] [PubMed]
  14. Chang, J.Y.; Mehran, R.J.; Feng, L.; Verma, V.; Liao, Z.; Welsh, J.W.; Lin, S.H.; O’Reilly, M.S.; Jeter, M.D.; Balter, P.A.; et al. Stereotactic ablative radiotherapy for operable stage I non-small-cell lung cancer (revised STARS): Long-term results of a single-arm, prospective trial with prespecified comparison to surgery. Lancet Oncol. 2021, 22, 1448–1457. [Google Scholar] [CrossRef]
  15. Chi, A.; Liao, Z.; Nguyen, N.P.; Xu, J.; Stea, B.; Komaki, R. Systemic review of the patterns of failure following stereotactic body radiation therapy in early-stage non-small-cell lung cancer: Clinical implications. Radiother. Oncol. 2010, 94, 1–11. [Google Scholar] [CrossRef] [PubMed]
  16. Senthi, S.; Lagerwaard, F.J.; Haasbeek, C.J.A.; Slotman, B.J.; Senan, S. Patterns of disease recurrence after stereotactic ablative radiotherapy for early stage non-small-cell lung cancer: A retrospective analysis. Lancet Oncol. 2012, 13, 802–809. [Google Scholar] [CrossRef]
  17. Shultz, D.B.; Trakul, N.; Abelson, J.A.; Murphy, J.D.; Maxim, P.G.; Le, Q.T.; Loo, B.W., Jr.; Diehn, M. Imaging features associated with disease progression after stereotactic ablative radiotherapy for stage I non-small-cell lung cancer. Clin. Lung Cancer 2014, 15, 294–301. [Google Scholar] [CrossRef]
  18. Palma, D.; Senan, S. Stereotactic radiation therapy: Changing treatment paradigms for stage I nonsmall cell lung cancer. Curr. Opin. Oncol. 2011, 23, 133–139. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Samper Ots, P.M.; Vallejo Ocaña, C.; Martin Martin, M.; Celada Álvarez, F.J.; Farga Albiol, D.; Almendros Blanco, P.; Hernandez Machancoses, A.; Rico Oses, M.; Flamarique Andueza, S.; Romero Ruperto, F.; et al. Stereotactic body radiotherapy for early-stage non-small cell lung cancer: A multicentre study by the Oncologic Group for the Study of Lung Cancer (Spanish Radiation Oncology Society). Clin. Transl. Oncol. 2022, 24, 342–349. [Google Scholar] [CrossRef] [PubMed]
  20. Ryuno, Y.; Abe, T.; Iino, M.; Saito, S.; Aoshika, T.; Oota, T.; Igari, M.; Hirai, R.; Kumazaki, Y.; Kaira, K.; et al. High-dose stereotactic body radiotherapy using CyberKnife® for stage I peripheral lung cancer: A single-center retrospective study. Radiat. Oncol. 2022, 17, 128. [Google Scholar] [CrossRef] [PubMed]
  21. Aerts, H.J.; Velazquez, E.R.; Leijenaar, R.T.; Parmar, C.; Grossmann, P.; Carvalho, S.; Bussink, J.; Monshouwer, R.; Haibe-Kains, B.; Rietveld, D.; et al. Decoding tumour phenotype by noninvasive imaging using a quantitative radiomics approach. Nat. Commun. 2014, 5, 4006. [Google Scholar] [CrossRef] [PubMed]
  22. Gillies, R.J.; Kinahan, P.E.; Hricak, H. Radiomics: Images Are More than Pictures, They Are Data. Radiology 2016, 278, 563–577. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Diehn, M.; Nardini, C.; Wang, D.S.; McGovern, S.; Jayaraman, M.; Liang, Y.; Aldape, K.; Cha, S.; Kuo, M.D. Identification of noninvasive imaging surrogates for brain tumor gene-expression modules. Proc. Natl. Acad. Sci. USA 2008, 105, 5213–5218. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Segal, E.; Sirlin, C.B.; Ooi, C.; Adler, A.S.; Gollub, J.; Chen, X.; Chan, B.K.; Matcuk, G.R.; Barry, C.T.; Chang, H.Y.; et al. Decoding global gene expression programs in liver cancer by noninvasive imaging. Nat. Biotechnol. 2007, 25, 675–680. [Google Scholar] [CrossRef] [PubMed]
  25. Zambirinis, C.P.; Midya, A.; Chakraborty, J.; Chou, J.F.; Zheng, J.; McIntyre, C.A.; Koszalka, M.A.; Wang, T.; Do, R.K.; Balachandran, V.P.; et al. Recurrence after Resection of Pancreatic Cancer: Can Radiomics Predict Patients at Greatest Risk of Liver Metastasis? Ann. Surg. Oncol. 2022, 29, 4962–4974. [Google Scholar] [CrossRef]
  26. Zhu, W.; Huang, X.; Qi, Q.; Wu, Z.; Min, X.; Zhou, A.; Xu, P. Artificial Neural Network-Based Ultrasound Radiomics Can Predict Large-Volume Lymph Node Metastasis in Clinical N0 Papillary Thyroid Carcinoma Patients. J. Oncol. 2022, 2022, 1–11. [Google Scholar] [CrossRef] [PubMed]
  27. Bortolotto, C.; Lancia, A.; Stelitano, C.; Montesano, M.; Merizzoli, E.; Agustoni, F.; Stella, G.; Preda, L.; Filippi, A.R. Radiomics features as predictive and prognostic biomarkers in NSCLC. Expert Rev. Anticancer Ther. 2021, 21, 257–266. [Google Scholar] [CrossRef] [PubMed]
  28. Dissaux, G.; Visvikis, D.; Da-Ano, R.; Pradier, O.; Chajon, E.; Barillot, I.; Duvergé, L.; Masson, I.; Abgral, R.; Santiago Ribeiro, M.-J.; et al. Pretreatment 18F-FDG PET/CT Radiomics Predict Local Recurrence in Patients Treated with Stereotactic Body Radiotherapy for Early-Stage Non–Small Cell Lung Cancer: A Multicentric Study. J. Nucl. Med. 2020, 61, 814–820. [Google Scholar] [CrossRef]
  29. Franceschini, D.; Cozzi, L.; De Rose, F.; Navarria, P.; Fogliata, A.; Franzese, C.; Pezzulla, D.; Tomatis, S.; Reggiori, G.; Scorsetti, M. A radiomic approach to predicting nodal relapse and disease-specific survival in patients treated with stereotactic body radiation therapy for early-stage non-small cell lung cancer. Strahlenther. Und Onkol. 2020, 196, 922–931. [Google Scholar] [CrossRef]
  30. Huynh, E.; Coroller, T.P.; Narayan, V.; Agrawal, V.; Hou, Y.; Romano, J.; Franco, I.; Mak, R.H.; Aerts, H.J. CT-based radiomic analysis of stereotactic body radiation therapy patients with lung cancer. Radiother. Oncol. 2016, 120, 258–266. [Google Scholar] [CrossRef] [PubMed]
  31. Mattonen, S.A.; Palma, D.A.; Haasbeek, C.J.A.; Senan, S.; Ward, A.D. Early prediction of tumor recurrence based on CT texture changes after stereotactic ablative radiotherapy (SABR) for lung cancer. Med. Phys. 2014, 41, 033502. [Google Scholar] [CrossRef] [PubMed]
  32. Oikonomou, A.; Khalvati, F.; Tyrrell, P.N.; Haider, M.A.; Tarique, U.; Jimenez-Juan, L.; Tjong, M.C.; Poon, I.; Eilaghi, A.; Ehrlich, L.; et al. Radiomics analysis at PET/CT contributes to prognosis of recurrence and survival in lung cancer treated with stereotactic body radiotherapy. Sci. Rep. 2018, 8, 4003. [Google Scholar] [CrossRef] [PubMed]
  33. Starkov, P.; Aguilera, T.A.; Golden, D.I.; Shultz, D.B.; Trakul, N.; Maxim, P.G.; Le, Q.T.; Loo, B.W.; Diehn, M.; Depeursinge, A.; et al. The use of texture-based radiomics CT analysis to predict outcomes in early-stage non-small cell lung cancer treated with stereotactic ablative radiotherapy. Br. J. Radiol. 2019, 92, 20180228. [Google Scholar] [CrossRef]
  34. Van Griethuysen, J.J.M.; Fedorov, A.; Parmar, C.; Hosny, A.; Aucoin, N.; Narayan, V.; Beets-Tan, R.G.H.; Fillion-Robin, J.-C.; Pieper, S.; Aerts, H.J.W.L. Computational Radiomics System to Decode the Radiographic Phenotype. Cancer Res. 2017, 77, e104–e107. [Google Scholar] [CrossRef] [Green Version]
  35. Ernani, V.; Appiah, A.K.; Marr, A.; Zhang, C.; Zhen, W.; Smith, L.M.; Ganti, A.K. Adjuvant Systemic Therapy in Patients with Early-Stage NSCLC Treated With Stereotactic Body Radiation Therapy. J. Thorac. Oncol. 2019, 14, 475–481. [Google Scholar] [CrossRef] [Green Version]
  36. Antonia, S.J.; Villegas, A.; Daniel, D.; Vicente, D.; Murakami, S.; Hui, R.; Kurata, T.; Chiappori, A.; Lee, K.H.; De Wit, M.; et al. Overall Survival with Durvalumab after Chemoradiotherapy in Stage III NSCLC. N. Engl. J. Med. 2018, 379, 2342–2350. [Google Scholar] [CrossRef]
  37. Spigel, D.R.; Faivre-Finn, C.; Gray, J.E.; Vicente, D.; Planchard, D.; Paz-Ares, L.; Vansteenkiste, J.F.; Garassino, M.C.; Hui, R.; Quantin, X.; et al. Five-Year Survival Outcomes from the PACIFIC Trial: Durvalumab After Chemoradiotherapy in Stage III Non–Small-Cell Lung Cancer. J. Clin. Oncol. 2022, 40, 1301–1311. [Google Scholar] [CrossRef]
  38. Chi, A.; Nguyen, N.P. Rationale for Combing Stereotactic Body Radiation Therapy with Immune Checkpoint Inhibitors in Medically Inoperable Early-Stage Non-Small Cell Lung Cancer. Cancers 2022, 14, 3144. [Google Scholar] [CrossRef]
  39. Yu, W.; Tang, C.; Hobbs, B.P.; Li, X.; Koay, E.J.; Wistuba, I.I.; Sepesi, B.; Behrens, C.; Rodriguez Canales, J.; Parra Cuentas, E.R.; et al. Development and Validation of a Predictive Radiomics Model for Clinical Outcomes in Stage I Non-small Cell Lung Cancer. Int. J. Radiat. Oncol. Biol. Phys. 2018, 102, 1090–1097. [Google Scholar] [CrossRef]
Figure 1. Impact of the radiomic factor named “LargeAreaEmphasis_LHH” on overall survival in the validation cohort. Patients were divided into two groups according to the median value of “LargeAreaEmphasis_LHH” (low value group and high value group). The univariate analysis using the log-rank test was conducted (p = 0.044, 5-year OS of 70.7% vs. 50.3%).
Figure 1. Impact of the radiomic factor named “LargeAreaEmphasis_LHH” on overall survival in the validation cohort. Patients were divided into two groups according to the median value of “LargeAreaEmphasis_LHH” (low value group and high value group). The univariate analysis using the log-rank test was conducted (p = 0.044, 5-year OS of 70.7% vs. 50.3%).
Cancers 14 03859 g001
Figure 2. Boxplots of the estimated OS time and the real OS time in the validation cohort. The estimated survival time was calculated by the formula from the multiple linear regression analysis. with the real survival time in the validation cohort. The average estimated OS time was 33.7 months, and the average real OS time was 37.8 months (p = 0.11 by the paired t-test).
Figure 2. Boxplots of the estimated OS time and the real OS time in the validation cohort. The estimated survival time was calculated by the formula from the multiple linear regression analysis. with the real survival time in the validation cohort. The average estimated OS time was 33.7 months, and the average real OS time was 37.8 months (p = 0.11 by the paired t-test).
Cancers 14 03859 g002
Table 1. Patient characteristics of the training and the validation cohort.
Table 1. Patient characteristics of the training and the validation cohort.
FactorTrainingValidationp Value
NRateNRate
N250108
Histology
Adenocarcinoma4718.8%1816.7%0.98
Squamous cell carcinoma104.0%43.7%
Others156.0%76.5%
None17871.2%7973.1%
Location
Central3915.6%1211.1%0.26
Peripheral21184.4%9688.9%
Site
Right upper lobe7228.8%4037.0%0.31
Right middle lobe208.0%76.5%
Right lower lobe5722.8%2119.4%
Left upper lobe6526.0%2018.5%
Left lower lobe3614.4%2018.5%
C/T ratio
≥50%20080.0%8982.4%0.66
<50%5020.0%1917.6%
Maximum tumor diameter
Mean (mm)20.520.80.85
SD10.59.3
Age
Mean (y.o.)77.177.80.47
SD8.98.1
KPS
Mean (%)89.088.50.52
SD5.86.4
SUV-max
Mean5.95.50.59
SD5.14.2
Total radiation dose
Mean (Gy)52.551.50.11
SD4.67.6
Sex
Male16967.6%7670.4%0.60
Female8132.4%3229.6%
RT technique
VMAT17168.4%7064.8%0.23
Fixed multi-port7530.0%3330.6%
Both41.6%54.6%
Salvage treatment
Systemic therapy156.0%98.3%0.45
Radiotherapy93.6%32.8%
Chemoradiotherapy31.2%00.0%
Surgery31.2%00.0%
Other treatment20.8%00.0%
None2911.6%1110.2%
Unknown10.4%32.8%
No recurrence18875.20%8275.9%
C/T―consolidation/tumor; KPS―Karnofsky performance status; SUV-max―maximum standard uptake value; RT―radiotherapy; VMAT―volumetric modulated arc therapy.
Table 2. Univariate analysis of background factors using the log-rank test.
Table 2. Univariate analysis of background factors using the log-rank test.
FactorOSLRFSPFS
5-Year
OS (%)
95%CIp Value5-Year
LRFS (%)
95%CIp Value5-Year
PFS (%)
95%CIp Value
KPS
Low value35.516.2–55.50.002863.640.4–79.80.001250.130.1–67.10.00053
High value61.551.5–70.0 79.971.6–85.9 65.956.4–73.8
SUV-max
Low value61.643.8–75.30.1179.764.7–88.90.009053.135.9–67.60.11
High value41.723.9–58.6 60.142.8–73.8 47.031.5–61.0
RT technique
VMAT60.548.3–70.60.5277.366.9–84.80.7258.847.5–68.50.26
Fixed multi-port57.042.5–69.1 76.262.7–85.3 70.255.8–80.7
Total radiation dose
Low value53.940.8–65.30.8070.357.6–79.90.07458.145.0–69.10.28
High value64.452.4–74.0 82.572.2–89.3 68.357.4–77.0
Age
Low value63.451.8–73.00.3379.168.7–86.30.7063.851.9–73.40.60
High value48.733.5–62.4 75.162.0–84.2 63.651.0–73.8
Histology
Adenocarcinoma64.639.4–81.40.2569.948.1–83.90.8768.747.5–82.70.73
Squamous cell carcinoma90.047.3–98.5 85.733.4–97.9 77.134.5–93.9
Others43.815.7–69.1 82.144.4–95.3 65.031.0–85.4
None54.943.8–64.7 78.769.7–85.3 60.850.1–69.9
Maximum tumor diameter
Low value66.654.2–76.40.04288.379.5–93.50.001768.656.3–78.10.12
High value48.034.6–60.1 64.851.7–75.1 56.944.3–67.7
Sex
Male52.140.8–62.10.03070.459.9–78.60.01159.948.9–69.20.26
Female68.251.5–80.1 90.379.0–95.6 69.454.7–80.2
Location
Central22.15.9–44.60.0004058.134.0–76.10.0006326.25.8–53.20.0011
Peripheral63.954.4–72.0 80.672.4–86.6 69.060.5–76.1
Site
Right upper lobe66.646.0–80.90.8874.257.5–85.10.4855.740.0–68.70.49
Right middle lobe48.217.5–73.7 76.648.8–90.5 64.530.4–85.1
Right lower lobe51.533.1–67.2 74.655.4–86.5 68.249.7–81.1
Left upper lobe62.445.3–75.5 85.367.9–93.7 70.252.2–82.6
Left lower lobe53.530.4–72.0 75.454.1–87.8 60.735.8–78.5
C/T ratio
≥50%57.847.6–66.60.6265.646.2–79.40.04252.634.1–68.10.11
<50%58.235.9–75.0 80.672.1–86.7 66.356.6–74.4
Salvage treatment
Systemic therapy40.410.6–69.40.005411.30.7–38.8<0.0001N/AN/A<0.0001
Radiotherapy23.71.0–63.8 N/AN/A N/AN/A
ChemoradiotherapyN/AN/A N/AN/A N/AN/A
Surgery66.75.4–94.5 66.75.4–94.5 N/AN/A
Other treatmentN/AN/A N/AN/A N/AN/A
NoneN/AN/A N/AN/A N/AN/A
UnknownN/AN/A N/AN/A N/AN/A
No recurrence65.955.6–74.3 99,294.7–99.9 100.0100.0–100.0
OS―overall survival; LRFS―local relapse-free survival; PFS―progression-free survival; KPS―Karnofsky performance status; SUV-max―maximum standard uptake value; VMAT―volumetric modulated arc therapy; N/A―not applicable.
Table 3. Multivariate analysis of radiomic factors using the Cox proportional hazards regression model with the stepwise method.
Table 3. Multivariate analysis of radiomic factors using the Cox proportional hazards regression model with the stepwise method.
FactorHazard RatioLower 95%CIUpper 95%CIp ValueCoefficient
OS
90 Percentile_HHH0.51970.32230.83810.007267−0.6545
LargeAreaEmphasis_LHH1.79901.12602.87500.0141200.5871
Mean_HHH1.97301.22403.17900.0052550.6795
Median_HLL1.98401.24003.17400.0042620.6852
LRFS
InverseVariance_HLL3.54501.645007.6380.00123101.2655
SmallDependenceHighGrayLevelEmphasis_HHH0.21330.092020.49420.0003138−1.5453
PFS
SmallDependenceHighGrayLevelEmphasis_HHH0.37150.20650.66860.0009571−0.9901
TotalEnergy_HHL2.46101.39104.35500.00198400.9006
JointEntropy_HLL0.49690.29430.83880.0088420−0.6994
OS―overall survival; LRFS―local relapse-free survival; PFS―progression-free survival.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sawayanagi, S.; Yamashita, H.; Nozawa, Y.; Takenaka, R.; Miki, Y.; Morishima, K.; Ueno, H.; Ohta, T.; Katano, A. Establishment of a Prediction Model for Overall Survival after Stereotactic Body Radiation Therapy for Primary Non-Small Cell Lung Cancer Using Radiomics Analysis. Cancers 2022, 14, 3859. https://doi.org/10.3390/cancers14163859

AMA Style

Sawayanagi S, Yamashita H, Nozawa Y, Takenaka R, Miki Y, Morishima K, Ueno H, Ohta T, Katano A. Establishment of a Prediction Model for Overall Survival after Stereotactic Body Radiation Therapy for Primary Non-Small Cell Lung Cancer Using Radiomics Analysis. Cancers. 2022; 14(16):3859. https://doi.org/10.3390/cancers14163859

Chicago/Turabian Style

Sawayanagi, Subaru, Hideomi Yamashita, Yuki Nozawa, Ryosuke Takenaka, Yosuke Miki, Kosuke Morishima, Hiroyuki Ueno, Takeshi Ohta, and Atsuto Katano. 2022. "Establishment of a Prediction Model for Overall Survival after Stereotactic Body Radiation Therapy for Primary Non-Small Cell Lung Cancer Using Radiomics Analysis" Cancers 14, no. 16: 3859. https://doi.org/10.3390/cancers14163859

APA Style

Sawayanagi, S., Yamashita, H., Nozawa, Y., Takenaka, R., Miki, Y., Morishima, K., Ueno, H., Ohta, T., & Katano, A. (2022). Establishment of a Prediction Model for Overall Survival after Stereotactic Body Radiation Therapy for Primary Non-Small Cell Lung Cancer Using Radiomics Analysis. Cancers, 14(16), 3859. https://doi.org/10.3390/cancers14163859

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop