Abstract
Background
Liver metastases limit survival in colorectal cancer. Earlier detection of (occult) metastatic disease may benefit treatment and survival.
Objective
The objective of this article is to evaluate the potential of whole-liver CT texture analysis of apparently disease-free liver parenchyma for discriminating between colorectal cancer (CRC) patients with and without hepatic metastases.
Methods
The primary staging CT examinations of 29 CRC patients were retrospectively analysed. Patients were divided into three groups: patients without liver metastases (
Results
Mean E1.5 and E2.5 for the whole liver in patients with synchronous metastases were significantly higher compared with the non-metastatic patients (
Conclusion
Texture analysis of the apparently non-diseased liver holds promise to differentiate between CRC patients with and without metastatic liver disease. Further research is required to determine whether these findings may be used to benefit the prediction of metachronous liver disease.
Introduction
The liver is the most common and often unique site of metastasis in patients presenting with colorectal cancer (CRC).1–4 Studies have reported that >20% of CRC patients present with synchronous metastasis at primary staging. 5 Furthermore, cumulative rates of metachronous liver metastases of 4.3% have been reported at one-year follow-up, increasing up to 14.5% at five-year follow-up. 3 It has been hypothesised that these metachronous liver metastases may in fact develop from occult liver metastases that are already present at the time of primary staging.4,6 The presence of occult liver metastases (as established at histopathological examination after liver surgery) is associated with an unfavourable prognostic outcome regarding overall survival, 7 and addition of chemotherapy may reduce the presence of micrometastases. 8 With conventional staging tools it has so far not been feasible to establish the presence of micrometastases within the liver. This could, however, be of clinical importance as it may help establish a more risk-adapted (neo)adjuvant treatment strategy.
Previous studies have mainly focused on methods targeting liver perfusion to detect occult liver metastases, based on the hypothesis that changes in the hemodynamics of the liver occur before focal metastases become visible on morphological imaging. Some promising results have been reported for measuring the Doppler Perfusion Index (DPI) on ultrasound9,10 or hepatic perfusion indices (HPI) and hepatic transit time (HTT) on computed tomography (CT) or magnetic resonance imaging (MRI).11–18 Challenges with these techniques, however, are that they are operator dependent (DPI)19,20 and require a relatively complex image acquisition and post-processing algorithm.
An alternative approach would be to evaluate the spatial heterogeneity of the liver by means of texture analyses, which can be readily obtained from routinely acquired clinical CT data without the need for extra imaging procedures. Texture analysis focuses on the distribution and relationships of grey-level values within CT images of the liver. Similar to overall changes in perfusion, it can be expected that the texture of the liver parenchyma will change when the liver becomes invaded with tumoural cells. Apart from morphologically visible metastases, the presence of tumour may also affect the spatial heterogeneity of the remaining – apparently non-diseased – part of the liver.21,22 Previous studies in CRC patients without metastases have demonstrated that texture parameters (entropy and uniformity) of the normal liver parenchyma were significantly correlated with the hepatic hemodynamic status and glucose metabolism23,24 and may also be used as a predictor of survival. 25 Furthermore, results from a study by Ganeshan et al. suggest that – in patients with known CRC liver metastases – texture parameters may be more sensitive to tumour-related changes in the apparently disease-free areas of the liver parenchyma than measurements of hepatic attenuation or perfusion. 22
In order to determine whether texture analyses may be used to assess occult metastatic disease and predict the development of metachronous liver disease, a first step should be to investigate whether there are any changes in the overall liver CT texture in patients with known metastases. A next step would be to see whether such changes may already be observed in patients who will develop metastases at a later stage in the course of their disease. Hence, the aim of our study was to assess whether any differences can be detected in the CT texture of the apparently disease-free liver parenchyma of patients without and patients with hepatic metastasis from CRC.
Materials and methods
Patients
This study retrospectively analysed 29 patients who were diagnosed with colon or rectal cancer at our institution between March 2007 and August 2011. All patients routinely underwent contrast-enhanced abdominal CT as a part of their primary staging work-up. Three groups of patients were considered: Group A (
CT acquisition
Contrast-enhanced CT including a portal venous phase scan of the liver was performed as part of the routine liver/abdomen CT imaging protocol for CRC at our institution using multi-slice CT equipment (Siemens Somatom Sensation 16 or Somatom Definition Flash, Siemens Healthcare, Erlangen, Germany; Toshiba Aquilion 16, Toshiba Medical Systems, Tokyo, Japan; Philips Brilliance 64 or Philips Gemini TF 64, Philips Medical Systems, Best, The Netherlands). The portal venous phase CT images used for texture analysis were obtained with a tube voltage of 120 kVp. The in-plane pixel resolution for the liver images varied from 0.55 to 0.85 mm2, with a slice thickness of 3 mm.
Image assessment
Images were transferred from the hospital’s picture archiving and communication system (PACS) to an offline workstation for texture analysis. Regions of interest (ROIs) covering the whole-liver volume were manually drawn on the portal phase CT images by a single reader (WvO). ROIs were drawn to cover the whole surface of the liver on each consecutive slice, excluding the border of the liver (to avoid partial volume effects), any visible lesions (metastases or benign focal liver lesions such as haemangiomas, or cysts), major portal and hepatic veins and the inferior vena cava. Because the border of the vena cava is often not clearly defined, we chose not to include the caudate lobe of the liver for the ROI delineation and texture analyses. A representative example of the manual ROI delineation is shown in Figure 1.
Representative example of regions of interest drawn on a portal venous phase computed tomography (CT) image of the liver (a). Liver borders, any visible focal liver lesions, major vessels and the caudate lobe were manually excluded. The images represent the different types of filtrations: no filtration (a), fine filtration (b), medium filtration (c) and coarse filtration (d).
Texture analysis
Texture analyses were performed using a dedicated script written in MATLAB (MathWorks Inc, Natick, MA, USA) by one of the authors (RSS). The mathematical analysis technique used for this study is described in detail in the appendix. It comprised two main stages: (a) image filtration, followed by (b) quantification of texture. For image filtration, a Laplacian of Gaussian (LoG) band-pass filter was applied to the ROI using sigma (σ) values of 0.5 (fine-scale filtration), 1.5 (medium-scale filtration) and 2.5 (coarse-scale filtration). Subsequently, the texture of the liver parenchyma was characterized by mean grey-level intensity (M), entropy (E) and uniformity (U).
Statistical analyses
Statistical analyses were performed using the Statistical Package for the Social Sciences (SPSS Inc, version 16.0, Chicago, IL, USA). Texture parameters were compared between group A (no metastases) and group B (synchronous metastases), and between group A and group C (metachronous metastases) using a Student’s
Results
Patient characteristics
Baseline characteristics of patients with colorectal liver metastasis
Group A: patients without liver metastasis; Group B: patients with synchronous liver metastasis at the time of primary staging; Group C: patients with metachronous liver metastasis within 18 months after primary staging. N/A: not applicable.
Median follow-up from the date of the primary staging computed tomography (CT).
Median interval between the primary staging CT and the detection of metachronous metastatic disease.
D: the maximal axial diameter of the largest metastasis.
Texture analyses
Mean (±standard deviation) of CT texture parameters of the apparently non-diseased liver parenchyma measured at different scale values (σ) on portal venous phase CT
Group A: patients without liver metastasis; Group B: patients with synchronous liver metastases at the time of primary staging; Group C: patients with metachronous liver metastases within 18 months following primary staging; Significant results are printed in bold. CT: computed tomography.
Diagnostic performance for assessment of synchronous liver metastases
Diagnostic performance for entropy and uniformity with or without filtrations in assessing synchronous colorectal liver metastasis by receiver operating characteristics (ROC) curve analyses
Data are percentages; absolute numbers are given in parentheses. AUC: area under the ROC curve; PPV: positive predictive value; NPV: negative predictive value; CI: confidence interval.
Discussion
The goal of this study was to assess whether any differences can be detected in the CT texture of the apparently disease-free liver parenchyma of patients without and patients with hepatic metastasis from CRC. Our results indicate that entropy and uniformity with medium to coarse filtrations were useful for detecting differences in spatial heterogeneity between the normal-appearing liver parenchyma of patients with and without liver metastases at the time of primary staging. AUCs for assessment of the presence of synchronous liver metastasis using these texture parameters ranged between 0.73 and 0.78.
Apparently, CT texture analysis has the potential to detect changes in the overall structure of the liver at a level that goes beyond the resolution of morphological CT imaging. Higher entropy and corresponding lower uniformity were found in the non-diseased part of the liver (i.e. excluding the metastatic lesions) of patients with synchronous metastatic disease as compared to those without. In theory, entropy is a measure of texture irregularity. It is almost the opposite of uniformity, which is a measure of the uniformity of grey-level distribution indicating tissue homogeneity. The parenchymal structure of the non-diseased part of the liver in patients with known liver metastases appears to be more heterogeneous than that of a patient without any known liver metastases. There are several hypotheses that may explain this phenomenon. First, in patients with proven focal metastatic lesions, additional micrometastases may be present in the apparently disease-free remaining liver, resulting in a more heterogeneous overall liver structure. Previous studies reported that additional intrahepatic micrometastases (i.e. vascular, biliary, lymphatic infiltration and satellite nodules) were present within the surgically resected specimens in 58%–81% of CRC patients who had not undergone neoadjuvant chemotherapy.8,26,27 Second, the texture measures of the liver parenchyma may be influenced by tumour-induced changes in liver perfusion since different pre-clinical and clinical studies have suggested that measurements of liver texture on CT may be associated with liver vascularity.23,24,28 Interestingly, our results also showed a subtle trend – similar to the synchronous metastases group – towards increased entropy and decreased uniformity values in the patients who had no metastases, but developed metachronous liver disease within two years after the primary staging CT on which the texture analyses were performed. Although results were not significant and our metachronous group is obviously too small to draw any meaningful conclusions, this is a potentially interesting finding that warrants further research as it may indicate that tumoural changes – whether due to the presence of micrometastases or perfusional changes – are already occurring in the liver preceding any visible morphological changes.
In contrast with our results, a previous study by Ganeshan et al. reported lower entropy and higher uniformity values in the liver parenchyma of patients with liver metastasis. 22 A contributing factor to these conflicting results could be that in this previous report a different imaging protocol was used, including both arterial and portal phase CT data which were derived from a dynamic scan protocol (consequently using lower contrast doses and higher injection rates). As a result, the portal venous phase will occur earlier than in a routine multi-phase CT imaging protocol as used in our study. In their study, significant results were mainly found in the arterial phase, while no significant results could be demonstrated in the portal venous phase (contradictory to our findings). Ganeshan reported that different contrast phases resulted in variations in entropy and uniformity, which may contribute to these inconsistent findings. There were also some differences in the study groups: Ganeshan et al. included patients during surveillance after primary CRC resection, while our study included CRC patients at the time of primary staging (prior to any treatment). This might result in different effects on hepatic haemodynamics caused by the presence vs. absence of the primary tumour.11,24 Moreover, in the study by Ganeshan et al. it was unknown from the archival study dataset whether the study patients had been exposed to any chemotherapeutic agents (which was known not to be the case in any of our study patients), which may have had a considerable effect on the texture measurements. Finally, in our study the whole volume of the liver was assessed rather than a single mid-liver axial section, as was the case in their previous study. Whole-volume measurements may be more representative of tissue heterogeneity and result in smaller measurement variations compared with single-slice analyses. On the other hand, the texture of specific hepatic segments in which liver metastases are located might be different from that of other segments. These subtle changes in more focal areas of the liver may be averaged out by choosing a whole-volume approach, an issue that requires further investigation. Altogether, larger study cohorts with uniform CT protocols and study populations are required to gain better insight into the complex relationship between texture measurements, treatment effects and underlying micro-architectural changes in the liver parenchyma.
Our study is limited by its retrospective nature and the single-reader design. Although we believe that interobserver variations in whole-volume delineation will likely be limited as has also previously been demonstrated in other cancer types, this is an issue that will need to be addressed by future studies. 29 Also, the total number of patients – particularly for the metachronous liver metastasis – is small. Further research is required to validate the value of CT texture analyses and particularly to prospectively investigate its potential for the prediction of the development of metachronous disease.
Clinical impact and conclusions
In CRC patients, the presence of intrahepatic micrometastases has been demonstrated to be associated with poorer prognosis. 7 Previous studies have shown that preoperative neoadjuvant chemotherapy may reduce the incidence of intrahepatic micrometastases in patients with colorectal liver metastasis. 8 Prediction of the presence of micrometastases by non-invasive imaging methods may therefore potentially impact the therapeutic strategy. For the patients with colorectal liver metastases who are scheduled to undergo hepatectomy, preoperative identification of those patients who would likely have micrometastases in addition to the visible liver metastases might allow a better stratification for preoperative neoadjuvant chemotherapy. For CRC patients without overt liver metastasis, accurate screening of micrometastases at the time of primary staging may allow personalized treatment strategies, offering chemotherapy for the high-risk patients, while omitting it for the lower-risk patients. Conventional imaging methods are incapable of detecting occult liver metastasis. CT texture analysis is a promising noninvasive approach that allows quantification of the spatial heterogeneity of the liver parenchyma. Texture data can be readily obtained from routine CT scans without the need for additional imaging procedures. Our results suggest that texture analysis holds promise to detect subtle changes in the ‘normal’ liver parenchyma of patients with known colorectal liver metastases that might reflect the presence of tumour-induced structural and/or haemodynamic changes. Furthermore, a similar subtle trend might be observed in patients who develop metachronous liver lesions later during the course of their disease, although our current results on this matter are too preliminary to draw any meaningful conclusions. These results warrant further investigation to evaluate whether CT texture analysis will have potential value for the detection of occult liver metastases and predict which patients are at risk to develop metachronous liver disease.
Footnotes
Funding
This research received no specific grant from any funding agency in the public, commercial, or not-for-profit sectors.
Conflict of interest
None declared.
