Circulating Tumor Cells Count and Morphological Features in Breast, Colorectal and Prostate Cancer

Abstract

Background

Presence of circulating tumor cells (CTC) in patients with metastatic breast, colorectal and prostate cancer is indicative for poor prognosis. An automated CTC (aCTC) algorithm developed previously to eliminate the variability in manual counting of CTC (mCTC) was used to extract morphological features. Here we validated the aCTC algorithm on CTC images from prostate, breast and colorectal cancer patients and investigated the role of quantitative morphological parameters.

Methodology

Stored images of samples from patients with prostate, breast and colorectal cancer, healthy controls, benign breast and colorectal tumors were obtained using the CellSearch system. Images were analyzed for the presence of aCTC and their morphological parameters measured and correlated with survival.

Results

Overall survival hazard ratio was not significantly different for aCTC and mCTC. The number of CTC correlated strongest with survival, whereas CTC size, roundness and apoptosis features reached significance in univariate analysis, but not in multivariate analysis. One aCTC/7.5 ml of blood was found in 7 of 204 healthy controls and 9 of 694 benign tumors. In one patient with benign tumor 2 and another 9 aCTC were detected.

Significance of the study

CTC can be identified and morphological features extracted by an algorithm on images stored by the CellSearch system and strongly correlate with clinical outcome in metastatic breast, colorectal and prostate cancer.

About the Authors:

Sjoerd T. Ligthart

Affiliation: Department of Medical Cell BioPhysics, MIRA Institute, University of Twente, Enschede, The Netherlands

Frank A. W. Coumans

Affiliation: Department of Medical Cell BioPhysics, MIRA Institute, University of Twente, Enschede, The Netherlands

Francois-Clement Bidard

Affiliation: Department of Medical Oncology, Institut Curie, Paris, France

Lieke H. J. Simkens

Affiliation: Department of Medical Oncology, Academic Medical Center, University of Amsterdam, Amsterdam, The Netherlands

Cornelis J. A. Punt

Affiliation: Department of Medical Oncology, Academic Medical Center, University of Amsterdam, Amsterdam, The Netherlands

Marco R. de Groot

Affiliation: Department of Internal Medicine, Medisch Spectrum Twente, Enschede, The Netherlands

Gerhardt Attard

Affiliation: The Royal Marsden NHS Foundation Trust, Sutton, Surrey, United Kingdom

Johann S. de Bono

Affiliation: The Royal Marsden NHS Foundation Trust, Sutton, Surrey, United Kingdom

Jean-Yves Pierga

Affiliation: Department of Medical Oncology, Institut Curie, Paris, France

Leon W. M. M. Terstappen

* E-mail: [email protected]

Affiliation: Department of Medical Cell BioPhysics, MIRA Institute, University of Twente, Enschede, The Netherlands

Introduction

Case studies showing the presence of tumor cells in blood of cancer patients have been reported for more than a century [1]–[4]. These circulating tumor cells (CTC) are extremely rare and technology to reliably enumerate CTC has only become available in recent years [5]. Studies conducted with the validated CellSearch system established the relation between the presence of CTC and poor outcome [6]–[10]. The ability to assess the presence of treatment targets on CTC demonstrates the potential for a real time liquid biopsy and has certainly spurred the interest in CTC [11]–[17].

The recent emergence of different technologies to identify CTC has been accompanied with a large range of CTC definitions urging the need for standardization [5], [18]–[31]. Although different CTC characteristics may indeed relate to clinical outcome [32], it is of utmost importance that counting and characterization of CTC can be done accurately, reproducibly, and is validated in controlled clinical studies. Studies have shown that even in the FDA cleared CellSearch system intra-reviewer, inter-reviewer, and inter-laboratory variability is substantial [33], [34]. Recently, we developed a CTC detection algorithm that counts CTC in images recorded by the CellSearch system [35]. This algorithm was not developed to copy human reviewers that assign events as CTC, but it used survival data of metastatic prostate cancer patients to arrive at a definition that optimally stratified the patients into groups with favorable and unfavorable survival [35]. This algorithm eliminates reviewer variability, is fast and decreases the cost of the CTC assay by elimination of the time a reviewer spends on reviewing the images. Furthermore, it delivers quantitative information about the objects it counts as CTC. CTC are morphologically very heterogeneous 5 and in case studies in breast and colorectal cancer it was found that breast cancer CTC are somewhat rounder than cells from colorectal samples [36] In this study, we validated the CTC algorithm on images of patients with metastatic colorectal and breast cancer and evaluated the correlation of quantitative morphological parameters with clinical outcome of patients with metastatic prostate, colorectal and breast cancer before and after initiation of a new line of therapy.

Materials and Methods

Patients

Blood samples enriched for CTC of patients before they received treatment (baseline) and at first follow-up after initiation of a new line of therapy from prospective multicenter CTC studies were used in this study. For metastatic breast cancer, fluorescence images of 283 blood samples processed on the CellSearch system from 179 patients could be imported from the IMMC-01 study and 482 image sets from 248 patients from the IC2006-04 study [7], [37]. For metastatic colorectal cancer, CellSearch images of 1691 image sets from 507 patients could be imported from the CAIRO-2 and IMMC-06 studies [6], [38]. For castration resistant prostate cancer (CRPC), CellSearch images of 370 image sets from 185 patients could be imported from the IMMC-38 study and 189 image sets from 100 patients could be imported from the Abiraterone study [8], [39], [40]. From the IMMC-26 study 317 image sets from 93 patients with benign colorectal disease and 200 image sets from 61 patients with benign breast disease were imported [41], [42]. In addition, 204 image sets from healthy controls from the IMMC-01 and IMMC-06 study were included in this study. In Table S1 the detailed characteristics of the patients enrolled in these studies are shown. For these studies all patients and healthy volunteers provided written informed consent. The studies were approved by the ethics commissions and institutional review boards of Institut Curie Paris, Royal Marsden, London, Radboud Nijmegen, Medisch Spectrum Twente and centers participating in the original studies [6], [7], [8], [37], [38], [39], [40], [41], [42].

Manual Enumeration of Circulating Tumor Cells

The CellSearch system (Veridex LLC, Huntingdon Valley, PA, USA) was used to enumerate CTC. The system consists of a CellTracks Autoprep® for sample preparation and a CellTracks Analyzer II® for sample analysis [5], [43] The CellTracks Autoprep immuno-magnetically enriches epithelial cells from 7.5 ml of blood targeting the Epithelial Cell Adhesion Molecule (EpCAM). The enriched sample is labeled with phycoerythrin-conjugated (PE) antibodies directed against cytokeratins (CK) 8, 18 and 19, an allophycocyanin-conjugated (APC) antibody against CD45 and the nuclear dye 4′,6-diamidino-2-phenylindole (DAPI). After preparation the enriched sample is transferred to a cartridge contained in a CellTracks MagNest. All ferrofluid labeled objects are pulled towards the imaging surface of the cartridge, which is placed on the CellTracks Analyzer II for image acquisition. The analyzer is a four-color semi-automated fluorescence microscope that captures digital images in four different fluorescent channels using a 10×/0.45 NA objective and a charge-coupled device camera with 6.7×6.7 µm2 sized pixels. For each cartridge, 144–180 4-layer tiff images of DAPI, FITC, PE, APC are saved per patient. After imaging, the system preselects objects that are positive for DNA and CK and shows the images of all fluorescence channels together with a DNA/CK overlay of all selected objects in a gallery. Finally, a trained operator selects objects from the gallery as manual CTC (mCTC) if they are DNA+CK+, CD45-, have a cell-like morphology, and are larger than 4×4 µm2.

Automated Enumeration of Circulating Tumor Cells

The recorded image sets from the patients were copied to central hard drives and analyzed by an algorithm developed in Matlab 2009a, using the DipImage toolbox (Delft University, Delft, The Netherlands). The algorithm was optimized using survival of 185 castration resistant prostate cancer patients as a training parameter and is described in more detail elsewhere [35]. In brief the optimal CTC definition was selected based on 1. High Cox Hazard ratio for borh baseline and follow-up samples; 2 Higher HR for folllow-up than baseline samples; 3 Low relative and absolute number of CTC in 68 control samples. The algorithm was then applied on 144–180 images from each patient. First objects in the CK-PE channel are selected by dynamic thresholding and their outline is saved [44]. Next, the outline of each object is used to perform measurements in the DNA-DAPI, CK-PE, and CD45-APC channels. Finally, the algorithm, optimized using survival data from metastatic prostate cancer patients, selects automated CTC (aCTC). The inclusion criteria for an aCTC were defined as a CK-PE standard deviation higher than 50 counts, a DNA-DAPI peak value of at least 170 counts, a CD45-APC peak value less than 60 counts and a size within 75–500 pixels (35–225 µm2). However, the upper aCTC size limit was increased to 2000 pixels (∼900 µm2) after the observation that aCTC sizes in metastatic breast cancer are larger [45]. For every patient image set, the included objects were added to arrive at a final aCTC count per patient. To compare the overall aCTC and mCTC enumeration, a frequency plot was created using the data from the three patient groups, benign tumors and healthy controls.

Additional Morphological Measurements on aCTC

We performed a number of measurements on the aCTC, including cell roundness, nuclear to cytokeratin (NCK) ratio, the presence of cell clusters, and the number of speckles within the cytoplasm. Cell roundness was defined as the ratio of the signal perimeter squared divided by 4π times the signal area (termed P2A in the literature). Thus, roundness equals one for a perfect circle. The NCK-ratio was defined as the area of the DNA-DAPI signal, divided by the area of the CK-PE signal. For this parameter, an extra segmentation step for DNA-DAPI was performed to find DNA-DAPI+ objects larger than 10 µm2 (22 pixels, the smallest DNA-DAPI nucleus area). If multiple DNA-DAPI+ objects were found close to the CK signal, the center-of-mass distances between the CK and DNA signals were determined. The object with (partial) overlap between CK and DNA and the smallest distance (no more than 17 µm, or half the size limit) was selected to determine the NCK-ratio. Clusters of cells were defined as the number of DNA-DAPI objects larger than 10 µm2 and located within each CK-PE signal. Finally, speckles within the CK-PE signals are indicative of damaged or apoptotic CTC [46]. The number of dot-like structures within the CK-PE signal was counted and used as a measure of apoptosis. Dot-like structures were counted by an algorithm developed for the detection of fluorescence in situ hybridization probes within a cell nucleus [47].

Identification and Morphological Measurements of Leukocytes

Due to a-specific binding of EpCAM-ferrofluid particles during the CTC enrichment, leukocytes are carried over and are present at the analysis surface of the sample cartridge. Leukocytes are identified by their DAPI staining, expression of CD45 and lack of cytokeratin expression. We applied the triangle threshold method to determine the outline of leukocytes in the CD45-APC channel [44]. Leukocytes were defined as objects that had CK-PE peak value lower than 100 counts, DNA-DAPI peak value of at least 100 counts, CD45-APC peak value of at least 100 counts and a size range between 50 and 2000 pixels (21–900 µm2). Morphological parameters were extracted as described for aCTC.

Statistical Analysis

SPSS 16 (IBM, New York, USA) was used for survival analysis. Least-squares linear regression analysis was performed to compare aCTC and mCTC, baseline and follow-up samples, and to compare different studies of the same type of cancer. Kaplan-Meier analysis was performed by dichotomizing the patient population on the clinically used cut-off value of five CTC for breast and prostate cancer, and three CTC for colorectal cancer. Beyond these cut-offs, other cut-off values of 1, 5, 10, and 100 for aCTC or mCTC were tested to dichotomize the patient population per study for univariate Cox hazard regression. Log-rank tests were used to compare survival between groups. Cox hazard ratios (HRs) with 95% confidence intervals (CI) were determined and collected in a table. Survival was defined as the time between baseline sample and death from any cause.

To determine whether the morphological parameters determined on aCTC provide additional information that is relevant for prognosis, multivariate Cox proportional hazards regression was performed (forward conditional, with pin 0.10 and pout 0.20), controlling for cancer type (colorectal, prostate, breast), with all morphological parameters and the aCTC count as continuous variables. The number of aCTC was 10log transformed before inclusion in the regression [48], [49]. Because the morphological parameters could not be determined on patients without aCTC, only patients with at least 1 aCTC were included. For parameters that are determined for each individual aCTC, the median value within a patient was used since this represents the predominant characteristic. The nonparametric Mann-Whitney U-test was used to compare morphological parameters of different types of cancer. The t-test could not be used because the parameters are not normally distributed.

Results

Frequency of aCTC Versus mCTC

aCTC and mCTC were enumerated in healthy controls, patients with benign disease and patients with metastatic colorectal, prostate and breast cancer before initiation of therapy. The frequency distribution of the CTC is shown in figure 1 and the 25th, 50th, and 75th percentiles are provided at the right hand side of the figure. Figure 2 shows the scatter plot of aCTC versus mCTC per study, as well as the linear regression coefficients. Detailed quantitative information about counted ranges, median, mean, SD, and linear regression coefficients are presented in Table S2. In summary, breast CTC counts ranged from 0–10111 aCTC (median 2/0 for baseline/follow-up) and 0–23618 mCTC (median 3/0 for baseline/follow-up); colorectal CTC counts from 0–448 (median 1/0) and 0–351 mCTC (median 0/0); prostate CTC counts ranged from 0–4970 aCTC (median 7/3) and 0–5925 mCTC (median 7/2) per 7.5 ml of blood. Benign patients ranged from 0–218 (median 0) aCTC and 0–12 mCTC (median 0). Finally, healthy controls had up to 1 aCTC and mCTC (both medians 0). Correlation R2 between mCTC and aCTC was 0.48 for 725 breast samples, 0.84 for 1729 colorectal samples, and 0.68 for 559 prostate samples.

[Figure omitted. See PDF.]

Figure 1. Scatter plot of mCTC and aCTC by cancer type.

The frequencies are shown of mCTC and aCTC found in healthy controls, patients with benign disease and patients with metastatic breast, prostate, and colorectal cancer before initiation of a new line of therapy. The dots in the scatter plot are spread for viewing purposes. The 25th, 50th, and 75th percentile are given by the black markers and in numerical values in the right part of the figure. At the start of each row, the percentage and number of patients that have 0 CTC are shown.

https://doi.org/10.1371/journal.pone.0067148.g001

[Figure omitted. See PDF.]

Figure 2. Scatter plot of the number of aCTC versus mCTC for every study.

Horizontal and vertical dashed lines were drawn to show clinically used cut-offs, which are ≥5 for breast and prostate, and ≥3 for colorectal cancer. At the top of each graph, the linear regression coefficients are shown.

https://doi.org/10.1371/journal.pone.0067148.g002

Measurements of aCTC Morphological Parameters

Recording of the outline of aCTC permits the extraction and quantification of morphological parameters from the identified aCTC. The distribution of the size, NCK ratio and roundness of aCTC in blood of patients with metastatic breast, colorectal, and prostate cancer before initiation of therapy and at first follow up after initiation of therapy is shown in figure 3. Examples for these three morphological parameters are shown at the top of the figure, with histograms for CTC from colorectal, breast and prostate cancer below each parameter. The aCTC size (µm2) distribution for colorectal, breast and prostate cancer is shown in panels A, B and C respectively. The NCK-ratio (range 0.5 to 5) for breast, colorectal, and prostate cancer CTC is shown in panels E, F and G respectively. aCTC appear in a large range of shapes. To illustrate this diversity in CTC shape, a measure of roundness (range 1 to 5) for colorectal, breast and prostate cancer CTC is shown in panels H, I, and J respectively. For comparison the size and roundness properties of leukocytes are shown in figure 3, panels D and K. Leukocytes have a narrow size distribution, as can be seen in panel D. The NCK-ratio is not shown as leukocytes do not express cytokeratin. R2 correlation between measurements before and after initiation of therapy is provided in each panel. aCTC morphological parameters were quantified by mean and standard deviation and are provided in table 1. To test whether the morphological parameters from table 1 were significantly different between cancer types, Mann-Whitney U-test p-values were derived using all image sets of each type of cancer. Derived p-values for median size of aCTC between colorectal and prostate was 0.90 and 0.35 for median NCK-ratio between breast and colorectal. The percentage apoptotic cells between breast and prostate had p = 0.18. For median roundness of aCTC p-values were 0.27, 0.008, and 0.15 for comparison between breast-colorectal, breast-prostate, and colorectal-prostate, respectively. For all other morphological parameters from table 1 p was less than 0.001 between all types of cancer.

[Figure omitted. See PDF.]

Figure 3. Histograms of morphological parameters of aCTC by cancer type.

Breast (panel A, E, H), colorectal (panel B, F, I), and prostate samples (panel C, G, J) are shown next to morphological parameters of leukocytes (panel D, K). Baseline samples (black bars) are shown next to first follow-up samples (white bars). Coefficients of determination of a linear regression between the time points are shown next to the bars. Examples of morphological parameter values of aCTC are given at the top of the figure, as well as examples of the outlines after segmentation. The objects seem smaller in the images than in the outlines. This is due to scaling of the image for print, the outlines show the area of the object that is distinguishable from noise, and thus represents the true size of the object. The 5 µm scale bar applies to all images if not stated differently. N/A: not applicable.

https://doi.org/10.1371/journal.pone.0067148.g003

[Figure omitted. See PDF.]

Table 1. Size, roundness and nuclear to cytokeratin ratio of CTC in breast, colorectal cancer and leukocytes present in the CTC enriched samples.

https://doi.org/10.1371/journal.pone.0067148.t001

The size of the objects can be influenced by the presence of multiple CTC representing CTC clusters. These clusters are counted as one CTC in both the manual identification as well as the automated identification of CTC. The presence of aCTC clusters was estimated by the number of nuclei that could be detected within the CK-PE area. Percentages of clustered cells per patient are shown in table 2 for patients with ≥5 aCTC for image sets of both baseline and follow-up. None of the patients showed a significant rise or drop (outside 95% CI) in the number of clusters after initiation of therapy. The presence of speckles of cytokeratin is thought to be associated with CTC undergoing apoptosis. As a measure of these speckled aCTC the percentage of all aCTC that have 3 or more dot-like structures in their cytokeratin signal is shown in table 2 for both baseline and follow-up samples. None of the patients had a significant rise or drop in the number of speckled aCTC between baseline and first follow-up measurement.

[Figure omitted. See PDF.]

Table 2. Clusters and speckled aCTC at baseline and follow-up samples from breast, colorectal and prostate cancer patients with ≥5 aCTC.

https://doi.org/10.1371/journal.pone.0067148.t002

mCTC, aCTC and Morphological Parameters Versus Survival in Metastatic Breast, Colorectal and Prostate Cancer

To illustrate the relationship between survival and CTC, Kaplan Meier plots were generated dividing the patients in favorable and unfavorable groups based on a threshold of 3 CTC for metastatic colorectal cancer and 5 CTC for metastatic prostate and breast cancer. The Kaplan-Meier plots shown in figure 4 are for aCTC and mCTC before initiation of therapy (panels A, C, E) and at first follow-up after initiation of therapy (panels B, C, F). The relationship between aCTC number and overall survival in metastatic breast colorectal, and prostate cancer is shown in table 3. To investigate the impact of aCTC morphological parameters on survival, univariate and multivariate analysis was performed using the measured morphological parameters as covariates next to the number of aCTC or mCTC. HRs were calculated for each covariate as a continuous variable. Table 4 shows the HRs, derived p-values, and chi-squared.

[Figure omitted. See PDF.]

Figure 4. Kaplan-Meier plots of baseline and follow-up samples by cancer type, and mCTC, or aCTC.

Metastatic breast (panel A and B), colorectal (panel C and D), and prostate cancer patients (panel E and F) are shown with favorable and unfavorable mCTC (black lines) and aCTC definition (grey lines). For every type of cancer, patients were pooled from the two studies of each type as shown in Table S1.

https://doi.org/10.1371/journal.pone.0067148.g004

[Figure omitted. See PDF.]

Table 3. Relation between aCTC number and overall survival in metastatic colorectal, breast and prostate cancer.

https://doi.org/10.1371/journal.pone.0067148.t003

[Figure omitted. See PDF.]

Table 4. Univariate analysis of variables for all follow-up samples used in this study.

https://doi.org/10.1371/journal.pone.0067148.t004

Next, all variables of aCTC were entered in a multivariate Cox regression analysis. mCTC were left out of this analysis, as their quantitative parameters are unknown. Hazard models were calculated using all first follow-up image sets for all cancer types together and for each type separately. When all image sets were combined, we controlled for differences in overall survival due to cancer type. Results are shown in table 5. The optimal model used the number of aCTC, and their median roundness per patient and had a Chi-square of 181 with 3 degrees of freedom (with the third degree the control variable for colorectal cancer). The model shows that more objects impose a bigger threat and that objects with a roundness close to one may be more dangerous to patients than objects that are elongated. When including only samples from one type of cancer, the median size of the aCTC, and the median NCK-ratio of the aCTC had a significant influence on survival in breast cancer samples. HRs smaller than one show that in breast cancer small objects represent greater hazard than large objects. Furthermore, objects with small DAPI signals with respect to their CK signal, i.e. a low NCK-ratio, are more dangerous for patients than cells with high NCK-ratio. Other morphological parameters did not improve the multivariate models.

[Figure omitted. See PDF.]

Table 5. Multivariate analyses of aCTC variables on all first follow-up samples, and for each type of cancer.

https://doi.org/10.1371/journal.pone.0067148.t005

Discussion

Recently, we showed that counting CTC by an automated algorithm is preferable to manual counting by a trained reviewer [35]. The algorithm performed similar to the human reviewer in terms of prognostic value: similar HRs were found. In addition, the algorithm takes no operator time and has 0% variability, against an inter-laboratory variability of 4% to 31% for mCTC (median 14%) reported previously by another group [33]. The aCTC algorithm was optimized to obtain the best separation in survival risk and not to obtain the best correlation with mCTC. The training set used image sets from a prospective metastatic prostate cancer study and was validated on an independent data set of prostate cancer patients included in the Abiraterone phase I/II trials [39], [40]. The aCTC algorithm with the highest prognostic value included objects with DNA and CK, without CD45, and which have a size between 34 µm2 and 224 µm2.

Here we applied the aCTC algorithm to image sets of patients with benign and malignant breast and colorectal cancer, as well as to image sets of healthy controls. The upper size limit was increased from 224 µm2 to 900 µm2, concordant with a recent study showing a larger size of CTC in breast cancer [45]. This change resulted in an aCTC increase by a factor of 1.3 for prostate cancer with less than 15% effect on the HR. The increase in the number of aCTC was a factor 1.4 for colorectal image sets and 1.9 for breast cancer image sets, and thus indicative for the larger average size of CTC in breast cancer. The frequency distribution of CTC illustrated in figure 1 shows that the aCTC and mCTC counts do not differ greatly. Separation of the results into the different studies however showed a clear discrepancy between the counts obtained with samples from the IMMC-01 breast cancer study as compared to all other studies (figure 2b). The cause of this discrepancy was traced back to the CD45 staining of objects, as the CD45 staining in the IMMC-01 study was much higher as compared to the IC2006–04 study. This resulted in the rejection of many objects in the IMMC-01 study. Lowering of the CD45 threshold eliminated the discrepancy between the number of aCTC in the IMMC-01 and other studies (data not shown). IMMC-01 was the first of all studies in this report, potential causes for the higher CD45-APC signal in IMMC-01 are a higher concentration of the CD45-APC fluorochrome, differences in sample preparation or differences in filter cubes giving rise to more crosstalk of CK-PE signal into the CD45-APC channel. Most image sets of IMMC-01, 06 and 38 were scanned on the CellSpotter Analyzer instead of the CellTracks Analyzer II used for the other studies, therefore differences between CellTracks and CellSpotter do not offer an explanation for the differences in CD45-APC signal. The IMMC-01 study was the first large clinical trial run with the CellSearch system and during this study blood samples were not automatically processed by the CellSearch Autoprep system. Kaplan-Meier analysis on the IC2006–04 study alone showed that the aCTC and mCTC had the same prognostic value. Therefore, future breast cancer blood samples, processed using current CellSearch equipment, can be accurately analysed with the algorithm. As different microscopes (CellSpotter Analyzer & CellTracks Analyzer II) were used in this study one might assume that other microscopes can be used for the analysis of the images with this algorithm. One however has to keep in mind that the same objective (10X, NA 0,45) was used. For successful application of another imaging system it is of great importance that the size in the sample corresponding to a pixel, and the typical intensity of a tumor cell will need to be matched. The control samples did not contain more than one aCTC or mCTC, showing that these definitions have very few false positives. The aCTC frequencies detected in patients with benign tumors were slightly less as compared to the mCTC frequencies, as were the numbers found in cancer patient blood samples. From Figure 4 and table 3 it is evident that identification of CTC by manual review and the automated algorithm performs equally well in discriminating between patients with favourable and unfavourable outcome. Except for the follow-up image sets of the colorectal cancer patients. The latter may be contributed to the low incidence of CTC for this group of patients.

Comparison of aCTC morphological parameters between cancer types shows a larger size, a relatively large nucleus to cytokeratin ratio and rounder CTC in breast cancer as compared to prostate and colorectal cancer. The distribution of the size of aCTC is exponential in CTC of both prostate and colorectal cancer, which is cut off at the lower end by the CTC algorithm, in contrast with leukocytes that display a more bell-shaped size distribution. These distributions shows that a wide range of CTC and CTC fragments are present in the blood of patients, as was reported earlier [46]. It was found in earlier work that the number of tumor micro particles (TMPs) is also prognostic for survival [32], it is therefore interesting to include these TMPs in the analysis, as they are present in abundant numbers. Especially in colorectal cancer, TMPs might reduce the large Poisson error that is currently present by using a clinical cut-off value of three mCTC. Clusters and speckled cells were found among aCTC, but not among leukocytes. A few exceptional colorectal samples contained a large number of apoptotic cells and clusters, but the presence of these cells or clusters had no detectable impact on survival. Breast cancer aCTC showed the highest number of speckled cells, while prostate cancer aCTC had more clusters, although the differences between cancer types are small. No significant rise or drop in the number of clusters or apoptotic cells was found between baseline and follow-up samples in any patient.

We performed continuous univariate and multivariate survival analysis to investigate the influence of aCTC and its morphological parameters on survival. CTC measurements at first-follow up after treatment were used for this analysis as the CTC “surviving” this treatment more strongly related to survival and might contain typical morphological features. Median roundness was significant (p<0.05) when we included all cancer types in the multivariate analysis, with rounder cells predicting poorer survival. When we performed multivariate analysis by cancer type, small cells with small nuclei compared to the cytoplasm were significantly (p<0.05) associated with poor survival only for breast cancer patients. Currently, we have no biological explanation for this phenomenon. However, research focused on separating normal blood cells from CTC on the basis of size should consider this observation. The association of rounder and/or smaller cells with poorer survival suggest that these maybe features of cancer stem cells, similar to the observations made for hematopoietic stem cells, which are clearly smaller and rounder as compared with the other hematopoietic progenitor cells [50]. If this is indeed the case size dependent methods for CTC enrichment will not be able to detect these putative cancer stem cells.

In conclusion, we have validated a CTC automated algorithm, which was trained on samples from metastatic prostate cancer patients, on other types of cancer samples and samples from patients with benign tumors and healthy controls. We show that differences exist in physical characteristics between CTC from different origin, although their impact on survival is limited. Standardized identification and counting of CTC is imperative to move the field forward and will enable the molecular characterization of those objects/cells that are indeed CTC.

Supporting Information

[Figure omitted. See PDF.]

Table S1.

Patient characteristics.

https://doi.org/10.1371/journal.pone.0067148.s001

(DOCX)

Table S2..

https://doi.org/10.1371/journal.pone.0067148.s002

(DOCX)

Author Contributions

Conceived and designed the experiments: SL FC FCB GA LT. Performed the experiments: SL FC. Analyzed the data: SL FC FCB GA LT. Contributed reagents/materials/analysis tools: SL FC FCB LS CP MdG GA JdB JYP. Wrote the paper: SL FC FCB GA LT.

Citation: Ligthart ST, Coumans FAW, Bidard F-C, Simkens LHJ, Punt CJA, de Groot MR, et al. (2013) Circulating Tumor Cells Count and Morphological Features in Breast, Colorectal and Prostate Cancer. PLoS ONE 8(6): e67148. https://doi.org/10.1371/journal.pone.0067148

References

1. Ashworth TR (1869) A case of cancer in which cells similar to those in the tumours were seen in the blood after death. Australian Medical Journal 14: 146–147.

2. Engell HC (1955) Cancer cells in the circulating blood; a clinical study on the occurrence of cancer cells in the peripheral blood and in venous blood draining the tumour area at operation. Acta chirurgica Scandinavica Supplementum 201: 1–70.

3. Myerowitz RL, Edwards PA, Sartiano GP (1977) Carcinocythemia (Carcinoma Cell Leukemia) due to metastatic carcinoma of breast- report of a case. Cancer 40: 3107–3111.

4. Yam LT, Janckila AJ (1987) Immunocytodiagnosis of carcinocythemia in disseminated breast cancer. Acta Cytologica 31: 68–72.

5. Allard WJ, Matera J, Miller MC, Repollet M, Connelly MC, et al. (2004) Tumor cells circulate in the peripheral blood of all major carcinomas but not in healthy subjects or patients with nonmalignant diseases. Clinical Cancer Research 10: 6897–6904.

6. Cohen SJ, Punt CJA, Iannotti N, Saidman BH, Sabbath KD, et al. (2008) Relationship of circulating tumor cells to tumor response, progression-free survival, and overall survival in patients with metastatic colorectal cancer. Journal of Clinical Oncology 26: 3213–3221.

7. Cristofanilli M, Budd GT, Ellis MJ, Stopeck A, Matera J, et al. (2004) Circulating tumor cells, disease progression, and survival in metastatic breast cancer. New England Journal of Medicine 351: 781–791.

8. de Bono JS, Scher HI, Montgomery RB, Parker C, Miller MC, et al. (2008) Circulating Tumor Cells Predict Survival Benefit from Treatment in Metastatic Castration-Resistant Prostate Cancer. Clinical Cancer Research 14: 6302–6309.

9. Krebs MG, Sloane R, Priest L, Lancashire L, Hou JM, et al. (2011) Evaluation and prognostic significance of circulating tumor cells in patients with non-small-cell lung cancer. Journal of clinical oncology 29: 1556–1563.

10. Matsusaka S, Chin K, Ogura M, Suenaga M, Shinozaki E, et al. (2010) Circulating tumor cells as a surrogate marker for determining response to chemotherapy in patients with advanced gastric cancer. Cancer science 101: 1067–1071.

11. Attard G, Swermenhuis JF, Olmos D, Reid AHM, Vickers E, et al. (2009) Characterization of ERG, AR and PTEN Gene Status in Circulating Tumor Cells from Patients with Castration-Resistant Prostate Cancer. Cancer Research 69: 2912–2918.

12. de Bono JS, Attard G, Adjei A, Pollak MN, Fong PC, et al. (2007) Potential applications for circulating tumor cells expressing the insulin-like growth factor-I receptor. Clinical Cancer Research 13: 3611–3616.

13. Hayes DF, Walker TM, Singh B, Vitetta ES, Uhr JW, et al. (2002) Monitoring expression of HER-2 on circulating epithelial cells in patients with advanced breast cancer. International Journal of Oncology 21: 1111–1117.

14. Meng S, Tripathy D, Shete S, Ashfaq R, Saboorian H, et al. (2006) uPAR and HER-2 gene status in individual breast cancer cells from blood and tissues. Proceedings of the National Academy of Sciences of the United States of America 103: 17361–17365.

15. Meng SD, Tripathy D, Shete S, Ashfaq R, Haley B, et al. (2004) HER-2 gene amplification can be acquired as breast cancer progresses. Proceedings of the National Academy of Sciences of the United States of America 101: 9393–9398.

16. Swennenhuis JF, Tibbe AG, Levink R, Sipkema RC, Terstappen LW (2009) Characterization of circulating tumor cells by fluorescence in situ hybridization. Cytometry Part A 75: 520–527.

17. van de Stolpe A, Pantel K, Sleijfer S, Terstappen LW, den Toonder JMJ (2011) Circulating Tumor Cell Isolation and Diagnostics: Toward Routine Clinical Use. Cancer Research 71: 5955–5960.

18. Racila E, Euhus D, Weiss AJ, Rao C, McConnell J, et al. (1998) Detection and characterization of carcinoma cells in the blood. Proceedings of the National Academy of Sciences of the United States of America 95: 4589–4594.

19. Alix-Panabieres C, Brouillet JP, Fabbro M, Yssel H, Rousset T, et al. (2005) Characterization and enumeration of cells secreting tumor markers in the peripheral blood of breast cancer patients. Journal of immunological methods 299: 177–188.

20. Kahn HJ, Presta A, Yang LY, Blondal J, Trudeau M, et al. (2004) Enumeration of circulating tumor cells in the blood of breast cancer patients after filtration enrichment: correlation with disease stage. Breast Cancer Research and Treatment 86: 237–247.

21. Krivacic RT, Ladanyi A, Curry DN, Hsieh HB, Kuhn P, et al. (2004) A rare-cell detector for cancer. Proceedings of the National Academy of Sciences of the United States of America 101: 10501–10504.

22. Nagrath S, Sequist LV, Maheswaran S, Bell DW, Irimia D, et al. (2007) Isolation of rare circulating tumour cells in cancer patients by microchip technology. Nature 450: 1235–1239.

23. Pachmann K, Camara O, Kavallaris A, Krauspe S, Malarski N, et al. (2008) Monitoring the response of circulating epithelial tumor cells to adjuvant chemotherapy in breast cancer allows detection of patients at risk of early relapse. Journal of Clinical Oncology 26: 1208–1215.

24. Vona G, Sabile A, Louha M, Sitruk V, Romana S, et al. (2000) Isolation by size of epithelial tumor cells - A new method for the immunomorphological and molecular characterization of circulating tumor cells. American Journal of Pathology 156: 57–63.

25. Gleghorn JP, Pratt ED, Denning D, Liu H, Bander NH, et al. (2010) Capture of circulating tumor cells from whole blood of prostate cancer patients using geometrically enhanced differential immunocapture (GEDI) and a prostate-specific antibody. Lab on a Chip 10: 27–29.

26. Hsieh HB, Marrinucci D, Bethel K, Curry DN, Humphrey M, et al. (2006) High speed detection of circulating tumor cells. Biosensors & Bioelectronics 21: 1893–1899.

27. Saliba AE, Saias L, Psychari E, Minc N, Simon D, et al. (2010) Microfluidic sorting and multimodal typing of cancer cells in self-assembled magnetic arrays. Proceedings Of The National Academy Of Sciences Of The United States Of America 107: 14524–14529.

28. Schwarzenbach H, Alix-Panabieres C, Muller I, Letang N, Vendrell JP, et al. (2009) Cell-free Tumor DNA in Blood Plasma As a Marker for Circulating Tumor Cells in Prostate Cancer. Clinical Cancer Research 15: 1032–1038.

29. Stakenborg T, Liu CX, Henry O, O'Sullivan CK, Fermer C, et al.. (2010) Lab-on-a-chip for the isolation and characterization of circulating tumor cells. 2010 Annual International Conference of the Ieee Engineering in Medicine and Biology Society. 292–294.

30. Talasaz AH, Powell AA, Huber DE, Berbee JG, Roh KH, et al. (2009) Isolating highly enriched populations of circulating epithelial cells and other rare cells from blood using a magnetic sweeper device. Proceedings Of The National Academy Of Sciences Of The United States Of America 106: 3970–3975.

31. Zheng SY, Lin HK, Lu B, Williams A, Datar R, et al. (2011) 3D microfilter device for viable circulating tumor cell (CTC) enrichment from blood. Biomedical Microdevices 13: 203–213.

32. Coumans FAW, Doggen CJM, Attard G, de Bono JS, Terstappen LWMM (2010) All circulating EpCAM+CK+CD45- objects predict overall survival in castration-resistant prostate cancer. Annals of Oncology 21: 1851–1857.

33. Kraan J, Sleijfer S, Strijbos MH, Ignatiadis M, Peeters D, et al. (2010) External quality assurance of circulating tumor cell enumeration using the CellSearch((R)) system: A feasibility study. Cytometry B Clin Cytom 80B: 112–118.

34. Tibbe AGJ, Miller MC, Terstappen L (2007) Statistical considerations for enumeration of circulating tumor cells. Cytometry Part A 71A: 154–162.

35. Ligthart ST, Coumans F, Attard G, Mulick Cassidy A, De Bono JS, et al.. (2011) Unbiased and Automated Identification of a Circulating Tumour Cell Definition that Associates with Overall Survival. Plos One 6.

36. Marrinucci D, Bethel K, Bruce RH, Curry DN, Hsieh B, et al. (2007) Case study of the morphologic variation of circulating tumor cells. Human Pathology 38: 514–519.

37. Pierga JY, Hajage D, Bachelot T, Delaloge S, Brain E, et al. (2011) High independent prognostic and predictive value of circulating tumor cells compared with serum tumor markers in a large prospective trial in first-line chemotherapy for metastatic breast cancer patients. Annals of Oncology 23: 618–624.

38. Tol J, Koopman M, Miller MC, Tibbe A, Cats A, et al. (2010) Circulating tumour cells early predict progression-free and overall survival in advanced colorectal cancer patients treated with chemotherapy and targeted agents. Annals of Oncology 21: 1006–1012.

39. Attard G, Reid AHM, A'Hern R, Parker C, Oommen NB, et al. (2009) Selective Inhibition of CYP17 With Abiraterone Acetate Is Highly Active in the Treatment of Castration-Resistant Prostate Cancer. Journal of Clinical Oncology 27: 3742–3748.

40. Attard G, Reid AHM, Yap TA, Raynaud F, Dowsett M, et al. (2008) Phase I clinical trial of a selective inhibitor of CYP17, abiraterone acetate, confirms that castration-resistant prostate cancer commonly remains hormone driven. Journal of Clinical Oncology 26: 4563–4571.

41. Franken B, De Groot MR, Mastboom WJB, Vermes I, Van der Palen J, et al. (2012) Circulating Tumor Cells, Disease Recurrence and Survival in Newly Diagnosed Breast Cancer. Breast Cancer Research 14: R133.

42. Scholten L, Terstappen LWMM, van der Palen J, Mastboom WJB, Tibbe AGJ, et al.. (2012) Circulating tumor cells as a possible prognostic tool in newly diagnosed nonmetastatic colorectal cancer? Journal of Clinical Oncology 30: suppl 4; abstr 395.

43. Kagan M, Howard D, Bendele T, Mayes J, Silvia J, et al. (2002) A sample preparation and analysis system for identification of circulating tumor cells. Journal of Clinical Ligand Assay 25: 104–110.

44. Zack GW, Rogers WE, Latt SA (1977) Automatic-Measurement Of Sister Chromatid Exchange Frequency. Journal Of Histochemistry & Cytochemistry 25: 741–753.

45. Ligthart ST, Bidard FC, Decraene C, Bachelot T, Delaloge S, et al. (2013) Unbiased quantitative assessment of Her-2 expression of circulating tumor cells in patients with metastatic and non metastatic breast cancer. Annals of Oncology 24: 1231–1238.

46. Larson CJ, Moreno JG, Pienta KJ, Gross S, Repollet M, et al. (2004) Apoptosis of circulating tumor cells in prostate cancer patients. Cytometry Part A 62A: 46–53.

47. Ligthart ST, Swennenhuis JF, Greve J, Terstappen LWMM (2012) FISH Probe Counting in Circulating Tumor Cells. In: Hamilton G. Cytokeratins- Tools in Oncology. New York: InTech. 158 p.

48. Scher HI, Jia XY, de Bono JS, Fleisher M, Pienta KJ, et al. (2009) Circulating tumour cells as prognostic markers in progressive, castration-resistant prostate cancer: a reanalysis of IMMC38 trial data. Lancet Oncology 10: 233–239.

49. Coumans F, Ligthart ST, Uhr J, Terstappen LWMM (2012) Challenges in the enumeration and phenotyping of CTC. Clinical Cancer Research 18: 5711–5718.

50. Terstappen LWMM, Huang S, Safford M, Lansdorp PM, Loken MR (1991) Sequential generations of hematopoietic colonies derived from single non lineage committed progenitor cells. Blood 77: 1218–1227.