ORIGINAL RESEARCH article

Front. Oncol., 03 May 2023

Sec. Thoracic Oncology

Volume 13 - 2023 | https://doi.org/10.3389/fonc.2023.1131816

Preoperative computed tomography-based tumoral radiomic features prediction for overall survival in resectable non-small cell lung cancer

  • 1. Department of Thoracic Surgery, The Second Affiliated Hospital of Harbin Medical University, Harbin, China

  • 2. Department of Radiology, The Second Affiliated Hospital of Harbin Medical University, Harbin, China

Abstract

Objectives:

The purpose of this study was to evaluate whether preoperative radiomics features could meliorate risk stratification for the overall survival (OS) of non-small cell lung cancer (NSCLC) patients.

Methods:

After rigorous screening, the 208 NSCLC patients without any pre-operative adjuvant therapy were eventually enrolled. We segmented the 3D volume of interest (VOI) based on malignant lesion of computed tomography (CT) imaging and extracted 1542 radiomics features. Interclass correlation coefficients (ICC) and LASSO Cox regression analysis were utilized to perform feature selection and radiomics model building. In the model evaluation phase, we carried out stratified analysis, receiver operating characteristic (ROC) curve, concordance index (C-index), and decision curve analysis (DCA). In addition, integrating the clinicopathological trait and radiomics score, we developed a nomogram to predict the OS at 1 year, 2 years, and 3 years, respectively.

Results:

Six radiomics features, including gradient_glcm_InverseVariance, logarithm_firstorder_Median, logarithm_firstorder_RobustMeanAbsoluteDeviation, square_gldm_LargeDependenceEmphasis, wavelet_HLL_firstorder_Kurtosis, and wavelet_LLL_firstorder_Maximum, were selected to construct the radiomics signature, whose areas under the curve (AUCs) for 3-year prediction reached 0.857 in the training set (n=146) and 0.871 in the testing set (n=62). The results of multivariate analysis revealed that the radiomics score, radiological sign, and N stage were independent prognostic factors in NSCLC. Moreover, compared with clinical factors and the separate radiomics model, the established nomogram exhibited a better performance in predicting 3-year OS.

Conclusions:

Our radiomics model may provide a promising non-invasive approach for preoperative risk stratification and personalized postoperative surveillance for resectable NSCLC patients.

Introduction

Non-small cell lung cancer (NSCLC), the most essential subtype of lung cancer, represents a prevalent malignant tumor with an unsatisfactory prognosis (1). In recent years, an expanding availabilities of targeted therapy and immune checkpoint inhibitors (ICIs), has been approved for lung cancer and improved the long-term prognosis of NSCLC patients. However, due to the secondary mutations and low response rate to ICIs, only a limited number of patients can benefit from those therapeutic approaches (2, 3). Thus, precision diagnose and comprehensive prognostic evaluation are essential steps when dealing with resectable NSCLC patients, in order to select the most appropriate treatment. The tumor node metastasis (TNM) staging system is still a classic evaluation approach which can assist in adjuvant therapy choices and predict the outcome of NSCLC (4). Nonetheless, patients with the same TNM stage typically manifest different clinical outcomes, which is largely attributed to tumor heterogeneity and anatomical factors (5). With the advent of multi-omics, evaluation of prognostic features based on multidisciplinary methods make personalized medicine possible.

Medical imaging has been long regarded as standard procedure for early screening, treatment decision-making, and postoperative surveillance of cancer patients. Computed tomography (CT) imaging, commonly stored in the form of Digital Imaging and Communications in Medicine (DICOM), can be conveniently obtained and utilized for quantitative assessment. Over the last decade, radiomics has emerged as a hot research field that provides massive high-dimensional feature space derived from raw imaging data by automatically high-throughput algorithm (6). There is growing evidence that quantitative parameters and features mined from functional and morphological images offer a new perspective for tumor phenotypes and microenvironment, which also have a significant complementary interrelation with other omics approaches such as genomics, hematology, and proteomics (79). Radiomics analyses based on intratumoral and peritumoral regions have been extensively used for exploring underlying biological process, predicting pathological characteristics, evaluating the drug treatment response, and assisting therapeutic decision-making in several human carcinomas (1014).

As a novelly emerging tool, radiomics provides a new direction of exploring intratumor heterogeneity and predictive markers using a noninvasive evaluation in lung cancer. Using three machine-learning (ML) classifiers derived from radiomics, Liu et al. estimated the benefit from Nivolumab treatment and progression probability in patients with stage IIIB/IV NSCLC (15). In terms of histological subtypes, an automatic deep-learning radiomics model showed satisfactory performance to distinguish lung squamous cell carcinoma (LUSC), lung adenocarcinoma (LUAD), and small cell lung cancer (SCLC) on CT images (16). However, there are limited studies investigating the contingent value of radiomics in improving prognosis stratification of resectable NSCLC patients.

Accordingly, the current study aimed to analyze the association of radiomics features with 3-year overall survival (OS) in enrolled patient cohort. Furthermore, the integration of radiomics model and clinicopathological traits was conducted to establish a comprehensive nomogram, which strengthen the prediction ability and may provide assistance in improving follow-up plans and individualized treatment in clinical practice.

Materials and methods

Patient data and study design

This study was approved by the Ethics and Scientific Committees of the Second Affiliated Hospital of Harbin Medical University (Approval Number: KY2022-144), due to the retrospective nature of this study, written informed consent for CT images was waived. The study design is illustrated in Figure 1.

Figure 1

493 consecutive NSCLC patients who underwent radical surgery (segmental resection, wedge resection, and lobectomy) at our institute from March 2015 to May 2019 were preliminarily included. All patients fulfilled the following inclusion criteria: (1) CT performed within 2 weeks before surgery; (2) Available CT images stored in DICOM format; and (3) Primary NSCLC confirmed by histopathology. The subsequent patient selection and exclusion criteria were visualized in Figure 1A. Consequently, we finally recruited 208 consecutive patients, including 167 cases with LUAD, 35 cases with LUSC, 3 cases with large cell carcinoma (LCC), 2 cases with adenosquamous carcinoma (ASC), and 1 case with carcinoid. According to the random allocation scheme and a ratio of 7:3, all patients were separated into two individual cohorts: 146 for training and 62 for testing. Patient grouping and corresponding clinicopathological traits are recapitulated in Table 1.

Table 1

CovariatesTraining set
(n=146)
Testing set
(n=62)
Total set
(n=208)
P-value
Age, no (%)≤ 6599(67.81)48(77.42)147(70.67)0.220
> 6547(32.19)14(22.58)61(29.33)
Gender, no (%)Female77(52.74)35(56.45)112(53.85)0.735
Male69(47.26)27(43.55)96(46.15)
Smoking history, no (%)No89(60.96)41(66.13)130(62.50)0.584
Yes57(39.04)21(33.87)78(37.50)
T stage, no (%)T1-2131(89.73)56(90.32)187(89.90)1.000
T3-415(10.27)6(9.68)21(10.10)
N stage, no (%)N0115(78.77)48(77.42)163(78.37)0.975
N1-231(21.23)14(22.58)45(21.63)
Pathologic TNM stage, no (%)Stage I101(69.18)46(74.19)147(70.67)0.754
Stage II21(14.38)7(11.29)28(13.46)
Stage III24(16.44)9(14.52)33(15.87)
Lateral location, no (%)Left60(41.10)30(48.39)90(43.27)0.413
Right86(58.90)32(51.61)118(56.73)
Lobe location, no (%)Mid-lower65(44.52)32(51.61)97(46.63)0.432
Upper81(55.48)30(48.39)111(53.37)
Location classification, no (%)Central29(19.86)7(11.29)36(17.31)0.196
Peripheral117(80.14)55(88.71)172(82.69)
Max diameter, no (%)≤3cm112(76.71)49(79.03)161(77.40)0.854
>3cm34(23.29)13(20.97)47(22.60)
Pathological type, no (%)LUAD112(76.71)55(88.71)167(80.29)0.138
LUSC29(19.86)6(9.68)35(16.83)
Other5(3.42)1(1.61)6(2.88)
Histological grade, no (%)Well24(16.44)8(12.9)32(15.38)0.644
Moderate44(30.14)15(24.19)59(28.37)
Poor50(34.25)26(41.94)76(36.54)
Unknown28(19.18)13(20.97)41(19.71)
Radiological sign, no (%)Cavitary appearance30(20.55)21(33.87)51(24.52)0.091
Part-solid appearance40(27.4)12(19.35)52(25)
Pure ground-glass appearance16(10.96)10(16.13)26(12.5)
Pure solid appearance60(41.1)19(30.65)79(37.98)

Patient grouping and corresponding clinicopathological traits.

LUAD, lung adenocarcinoma; LUSC, lung squamous cell carcinoma.

Follow-up

After radical surgery, all enrolled subjects were followed up by outpatient review or telephone every 3 months for the first year and every 6 months thereafter. We applied 3-year OS as the primary study endpoint, which is construed as the time between the operation and the date of all-cause death.

CT image acquisition and pre-processing

Helical CT images of all enrolled subjects were acquired by 64-channel CT (Discovery 750, GE Healthcare, Milwaukee, USA) and 256-channel CT (Revolution CT, GE Healthcare, Waukesha, WI, USA). Detailed scanning parameters were as follows: tube voltage, 120 kV; tube current, 100-250 mAs; slice thickness, 0.625-5 mm; field of view (FOV), 350-400mm; 512 x 512 matrix; and reconstructed slice thickness, 0.625-3mm. Filtered back projection (FBP) and adaptive statistical iterative reconstruction (ASIR) level 40% were utilized to reconstruct all enrolled CT images. A standard kernel was also used in the reconstruction procedure.

Due to the potential differences in specifications caused by distinct reconstruction slice thickness and voxel spacing, we performed image pre-processing before radiomics feature extracting. Specifically, we chose to resample the raw data to a standard voxel spacing of 1x1x1 mm3 by near interpolation algorithm (17). Moreover, we adopted 25 HU as the fixed bin width to perform gray level discretization (18).

Tumor segmentation and radiomics feature extraction

CT images stored in DICOM format were loaded to 3D slicer (Version 4.1.1) software to perform tumor segmentation. As an open source platform, 3D slicer software has superiority in image interactive segmentation and medical raw data processing (19). A thoracic surgeon with 5 years of experience first manually delineated all lesions slice by slice in the lung window of CT images, the volumes of interest (VOIs) were then inspected and revised by an experienced radiologist. They were blinded to the clinical outcome and medical records of enrolled patients, in parallel, intratumoral or tumor-adjacent vessels, bronchi, and air spaces were carefully avoided during delineation of the tumor volume. To improve the reproducibility of radiomics features, the VOIs of 50 randomly selected CT images were repetitively delineated by another experienced thoracic surgeon, which were used for subsequent interclass correlation coefficient (ICC) analysis.

In our study, 1542 radiomics features were extracted from each manually-defined VOI using “pyradiomics” python package, the detailed definitions of which are congruity with Imaging Biomarker Standardization Initiative (IBSI). The extracted radiomics features can be assigned to the following categories: (1) Gray level dependence matrix-based features (GLDM); (2) neighboring gray tone difference matrix-based features (NGTDM); (3) gray level co-occurrence matrix-based features (GLCM); (4) gray level run-length matrix-based features (GLRLM); (5) gray level size zone matrix-based features (GLSZM); (6) first-order statistics features; (7) shape-based features, and (8) transformed features: features extracted from images pre-processed with several built-in filters including laplacian of gaussian (LoG), wavelet, logarithm, square, square root, and gradient. Detailed mathematic definitions and feature explanations can be acquired from the previous literature (20, 21).

Feature selection and radiomics model construction

To assess the consistency and robustness of extracted radiomics features between the two thoracic surgeons, we calculated the interclass correlation coefficient (ICC) value for each radiomics feature. The radiomics features with ICC values > 0.75 were retained for deeper analysis. Next, we performed data pre-processing using Z-score transformation by “caret” R package. Least absolute shrinkage and selection operator (LASSO) Cox regression analysis were utilized for further feature selection and model construction. In the training set, LASSO with ten-fold cross-validation effectively avoided over-fitting and identified optimal features with nonzero coefficients by applying a constraint on the model hyperparameter (λ). Finally, we obtained a linear combination of optimal radiomics features weighted with the regression coefficients during multiple computing. All the selected features for model construction are detailed in Table 2.

Table 2

Selected radiomics featureRadiomics groupFilter associatedRegression coefficient
R1InverseVarianceGLCMgradient0.068
R2MedianFirst orderlogarithm0.096
R3RobustMeanAbsoluteDeviationFirst orderlogarithm0.181
R4LargeDependenceEmphasisGLDMsquare0.375
R5KurtosisFirst orderwavelet_HLL0.016
R6MaximumFirst orderwavelet_LLL0.269

Six selected radiomics features.

GLCM, gray level co-occurrence matrix; GLDM, Gray level dependence matrix.

Radiomics model validation and evaluation

The prognostic efficacy of radiomics model established in the training set was subsequently verified in the testing set and the total enrolled subjects. Patients in these three sets were respectively dichotomized into low- and high-risk groups based on the median radiomics score threshold. To appraise diversities in the overall survival between two cohorts with different radiomics score, R packages survival and survMiner were used to perform Kaplan-Meier survival analysis. We utilized receiver operating characteristic (ROC) curve analysis and calculated the area under the curve (AUC) to access the sensitivity and specificity of the prognostic model. The R packages survminer, survival, and timeROC were employed to enable this process. Moreover, Harrell’s concordance index (C-index) was applied to estimate the predictive accuracy of the signature. AUC and C-index both ranges from 0.5 (poor predictive performance) to 1 (perfect predictive performance).

A nomogram comprising of radiomics score and clinicopathological parameters (age, gender, smoking history, pathological type, lobe location, location classification, lateral location, max diameter, T stage, N stage, and TNM stage) was set up to predict the 1-, 2-, and 3-year OS of resectable NSCLC patients. Correction curve analysis, ROC curve analysis, and C-index were applied to evaluate the predictive performance of the nomogram. In addition, we performed decision curve analysis (DCA) to compare the clinical benefit of the classical TNM stage, established radiomics model, and the nomogram by quantifying the net benefits (22).

Statistical analyses

All statistical analyses were carried out by R software (version 4.1.2) and Python software (version 3.4.3). The statistically significant threshold was set to p value < 0.05 (*p < 0.05, **p < 0.01, ***p < 0.001).

Results

Construction of the radiomics model

All enrolled 208 subjects were randomly divided into training set (n=146) and testing set (n=62) with a ratio of 7:3, and there was no statistical difference in terms of clinical features (Table 1). The training set was applied to create the prognostic radiomics model, we then used the testing set and the total set to validate the constructed model.

According to the 1542 radiomics features extracted by two experienced thoracic surgeons, The ICC analysis indicated good consistency in 1178 (76.39%) radiomics features, moderate consistency in 266 (17.25%) radiomics features, and poor consistency in 98 (6.36%) radiomics features. We ultimately included the 1178 most stable features in the subsequent model construction. LASSO Cox regression analysis is widely used in multiple regression analysis, which not only optimize the selection of characteristics with a deficient correlation and prominent predicted value from high-dimensional data, but also improve the forecast accuracy. We performed tenfold cross-validation to select the minimal penalty term (λ) (Figures 2A, B).

Figure 2

Afterwards, we established a prognostic radiomics model implicating six selected features for resectable NSCLC patients. The formula was constructed as follows: Radiomics score = 0.068 × gradient_glcm_InverseVariance + 0.096 × logarithm_firstorder_Median + 0.181 × logarithm_firstorder_RobustMeanAbsoluteDeviation + 0.375 × square_gldm_LargeDependenceEmphasis + 0.016 × wavelet_HLL_firstorder_Kurtosis + 0.269 × wavelet_LLL_firstorder_Maximum.

Survival analyses in the training set, testing set, and the entire NSCLC set

We first dichotomize 146 NSCLC patients of training set into low- and high-risk groups based on the median radiomics score. Kaplan-Meier analysis showed significantly distinct prognoses between the two risk groups (Figure 2C). In addition, the distributions of the radiomics score, survival time, and survival status of each NSCLC patient in the training set were shown in Figure 2D. To further verify the prognostic performance of the constructed radiomics model, we acquired the radiomics score of each NSCLC patient in the testing set and the total NSCLC set using the radiomics score formula. Employing similar segmentation method, we divided the testing set and the total NSCLC set into subgroups with different risk level. The Kaplan-Meier analysis demonstrated that the patients in the high-risk group had significantly poorer outcome (Figures 2F, 3A). The distribution patterns of the radiomics score, survival time, and survival status of both sets were illustrated in Figures 2G, 3B.

Figure 3

Performance evaluation of the radiomics model

ROC curve analysis indicated that the AUC of established radiomics model reached 0.857 and 0.871 at 3 years in the training set and the testing set (Figures 2E, H). Notably, the AUC of our model in total NSCLC set reached 0.792, 0.839, and 0.861 at 1 year, 2 years, and 3 years, respectively (Figure 3C). When compared with other clinicopathological variables in terms of ROC analysis, the radiomics model still had an advantage in evaluating precision and sensitivity (Figure 3D). Furthermore, The C-index curve showed promising predictive accuracy of the radiomics model from another dimension (Figure 3E).

To assess the independent prognostic value of the radiomics model, we incorporated age, gender, smoking history, pathological type, lobe location, location classification, lateral location, max diameter, T stage, N stage, TNM stage, radiological sign, and radiomics score into the univariate and multivariate Cox regression analyses. The results suggested that the established radiomics model, N stage, and radiological sign could act as independent predictors for overall survival in resectable NSCLC patients (Figure 4, p < 0.001, p < 0.05, and p < 0.05).

Figure 4

To examine the applicability of our radiomics model in subgroups stratified by different clinicopathological traits, the total 208 NSCLC subjects were dichotomized into disparate subgroups by age, gender, smoking history, histological grade, lobe location, location classification, lateral location, max diameter, T stage, N stage, TNM stage, and radiological sign. In the overwhelming majority of subgroups, our radiomics model could accurately distinguish those high-risk patients with worse outcome (Figure 5 and Supplementary Figure S1).

Figure 5

Correlation between constructed radiomics model and clinicopathological features

We next performed correlation analysis to dig deeper connections between radiomics score and the clinicopathological traits obtained from the electronic medical record. The results revealed that higher radiomics score with worse survival was significantly associated with age > 65 (p < 0.05), male (p < 0.001), smoking history (p < 0.001), T2-4 (p < 0.001), N1-2 (p < 0.001), TNM stage II-III (p < 0.001), max diameter > 3cm (p < 0.001), poor differentiation grade (p < 0.001), central-type NSCLC (p < 0.001), LUSC and other pathological types (p < 0.01), middle and lower lobes (p < 0.01), and pure solid appearance (p < 0.001), which is similar to previous experiences acquired from clinical practice and may provide additional clues to clinical management of NSCLC (Figure 6A and Supplementary Figure S2).

Figure 6

Development and assessment of the prognostic nomogram

In order to furnish a comprehensive prognostic tool to predict the survivability of resectable NSCLC patients at 1, 2, and 3 years, we constructed a nomogram integrating the radiomics model and clinicopathological parameters (Figure 6B). In terms of predicting OS at 1 year, 2 years, and 3 years, calibration plot displayed that there is a decent consistency between the prediction curve and the ideal curve (Figure 7A). The AUC values of the prognostic nomogram were 0.863, 0.870, and 0.898 at 1-, 2−, and 3−year, respectively, which exhibited better performance than radiomics model and other clinicopathological parameters (Figures 7B, C). Moreover, the C-index values of the comprehensive nomogram and established radiomics model reached 0.854 and 0.814 at 3 years, respectively. Clinical utility is commonly utilized to measure the practical clinical value of prognostic models. Subsequent DCA curves revealed that the appropriate combination of radiomics model and the comprehensive nomogram may bring significantly more benefit than TNM staging system in clinical work (Figure 7D).

Figure 7

Discussion

In this retrospective study, we explored the potential value of radiomics features in predicting 3-year OS for resectable NSCLC patients varied from pathologic stage IA to stage IIIA. The established radiomics model exhibited good prediction accuracy with an AUC of 0.857 in the training set (n = 146) and an AUC of 0.871 in the testing set (n = 62). The stratified analysis indicated that our radiomics model could distinguish those patients with worse prognosis in the vast majority of subgroups. Subsequently, a comprehensive nomogram incorporating radiomics model and clinical parameters further enhanced the prognostic performance of the single radiomics model with an AUC of 0.898. Thus, this radiomics classifier could be an advantageous noninvasive biomarker in the whole-course clinical management of resectable NSCLC patients.

The generic workflow of NSCLC biomarkers comprises various prognostic models reliant on clinical elements including lserum tumor markers, specific gene expression, next-generation sequencing, and circulating tumor DNA (ctDNA) (2325). For instance, lu et al. constructed the tumor mutation index (TMI) model based on ctDNA sequencing to predict OS and recognize NSCLC patients who may respond well to monotherapy with atezolizumab or docetaxel (26). To improve NSCLC treatment in terms of chemotherapy and immunotherapy, Guo et al. developed a 7-gene predictive signature using qRT-PCR assays based on 337 snap-frozen NSCLC tissues. Patients receiving adjuvant chemotherapy were precisely identified with significantly better disease-specific survival in the predicted benefit group via the multi-gene prognostic signature (27). However, the few available predictive models based on quantitative gene expression levels in the clinical practice have usual limitations such as invasive procedures, time-consuming, cost-effectiveness, and some degree of interference to clinical workflow.

In contrast, radiomics have exhibited a bright prospect for prognosis, diagnosis, and treatment response prediction, as well as long-term health surveillance of NSCLC treatment in a non-invasive modality. Several studies focused on exploring potential combined signature based on radiomics features for NSCLC patients at specific pathologic stage. Huang et al. extracted 132 texture features from CT images of early stage NSCLC (stage I or II) and obtained a better performance for disease-free survival (DFS) prediction (C-index = 0.72) when incorporating the radiomics signature into a comprehensive nomogram (28). Xie et al. enrolled 554 candidates with resected stage I LUAD from three multicenter cohorts and further recognized potential subjects who may benefit from adjuvant chemotherapy (29). Moreover, a number of studies mined high-dimensional clues from functional and metabolic images of 18F-FDG-PET/CT, aiming to improve clinical decision of epidermal growth factor receptor (EGFR) tyrosine kinase inhibitors (TKIs) or immune checkpoint inhibitors (ICIs) treatments for specific NSCLC patient populations (3034). In this study, we restricted the entrance criteria to preoperative CT scans, since CT scans serve as the major approaches for lung cancer screening and whole-course monitoring in a real-world clinical environment, especially in relatively early stages of NSCLC.

Six screened radiomics features used in our prediction model consisted of four first-order-based features (Maximum, Median, Kurtosis, and RobustMeanAbsoluteDeviation), one GLDM-based feature (LargeDependenceEmphasis), and one GLCM-based feature (InverseVariance), which were extracted from CT images pre-processed with built-in filters including wavelet_HLL, wavelet_LLL, logarithm, square, and gradient. First-order statistics features are generally acquired by measuring the gray values of region of interest (ROI) cropping, which reflects the intratumoral distribution of grayscale intensity (20). Previous literature indicated that two First-order statistics features may serve as radiomics predictors for identifying invasive phenotype of solitary pulmonary nodule (35). Moreover, GLDM-based features manifest the intrinsic grayscale associations of central voxel with neighboring voxels, which may be a reflection of heterogeneity and homogeneity of tumors. Padmakumari et al. demonstrated that LargeDependenceEmphasis, a radiomics feature from GLDM, can exhibit robust performance in discriminating lung cancer from tuberculosis with an AUC of 0.92 (36). GLCM describes the integrated information about spatial correlation characteristics of pixel pairs in terms of the pattern of grayscale arrangement, direction, distance, and gray value (37). Notably, textural features derived from GLCM have been demonstrated to have pathological association and can be applied to the diagnosis of malignant lesion in breast cancer (38). Granata et al. identified Correlation from GLCM as a reliable predictor for recognizing tumor recurrence in colorectal liver metastases patients (39).

In our study, we identified the established radiomics model, N stage, and radiological sign as independent predicting indicators for OS in in resectable NSCLC patients. Indeed, several studies have confirmed the accuracy and reliability of radiomics signature in predicting prognosis in NSCLC (4042). Yang et al. incorporated the radiomics signatures and four clinicopathological features (N stage, T stage, age, and sex) to construct a comprehensive nomogram for survival prediction in NSCLC patients at stage I/II, the performance of which was measured by a C-index of 0.710 (41). The radiomics model and corresponding nomogram in this study showed more satisfactory performance with C-index values of 0.814 and 0.854, respectively. Subsequent DCA and calibration analyses further supported their clinical utility. In the current clinical practice, adjuvant chemotherapy after surgery is not recommended for patients with pathologic stage IA and there is a controversial debate regarding its potential benefit for stage IB (43). In the model evaluation section, our radiomics model seemed to be more robust than traditional TNM staging system. In stratified analyses, the established radiomics model could dichotomize participants with pathologic stage I into high- and low-risk groups by Kaplan-Meier method, which would be helpful to performed personalized treatment interventions on these high-risk patients with worse prognosis.

Some limitations of this study have to be acknowledged. First, due to the relatively small sample size from single center and the retrospective nature, potential selection bias may obstruct the robustness and generalizability of our radiomics model. Therefore, it is necessary to recruit more subjects and perform multicentric external verification in future research. Second, the 3-year follow-up period of included subjects was relatively short, we will conduct the remaining follow-up until 5 years in the next work. Third, we delineated all lesions manually, which was laborious and time-consuming. Automatic delineation based on deep learning method is worth further study to improve the workflow of radiomics in busy clinical practice.

In summary, the current study proposed a novel non-invasive approach based on preoperative CT scans that can predict OS in patients with NSCLC after radical surgery, which may provide clues to help clinicians improve clinical decisions and guide personalized treatment. However, further external validation is warranted before its widespread application in clinical practice.

Statements

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving human participants were reviewed and approved by The Ethics and Scientific Committees of the Second Affiliated Hospital of Harbin Medical University. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author contributions

BP conceived the study, conducted radiomics analysis, and wrote the preliminary draft. KW and RX participated in the radiomics analysis and revised the manuscript. CG, TL, YL, CW, and XC contributed to study design and data preprocess. YW contributed to tumor segmentation and provided professional suggestions for statistical analysis. ZS, JS, and CX contributed to data collection. The manuscript was revised and approved by LZ. All authors contributed to the article and approved the submitted version.

Acknowledgments

We sincerely thank professor Huijie Jiang (Department of Radiology, The Second Affiliated Hospital of Harbin Medical University, China), for her precious suggestions to the work. We also thank all participants for their contribution to the study.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2023.1131816/full#supplementary-material

References

  • 1

    MullerDCJohanssonMBrennanP. Lung cancer risk prediction model incorporating lung function: development and validation in the UK biobank prospective cohort study. J Clin Oncol (2017) 35(8):861–9. doi: 10.1200/JCO.2016.69.2467

  • 2

    HeJHuYHuMLiB. Development of PD-1/PD-L1 pathway in tumor immune microenvironment and treatment for non-small cell lung cancer. Sci Rep (2015) 5:13110. doi: 10.1038/srep13110

  • 3

    PascheBGrantS. Non-small cell lung cancer and precision medicine: a model for the incorporation of genomic features into clinical trial design. JAMA (2014) 311(19):1975–6. doi: 10.1001/jama.2014.3742

  • 4

    GroomePABolejackVCrowleyJJKennedyCKrasnikMSobinLHet al. The IASLC lung cancer staging project: validation of the proposals for revision of the T, n, and m descriptors and consequent stage groupings in the forthcoming (seventh) edition of the TNM classification of malignant tumours. J Thorac Oncol (2007) 2(8):694705. doi: 10.1097/JTO.0b013e31812d05d5

  • 5

    Rami-PortaRAsamuraHGoldstrawP. Predicting the prognosis of lung cancer: the evolution of tumor, node and metastasis in the molecular age-challenges and opportunities. Transl Lung Cancer Res (2015) 4(4):415–23. doi: 10.3978/j.issn.2218-6751.2015.07.11

  • 6

    TunaliIGilliesRJSchabathMB. Application of radiomics and artificial intelligence for lung cancer precision medicine. Cold Spring Harb Perspect Med (2021) 11(8):a039537. doi: 10.1101/cshperspect.a039537

  • 7

    GilliesRJKinahanPEHricakH. Radiomics: images are more than pictures, they are data. Radiology (2016) 278(2):563–77. doi: 10.1148/radiol.2015151169

  • 8

    WeiWWangKLiuZTianKWangLDuJet al. Radiomic signature: a novel magnetic resonance imaging-based prognostic biomarker in patients with skull base chordoma. Radiother Oncol (2019) 141:239–46. doi: 10.1016/j.radonc.2019.10.002

  • 9

    ThawaniRMcLaneMBeigNGhoseSPrasannaPVelchetiVet al. Radiomics and radiogenomics in lung cancer: a review for the clinician. Lung Cancer (2018) 115:3441. doi: 10.1016/j.lungcan.2017.10.015

  • 10

    WesdorpNJHellingmanTJansmaEPvan WaesbergheJTMBoellaardRPuntCJAet al. Advanced analytics and artificial intelligence in gastrointestinal cancer: a systematic review of radiomics predicting response to treatment. Eur J Nucl Med Mol Imaging (2021) 48(6):1785–94. doi: 10.1007/s00259-020-05142-w

  • 11

    ZhongYSheYDengJChenSWangTYangMet al. Deep learning for prediction of N2 metastasis and survival for clinical stage I non-small cell lung cancer. Radiology (2022) 302(1):200–11. doi: 10.1148/radiol.2021210902

  • 12

    LiYZhangYFangQZhangXHouPWuHet al. Radiomics analysis of [18F]FDG PET/CT for microvascular invasion and prognosis prediction in very-early- and early-stage hepatocellular carcinoma. Eur J Nucl Med Mol Imaging (2021) 48(8):2599–614. doi: 10.1007/s00259-020-05119-9

  • 13

    FangJZhangBWangSJinYWangFDingYet al. Association of MRI-derived radiomic biomarker with disease-free survival in patients with early-stage cervical cancer. Theranostics (2020) 10(5):2284–92. doi: 10.7150/thno.37429

  • 14

    KirienkoMSolliniMCorbettaMVoulazEGozziNInterlenghiMet al. Radiomics and gene expression profile to characterise the disease and predict outcome in patients with lung cancer. Eur J Nucl Med Mol Imaging (2021) 48(11):3643–55. doi: 10.1007/s00259-021-05371-7

  • 15

    LiuCGongJYuHLiuQWangSWangJ. A CT-based radiomics approach to predict nivolumab response in advanced non-Small-Cell lung cancer. Front Oncol (2021) 11:544339. doi: 10.3389/fonc.2021.544339

  • 16

    GuoYSongQJiangMGuoYXuPZhangYet al. Histological subtypes classification of lung cancers on CT images using 3D deep learning and radiomics. Acad Radiol (2021) 28(9):e258–66. doi: 10.1016/j.acra.2020.06.010

  • 17

    ZhaoHLiWLyuPZhangXLiuHLiangPet al. TCGA-TCIA-Based CT radiomics study for noninvasively predicting Epstein-Barr virus status in gastric cancer. AJR Am J Roentgenol (2021) 217(1):124–34. doi: 10.2214/AJR.20.23534

  • 18

    LarueRTHMvan TimmerenJEde JongEECFelicianiGLeijenaarRTHSchreursWMJet al. Influence of gray level discretization on radiomic feature stability for different CT scanners, tube currents and slice thicknesses: a comprehensive phantom study. Acta Oncol (2017) 56(11):1544–53. doi: 10.1080/0284186X.2017.1351624

  • 19

    FedorovABeichelRKalpathy-CramerJFinetJFillion-RobinJCPujolSet al. 3D slicer as an image computing platform for the quantitative imaging network. Magn Reson Imaging (2012) 30(9):1323–41. doi: 10.1016/j.mri.2012.05.001

  • 20

    MayerhoeferMEMaterkaALangsGHäggströmISzczypińskiPGibbsPet al. Introduction to radiomics. J Nucl Med (2020) 61(4):488–95. doi: 10.2967/jnumed.118.222893

  • 21

    van GriethuysenJJMFedorovAParmarCHosnyAAucoinNNarayanVet al. Computational radiomics system to decode the radiographic phenotype. Cancer Res (2017) 77(21):e104–7. doi: 10.1158/0008-5472.CAN-17-0339

  • 22

    FitzgeraldMSavilleBRLewisRJ. Decision curve analysis. JAMA (2015) 313(4):409–10. doi: 10.1001/jama.2015.37

  • 23

    VaidyaPBeraKGuptaAWangXCorredorGFuPet al. CT derived radiomic score for predicting the added benefit of adjuvant chemotherapy following surgery in stage I, II resectable non-small cell lung cancer: a retrospective multi-cohort study for outcome prediction. Lancet Digit Health (2020) 2(3):e116–28. doi: 10.1016/s2589-7500(20)30002-9

  • 24

    LiBCuiYDiehnMLiR. Development and validation of an individualized immune prognostic signature in early-stage nonsquamous non-small cell lung cancer. JAMA Oncol (2017) 3(11):1529–37. doi: 10.1001/jamaoncol.2017.1609

  • 25

    DevarakondaSRotoloFTsaoMSLancIBrambillaEMasoodAet al. Tumor mutation burden as a biomarker in resected non-Small-Cell lung cancer. J Clin Oncol (2018) 36(30):29953006. doi: 10.1200/JCO.2018.78.1963

  • 26

    LuJWuJLouYShiQXuJZhangLet al. Blood-based tumour mutation index act as prognostic predictor for immunotherapy and chemotherapy in non-small cell lung cancer patients. biomark Res (2022) 10(1):55. doi: 10.1186/s40364-022-00400-5

  • 27

    GuoNLDowlatiARaeseRADongCChenGBeerDGet al. A predictive 7-gene assay and prognostic protein biomarkers for non-small cell lung cancer. EBioMedicine (2018) 32:102–10. doi: 10.1016/j.ebiom.2018.05.025

  • 28

    HuangYLiuZHeLChenXPanDMaZet al. Radiomics signature: a potential biomarker for the prediction of disease-free survival in early-stage (I or II) non-small cell lung cancer. Radiology (2016) 281(3):947–57. doi: 10.1148/radiol.2016152234

  • 29

    XieDWangTTHuangSJDengJJRenYJYangYet al. Radiomics nomogram for prediction disease-free survival and adjuvant chemotherapy benefits in patients with resected stage I lung adenocarcinoma. Transl Lung Cancer Res (2020) 9(4):1112–23. doi: 10.21037/tlcr-19-577

  • 30

    MuWJiangLZhangJShiYGrayJETunaliIet al. Non-invasive decision support for NSCLC treatment using PET/CT radiomics. Nat Commun (2020) 11(1):5228. doi: 10.1038/s41467-020-19116-x

  • 31

    MuWTunaliIGrayJEQiJSchabathMBGilliesRJ. Radiomics of 18F-FDG PET/CT images predicts clinical benefit of advanced NSCLC patients to checkpoint blockade immunotherapy. Eur J Nucl Med Mol Imaging (2020) 47(5):1168–82. doi: 10.1007/s00259-019-04625-9

  • 32

    ZhangJZhaoXZhaoYZhangJZhangZWangJet al. Value of pre-therapy 18F-FDG PET/CT radiomics in predicting EGFR mutation status in patients with non-small cell lung cancer. Eur J Nucl Med Mol Imaging (2020) 47(5):1137–46. doi: 10.1007/s00259-019-04592-1

  • 33

    MuWJiangLShiYTunaliIGrayJEKatsoulakisEet al. Non-invasive measurement of PD-L1 status and prediction of immunotherapy response using deep learning of PET/CT images. J Immunother Cancer (2021) 9(6):e002118. doi: 10.1136/jitc-2020-002118

  • 34

    MonacoLDe BernardiEBonoFCortinovisDCrivellaroCEliseiFet al. The "digital biopsy" in non-small cell lung cancer (NSCLC): a pilot study to predict the PD-L1 status from radiomics features of [18F]FDG PET/CT. Eur J Nucl Med Mol Imaging (2022) 49(10):3401–11. doi: 10.1007/s00259-022-05783-z

  • 35

    WengQZhouLWangHHuiJChenMPangPet al. A radiomics model for determining the invasiveness of solitary pulmonary nodules that manifest as part-solid nodules. Clin Radiol (2019) 74(12):933–43. doi: 10.1016/j.crad.2019.07.026

  • 36

    Thattaamuriyil PadmakumariLGuidoGCarusoDNacciIDel GaudioAZerunianMet al. The role of chest CT radiomics in diagnosis of lung cancer or tuberculosis: a pilot study. Diagnostics (Basel) (2022) 12(3):739. doi: 10.3390/diagnostics12030739

  • 37

    QiuLZhangXMaoHFangXDingWZhaoLet al. Comparison of comprehensive morphological and radiomics features of subsolid pulmonary nodules to distinguish minimally invasive adenocarcinomas and invasive adenocarcinomas in CT scan. Front Oncol (2022) 11:691112. doi: 10.3389/fonc.2021.691112

  • 38

    LiXGuindaniMNgCSHobbsBP. Spatial Bayesian modeling of GLCM with application to malignant lesion characterization. J Appl Stat (2018) 46(2):230–46. doi: 10.1080/02664763.2018.1473348

  • 39

    GranataVFuscoRDe MuzioFCutoloCSetolaSVGrassiRet al. Radiomics textural features by MR imaging to assess clinical outcomes following liver resection in colorectal liver metastases. Radiol Med (2022) 127(5):461–70. doi: 10.1007/s11547-022-01477-6

  • 40

    WangTDengJSheYZhangLWangBRenYet al. Radiomics signature predicts the recurrence-free survival in stage I non-small cell lung cancer. Ann Thorac Surg (2020) 109(6):1741–9. doi: 10.1016/j.athoracsur.2020.01.010

  • 41

    YangLYangJZhouXHuangLZhaoWWangTet al. Development of a radiomics nomogram based on the 2D and 3D CT features to predict the survival of non-small cell lung cancer patients. Eur Radiol (2019) 29(5):2196–206. doi: 10.1007/s00330-018-5770-y

  • 42

    WangXDuanHLiXYeXHuangGNieS. A prognostic analysis method for non-small cell lung cancer based on the computed tomography radiomics. Phys Med Biol (2020) 65(4):045006. doi: 10.1088/1361-6560/ab6e51

  • 43

    Artal CortésÁCalera UrquizuLHernando CuberoJ. Adjuvant chemotherapy in non-small cell lung cancer: state-of-the-art. Transl Lung Cancer Res (2015) 4(2):191–7. doi: 10.3978/j.issn.2218-6751.2014.06.01

Summary

Keywords

computed tomography, radiomics, non-small cell lung cancer, prognosis, LASSO

Citation

Peng B, Wang K, Xu R, Guo C, Lu T, Li Y, Wang Y, Wang C, Chang X, Shen Z, Shi J, Xu C and Zhang L (2023) Preoperative computed tomography-based tumoral radiomic features prediction for overall survival in resectable non-small cell lung cancer. Front. Oncol. 13:1131816. doi: 10.3389/fonc.2023.1131816

Received

26 December 2022

Accepted

24 April 2023

Published

03 May 2023

Volume

13 - 2023

Edited by

Hongda Liu, Nanjing Medical University, China

Reviewed by

Nika Guberina, University Duisburg-Essen, Germany; Burcin Celik, Ondokuz Mayıs University, Türkiye

Updates

Copyright

*Correspondence: Linyou Zhang,

†These authors have contributed equally to this work

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics