- Research article
- Open Access
CT texture analysis: a potential tool for predicting the Fuhrman grade of clear-cell renal carcinoma
Cancer Imagingvolume 19, Article number: 6 (2019)
The purpose of this study was to analyze the image heterogeneity of clear-cell renal-cell carcinoma (ccRCC) by computer tomography texture analysis and to provide new objective quantitative imaging parameters for the pre-operative prediction of Fuhrman-grade ccRCC.
A retrospective analysis of 131 cases of ccRCCs was performed by manually depicting tumor areas. Then, histogram-based texture parameters were calculated. The texture-feature values between Fuhrman low- (Grade I-II) and high-grade (Grade III-IV) ccRCCs were compared by two independent sample t-tests (False Discovery Rate correction), and receiver operating characteristic curve (ROC) was used to evaluate the efficacy of using texture features to predict Fuhrman high- and low-grade ccRCCs.
There were no statistical differences for any texture parameters without filtering (p > 0.05). There was a statistically significant difference between the entropy (fine) of the corticomedullary phase and the entropy (fine and coarse) of the nephrographic phase after Laplace of Gaussian filtering. The area under the ROC of the entropy was between 0.74 and 0.83.
Computer tomography texture features can predict the Fuhrman grading of ccRCC pre-operatively, with entropy being the most important imaging marker for clinical application.
Renal carcinoma is the most common adult malignant epithelial tumor of the kidney . The most common pathological type is clear-cell renal cell carcinoma (ccRCC), which accounts for 60–85% of renal carcinoma . Fuhrman grade is currently the most common classification method for ccRCC. It classifies renal carcinoma into four different pathological grades according to the size, shape, staining, and presence or absence of nucleoli of the nuclei of cancer cells. Among them, Grades I–II are low and Grades III–IV are high . Current studies have shown that the grade of Fuhrman classification is closely related to the growth rate of tumors and the prognosis of patients [4,5,6]. High-grade tumors have higher invasive capacities, higher possibility of metastasis, and poor prognosis .
Fuhrman grade is an independent prognostic factor . Currently Nephron-sparing surgery has been used more and more in the treatment of T1-stage ccRCC. But T1 tumors with a high Fuhrman grade had a higher malignant potential for subsequent tumor recurrence . Therefore, pre-operative confirmation of the Fuhrman grading of ccRCC is critical for the selection of surgical options and prognosis.
Computer tomography (CT) texture analysis (TA) is a method used to quantitatively analyze spatial heterogeneity of lesions on CT images . It is filtered by a Laplace of Gaussian (LoG) operation to reduce photon noise while highlighting the changes in the characteristics of the grayscale images [11, 12]. It has been applied to the identification of renal tumors , intraductal papillary mucinous neoplasm classification of the pancreas , and benign and malignant prediction of mediastinal lymph nodes .
This study aimed to explore the correlation between the texture features of ccRCC and Fuhrman grading by CTTA and to set up objective and effective quantitative parameters of imaging to evaluate the heterogeneity of ccRCC, and further help to predict the pre-operative Fuhrman grading and improve the prognosis of the patients.
This retrospective study was approved by the Institutional Review Board of First Affiliated Hospital of College of Medical Science, Zhejiang University, and the requirement for informed consent was waived.
The study retrospectively analyzed 131 patients with ccRCC from June 2013 to December 2017. The inclusion criteria were the following: (1) patients underwent pre-operative renal enhanced CT scans, including noncontrast phase, corticomedullary phase, and nephrographic phase; (2) the patients were scanned using the same scanning protocol on the same CT scaner; (3) the lesions showed at least seven slices on the CT axial images; (4) tumors were surgically removed or confirmed as ccRCC via percutaneous biopsy; (5) cases were not treated before CT scan. Exclusion criteria were as follows: (1) cases showing Fuhrman-grade ambiguity, such as undetermined between Grade II and III, and (2) CT images having obvious artifacts.
According to the Fuhrman classification, we divided the cases into a low-grade group (Grade I–II) and a high-grade group (Grade III–IV). There were 77 patients in the low-grade group, 59 males, and 16 females, aged 25–81 (53.82 ± 15.88) years; there were 54 high-grade patients, 41 males and 13 females, aged 48–77 (61.59 ± 11.52) years. All cases in this study were unilateral lesions.
All cases were scanned using a GE Lightspeed 64-row CT machine with a tube voltage of 120 kV and a tube current of 150–350 mA, and the non-ionic contrast agent (iodine content 300 mg/mL) was infused into the peripheral vein at an infusion rate of 3 mL/s and an infusion dose of 80–100 mL. The scan ranged from the adrenal gland to the inferior pole of the kidney with a thickness of 3 mm per layer. If the tumor had a large diameter, the scan was expanded to include the entire tumor. The corticomedullary phase was delayed by 30 s and the nephrographic phase was delayed by 90 s. We manually traced the region of interest (ROI) based on the boundary of the tumor on the CT images of the corticomedullary and nephrographic phases. The boundary of the lesion could not be accurately identified on the noncontrast phase images, which was not used in this study. We used two radiologists with 10 years of work experience to do ROI delineation, and two radiologists segmentation results were superimposed to generate the final ROI. Each ROI was validated by a senior radiologist with 20 years of experience. ROI sketching are completed in the self-compiled MATLAB (version R2013b, MathWorks, USA) program.
The data processing of texture analysis was also based on MATLAB, and the algorithm was basically consistent with the current published partial CT texture-analysis research [16,17,18,19]. The texture analysis included the following two processes: (1) filtering of the original images using LoG filtering. In this study, we used sigma (s) values of 1.0 (fine-scale filtration) and 2.5 (coarse-scale filtration), and the original unfiltered CT images were preserved (Fig.1). (2) Quantification of image texture. Currently, there are several texture parameters. In this study, only the most representative first-order histogram textures were calculated, including mean gray intensity, entropy, standard deviation, skewness, and kurtosis. The mean gray intensity is the average of the pixel intensities, and the entropy is the amount of information contained in the aggregated features of the grayscale distribution in the images, reflecting the complexity of the images. The standard deviation is the degree of dispersion of individual pixel values in the grayscale distribution of the images. Skewness is a quantitative indicator of histogram asymmetry, and kurtosis is a quantitative indicator of histogram spikes.
This study used version 3.3.2 of R software (http://www.R-project.org) for statistical analysis of the data. The texture parameters of the Fuhrman low−/high-grade groups were compared using a two-sample t-test. In order to avoid false positives after multiple comparisons, we also performed multiple comparison corrections based on False Discovery Rate (FDR), with a statistical difference of p < 0.05. Subsequently, the receiver operating characteristic (ROC) curve analysis was performed on the statistically significant texture parameters, and the area under the ROC curve (AUC) and sensitivity and specificity values were calculated to evaluate the effectiveness of each parameter in distinguishing the Fuhrman grading.
Table 1 shows that no texture parameters of the two groups were statistically different (p > 0.05) without the addition of filtering. The mean gray intensity and skewness showed no significant differences after LoG filtering (p > 0.05). There was a statistically significant difference in standard deviation of the corticomedullary phase (coarse filtration) (p < 0.05), and there were significant differences in kurtosis (coarse filtration) of the corticomedullary and nephrographic phases (p < 0.05). There were statistically significant differences in all entropy after LoG filtering (p < 0.01). However, entropy (fine filtration) of the corticomedullary phase and entropy (fine and coarse filtration) of the nephrographic phase showed significant statistical differences after FDR corrections.
Table 2 lists the AUC, the optimal threshold, and the corresponding accuracy after FDR correction of the entropy-predicted Fuhrman low/high grades. The corticomedullary phase entropy (fine filtration) had an AUC of 0.74, sensitivity of 0.76, specificity of 0.65 and cutoff value of 4.66. The entropy (fine filtration) of the nephrographic phase had an AUC of 0.80, sensitivity of 0.95, specificity of 0.54 and cutoff value of 4.27. The nephrographic phase entropy (coarse filtration) had an AUC of 0.83, sensitivity of 0.82, specificity of 0.77 and cutoff value of 2.55(Fig.2).
There is an important relationship between the Fuhrman classification of ccRCC and prognosis, so several non-invasive methods have been used to predict the Fuhrman grading of ccRCC. MR(Magnetic Resonance)has functional imaging studies based on diffusion and perfusion [20, 21], but MR is costly, which leads to lower popularity. Wang et al. used the RENAL nephrometry score based on CT images for the prediction . There have also been a large number of CT-based semi-quantitative and quantitative studies [23, 24], which indicates that CT is a convenient and effective method for predicting the Fuhrman classification of ccRCC. In this study, we further explored the CT data. After applying the LoG filter to pre-process the images, entropy was found to be a very important factor in predicting the texture parameters of Fuhrman grading.
The results suggest that there was no statistically significant difference in any texture parameters of the unfiltered images. Previously, Huhdanpaa et al. also performed Fuhrman-grade prediction based on the histogram parameters of ordinary CT images. The results show that only the inter-quartile range of the nephrographic phase has significant statistical significance. There was no statistical difference in mean, standard deviation, skewness, and kurtosis , which was consistent with our unfiltered results. Further, we carried out LoG filtering and found that some texture parameters showed significant differences after LoG filtering. LoG filtering is an advanced image-filtering method that combines Laplacian filtering and Gaussian filtering. Laplacian filtering can highlight the grayscale mutation region in the image and enhance the grayscale contrast. Gaussian operators can suppress the noise brought about by the Laplacian operator [11, 26]. The low filter value corresponds to the fine texture features and the high filter value corresponds to the coarse texture features. Meghan and others also found that only LoG filtered texture parameters were significantly correlated with cirrhosis grade , and our results once again confirm that LoG filtering can improve the ability to detect disease heterogeneity.
Heterogeneity is an important feature of malignant tumors and is closely related to the adverse biological processes of tumors. CTTA is a technique for effectively assessing tumor heterogeneity [12, 28]. Zhu et al. retrospectively evaluated 255 cases of ccRCC and found that low enhancement of medullary tumors was an independent factor in predicting high-grade tumors, but their experiments required higher ROI and extensive experience of the physician in measurements to avoid areas of obvious necrosis, large blood vessels, and calcification. In addition, they only selected one slice of the tumor imaging, but not the entire tumor . The method of Zhu et al. relied too much on personal experience, and the consistency between the measures was difficult to guarantee. Pichler et al. believed necrosis to be an independent prognostic indicator of ccRCC, so that avoiding necrotic regions almost completely ignored the heterogeneity of ccRCC, which was an important feature of tumors . Hebert et al. performed Fuhrman-grade analysis based on single slices and the entire lesion, and they found that the pathological grade was not associated with the enhancement parameters of a single slice, while the enhancement parameters of the volume measurement were related to the grade . Choosing the suitable slice is important when only analyzing one slice of image. However, the result of a single slice cannot fully reflect the heterogeneity of the entire tumor. Therefore, this study used volumetric measurements, which were relatively cumbersome, but the assessment of heterogeneity was very important.
Our study aimed to reveal a group of common first-order parameters to reflect tumor heterogeneity, including standard deviation, mean gray intensity, skewness, kurtosis, and entropy. It is worth mentioning that we have introduced entropy, which is an abstract quantitative indicator of texture irregularity and chaos, reflecting the degree of disorder of the images . All of these parameters can reflect the overall texture features of the tumor to some extent, and are used to objectively and quantitatively analyze the heterogeneity of ccRCC. The final result verifies our assumption that entropy is indeed an independent and excellent texture parameter. However, counterintuitively, the average entropy of low-grade tumors is higher than that of high-grade tumors; that is, the higher heterogeneity was found in low-grade tumors. However, this result is consistent with the results of Heidi et al. and Kousei et al [24, 31]. We speculate that it may be due to the fact that the micro-vessels of the lower-grade tumors are more abundant and the degree of intra-tumor intensification is higher. In addition, the concept of entropy is relatively abstract. The larger entropy value indicates a more random density distribution in the lesions, and the smaller entropy value indicates a relatively uniform density distribution. High-grade tumors, due to their relatively large liquefaction necrosis range, will result in reduced entropy.
In this study, there was a statistically significant difference in the standard deviation (coarse filtration) of the corticomedullary phase and the kurtosis (coarse filtration) of the corticomedullary and nephrographic phases before FDR correction, but they were not corrected by multiple comparisons. Chalkidou pointed out that there were serious false positives in the current CCTA studies . Many studies have not performed multiple comparison corrections, resulting in poor reproducibility, and, accordingly, the results cannot be reproduced. In order to reduce the typeIerror, the study performed a strict FDR correction on the results. Surprisingly, entropy is statistically significant even after FDR correction, and in subsequent ROC curve analysis the entropy also showed a good degree of discrimination in tumor grades. The entropy calculation based on ordinary CT images is simple and convenient, and can objectively quantify the heterogeneity of tumors, and has the prospect of clinical transformation.
The limitations of this study were as follows. (1) This study was a single-center retrospective study with no external data validation. (2) The sample size was still relatively small. In fact, the number of low-grade Fuhrman tumors is very large, while the number of high-grade tumors is relatively small, and the number of cases between groups is easily imbalanced. (3) Since the tumor boundary is manually drawn, the interference of the volume effect cannot be completely avoided. (4) The WHO/International Society of Urologic Pathology grading system will replace Fuhrman grading . In the future, we plan to collect cases based on new grading.
This results of this study show that LoG filtering can reduce image noise while highlighting the degree of unevenness inside a tumor, which is a very effective image pre-processing method. The entropy value can objectively quantify tumor heterogeneity, and the calculation is simple and convenient. It can be directly applied to conventional CT to assist clinicians in predicting Fuhrman grade.
The area under the curve
Clear-cell renal-cell carcinoma
False Discovery Rate
Laplace of Gaussian
Receiver operating characteristic
Region of interest
Capitanio U, Montorsi F. Renal cancer. Lancet. 2015;387:894.
Patard JJ, Leray E, Riouxleclercq N, Cindolo L, Ficarra V, Zisman A, et al. Prognostic value of histologic subtypes in renal cell carcinoma: a multicenter experience. J Clin Oncol. 2005;175:481–2.
Fuhrman SA, Lasky LC, Limas C. Prognostic significance of morphologic parameters in renal cell carcinoma. Int J Clin Pract. 2004;58:333–6.
Minardi D, Lucarini G, Mazzucchelli R, Milanese G, Natali D, Galosi AB, et al. Prognostic role of Fuhrman grade and vascular endothelial growth factor in pT1a clear cell carcinoma in partial nephrectomy specimens. J Urol. 2005;174:1208–12.
Li XS, Yao L, Gong K, Yu W, He Q, Zhou LQ, et al. Growth pattern of renal cell carcinoma (RCC) in patients with delayed surgical intervention. J Cancer Res Clin Oncol. 2012;138:269–74.
Sahni VA, Silverman SG. Imaging management of incidentally detected small renal masses. Semin Interv Radiol. 2014;31:009–19.
Donat SM, Diaz M, Bishoff JT, Coleman JA, Dahm P, Derweesh IH, et al. Follow-up for clinically localized renal neoplasms: AUA guideline. J Urol. 2013;190:407–16.
Ljungberg B, Cowan NC, Hanbury DC, Hora M, Kuczyk MA, Merseburger AS, et al. EAU guidelines on renal cell carcinoma: the 2010 update. Eur Urol. 2010;58:398–406.
Suzuki K, Mizuno R, Mikami S, Tanaka N, Kanao K, Kikuchi E, et al. Prognostic significance of high nuclear grade in patients with pathologic T1a renal cell carcinoma. Jpn J Clin Oncol. 2012;42:831–5.
Lubner MG, Smith AD, Sandrasegaran K, Sahani DV, Pickhardt PJ. CT texture analysis: definitions, applications, biologic correlates, and challenges. Radiographics. 2017;37:1483–503.
Miles KA, Ganeshan B, Hayball MP. CT texture analysis using the filtration-histogram method: what do the measurements mean? Cancer imaging the official publication of the international Cancer imaging. Society. 2013;13:400.
Ganeshan B, Miles KA. Quantifying tumour heterogeneity with CT. Cancer imaging the official publication of the international cancer imaging. Society. 2013;13:140–9.
Raman SP, Chen Y, Schroeder JL, Huang P, Fishman EK. CT texture analysis of renal masses: pilot study using random forest classification for prediction of pathology. Acad Radiol. 2014;21:1587–96.
Hanania AN, Bantis LE, Feng Z, Wang H, Tamm EP, Katz MH, et al. Quantitative imaging to evaluate malignant potential of IPMNs. Oncotarget. 2016;7(52):85776–84.
Andersen MB, Harders SW, Ganeshan B, Thygesen J, Torp Madsen HH. Rasmussen F. CT texture analysis can help differentiate between malignant and benign lymph nodes in the mediastinum in patients suspected for lung cancer. Acta Radiol. 2016;3:669–76.
Rao SX, Lambregts DM, Schnerr RS, Beckers RC, Maas M, Albarello F, et al. CT texture analysis in colorectal liver metastases: a better way than size and volume measurements to assess response to chemotherapy? United European Gastroenterol. 2016;4:257–263.
Ba-Ssalamah A, Muin D, Schernthaner R, Kulinna-Cosentini C, Bastati N, Stift J, et al. Texture-based classification of different gastric tumors at contrast-enhanced CT. Eur J Radiol. 2013;82:E537–E43.
Ganeshan B, Goh V, Mandeville HC, Ng QS, Hoskin PJ, Miles KA. Non-small cell lung cancer: histopathologic correlates for texture parameters at CT. Radiology. 2013;266:326–36.
Rao SX, Lambregts DM, Schnerr RS, Van OW, van Nijnatten TJ, Martens MH, et al. Whole-liver CT texture analysis in colorectal cancer: does the presence of liver metastases affect the texture of the remaining liver? United European Gastroenterol J. 2014;2:530–8.
Rosenkrantz AB, Niver BE, Fitzgerald EF, Babb JS, Chandarana H, Melamed J. Utility of the apparent diffusion coefficient for distinguishing clear cell renal cell carcinoma of low and high nuclear grade. AJR Am J Roentgenol. 2010;185:2099.
Lanzman RS, Robson PM, Sun MR, Patel AD, Mentore K, Wagner AA, et al. Arterial spin-labeling MR imaging of renal masses: correlation with histopathologic findings. Radiology. 2012;265:799.
Wang HK, Zhu Y, Yao XD, Zhang SL, Dai B, Zhang HL, et al. External validation of a nomogram using RENAL nephrometry score to predict high grade renal cell carcinoma. J Urol. 2012;187:1555–60.
Coy H, Young JR, Douek ML, Pantuck A, Brown MS, Sayre J, et al. Association of qualitative and quantitative imaging features on multiphasic multidetector CT with tumor grade in clear cell renalcell carcinoma. Abdom Radiol. 2019;44:180.
Kousei I, Leite LV, Pakalniskis MG, Lee DK, Holanda DG, Kuehn DM. Tumor grade of clear cell renal cell carcinoma assessed by contrast-enhanced computed tomography. Springerplus. 2014;3:1–7.
Huhdanpaa H, Hwang D, Cen S, Quinn B, Nayyar M, Zhang X, et al. CT prediction of the Fuhrman grade of clear cell renal cell carcinoma (RCC): towards the development of computer-assisted diagnostic method. Abdom Imaging. 2015;40:3168–74.
Shen Q, Shan Y, Hu Z, Chen W, Yang B, Han J, et al. Quantitative parameters of CT texture analysis as potential markersfor early prediction of spontaneous intracranial hemorrhage enlargement. Eur Radiol. 2018;28:4389–96.
Lubner MG, Malecki K, Kloke J, Ganeshan B, Pickhardt PJ. Texture analysis of the liver at MDCT for assessing hepatic fibrosis. Abdom Radiol. 2017;42:1–10.
Davnall F, Yip CSP, Ljungqvist G, Selmi M, Ng F, Sanghera B, et al. Assessment of tumor heterogeneity: an emerging imaging tool for clinical practice? Insights Into Imaging. 2012;3:573–89.
Zhu YH, Wang X, Zhang J, Chen YH, Kong W, Huang YR. Low enhancement on multiphase contrast-enhanced CT images: an independent predictor of the presence of high tumor grade of clear cell renal cell carcinoma. AJR Am J Roentgenol. 2014;203:295–300.
Pichler M, Hutterer GC, Chromecki TF, Jesche J, Kampelkettner K, Rehak P, et al. Histologic tumor necrosis is an independent prognostic Indicator for clear cell and papillary renal cell carcinoma. Am J Clin Pathol. 2012;137:283.
Vargas HA, Delaney HG, Delappe EM, Wang Y, Zheng J, Moskowitz CS, et al. Multiphasic contrast-enhanced MRI: single-slice versus volumetric quantification of tumor enhancement for the assessment of renal clear-cell carcinoma fuhrman grade. J Magn Reson Imaging Jmri. 2013;37:1160–7.
Chalkidou A, O'Doherty MJ, Marsden PK. False discovery rates in PET and CT studies with texture features: a systematic review. PLoS One. 2015;10:e0124165.
Srigley JR, Delahunt B, Eble JN, Egevad L, Epstein JI, Grignon D, et al. The International Society of Urological Pathology (ISUP) Vancouver classification of renal neoplasia. Am J Surg Pathol. 2013;37:1469–89.
This study has received funding by the Department of Health of Zhejiang Province, China (No. 2017KY051).
This study has received funding by the Department of Health of Zhejiang Province, China (No. 2018KY582).
This study has received funding by the Department of Health of Zhejiang Province, China (No. 2019KY551).
Availability of data and materials
Ethics approval and consent to participate
This retrospective study was approved by the Institutional Review Board of First Affiliated Hospital of College of Medical Science, Zhejiang University, and the requirement for informed consent was waived.
Consent for publication
The authors declare that they have no competing interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.