- 1Department of Imaging, Nanfang Hospital, Southern Medical University, Guangzhou, China
- 2PET/CT Center, Department of Nuclear Medicine, First Affiliated Hospital of Shenzhen University, Shenzhen Second People’s Hospital, Shenzhen, China
- 3College of Information Engineering, Shenzhen University, Shenzhen, China
Background: Many patients experience recurrence of renal cell carcinoma (RCC) after radical and partial nephrectomy. Radiomics nomogram is a newly used noninvasive tool that could predict tumor phenotypes.
Objective: To investigate Radiomics Features (RFs) associated with progression-free survival (PFS) of RCC, assessing its incremental value over clinical factors, and to develop a visual nomogram in order to provide reference for individualized treatment.
Methods: The RFs and clinicopathological data of 175 patients (125 in the training set and 50 in the validation set) with clear cell RCC (ccRCC) were retrospectively analyzed. In the training set, RFs were extracted from multiphase enhanced CT tumor volume and selected using the stability LASSO feature selection algorithm. A radiomics nomogram final model was developed that incorporated the RFs weighted sum and selected clinical predictors based on the multivariate Cox proportional hazard regression. The performances of a clinical variables-only model, RFs-only model, and the final model were compared by receiver operator characteristic (ROC) analysis and DeLong test. Nomogram performance was determined and validated with respect to its discrimination, calibration, reclassification, and clinical usefulness.
Results: The radiomics nomogram included age, clinical stage, KPS score, and RFs weighted sum, which consisted of 6 selected RFs. The final model showed good discrimination, with a C-index of 0.836 and 0.706 in training and validation, and good calibration. In the training set, the C-index of the final model was significantly larger than the clinical-only model (DeLong test, p = 0.008). From the clinical variables-only model to the final model, the reclassification of net reclassification improvement was 18.03%, and the integrated discrimination improvement was 19.08%. Decision curve analysis demonstrated the clinical usefulness of the radiomics nomogram.
Conclusion: The CT-based RF is an improvement factor for clinical variables-only model. The radiomics nomogram provides individualized risk assessment of postoperative PFS for patients with RCC.
Introduction
Renal cell carcinoma (RCC) is a malignant tumor originating from the proximal tubular epithelial system of renal parenchyma, and accounts for about 85% of all adult renal malignant tumors. It is estimated that, in 2020, there were 431,288 new cases of RCC worldwide, resulting in 179,368 deaths, and accounting for 2.2% and 1.8% of global new cancer morbidity and mortality, respectively. In addition, the incidence of renal cancer is increasing yearly (1). The clear cell renal cell carcinoma (ccRCC) is the most common subtype accounting for about 75% of all RCC (2), and is associated with high invasion and poor prognosis (3, 4). According to the AJCC Tumor Classification Criteria eighth Edition (2017) (5), surgery is the preferred treatment for patients with stage I–III RCC, and is associated with a 5-year survival rate of 71% to 91% (6). However, approximately 20% to 30% of patients will relapse after surgery (7). If we can predict these patients with high risk of recurrence before surgery, and give them targeted treatment and close follow-up, it will be very helpful to improve the prognosis of these patients.
Traditional radiotherapy and chemotherapy have poor efficacy for RCC, and there are no effective adjuvant therapies for RCC. A recent clinical trial showed that a subset of patients with more aggressive disease could benefit from targeted therapy after surgery (8). According to National Comprehensive Cancer Network and European Association of Urology guidelines, adjuvant therapy can reduce the recurrence rate of stage III ccRCC, which is associated with a high recurrence risk. However, about 50% of patients in this high-risk subgroup still do not have postoperative recurrence, and do not need to receive expensive adjuvant targeted therapy. Therefore, there is a need to develop prognostic factors to identify patients who will and will not benefit from adjuvant targeted therapy.
There are currently no markers for a diagnosis of RCC. Tumor stage and pathological nuclear grade are the most important prognostic factors. Nevertheless, distinct outcomes are demonstrated in patients with equivalent tumor-node-metastasis stage and pathological grade; they cannot fully address the issue of individualized treatment for patients with different recurrence risks (9, 10).
The field of artificial intelligence and radiomics has developed rapidly in recent years. Many studies have demonstrated that radiomics can be used to assess the heterogeneity of tumors, thus providing clinicians with more accurate prognostic information to inform treatment decisions (11–13). It has been increasingly reported that radiomics can be used for differentiating benign and malignant renal tumors, as well as discriminating high and low Fuhrman nuclear ccRCC (14, 15). However, to the best of our knowledge, no study has evaluated radiomics for its ability to predict the aggressive potential of ccRCC.
The purpose of this study was to investigate Radiomics Features (RFs) associated with progression-free survival (PFS) of RCC, assessing its incremental value over clinical factors, and to develop a visual nomogram in order to provide reference for individualized treatment and prognosis evaluation of RCC.
Materials and Methods
Patients
This retrospective study was approved by the Ethics Committee of Southern Medical University, and because of the retrospective nature of the analysis, the requirement of informed patient consent was waived. The data of patients with RCC who were treated in the Department of Urology, Nanfang Hospital, Southern Medical University, from March 2011 to March 2016 were retrospectively collected. Patients were randomly divided into a training set and validation set in a ratio of approximately 7:3. Clinical, pathological, and surgical data collected included age, sex, symptoms (low back pain, hematuria, emaciation, low-grade fever, cough, abdominal mass, and paraneoplastic syndrome), the interval from diagnosis to treatment, Karnofsky performance status (KPS) score, hemoglobin level, serum calcium level, neutrophil count, platelet count, maximum tumor size, stage, and pathological subtype, WHO/ISUP nuclear grade, overall tumor size, growth pattern, necrosis and calcification, treatment method, and adjuvant treatment. Inclusion criteria for the study were as follows: (1) Histological subtype ccRCC; (2) Clinical stage I–III (stage I–II, T1-2N0M0 and stage III, T1-2N1M0, T3N0-1M0); (3) Received radical nephrectomy or partial nephrectomy; and (4) Complete enhanced computed tomography (CT) imaging data, containing non-contrast phase (NCP), cortico-medullary phase (CMP), nephrographic phase (NP), and excretory phase (EP). Patients with a complete cystic renal tumor, positive tumor margins, and inadequate CT images were excluded.
CT Parameters
A Siemens 64-slice (Somatom Definition CT scanner; Siemens Medical Solutions Company, Malvern, PA, USA) CT scanners were used. A supine CT scan was performed from the diaphragmatic apex to the lower poles of both kidneys with breath holding using the following parameters: CT tube voltage = 120 kV, tube current = 150–320 mA, layer thickness = 5 mm, layer spacing = 5 mm, field of view = 360 mm, matrix = 512 × 512. After obtaining unenhanced images, Omnipaque (GE Healthcare) was injected into the anterior elbow vein with a high-pressure syringe at a dose of 2 ml/kg and an injection rate of 2.5 ml/s, with a maximum dose of 160 ml. Enhanced CT scanning was started at 25–30 s, 75–80 s, and 180–200 s.
Image Segmentation
Image segmentation was performed by 2 radiologists with 5 and 8 years, respectively, of abdominal imaging diagnosis experience using ITK-snap software. CT images were obtained from the PACS system, with a window width of 300–400 HU, a window level of 45–65 HU, and a slice thickness of 5 mm. The volume of the tumor was selected as the region of interest (ROI), and the edges were kept about 1 mm away from tumor edges to reduce interference from adjacent tissues (such as fat or normal renal tissue). Based on the threshold, areas with CT values less than −55 HU and greater than 350 HU pixels were filtered out. Intra-group and inter-group correlation coefficients (ICCs) were used to ensure stability and repeatability. First, 40 images were randomly segmented by the 2 physicians to assess reproducibility between groups. A week later, Doctor A repeated the same procedure to assess the reproducibility within the group. The results showed that ICC > 0.80 between groups and within groups; it means that the image segmentation was consistent, and the remaining image segmentation was performed by Dr. A.
Radiomics Feature Extraction and Selection
RFs were extracted, preprocessed, and filtered from segmented images using the PyRadiomics computing platform. First, the original CT images and 3D segmented images were imported into the platform for loading. Then, the image was preprocessed based on the Simple ITK software package embedded in the platform to ensure that the isotropic voxels of texture feature and shape feature are equidistant from adjacent bodies in all directions. Then, the preprocessed image was filtered based on the platform built-in filter, including PyWavelets and Simple ITK for wavelet filter and logarithmic filter, and NumPy for the remaining filters. Finally, four custom feature extraction methods based on the platform were used for radiomics feature extraction.
In the training set, feature selection was performed based on the stability LASSO algorithm. It performed through 100 times hierarchical 5-fold cross-validation. The LassoCV method used automatically found the optimal penalty coefficient α value through k-fold cross-verification. The criterion of α selection is to minimize generalization error. RFs were extracted in the training set with 5-fold cross-validations, and RFs with R2 > 0.8 were retained in the 5-fold cross-validation test set. The top 20% RFs of the 100 times 5-fold cross-validation were selected as the final features (Figure 1).
Development of Clinical Variable-Only, RFs-Only, and the Final Model
Multivariate Cox proportional hazard regression models were used to estimate the correlation coefficients of selected RFs to PFS, and to calculate the RFs-weighted sum as an independent variable in both training and validation sets. In the training set, univariate and multivariate Cox proportional hazard regression models were used to investigate factors associated with PFS. Independent variables with a value of p < 0.05 in univariate results were entered into a multivariate model, and variables that were significant in the multivariate model were considered factors associated with PFS, and estimated hazard ratios (HRs) were calculated. The final associated factors were used to build a multivariate Cox regression model. The optimal cutoff value of continuous variables was determined by the maximally selected rank statistics from the “maxstat” R package and used to predict PFS. Kaplan-Meier function was used to predict PFS for categorical variables.
Assessment of the Performance of Different Models
The probabilities were used as an independent continuous variable in receiver operating characteristic (ROC) curve analysis, and the C-index (area under curve the ROC curve [AUC]) and likelihood parameters were determined. A cutoff value was determined by maximizing the Youden index and used to predict PFS. ROC analysis was used in the validation set to assess the final model’s diagnostic effectiveness.
To compare the performance of a clinical variable-only model, RFs-only model, and the final model, ROC analysis and the DeLong test were used. The Hosmer–Lemeshow test was used to check the calibration. Decision curve analysis was used to observe the net benefits. Finally, index net reclassification improvement (NRI) and integrated discrimination improvement (IDI) were used to calculate the increment from the clinical variables-only model to the final model.
Statistical Analysis
Continuous data were presented as mean ± standard deviation, and categorical data were presented as number and percentage (%). For comparisons of means between groups, Student’s independent t-test or Mann–Whitney U test was used, depending on normality assumption. Categorical data were tested using chi-square test or Fisher’s exact text (if an expected value ≤ 5 was found). In all analysis, a 2-tailed value of p < 0.05 was considered to indicate statistical significance. Statistical analyses were performed using IBM SPSS version 25 software (IBM Corporation, Somers, New York). A nomogram was established using the associated factors in the training set with R statistical software (version 3.5.2) and the “rms” package. The surv_cutpoint function in the “surviminer” package finds the best cutoff value for a continuous variable. The “Predict ABEL” package is used to calculate the NRI and IDI. The decision curve analysis was also performed with R software and “rmda” package.
Results
Patient Clinical Characteristics and RFs of the Training and Validation Sets
A total of 175 patients were included in the analysis, with 125 in the training set and 50 in the validation set. The mean age of the training set was 52.31 ± 14.51 years, and that of the validation was 52.06 ± 13.19 years. The male:female ratio of the training set was 2.05:1, and that of the validation set was 1.63:1. The mean PFS of the training set was 55.83 ± 22.25 months, and that of the validation set was 61.06 ± 20.29 months. In the training set and validation set, 32.80% and 24.00% of patients experienced disease progression, respectively. Patient characteristics are summarized in Table 1, all clinical characteristics, follow-up results, and RFs of the training set and validation set were comparable (all, p > 0.05).
A total of 107 RFs were extracted from the 3D multiphase CT images of each phase of each patient. The RF features were categorized as follows: (1) first-order statistics, (2) shape-based features, and (3) texture features. Of the 107 RFs, there were18 first-order statistics features, 14 shape-based features, 24 gray-level cooccurrence matrix (GLCM) features, 16 gray-level size zone matrix (GLSZM) features, 16 gray-level run length matrix (GLRLM) features, 14 gray-level dependence matrix (GLDM) features, and 5 neighboring gray-tone difference matrix (NGTDM) features. A total of 428 (4×107) RFs were extracted from the 4 phase CT images. There were 6 RFs final selected, including X1, X2, X3, X4, X5, and X6. Table 2 shows RF’s designation, phase, abbreviation, classification, and description. A multivariate Cox regression model for PFS of the training set was established to integrate the RFs indices into a single index, the RFs-weighted sum, using the formula: RFs weighted sum = 5.5056 × 10-5 × X1 + 0.0013 × X2 + 0.0028 × X3 + 0.0053 × X4 − 10.4363 × X5 + 0.0048 × X6.
Predictive Model of PFS Using the Training Set
The univariate and multivariate Cox regression analyses results of the relations of independent variables to PFS in the training set are shown in Table 3. Variables significant in univariate results were entered into the multivariate model. Because of the high correlation between T stage and clinical stage (r = 0.63, p < 0.001), T stage was excluded in the final model. Thus, the final model for PFS was established using age, clinical stage, KPS score, and RFs-weighted sum (Table 4).
The results showed that patients with higher age, higher clinical stage, KPS score < 80, and larger RFs-weighted sum were more likely to have disease progression. The surv_cutpoint functions for age and RFs-weighted sum are shown in Figures 2A, B, the cutoff value of RFs-weighted sum and age was 48 and −0.73, and significant differences were found between age ≥ 48 and < 48 (p = 0.013) and between RFs-weighted sum ≥ −0.73 and < −0.73 (p < 0.001).The Kaplan–Meier survival functions for clinical stage and KPS score are shown in Figures 3A, B, and significant differences were found between clinical stage I–II and III (p < 0.001), and between KPS score ≥ 80 and < 80 (p = 0.006).
Figure 2 Surv_cutpoint function and survival analysis of PFS in the training set. (A) RFs-weighted sum. (B) Age.
Figure 3 Kaplan–Meier survival analysis of PFS in the training set. (A) Clinical stage. (B) KPS score.
ROC Analysis and Nomogram
Table 5 and Figure 4 show the results of ROC analysis of the final models of the training set and validation set. The C-index of training and validation models was 0.836 and 0.706, respectively. The Hosmer–Lemeshow test indicated that the final models of the training set (chi-square = 15.05, p = 0.058) and validation set (chi-square = 13.84, p = 0.086) were acceptable. A nomogram of the final model was established for clinical use (Figure 5), and included risk estimations of PFS, and 1-, 3-, and 5-year survival.
Figure 5 A nomogram for PFS was established that included age, clinical stage, KPS score, and RFs-weighted sum.
Comparisons of Clinical Variables-Only, RFs-Only, and Final Models
Three models were compared to investigate the importance of the RFs indices. The models included a clinical variables-only model (age, clinical stage, and KPS score), RFs-only model (RFs-weighted sum), and final model (age, clinical stage, KPS score, and RFs-weighted sum). The ROC analysis of the models is shown in Figure 6. It was found that the C-index of the final model was significantly larger than the clinical-only model (DeLong test, p = 0.008), but not significant compared to the RFs-only model (p = 0.402). However, no significance was found among the paired comparisons of models in the validation set (all p > 0.05).
Figure 6 ROC results of the clinical variables-only model, RFs only model, and final model of the training set (A) and validation set (B).
The decision curve analysis is shown in Figure 7. The net benefit was higher in the final model than in the clinical variables-only model all the time and also higher than that of the RFs-only model after a 0.4 risk-threshold. The NRI from the clinical variables-only model to the final model was 18.03% (Z = 1.80, p = 0.072), marginally significant; and the IDI was 19.08% (Z = 6.39, p < 0.001).
Figure 7 Decision curve analysis results of the clinical variables-only model, RFs only model, and final model.
Discussion
In this retrospective analysis, we developed a radiomics nomogram that incorporates three clinical factors and RFs-weighted sum for noninvasive, individualized prediction of PFS in patients with clinical stage I–III ccRCC, which can enable physicians to select reasonable treatment tactics and individualized monitoring protocols to improve clinical outcomes. To the best of our knowledge, this is the first prediction model developed to predict PFS of resectable ccRCC using CT-based radiomics. Through cross-validation and calibration, the RFs selection of this study ensures reliability and avoids over-fitting of the model (16). There were 41 progressive cases in the training set used for modeling, and the nomogram contains 4 factors. The variable selection is consistent with the 10–15 EPV (Event per Variable) criteria proposed by Peduzzi et al. (17), indicating that the model is reliable. The proposed radiomics nomogram demonstrated favorable discrimination in both the training set (C-index, 0.836) and validation set (C-index, 0.706), with high sensitivity, specificity, and accuracy.
Many factors, including clinical, anatomical, pathological, and molecular factors and treatment methods, are related to the prognosis of RCC. Early treatment, being asymptomatic, and a higher KPS score (>80) are very important for prolonging the survival of patients with RCC (18). Inflammatory cells (19–21), the TNM classification, and the histological factors included tumor nuclear grade, subtype, sarcomatoid features, microvascular invasion, tumor necrosis, and collection system invasion play a very important role in the tumor prognosis (22, 23). Many molecular factors such as CAIX, VEGF, HIF, Ki67, p53, p21, cell cycle PTEN, E-cadherin, CD44, and CXCR4, as well as other cell cycle and proliferation markers, may be associated with the prognosis of RCC (24–27). Partial nephrectomy is associated with improved survival of early RCC (28). Park et al. (29) reviewed preoperative laboratory data in 747 RCC patients and revealed that clinical information supporting aggressive ccRCC included an older age, larger size, lower hemoglobin, albumin, and calcium, as well as higher platelet and neutrophil. However, few radiologic parameters have been reported as prognostic factors of ccRCC in contrast to pathological markers. We enrolled these variables in this study; univariate Cox regression analysis showed that age, KPS score, neutrophils, tumor size, T stage, N stage, clinical stage, and calcification were associated with PFS, which was consistent with previous studies. However, affected by the RFs weighted sum, the clinical factors of the final model were only age, KPS score, and clinical stage. While RCC can be seen in all age groups, the median age of onset is 64 years, with a high incidence ranging from 50 to 70 years old. The mean age of patients in this study was 52 years, which is lower than that reported in the literature, and this may be related to the small sample size, patient race, and the relatively early disease stage.
Various scoring systems have been developed to predict the risk of postoperative recurrence in patients with RCC. The response evaluation criteria in solid tumors (RECIST) is the most commonly used prognostic evaluation method for tumors (30). However, it is impossible to predict the treatment effect before treatment. Other systems include the University of California, Los Angeles Integrated Staging System (UISS) (31), the Stage Size Grade and Necrosis (SSIGN) model (32), the Leibovich scoring system (33), the Kattan Nomogram (34), and the Karakiewic prognostic model (35). However, some of the parameters used in the model such as tumor necrosis and clinical presentation are subject to inter-observer variability. Different observation end points of different models result in different accuracy results in different studies, and ethnic differences and tumor diversity also limit the use of some systems; in various external validation samples, the results are not consistent (36, 37). Hence, further research and validation are needed.
RFs contain information about tumor heterogeneity and can reflect tumor phenotypes. Our study has filled a gap in the literature on PFS risk of stage I–III RCC in the setting of radiomics. In the recent literature, Radiomics nomogram has demonstrated excellent efficacy in differential diagnosis, nuclear grading, prognosis, and gene expression of RCC (38–43). Among the 6 RFs selected in this study, there were 3 features from the corticomedullary phase, suggesting that the corticomedullary phase may contain more abundant information to predict PFS. The results showed that RFs-weighted sum was an important factor that improved the diagnostic efficiency of the clinical variables-only model. The decision curve analysis revealed that using the radiomics nomogram to predict PFS in patients with stage I–III ccRCC presents more notable benefits than solely relying on clinical variables-only model.
There are several limitations to our study. First, owing to the limitation of the retrospective study and small number of cases, the follow-up time we used was at least 5 years. It would be more interesting to enroll patients without recurrence evidence for more than 5 or 10 years. Second, as a single-center study, the patient population was relatively homogeneous and small. During the 5-year recruiting period, there a large proportion of patients with stage I–II in this study (training set: 84%, verification set: 90%). A large-scale independent prospective multicenter study is needed to evaluate the generalizability of the results, and further work would focus on it.
In conclusion, this study presented a CT-based radiomics nomogram that showed satisfactory performance in predicting PFS in patients with stage I–III ccRCC, as a non-invasive and quantitative method that can be used as an efficient tool to complement individualized treatment.
Data Availability Statement
The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding author.
Ethics Statement
The studies involving human participants were reviewed and approved by the Ethics Committee of Southern 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
HZ, FY, and GW contributed to conception and design of the study. MC, LY, AQ, WC, and SY organized the database. HZ and FY performed the statistical analysis. HZ wrote the first draft of the manuscript. FY wrote sections of the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.
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.
References
1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. 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. doi: 10.3322/caac.21660
2. Lipworth L, Morgans AK, Edwards TL, Barocas DA, Chang SS, Herrell SD, et al. Renal Cell Cancer Histological Subtype Distribution Differs by Race and Sex. BJU Int (2016) 117(2):260–5. doi: 10.1111/bju.12950
3. Keegan KA, Schupp CW, Chamie K, Hellenthal NJ, Evans CP, Koppie TM. Histopathology of Surgically Treated Renal Cell Carcinoma: Survival Differences by Subtype and Stage. J Urol (2012) 188(2):391–7. doi: 10.1016/j.juro.2012.04.006
4. Capitanio U, Cloutier V, Zini L, Isbarn H, Jeldres C, Shariat SF, et al. A Critical Assessment of the Prognostic Value of Clear Cell, Papillary and Chromophobe Histological Subtypes in Renal Cell Carcinoma: A Population-Based Study. Bju Int (2009) 103(11):1496–500. doi: 10.1111/j.1464-410X.2008.08259.x
5. Paner GP, Stadler WM, Hansel DE, Montironi R, Lin DW, Amin MB. Updates in the Eighth Edition of the Tumor-Node-Metastasis Staging Classification for Urologic Cancers. Eur Urol (2018) 73(4):560–9. doi: 10.1016/j.eururo.2017.12.018
6. Tsui KH, Shvarts O, Smith RB, Figlin RA, deKernion JB, Belldegrun A. Prognostic Indicators for Renal Cell Carcinoma: A Multivariate Analysis of 643 Patients Using the Revised 1997 TNM Staging Criteria. J Urol (2000) 163(4):1090–5; quiz 1295. doi: 10.1016/s0022-5347(05)67699-9
7. Stewart SB, Thompson RH, Psutka SP, Cheville JC, Lohse CM, Boorjian SA, et al. Evaluation of the National Comprehensive Cancer Network and American Urological Association Renal Cell Carcinoma Surveillance Guidelines. J Clin Oncol (2014) 32(36):4059–65. doi: 10.1200/JCO.2014.56.5416
8. Ravaud A, Motzer RJ, Pandha HS, George DJ, Pantuck AJ, Patel A, et al. Adjuvant Sunitinib in High-Risk Renal-Cell Carcinoma After Nephrectomy. N Engl J Med (2016) 375(23):2246–54. doi: 10.1056/NEJMoa1611406
9. Amin MB, Greene FL, Edge SB, Compton CC, Gershenwald JE, Brookland RK, et al. The Eighth Edition AJCC Cancer Staging Manual: Continuing to Build a Bridge From a Population-Based to a More “Personalized” Approach to Cancer Staging. CA Cancer J Clin (2017) 67(2):93–9. doi: 10.3322/caac.21388
10. Martínez-Salamanca JI, Huang WC, Millán I, Bertini R, Bianco FJ, Carballido JA, et al. Prognostic Impact of the 2009 UICC/AJCC TNM Staging System for Renal Cell Carcinoma With Venous Extension. Eur Urol (2011) 57(1):120–7. doi: 10.1016/j.eururo.2010.10.001
11. van Griethuysen JJM, Fedorov A, Parmar C, Hosny A, Aucoin N, Narayan V, et al. Computational Radiomics System to Decode the Radiographic Phenotype. Cancer Res (2017) 77(21):e104–7. doi: 10.1158/0008-5472.CAN-17-0339
12. Ganeshan B, Nathan P, Juttla JK, Vinayan A, Miles KA. Assessment of Response to Tyrosine Kinase Inhibitors in Metastatic Renal Cell Cancer: CT Texture as a Predictive Biomarker. Radiology (2011) 261(1):165–71. doi: 10.1148/radiol.11110264
13. Lubner MG, Stabo N, Abel EJ, Del Rio AM, Pickhardt PJ. CT Textural Analysis of Large Primary Renal Cell Carcinomas: Pretreatment Tumor Heterogeneity Correlates With Histologic Findings and Clinical Outcomes. AJR Am J Roentgenol (2016) 207(1):96–105. doi: 10.2214/AJR.15.15451
14. Bektas CT, Kocak B, Yardimci AH, Turkcanoglu MH, Yucetas U, Koca SB, et al. Clear Cell Renal Cell Carcinoma: Machine Learning-Based Quantitative Computed Tomography Texture Analysis for Prediction of Fuhrman Nuclear Grade. Eur Radiol (2019) 29(3):1153–63. doi: 10.1007/s00330-018-5698-2
15. Feng Z, Rong P, Cao P, Zhou Q, Zhu W, Yan Z, et al. Machine Learning-Based Quantitative Texture Analysis of CT Images of Small Renal Masses: Differentiation of Angiomyolipoma Without Visible Fat From Renal Cell Carcinoma. Eur Radiol (2018) 28(4):1625–33. doi: 10.1007/s00330-017-5118-z
16. Babyak MA. What You See may Not be What You Get: A Brief, Nontechnical Introduction to Overfitting in Regression-Type Models. Psychosom Med (2004) 66(3):411–21. doi: 10.1097/01.psy.0000127692.23278.a9
17. Peduzzi P, Concato J, Kemper E, Holford TR, Feinstein AR. A Simulation Study of the Number of Events Per Variable in Logistic Regression Analysis. J Clin Epidemiol (1996) 49(12):1373–9. doi: 10.1016/s0895-4356(96)00236-3
18. Van Poppel H, Da Pozzo L, Albrecht W, Matveev V, Bono A, Borkowski A, et al. A Prospective, Randomised EORTC Intergroup Phase 3 Study Comparing the Oncologic Outcome of Elective Nephron-Sparing Surgery and Radical Nephrectomy for Low-Stage Renal Cell Carcinoma. Eur Urol (2011) 59(4):543–52. doi: 10.1016/j.eururo.2010.12.013
19. Grivennikov SI, Greten FR, Karin M. Immunity, Inflammation, and Cancer. Cell (2010) 140(6):883–99. doi: 10.1016/j.cell.2010.01.025
20. Szczerba BM, Castro-Giner F, Vetter M, Krol I, Gkountela S, Landin J, et al. Neutrophils Escort Circulating Tumour Cells to Enable Cell Cycle Progression. Nature (2019) 566(7745):553–7. doi: 10.1038/s41586-019-0915-y
21. Bambace NM, Holmes CE. The Platelet Contribution to Cancer Progression. J Thromb Haemost (2011) 9(2):237–49. doi: 10.1111/j.1538-7836.2010.04131.x
22. Zhang L, Zha Z, Qu W, Zhao H, Yuan J, Feng Y, et al. Tumor Necrosis as a Prognostic Variable for the Clinical Outcome in Patients With Renal Cell Carcinoma: A Systematic Review and Meta-Analysis. BMC Cancer (2018) 18(1):870. doi: 10.1186/s12885-018-4773-z
23. Fuhrman SA, Lasky LC, Limas C. Prognostic Significance of Morphologic Parameters in Renal Cell Carcinoma. Am J Surg Pathol (1982) 6(7):655–63. doi: 10.1097/00000478-198210000-00007
24. Sim SH, Messenger MP, Gregory WM, Wind TC, Vasudev NS, Cartledge J, et al. Prognostic Utility of Pre-Operative Circulating Osteopontin, Carbonic Anhydrase IX and CRP in Renal Cell Carcinoma. Br J Cancer (2012) 107(7):1131–7. doi: 10.1038/bjc.2012.360
25. Choueiri TK, Pal SK, McDermott DF, Morrissey S, Ferguson KC, Holland J, et al. A Phase I Study of Cabozantinib (XL184) in Patients With Renal Cell Cancer. Ann Oncol (2014) 25(8):1603–8. doi: 10.1093/annonc/mdu184
26. Choueiri TK, Escudier B, Powles T, Mainwaring PN, Rini BI, Donskov F, et al. Cabozantinib Versus Everolimus in Advanced Renal-Cell Carcinoma. N Engl J Med (2015) 373(19):1814–23. doi: 10.1056/NEJMoa1510016
27. Minardi D, Lucarini G, Milanese G, Di Primio R, Montironi R, Muzzonigro G. Loss of Nuclear BAP1 Protein Expression Is a Marker of Poor Prognosis in Patients With Clear Cell Renal Cell Carcinoma. Urol Oncol (2016) 34(8):338. e11–8. doi: 10.1016/j.urolonc.2016.03.006
28. Boorjian S. Commentary on “Long-Term Survival Following Partial vs. Radical Nephrectomy Among Older Patients With Early-Stage Kidney Cancer”. Tan HJ, Norton EC, Ye Z, Hafez KS, Gore JL, Miller DC. Urol Oncol (2012) 31 (1):126–27. doi: 10.1016/j.urolonc.2012.11.004
29. Park YH, Baik KD, Lee YJ, Ku JH, Kim HH, Kwak C. Late Recurrence of Renal Cell Carcinoma >5 Years After Surgery: Clinicopathological Characteristics and Prognosis. BJU Int (2012) 110(11 Pt B):E553–8. doi: 10.1111/j.1464-410X.2012.11246.x
30. Hillman SL, An MW, O’Connell MJ, Goldberg RM, Schaefer P, Buckner JC, et al. Evaluation of the Optimal Number of Lesions Needed for Tumor Evaluation Using the Response Evaluation Criteria in Solid Tumors: A North Central Cancer Treatment Group Investigation. J Clin Oncol (2009) 27(19):3205–10. doi: 10.1200/JCO.2008.18.3269
31. Patard JJ, Kim HL, Lam JS, Dorey FJ, Pantuck AJ, Zisman A, et al. Use of the University of California Los Angeles Integrated Staging System to Predict Survival in Renal Cell Carcinoma: An International Multicenter Study. J Clin Oncol (2004) 22(16):3316–22. doi: 10.1200/JCO.2004.09.104
32. Sorbellini M, Kattan MW, Snyder ME, Reuter V, Motzer R, Goetzl M, et al. A Postoperative Prognostic Nomogram Predicting Recurrence for Patients With Conventional Clear Cell Renal Cell Carcinoma. J Urol (2005) 173(1):48–51. doi: 10.1097/01.ju.0000148261.19532.2c
33. Leibovich BC, Blute ML, Cheville JC, Lohse CM, Frank I, Kwon ED, et al. Prediction of Progression After Radical Nephrectomy for Patients With Clear Cell Renal Cell Carcinoma: A Stratification Tool for Prospective Clinical Trials. Cancer (2003) 97(7):1663–71. doi: 10.1002/cncr.11234
34. Kattan MW, Reuter V, Motzer RJ, Katz J, Russo P. A Postoperative Prognostic Nomogram for Renal Cell Carcinoma. J Urol (2001) 166(1):63–7. doi: 10.1097/00005392-200107000-00015
35. Karakiewicz PI, Briganti A, Chun FK, Trinh QD, Perrotte P, Ficarra V, et al. Multi-Institutional Validation of a New Renal Cancer-Specific Survival Nomogram. J Clin Oncol (2007) 25(11):1316–22. doi: 10.1200/JCO.2006.06.1218
36. Beisland C, Gudbrandsdottir G, Reisæter LA, Bostad L, Wentzel-Larsen T, Hjelle KM. Contemporary External Validation of the Leibovich Model for Prediction of Progression After Radical Surgery for Clear Cell Renal Cell Carcinoma. Scand J Urol (2015) 49(3):205–10. doi: 10.3109/21681805.2014.980844
37. Pichler M, Hutterer GC, Chromecki TF, Jesche J, Kampel-Kettner K, Rehak P, et al. External Validation of the Leibovich Prognosis Score for Nonmetastatic Clear Cell Renal Cell Carcinoma at a Single European Center Applying Routine Pathology. J Urol (2011) 186(5):1773–7. doi: 10.1016/j.juro.2011.07.034
38. Nie P, Yang G, Wang Z, Yan L, Miao W, Hao D, et al. A CT-Based Radiomics Nomogram for Differentiation of Renal Angiomyolipoma Without Visible Fat From Homogeneous Clear Cell Renal Cell Carcinoma. Eur Radiol (2020) 30(2):1274–84. doi: 10.1007/s00330-019-06427-x
39. Kang B, Sun C, Gu H, Yang S, Yuan X, Ji C, et al. T1 Stage Clear Cell Renal Cell Carcinoma: A CT-Based Radiomics Nomogram to Estimate the Risk of Recurrence and Metastasis. Front Oncol (2020) 10:579619. doi: 10.3389/fonc.2020.579619
40. Zheng Z, Chen Z, Xie Y, Zhong Q, Xie W. Development and Validation of a CT-Based Nomogram for Preoperative Prediction of Clear Cell Renal Cell Carcinoma Grades. Eur Radiol (2021) 31(8):6078–86. doi: 10.1007/s00330-020-07667-y
41. Huang Y, Zeng H, Chen L, Luo Y, Ma X, Zhao Y. Exploration of an Integrative Prognostic Model of Radiogenomics Features With Underlying Gene Expression Patterns in Clear Cell Renal Cell Carcinoma. Front Oncol (2021) 11:640881. doi: 10.3389/fonc.2021.640881
42. Zhang H, Yin F, Chen M, Qi A, Lai Z, Yang L, et al. A Reliable Prediction Model for Renal Cell Carcinoma Subtype Based on Radiomic Features From 3d Multiphase Enhanced CT Images. J Oncol (2021) 2021:6595212. doi: 10.1155/2021/6595212
Keywords: renal cell carcinoma (RCC), Radiomics, CT, progression-free survival (PFS), predict model, artificial intelligence
Citation: Zhang H, Yin F, Chen M, Yang L, Qi A, Cui W, Yang S and Wen G (2022) Development and Validation of a CT-Based Radiomics Nomogram for Predicting Postoperative Progression-Free Survival in Stage I–III Renal Cell Carcinoma. Front. Oncol. 11:742547. doi: 10.3389/fonc.2021.742547
Received: 16 July 2021; Accepted: 29 November 2021;
Published: 27 January 2022.
Edited by:
Giorgio Russo, Institute of Bioimaging and Molecular Physiology (CNR), ItalyReviewed by:
Yuanshen Zhao, Chinese Academy of Sciences (CAS), ChinaArie Carneiro, Albert Einstein Israelite Hospital, Brazil
Dapeng Hao, The Affiliated Hospital of Qingdao University, China
Copyright © 2022 Zhang, Yin, Chen, Yang, Qi, Cui, Yang and Wen. 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: Ge Wen, bTEzMzYwMDIyMTY2QDE2My5jb20=
†These authors have contributed equally to this work