[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Skip to main content

CT-based deep learning radiomics biomarker for programmed cell death ligand 1 expression in non-small cell lung cancer

Abstract

Background

Programmed cell death ligand 1 (PD-L1), as a reliable predictive biomarker, plays an important role in guiding immunotherapy of lung cancer. To investigate the value of CT-based deep learning radiomics signature to predict PD-L1 expression in non-small cell lung cancers(NSCLCs).

Methods

259 consecutive patients with pathological confirmed NSCLCs were retrospectively collected and divided into the training cohort and validation cohort according to the chronological order. The univariate and multivariate analyses were used to build the clinical model. Radiomics and deep learning features were extracted from preoperative non-contrast CT images. After feature selection, Radiomics score (Rad-score) and deep learning radiomics score (DLR-score) were calculated through a linear combination of the selected features and their coefficients. Predictive performance for PD-L1 expression was evaluated via the area under the curve (AUC) of receiver operating characteristic, the calibration curves, and the decision curve analysis.

Results

The clinical model based on Cytokeratin 19 fragment and lobulated shape obtained an AUC of 0.767(95% CI: 0.673–0.860) in the training cohort and 0.604 (95% CI:0.477–0.731) in the validation cohort. 11 radiomics features and 15 deep learning features were selected by LASSO regression. AUCs of the Rad-score were 0.849 (95%CI: 0.783–0.914) and 0.717 (95%CI: 0.607–0.826) in the training cohort and validation cohort, respectively. AUCs of DLR-score were 0.938 (95%CI: 0.899–0.977) and 0.818(95%CI:0.727–0.910) in the training cohort and validation cohort, respectively. AUCs of the DLR-score were significantly higher than those of the Rad-score and the clinical model.

Conclusion

The CT-based deep learning radiomics signature could achieve clinically acceptable predictive performance for PD-L1 expression, which showed potential to be a surrogate imaging biomarker or a complement of immunohistochemistry assessment.

Peer Review reports

Background

Lung cancer is the second most commonly diagnosed cancer worldwide [1, 2]. Although early-stage detection through low-dose CT screening and mini-invasive video-assisted thoracoscopic surgery have improved patients’ survival and life quality greatly, lung cancer remains the leading cause of cancer-related death due to approximately 80% of lung cancers being diagnosed at advanced stage, which are unresectable and systemic chemotherapy is the only option [3, 4]. Only 20-40% of patients response to the standard platinum-based chemotherapy [5]. With the development of immunotherapy, the management of lung cancers has evolved enormously recently. Non-small cell lung cancer (NSCLC) accounts for about 85% in all lung cancers, including the most common subtypes such as lung adenocarcinoma and lung squamous cell carcinoma [6]. The first-generation antibody-based immunotherapy, which targets at blocking the receptor and/or ligand interactions of molecules, such as programmed cell death protein 1 (PD-1) and its ligand (PD-L1) or cytotoxic T lymphocyte antigen-4 can modulate antitumor responses, had shown remarkably response durable in NSCLCs [7]. Unfortunately, only 17-21% of patients with NSCLCs demonstrated a response to anti–PD-1 or PD-L1 therapy [8, 9]. Therefore, acknowledgment of which patients would benefit from immune checkpoint inhibitors (ICIs) is needed in NSCLCs treatment strategies. Many studies revealed that tumor mutational burden and PD-L1 expression were independent predictive factors for the response of ICIs [10, 11]. However, acknowledging these biomarkers requires invasive procedures to obtain tumor tissue specimens for gene sequencing or immunohistochemistry (IHC) staining, which are time-consuming and expensive. Furthermore, obtaining tissue specimens is difficult, even impossible in most clinical scenarios for patients with advanced NSCLC.

As a noninvasive technique, CT has been widely implemented in the diagnosis, staging, treatment planning, and response assessment of NSCLCs through radiologists’ visual interpreting, which uses only a few metrics of imaging. The development of computer science and artificial intelligent results into the emergence of radiomics, which extracts high-dimensional features from medical imaging data to decode imaging phenotype to achieve comprehensive clinical goals [12]. In recent years, convolutional neural network (CNN) with multiple network structures has been widely used in radiological tumor research. They can extract a large number of useful deep learning (DL) features for tumor grade prediction, lymph node metastasis prediction and risk prognosis prediction [13,14,15]. However, the construction of CNN often requires a large number of samples, and most medical studies often have a small sample size. Therefore, transfer learning is widely used in the field of medical deep learning, which can alleviate the limitation of small data sets [16]. Transfer learning involves the use of pre-trained neural networks on other images and allows existing training models to be applied to unsolved problems, thus greatly reducing the need for a large amount of training data. The purpose of this study was to develop deep learning radiomics signature as a surrogate imaging biomarker for PD-L1 expression of NSCLC, using transfer learning to extract features from CT images, in order to provide decision-making support for selecting patients who would benefit from ICIs treatment.

Materials and methods

Patients

This study was approved by the Local Hospital Ethics Committee of our hospital with the waiver of informed consent due to the retrospective nature. We searched our hospital’s database for patients who had pathological diagnosis of NSCLC via biopsy or surgery receiving CT examinations being performed within 3 months before biopsy or surgery and IHC staining for PD-L1 expression. Exclusion criteria were as follows: (1) insufficient image quality for nodule segmenting (n = 5); (2) undergoing anti-tumor therapy (radiotherapy, chemotherapy or chemoradiotherapy) before biopsy or surgery (n = 7); (3) Radiomics feature or deep learning feature extraction failed (n = 30). At first, 301 deemed eligible patients were identified. After excluding 42 patients, the final study cohort included 259 patients.who were divided into a training cohort and a validation cohort at a ratio of 7:3 according to the chronological order (Fig. 1). Patients before October 2021 were assigned to training sets(PD-L1 negative n = 131, PD-L1 positive n = 32), and patients after that wereassigned to validation sets (PD-L1 negative n = 67, PD-L1 positive n = 29). There was no significant difference in the distribution of PD-L1 expression in the training set and the validation set (P = 0.053), which could be used for model establishment and validation.

Fig. 1
figure 1

Flowchart shows inclusion and exclusion criteria for the study

CT acquisition and interpretation

Preoperative chest CT examinations were performed on CT scanners (Brilliance iCT, Philips Medical Systems). The acquisition parameters were as follows: 0.625 mm x 128 of a collimation; 120 kVp of tuber voltage; automatic exposure control (AEC) of tube current; a reconstruction slice thickness of 1.5 mm and a gap of 1 mm; field of view of 350 × 350 mm; matrix of 512 × 512.

Two experienced radiologists (2 years and 5 years chest CT interpretation experience), who were blinded to the clinical and PD-L1 expression, interpreted the thin slice CT images in the lung window setting to obtain the semantic features of nodules on the PACS (Picture Archiving and Communication Systems). In the semantic description of nodules, CT characteristics as following were included: vacuolar sign, cavity, pleural thickening, pleural indentation, hilar adenopathy, mediastinal adenopathy, vessel convergence, location, lobulated shape, spiculation, airbronchial sign, types of nodules and size.

Feature extraction, feature selection and signature construction

The thin slice images in Digital Imaging and Communication in Medicine (DICOM) format derived from PACS were transferred to ITK-SNAP 3.8.0 (http://www.itksnap.org) of a personal computer for tumor segmentation and radiomics feature extraction. The regions of interest (ROIs) of tumor were manually drew using ITK-SNAP software on each thin slice to convert a three-dimensional volume of interest (3D-VOIs). Two radiologists with 2 years and 5 years of experience in thoracic CT interpretation segmented tumors in 30 randomly selected patients independently. The radiologist with 2 years of experience segmented all the nodules manually. 1834 radiomics features were extracted from the 3D-VOIs using the open-source software package Pyradiomics (https://github.com/Radiomics/pyradiomics).

A pre-trained CNN, ResNet 50, was used for transfer learning to extract deep learning features from thin CT images of NSCLCs. First, the image with the largest tumor area per patient was selected and the grayscale values were normalized into the range [− 1,1] using a min-max transformation. Then each cropped subregion image was resized to 224 × 224, and the resulting image was used as model input [17, 18].

Z-Score normalization was implemented to reduce the influence of features’ scales. Intraclass correlation coefficients (ICC) of interobservers were implemented to roll out those radiomics or deep learning features with low repeatability (ICC ≤ 0.75). Pearson correlation was performed to exclude the radiomics or deep learning features with high correlations (r > 0.90). Then the retained features were introduced into the least absolute shrinkage and selection operator (LASSO) regression with 5-fold cross-validation to select radiomics or deep learning features which were strong association with PD-L1. Radiomics score (Rad-score) and deep learning radiomics score (DLR-score) for each patient were calculated through a linear combination of the selected features weighted by their coefficient. Feature selection procedure was implemented on both radiomics features and deep learning features in the training cohort. The signatures trained on the training cohort were applied to the validation cohort for testing in independent cases. The workflow of model building is shown in Fig. 2.

Fig. 2
figure 2

The models workflow

PD-L1 testing

PD-L1 expression of NSCLCs was assessed using immunohistochemical staining and reported as tumor proportion score (TPS). TPS is defined as the percentage of tumor cells stained with PD-L1 membrane of any intensity. The PD-L1 expression was dichotomized according to TPS level (TPS < 1% is negative, TPS ≥ 1% is positive). Finally, there were 198 PD-L1 negative cases and 61 PDL1 positive cases.

Statistical analysis and Model Development

For clinical metrics, univariate analysis (t-test, Mann-Whitney U rank test, χ2 test, or Fisher’s precise probability test) was used to select those which were related to PD-L1 expression, and then introduced into a multivariate logistic regression model. The predictive performances of the clinical model, the Rad-score and the DLR-score were assessed by area under the curve (AUC) of receiver operating characteristic (ROC) which were compared using the DeLong method. The calibration effectiveness, the goodness of fit, the net benefit and the clinical effectiveness of the better model were evaluated using the calibration curve, the Hosmer-Lemeshow and the decision curve. ROC analysis of models was performed to obtain the optimal cut-off value [19, 20].

Results

Patient characteristics and clinical model

Demographic characteristics of all patients were showed in Table 1. At univariable analysis, Cytokeratin 19 fragment (P = 0.047) and lobulated shape(P = 0.014) were related to PD-L1 with statistical significance (Table 2) in training cohort. Multivariate logistic regression revealed that Cytokeratin 19 fragment and lobulated shape were independent risk factors for PD-L1 (Table 2). A clinical model based on these two clinical features yielded predictive performance with an AUC of 0.767(95% CI: 0.673–0.860) for PD-L1 in the training cohort and 0.604 (95% CI:0.477–0.731) in validation cohort.

Table 1 Demographic and CT characteristics of all patients
Table 2 The univariate and multivariate logistic regression of PD-L1 expression based on clinical and CT characteristics in training cohort

Radiomics score and deep learning Radiomics score building and evaluation

After ICC and pearson analysis, 205 radiomics features were introduced into the LASSO model. 11 features were selected to establish the Rad-score as [21]: Rad-score = 0.19631901840490795 + 0.002025*exponential_glszm_SmallAreaEmphasis + 0.022456*lbp_3D_k_ngtdm_Busyness + 0.030649*lbp_3D_m1_glszm_SizeZoneNonUniformity + 0.028801*lbp_3D_m1_glszm_SmallAreaLowGrayLevelEmphasis + 0.009714*lbp_3D_m2_firstorder_Median-0.032295*log_sigma_2_0_mm_3D_glszm_Zone%-0.028661*squareroot_firstorder_Skewness + 0.024971*wavelet_HHL_glszm_ZoneVariance + 0.020694*wavelet_HLL_glcm_ClusterProminence-0.002185*wavelet_LHL_firstorder_Median + 0.019159 * wavelet_LLL_firstorder_Kurtosis. In deep learning, 2048 DL radiomics features were extracted. After ICC and Pearson analysis, 7 features were found to be excluded. Then LASSO regression revealed that 19 DL features were strongly associated with PD-L1. Following a rule of thumb, we have retained the 15 characteristics that are most relevant to PD-L1 [22]. The DLR-score was calculated as [17]: DLR-score = 0.19631901840490795- 0.029073 * DL_85 -0.008362 * DL_87 + 0.019473 * DL_148 -0.001373 * DL_263 -0.034845 * DL_304 -0.030973 * DL_414 -0.005602 * DL_456 -0.002396 * DL_642 -0.004188 * DL_780 + 0.009646 * DL_819 + 0.018088 * DL_899 -0.013829 * DL_1011 -0.006451 * DL_1038-0.012642 * DL_1252 -0.008234 * DL_1681. AUCs of Rad-score for PD-L1 expression were 0.849 (95%CI: 0.783–0.914) and 0.717 (95%CI: 0.607–0.826) in the training cohort and validation cohort, respectively. AUCs of DLR-score for PD-L1 expression were 0.938(95%CI: 0.899–0.977) and 0.818(95%CI:0.727–0.910) in the training cohort and validation cohort, respectively. The predictive performance of the DLR-score was significantly higher than that of the clinical model and Rad-score in both cohorts (Fig. 3). The results of the ROC analysis of each model are shown in Table 3. The Youden index of the ROC curve determined the optimal cut-off value. The optimal cut-off value derived from DLR-score was 0.246, with sensitivity, specificity, PPV, NPV and accuracy, of 90.6%, 87.8%, 64.4%, 97.5%, and 88.3%, respectively. The comparison of clinical and CT features between high (more than cut-off value) and low ( less than cut-off value) DLR-score was shown in Table 4. Intending to examine the interpretability of DL features, we also visualize the network by applying a gradient-weighted class activation mapping that can provide a rough positioning map to highlight important areas of the classification target (Figs. 4 and 5). The calibration curves of the radiomics model showed a good calibration effect on the predictive efficacy of PD-L1 expression in the training cohort (Fig. 6A), and the Hosmer-Lemeshow test showed nonsignificant statistic in the training cohort (P = 0.141), indicating that there was no significant difference between prediction and pathology result. The decision curve analysis shows that radiomics model achieved a high net benefit at most probability thresholds, indicating that DLR-score could achieve excellent clinical effectiveness when the probability of threshold is approximately between 0% and about 80% in the training cohort (Fig. 6B).

Table 3 The performance of deep transfer learning radiomics signature to predict PD-L1 expression
Table 4 The comparison of clinical and CT features between high DLRad-score and low DLRad-score
Fig. 3
figure 3

The ROC curves of clinical model, Radscore, and DLR-score in training cohort (A) and validation cohort (B), respectively

Fig. 4
figure 4

A 51-year-old man with NSCLC. (A) Axial CT image shows a solid nodule of the left upper lobe. DLR-score 0.839. (B) The Photomicrograph shows a positive expression of PD-L1(IHC; x400). (C) Grad-CAM visualization. Grad-CAM, gradient-weighted class activation mapping

Fig. 5
figure 5

A 78-year-old man with NSCLC. (A) Axial CT image shows a ground glass nodule of the left lower lobe, DLR-score 0.141. (B) The Photomicrograph shows a negative expression of PD-L1(IHC; x400). (C) Grad-CAM visualization. Grad-CAM, gradient-weighted class activation mapping

Fig. 6
figure 6

The calibration curves of the DLR-score in the training cohort (A). The decision curve analyses of the three models in training cohort (B). The DLR-score can achieve a greater net effect than the other two models

Discussion

In this study, we probed whether deep learning signature derived from preoperative CT could be used to predict PD-L1 expression. We developed and validated three models, clinical, radiomics, and DLR-score, for predicting PD-L1 expression by quantitative analysis of CT images of NSCLCs. In the training and validation cohorts, the DLR-score showed the best predictive performance compared to other models. The AUCs of DLR-score were 0.938(95%CI: 0.899–0.977) in the training cohort and 0.818(95%CI:0.727–0.910) in the validation cohort. DCA showed that the DLR-score can improve the predictive performance of PD-L1 expression. The high predictive performance of the DLR-score showed the possibility to be a noninvasive surrogate biomarker for PD-L1 facilitating the selection of patients who would benefit from ICIs treatment.

PD-L1 plays an important role in guiding immunotherapy of lung cancer. When the PD-L1 of tumor cells binds to the PD-1 on the surface of immune cells, a negative immune response will occur, resulting in the escape of tumor cells and promoting the occurrence, development and metastasis of tumors [23]. PD-1/PD-L1 inhibitors kill tumor cells by blocking the binding of the PD-1/PD-L1 pathway, relieving the negative regulation of immune cells and preventing immune escape [24]. PD-1/PD-L1 inhibitors were included in the first category of recommendations in NSCLC’s NCCN guidelines [25, 26]. When stratified analysis of PD-L1 in tumor tissue, it was found that PD-1/PD-L1 inhibitor was more effective in patients with PD-L1 positive [27,28,29]. Therefore, there is an urgent need for screening patients who might be most likely to benefit from ICIs treatment. Recent studies showed that PD-L1 expression, tumor mutational burden and tumor immune microenvironment can be used as predictive biomarkers to predict the response of ICIs treatment, but these require not only invasive procedures to obtain tissue specimens, but also time-consuming and expensive laboratory tests [7, 30]. Therefore, these biomarkers are not broadly available in clinical scenarios, especially in patients with advanced NSCLCs for tissue specimens are difficult to obtain sometimes, even impossible. Furthermore, malignancies are heterogenetic, and tissue specimens, especially those through biopsy, may harbor sample errors. Therefore, clinicians still confront the challenge to choose suitable patients for ICIs treatment.

Radiomics is a data-driven discipline based on widely available imaging data that can be used to improve diagnosis, prognosis, and clinical decision support [31]. In order to provide decision support of ICIs treatment, several studies have investigated the association between radiomics signature with PD-L1 expression and tumor immune microenvironment in several kinds of solid tumors [32,33,34]. Regarding PD-L1 expression in NSCLC, Jiang et al [35] derived radiomics signatures from CT, PET, and PET/CT, which achieved predictive performance to identify PD-L1 expression over 1% with AUC of 0.86, 0.62, and 0.85, respectively. Using the same algorithm, Sun et al [36] reported preoperative CT-derived radiomics signature obtained AUCs of 0.786 and 0.807 in the training and validation cohort, respectively. When combined with clinicopathological features, the predictive performance increased to 0.829 and 0.848, respectively. However, the authors did not evaluate whether the difference was statistically significant. Our present study achieved a similar predictive performance with these two studies. In our study, the AUC of our clinical model and the radiomics score were 0.767 (95%CI: 0.673–0.860) and 0.849 (95% CI: 0.783–0.914), respectively, which was similar with the results of the above two studies.

Deep learning has made significant progress in the field of medical image analysis by mining high-throughput information from medical images for recognition of images or prediction of gene expression [37,38,39,40]. Wang et al. developed an end-to-end deep learning model based on CT images, which obtained AUCs of 0.85 (95% CI 0.83–0.88) in the primary cohort and 0.81 (95% CI 0.79–0.83) in the independent validation cohort to predict EGFR mutation status in lung cancers, respectively [38]. Chen et al. developed a deep learning-based method for the automatic segmentation of meningiomas from multiparametric MR images, and the AUC of the radiomics model with automatic segmentation was comparable to the AUC of the manual segmentation model in the internal (0.95 vs. 0.93, p = 0.176) and external (0.88 vs. 0.91, p = 0.419) test cohort [37]. In this study, a pre-trained CNN, ResNet 50, was implemented to extract 2048 deep learning features from CT images of NSCLCs, and 15 features were found to be strongly associated with PD-L1 expression. Compared to the clinical model and radiomics model, the DLR-score demonstrated the highest predictive performance for PD-L1 expression with AUCs of 0.938(95%CI: 0.899–0.977) and 0.818(95%CI:0.727–0.910) in the training cohort and the validation cohort, respectively. The optimal cut-off value derived from DLR-score was 0.246 achieved sensitivity, specificity, PPV, NPV, and accuracy, of 90.6%, 87.8%, 64.4%, 97.5%, and 88.3%, respectively, which were more improved than the other two models. The decision curve analysis also showed that the clinical net benefit of the DLR-score was higher than that of clinical models and radiomics score, both in the training and validation cohort.

Several limitations of this study need to be acknowledged. First, this was a retrospective study with a small sample size and no external validation cohort. Second, the imbalance distribution between PD-L1 positive and negative expression may impact on the predictive performance of the model. Third, to avoid overfitting, transfer learning often requires a large sample size, and the sample size in this study was clearly not sufficient for 3D analysis, so we used the image of the largest tumor area, rather than using 3D whole tumor volume to extract DL features [17]. However, this approach is time-saving and may be more clinically appropriate. Finally, owing to surgical confirmed cases at an early stage and a short follow-up period, only several patients received ICIs treatment after surgery due to recurrence or metastasis at present. Therefore, the predictive performance of the CT-based deep learning radiomics signature for treatment response was not evaluated at this study.

Conclusion

In conclusion, this study developed clinical, radiomics and deep learning models to predict PD-L1 expression in NSCLCs non-invasively. It showed the CT-based deep learning radiomics model could achieve clinically acceptable predictive performance in both training and validation cohorts. The deep learning radiomics signature could offer a surrogate imaging biomarker or a complement for IHC analysis, which could facilitate clinical decision support in identifying NSCLC patients who are likely to benefit from ICIs treatment.

Data availability

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

Abbreviations

3D-VOIs:

Three-dimensional volume of interest

AUC:

Area under the curve

CNN:

Convolutional neural network

DICOM:

Digital Imaging and Communications in Medicine

DL:

Deep learning

DLR-score:

Deep learning radiomics score

ICCs:

Intraclass correlation coefficients

ICIs:

Immune checkpoint inhibitors

IHC:

Immunohistochemistry

LASSO:

Least absolute shrinkage and selection operator

NSCLC:

Non-small cell lung cancer

PACS:

Picture archiving and communication system

PD-1:

Programmed cell death protein 1

PD-L1:

Programmed cell death ligand 1

Rad-score:

Radiomics score

ROI:

Regions of interest

TPS:

Tumor proportion score

References

  1. Xia C, Dong X, Li H, Cao M, Sun D, He S, Yang F, Yan X, Zhang S, Li N, et al. Cancer statistics in China and United States, 2022: profiles, trends, and determinants. Chin Med J (Engl). 2022;135(5):584–90.

    Article  PubMed  Google Scholar 

  2. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, Bray F. Global Cancer statistics 2020: GLOBOCAN estimates of incidence and Mortality Worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209–49.

    Article  PubMed  Google Scholar 

  3. Sands J, Tammemagi MC, Couraud S, Baldwin DR, Borondy-Kitts A, Yankelevitz D, Lewis J, Grannis F, Kauczor HU, von Stackelberg O, et al. Lung screening benefits and challenges: a review of the data and outline for implementation. J Thorac Oncol. 2021;16(1):37–53.

    Article  PubMed  Google Scholar 

  4. Houston KA, Henley SJ, Li J, White MC, Richards TB. Patterns in lung cancer incidence rates and trends by histologic type in the United States, 2004–2009. Lung Cancer. 2014;86(1):22–8.

    Article  PubMed  Google Scholar 

  5. Santabarbara G, Maione P, Rossi A, Palazzolo G, Gridelli C. The role of pembrolizumab in the treatment of advanced non-small cell lung cancer. Ann Transl Med. 2016;4(11):215.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Herbst RS, Morgensztern D, Boshoff C. The biology and management of non-small cell lung cancer. Nature. 2018;553(7689):446–54.

    Article  CAS  PubMed  Google Scholar 

  7. Rizvi H, Sanchez-Vega F, La K, Chatila W, Jonsson P, Halpenny D, Plodkowski A, Long N, Sauter JL, Rekhtman N, et al. Molecular determinants of response to Anti-programmed Cell death (PD)-1 and anti-programmed death-ligand 1 (PD-L1) blockade in patients with non-small-cell Lung Cancer profiled with targeted next-generation sequencing. J Clin Oncol. 2018;36(7):633–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Topalian SL, Hodi FS, Brahmer JR, Gettinger SN, Smith DC, McDermott DF, Powderly JD, Carvajal RD, Sosman JA, Atkins MB, et al. Safety, activity, and immune correlates of anti-PD-1 antibody in cancer. N Engl J Med. 2012;366(26):2443–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Garon EB, Ciuleanu TE, Arrieta O, Prabhash K, Syrigos KN, Goksel T, Park K, Gorbunova V, Kowalyszyn RD, Pikiel J, et al. Ramucirumab plus Docetaxel versus placebo plus docetaxel for second-line treatment of stage IV non-small-cell lung cancer after disease progression on platinum-based therapy (REVEL): a multicentre, double-blind, randomised phase 3 trial. Lancet. 2014;384(9944):665–73.

    Article  CAS  PubMed  Google Scholar 

  10. Kluger HM, Zito CR, Turcu G, Baine MK, Zhang H, Adeniran A, Sznol M, Rimm DL, Kluger Y, Chen L, et al. PD-L1 studies across Tumor types, its Differential expression and predictive value in patients treated with Immune Checkpoint inhibitors. Clin Cancer Res. 2017;23(15):4270–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Shi WJ, Zhao W. Biomarkers or factors for predicting the efficacy and adverse effects of immune checkpoint inhibitors in lung cancer: achievements and prospective. Chin Med J (Engl). 2020;133(20):2466–75.

    Article  CAS  PubMed  Google Scholar 

  12. Aerts HJ, Velazquez ER, Leijenaar RT, Parmar C, Grossmann P, Carvalho S, Bussink J, Monshouwer R, Haibe-Kains B, Rietveld D, et al. Decoding tumour phenotype by noninvasive imaging using a quantitative radiomics approach. Nat Commun. 2014;5:4006.

    Article  CAS  PubMed  Google Scholar 

  13. Liu Y, Balagurunathan Y, Atwater T, Antic S, Li Q, Walker RC, Smith GT, Massion PP, Schabath MB, Gillies RJ. Radiological image traits Predictive of Cancer Status in Pulmonary nodules. Clin Cancer Res. 2017;23(6):1442–9.

    Article  PubMed  Google Scholar 

  14. Ehteshami Bejnordi B, Veta M, van Johannes P, van Ginneken B, Karssemeijer N, Litjens G, van der Laak J, the, Hermsen CC, Manson M et al. QF : Diagnostic Assessment of Deep Learning Algorithms for Detection of Lymph Node Metastases in Women With Breast Cancer. JAMA 2017, 318(22):2199–2210.

  15. Ardila D, Kiraly AP, Bharadwaj S, Choi B, Reicher JJ, Peng L, Tse D, Etemadi M, Ye W, Corrado G, et al. End-to-end lung cancer screening with three-dimensional deep learning on low-dose chest computed tomography. Nat Med. 2019;25(6):954–61.

    Article  CAS  PubMed  Google Scholar 

  16. Christodoulidis S, Anthimopoulos M, Ebner L, Christe A, Mougiakakou S. Multisource transfer learning with convolutional neural networks for lung pattern analysis. IEEE J Biomed Health Inf. 2017;21(1):76–84.

    Article  Google Scholar 

  17. Zheng Y-m, Che J-y, Yuan M-g, Wu Z-j, Pang J, Zhou R-z, Li X-l, Dong C. A CT-Based Deep Learning Radiomics Nomogram to predict histological grades of Head and Neck squamous cell carcinoma. Acad Radiol. 2023;30(8):1591–9.

    Article  PubMed  Google Scholar 

  18. Zeng Q, Li H, Zhu Y, Feng Z, Shu X, Wu A, Luo L, Cao Y, Tu Y, Xiong J et al. Development and validation of a predictive model combining clinical, radiomics, and deep transfer learning features for lymph node metastasis in early gastric cancer. Front Med 2022, 9.

  19. Zhang Y, Ko C-C, Chen J-H, Chang K-T, Chen T-Y, Lim S-W, Tsui Y-K, Su M-Y. Radiomics Approach for Prediction of Recurrence in Non-functioning Pituitary Macroadenomas. Front Oncol. 2020;10:590083.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Zhu F, Yang C, Xia Y, Wang J, Zou J, Zhao L, Zhao Z. CT-based radiomics models may predict the early efficacy of microwave ablation in malignant lung tumors. Cancer Imaging: Official Publication Int Cancer Imaging Soc. 2023;23(1):60.

    Article  Google Scholar 

  21. Wang Y, Bi Q, Deng Y, Yang Z, Song Y, Wu Y, Wu K. Development and validation of an MRI-based Radiomics Nomogram for Assessing Deep Myometrial Invasion in Early Stage Endometrial Adenocarcinoma. Acad Radiol. 2023;30(4):668–79.

    Article  PubMed  Google Scholar 

  22. Halligan S, Menu Y, Mallett S. Why did European Radiology reject my radiomic biomarker paper? How to correctly evaluate imaging biomarkers in a clinical setting. Eur Radiol. 2021;31(12):9361–8.

    Article  PubMed  PubMed Central  Google Scholar 

  23. He J, Hu Y, Hu M, Li B. Development of PD-1/PD-L1 pathway in Tumor Immune Microenvironment and Treatment for Non-small Cell Lung Cancer. Sci Rep. 2015;5:13110.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Steven A, Fisher SA, Robinson BW. Immunotherapy for lung cancer. Respirology. 2016;21(5):821–33.

    Article  PubMed  Google Scholar 

  25. Herbst RS, Giaccone G, de Marinis F, Reinmuth N, Vergnenegre A, Barrios CH, Morise M, Felip E, Andric Z, Geater S, et al. Atezolizumab for First-Line treatment of PD-L1-Selected patients with NSCLC. N Engl J Med. 2020;383(14):1328–39.

    Article  CAS  PubMed  Google Scholar 

  26. Gadgeel S, Rodriguez-Abreu D, Speranza G, Esteban E, Felip E, Domine M, Hui R, Hochmair MJ, Clingan P, Powell SF, et al. Updated analysis from KEYNOTE-189: Pembrolizumab or Placebo Plus Pemetrexed and Platinum for previously untreated metastatic nonsquamous non-small-cell Lung Cancer. J Clin Oncol. 2020;38(14):1505–17.

    Article  CAS  PubMed  Google Scholar 

  27. Rittmeyer A, Barlesi F, Waterkamp D, Park K, Ciardiello F, von Pawel J, Gadgeel SM, Hida T, Kowalski DM, Dols MC, et al. Atezolizumab versus Docetaxel in patients with previously treated non-small-cell lung cancer (OAK): a phase 3, open-label, multicentre randomised controlled trial. Lancet. 2017;389(10066):255–65.

    Article  PubMed  Google Scholar 

  28. Herbst RS, Baas P, Kim DW, Felip E, Perez-Gracia JL, Han JY, Molina J, Kim JH, Arvis CD, Ahn MJ, et al. Pembrolizumab versus Docetaxel for previously treated, PD-L1-positive, advanced non-small-cell lung cancer (KEYNOTE-010): a randomised controlled trial. Lancet. 2016;387(10027):1540–50.

    Article  CAS  PubMed  Google Scholar 

  29. Borghaei H, Paz-Ares L, Horn L, Spigel DR, Steins M, Ready NE, Chow LQ, Vokes EE, Felip E, Holgado E, et al. Nivolumab versus Docetaxel in Advanced Nonsquamous Non-small-cell Lung Cancer. N Engl J Med. 2015;373(17):1627–39.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. McGranahan N, Furness AJ, Rosenthal R, Ramskov S, Lyngaa R, Saini SK, Jamal-Hanjani M, Wilson GA, Birkbak NJ, Hiley CT, et al. Clonal neoantigens elicit T cell immunoreactivity and sensitivity to immune checkpoint blockade. Science. 2016;351(6280):1463–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Shur JD, Doran SJ, Kumar S, Ap Dafydd D, Downey K, O’Connor JPB, Papanikolaou N, Messiou C, Koh DM, Orton MR. Radiomics in Oncology: a practical guide. Radiographics. 2021;41(6):1717–32.

    Article  PubMed  Google Scholar 

  32. Wen Q, Yang Z, Zhu J, Qiu Q, Dai H, Feng A, Xing L. Pretreatment CT-Based Radiomics Signature as a potential imaging Biomarker for Predicting the expression of PD-L1 and CD8 + TILs in ESCC. Onco Targets Ther. 2020;13:12003–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Sun R, Limkin EJ, Vakalopoulou M, Dercle L, Champiat S, Han SR, Verlingue L, Brandao D, Lancia A, Ammari S, et al. A radiomics approach to assess tumour-infiltrating CD8 cells and response to anti-PD-1 or anti-PD-L1 immunotherapy: an imaging biomarker, retrospective multicohort study. Lancet Oncol. 2018;19(9):1180–91.

    Article  CAS  PubMed  Google Scholar 

  34. Tian Y, Komolafe TE, Zheng J, Zhou G, Chen T, Zhou B, Yang X. Assessing PD-L1 expression level via Preoperative MRI in HCC based on integrating Deep Learning and Radiomics features. Diagnostics (Basel). 2021;11(10):1875.

    Article  CAS  PubMed  Google Scholar 

  35. Jiang M, Sun D, Guo Y, Guo Y, Xiao J, Wang L, Yao X. Assessing PD-L1 expression level by Radiomic features from PET/CT in Nonsmall Cell Lung Cancer patients: an initial result. Acad Radiol. 2020;27(2):171–9.

    Article  PubMed  Google Scholar 

  36. Sun Z, Hu S, Ge Y, Wang J, Duan S, Song J, Hu C, Li Y. Radiomics study for predicting the expression of PD-L1 in non-small cell lung cancer based on CT images and clinicopathologic features. J Xray Sci Technol. 2020;28(3):449–59.

    CAS  PubMed  Google Scholar 

  37. Chen H, Li S, Zhang Y, Liu L, Lv X, Yi Y, Ruan G, Ke C, Feng Y. Deep learning-based automatic segmentation of meningioma from multiparametric MRI for preoperative meningioma differentiation using radiomic features: a multicentre study. Eur Radiol. 2022;32(10):7248–59.

    Article  PubMed  Google Scholar 

  38. Wang S, Shi J, Ye Z, Dong D, Yu D, Zhou M, Liu Y, Gevaert O, Wang K, Zhu Y et al. Predicting EGFR mutation status in lung adenocarcinoma on computed tomography image using deep learning. Eur Respir J 2019, 53(3).

  39. Tian P, He B, Mu W, Liu K, Liu L, Zeng H, Liu Y, Jiang L, Zhou P, Huang Z, et al. Assessing PD-L1 expression in non-small cell lung cancer and predicting responses to immune checkpoint inhibitors using deep learning on computed tomography images. Theranostics. 2021;11(5):2098–107.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Zhu Y, Man C, Gong L, Dong D, Yu X, Wang S, Fang M, Wang S, Fang X, Chen X, et al. A deep learning radiomics model for preoperative grading in meningioma. Eur J Radiol. 2019;116:128–34.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

Thanks to all authors for their contribution to this article.

Funding

This study has received funding by the National Natural Science Foundation of China (82172026) (Jingshan Gong).

Author information

Authors and Affiliations

Authors

Contributions

Ting Xu: Formal analysis, Writing–original draft; Xiaowen Liu: Conceptualization; Jingshan Gong: Funding acquisition, Methodology, Writing – review & editing; Yaxi Chen: Data curation; Changsi Jiang: Data curation; Shuxing Wang: Data curation.

Corresponding author

Correspondence to Jingshan Gong.

Ethics declarations

Ethical approval

This study was approved by the Shenzhen People’s Hospital Ethics Committee with waiver of informed consent due to the nature of retrospective study (NO. LL-KY-2021058). All methods were conducted ethically in accordance with the World Medical Association Declaration of Helsinki.

Consent for publication

Not applicable.

Competing interests

The authors of this manuscript declare no relationships with any companies, whose products or services may be related to the subject matter of the article.

Additional information

Publisher’s Note

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

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/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Xu, T., Liu, X., Chen, Y. et al. CT-based deep learning radiomics biomarker for programmed cell death ligand 1 expression in non-small cell lung cancer. BMC Med Imaging 24, 196 (2024). https://doi.org/10.1186/s12880-024-01380-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12880-024-01380-8

Keywords