Skip to main content
Download PDF

CT radiomics-based long-term survival prediction for locally advanced non-small cell lung cancer patients treated with concurrent chemoradiotherapy using features from tumor and tumor organismal environment

Radiation Oncology volume 17, Article number: 184 (2022) Cite this article

Abstract

Background

Definitive concurrent chemoradiotherapy (CCRT) is the standard treatment for locally advanced non-small cell lung cancer (LANSCLC) patients, but the treatment response and survival outcomes varied among these patients. We aimed to identify pretreatment computed tomography-based radiomics features extracted from tumor and tumor organismal environment (TOE) for long-term survival prediction in these patients treated with CCRT.

Methods

A total of 298 eligible patients were randomly assigned into the training cohort and validation cohort with a ratio 2:1. An integrated feature selection and model training approach using support vector machine combined with genetic algorithm was performed to predict 3-year overall survival (OS). Patients were stratified into the high-risk and low-risk group based on the predicted survival status. Pulmonary function test and blood gas analysis indicators were associated with radiomic features. Dynamic changes of peripheral blood lymphocytes counts before and after CCRT had been documented.

Results

Nine features including 5 tumor-related features and 4 pulmonary features were selected in the predictive model. The areas under the receiver operating characteristic curve for the training and validation cohort were 0.965 and 0.869, and were reduced by 0.179 and 0.223 when all pulmonary features were excluded. Based on radiomics-derived stratification, the low-risk group yielded better 3-year OS (68.4% vs. 3.3%, p < 0.001) than the high-risk group. Patients in the low-risk group had better baseline FEV1/FVC% (96.3% vs. 85.9%, p = 0.046), less Grade ≥ 3 lymphopenia during CCRT (63.2% vs. 83.3%, p = 0.031), better recovery of lymphopenia from CCRT (71.4% vs. 27.8%, p < 0.001), lower incidence of Grade ≥ 2 radiation-induced pneumonitis (31.6% vs. 53.3%, p = 0.040), superior tumor remission (84.2% vs. 66.7%, p = 0.003).

Conclusion

Pretreatment radiomics features from tumor and TOE could boost the long-term survival forecast accuracy in LANSCLC patients, and the predictive results could be utilized as an effective indicator for survival risk stratification. Low-risk patients might benefit more from radical CCRT and further adjuvant immunotherapy.

Trial registration:

retrospectively registered.

Introduction

Definitive concurrent chemoradiotherapy (CCRT) is the standard treatment for patients with unresectable locally advanced non-small cell lung cancer (LANSCLC). In the past two decades, concomitant regimens achieved promising tumor local control and long-term survival. With improved outcome, the maintenance of an adequate pulmonary function is essential to ensure acceptable quality of life and adjuvant immunotherapy. However, many patients with LANSCLC are diagnosed with pre-existing lung comorbidities, which significantly increases the risk for radiation-induced lung toxicity (RILT) [1, 2].

Most existing RILT prediction models largely focused on clinical prognostic factors (CPFs) and dose-volume histogram parameters [3,4,5], but remained insufficient. Recently, machine learning methods have been reported to improve the capacity of the predictive modelling [6,7,8,9], compared with logistic regression widely used in normal tissue complication probability model.

Moreover, radiomics analysis, attempting to identify computational biomarkers potentially hidden within high-throughput imaging data [10, 11], has been demonstrated the added predictive value for overall survival (OS) [12,13,14] or RILT [8, 9]. However, most of them rely on the radiomic information from tumor or its surrounding peritumoral region, few studies have been designed based on the radiomics analysis of tumor organismal environment (TOE).

Similar to other published reports [15, 16], our previous study [17] indicated that pulmonary function test (PFT) was significantly related to patients’ long-term survival. However, it failed to predict progression-free survival (PFS). Even though patients with worse FEV1/FVC% or DLCO% showed a high objective response rate (ORR) to CCRT, their survival outcomes were still poor, hinting that TOE, the status of lungs in the case of LANSCLC might play an indispensable role in the prognostic prediction after CCRT. As some patients could not tolerate well with PFT, radiomics analysis using machine learning method might be an effective technique to investigate the relationship of tumor and TOE, due to its accessibility.

In this study, we utilized computed tomography (CT) images before CCRT to develop an image-based machine learning framework to analyze the relationship of primary lung tumor and bilateral lungs for long-term survival prediction in LANSCLC. To balance the training accuracy and predictive capability using relative small number of patient samples, an integrated feature selection and model training (IFSMT) approach was developed to extract the most critical quantitative radiomic features from both tumor and lungs. A radiomic-based risk stratification was built to distinguish high-risk and low-risk patients and provided evidence for clinical decision making.

Methods

Study population

Consecutive patients irradiated for lung cancer from September 2011 to April 2019 in our institution were retrospectively screened. Inclusion criteria included: (1) histologically confirmed NSCLC; (2) unresectable stage III disease (AJCC/UICC 8th staging criteria) proven by chest and upper abdominal CT, brain magnetic resonance imaging (MRI), bone scan and/or positron emission tomography-computed tomography (PET-CT); (3) definitive radiotherapy with concurrent chemotherapy was administered; (4) stay followed-up no less than 6 months since the start of radiotherapy (unless death or disease progression was documented); (5) complete clinical records. Patients that met the inclusion criteria were randomly assigned into the training and validation cohort, with the numbers at a ratio 2:1.

Planning CT image acquisition

The four-dimensional (4D) planning CT scan was performed 1–2 weeks prior to treatment, using multiple CT simulation positioning machines with varied parameter settings in our institution (detailed in Additional File 1). Ten phases of the breathing cycle were reconstructed, including: 0%, 10%, 20%, 30%, 40%, 50%, 60%, 70%, 80%, and 90%. The segmentation and radiomics were then performed on the 20% phase (middle exhale phase) with a consistent mediastinum/lung window level setting.

Radiotherapy and concurrent chemotherapy

Patients were positioned supine and immobilized in a vacuum pad. They were scanned from the Atlas to the second lumbar vertebra level with 0.3-0.5 cm thickness slices to obtain the stimulation CT images. The respiration motion was recorded by performing 4DCT scanning. The maximum intensity projection images were reconstructed using the images collected in 10 phases of respiratory cycle. Gross tumor volume (GTV) was delineated to cover the tumor and involved regional nodes visible on each phase of the 4DCT. The total volumes of GTVs across the 10 respiratory phases CT composed the internal target volume (ITV). Planning target volumes (PTVs) were created by expanding GTV and clinical target volume with 6 mm. Lungs were delineated according to the atlases for organs at risk (OARs) in thoracic radiation therapy [18], but GTV was excluded from the lung delineation. A dose of 60-76 Gy was prescribed to PTV-GTV in 22–33 fractions, with 2-3 Gy per fraction performed once daily, using intensity modulated radiation therapy technique. The dose constraints for OARs were: V20 < 35% for lungs; mean lung dose < 19 Gy; maximun dose (Dmax) of esophagus < 105% prescription dose; Dmax of spinal cord < 46 Gy; V30 < 40% for heart.

All patients received platinum-based double agents weekly or every three weeks. The regimens included docetaxel/paclitaxel/etopside/pemetrexed plus platinum.

Evaluation and follow-up

The baseline characteristics of each patient before entry were reviewed attentively and extracted from their medical records, including blood tests, PFT, blood gas analysis (BGA) and radiologic tests. All included patients received regular radiologic follow-up, including chest and upper abdominal CT and brain MRI performed every 3 ~ 6 months in the first 2 years, and every 6 ~ 12 months thereafter. PET-CT, bone scan, and biopsy were recommended if clinically required. The responses to CCRT were first assessed by an independent radiation oncologist and confirmed by a senior physician at 4 ~ 6 weeks post CCRT, based on Response Evaluation Criteria in Solid Tumors 1.1. Another senior radiologist was consulted for disagreement. Therapeutic toxicities were graded and recorded according to Common Terminology Criteria for Adverse Events 4.0.

OS modelling procedures

The whole procedures were illustrated in Fig. 1. For both cohorts of patients, the regions of interest (ROIs) corresponding to GTV and lungs were delineated by an auto-contouring software tool CezanneDraw™ v1.0 (Homology Medical, Ningbo, China, 2020) using the CT slices and manually modified by radiation oncologists if necessary. One 3D bounding box was fitted for each ROI. And inside the bounding box, the CT values of the ROI voxels were retained while the values of other voxels were marked by zero. CT values of voxels in each bounding box were then interpolated to a resolution of 1 mm×1 mm×5 mm and resampled into 400 discrete values (called bins) with absolute discretization from − 1000 to 3000 Hounsfield units, leading to a fixed bin size of 10 Hounsfield units.

Fig. 1
figure 1

Schematic overview of the integrated feature selection and model training (IFSMT) approach. IFSMT approach consisted of five steps: (1) volumetric data pre-processing; (2) delineation; (3) feature extraction; (4) integrated feature selection and model training; (5) model validation using leave-one-out cross-validation (LOOCV).

Full size image

A total of 92 tumor-related and lung-related features were then computed for both ROIs and used as the input feature pool for the machine learning framework by the LIFEx software (version 3.44) [19]. The imaging-based features covered two categories of texture features and first order features. The texture features consisted of four sub-categories of matrix based texture features. These matrices included the grey-level co-occurrence matrix (GLCM), neighborhood grey-level different matrix (NGLDM), grey-level run length matrix (GLRLM) and grey-level zone length matrix (GLZLM). The first order features included indices from shape, indices from histogram and conventional indices.

The machine learning based classification method used to predict the two-class 3-year survival status for each individual patient was support vector machine (SVM) [20]. The SVM mapped the features of training data into a high-dimensional feature space through a kernel function and utilizes a hyper-plane to optimally separate the training data points into two categories. To reduce the possibility of overfitting, only a subset of features from the feature pool could be selected for the input of SVM. In this study, the IFSMT approach was developed to maximize the fitting accuracy and minimize the overfitting potential. This posteriori approach applied the genetic algorithm (GA) for the feature selection, which was illustrated in Fig. 2 and Additional File 2. A chromosome represents a feature template working with SVM of certain configuration for diagnosing purpose. The SVM is implemented in leave-one-out cross-validation (LOOCV) fashion to score a chromosome. In each generation, the chromosomes of higher scores may go through mutation, partially changing feature encoding, and crossover, partially exchanging feature encoding, to make new ones to replace those of lower scores. Collect the chromosome of best score from each generation into a group. And the best one in the group is the result of the model. Manual reconfiguration of SVM is not included in the model.

Fig. 2
figure 2

Schematic overview of the genetic algorithm (GA) in the integrated feature selection and model training (IFSMT) approach. A chromosome is scored with LOOCV-SVM. The chromosomes of higher scores may go through mutation and crossover to make new ones to replace those of lower scores. Collect the chromosome of best score from each generation into a group. And the best one in the group is the result of the model. Abbreviation: LOOCV, leave-one-out cross-validation; SVM, support vector machine

Full size image

Once the optimal set of features was determined, the SVM models were trained again on the training cohort. In this study, after extensive experimental comparisons, the linear kernel was chosen for SVM and optimal hyper parameters of the SVM (C, ε and γ) were determined through exhaustive search in the parametric space. Receiver operating characteristics (ROC) curves were obtained by varying threshold of the decision variable, the signed distance to decision hyper-plane. Area under curve (AUC) for each ROC was calculated for training cohort. The trained models were then used to predict the survival status for each individual patient in the validation cohort, and ROCs and their corresponding AUCs were also calculated. All the above feature selection and machine learning approaches were implemented on the cloud-based clinical data service platform iRAAS® v2.0 (Homology Medical, Ningbo, China, 2020).

To assess the importance of each selected feature to the accurate prediction of the clinical outcome, a one-by-one feature evaluation procedure was designed. This procedure tested the importance of each feature by deleting each feature from the selected feature set and calculating the reduction of the AUC for the model trained with the original selected features except this specific feature. This reduction of model performance was used as the importance weight (IW) of this feature. All the selected features were then sorted according to their IWs. To further assess the importance organismal features, the AUC for the model trained with the original selected features excluding all the lung-related features were also calculated.

Statistical methods

OS was defined as the time from radiotherapy start to the last follow up, which ended at November 30th, 2021, or death. A t-test was used to determine if there was significant difference between the means of continuous variables, while Fisher’s exact test was performed to reveal the difference in distribution between two groups of categories variables. The association between radiomic features and PFT/BGA indicators was examined using Pearson’s correlation coefficient. A p-value < 0.05 (two-sided) were considered as statistically significant. Missing data were excluded from the statistical analysis. Statistics were performed using SPSS 22.0 (IBM, Chicago, IL, USA).

To report the model fitting accuracy and the prediction capability, the true positive rate (TPR), true negative rate (TNR), F1 score, overall prediction accuracy, average prediction accuracy for the training cohort and validation cohort were calculated based on the SVM model. Herein, death is marked as the positive. The overall prediction accuracy was expressed as the number correctly predicted patients / the number of all patients; and the average prediction accuracy = (TPR + TNR)/2.

To assess the prognostic value of the survival status model, the predicted 3-year survival status was adopted respectively as the clinical risk estimator to stratify the patients into the high-risk and low-risk groups. Patients with negative predicted survival status were classified into the low-risk group and the others with positive predicted survival status into the high-risk group. Kaplan-Meier curves for both groups were displayed to illustrate its effectiveness and log-rank test was performed.

Results

Patient characteristics

A total of 298 LANSCLC patients were included for analysis, with 200 in the training cohort and 98 in the validation cohort. The baseline and treatment-related characteristics were comparable between these two cohorts (Additional File 3). There were 57 females and 241 males in the whole cohort, with the median age of 59 years (range, 28–81 years). Squamous cell carcinoma was the predominant histologic type both in the training (46.5%) and validation (62.2%) cohorts.

OS modelling

With the median follow-up of 27.7 (range, 4.0 ~ 122.7) months for all and 67.0 months (range, 36.2 ~ 122.7 months) for event-free patients, our cohort demonstrated the estimated median OS of 27.6 (95% confidence interval (CI), 22.3 ~ 33.0) months, and the 3-year OS rate was 43.0% (95%CI, 37.3%~48.7%).

As shown in Table 1, the overall prediction accuracy for 3-year survival status was 92.50% and 85.71%, and the AUC of the ROC was 0.965 and 0.869, respectively, in the training and validation cohort.

Table 1 OS Training and validation accuracy using all selected features or without pulmonary features
Full size table

Stratification of patients in the validation cohort with machine learning model

In the validation cohort, 60 (61.2%) of 98 patients were stratified into the high-risk group and 38 (44.1%) into the low-risk group. CCRT was more successful in patients in the low-risk group than those in the high-risk group. The ORR was 84.2% (32/38) and 66.7% (40/60) in the low-risk and high-risk group, respectively (p = 0.003) (Additional File 4). And the low-risk group yielded better 3-year OS (68.4% versus 3.3%, p < 0.001, log-rank) than the high-risk group (Fig. 3B). What’s more, the rate of Grade ≥ 2 pneumonitis was 31.6% (12/38), versus 53.3% (32/60) (p = 0.040) in the low-risk and high-risk group. The typical presentation of two patients in the low-risk and high-risk group was illustrated in Fig. 4.

Fig. 3
figure 3

Kaplan-Meier curves for the training and validation cohort, with all selected features (a, b), and without pulmonary features (c, d), respectively

Full size image
Fig. 4
figure 4

Two patients from the high-risk (A) and low-risk group (B). From the series CT images, there were discernible distinction observed in tumor and pulmonary status between the two cases. High-risk patient had heterogeneous primary lung tumor and chronic obstructive pneumonia (A-1), while low-risk patient had relatively homogeneous primary lung tumor and better pulmonary condition (B-1). Significant tumor remission was achieved in low-risk patient without obvious radiation pneumonitis after CCRT (B-2), while high-risk patient had stable disease and developed Grade 2 radiation pneumonitis in bilateral lungs (A-2). The GLZLM matrices for tumor, GLRLM matrices along + x axis for lungs, and histograms of HU values for tumor (A-3, B-3) were displayed. It was found that the short homogeneous runs and the non-uniformity of the grey-levels (CT value) were emphasized for high-risk patient compared to those of low-risk patient

Full size image

Correlation of selected radiomic features to the model performance

A total of 9 features were selected in the proposed model, including 5 tumor-related features and 4 lung-related features. In Table 2, the IW of each selected feature for both training and validation cohorts were listed in the order from high to low. The imaging features from lungs ranked at 2nd, 4th, 5th and 8th in the all 9 features in the training cohort, and 1st, 3rd, 6th, and 8th in the validation cohort. When all pulmonary features were excluded from the selected feature set, the AUCs for the training and validation cohorts were reduced by 0.179 and 0.223, respectively (Fig. 5). Figure 4 showed two patients in the low-risk and high-risk groups.

Table 2 Selected features and their importance rank for the training and validation cohort
Full size table
Fig. 5
figure 5

The receiver operating characteristics (ROC) curves for the training and validation cohort, with all selected features (a, b), and without pulmonary features (c, d), respectively

Full size image

Correlation of radiomic features to the PFT/BGA indicators

Pearson correlation analysis (Fig. 6A) demonstrated that FEV1/FVC% had modest correlation with three pulmonary features (SHAPE_Volume_mL, GLRLM_LRE and GLRLM_RP) (all Pearson correlation >|0.45|), and mild correlation with CONV_SUVstd of lungs and GLZLM_GLNUz of tumor (all Pearson correlation >|0.25|). Other PFT/BGA indicators and radiomic features were not well correlated.

Patients in the low-risk group had better baseline FEV1/FVC% (median, 96.3% vs. 85.9%, p = 0.046) compared with those in the high-risk group (Fig. 6B). Kaplan-Meier analysis indicated that better baseline FEV1/FVC% (p = 0.006) and SaO2 (p = 0.039) could exhibit superior OS, DLCO% (p = 0.063) had a tendency to be associated with OS, however, pO2 (p = 0.110) and AaDO2 (p = 0.299) failed to predict OS (Additional File 5).

Fig. 6
figure 6

Pearson correlation coefficient heatmap for radiomic features and PFT/BGA indicators (A), and the distribution of PFT/BGA indicators between the low-risk and high-risk group (B-E). A. The FEV1/FVC% had modest correlation with three pulmonary features (SHAPE_Volume_mL, GLRLM_LRE and GLRLM_RP) (all Pearson correlation >|0.45|), and mild correlation with CONV_SUVstd of lungs and GLZLM_GLNUz of tumor (all Pearson correlation >|0.25|). Other PFT/BGA indicators and radiomic features were not well correlated. B-E. Patients in the low-risk group had better baseline FEV1/FVC% (median, 96.3% vs. 85.9%, p = 0.046) compared with those in the high-risk group. No significant difference of DLCO% (median, 84.3% vs. 77.5%, p = 0.136), pO2 (median, 86.0 vs. 86.0 mmHg, p > 0.999), AaDO2 (median, 18.0 vs. 18.0 mmHg, p > 0.999) and SaO2 (median, 96.8% vs. 96.5%, p = 0.634) was found. Abbreviation: PFT, pulmonary function test; BGA, blood gas analysis

Full size image

Dynamic changes of lymphocyte counts before and after CCRT

Although there was no significant difference in lymphocyte counts before CCRT (median, 1650 vs. 1650 cells/mm3, p > 0.99) between the low-risk and high-risk group (Additional File 6), patients in the low-risk group had less Grade ≥ 3 lymphopenia (63.2% vs. 83.3%, p = 0.031) during CCRT, and more patients in the low-risk group could recover to normal level (≥ 1000 cells/mm3) at 4 ~ 6 weeks post CCRT (71.4% vs. 27.8%, p < 0.001). Kaplan–Meier curves demonstrated that better recovery to normal level (≥ 1000 cells/mm3) at 4 ~ 6 weeks post CCRT (3-year OS rate, 47.5% versus 14.3%, p = 0.001) was a prognostic factor of OS (Additional File 5).

Discussions

Application of radiomics to the long-term survival prediction for LANSCLC after CCRT is a reasonable extension under the background of the field-wide adoption of machine learning methods. Other than previous works focused on the features from tumor and peritumoral tissue, the relationship between tumor and TOE is increasingly attached importance. Significant association was found between pulmonary function and radiomic features extracted from the lungs of CT images [21,22,23]. In current study, the long-term survival forecast accuracy of LANCLC patients after CCRT was demonstrated to be boosted by integrating primary tumor characteristics and pulmonary features from pretreatment CT images. Based on the CT-based predictive model, patients could be precisely stratified into the low-risk and high-risk group before treatment, which should be considered in individualized treatment decision-making process.

From the importance rank of the selected features, it could be confirmed that two features from tumor, GLRLM_SRE and GLZLM_GLNUz which represent the inhomogeneity of CT images [19], remained important factors determining OS, which were consistent with published literatures [24, 25]. Meanwhile, the ranking of pulmonary features underlined their indispensable role in the OS forecast. Our results of the significant difference between fitting and prediction accuracies with and without pulmonary features in model performance further support this finding, implying that the TOE, herein the pulmonary environment, might have a significant impact in LANSCLC patients with large tumor burden and limited pulmonary function. Accordingly, the relatively longer OS for patients with healthier pulmonary status could possibly contribute to their more tolerance to radical CCRT and less incidence of severe lung toxicities.

PFT have been reported to predict the risk of RILT after CCRT [26,27,28,29]. Our previous work showed that FEV1/FVC% and DLCO% were prognostic factors for long-term survival but not for PFS [17], implying that long-term survival outcomes might not be achievable due to detriment of pulmonary function even though patients had good early response to CCRT. To further interpret the underlying role of these selected radiomic lung features, the correlation between radiomic features and PFT/BGA indicators were explored in depth and it was confirmed that FEV1/FVC% was well correlated with radiomic pulmonary features. This correlation between the pulmonary ventilation function and selected radiomic pulmonary features for OS prediction reaffirms the findings in Occhipinti et al.’s study that the changes in lung function, such as bronchial thickening and honeycombing, can be mechanistically explained based on morphological CT features [23]. And it might additionally imply that the tumor not only interacts with cells in its immediate vicinity, but also communicates with the entire host organ [30], just as suggested by a prior study [31] that the tumor and TOE could possibly interact in a bi-directional way.

In the aspect of methodology, the machine learning framework in this study used SVM combined with the proposed IFSMT approach to iteratively select features using GA and improve the accuracy of the prediction model. Our avoidance of topical deep-learning frameworks, such as deep convolutional neural network, is due to the intrinsic weaknesses of overfitting and blackbox for these frameworks. To ease the problem of overfitting, the deep-learning frameworks are more suitable for the learning tasks armed with big data as learning samples. However, the number of patients in current study for model training was relatively small, which intensively restricts the application of deep-learning frameworks which may have millions of parameters and thousands of decision making variables. The SVM is equivalent to an optimized three-layer neural network with only one hidden layer. This simplified neural network architecture substantially reduces the potential of overfitting. Additionally, in contrast with the problem of blackbox for deep learning framework, the features used in modeling are explicitly created and selected with the IFSMT approach. Therefore, each feature had an explicit clinical or physical meaning relevant to image of a specific ROI, which made it easy to apprehend the behind-the-scene mechanism of the survival status prediction and directly related the comprehensible clinical and image oriented indices to the clinical outcome. The effectiveness of IFSMT approach had been demonstrated by high AUC values achieved for the survival status prediction.

The most recent work on prognostic model for the survival outcome for NSCLC patients treated with CCRT demonstrated that pretreatment CT texture features provided prognostic information beyond CPFs [12]. However, it didn’t provide the result in terms of AUC or employ the validation cohort. In another predictive model conducted by Dehing-Oberije C et al. [32], which used CPF indices only, the AUC was 0.74 for the training cohort, 0.75 and 0.76 for the two separate validation cohorts. The improvement of model performance by imaging features in current study is discernible with the AUCs of 0.965 and 0.869 for the training and validation cohort, which could be attributed to inclusion of the image-based pulmonary features.

What’s more, the predictive OS results using imaging features in our study with machine learning could be utilized as an effective indicator for the survival risk stratification of these patients, which could potentially individualize CCRT regimen and adjuvant treatment from the perspective of personalized medicine. For example, immunotherapy has evolved into a standard adjuvant treatment option for LANSCLC patients treated with definitive CCRT. Based on the promising results of the phase III PACIFIC study [33, 34], adjuvant immunotherapy resulted in a significant prolonged PFS and OS for those patients. To be noticed, the most common grade 3 or 4 adverse event in the durvalumab arm was pneumonia (4.4%), followed by pneumonitis or radiation pneumonitis (3.4%), and Asian patients seemed to have a higher rate of any grade pneumonitis (73.6%) and severe pneumonitis (5.6%) [35]. Thus, based on the survival risk stratification of LANSCLC patients in this study, low-risk patients might have several potential advantages for adjuvant immunotherapy: (1) supporting role of better pulmonary function and quality of life; (2) superior tumor remission with less pulmonary toxicities; (3) less severe lymphopenia during CCRT and better recovery of lymphopenia from CCRT. However, for high-risk LANSCLC patients who had worse baseline FEV1/FVC%, higher rate of Grade ≥ 3 lymphopenia during CCRT, worse recovery of lymphopenia from CCRT, and higher incidence of radiation-induced pneumonitis, radical CCRT or further adjuvant immunotherapy might not be feasible because of poor organ functions and high probability of severe complications. Therefore, pretreatment radiomics-based risk stratification of LANSCLC patients using features from tumor and TOE could provide direct evidences to effectively support the treatment decision making.

It should also be noted that there were a few limitations in this study. First, the absence of external validation was the major disadvantage. Nevertheless, multiple CT simulation machines were available in our institution (Additional File 1). The high AUC values were generated from these different scanners with varied parameter settings, demonstrating the great robustness of our model. Besides, Zhao et al. considered that radiomic features in lung cancer were reproducible over a wide range of imaging settings [36]. Multicenter validations with larger samples are warranted for the ultimate application of this model clinically. Second, there might be some variability in multiple observer delineations in our study. E et al. reported that although the ROIs delineation tended to be different between individual experts, an overall high AUC value could still be achieved [37]. Third, we focused only on the radiomic analysis of pretreatment planning CT in this study, and other imaging modalities, such as PET-CT [38] and MRI, still need to be investigated as to whether they could also yield complementary information which would facilitate more accurate predictive models.

Conclusion

Pretreatment CT-based radiomics features from tumor and TOE could improve the long-term survival forecast accuracy in LANSCLC patients treated with CCRT using machine learning. The predictive results could be utilized as an effective indicator for the stratification of these patients into the low-risk and high-risk groups. It was further confirmed that patients in the low-risk group had better baseline FEV1/FVC%, less severe lymphopenia during CCRT, better recovery of lymphopenia from CCRT, lower incidence of radiation-induced pneumonitis, superior tumor remission and long-term survival, which might suggest more benefit for these patients from radical CCRT or further adjuvant immunotherapy.

Data Availability

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

Abbreviations

CI:

confidence interval

CCRT:

concurrent chemoradiotherapy

LANSCLC:

locally advanced non-small cell lung cancer

RILT:

radiation-induced lung toxicity

CPF:

clinical prognostic factor OS overall survival

TOE:

tumor organismal environment

PFT:

pulmonary function test

PFS:

progression-free survival

ORR:

objective response rate

CT:

computed tomography IFSMT integrated feature selection and model training

MRI:

magnetic resonance imaging

PET-CT:

positron emission tomography-computed tomography

4D:

four-dimensional

GTV:

gross tumor volume

ITV:

internal target volume

PTV:

planning target volume

OAR:

organ at risk

BGA:

blood gas analysis

ROI:

regions of interest

GLCM:

grey-level co-occurrence matrix

NGLDM:

neighborhood grey-level different matrix

GLRLM:

grey-level run length matrix

GLZLM:

grey-level zone length matrix

SVM:

support vector machine

GA:

genetic algorithm

LOOCV:

leave-one-out cross-validation

ROC:

receiver operating characteristics

AUC:

area under curve

IW:

importance weight TPR true positive rate

TNR:

true negative rate

References

  1. Zhou Z, Song X, Wu A, Liu H, Wu H, Wu Q, et al. Pulmonary emphysema is a risk factor for radiation pneumonitis in NSCLC patients with squamous cell carcinoma after thoracic radiation therapy. Sci Rep. 2017;7(1):2748.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Li F, Zhou Z, Wu A, Cai Y, Wu H, Chen M, et al. Preexisting radiological interstitial lung abnormalities are a risk factor for severe radiation pneumonitis in patients with small-cell lung cancer after thoracic radiation therapy. Radiat Oncol. 2018;13(1):82.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Roach M 3rd, Gandara DR, Yuo HS, Swift PS, Kroll S, Shrieve DC, et al. Radiation pneumonitis following combined modality therapy for lung cancer: analysis of prognostic factors. J Clin Oncol. 1995;13(10):2606–12.

    Article  Google Scholar 

  4. Bradley JD, Hope A, El Naqa I, Apte A, Lindsay PE, Bosch W, et al. A nomogram to predict radiation pneumonitis, derived from a combined analysis of RTOG 9311 and institutional data. Int J Radiat Oncol Biol Phys. 2007;69(4):985–92.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Hope AJ, Lindsay PE, El Naqa I, Alaly JR, Vicic M, Bradley JD, et al. Modeling radiation pneumonitis risk with clinical, dosimetric, and spatial parameters. Int J Radiat Oncol Biol Phys. 2006;65(1):112–24.

    Article  PubMed  Google Scholar 

  6. Valdes G, Solberg TD, Heskel M, Ungar L, Simone CB 2. Using machine learning to predict radiation pneumonitis in patients with stage I non-small cell lung cancer treated with stereotactic body radiation therapy. Phys Med Biol. 2016;61(16):6105–20. nd. .

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Yakar M, Etiz D, Metintas M, Ak G, Celik O. Prediction of Radiation Pneumonitis With Machine Learning in Stage III Lung Cancer: A Pilot Study. Technol Cancer Res Treat. 2021;20:15330338211016373.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Bourbonne V, Da-Ano R, Jaouen V, Lucia F, Dissaux G, Bert J, et al. Radiomics analysis of 3D dose distributions to predict toxicity of radiotherapy for lung cancer. Radiother Oncol. 2021;155:144–50.

    Article  CAS  PubMed  Google Scholar 

  9. Puttanawarut C, Sirirutbunkajorn N, Tawong N, Jiarpinitnun C, Khachonkham S, Pattaranutaporn P, et al. Radiomic and Dosiomic Features for the Prediction of Radiation Pneumonitis Across Esophageal Cancer and Lung Cancer. Front Oncol. 2022;12:768152.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Lee G, Lee HY, Park H, Schiebler ML, van Beek EJR, Ohno Y, et al. Radiomics and its emerging role in lung cancer research, imaging biomarkers and clinical management: State of the art. Eur J Radiol. 2017;86:297–307.

    Article  PubMed  Google Scholar 

  11. Thawani R, McLane M, Beig N, Ghose S, Prasanna P, Velcheti V, et al. Radiomics and radiogenomics in lung cancer: A review for the clinician. Lung cancer (Amsterdam Netherlands). 2018;115:34–41.

    Article  Google Scholar 

  12. Fried DV, Tucker SL, Zhou S, Liao Z, Mawlawi O, Ibbott G, et al. Prognostic value and reproducibility of pretreatment CT texture features in stage III non-small cell lung cancer. Int J Radiat Oncol Biol Phys. 2014;90(4):834–42.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Sun W, Jiang M, Dang J, Chang P, Yin FF. Effect of machine learning methods on predicting NSCLC overall survival time based on Radiomics analysis. Radiation Oncol (London England). 2018;13(1):197.

    Article  Google Scholar 

  14. Lu MT, Ivanov A, Mayrhofer T, Hosny A, Aerts H, Hoffmann U. Deep Learning to Assess Long-term Mortality From Chest Radiographs. JAMA Netw open. 2019;2(7):e197416.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Duijm M, van der Voort van Zyp NC, Granton PV, van de Vaart P, Mast ME, Oomen-de Hoop E, et al. Prognostic factors of local control and disease free survival in centrally located non-small cell lung cancer treated with stereotactic body radiation therapy. Acta Oncol. 2020;59(7):809–17.

    Article  PubMed  Google Scholar 

  16. Kang J, Ning MS, Feng H, Li H, Bahig H, Brooks ED, et al. Predicting 5-Year Progression and Survival Outcomes for Early Stage Non-small Cell Lung Cancer Treated with Stereotactic Ablative Radiation Therapy: Development and Validation of Robust Prognostic Nomograms. Int J Radiat Oncol Biol Phys. 2020;106(1):90–9.

    Article  CAS  PubMed  Google Scholar 

  17. Qiu B, Xiong M, Luo Y, Li Q, Chen N, Chen L, et al. Hypofractionated Intensity Modulated Radiation Therapy With Concurrent Chemotherapy in Locally Advanced Non-Small Cell Lung Cancer: A Phase II Prospective Clinical Trial (GASTO1011). Pract Radiat Oncol. 2021;11(5):374–83.

    Article  PubMed  Google Scholar 

  18. Kong FM, Ritter T, Quint DJ, Senan S, Gaspar LE, Komaki RU, et al. Consideration of dose limits for organs at risk of thoracic radiotherapy: atlas for lung, proximal bronchial tree, esophagus, spinal cord, ribs, and brachial plexus. Int J Radiat Oncol Biol Phys. 2011;81(5):1442–57.

    Article  PubMed  Google Scholar 

  19. Nioche C, Orlhac F, Boughdad S, Reuzé S, Goya-Outi J, Robert C, et al. LIFEx: A Freeware for Radiomic Feature Calculation in Multimodality Imaging to Accelerate Advances in the Characterization of Tumor Heterogeneity. Cancer Res. 2018;78(16):4786–9.

    Article  CAS  PubMed  Google Scholar 

  20. Smola AJ, Scholkopf B. A tutorial on support vector regression. Stat Comput. 2004;14(3):199–222.

    Article  Google Scholar 

  21. Wu MT, Chang JM, Chiang AA, Lu JY, Hsu HK, Hsu WH, et al. Use of quantitative CT to predict postoperative lung function in patients with lung cancer. Radiology. 1994;191(1):257–62.

    Article  CAS  PubMed  Google Scholar 

  22. Lafata KJ, Zhou Z, Liu JG, Hong J, Kelsey CR, Yin FF. An Exploratory Radiomics Approach to Quantifying Pulmonary Function in CT Images. Sci Rep. 2019;9(1):11509.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Occhipinti M, Paoletti M, Bartholmai BJ, Rajagopalan S, Karwoski RA, Nardi C, et al. Spirometric assessment of emphysema presence and severity as measured by quantitative CT and CT-based radiomics in COPD. Respir Res. 2019;20(1):101.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Ganeshan B, Panayiotou E, Burnand K, Dizdarevic S, Miles K. Tumour heterogeneity in non-small cell lung carcinoma assessed by CT texture analysis: a potential marker of survival. Eur Radiol. 2012;22(4):796–802.

    Article  PubMed  Google Scholar 

  25. Bradley JD, Ieumwananonthachai N, Purdy JA, Wasserman TH, Lockett MA, Graham MV, et al. Gross tumor volume, critical prognostic factor in patients treated with three-dimensional conformal radiation therapy for non-small-cell lung carcinoma. Int J Radiat Oncol Biol Phys. 2002;52(1):49–57.

    Article  PubMed  Google Scholar 

  26. Lopez Guerra JL, Gomez D, Zhuang Y, Levy LB, Eapen G, Liu H, et al. Change in diffusing capacity after radiation as an objective measure for grading radiation pneumonitis in patients treated for non-small-cell lung cancer. Int J Radiat Oncol Biol Phys. 2012;83(5):1573–9.

    Article  PubMed  Google Scholar 

  27. Park YH, Kim JS. Predictors of radiation pneumonitis and pulmonary function changes after concurrent chemoradiotherapy of non-small cell lung cancer. Radiation Oncol J. 2013;31(1):34–40.

    Article  Google Scholar 

  28. Torre-Bouscoulet L, Muñoz-Montaño WR, Martínez-Briseño D, Lozano-Ruiz FJ, Fernández-Plata R, Beck-Magaña JA, et al. Abnormal pulmonary function tests predict the development of radiation-induced pneumonitis in advanced non-small cell lung Cancer. Respir Res. 2018;19(1):72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Zhou Y, Yan T, Zhou X, Cao P, Luo C, Zhou L, et al. Acute severe radiation pneumonitis among non-small cell lung cancer (NSCLC) patients with moderate pulmonary dysfunction receiving definitive concurrent chemoradiotherapy: Impact of pre-treatment pulmonary function parameters. Strahlentherapie und Onkologie: Organ der Deutschen Rontgengesellschaft [et al]. 2020;196(6):505–14.

    Article  Google Scholar 

  30. Dieterich LC, Bikfalvi A. The tumor organismal environment: Role in tumor development and cancer immunotherapy. Sem Cancer Biol. 2020;65:197–206.

    Article  Google Scholar 

  31. Khorrami M, Prasanna P, Gupta A, Patil P, Velu PD, Thawani R, et al. Changes in CT Radiomic Features Associated with Lymphocyte Distribution Predict Overall Survival and Response to Immunotherapy in Non-Small Cell Lung Cancer. Cancer Immunol Res. 2020;8(1):108–19.

    Article  CAS  PubMed  Google Scholar 

  32. Dehing-Oberije C, Yu S, De Ruysscher D, Meersschout S, Van Beek K, Lievens Y, et al. Development and external validation of prognostic model for 2-year survival of non-small-cell lung cancer patients treated with chemoradiotherapy. Int J Radiat Oncol Biol Phys. 2009;74(2):355–62.

    Article  PubMed  Google Scholar 

  33. Antonia SJ, Villegas A, Daniel D, Vicente D, Murakami S, Hui R, et al. Durvalumab after Chemoradiotherapy in Stage III Non-Small-Cell Lung Cancer. N Engl J Med. 2017;377(20):1919–29.

    Article  CAS  PubMed  Google Scholar 

  34. Antonia SJ, Villegas A, Daniel D, Vicente D, Murakami S, Hui R, et al. Overall Survival with Durvalumab after Chemoradiotherapy in Stage III NSCLC. N Engl J Med. 2018;379(24):2342–50.

    Article  CAS  PubMed  Google Scholar 

  35. Peng L, Wu YL. Immunotherapy in the Asiatic population: any differences from Caucasian population? J Thorac disease. 2018;10(Suppl 13):1482-s93.

    Google Scholar 

  36. Zhao B, Tan Y, Tsai WY, Schwartz LH, Lu L. Exploring Variability in CT Characterization of Tumors: A Preliminary Phantom Study. Translational Oncol. 2014;7(1):88–93.

    Article  Google Scholar 

  37. Lu EL, Li L, Yang L, Schwartz H, Zhao LH. B. Radiomics for Classification of Lung Cancer Histological Subtypes Based on Nonenhanced Computed Tomography. Acad Radiol. 2019;26(9):1245–52.

    Article  PubMed  Google Scholar 

  38. Cook GJ, Yip C, Siddique M, Goh V, Chicklore S, Roy A, et al. Are pretreatment 18F-FDG PET tumor textural features in non-small cell lung cancer associated with response and survival after chemoradiotherapy? Journal of nuclear medicine: official publication. Soc Nuclear Med. 2013;54(1):19–26.

    Article  Google Scholar 

Download references

Funding

This work was supported by National Natural Science Foundation of China [Grant Number 82073328], National Key R&D Program of China [Grant Number 2018YFC0116800] and Suzhou Municipal Science and Technology Program [Grant Number SYS2018009]. The funders had no role in study design; in the collection, analysis and interpretation of data; in the writing of the report; and in the decision to submit the article for publication.

Author information

Author notes

  1. Nai-Bin Chen, Mai Xiong and Rui Zhou contributed equally as co-first authors.

Authors and Affiliations

  1. Department of Radiation Oncology, State Key Laboratory of Oncology in South China, Collaborative Innovation Center for Cancer Medicine, Sun Yat-sen University Cancer Center, No.651 Dongfeng Road East, 510060, Guangzhou, China

    Nai-Bin Chen, Rui Zhou, Bo Qiu, Chu Chu, Qi-Wen Li, Bin Wang, Jin-Yu Guo & Hui Liu

  2. Department of Imaging Diagnosis and Interventional Center, State Key Laboratory of Oncology in South China, Collaborative Innovation Center for Cancer Medicine, Sun Yat-sen University Cancer Center, No.651 Dongfeng Road East, 510060, Guangzhou, Guangdong, China

    Kang-Qiang Peng & Chuan-Miao Xie

  3. Department of Cardiac Surgery, The First Affiliated Hospital of Sun Yat-sen University, Guangzhou, Guangdong, China

    Mai Xiong

  4. Department of Pulmonary and Critical Care Medicine, The First Affiliated Hospital of Sun Yat-sen University, Guangzhou, Guangdong, China

    Yi-Feng Luo

  5. Homology Medical Technologies Inc., Ningbo, Zhejiang, China

    Yin Zhou & Hai-Hang Jiang

  6. Guangzhou Xinhua University, Guangzhou, Guangdong, China

    Su Zhou

Authors
  1. Nai-Bin Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Mai Xiong
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Rui Zhou
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Yin Zhou
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Bo Qiu
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Yi-Feng Luo
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Su Zhou
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Chu Chu
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Qi-Wen Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Bin Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Hai-Hang Jiang
    View author publications

    You can also search for this author in PubMed Google Scholar

  12. Jin-Yu Guo
    View author publications

    You can also search for this author in PubMed Google Scholar

  13. Kang-Qiang Peng
    View author publications

    You can also search for this author in PubMed Google Scholar

  14. Chuan-Miao Xie
    View author publications

    You can also search for this author in PubMed Google Scholar

  15. Hui Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

HL and CMX contributed to the conception and design of the study. NBC, MX, RZ, BQ, YFL, SZ, CC, QWL, BW, JYG and KQP contributed to acquisition of data. NBC, MX, RZ, YZ, BQ, YFL and HHJ analyzed and interpreted data. NBC, MX and RZ drafted the article. HL and CMX revised it critically and made the final approval of the version to be submitted. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Chuan-Miao Xie or Hui Liu.

Ethics declarations

Ethics approval and consent to participate

This study was conducted based on the ethical standards of the Declaration of Helsinki. It was reviewed and approved by the review board of Sun Yat-sen University Cancer Center. Since it was an anonymous retrospective study, we required and were granted a waiver of written informed-consent.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chen, NB., Xiong, M., Zhou, R. et al. CT radiomics-based long-term survival prediction for locally advanced non-small cell lung cancer patients treated with concurrent chemoradiotherapy using features from tumor and tumor organismal environment. Radiat Oncol 17, 184 (2022). https://doi.org/10.1186/s13014-022-02136-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13014-022-02136-w

Keywords

Radiation Oncology

ISSN: 1748-717X

Contact us