1. Introduction
Pancreatic cancer is among the top cancers with high mortality affecting over 200,000 deaths every year, worldwide [
1,
2]. Unfortunately, a large proportion of patients are identified at an advanced stage with poor prognosis, especially for metastatic pancreatic cancer, for which the liver is the most common site, accounting for 37–41.9% of the initially diagnosed cases, with a 5-year overall survival (OS) rate of 2%, and a median life expectancy of less than 1 year [
3,
4]. FOLFIRINOX and gemcitabine-based chemotherapy regimens are common treatments for metastatic pancreatic cancer. It was observed that FOLFIRINOX was superior to gemcitabine alone in progression-free survival (PFS), response, and OS in the treatment of metastatic pancreatic cancer, and another phase 3 trial confirmed the superiority of gemcitabine plus nanoparticle albumin-bound paclitaxel (nab-paclitaxel) as well [
5,
6]. Assessment of treatment response to chemotherapy is important, especially for the non-responders, which might provide these patients opportunities to find more appropriate treatment plans in time.
The response evaluation criteria in solid tumors (RECIST) as a frequently used tool for the assessment of tumor response only takes size change of the lesions into account [
7]. Choi et al. developed comprehensive criteria incorporating changes in both tumor size and attenuation, which defines partial response (PR) as a decrease of >10% in tumor size or a decrease of >15% in tumor attenuation. These criteria define progressive disease (PD) as a tumor size increase of >10% without meeting the PR criteria [
8]. The definition of PR according to the modified Choi criteria is a 15% reduction in enhancement and a 10% reduction in size [
9]. No studies have been conducted to assess the response to chemotherapy in pancreatic cancer with liver metastasis using these criteria.
Analysis of tumor lesion heterogeneity will reveal vital information concerning response to treatment [
10,
11]. CT texture analysis is an emerging technique used to process images and hence help to characterize lesion heterogeneity. This technique analyzes the relationship and distribution of pixel gray levels within a lesion and reveals spatial variations among individual gray patterns or levels [
10,
12,
13]. Extracted features include kurtosis, skewness, mean of positive pixels (MPP), entropy, standard deviation (SD), and mean gray-level intensity (MI) [
14]. It has been proved that CT texture analysis played an important role in the assessment of tumor response to various treatments and had the potential to establish a more precise assessment criteria for response evaluation [
15,
16,
17,
18]. However, texture analysis as a non-invasive tool for the assessment of chemotherapy treatment response in pancreatic cancer with liver metastasis has not been investigated before.
In this study, we assessed alterations in CT texture of metastatic liver lesions following gemcitabine-based chemotherapy in pancreatic cancer patients and then compared the effectiveness of this approach with existing assessment criteria in evaluating treatment response and determining time to progression (TTP).
2. Materials and Methods
2.1. Study Population
Prior to the study, approval of this study was waived by the institutional review board of Peking Union Medical College Hospital for its retrospective nature. Patients with pancreatic cancer with liver metastasis undergoing chemotherapy between September 2014 and October 2018 were identified from our institutional electronic medical database. Inclusion criteria: (a) They were pathologically diagnosed to be pancreatic cancer with liver metastasis; (b) they were chemotherapy naïve and received gemcitabine-based chemotherapy (gemcitabine plus an oral fluoropyrimidine anticancer agent, S-1) as first-line treatment; (c) baseline contrast-enhanced CT had been done in 2 weeks prior to treatment initiation; (d) contrast-enhanced CT had been done following treatment to monitor response. A total of 168 patients met the inclusion criteria. Patients were excluded if they had previously undergone radiofrequency or microwave ablation of the liver metastatic lesions (n = 17), had infiltrative HCC without any accurately delineable lesion (n = 20), died during the first 15 days after treatment initiation (n = 3), and were without a regular documented follow-up after chemotherapy initiation or until disease progression (n = 18).
2.2. Follow-Up and Endpoints
The follow-up (radiological, clinical, and biological) was carried out every 2 months as outlined by the institutional protocol. During follow-up, radiological assessment was performed with contrast-enhanced CT scans of the pelvis, abdomen, and thorax. TTP was defined as the time from the start of chemotherapy to radiologic progression, which was the chosen end point detailed by RECIST 1.1. Participants who did not exhibit radiologic progression at the end of the follow-up were censored.
2.3. Protocol of CT Examination
All images were acquired in the Department of Radiology at our hospital. All the scans were done on 128-detector CT scanners (Siemens SOMATOM Definition Flash, Siemens Healthcare, Forchheim, Germany). The scanning parameters were as follows: tube voltage, 120 kVp; tube current, 150 mAs (thorax) or 200 mAs (abdomen and pelvis) with dose modulation; gantry rotation time, 0.5 s; table increment 46.8 mm per rotation; matrix 512 × 512. Images were routinely reconstructed with 5.0 mm slice thickness and 5.0 mm intervals. Non-ionic contrast media (Ultravist, 370 mg of iodine per milliliter, Schering, Berlin, Germany) were injected with 1.5 mL per kilogram of body weight at a rate of 3.0 mL/s using an automatic power injector. Bolus tracking was applied, and the thorax and upper abdomen (from the supraclavicular fossa to the iliac crest) were imaged in the arterial phase of enhancement (25-s delay), and the abdomen and pelvis (from the dome of the diaphragm to the pubis) were imaged in the portal venous phase of enhancement (70-s delay), according to the protocol in our institution.
2.4. Texture Analysis
All images were analyzed by two radiologists (20 years and 5 years of abdominal imaging experience, respectively) who were blinded to the clinical outcome of the participants. Measurements were made on metastatic liver lesions according to the protocol provided by previous study [
10], and tumor dimension and attenuation were measured using an image viewing software (Centricity, GE Medical Systems, Milwaukee, WI, USA). Texture analysis was carried out on the lesions using a commercial software (TexRAD, TexRAD Ltd., Cambridge, UK) by drawing a region of interest (ROI) around the peripheral margin of the lesions on CT images showing the largest area of the metastatic lesions. Various texture parameters were quantified prior and post to the application of spatial band-pass filters. Filter values on a scale of 0 to 6 were adjusted to select an appropriate scale. Where 0 = absence of filtration (spatial scaling factor (SSF) = 0), 2 = fine texture, 3 and 4 indicate different degrees of medium texture and 5 and 6 indicate different degrees of coarse texture. This analysis yielded derived CT texture images showing imaging characteristics at various spatial scales in the ROI (
Figure 1).
2.5. Statistical Analysis
Quantitative data were shown as mean ± SD or medians. Categorical variables were presented as percentages. Mann-Whitney U test, Fisher exact, or Chi-square test were utilized appropriately to compare proportions and ranks of variables between training and validation cohorts. Texture parameters extracted from liver metastatic sites at baseline CT and at CT after two cycles of treatment, and the percentage change from the baseline value were all statistically analyzed. The association between texture parameters, RECIST 1.1, Choi, and modified Choi response criteria after two cycles of treatment and measured TTP were evaluated with the Kaplan-Meier analysis. Receiver operating characteristic (ROC) analysis was performed to identify threshold values for texture parameters. Kaplan-Meier analysis of participants with values below or above the thresholds were done to show the proportion of patients who did not have disease progression at any given time. Nonparametric log-rank test was performed to examine the differences between Kaplan-Meier curves for texture parameters, RECIST 1.1, Choi, and modified Choi criteria. Independent predictors of TTP among clinical and texture parameters were identified by performing multivariate analysis. All statistical analysis was done on SPSS 22.0 (IBM, New York, NY, USA).
p < 0.05 indicated a significant difference, however, when facing multiple comparisons, a stepwise Holm–Bonferroni procedure was performed to reduce the potential for type I errors [
19].
3. Results
In this study, we finally analyzed 110 pancreatic cancer patients with liver metastasis. These patients were randomly divided into two cohorts with 60 participants (mean age of 61.1 ± 8.6 years; range, 42–76 years) constituted a training cohort; while 50 participants (mean age, 56.1 ± 10.3 years; range, 31–73 years) constituted the validation cohort. The demographics and clinical features did not show significant difference between the validation and training cohorts (
Table 1). In the training cohort, 10% of the patients (6/60) exhibited more than three liver target lesions, while in the validation cohort, 24% of patients (12/50) had >three target lesions. The initial treatment responses in the training cohort were shown in
Table 2, categorized using RECIST 1.1, Choi, and modified Choi criteria. Follow-up imaging revealed that 42 of 60 patients in the training cohort and 32 of 50 patients in the validation cohort showed progression at the end of the follow-up. The median TTP was 178 (range, 36–1441) days in the validation cohort and 241 (range, 43–1164) days in the training cohort.
The median values with their ranges for texture parameters for absolute scale values at baseline and after two treatment cycles, as well as corresponding percentage change, were shown in
Table 3.
Predictive ability of texture parameters at different SSF were shown in
Table 4,
Table 5 and
Table 6. With fine texture (SSF = 2), baseline MI and percentage change in both SD and kurtosis, were found to be significantly correlated with TTP. With medium texture (SSF = 4), baseline kurtosis, and percentage change in MPP, were found to be significant predictors. Baseline kurtosis at SSF3, as well as percentage changes in entropy and MPP at SSF5, were also found to be significantly associated with TTP. Kaplan-Meier analysis of the proportion of patients without disease progression revealed significant differences for the above-mentioned texture parameters, which were better than those for RECIST 1.1, Choi-, and modified Choi-defined responses after chemotherapy (
p value < 0.05 vs.
p value = 0.398,
p value = 0.142, and
p value = 0.536, respectively) (
Figure 2).
After the performance of multivariate Cox regression analysis, baseline MI at fine texture scale (SSF2,
p value = 0.028), baseline kurtosis at medium texture scales (SSF3,
p value = 0.032; SSF4,
p value = 0.049), percentage change in kurtosis and SD at fine texture scale (SSF2,
p value = 0.033 and
p value = 0.016, respectively), and percentage change in MPP at medium texture scale (SSF4,
p value = 0.016) were found to be independent prognostic factors for TTP (
Table 7).
However, these initially obtained results from training cohort needed further validation in an independent cohort. Optimal thresholds of above-mentioned independent texture parameters were acquired using ROC, and they were used to dichotomize texture parameters in the validation cohort. After application of Kaplan-Meier analysis, percentage change in SD at fine (SSF2) texture scale was confirmed as the only significant predictor of TTP (
p value = 0.019) (
Table 8). Percentage change in SD exceeding −24% at SSF2 was significantly associated with longer TTP following gemcitabine-based chemotherapy (
Figure 3).
4. Discussion
In this study, we found that texture features of liver metastatic lesions on contrast-enhanced CT images could accurately reflect treatment response and predict TTP in pancreatic cancer patients treated with chemotherapy, and these texture parameters showed great advantages over traditional assessment criteria. Importantly, we demonstrated that the percentage change in SD at fine texture scale had the ability to serve as an independent predictive biomarker, and its correlated optimal threshold value was validated in an independent cohort of patients.
RECIST 1.1 as a commonly used imaging-defined response assessment tool, has been widely used in the evaluation of various cancers. However, accurate assessment of clinical benefits from chemotherapy using RECIST 1.1 was constrained by the fact that many patients might exhibit tumor shrinkage to a degree that did not meet the criteria of RECIST 1.1, but still exhibited significant prolonged survival [
10]. Tumor necrosis might occur in response to therapy without appreciable effects on tumor size, despite reduced tumor vascularization and tumor attenuation [
10,
20,
21]. Therefore, response evaluation by RECIST 1.1 may significantly underestimate the number of patients experiencing treatment benefits [
7]. Tumor attenuation was previously used to assess treatment response in other cancers [
22]. Choi and modified Choi criteria, which included changes in tumor enhancement, were shown to be superior to RECIST 1.1 in their capacity to predict drug efficacy for various cancers [
8,
9,
23,
24,
25,
26,
27,
28]. Nevertheless, all these assessment tools did not reveal significant differences in Kaplan-Meier curves of the proportion of patients without disease progression, suggesting that existing assessment criteria were not sufficient to reflect underlying alterations in tumor heterogeneity. Therefore, efforts should be made to explore new assessment tools or biomarkers for the identification of pancreatic cancer patients with liver metastasis who might benefit from chemotherapy. Such biomarkers will inform the use of alternative therapies on possible non-responding patients earlier and help better understand disease heterogeneity.
Heterogeneity of neovascularization reflects conditions of hypoxia inside tumors, which is closely linked to high risk of invasion, metastasis, immunosuppression, and unfavorable response to chemotherapy [
29]. Tumor response to a certain treatment could be obviously affected by its heterogenous blood supply. In this case, poor vascularization compromises drug availability at tumor sites. The emerging technique of CT texture analysis permits quantitative extraction of tumor heterogeneity, which enables assessment of tumor neoangiogenesis and hypoxia. Multiple previous studies have characterized the relationship between CT texture parameters and histopathological indicators of hypoxia or angiogenesis in various cancers [
30,
31]. CT texture analysis may also predict survival and response to treatment in cancer patients [
10,
32]. To our knowledge, the use of texture parameters as tumor response assessment criteria and predictive biomarkers in pancreatic cancer patients with liver metastasis undergoing chemotherapy has not been explored before.
Our results demonstrated the independent predictive significance of percentage change in liver metastatic lesion SD for TTP in pancreatic cancer patients with liver metastasis. The identified threshold values found by the training cohort were validated by another independent cohort, which increased the generalizability of the present findings [
33]. SD is a measure of data dispersion from the mean and it increases in proportion to the square root of the number of objects highlighted and their mean intensity difference relative to background (i.e., dark and bright objects are both positive) [
34]. Different from RECIST 1.1, Choi, and modified Choi criteria, the percentage change in texture SD (SSF = 2) of liver metastatic lesions was an independent predictor of TTP. By applying a threshold of −24% for SD at a fine scale (SSF = 2), the performance of the obtained Kaplan-Meier curves of proportion of patients without disease progression was notably different from and better than those obtained by RECIST 1.1, Choi, and modified Choi criteria. This revealed that SD as one of the indicators of tumor heterogeneity provided more details in the level of pixels or gray scales other than traditional enhancement and size change, thereby enabling complementary evaluation. This observation was supported by multiple recent studies involving unresectable pancreatic cancer and pancreatic head cancer [
28]. Heterogeneity at fine texture scale was proposed to reflect lesion vascular permeability [
35], which suggested that changes in vascular permeability caused by chemotherapy might explain the survival outcomes of pancreatic cancer patients with liver metastasis we observed in this study.
Due to its retrospective nature, this study was constrained by the limited number of participants, thus a prospective study with lager population needed to be conducted in the future to verify our conclusion. Therefore, it was still uncertain that CT texture of liver metastases could modify the clinical management of such patients, especially when multiple liver lesions or large, necrotic lesions presented. Ahn et al. [
15] conducted a study to determine whether CT texture analysis of liver metastasis from colorectal cancer is predictive of therapeutic response after cytotoxic chemotherapy, and they found that the lower skewness and narrower SD showed good performance. Nakanishi et al. [
36] also conducted a study focusing on the texture analsyis of liver metastasis, and they aimed to develop a radiomics-based prediction model for the response of colorectal liver metastases to oxaliplatin-based chemotherapy. Overall, 126 liver metastases were analyzed and multivariate analysis revealed that high radiomics scores was independently associated with good response. These above-mentioned studies might give us some inspiration for CT texture of liver metastases in response evaluation of pancreatic cancer. Another limitation was that only the largest cross-sectional area of the lesions, not the whole lesions, were analyzed, which therefore might not reflect the characteristics of liver metastatic lesions comprehensively. A 3D texture analysis software should be applied to evaluate the whole lesion heterogeneity and thus reproducibility of the study might be improved. Finally, selection bias could not be avoided as infiltrative metastatic lesions whose boundaries could not be clearly defined, were excluded from analysis.