Multi-level evidence that circulating CK18 is a biomarker of tumour burden in colorectal cancer.

BACKGROUND: Circulating total cytokeratin 18 (tCK18) and/or caspase cleaved cytokeratin 18 (cCK18) (measured by M65 and M30 enzyme-linked immunosorbent assays (ELISAs), respectively) are used as pharmacodynamic (PD) biomarkers of epithelial cell death in clinical trials. Having validated these ELISAs, we assessed their utility in colorectal cancer (CRC).
METHODS: We applied the assays in several settings: 53 controls; 97 patients undergoing surgery and 74 patients with metastatic CRC undergoing chemotherapy (55 first line; 56 patients with repeated sampling through chemotherapy). Prognostic significance was evaluated using Kaplan-Meier life tables and Cox models; PD utility was assessed by analysis of repeated measures.
RESULTS: Median cCK18 and tCK18 levels were elevated in patients with cancer (both P=0.0001), and among cancer patients, there were increasing trends from early to advanced stages (both P(trends)=0.0001). Increasing tCK18 predicted for reduced survival after surgery with curative intent (adjusted hazard ratio (HR) for doubling in concentration 1.77, 95% CI: 1.04, 3.01) and after first-line chemotherapy in metastatic disease (adjusted HR per doubling in concentration=1.78, 95% CI: 1.37, 2.30). In patients with progressive disease during chemotherapy, repeated sampling revealed profiles with high baselines and progressive upwardly increases after cycle 1.
CONCLUSION: This study provides evidence for cytokeratin 18 (CK18) as a prognostic and PD biomarker in patients with CRC and supports continued deployment of circulating CK18 in biomarker-enhanced trials.

Colorectal cancer (CRC) is the second most common cancer in Europe accounting for 183 000 new cancers per year. The integration of mechanism-based therapies (MBTs, e.g., cetuximab or bevacizumab) with chemotherapy in CRC has expanded treatment options, resulting in survival gains in some patients (Cunningham ). With the introduction of MBTs, there is a need to develop biomarkers to prioritise new agent development (predictive and prognostic biomarkers), optimise dosing schedules (pharmacodynamic (PD) biomarkers), and rationalise combination selection (Cummings ).

Circulating products of epithelial cell death have been used as PD biomarkers for MBT responses both in animal models (Olofsson ) and in clinical trials (Cummings ). The M30 apoptosense and M65 enzyme-linked immunosorbent assays (ELISAs) detect circulating cytokeratin 18 (CK18) fragments: the M30 ELISA selectively detects caspase cleaved CK18 (cCK18) (apoptosis), whereas the M65 assay detects cCK18 and uncleaved CK18 (total cell death, e.g., apoptosis and necrosis) (Cummings , 2006). We previously validated these assays (Cummings , 2006) and determined they are ‘fit for purpose’ (Greystoke ). To-date, we have reported their potential utility in testicular cancer (de Haas ), pancreatic cancer (Dive ), small cell lung cancer (SCLC) (Hou ; Dean ), and lymphoma (Greystoke ), and evaluated their performance in early clinical trials (Dean ; Gandhi ). A few studies with low sample numbers have evaluated serum cCK18 and total CK18 (tCK18) in patients with gastrointestinal cancers (Scott ) and CRC (Ausch , 2009b).

Previous studies in SCLC and testicular cancer suggested that elevated baseline tCK18 levels are associated with poor prognosis, and that in these chemo-sensitive tumours, response is accompanied by falling tCK18 levels (de Haas ; Hou ). In a colorectal xenograft model, tCK18 correlated with tumour burden, whereas measurement of cCK18 correlated with treatment-related tumour death (Cummings ). Based on these observations, we hypothesised that (i) elevated levels of tCK18 at baseline are associated with poor prognosis, (ii) circulating cCK18 levels are associated with tumour apoptosis following therapy, and (iii) circulating tCK18 correlates with tumour burden in CRC.

Patients and methods

Study design and patients

We used a mixed cross-sectional prospective design to assess the utility of serum cCK18 and tCK18 across multiple clinical settings in 224 subjects (Figure 1): (i) 53 individuals undergoing screening colonoscopy (controls); (ii) 97 undergoing surgery (Renehan ); (iii) 74 patients undergoing palliative chemotherapy. Survival analysis was restricted to 66 surgery patients who underwent curative resection and 55 patients undergoing first-line chemotherapy. In addition, 56 patients underwent repeated sampling throughout chemotherapy.

Figure 1

Flow diagram of analyses – a mixed cross-sectional prospective design. Abbreviations: CRC=colorectal cancer; FU=follow-up.

All treatments were carried out at the Christie NHS Foundation Trust, Manchester, UK. The study was carried out according to the Good Clinical Practice guidelines, and following local ethics committee approval.

Sampling and assays

Blood was collected in tubes containing a silica clot activator and centrifuged at 2000 g for 10 min to obtain serum. The sample collection, processing, storage, and measurement of cCK18 and tCK18 were as previously described (Greystoke ) (see Supplementary Material).

Baseline samples were collected 30 min prior to colonoscopy (group I), the day before surgery (group II), or the first day of chemotherapy (group III). Venous mesenteric sampling was performed intra-operatively prior to ligation of the pedicle drawing blood from the main draining vessel, for example, the inferior mesenteric vein for recto-sigmoid tumours. Repeated sampling in group III patients varied depending on clinical practice, but were typically collected on days 1, 3, 8, 15, 21, 28, 35, 42, 49, and 56, and thereafter determined by clinical indication.

Statistical analysis

CK18 biomarker distributions in patients with CRC were highly skewed (see Supplementary Table 1). Thus, comparisons were performed using Mann–Whitney U- and Kruskal–Wallis tests; paired samples were compared using Wilcoxan rank tests; and nonparametric tests for trend across ordered groups (Cuzick, 1985).

For prognostic analyses, we followed the REMARK guidelines (McShane ). Survival curves were derived using Kaplan–Meier analyses and comparisons performed using log-rank tests. Cox models were applied to assess whether cCK18 and tCK18 were independent prognostic factors for survival. We used a seven-step approach advocated by Hosmer and Lemeshow (2000), forcing cCK18 and tCK18 data separately into the final model. Biomarkers were handled as both continuous and categorical variables in separate models – the former, to test for independent prognostic significance; the latter, to allow extension of the univariate survival analyses and explore for potential clinically relevant cut points. For the continuous data, to reduce the right skewness of the distributions of the CK18 concentrations, the base 2 logarithmic transformation was used, which leads to a convenient interpretation: the hazard ratio (HR) associated with a change of one unit on the log2 scale corresponds to the HR associated with a doubling in CK18 level on the original scale. For all models, we excluded variables if the partial likelihood test P-value>0.05. In initial models for surgically treated patients, we considered age, sex, stage, and carcinoembryonic antigen (CEA: as tertiles); in initial models for patients with metastatic disease, we considered age, sex, presence of liver metastases, treatment regimens, CEA, alkaline phosphatise, aspartate transaminase, and lactate dehydrogenase as continuous variables (after natural log transformation). As these analytes correlated with each other and with cCK18 and tCK18 (see Supplementary Table 2), models were tested for co-linearity. Predictive accuracy and calibration were tested, respectively, using a concordance index (C-statistic; values ranging from 0.5 to 1.0) and the calibration statistic, analogous to the Hosmer–Lemeshow goodness-of-fit test (with values ⩽20 indicating good agreement).

In patients undergoing chemotherapy, we assessed PD response using analyses of repeated measures but for pragmatic purposes, we limited analyses to regression – joint mixed effects modelling may better test for differences between biomarkers, but the complexity of such modelling was beyond the scope of this paper (Steele, 2008). Instead, each patient was handled as a dummy variable to take account of intra-individual correlation. We investigated potential nonlinear relationships between the biomarkers and time on treatment using restricted cubic splines with three knots at 3, 14, and 21 days, chosen to represent immediate changes, and two subsequent time points at which future biomarker-directed trials may inform early decisions on continuation of trial drug. Carcinoembryonic antigen was checked at the start of every cycle – a 20% rise after cycle 2 was considered clinically relevant, while a rise in tCK18 >12% from cycle 2 onwards was considered a positive signal; larger than the documented intra-person variability in the assay (Cummings ). Response to therapy was assessed by CT scan every 2 months during therapy. Patients were categorised into progressive disease or clinical benefit (either stable disease or partial remission) based on RECIST 1.0 (Therasse ). All analyses were performed using STATA version 11.1 (College Station, TX, USA).

Results

Baseline characteristics

Baseline characteristics are shown in Table 1. The proportion of men was greater in subjects receiving chemotherapy (78.2%) compared with controls (54.7%) and surgically treated subjects (51.6% P=0.004); patients with cancer were older (P=0.0004) and had a lower BMI (P=0.02) than controls. However, in the control group serum, cCK18 and tCK18 were uninfluenced by gender, age, and BMI (see Supplementary Tables 3–5). As might be expected, median CEA levels were elevated in patients receiving palliative chemotherapy compared with surgically treated patients (P=0.0002).

Table 1

Baseline characteristics for groups I, II, and III

	Group I Normal colonoscopy	Group II CRC undergoing surgery	Group III First-line metastatic CRC
No. of patients	53	97	55
Men : women	29 : 24	50 : 47	43 : 12
Median age (IQR) years	60.2 (58.0–62.8)	68.1 (55.3–75.3)	65.0 (57.0–72.0)
Performance status
0 (%)	NA	65 (67)	25 (45)
1 (%)	NA	20 (21)	20 (36)
2 (%)	NA	12 (12)	10 (18)
Tumour stage a
I and II (%)	NA	53 (55)	NA
III and IV (%)	NA	44 (45)	NA
Distant metastases – no liver (%)	NA	NA	21 (38)
Distant metastases – liver (%)	NA	NA	34 (62)
Surgical treatment
Curative segmental resection (%)	NA	50 (52)b	NA
Curative – exenterative surgery (%)	NA	4 (4)c	NA
Curative segmental resection with metastases (%)	NA	12 (12)d	NA
Palliative (%)	NA	31 (32)	NA
Chemotherapy
Ox+(5-FU/Cap) (%)	NA	NA	28 (51)
Ir+(5-FU/Cap) (%)	NA	NA	6 (11)
Cap (%)	NA	NA	11 (20)
Others (%)	NA	NA	10 (18)
Median serum CEA (IQR) ng ml−1	NA	7 (3–51)	18.5 (3–95)

Abbreviations: Cap=capecitabine; CEA=carcinoembryonic antigen; CRC=colorectal cancer; FU=fluorouracil; Ir=irinotecan; IQR=inter-quartile range; NA=not applicable; NR=not recorded; Ox=oxaliplatin.

TNM staging 5th edition.

20 anterior resections, 10 abdomino-perineal resections, 2 Hartmann resections, 7 left hemi-colectomies, 6 right hemi-colectomies, 3 subtotal colectomies, and 2 local excisions.

Three total pelvic clearances and one cytoreductive surgery.

Three anterior resections; one abdomino-perineal resection; one Hartmann resection; four left hemi-colectomies; two right hemi-colectomies, and one total pelvic clearance.

Comparisons of cancer patients vs controls

Concentrations were significantly increased in patients with cancer for tCK18 (median (inter-quartile range): 410 (287, 942) vs 185 (151, 235) U l−1, P=0.0001) and cCK18 (162 (106, 297) vs 113 (168, 89) U l−1, P=0.0001) compared with those in controls (see Supplementary Table 6). Among patients with cancers, there were trends of increasing median concentrations from stage I and II (surgery), through stage III and locally advanced (surgery), to patients with metastatic disease for both biomarkers (both Ptrends=0.0001) (Figure 2).

Figure 2

(A) Serum tCK18 concentrations by disease groups of increasing stage. (B) Serum cCK18 concentrations by disease groups of increasing stage. The P-values at the top of each figure represent comparisons of medians for all patients with cancer vs controls. The tests for trends used the Cuzick test for nonparametric distributions of ordered groups. Pairwise comparisons (Kruskal–Wallis tests) were also performed as follows: for tCK18, group IIa vs I, P=0.0001; IIb vs IIa, P=0.0001; III vs IIb, P=0.486; and for cCK18, IIa vs I, P=0.615; IIb vs IIa, P=0.0001; III vs IIb, P=0.084. Note: log scale on y axis. Abbreviation: IQR=inter-quartile range.

Among 16 patients undergoing resection, median levels from tumour-draining mesenteric veins were higher than those in peripheral blood for tCK18 (medians: 446 vs 308 U l−1, Ppaired=0.003) and to a lesser extent cCK18 (189 vs 98 U l−1, Ppaired=0.03) (see Supplementary Figure 1A). Median levels of tCK18 reduced after surgery in 11 patients undergoing resection (medians: 310 vs 268 U l−1, P=0.07), but not cCK18 (P=0.37) (see Supplementary Figure 1B).

Prognostic modelling

In 66 patients who underwent curative surgery, increasing tertiles of tCK18 and cCK18 were associated with significantly reduced cancer-specific survival (log rank: P=0.0284 and P=0.0097, respectively) (Figure 3A). For the 55 patients receiving palliative first-line chemotherapy, increasing biomarker levels were associated with reduced overall survival for tCK18 (P=0.0017), but less so for cCK18 (P=0.0690) (Figure 3B).

Figure 3

(A) Cancer-specific survival in patients with CRC undergoing surgical resection with curative intent. (B) Overall survival in patients undergoing first-line conventional chemotherapy for metastatic CRC.

We tested these associations further using multivariate Cox proportional hazards models. In surgically treated patients, the models with tCK18 and cCK18 as continuous variables (log transformed to base 2) were borderline significant, but in general, the prognostic associations of tCK18 and cCK18 were attenuated after adjustment included stage (Table 2). However, in patients with metastatic CRC, increasing concentrations of the biomarkers were independently significant for tCK18 (HR per doubling in concentration=1.78, 95% CI: 1.37, 2.30, P<0.001) and cCK18 (HR per doubling in concentration=1.93, 95% CI: 1.32, 2.81, P=0.001) (Table 3).

Table 2

Multivariate Cox models of total CK18 (tCK18) and caspase cleaved CK18 (cCK18) prognostic influence on cancer-specific survival in patients with colorectal cancer undergoing surgical resection (n=66)

	Categories (U l−1)	Hazard ratio (95% CIs)	P continuous	P vs referent
tCK18 model a
Continuousb	Per unit log2	1.77 (1.04–3.01)	0.036
Categorical
Tertile 1 (n=22)	154–242	1.00
Tertile 2 (n=22)	243–430	1.47 (0.34–6.31)		0.601
Tertile 3 (n=22)	431–3608	2.82 (0.72–11.07)		0.136
cCK18 model a
Continuousb	Per unit log2	1.84 (0.99–3.39)	0.053
Categorical
Tertile 1 (n=22)	45–95	1.00
Tertile 2 (n=22)	96–148	3.41 (0.69–17.01)		0.134
Tertile 3 (n=22)	149–920	4.72 (0.94–23.75)		0.060

Abbreviation: CI=confidence interval.

The calibration and performance characteristics for cancer-specific survival model for M65 were: Harrell’s C-statistic=0.5557 and likelihood ratio for goodness-of-fit=11.924; for M30, C=0.7684 and likelihood ratio for goodness-of-fit=13.946.

All final models included were adjustment for age, sex, stage, and tCK18 and cCK18, respectively.

Interpreted as change in hazard ratio associated with a doubling in CK18 level on the original scale.

Table 3

Multivariate Cox models of total CK18 (tCK18) and caspase cleaved CK18 (cCK18) prognostic influence on overall survival in patients undergoing first-line conventional chemotherapy for metastatic colorectal cancer (n=55)

	Categories (U l−1)	Hazard ratio (95% CIs)	P continuous	P vs referent
tCK18 model a
Continuousb	Per unit log2	1.78 (1.37–2.30)	<0.001
Categorical
Tertile 1 (n=18)	215–370	1.00
Tertile 2 (n=19)	371–1189	3.48 (1.25–9.70)		0.017
Tertile 3 (n=18)	1190–7098	5.42 (2.04–14.3)		0.001
cCK18 model a
Continuousb	Per unit log2	1.93 (1.32–2.81)	0.001
Categorical
Tertile 1 (n=17)	64–179	1.00
Tertile 2 (n=20)	180–409	1.50 (0.62–3.63)		0.367
Tertile 3 (n=18)	410–2388	3.18 (1.27–7.94)		0.013

Abbreviation: CI=confidence interval.

The calibration and performance characteristics for overall survival model for tCK18 were: Harrell’s C-statistic=0.6949 and likelihood ratio for goodness-of-fit=0.867; for cCK18, C=0.6384 and likelihood ratio for goodness-of-fit=3.489.

All final models included were adjustment for age, sex, stage, and tCK18 and cCK18, respectively.

Interpreted as change in hazard ratio associated with a doubling in CK18 level on the original scale.

Longitudinal chemotherapy-related effects

Fifty-six patients receiving chemotherapy underwent repeated sampling (37 patients first-line and 19 patients second- or third-line chemotherapy). Patients with disease progression on chemotherapy had higher median baseline tCK18 (999 vs 389 U l−1, P=0.004) and cCK18 (383 vs 225 U l−1, P=0.019) levels than patients experiencing clinical benefit. In 24 patients with progressive disease, there was an initial decrease in biomarker levels during cycle 1 followed by a progressive upward trend thereafter (Figure 4). In 32 patients with clinical benefit, an initial decrease in biomarker was observed in tCK18 only with no subsequent change.

Figure 4

PD monitoring of tCK18 and cCK18 in 56 patients with metastatic CRC undergoing conventional chemotherapy, categorised by treatment response. The y axis is percentage change in biomarker. Curves generated using cubic splines knotted at days 3, 14, and 21. Shaded area represents 95% CI.

We noted examples where CK18 was more informative than CEA. Thus, of 24 patients with progressive disease on therapy, 9 patients had clinically relevant increases in CEA prior to tCK18 (median 21 days; range 7–28), 3 patients had simultaneous increases in CEA and tCK18, 5 patients had increases in tCK18 prior to CEA (median 21 days; range 7–48), and 7 patients had increases in tCK18 (at median 49 days post start of therapy; range 21–147), but never had clinically relevant CEA rises. Additionally of the 74 patients with metastatic disease who received chemotherapy, 21 subjects had undetectable CEA or levels at the lower limit of quantification.

Finally, we explored whether changes in tCK18 and cCK18 reflect treatment-related toxicity. Using aspirate transaminase levels as a surrogate of hepatic dysfunction, we found no patients with protracted elevated serum concentrations (see Supplementary Figure 2).

Discussion

This study assessed the prognostic and PD utility of circulating CK18 in a range of clinical settings in CRC. The study demonstrated that serum tCK18 and cCK18 levels were elevated in patients with cancer compared with controls, and among patients with cancers, there were trends of increasing median concentrations from early to advanced stages. Levels were higher in tumour-draining venous samples compared with peripheral sampling; and levels fall following tumour resection. Collectively, these observations support the hypothesis that these biomarkers are tumour-derived.

The utility of tCK18 and cCK18 as a diagnostic biomarker of malignancy is limited by the overlap with healthy controls. However, elevated levels of tCK18 and cCK18 prior to treatment were associated with poorer prognosis in both the surgical and metastatic setting, though in the surgical setting much of this prognostic impact was accounted for by associations with stage.

We explored the utility of serum tCK18 and cCK18 as a PD biomarker in patients receiving chemotherapy for metastatic CRC and found that patients with progressive disease on therapy have elevated serum CK18 pre-treatment, tend to show a reduction in these biomarkers at 14–21 days, but then levels progressively increase, returning to and exceeding pre-treatment levels. This is in contrast to patients with clinical benefit from therapy where pre-treatment levels are lower, and after commencement of therapy, either decrease and then flatten (tCK18) or remain constant throughout (cCK18). Collectively, these findings suggest the potential utility of circulating CK18 as a surrogate response biomarker in CRC (Cancer Research UK Glossary, 2011).

Study strengths and limitations

This study has the advantage that we evaluated previously validated biomarkers where profiles had already been examined in relevant xenograft models (Cummings ). We performed prognostic modelling of tCK18- and cCK18-derived data in both surgically treated patients and patients with metastatic disease, and showed different patterns – perhaps, not surprisingly, in patients with metastatic disease where biomarkers levels were most elevated. We then prospectively evaluated changes in serum tCK18 and cCK18 throughout chemotherapy and beyond allowing us to explore longitudinal patterns. This approach will be required to inform future biomarker-adaptive trial designs (Wu ). Our study demonstrated the feasibility of undertaking biomarker dynamic studies in ‘real-life’ clinical practice.

The study findings should be interpreted in light of some limitations. First, patient recruitment was in ‘standard of care’ settings such that treatment selection and follow-up protocols were not precisely controlled. We addressed these weaknesses with multiple adjustments in our multivariate analyses. Although sample sizes per setting were relatively small, by evaluating a number of clinical scenarios, we were able to test a singular hypothesis; that circulating CK18 is tumour-derived.

Context of other studies

Elevated levels of circulating tCK18 and cCK18 have been reported in CRC (Koelink ; Ausch , 2009b; Brandt ). Similar to the present study, Koelink showed in 49 patients with surgically treated CRC that both biomarkers correlated with stage and tumour volume, but reported significantly independent increased risk of poor outcome in their multivariate analysis for cCK18 but not for tCK18. Ausch suggested that cCK18 may be a more accurate reflection of tumour burden than tCK18 in patients with surgically resected CRC. Reductions in CK18 levels measured by both tCK18 and cCK18 post surgical resection were noted (Koelink ; Ausch ).

The fall in tCK18 seen following chemotherapy in CRC agrees with data previously published in testicular cancer (de Haas ) and SCLC (Hou ). However, in the former, there were no patients who progressed during chemotherapy; in the latter, samples were taken only throughout the first cycle from patients responding to chemotherapy. Both papers commented that they were unable to distinguish the profile of response from disease progression, as was done in the present study.

Other literature investigating these assays in CRC has not reported on the kinetics of throughout chemotherapy. In the 36 patients with CRC evaluated by Scott , they only took samples on day 2. In a mixed study of 35 patients with gastrointestinal cancers (25 with CRC), it was suggested that rising cCK18 levels at the end of the first cycle were associated with clinical benefit (Brandt ), and contrasts with data presented here.

Plausible mechanisms and implications

Initial evaluations of CK18 assays lay in the early assessment of treatment-related tumour death and subsequent prediction of response to therapy (Dean ; Gandhi ). However, in patients treated with chemotherapy, early CK18 changes may be unhelpful, due to overlap with host toxicity (Greystoke ). Although we did not observe any significant changes within the first 48 h, we have shown that for many patients there is a decrease in tCK18, and to a lesser extent cCK18, from 1 week after chemotherapy. We speculate that there is an initial effect of chemotherapy on the population of cells that are chemo-sensitive leading to an initial reduction in overall tumour burden, with the later increases in circulating CK18 reflecting subsequent growth in the population of chemo-resistant cells. If this hypothesis is correct, then these biomarkers could be used as surrogates of treatment response to monitor the development of chemo-resistance and lead to early changes in therapy.

This study highlights the potential utility of these biomarkers to guide care in patients with CRC. In particular, tCK18 might guide drug development in prioritising agents and combinations. On the other hand, cCK18 might have limited utility in monitoring chemotherapy because of the overlap between host toxicity and tumour cell death, but might be most useful in the development of drugs directly targeting the apoptotic pathways (Dean ).

Unanswered questions and future research

Carcinoembryonic antigen is commonly used during chemotherapy and subsequent follow-up in patients with CRC. However, its utility in making early decisions on the effectiveness of therapy is limited by two factors: (i) in approximately a third of patients circulating CEA is undetectable; (ii) in patients with detectable CEA, rising levels can be seen following the first cycle of therapy despite subsequent response (Sorbye and Dahl, 2004). It is likely that the information obtained by measurement of CK18 will be used in conjunction with CEA for clinical decision making, but the optimum method for interpreting these in combination is yet to be determined.

In the future, many MBT agents will be administered in combination with conventional chemotherapy in patients with CRC (Cunningham ). Therefore, knowledge of the behaviour of cell death biomarkers following chemotherapy alone will be vital in the interpretation of biomarker profiles in these trials. This study supports the continued deployment of circulating CK18 into biomarker-enhanced trials, particularly in patients undergoing chemotherapy for metastatic CRC.