Skip to main content

ORIGINAL RESEARCH article

Front. Oncol., 18 February 2020
Sec. Cancer Imaging and Image-directed Interventions

CT-Based Radiomic Signature as a Prognostic Factor in Stage IV ALK-Positive Non-small-cell Lung Cancer Treated With TKI Crizotinib: A Proof-of-Concept Study

\nHailin Li,&#x;Hailin Li1,2Rui Zhang&#x;Rui Zhang1Siwen Wang,&#x;Siwen Wang2,3Mengjie Fang,Mengjie Fang2,3Yongbei Zhu,Yongbei Zhu2,4Zhenhua Hu,
Zhenhua Hu2,3*Di Dong,
Di Dong2,3*Jingyun Shi
Jingyun Shi5*Jie Tian,
Jie Tian2,4*
  • 1School of Automation, Harbin University of Science and Technology, Harbin, China
  • 2CAS Key Laboratory of Molecular Imaging, Institute of Automation, Chinese Academy of Sciences, Beijing, China
  • 3School of Artificial Intelligence, University of Chinese Academy of Sciences, Beijing, China
  • 4Beijing Advanced Innovation Center for Big Data-Based Precision Medicine, School of Medicine, Beihang University, Beijing, China
  • 5Department of Radiology, Shanghai Pulmonary Hospital, Tongji University School of Medicine, Shanghai, China

Objectives: To identify a computed tomography (CT)-based radiomic signature for predicting progression-free survival (PFS) in stage IV anaplastic lymphoma kinase (ALK)-positive non-small-cell lung cancer (NSCLC) patients treated with tyrosine kinase inhibitor (TKI) crizotinib.

Materials and Methods: This retrospective proof-of-concept study included a cohort of 63 stage IV ALK-positive NSCLC patients who had received TKI crizotinib therapy for model construction and validation. Another independent cohort including 105 stage IV EGFR-positive NSCLC patients was also used for external validation in EGFR-TKI treatment. We initially extracted 481 quantitative three-dimensional features derived from manually segmented tumor volumes of interest. Pearson's correlation analysis along with the least absolute shrinkage and selection operator (LASSO) penalized Cox proportional hazards regression was successively performed to select critical radiomic features. A CT-based radiomic signature for PFS prediction was obtained using multivariate Cox regression. The performance evaluation of the radiomic signature was conducted using the concordance index (C-index), time-dependent receiver operating characteristic (ROC) analysis, and Kaplan–Meier survival analysis.

Results: A radiomic signature containing three features showed significant prognostic performance for ALK-positive NSCLC patients in both the training cohort (C-index, 0.744; time-dependent AUC, 0.895) and the validation cohort (C-index, 0.717; time-dependent AUC, 0.824). The radiomic signature could significantly risk-stratify ALK-positive NSCLC patients (hazard ratio, 2.181; P < 0.001) and outperformed other prognostic factors. However, no significant association with PFS was captured for the radiomic signature in the EGFR-positive NSCLC cohort (log-rank tests, P = 0.41).

Conclusions: The CT-based radiomic features can capture valuable information regarding the tumor phenotype. The proposed radiomic signature was found to be an effective prognostic factor in stage IV ALK mutated nonsynchronous nodules in NSCLC patients treated with a TKI.

Introduction

Non-small-cell lung cancer (NSCLC) accounts for 85% in lung cancer (1), making it the leading cause of cancer-related mortality across the globe (2). An approximate estimate of 5% NSCLC patients are found to have rearrangements in the anaplastic lymphoma kinase (ALK) gene. Although the incidence of ALK seems relatively low among NSCLC patients, a number of nearly 40,000 ALK-positive cases occur annually worldwide (3), and about half of them come from Asia.

Crizotinib, which is regarded as the first-line treatment for stage IV NSCLC patients with ALK gene rearrangements according to the clinical practice guidelines given by the American Society of Clinical Oncology (ASCO) (4), has a proven therapeutical efficacy. And those who respond to crizotinib may show rapid improvement in symptoms, including cough, dyspnea, and pain (5, 6).

In spite of the potency of crizotinib, however, some patients still suffered from the disease progression, and the median survival time just ranged from 7 months to 1 year (4, 7, 8). Some of the reported mechanisms underlying the acquired resistance to crizotinib in ALK translocated cancers involve secondary mutations, ALK amplification, KIT amplification, and autophosphorylation of epidermal growth factor receptor (EGFR) (913). Therefore, at the genetic level, resistance to tyrosine kinase inhibitor (TKI) is hard to detect (14).

In patients with ALK rearrangements, when disease progresses after crizotinib therapy, other ALK inhibitor (alectinib or brigatinib or ceritinib) or local therapy (e.g., stereotactic ablative radiotherapy or surgery) is recommended by the National Comprehensive Cancer Network (NCCN) guidelines (15). Thus, if the progression on crizotinib can be predicted, other ALK inhibitor or local therapy can be considered in the initial systemic therapy strategy, which will help the individualized treatment decision-making. Therefore, finding a prognostic marker to predict the development of TKI resistance of ALK-positive patients will be of great significance.

Using automatic feature extraction algorithms, radiomics is capable of converting embedded information in medical images into mineable data (16, 17), which has been widely applied in the prediction of preoperative distant metastasis, histologic subtype classification, and so on (1820). Prognosis based on radiomics is gaining popularity as associations between radiomic features and the underlying genomic patterns emerge in various cancers (2124). On one hand, a latest study showed that a subset of radiomic features were able to consistently capture texture information about the underlying tissue histology, but some of them were incapable to be observed at the purely human level (25, 26). On the other hand, computed tomography (CT) imaging modality can quantify the tissue density and is widely used in radiology researches.

Gene mutations in adenocarcinoma include EGFR mutation, kirsten rat sarcoma viral oncogene (KRAS) mutation, ALK mutation or c-ros oncogene 1 receptor kinase (ROS1) rearrangement, etc. EGFR mutations occupy the majority of gene mutations in adenocarcinoma. Previous studies ever used radiomics to predict progression-free survival (PFS) in stage IV EGFR-mutant NSCLC patients with EGFR-TKI therapy (24). However, as another mutation independent of EGFR, no studies have yet been conducted to investigate the potential of radiomics in the prediction of ALK inhibitor-related PFS in ALK rearranged patients with stage IV NSCLC.

Consequently, this proof-of-concept study aimed to launch and verify a radiomic signature to predict PFS in stage IV ALK-positive NSCLC patients treated with TKI crizotinib, and further investigate the association of this radiomic signature with PFS in EGFR-mutant NSCLC patients with EGFR-TKI therapy.

Materials and Methods

Patients

From January 2012 to January 2017, a total of 63 stage IV ALK-positive NSCLC patients (33 males and 30 females) and 105 stage IV EGFR-positive NSCLC patients (42 males and 63 females) were enrolled. This study was examined and approved by the review committee as a retrospective research. Confidential information for each patient has been hidden. The inclusion criteria for ALK-positive patients were as follows: (1) histopathologically confirmed stage IV NSCLC according to the TNM classification system of the American Joint Committee on Cancer (AJCC); (2) available test results for ALK mutation status; and (3) underwent contrast-enhanced CT scans 2 weeks before the crizotinib therapy. Patients were ruled out for any of the following reasons: (1) underwent other antitumor therapies such as resection for local advanced or metastatic disease; and (2) death for other reasons during follow-up visits. A rapid prescreening for ALK mutation was done by immunohistochemistry (IHC). Positive findings were confirmed subsequently by fluorescence in situ hybridization (FISH) analysis (27, 28). The inclusion criteria for EGFR-positive patients were as follows: (1) histopathologically confirmed stage IV NSCLC according to the TNM classification system of the AJCC; (2) available data for the EGFR mutation status; and (3) underwent contrast-enhanced CT scans 2 weeks before the EGFR-TKI treatment. Patients were ruled out for any of the following reasons: (1) underwent other antitumor therapies such as resection for locally advanced or metastatic disease; and (2) death for other reasons during follow-up visits. Note that patients with synchronous nodules were not included in the study.

PFS was counted from the start of treatment with TKI to the confirmation of disease progression or death. Basic clinical characteristics including sex, age, and smoking status were recorded. The clinical characteristics and follow-up information for patients are illustrated in Table 1.

TABLE 1
www.frontiersin.org

Table 1. Baseline demographic and clinicopathologic characteristics of patients in the training and validation cohorts.

CT Image Acquisition and Tumor Segmentation

CT scans were examined after a 60-s delay following an intravenous injection of 100 ml Iopromide (Ultravist-300; Bayer Schering Pharma, Berlin Germany) at a rate of 3 ml/s for enhancement. The scanning parameters were 120 kV, 160 mAs, 0.6 s rotation time, and a matrix of 512×512. Portal venous phase CT images were reconstructed using a slice thickness of 1 mm.

A radiologist engaged in the chest CT interpretation for more than 7 years segmented the primary tumors of all the patients slice by slice. Three-dimensional (3D) manual segmentation was performed using the ITK-SNAP software (www.itksnap.org, version 3.6.0).

Statistical Analysis

The designed workflow of the study is shown in Figure 1. All statistical tests were conducted using the R software (www.r-project.org, version 3.5.1). A summary of R packages adopted in the study can be found in Supplementary Appendix 1. To ensure that there were enough cases left in the validation cohort for model evaluation, we initially separated the dataset into the training cohort and the validation cohort at a 1:1 ratio. Besides, we further performed 10 random cohort allocations and conducted feature selections in the training cohorts to test the stability of the radiomic features to be selected. EGFR-positive NSCLC patients were used as another independent validation cohort. Features were initially normalized. Radiomic feature selection and radiomic signature construction were carried out based on the training cohort, whereas the validation cohorts were only for performance verification.

FIGURE 1
www.frontiersin.org

Figure 1. Study design. (A) Tumor segmentation. Shown are typical CT images of lung cancer patients with tumor contours and three-dimensional visualizations. (B) Radiomic feature extraction. Four types of radiomic features were extracted from VOIs. (C) Feature selection and statistical analysis. LASSO penalized Cox proportional hazards regression was adopted to select critical features. CT, computed tomography; VOI, volume of interest; LASSO, least absolute shrinkage and selection operator.

Radiomic Feature Extraction

Within each tumor volume of interest (VOI), 481 radiomic features across four feature types were extracted from 3D CT images. These characteristics embodied first-order, shape, texture, and wavelet features. Feature extraction procedures were conducted within the VOIs with in-house algorithms implemented in Matlab R2018b (The MathWorks, Natick, MA, USA).

First-order features characterized the distribution of intensity values, according to commonly used and basic metrics such as energy, entropy, kurtosis, etc. 3D shape descriptors were included to describe the shape of the tumor VOI, which was independent of the gray-level intensity distribution. Volume, surface area, sphericity, compactness, and other features were embodied. Gray-level co-occurrence matrix (GLCM) (29) described the second-order joint probability function of an image region, and gray-level run-length matrix (GLRLM) (30, 31) calculated texture features that represented voxel intensities of spatial distribution or patterns. Wavelet features effectively decoupled information via decomposing the original image. Details about feature extraction algorithms and parameters are given in Supplementary Appendix 2.

Feature Discovery and Radiomic Signature Building

Based on Harrell's guideline, the number of features must be reduced to <1/10 of the events in the multivariate Cox regression (32). Therefore, Pearson correlation analysis and the least absolute shrinkage and selection operator (LASSO) penalized Cox proportional hazards regression were successively performed to reduce the number of radiomic features and select critical ones. In the training cohort, the Pearson correlation coefficient matrix of all radiomic features was calculated. An absolute value of correlation coefficient >0.8 indicated a strong correlation between the two features, one of which was removed then. The LASSO penalized Cox proportional hazards regression was adopted to select critical features with nonzero weights. According to minimizing the partial likelihood deviance, an optimal value of λ was adopted via the 10-fold cross-validation in the training cohort. Finally, the radiomic signature was built by fitting the selected critical features using multivariate Cox regression.

Validation of Radiomic Signature

Associations between the radiomic signature and PFS were evaluated in both ALK-mutant cohorts. Based on the median radiomic signature value generated in the training cohort, patients were classified into high and low radiomic signature value groups, which were expected to represent rapid and slow progression treated with crizotinib, respectively. When not all the subjects continue to be followed up in the study, Kaplan–Meier survival curve is usually one of the best choices to be used to measure the fraction of subjects living for a certain amount of time after treatment. At the same time, the log-rank test is always used to obtain the statistical significance level in survival difference between two groups. In this study, the potential association between the radiomic signature and PFS was assessed by Kaplan–Meier survival curves as well as log-rank tests (G-rho rank test, rho = 1). The discriminative ability of the radiomic signature was assessed by time-dependent receiver operating characteristic (ROC) analysis and the Harrell's concordance index (C-index). Due to the relatively short survival time range of stage IV NSCLC patients in this study, time-dependent ROCs were plotted in both cohorts at 6 months. Furthermore, to identify whether our model had the prognostic value in other treatments, we constructed an EGFR-mutant cohort to further validate the radiomic signature by Kaplan–Meier survival curves along with log-rank tests.

Performance Assessment and Comparison of Prognostic Factors

To further verify the prognostic performance of the radiomic signature and test the prognostic ability of clinical characteristics, we established a clinical model. Aiming to investigate whether the combination of clinical characteristics and radiomic signature could further improve the model performance, we also established a combined model via multivariate Cox proportional hazards regression. C-index statistics were calculated to assess the model performance. Herein, the outputs of the two models could be deemed two prognostic factors for stage IV ALK-positive NSCLC patients treated with TKI crizotinib.

Moreover, to find the optimal prognostic factor, the association of each prognostic factor with PFS, including each clinical characteristic, the radiomic signature, the clinical model (constructed using significant clinical characteristics), and the combined model, was analyzed and compared with respect to the hazard ratio (HR) and corresponding P-values in the validation cohort, respectively. P-values are computed using the likelihood ratio test to measure whether the HR is different from 1. The forest plots of HR were conducted.

Results

Clinical Characteristics and PFS

In the ALK-mutant cohorts, a total of 63 patients who met the inclusion and exclusion requirements were enrolled, all of whom were stage IV ALK-mutant NSCLC treated with TKI crizotinib. The last follow-up was on April 30, 2016. Throughout the follow-up period, 43 patients experienced progression (21 patients in the training cohort, 22 patients in the validation cohort), and 20 patients were lost to follow-up (11 patients in the training cohort, 9 patients in the validation cohort). The PFS was balanced between the training cohort (median PFS, 248 days) and the validation cohort (median PFS, 265 days). Furthermore, there were no differences between the two cohorts with respect to baseline clinical characteristics including age, sex, and smoking status by chi-square tests (Table 1).

In the EGFR-mutant cohort, we included 105 patients who met the enrollment requirements, all of whom were stage IV EGFR-mutant NSCLC treated with TKI treatment. Throughout the follow-up period, the median PFS was 238 days.

Radiomic Signature Building

In the training cohort, 34 features were retained after Pearson's correlation analysis. Then, cluster shade and inverse variance calculated from GLCM were selected by the LASSO Cox regression. The value of parameter λ was 0.1876 (Figure 2). The three selected radiomic features stayed stable in the training cohorts from 10 random cohort allocations (Supplementary Appendix 3). Finally, the formula for the radiomic signature generated by multivariate Cox regression was shown below:

Radiomic Signature = 0.312×first order minimum×cluster                                               shade+0.602×short run high gray                                               level emphasis (SRHGLE)
FIGURE 2
www.frontiersin.org

Figure 2. Feature selection using the LASSO Cox model. (A) In the LASSO Cox model, the minimum standard is adopted to obtain the value of the super parameter λ by 10-fold cross-validation. The λ value was confirmed as 0.1876. (B) Shown here is a coefficient sectional view plotted against the log (λ) magnitude. Based on 10-fold cross-validation, the optimal λ corresponding to three nonzero coefficients were obtained where the vertical line was drawn. LASSO, least absolute shrinkage and selection operator.

Validation of Radiomic Signature

In the ALK-mutant patients, the radiomic signature was relevant to PFS in the training cohort [HR, 2.718; 95% confidence interval (CI): 1.599–4.621; P < 0.01] and the validation cohort (HR, 2.181; 95% CI: 1.370–3.471; P < 0.01). The C-index yielded from the radiomic signature was 0.744 (95% CI: 0.678–0.809) in the training cohort and 0.717 (95% CI: 0.614–0.821) in the validation cohort. Moreover, in the 10 random validation cohorts, the mean of the C-index was 0.709 (range: 0.648–0.759). Time-dependent ROC curves achieved a training area under the curve (AUC) of 0.895 at 6 months and a validation AUC of 0.824 (Figure 3).

FIGURE 3
www.frontiersin.org

Figure 3. Time-dependent ROC curves of the radiomic signature in the training cohort (A) and validation cohort (B) at 6 months. AUCs were used to assess the prognostic accuracy in both cohorts. ROC, receiver operating characteristics; AUC, area under the curve.

Based on the cutoff value of the median radiomic signature in the training cohort, we allocated patients into high and low signature value groups. Thereinto, the high signature value group (median survival: training cohort, 107 days; validation cohort, 208 days) had a shorter median PFS than the low signature value group (median survival: training cohort, 377 days; validation cohort, 279 days). In both cohorts, the PFS differed significantly between the high and low signature value groups based on Kaplan–Meier survival curves (Figure 4) with log-rank P-values of 0.00035 and 0.0085, respectively.

FIGURE 4
www.frontiersin.org

Figure 4. Kaplan–Meier survival curves of ALK-positive NSCLC patients. The P-values were calculated using log-rank tests. Kaplan–Meier survival analysis for patients stratified by median radiomic signature in the training (A) and validation (B) cohorts showed a significant association between the radiomic signature and PFS. ALK, anaplastic lymphoma kinase; NSCLC, non-small-cell lung cancer; PFS, progression-free survival.

In the EGFR-mutant cohort, however, the radiomic signature failed to significantly risk-stratify NSCLC patients by Kaplan–Meier survival curves with a log-rank P-value of 0.41 (Figure 5).

FIGURE 5
www.frontiersin.org

Figure 5. Kaplan–Meier survival curve of the radiomic signature for EGFR-positive NSCLC patients. The P-values were calculated using log-rank test. EGFR, epidermal growth factor receptor; NSCLC, non-small-cell lung cancer; PFS, progression-free survival.

Performance Analysis and Comparison of Prognostic Factors

Since the smoking status and age were found to be significantly correlated with PFS in the training cohort (P = 0.045 and P = 0.028, respectively), a clinical model was built by the two characteristics. The C-index yielded from the clinical model was 0.638 (95% CI: 0.498–0.778) in the training cohort and 0.525 (95% CI: 0.416–0.634) in the validation cohort. A trial was made to feed the smoking status, age, and the radiomic signature into a multivariate Cox regression, the result of which showed that adding the clinical characteristics did not result in an extra benefit compared with the radiomic signature alone, with a C-index of 0.789 (95% CI: 0.694–0.885) in the training cohort and 0.655 (95% CI: 0.614–0.821) in the validation cohort. We also conducted the feature selection on both radiomic features and clinical variables; only the radiomic features were reserved in the model.

The forest plots showed that the radiomic signature had the significant correlation with PFS in the validation cohort (P < 0.001). The combined model, which was built by the clinical characteristics and radiomic signature, also showed a significant correlation with PFS, but it had a much higher P-value than the radiomic signature alone (Figure 6). Thereinto, the radiomic signature also demonstrated superior prognostic efficacy to the clinical characteristics (sex, age, and smoking status), the clinical model (constructed using smoking status and age), and the combined model.

FIGURE 6
www.frontiersin.org

Figure 6. The forest plots of HR for each clinical characteristic, each selected radiomic feature, the radiomic signature, the clinical model, and the combined model in the validation cohort. The P-value was computed using the likelihood ratio test, and a *P < 0.05 was considered statistically significant. HR, hazard ratio.

Discussion

A CT-based radiomic signature was identified as an independent prognostic factor to predict PFS in stage IV ALK-positive NSCLC patients treated with TKI crizotinib. The radiomic signature managed to risk-stratify the patients into groups of rapid progression and slow progression. Our results showed no significant correlation between PFS and the radiomic signature in EGFR-mutant NSCLC patients, indicating the specific stratification ability of the radiomic signature in predicting ALK-TKI resistance. In addition, adding the clinical characteristics did not result in an extra benefit compared with the radiomic signature alone. The clinical model built by clinical characteristics failed to predict the prognosis as well, which demonstrated that the three clinical characteristics (sex, age, and smoking status) could not predict the prognosis in ALK-positive NSCLC treated with TKI crizotinib.

To our knowledge, the correlations between CT-based radiomic features and PFS in ALK-positive NSCLC patients treated with crizotinib have not been evaluated. The patients with a shorter PFS would particularly benefit if we could predict their response to crizotinib therapy and offer prompt guidance for alternative treatment strategies. For stage IV patients with rapidly progressed tumors, a platinum-based regimen is helpful (33). There is evidence showing that early palliative care together with standard care improved the overall mood, quality of life, and even survival in NSCLC patients, although these patients received less aggressive treatments compared to patients receiving standard care alone (34). Several ongoing studies are evaluating treatment options for ALK-positive patients. Crizotinib therapy can be continued for ALK-positive patients showing disease progression as long as they do not have multiple systemic symptomatic lesions (35). Ceritinib was suggested as a follow-up treatment by NCCN Panel for ALK-positive NSCLC patients who progressed after crizotinib therapy. Patients with ALK-positive metastatic NSCLC who have progressed or cannot tolerate crizotinib can also be treated with an oral TKI alectinib, which is approved by the FDA (36). Our novel radiomic signature showed potential in identifying patients with high and low risk of disease progression, and hence could guide the course of treatment. For patients who are prone to crizotinib resistance, it can be compared between the curative effects of targeted agents after chemotherapy and the first-line ALK inhibitor. Besides, although a new generation of TKIs are being promoted, the size of clinical data available for research is still very small. More researches on other TKI resistance may be inspired by the outcomes of this study in the future.

Radiomics, a novel and noninvasive method, can capture information from the entire tumor quantitatively and remove the obstacles hindering the prediction of PFS in patients with the same stage disease. In the process of building a radiomic signature, 481 features were reduced to two components using Pearson's correlation analysis and the LASSO Cox regression. Pearson's correlation analysis could greatly reduce the redundancy of the data by eliminating highly correlated features. The LASSO method is ideal for analyzing cases when the sample size is smaller than the feature dimension, and it generally gives accurate results and avoids overfitting (37, 38). Thus, the critical features obtained were combined as a radiomic signature, revealing sufficient discrimination in both cohorts. The selected radiomic features including cluster shade and SRHGLE are both texture features. Cluster shade measures the skewness and the uniformity of the GLCM. The higher the cluster shade, the greater the asymmetry about the mean. SRHGLE represents the measurement of the joint distribution of higher gray-level shorter run lengths of the GLRLM. They may reflect tumor heterogeneity in two ways, and a higher value of cluster shade or SRHGLE may indicate higher tumor heterogeneity. The interpretations between the two texture features and tumor heterogeneity need further investigation. Besides, in terms of the multivariate Cox regression formula, more rapid resistance to crizotinib is more likely to be caused by tumors with high heterogeneity.

This study still has some limitations, given the fact that the study is retrospective. Since the sample size was small and only Asian data sets were included, this proof-of-concept study needs further verification on large-scale cohorts. We excluded patients with synchronous nodules, which are commonly seen among NSCLC patients; further study on synchronous nodules should be performed in the future. In the EGFR cohort, not only one kind of crizotinib was used, the influence of which should be further studied. In addition, further prospective trials are in desperate need to be carried out to assure the predictive efficacy of the CT-based radiomic signature in TKI resistance.

In conclusion, this study discussed the correlation of a radiomic signature with PFS and its predictive potential in stage IV ALK-positive NSCLC patients treated with crizotinib. Besides, it can predict resistance to the ALK inhibitor, thereby enabling the effective prognostic prediction and stratification, which can improve the individualized treatment with ALK inhibitors.

Data Availability Statement

The datasets for this article are not publicly available because of patient information privacy. Requests to access the datasets should be directed to Prof. Jie Tian, amllLnRpYW4mI3gwMDA0MDtpYS5hYy5jbg==.

Ethics Statement

The studies involving human participants were reviewed and approved by the review committee of Shanghai Pulmonary Hospital. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author Contributions

JT, JS, and DD conceived and launched this study. HL, SW, and RZ designed the medical and statistical analyses. JS, ZH, and DD collected cases and implemented the control of image quality and clinical diagnosis. HL, MF, and YZ analyzed the data and carried out statistical experiments. HL, SW, MF, and JS provided result interpretation. HL, SW, and RZ wrote the first draft of this manuscript. DD, ZH, JT, and JS revised and edited the final version. All the authors reviewed and approved the manuscript.

Funding

This work was supported by the National Key R&D Program of China (2017YFC1308700, 2017YFA0205200, 2017YFC1309100, 2016YFC0102600, 2016YFA0100900, 2016YFA0100902), the National Natural Science Foundation of China (91959130, 81971776, 81771924, 81501616, 61622117, 81671759, 81227901, 81671851, 81527805), the Beijing Natural Science Foundation (L182061, JQ19027), the Bureau of International Cooperation of Chinese Academy of Sciences (173211KYSB20160053), Beijing Nova Program (Z181100006218046), the Scientific Instrument Developing Project of the Chinese Academy of Sciences (YZ201502, YZ201672), and the Youth Innovation Promotion Association CAS (2017175).

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.

Supplementary Material

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

Abbreviations

CT, computed tomography; PFS, progression-free survival; ALK, anaplastic lymphoma kinase; NSCLC, non-small-cell lung cancer; TKI, tyrosine kinase inhibitor; LASSO, least absolute shrinkage and selection operator; C-index, concordance index; ASCO, American Society of Clinical Oncology; EGFR, epidermal growth factor receptor; NCCN, National Comprehensive Cancer Network; ICC, interclass correlation coefficient; VOI, volume of interest; GLCM, gray-level co-occurrence matrix; GLRLM, gray-level run-length matrix; ROC, receiver operating characteristic; HR, hazard ratio; CI, confidence interval; AUC, area under the curve.

References

1. Molina JR, Yang P, Cassivi SD, Schild SE, Adjei AA, editors. Non-small cell lung cancer: epidemiology, risk factors, treatment, and survivorship. Mayo Clin Proc. (2008) 83:584–94. doi: 10.4065/83.5.584

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. (2018) 68:394–424. doi: 10.3322/caac.21492

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Shaw AT, Engelman JA. ALK in lung cancer: past, present, and future. J Clin Oncol. (2013) 31:1105–11. doi: 10.1200/JCO.2012.44.5353

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Shaw AT, Gandhi L, Gadgeel S, Riely GJ, Cetnar J, West H, et al. Alectinib in ALK-positive, crizotinib-resistant, non-small-cell lung cancer: a single-group, multicentre, phase 2 trial. Lancet Oncol. (2016) 17:234–42. doi: 10.1016/S1470-2045(15)00488-X

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Hallberg B, Palmer RH. Mechanistic insight into ALK receptor tyrosine kinase in human cancer biology. Nat Rev Cancer. (2013) 13:685. doi: 10.1038/nrc3580

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Solomon BJ, Mok T, Kim D-W, Wu Y-L, Nakagawa K, Mekhail T, et al. First-Line Crizotinib versus chemotherapy in ALK-positive lung cancer. N Engl J Med. (2014) 371:2167–77. doi: 10.1056/NEJMoa1408440

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Bang Y-J. Treatment of ALK-positive non–small cell lung cancer. Arch Pathol Lab Med. (2012) 136:1201–4. doi: 10.5858/arpa.2012-0246-RA

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Choi YL, Soda M, Yamashita Y, Ueno T, Takashima J, Nakajima T, et al. EML4-ALK mutations in lung cancer that confer resistance to ALK inhibitors. N Engl J Med. (2010) 363:1734–9. doi: 10.1056/NEJMoa1007478

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Friboulet L, Li N, Katayama R, Lee CC, Gainor JF, Crystal AS, et al. The ALK inhibitor ceritinib overcomes crizotinib resistance in non-small cell lung cancer. Cancer Discov. (2014) 4:662–73. doi: 10.1158/2159-8290.CD-13-0846

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Katayama R, Khan TM, Benes C, Lifshits E, Ebi H, Rivera VM, et al. Therapeutic strategies to overcome crizotinib resistance in non-small cell lung cancers harboring the fusion oncogene EML4-ALK. Proc Natl Acad Sci USA. (2011) 108:7535–40. doi: 10.1073/pnas.1019559108

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Sasaki T, Okuda K, Zheng W, Butrynski J, Capelletti M, Wang L, et al. The neuroblastoma-associated F1174L ALK mutation causes resistance to an ALK kinase inhibitor in ALK-translocated cancers. Cancer Res. (2010) 70:10038–43. doi: 10.1158/0008-5472.CAN-10-2956

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Doebele RC, Pilling AB, Aisner DL, Kutateladze TG, Le AT, Weickhardt AJ, et al. Mechanisms of resistance to crizotinib in patients with ALK gene rearranged non-small cell lung cancer. Clin Cancer Res. (2012) 18:1472–82. doi: 10.1158/1078-0432.CCR-11-2906

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Katayama R, Shaw AT, Khan TM, Mino-Kenudson M, Solomon BJ, Halmos B, et al. Mechanisms of acquired crizotinib resistance in ALK-rearranged lung cancers. Sci Transl Med. (2012) 4:120ra17. doi: 10.1126/scitranslmed.3003316

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Qian M, Zhu B, Wang X, Liebman M. Drug resistance in ALK-positive non-small cell lung cancer patients. Semin Cell Dev Biol. (2017) 64:150–7. doi: 10.1016/j.semcdb.2016.09.016

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Demetri GD, Benjamin RS, Blanke CD, Blay J-Y, Casali P, Choi H, et al. NCCN Task Force report: management of patients with gastrointestinal stromal tumor (GIST)–update of the NCCN clinical practice guidelines. J Natl Comprehen Cancer Netw. (2007) 5:S1–29; quiz S30. doi: 10.6004/jnccn.2010.0116

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Aerts HJ, Velazquez ER, Leijenaar RT, Parmar C, Grossmann P, Carvalho S, et al. Decoding tumour phenotype by noninvasive imaging using a quantitative radiomics approach. Nat Commun. (2014) 5:4006. doi: 10.1038/ncomms5644

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Gillies RJ, Kinahan PE, Hricak H. Radiomics: images are more than pictures, they are data. Radiology. (2015) 278:563–77. doi: 10.1148/radiol.2015151169

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Dong D, Tang L, Li Z-Y, Fang M-J, Gao J-B, Shan X-H, et al. Development and validation of an individualized nomogram to identify occult peritoneal metastasis in patients with advanced gastric cancer. Ann Oncol. (2019) 30:431–8. doi: 10.1093/annonc/mdz001

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Zhu X, Dong D, Chen Z, Fang M, Zhang L, Song J, et al. Radiomic signature as a diagnostic factor for histologic subtype classification of non-small cell lung cancer. Eur Radiol. (2018) 28:2772–8. doi: 10.1007/s00330-017-5221-1

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Bi WL, Hosny A, Schabath MB, Giger ML, Birkbak NJ, Mehrtash A, et al. Artificial intelligence in cancer imaging: clinical challenges and applications. CA Cancer J Clin. (2019) 69:127–57. doi: 10.3322/caac.21552

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Pickles MD, Lowry M, Gibbs P. Pretreatment prognostic value of dynamic contrast-enhanced magnetic resonance imaging vascular, texture, shape, and size parameters compared with traditional survival indicators obtained from locally advanced breast cancer patients. Investig Radiol. (2016) 51:177–85. doi: 10.1097/RLI.0000000000000222

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Leijenaar RT, Carvalho S, Hoebers FJ, Aerts HJ, van Elmpt WJ, Huang SH, et al. External validation of a prognostic CT-based radiomic signature in oropharyngeal squamous cell carcinoma. Acta Oncol. (2015) 54:1423–9. doi: 10.3109/0284186X.2015.1061214

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Wang S, Shi J, Ye Z, Dong D, Yu D, Zhou M, et al. Predicting EGFR mutation status in lung adenocarcinoma on computed tomography image using deep learning. Eur Respir J. (2019) 53:1800986. doi: 10.1183/13993003.00986-2018

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Song J, Shi J, Dong D, Fang M, Zhong W, Wang K, et al. A new approach to predict progression-free survival in stage IV EGFR-mutant NSCLC patients with EGFR-TKI therapy. Clin Cancer Res. (2018) 24:3583–92. doi: 10.1158/1078-0432.CCR-17-2507

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Bobholz S, Lowman A, Barrington A, Brehler M, McGarry S, Cochran EJ, et al. Radiomic features of multi-parametric MRI present stable associations with analogous histological features in brain cancer patients. arXiv[Preprint].arXiv:191207567. (2019).

Google Scholar

26. Stoyanova R, Chinea F, Kwon D, Reis IM, Tschudi Y, Parra NA, et al. An automated multiparametric MRI quantitative imaging prostate habitat risk scoring system for defining external beam radiation therapy boost volumes. Int J Radiat Oncol Biol Phys. (2018) 102:821–9. doi: 10.1016/j.ijrobp.2018.06.003

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Lindeman NI, Cagle PT, Beasley MB, Chitale DA, Dacic S, Giaccone G, et al. Molecular testing guideline for selection of lung cancer patients for EGFR and ALK tyrosine kinase inhibitors: guideline from the College of American Pathologists, International Association for the Study of Lung Cancer, and Association for Molecular Pathology. J Mol Diagn. (2013) 8:823–59. doi: 10.1097/JTO.0b013e318290868f

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Mino-Kenudson M, Chirieac LR, Law K, Hornick JL, Lindeman N, Mark EJ, et al. A novel, highly sensitive antibody allows for the routine detection of ALK-rearranged lung adenocarcinomas by standard immunohistochemistry. Clin Cancer Res. (2010) 16:1561–71. doi: 10.1158/1078-0432.CCR-09-2845

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Haralick RM, Shanmugam K. Textural features for image classification. IEEE Trans Syst Man Cybernet. (1973) 3:610–21. doi: 10.1109/TSMC.1973.4309314

CrossRef Full Text | Google Scholar

30. Galloway MM. Texture analysis using gray level run lengths. Comput Graph Image Proc. (1975) 4:172–9. doi: 10.1016/S0146-664X(75)80008-6

CrossRef Full Text

31. Dasarathy BV, Holder EB. Image characterizations based on joint gray level—run length distributions. Pattern Recogn Lett. (1991) 12:497–502. doi: 10.1016/0167-8655(91)80014-2

CrossRef Full Text | Google Scholar

32. Harrell FE Jr. Regression Modeling Strategies: With Applications to Linear Models, Logistic and Ordinal Regression, and Survival Analysis. Cham: Springer International Publishing (2015).

Google Scholar

33. Azzoli CG, Baker S Jr, Temin S, Pao W, Aliff T, Brahmer J, et al. American Society of Clinical Oncology Clinical Practice Guideline update on chemotherapy for stage IV non-small-cell lung cancer. J Clin Oncol. (2009) 27:6251–66. doi: 10.1200/JCO.2009.23.5622

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Ettinger DS, Wood DE, Aisner DL, Akerley W, Bauman J, Chirieac LR, et al. Non–small cell lung cancer, version 5.2017, NCCN clinical practice guidelines in oncology. J Natl Comprehen Cancer Netw. (2017) 15:504–35. doi: 10.6004/jnccn.2017.0050

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Shaw AT, Kim D-W, Nakagawa K, Seto T, Crinó L, Ahn M-J, et al. Crizotinib versus chemotherapy in advanced ALK-positive lung cancer. N Engl J Med. (2013) 368:2385–94. doi: 10.1056/NEJMoa1214886

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Larkins E, Blumenthal GM, Chen H, He K, Agarwal R, Gieser G, et al. FDA approval: alectinib for the treatment of metastatic, ALK-positive non-small cell lung cancer following crizotinib. Clin Cancer Res. (2016) 22:5171–6. doi: 10.1158/1078-0432.CCR-16-1293

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Gui J, Li H. Penalized cox regression analysis in the high-dimensional and low-sample size settings, with applications to microarray gene expression data. Bioinformatics. (2005) 21:3001–8. doi: 10.1093/bioinformatics/bti422

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Hepp T, Schmid M, Gefeller O, Waldmann E, Mayr A. Approaches to regularized regression – a comparison between gradient boosting and the Lasso. Methods Inf Med. (2016) 55:422–30. doi: 10.3414/ME16-01-0033

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: computed tomography, radiomics, non-small-cell lung cancer, tyrosine kinase inhibitor resistance, anaplastic lymphoma kinase

Citation: Li H, Zhang R, Wang S, Fang M, Zhu Y, Hu Z, Dong D, Shi J and Tian J (2020) CT-Based Radiomic Signature as a Prognostic Factor in Stage IV ALK-Positive Non-small-cell Lung Cancer Treated With TKI Crizotinib: A Proof-of-Concept Study. Front. Oncol. 10:57. doi: 10.3389/fonc.2020.00057

Received: 11 October 2019; Accepted: 13 January 2020;
Published: 18 February 2020.

Edited by:

Peng Huang, Johns Hopkins University, United States

Reviewed by:

Zhongxiang Ding, Hangzhou First People's Hospital, China
Seyedmehdi Payabvash, Yale University, United States
Jun Dang, First Affiliated Hospital of Chongqing Medical University, China

Copyright © 2020 Li, Zhang, Wang, Fang, Zhu, Hu, Dong, Shi and Tian. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Zhenhua Hu, emhlbmh1YS5odSYjeDAwMDQwO2lhLmFjLmNu; Di Dong, ZGkuZG9uZyYjeDAwMDQwO2lhLmFjLmNu; Jingyun Shi, c2hpamluZ3l1bjg5MTc5JiN4MDAwNDA7MTI2LmNvbQ==; Jie Tian, amllLnRpYW4mI3gwMDA0MDtpYS5hYy5jbg==

These authors have contributed equally to this work and share first authorship

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.