Impact of hepatocellular carcinoma heterogeneity on computed tomography as a prognostic indicator
We assessed the relationship between the heterogeneity of HCC on preoperative non-contrast-enhanced CT and patient prognosis. The heterogeneity of CT images from 122 patients was assessed and texture feature parameters such as mean, standard deviation (SD), entropy, mean of the positive pixels (MPP), skewness, and kurtosis were obtained using filtration. The relationship between CT texture features and 5-year overall survival (OS) or disease-free survival (DFS) was assessed. Multivariate regression analysis was performed to evaluate the independence of texture feature from clinical or pathological parameters. The Kaplan-Meier curves for OS or DFS was significantly different between patient groups dichotomized by cut-off values for all CT texture parameters with filtration at at least one filter level. Multivariate regression analysis showed the independence of most CT texture parameters on clinical and pathological parameters for OS with filtration at at least one filter level and without filtration except kurtosis. SD, entropy, and MPP with coarse filter, and skewness without filtration showed a significant correlation for DFS. CT texture features of non-contrast-enhanced CT images showed a relationship with HCC prognosis. Multivariate regression analysis showed the possibility of CT texture feature increase the prognostic prediction of HCC by clinical and pathological information.
The incidence of hepatocellular carcinoma (HCC) has increased due to the spread of hepatitis, and HCC is now the second leading cause of cancer-related death worldwide1. Because HCC has a high potential for vascular invasion, metastasis, and recurrence, its prognosis is poor. Among several therapies used to treat HCC, hepatic resection is considered the mainstay of curative therapy for localized HCC2,3. However, due to its high rate of recurrence after resection, preoperative prognostic prediction of HCC is important for the appropriate management of patients.
Heterogeneity is widely recognized as a feature of malignancy associated with adverse tumor biology, and it is assessed on medical images using a texture-analysis technique4. On computed tomography (CT), photon noise, caused by fluctuations in the number of photons in incident X-rays, obscures biologic heterogeneity4. CT texture analysis using filters to reduce photon noise is used for a variety of purposes, including assessment of malignant tumor prognosis5.
CT is widely used for the detection and assessment of HCC and is an indispensable preoperative examination. If information related to HCC prognosis is available from non-contrast-enhanced CT images, preoperative CT can provide additional meaningful data. The objective of this study was to assess the relationship between HCC heterogeneity and prognosis after the resection of HCC using non-contrast-enhanced CT texture analysis with a filter technique.
This study was approved by the Research Ethics Committee of our institution as a retrospective data analysis for medical imaging-based diagnoses, and the requirement for informed consent was waived by the Committee (26-73-1126). Additionally, all experiments and methods in this study were performed in accordance with the Declaration of Helsinki.
A total of 169 consecutive patients who underwent initial hepatectomy for HCC between January 2004 and September 2009 were included in this study. This study was approved by the Research Ethics Committee of the institution. A total of 16 patients were excluded due to the unavailability of CT images on the picture archiving system. We excluded a further 28 patients who underwent transcatheter arterial embolization (TACE) prior to CT, because the high density of deposited iodized oil may affect the texture features. A total of three patients with lesions too small (<1 cm) to evaluate on CT were additionally excluded. Finally, the study group consisted of 122 patients: 92 males (mean age, 63.5 years; range, 25–82 years) and 30 females (mean age, 69.2 years; range, 53–81 years). A total of 37 patients underwent TACE between CT and surgery (Table 1). Clinical patient’s information such as Child-Pugh score, serum levels of alpha-fetoprotein (AFP), tumor differentiation, pathological stage (pStage), or presence/absence of venous invasion was examined using database of surgical records.
CT studies were performed using a 4-, 16-, or 64-detector row CT scanner (Light Speed QX/I, LightSpeed Ultra, LightSpeed VCT, GE Medical Systems, Milwaukee, WI, USA) or a 16-detector row CT scanner (Aquilion 16, Toshiba Medical Systems Corp., Tochigi, Japan). The same clinical protocol was used: 120 kV, 180–280 mAs depending on body habitus, a matrix of 512, and a field of view of 350–400 mm. Non-contrast CT was performed using a 5-mm contiguous axial section to encompass the whole liver. Non-contrast CT was mostly performed as part of the contrast-enhanced study.
The location of the HCC was defined on unenhanced CT according to the surgical records, and the slice with the largest lesion was selected by one of the research conductors (a radiologist with 22 years of experience with abdominal CT). Images were loaded onto a workstation for further texture analysis. Texture analysis was performed by a single observer (a radiologist with 12 years of experience with abdominal CT) who was blinded to the clinical outcome, using TexRAD (TexRAD Ltd., www.texrad.com part of Feedback Plc., Cambridge, UK), a proprietary research software algorithm developed to visualize and quantify the texture properties of tissues from medical imaging scans. The region of interest (ROI) was initially delineated around the tumor periphery section by an observer and refined by excluding areas of air using a thresholding procedure that removed any pixels with attenuation values below −50 Hounsfield unit from the analysis. The assessment of texture feature comprised an initial filtration step in which a Laplacian of Gaussian spatial band-pass filter was used to selectively extract features of different sizes and intensity variations, at three different frequency scales: fine (features approximately 2 mm in width), medium (features approximately 4 mm in width) and coarse (features approximately 6 mm in width) (Fig. 1). The heterogeneity within the ROI was assessed with and without filtration, and texture parameters were calculated as follows: the mean grey-level intensity (mean), variation/dispersion from the mean grey-level intensity (standard deviation, SD), the irregularity or complexity of the grey signal (entropy), the average intensity of the positive grey-level signal pixel values within the ROI (mean of the positive pixels, MPP), asymmetry of the distribution (skewness), and pointiness or peakedness of the distribution (kurtosis)4,5. These parameters derived from the histogram analysis. SD, kurtosis and skewness describe the shape of the histogram representing the gray-level variation, asymmetry and peak with in the ROI, respectively. Entropy is a measure of texture irregularity and defined by following equation where l is the pixel levels (between l = 1 to k) in ROI, and p(l) is the probability of the occurrence of that pixel level.
The relationship between the CT texture features and 5-year overall survival (OS) or disease-free survival (DFS) was assessed by Kaplan-Meier analyses. Patients were dichotomized according to the best cut-off values, which were calculated by the log-rank test to classify the outcome of OS or DFS.
The multivariate Cox proportional hazard regression analysis using likelihood ratio test was performed to evaluate the independence of CT texture feature from clinical and pathological parameters, such as the Child-Pugh score, serum levels of AFP, tumor histological differentiation, pStage, or presence/absence of venous invasion, on OS or DFS. The correlations between CT texture features and clinical parameters were assessed: presence of HBs-Ag or HCV-Ab using paired-t-test; Child-Pugh classification, differentiation of HCC, or TNM stage using Spearman’s correlation; AFP using Pearson’s correlation. A P-value < 0.05 was considered to indicate statistical significance. The correlations of the CT texture features were assessed using the Spearman rank correlation coefficient. All statistical analyses were performed by a statistician using SPSS (version 18.0; IBM Corp, Armonk, NY, USA) statistical software.
Of the total number of patients, 12 were stage I, 60 stage II, 35 stage III, 12 stage IVA, and 3 stage IVB (Table 1). A total of 60 patients had a Child-Pugh score of 5, 46 a score of 6, 13 a score of 7, 2 a score of 8, and 1 a score of 9.Patients with the Child-Pugh score of 8–9 had inherited constitutional jaundice. AFP levels above 400 ng/ml were observed in 22 (18%) patients and venous invasion in 38 (31%) patients. A total of 77 (63.1%) patients had a single HCC, 24 (19.7%) patients had two HCCs, and 21 (17.2%) patients had three or more HCCs. The histopathologic differentiation of the tumors was well in 26 (21.3%) patients, moderate in 82 (67.2%) patients, poor in 10 (8.2%) patients and necrosis in 4 (3.3%) patients. Sixty-two (50.8%) patients died within 5 years, and the median OS was 50.6 months. A total of 98 (80.3%) patients suffered HCC relapse within 5 years, and the median DFS was 19.8 months.
The Kaplan-Meier curves for OS or DFS was significantly different between patient groups dichotomized by cut-off values for all CT texture parameters with filtration at at least single filter level (Table 2, Fig. 2). For SD and entropy, OS was significantly different without filtration or with filtration at each filter level (fine, medium and coarse). OS was different without filtration or with filtration (medium) for skewness. For mean, MPP and kurtosis, OS was not different without filtration between patient groups, but it was different with filtration (fine, fine and medium and coarse, and fine and medium, respectively). DFS was significantly different for skewness and kurtosis without filtration or with filtration (fine and medium, fine and coarse, respectively). DFS was different for mean, SD, entropy and MPP with filtration at single filter level (fine, coarse, coarse, and coarse, respectively).
Among the CT texture features and clinical and pathological parameters for OS, multivariate Cox proportional hazard regression analysis using likelihood test showed a significant effect of SD, entropy and MPP without filtration or with filtration at each filter level. For mean and skewness, a significant effect was found without filtration or filtration at single filter lever (fine, medium, respectively) (Table 2). SD, entropy and MPP at coarse filter, and skewness without filtration showed a significant effect for DFS. Thus, the possibility of CT texture feature increase the prognostic prediction of HCC by clinical and pathological information was shown. The correlation between CT texture features and clinical parameters (presence of HBs- Ag or HCV-Ab, Child-Pugh classification, differentiation of HCC, TNM stage or AFP) was not found (Table 3).
Strong positive correlations were found between the SD and entropy at each filter level and without filtration, and the mean and MPP, mean and SD, SD and MPP, and entropy and MPP at each filtration (Fig. 3). The skewness or kurtosis was not strongly related to mean, SD, entropy or MPP.
In this study, we assessed the relationship between non-contrast-enhanced CT texture features and the prognosis of HCC. When patients were dichotomized using the cut-off calculated using the log-rank test, the Kaplan-Meier curves for OS or DFS was significantly different for CT texture feature parameters. Using multivariate Cox proportional hazard regression analysis, CT texture features showed a significant prognostic effect for OS or DFS similar on the clinical parameters.
The heterogeneity of a tumor is thought to be associated with variations in genomic subtype, gene expression, angiogenic factors, and the tumoral microenviroment6. This may explain the difficulties encountered in the validation of oncology biomarkers owing to sampling bias and help predict therapeutic resistance7. The assessment of heterogeneity has been applied to classifying pulmonary nodules8,9,10,11, assessing disease severity12,13,14, and determining the prognoses of colorectal carcinoma15, non-small cell lung cancer16, and esophageal carcinoma17.
We assessed the heterogeneity of HCC using texture analysis to determine the prognosis of HCC. Similar to previous reports15,16,17, we determined a relationship between the prognosis of HCC and the texture features of preoperative non-contrast-enhanced CT images. To reduce photon noise, which may mask any underlying biological heterogeneity, a Laplacian of Gaussian filter using an edge-detecting filter with Gaussian smoothing to remove high-frequency image photon noise was adopted. This filter was also adopted in the previous studies of tumor prognosis using CT15,16,17. The Kaplan-Meier curves for OS or DFS was significantly different between patient groups dichotomized by cut-off values for all CT texture parameters with filtration at at least single filter level. Without filtration, OS or DFS was not significantly different for some CT texture parameters. These results may show the CT texture analysis with filtration delineated the biological heterogeneity which related to the prognosis more in detail than without filtration. Texture analysis with filtration may be useful for developing preoperative strategies for HCC and personalized medicine.
SD, entropy, and MPP showed similar results; OS was different in the Kaplan-Meier curves at each filter level, and DFS was different at the coarse filter. The Spearman rank correlation coefficient revealed a high positive correlation among these CT texture parameters. On the other hand, skewness and kurtosis showed different trends from SD, entropy, and MPP. The skewness and kurtosis are a higher order of statics, whereas SD, entropy, and MPP are first or second order statics4. This may explain that skewness and kurtosis behave differently from other texture features. We applied three levels for filtration. Among them, coarse filter of SD, entropy and MPP showed significant difference among dichotomized patient groups for both of OS and DFS. Focusing on SD, entropy, and MPP, the coarse filter was useful to predict both of OS and DFS of HCC patients.
For SD, entropy, and MPP, DFS did not differ between dichotomized patients groups based on cut-off values at lower filters (fine and medium) nor without filtration. This may show that the effect of photon noise to obscure heterogeneity at lower filter level or without filtration is more prominent than at coarse filter. Meanwhile, OS was different at each filter level. The reason to explain this difference between DFS and OS is not clear; however, the higher ratio (80.3%) for HCC relapse than the ratio of deaths (50.8%) within 5 years could be one of the causes. The incidence of recurrence remains after hepatectomy, with a 5-year recurrence rate of approximately 80%18,19,20, and the data of this study are consistent with previous reports.
The clinical or pathological prognostic factors for HCC were the Child-Pugh score, AFP serum levels, pathological differentiation, pStage, or presence/absence of venous invasion21,22,23. The independence of texture features with respect to these parameters was evaluated by multivariate Cox proportional hazard regression analysis. We found independence of most CT texture parameters except kurtosis on OS, and a few CT texture parameters (SD, entropy and MPP with coarse filter, and skewness at without filtration) showed independence on DFS. These results were similar to results of Log-rank test. This indicates CT texture feature may increase the prognostic prediction of HCC by clinical and pathological information. For the prevention of HCC recurrence, several therapies, such as transarterial radioactive iodine, adoptive immunotherapy, use of retinoid, interferon and vitamin K2, are reported24. The information of the prognostic prediction of HCC is considered useful for the selection of patients undergoing these prophylactic therapies. Also, frequent examination for recurrence check may be performed on patients at high risk of recurrence.
We demonstrated that texture analysis using non-contrast-enhanced CT images show a relationship with HCC prognosis. CT assessment of preoperative HCC requires a contrast agent, except in cases with contraindications, such as an allergy to the iodinated contrast agent, asthma, or renal failure. Usually, the information gained from non-contrast-enhanced CT imaging for HCC is limited, and thus enhanced CT images are regarded as more useful for forming preoperative strategies for HCC. In this study, we showed the significance of non-contrast-enhanced CT images in…