Next Article in Journal
Different Calculation Strategies Are Congruent in Determining Chemotherapy Resistance of Brain Tumors In Vitro
Next Article in Special Issue
Detecting Resistance to Therapeutic ALK Inhibitors in Tumor Tissue and Liquid Biopsy Markers: An Update to a Clinical Routine Practice
Previous Article in Journal
Pathogenic Genome Signatures That Damage Motor Neurons in Amyotrophic Lateral Sclerosis
Previous Article in Special Issue
Liquid Biopsy in Small Cell Lung Cancer—A Route to Improved Clinical Care?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Endothelium-Derived Extracellular Vesicles Associate with Poor Prognosis in Metastatic Colorectal Cancer

by
Afroditi Nanou
1,*,
Linda Mol
2,
Frank A. W. Coumans
1,
Miriam Koopman
3,
Cornelis J. A. Punt
4 and
Leon W. M. M. Terstappen
1,*
1
Department of Medical Cell BioPhysics, University of Twente, 7522ND Enschede, The Netherlands
2
Netherlands Comprehensive Cancer Organization, 6533AA Nijmegen, The Netherlands
3
Department of Medical Oncology, University Medical Center Utrecht, 3584CS Utrecht, The Netherlands
4
Julius Center for Health Sciences and Primary Care, Department of Epidemiology, University Medical Center Utrecht, 3584CG Utrecht, The Netherlands
*
Authors to whom correspondence should be addressed.
Cells 2020, 9(12), 2688; https://doi.org/10.3390/cells9122688
Submission received: 23 October 2020 / Revised: 1 December 2020 / Accepted: 11 December 2020 / Published: 15 December 2020
(This article belongs to the Special Issue Liquid Biopsy)

Abstract

:
Elevated, tumor-derived extracellular vesicle (tdEV) and circulating tumor cell (CTC) loads in metastatic cancer are associated with poor clinical outcome. Herein, we investigate whether endothelium-derived extracellular vesicles (edEVs) can be detected in the blood of metastatic colorectal cancer (mCRC) patients, and whether those vesicles associate with prognosis. The open-source ACCEPT (Automated CTC Classification, Enumeration, and Phenotyping) software was used to enumerate edEVs, tdEVs, and other objects from digitally stored CellSearch images acquired after CTC and circulating endothelial cell (CEC) enrichment from the blood of 395 mCRC patients before the initiation of a new therapy. Patients had participated in the prospective phase III CAIRO2 study. The presence of edEVs was found 5- to 10-fold higher than CECs. The hazard ratio (HR) (95% CI) of progression-free survival (PFS) for increased CTCs (≥3 in 7.5 mL), tdEVs (≥40 in 7.5 mL), and edEVs (≥287 in 4.0 mL.) was 1.4 (1.1–1.9), 2.0 (1.5–2.6), and 1.7 (1.2–2.5), respectively. The HR of Overall Survival (OS) for increased CTCs, tdEVs and edEVs was 2.2 (1.7–3.0), 2.7 (2.0–3.5), and 2.1 (1.5–2.8), respectively. There was no cut-off value for CECs, leading to a dichotomization of patients with a significant HR. Only tdEVs remained a significant predictor of OS in the final multivariable model.

1. Introduction

The necessity of non-invasive, disease-specific, and reliable biomarkers in different cancers is becoming urgent for the administration of targeted therapies and the evaluation of their efficacy in a timely fashion. In approximately 30% of advanced colorectal cancer patients, three or more circulating tumor cells (CTCs) in 7.5 mL of blood are detected, and these patients have a significantly shorter progression-free (PFS) and overall survival (OS) [1,2]. Therapy monitoring through increases or decreases in the number of CTCs is therefore only meaningful in these 30% of patients. Recently, we have shown that tumor-derived extracellular vesicles (tdEVs) are co-isolated with the CTCs and can be automatically enumerated in fluorescence images [3]. These tdEVs are present in approximately 20× higher frequencies compared to CTCs in patients with advanced prostate, breast, non-small-cell lung, and colorectal cancer, and tdEVs’ presence is also strongly associated with a poor clinical outcome [3,4,5]. The presence of circulating endothelial cells (CECs) in the blood of cancer patients has also been investigated, and has been found in higher frequencies in the blood of patients with advanced cancer compared to healthy donors [6]. The release of a portion of the detected endothelial cells can be attributed to their detachment from the blood vessel wall during the blood draw. As this varies between collected blood samples, it has a significant influence on the actual CEC counts and their reliability with regard to following disease processes [6]. We postulated that if endothelial cells produce extracellular vesicles in a similar way as tumor cells, we might be able to detect them and their count would most likely be more reliable, since their formation would be independent of the blood draw and the damage of the vasculature. For this study, we used a digitally stored set of fluorescence images generated from the blood samples of 395 advanced colorectal cancer patients receiving first-line chemotherapy and targeted agents before the initiation of a new therapy. The patients had participated in the phase III CAIRO2 study. The blood samples were processed with the CellSearch CTC kit for CTC and tdEV enumeration, as well as the CellSearch Endothelial Cell kit for CEC and edEV enumeration. The CTC and CEC data have been reported earlier [2,7]. We used the open-source ACCEPT (Automated CTC Classification, Enumeration, and Phenotyping) software to enumerate tdEVs and explore the presence of endothelium-derived extracellular vesicles (edEVs).

2. Materials and Methods

2.1. Patients

All patients included in the present analysis had participated in the prospective phase III CAIRO2 study (ClinicalTrials.gov identifier: NCT00208546) of the Dutch Colorectal Cancer Group [8]. All subjects had given their informed consent for inclusion before they participated in the study. The study was conducted in accordance with the Declaration of Helsinki, and the protocol was approved by the Committee on Research Involving Human Subjects Arnhem-Nijmegen (project identification code: CMO 2005/076).
Patients were randomly assigned to receive first-line treatment with capecitabine, oxaliplatin, and bevacizumab, or the same schedule with the addition of weekly cetuximab. Tumor response was assessed every 9 weeks using computed tomography (CT) imaging, and evaluated according to the Response Evaluation Criteria in Solid Tumors (RECIST) [9]. Details of the study have been reported previously [8,10].

2.2. Isolation of CTCs, tdEVs, CECs, and edEVs

The CellSearch CTC assay (Menarini Silicon Biosystems, Huntingdon Valley, PA, USA) was used to immunomagnetically isolate CTCs and tdEVs, based on their epithelial cell adhesion molecule (EpCAM) expression from 7.5 mL of blood collected in CellSave (Menarini Silicon Biosystems, Huntingdon Valley, PA, USA) tubes. The positively enriched objects were stained with the nuclear dye 4′,6-diamidino-2-phenylindole (DAPI), allophycocyanin-conjugated monoclonal antibodies against the leukocyte-specific cluster of differentiation 45 (CD45-APC), and phycoerythrin-conjugated monoclonal antibodies against cytokeratins 8, 18, and 19 (CK-PE). The suspension was placed in a cartridge contained within a magnet, as previously described [11].
The CellSearch Endothelial Cell assay (Menarini) was used to immunomagnetically enrich for CECs and edEVs expressing CD146 from 4.0 mL of CellSave blood. The positively enriched objects were stained with the nuclear dye DAPI and the fluorescently labeled monoclonal antibodies against CD45 (CD45-APC) and endoglin (CD105-PE). Image acquisition of both EpCAM- and CD146-enriched and stained objects within cartridges was performed on the CellTracks Analyzer II, a semi-automated fluorescence microscope equipped with computer-controlled X, Y, and Z stages, a 10× objective with 0.45 numerical aperture, a Mercury Arc lamp, a 12-bit charge-coupled device (CCD) camera, and filter cubes for DAPI, PE, APC, and fluorescein isothiocyanate (FITC). Typically, 175 images per fluorescence channel are acquired to cover the surface of the whole cartridge [11].
From this multicenter, phase III trial study, CTCs were enumerated before the initiation of therapy in 467 of the 755 advanced colorectal cancer patients and CECs in 473 of the 755 patients using the CellSearch system, and results have been reported earlier [2,7]. Fluorescence images of samples processed with the CTC and CEC kits were stored on digital versatile discs and transferred to hard disks in 450 of the 457 CTC samples and 395 of the 473 CEC samples. We re-analyzed the stored images of 395 patients for both CTC and CEC samples to identify the leukocytes, nucleated cells, and subpopulations of extracellular vesicles (EVs) co-isolated with both CTC and CEC kits.

2.3. Enumeration of tdEVs, edEVs, ldEVs, and Nucleated Cells by ACCEPT

For this study, the original CTC and CEC scores (manual counts) were used. To enumerate tdEVs, edEVs, leukocyte-derived extracellular vesicles (ldEVs), leukocytes, and nucleated cells, the digitally stored fluorescence image files were re-analyzed with the open-source Automated CTC Classification, Enumeration, and Phenotyping (ACCEPT) software (http://github.com/LeonieZ/ACCEPT), using the “Full Detection” processor. For the samples processed with the CEC kit, CD105+ leukocytes and CD105+ ldEVs were also identified and enumerated. The gates for each of the abovementioned populations in the samples processed with the CTC and CEC kits are provided in the Table S1.

2.4. Statistical Analysis

Statistical analysis was performed using the IBM SPSS platform for Windows, version 23.0 (SPSS Inc., Chicago, IL, United States) and MedCalc for Windows, version 18.0 (MedCalc Software, Ostend, Belgium). The primary objective was to assess the prognostic value of tdEVs and edEVs in mCRC patients. Eligible patients were assessable for these analyses if both a CTC and CEC sample was available prior to the start of a therapy. Patients were prospectively divided into patients with favorable and unfavorable CTC counts (<3 and ≥3, respectively). For tdEVs, the value of 40 was used as a cut-off, based on our previous results supporting the additional dichotomization of patients with favorable CTC counts using that value [5]. Overall survival (OS) was defined as the elapsed time in months between the baseline blood draw date and the date of death or last follow-up. Progression-free survival (PFS) was defined as the elapsed time in weeks between the baseline blood draw date and the date of progression or death, whichever occurred first. Patients alive at the end of the study or lost during the follow-up were censored. Median PFS and OS were defined by Kaplan–Meier survival curves. Kaplan–Meier survival curves of favorable and unfavorable groups were compared using a non-parametric log–rank test.
The open-source web application Cutoff Finder (http://molpath.charite.de/cutoff/) was used to calculate the hazard ratios (HRs) for PFS and OS, including 95% confidence intervals (CIs), over a wide range of cut-off values for CECs and edEVs of the 395 patients before the initiation of a new treatment. The optimal cut-off value of edEVs was defined as the point with the most significant split (minimum p-value, log–rank). Cutoff Finder uses the R code to provide optimization and visualization tools for cut-off determination [12].
Univariable Cox proportional hazards regression analysis was performed to determine the HR for PFS and OS, with 95% CI for each continuous and categorical variable. Continuous complete blood count (CBC)-based parameters, counts of CECs, CTCs, leukocytes, nucleated events, and EVs, isolated with the CEC and CTC kits, were log-transformed to improve the model fit. In case of the isolated events with the CTC and CEC kits, a constant value of 0.1 was added to deal with 0 values. A final multivariable Cox model was fit, including significant variables from the univariable Cox proportional hazards regression analysis. Due to correlation between some of the included variables, the final model was selected using backward stepwise elimination (pin = 0.05 and pout = 0.10).
The non-parametric, two-tailed Spearman’s Rho test was used to assess the strength of association between the available prognostic clinical parameters of the patients and the investigating biomarkers, namely CTCs, tdEVs, CECs, and edEVs. The same test was used to evaluate the strength of association between CBC-based parameters and the different classes of objects isolated with the CellSearch using the CTC and CEC kits. The non-parametric Wilcoxon signed ranks test was used to compare the distributions of matched leukocytes, ldEVs and nucleated events isolated with the different kits. Receiver operating characteristic (ROC) curves were used to assess the performance of a test using solely edEVs, CTCs, tdEVs, or a combination of edEVs with either CTCs or tdEVs. The classification variable was the shorter patient survival compared to the median OS of the patient cohort used. The areas under the curve (AUCs) were compared using the non-parametric DeLong approach [13].

3. Results

3.1. Identification of tdEVs, edEVs, ldEVs, Leukocytes, and Nucleated Cells by ACCEPT

Digitally stored fluorescence images, which were originated from the CellTracks Anayzer II, were analyzed by ACCEPT to identify tdEVs, edEVs, ldEVs, leukocytes, and nucleated cells, using the gates provided in Table S1. Figure 1 shows an example of a scatter plot for tdEVs (panel A), edEVs (panel B), ldEVs (panel C), and leukocytes (panel D). The ldEVs and leukocytes presented in this example were co-isolated using the CTC kit. For tdEVs, the CK-PE mean intensity was plotted against the DAPI mean intensity; for edEVs, the CD105-PE mean intensity was plotted against the DAPI mean intensity. For ldEVs and leukocytes, the CD45-APC mean intensity was plotted against the DAPI mean intensity. Next to each scatter plot, four images of objects falling in each gate are shown, and the position of each image in the scatter plot is indicated with numbers 1 through 4. Nucleated cells for which the lineage of origin could not be determined, due to a lack of CD45 and cytokeratin (CK)/endoglin (CD105) expression, fall on the y-axis of panels A–D.

3.2. Estimated Size Distribution of edEVs and tdEVs Derived from Fluorescence Images

The size of edEVs and tdEVs was estimated from the fluorescence images. More specifically, the PE area of each thumbnail falling in the edEV and tdEV gates (corresponding to CD105 and CKs, respectively) was measured and exported by ACCEPT. Subsequently, the size distribution plots of the herein-reported edEV and tdEV populations were constructed (Figure 2a,c). The two size distributions, even if they seem to be very similar, were significantly different (Mann–Whitney U test, p < 0.001), with tdEVs being slightly larger overall. The median radius of tdEVs detected in the fluorescence image datasets of 395 patient samples was 2.7 μm (interquartile range (IQR) = 1.5), with a minimum of 1.0 and a maximum of 6.9 μm. The median radius of edEVs was 2.5 μm (IQR = 1.2), with a minimum of 0.9 and maximum of 6.9 μm. Since some patient samples contributed more to these size distributions than others, as they contained more EVs, the respective sampling distributions were plotted, showing the distribution of the mean size of edEVs and tdEVs (Figure 2b,d). 80% of patient samples had an average edEV and tdEV radius smaller than 3 μm. The sampling distributions of the minimum and maximum sizes of edEVs and tdEVs is shown in Figure S1.
It is important to point out that even if the available fluorescence images allow for the detection of EVs, conclusions regarding their size cannot be unquestionable, due to the large pixel size (0.64 × 0.64 μm2) and the consequent lost details in the imaged objects.

3.3. Frequencies of CTCs, tdEVs, CECs, and edEVs in Colorectal Cancer Patients before the Initiation of a New Therapy

CTCs, tdEVs, CECs, and edEVs were enumerated in 395 mCRC patients after immunomagnetic enrichment using the CellSearch system. The results are summarized in Figure 3. Counts of CTCs ranged from 0 to 312 (median = 0, interquartile range (IQR) = 3), of tdEVs from 0 to 1268 (median = 30, IQR = 71), of CECs from 0 to 1250 (median = 27, IQR = 45), and of edEVs from 1 to 1504 (median = 139, IQR = 156).
In addition to tdEVs and edEVs, nucleated cells of unknown cell lineage of origin, CD45-expressing leukocytes, and ldEVs were identified and enumerated from the image sets generated after processing the blood samples with the CTC and CEC kits. In case of samples processed with the CEC kit, a portion of isolated leukocytes co-expressed CD45 and CD105, and were most likely monocytes in their differentiation transition to macrophages [14,15]. To facilitate the comparison of the frequencies of the co-isolated events between the two different kits used, all counts were normalized to the initial blood volume processed (7.5 mL for samples processed with the CTC kit and 4.0 mL for samples processed with the CEC kit), and the results are shown in Figure S2. Noteworthy is the observation that in contrast to tdEVs and edEVs, which appear in significantly higher frequencies than CTCs and CECs, respectively (p < 0.01, Wilcoxon signed ranks test), the number of ldEVs is lower than leukocytes in samples processed with either the CTC or CEC kits. The number of leukocytes and nucleated events isolated with the CEC kit are significantly higher when compared to the number of leukocytes and nucleated events isolated with the CTC kit (p < 0.01, Wilcoxon signed ranks test). In the case of the leukocytes, this observation can be attributed to the expression of CD146 on a subset of activated T-lymphocytes and natural killer (NK) cells, which are enriched with the CD146 ferrofluid aimed to enrich for endothelial cells [16,17,18]. In order to evaluate whether the leukocyte carryover is relative to the total leukocytes of the patients, as counted by a complete blood count (CBC), the Spearman’s Rho correlation coefficients ρ between all different CBC-based parameters and CTC and CEC isolated events were estimated. The results are summarized in Table S2. Even if many correlations were significant, most of them were very weak. Only CTCs were strongly associated with tdEVs.

3.4. Association of CTCs, tdEVs, CECs, and edEVs with Progression-Free and Overall Survival of Colorectal Cancer Patients before the Initiation of a New Therapy

To determine whether CECs and edEVs associate with PFS and OS, the HRs (95% CI) were calculated over a wide range of cut-off values for the CECs and edEVs of 395 patients, as illustrated in Figure S3. Only 1.8% of the possible CEC cut-off values can significantly dichotomize patients based on the HR for PFS (Figure S3A), and none of these values can result in a patient dichotomization with a significant HR for OS. On the contrary (Figure S3B), significant HRs were obtained over a wide range of edEV values (76.9% and 69.8% of all possible cutoff points, based on their PFS and OS, respectively). The optimal cut-off value of edEVs was defined as the one leading to the dichotomization of the patients with the most significant split (minimum p-value, log–rank) based on their HR for OS, which was found to be 287.
Association of CTCs with PFS and OS and of tdEVs with OS has been reported previously [1,2,5]. Here, we demonstrate the association of PFS (Figure 4) and OS (Figure 5) of 395 mCRC patients with CTCs (a), tdEVs (b), CECs (c), and edEVs (d), generating the respective Kaplan–Meier survival plots. Patients were dichotomized using a cut-off of three or more CTCs, 40 or more tdEVs, 66 or more CECs (vertical line, Figure S3A), and 287 or more edEVs (vertical line, Figure S3B). It is clear that next to CTCs and tdEVs, elevated edEVs (but not CECs) are associated with relatively worse PFS and OS. To evaluate whether edEVs could aid in the prognostication of mCRC patients with favorable and unfavorable CTC and tdEV counts, we generated ROC curves for CTCs, tdEVs, and edEVs, as well as the combined use of CTCs or tdEVs with edEVs (Figure S4), using as a classification variable the survival of patients dichotomized by the median OS. The AUC of the combined use of CTCs or tdEVs and edEVs was significantly higher (p < 0.05; DeLong) compared to the AUC of only one biomarker. Kaplan–Meier survival curves of PFS and OS were generated, stratifying patients into four different risk groups after combining their edEV and CTC (Figure 4e and Figure 5e) or edEV and tdEV (Figure 4f and Figure 5f) counts. Importantly, patients with unfavorable edEV counts and favorable CTC or tdEV counts had significantly shorter PFS and OS compared to patients with both favorable edEV and CTC or tdEV counts. In case of patients with unfavorable CTC or tdEV counts, edEV counts can further stratify patients only with respect to OS.

3.5. Univariable and Multivariable Associations between Potential Risk Factors and Clinical Outcome (PFS and OS) of mCRC Patients

All parameters, including patient characteristics (age; gender; Eastern Cooperative Oncology Group (ECOG) performance status; prior adjuvant therapy; number of metastatic sites; primary tumor in situ; sidedness of primary tumor; presence of KRAS, NRAS, and BRAF mutations; treatment arm), baseline serum lactate dehydrogenase (LDH), alkaline phosphatase (ALP), CBC-based parameters (absolute leukocyte counts, absolute platelet counts, hemoglobin), and all different classes isolated with the CTC and CEC kits (CTCs, tdEVs, CECs, edEVs, leukocytes, nucleated events, ldEVs, CD105+ leukocytes, and CD105+ ldEVs) were evaluated as potential risk factors of PFS and OS, using univariable Cox proportional hazards regression analysis. Next to the continuous log-transformed CTCs, tdEVs, CECs, and edEVs, the dichotomized CTCs, tdEVs, CECs and edEVs were also included in the analysis and are shown in Table 1. Surprisingly, leukocytes present in the cartridges after EpCAM enrichment showed a significant association with PFS and OS; the nucleated events after EpCAM enrichment were also significantly associated with OS (not PFS). In contrast to EpCAM-enriched leukocytes and nucleated events, the respective populations after CD146 enrichment were not significantly associated with worse PFS or OS.
For the multivariable analysis, the significant predictors from the univariable analysis were included. The number of metastatic sites, presence of BRAF mutation, and hemoglobin remained significant predictors of PFS in the final multivariable model (Table 2). The number of metastatic sites, presence of BRAF mutation, right-sidedness of the primary tumor, platelets (log-transformed), and tdEVs (log-transformed) remained significant predictors of OS in the final multivariable model (Table 2).

3.6. Association of Clinical Parameters with tdEVs, CTCs, CECs and edEVs

We evaluated the strength of association between clinical parameters and the biomarkers under investigation, namely CTCs, tdEVs, CECs, and edEVs, using the non-parametric Spearman’s Rho test. The results are summarized in Table S3. CECs were correlated significantly only with hemoglobin and prior adjuvant therapy, but the strength of that association was very weak (ρ < 0.2). Even if many clinical parameters were significantly correlated with edEVs, tdEVs, and CTCs, most associations were (very) weak (ρ < 0.4). Interestingly, CTCs and tdEVs had a significant association of moderate strength (0.4 < ρ < 0.6) with the LDH and ALP serum values.

4. Discussion

Liquid biopsies have attracted the interest of many research groups thanks to their potential to improve the disease management of cancer patients in a non-invasive and timely fashion. Circulating tumor cells have emerged as significant predictors of overall survival of castration-resistant prostate, metastatic breast, metastatic colorectal, and non-small-cell lung cancer patients [1,4,19,20]. Importantly, their phenotypic and genotypic characterization can shed light into new therapeutic targets and predict the treatments that will lead to the recession of specific patient groups [21,22]. However, the scarceness of CTCs in the blood, their frequently apoptotic status, and the leukocyte carryover after the chosen enrichment technique impedes their downstream genetic analysis [23,24]. Therefore, complementary research is being conducted towards the discovery of additional biomarkers predictive of treatment benefit.
Large, tumor-derived extracellular vesicles isolated together with CTCs from blood samples of metastatic prostate, breast, colorectal, and non-small cell lung cancer (NSCLC) patients are present in significantly higher frequencies than CTCs, while having an equivalent prognostic power to them [3,4,5]. Furthermore, they seem to better reflect the phenotype of the primary tumor compared to CTCs [25]. Studies of Di Vizio et al. support that large oncosomes with similar characteristics as the tdEVs we have reported in the present study (in regards to their size and expression of CKs and EpCAM), are released by non-apoptotic and highly aggressive amoeboid tumor cells [26]. Interestingly, we observed that CTCs and tdEVs have a significant association of moderate strength with the serum LDH and ALP values of mCRC patients. Elevated LDH serum levels of cancer patients have been linked to tumor hypoxia and worse clinical outcome [27,28,29,30]. It was previously shown that hypoxia stimulates the upregulation of tumor progression genes in cancer cell lines [31]. In alignment with these findings, Kallergi et al. demonstrated the expression of hypoxia-inducible factor-1α in around 60% of the total CTCs of around 75% of breast cancer patients [32]. Furthermore, increased ALP is an indicator of tumor invasiveness, since it is associated with bone and liver metastases. All these observations together support the hypothesis that hypoxia is a driving force of tumor aggressiveness, manifested by increased CTCs and tdEVs in the blood circulation of cancer patients.
Circulating endothelial cells hold great promise, since their elevated counts could be an indicator of tumor angiogenesis or vascular damage due to the tumor growth; thus, they could predict patient response to anti-angiogenic treatments. A variety of assays for CEC enumeration have been reported; however, the lack of a consensus on the CEC phenotype and of standardized techniques for CEC isolation and detection hinders the comparison of results between different studies and the drawing of incontrovertible conclusions [33]. The CEC assay that we used in the present study, immunomagnetically enriched CD146+ CECs from 4 mL of whole blood and identified CECs based on presence of nucleus, expression of CD105, and lack of CD45. CEC frequencies in blood of cancer patients are elevated compared to the CEC frequencies in healthy donors [6]. However, no association with clinical outcome could be observed in metastatic prostate cancer [34] and mCRC [7]. Our results confirm these findings for all possible cut-off values of CECs. One explanation is that CEC counts are biased due to the endothelial cells being detached during the blood draw [6]. A solution to eliminate that bias could be the use of antibodies targeting exclusively cancer-associated CECs [35]. However, in a study in which CECs were defined as DAPI+, CD34+, CD146+, and CD45−, and also expressed the immune check molecule CD276, which is frequently associated with cancer, no association could be found with poor outcome or treatment responses [36].
As tdEVs can be detected in blood samples enriched for CTCs, we asked ourselves whether edEVs could also be detected in blood samples enriched for CECs, and whether their presence is associated with clinical outcome. In contrast to CEC counts, which are influenced by a portion of endothelial cells originating from the vascular wall due to the blood draw, the edEV count is expected to be less biased.
Previous studies have reported the presence of edEVs in different pathophysiological conditions, as analyzed mainly with flow cytometry [37,38,39]. However, the need for a standardized isolation technique and a consensus on their definition has been stressed in several studies [40,41]. Duval et al. compared different antibody combinations using flow cytometry, and concluded that the CD146+, CD105+, and CD45- definition results in the most efficient detection of the edEV population [37]. Consistent with their probing, in our study CD146 was used to immunomagnetically enrich for CECs and edEVs, CD105 was used as an inclusion marker for their detection, and CD45 was an exclusion marker to discriminate them from the leukocyte carryover. We were able to detect edEVs, but the large pixel size of the acquired fluorescence images (0.64 × 0.64 μm2) let us estimate rather than accurately determine their size distribution (Figure 2). Since the smallest detected EVs in the present study, based on the immunofluorescence images, had a radius of 1 μm, questions arise regarding their biogenesis. What portion of the CellSearch isolated edEVs (and tdEVs) are apoptotic bodies and apoptotic microvesicles [42] remains to be addressed in a prospective study by the inclusion of an additional marker for apoptosis in the CellSearch runs, such as the apotracker or annexin V. A better classification of the herein-reported EVs into microvesicles and apoptotic bodies could lead to more selective biomarkers for the evaluation of the clinical outcome of the patients and their response to treatments. Moreover, the detection and further classification of edEVs and tdEVs in the plasma fraction should be explored further, as it is likely that the majority of edEVs and tdEVs will reside in this fraction. Whether or not they will have a similar relation with clinical outcome or encompass a different class of edEVs and tdEVs will have to be investigated.
Importantly, when plotting all possible cut-off values of edEV counts versus the HR, more than half of the values could significantly dichotomize the patients into more and less favorable groups of PFS (76.9% of cut-off values) and OS (69.8% of cut-off values) (Figure S3). This was into contrast with CECs, for which only 1.8% of possible cut-off values were significant in the case of only PFS. The cutoff value of 287 edEVs in 4 mL of peripheral blood was chosen as the value leading to the most significant split of the patients. The group of patients with elevated edEVs experienced significantly shorter median PFS and OS (Figure 3 and Figure 4) compared to the patients with <287 edEVs.
The weak association between edEVs and tdEVs (Table S3) suggests partially independent underlying biological functions of edEVs and tdEVs in the metastatic process of colorectal cancer. There is not much data on the biological functions of large edEVs. However, Hristov et al. reported that edEVs originated from apoptotic endothelial cells enhance the differentiation of endothelial progenitors [43]. In case of cancer, that differentiation could lead to increased angiogenesis at the tumor sites. In the present study, patients of both treatment arms received the angiogenesis inhibitor bevacizumab, which does not allow us to determine a possible effect of the drug on the edEV release. No difference in the edEV counts could be detected in patients in which the epidermal growth factor receptor (EGFR) inhibitor cetuximab was administered or not (data not shown).
Last but not least, we unexpectedly noticed that mCRC patients with increased counts of leukocytes co-isolated after EpCAM enrichment had significantly worse OS and PFS (Table 1). In contrast, no association with clinical outcome could be found between the carried-over leukocytes using the CD146 enrichment (Table 1). Previously, we reported the presence of a significantly larger number of leukocytes and nucleated cells of unknown origin after EpCAM enrichment in blood of NSCLC patients, compared to healthy donors [43]. Up to now, leukocytes detected in the CellSearch cartridges were considered to constitute a non-specific carryover. The weak correlation between the leukocyte counts in the EpCAM-enriched samples and the leukocyte counts from the available complete blood counts (Table S1) rejects a plausible explanation that the leukocyte carryover reflects the total leukocyte count of patients, with the latter known to be associated with patients’ clinical outcome in bladder, prostate, and non-small-cell lung cancers [44,45,46,47], as well as in colorectal cancer (Table 1). That observation raised questions about the subtype of leukocytes isolated with the CTC kit, which either overexpress Fcγ receptors that bind to the IgG immunoglobulins of the anti-EpCAM ferrofluid, or express EpCAM. Fc receptors are expressed on all hematopoietic cells, playing a key role in immune modulation, having both inhibitory and activating functions [48,49]. The increased expression of Fcγ receptors on the monocytes of gastric cancer patients has been previously reported [50]. An alternative plausible explanation would be the expression of EpCAM on a subset of leukocytes and their specific isolation using the CTC kit. EpCAM belongs in the CAM family, serving a diversity of cell functions [51]. Its expression is known to be restricted only in epithelial cells in healthy humans. De novo expression of EpCAM has been reported in mature human hepatocytes during various inflammatory liver diseases and liver regeneration and repair [52,53]. Furthermore, EpCAM is expressed also by a subset of thymocytes, as well as B and T lymphocytes in mice, facilitating the cells’ infiltration in inflammation sites [54,55]. Similarly, the upregulation of integrins, which comprise a CAM family, by effector T lymphocytes also facilitate their targeted access to inflammatory sites [56]. However, to our knowledge, it has not been investigated whereas EpCAM is expressed in any subsets of leukocytes in cancer patients that would facilitate the leukocytes’ infiltration in the tumor. The EpCAM expression on the surfaces of the leukocyte could emanate from (a) the activation of originally silent molecular pathways, (b) the uptake of tumor-derived extracellular vesicles by cell–EV membrane fusion [57], or (c) the internalization and expression of functional EpCAM mRNA by leukocytes [58]. The two latter suggestions are supported by the findings that tdEVs accumulate in all leukocyte subpopulations in murine models, altering the phenotypes of the recipient cells [59,60].

5. Conclusions

Endothelium-derived extracellular vesicles, isolated based on their CD146 expression and detected based on their CD105 expression, are significant predictors of PFS and OS in metastatic colorectal cancer, in contrast to circulating endothelial cells. Interestingly, the EpCAM-based co-isolated leukocytes constitute a significant predictor of OS and PFS, suggesting further investigation on the phenotype of that population.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4409/9/12/2688/s1, Figure S1: Sampling distributions of mininum edEV size (a), maximum edEV size (b), minimum tdEV size (c), and maximum tdEV size (d). The black lines in each graph correspond to the respective cumulative distributions. Figure S2: Frequencies of different populations co-isolated with the CTC (black dots) and CEC kits (grey dots), normalized to 1 mL of blood for comparison between the kits. Figure S3: Cut-off optimization of the baseline values for CECs (A) and edEVs (B) in mCRC patients. For each possible cut-off, CECs and edEVs were correlated with PFS (top) or overall survival (OS) (bottom). The hazard ratio (HR), including 95% CI, was plotted depending on the cutoff. The vertical line indicates the cut-off that results in the most significant correlation with OS. The value distribution of CECs and edEVs is shown as a rug plot at the bottom of the respective figure. Figure S4: Receiver operating characteristic (ROC) curves treating survival time dichotomized by the median OS time of the patient cohort as the classification variable. The addition of edEVs to CTCs (A) or tdEVs (B) results in a significantly (p < 0.05) higher area under the curve (AUC) compared to solely CTCs or tdEVs. Table S1: ACCEPT gates for the automated enumeration of tdEVs, edEVs, ldEVs, leukocytes, and nucleated events. The exponents a and b correspond to the CTC and CEC kit, respectively. Table S2: Correlation between CBC-based parameters, CTC and CEC kit-isolated objects of 395 mCRC patients, using the Spearman’s Rho correlation coefficient (ρ). Correlation is considered to be weak at ρ < 0.4 (in grey), moderate at 0.4 ≤ ρ <0.6 (in black), and strong at 0.6 ≤ ρ ≤0.8 (in bold black). * Indicates significance at 0.05 level (two-tailed); ** indicates significance at 0.01 level (two-tailed). Table S3: Correlation between clinical parameters and CTCs, tdEVs, CECs, and edEVs of 395 mCRC patients using the Spearman’s Rho correlation coefficient (ρ). Correlation is considered to be weak at ρ < 0.4 (in grey), moderate at 0.4 ≤ ρ <0.6 (in black), and strong at 0.6 ≤ ρ ≤0.8 (in bold black). ** Indicates significance at 0.01 level (two-tailed); * indicates significance at 0.05 level (two-tailed).

Author Contributions

Conceptualization, A.N. and L.W.M.M.T.; methodology, A.N. and L.W.M.M.T.; validation, A.N.; formal analysis, A.N.; investigation, A.N.; resources, L.M., F.A.W.C., M.K., C.J.A.P., and L.W.M.M.T.; writing—original draft preparation, A.N. and L.W.M.M.T.; writing—review and editing, A.N., L.M., F.A.W.C., M.K., C.J.A.P., and L.W.M.M.T.; supervision, L.W.M.M.T.; project administration, L.W.M.M.T.; funding acquisition, L.W.M.M.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Netherlands Organization for Scientific Research (NWO) Applied and Engineering Science Cancer-ID project no. 14190.

Acknowledgments

The authors acknowledge all patients and healthy donors, the blood samples of whom were used for the accomplishment of the present study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cohen, S.J.; Punt, C.J.; Iannotti, N.; Saidman, B.H.; Sabbath, K.D.; Gabrail, N.Y.; Picus, J.; Morse, M.; Mitchell, E.; Miller, M.C.; et al. Relationship of circulating tumor cells to tumor response, progression-free survival, and overall survival in patients with metastatic colorectal cancer. J. Clin. Oncol. 2008, 26, 3213–3221. [Google Scholar] [CrossRef] [PubMed]
  2. Tol, J.; Koopman, M.; Miller, M.C.; Tibbe, A.; Cats, A.; Creemers, G.J.M.; Vos, A.H.; Nagtegaal, I.D.; Terstappen, L.W.M.M.; Punt, C.J.A. Circulating tumour cells early predict progression-free and overall survival in advanced colorectal cancer patients treated with chemotherapy and targeted agents. Ann. Oncol. Off. J. Eur. Soc. Med Oncol. 2010, 21, 1006–1012. [Google Scholar] [CrossRef]
  3. Nanou, A.; Coumans, F.A.W.; van Dalum, G.; Zeune, L.L.; Dolling, D.; Onstenk, W.; Crespo, M.; Fontes, M.S.; Rescigno, P.; Fowler, G.; et al. Circulating tumor cells, tumor-derived extracellular vesicles and plasma cytokeratins in castration-resistant prostate cancer patients. Oncotarget 2018, 9, 19283–19293. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. De Wit, S.; Rossi, E.; Weber, S.; Tamminga, M.; Manicone, M.; Swennenhuis, J.F.; Groothuis-Oudshoorn, C.G.M.; Vidotto, R.; Facchinetti, A.; Zeune, L.L.; et al. Single tube liquid biopsy for advanced non-small cell lung cancer. Int. J. Cancer 2018. [Google Scholar] [CrossRef] [Green Version]
  5. Nanou, A.; Miller, M.C.; Zeune, L.L.; de Wit, S.; Punt, C.J.A.; Groen, H.J.M.; Hayes, D.F.; de Bono, J.S.; Terstappen, L. Tumour-derived extracellular vesicles in blood of metastatic cancer patients associate with overall survival. Br. J. Cancer 2020. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Rowand, J.L.; Martin, G.; Doyle, G.V.; Miller, M.C.; Pierce, M.S.; Connelly, M.C.; Rao, C.; Terstappen, L.W. Endothelial cells in peripheral blood of healthy subjects and patients with metastatic carcinomas. Cytometry A 2007, 71, 105–113. [Google Scholar] [CrossRef]
  7. Simkens, L.H.; Tol, J.; Terstappen, L.W.; Teerenstra, S.; Punt, C.J.; Nagtegaal, I.D. The predictive and prognostic value of circulating endothelial cells in advanced colorectal cancer patients receiving first-line chemotherapy and bevacizumab. Ann. Oncol. 2010, 21, 2447–2448. [Google Scholar] [CrossRef]
  8. Tol, J.; Koopman, M.; Cats, A.; Rodenburg, C.J.; Creemers, G.J.M.; Schrama, J.G.; Erdkamp, F.L.G.; Vos, A.H.; van Groeningen, C.J.; Sinnige, H.A.M.; et al. Chemotherapy, bevacizumab, and cetuximab in metastatic colorectal cancer. N. Engl. J. Med. 2009, 360, 563–572. [Google Scholar] [CrossRef] [Green Version]
  9. Therasse, P.; Arbuck, S.G.; Eisenhauer, E.A.; Wanders, J.; Kaplan, R.S.; Rubinstein, L.; Verweij, J.; Van Glabbeke, M.; van Oosterom, A.T.; Christian, M.C.; et al. New guidelines to evaluate the response to treatment in solid tumors. European organization for research and treatment of cancer, national cancer institute of the United States, national cancer institute of Canada. J. Natl. Cancer Inst. 2000, 92, 205–216. [Google Scholar] [CrossRef] [Green Version]
  10. Tol, J.; Koopman, M.; Rodenburg, C.J.; Cats, A.; Creemers, G.J.; Schrama, J.G.; Erdkamp, F.L.G.; Vos, A.H.; Mol, L.; Antonini, N.F.; et al. A randomised phase III study on capecitabine, oxaliplatin and bevacizumab with or without cetuximab in first-line advanced colorectal cancer, the CAIRO2 study of the Dutch Colorectal Cancer Group (DCCG). An interim analysis of toxicity. Ann. Oncol. J. Eur. Soc. Med Oncol. 2008, 19, 734–738. [Google Scholar] [CrossRef]
  11. Coumans, F.; Terstappen, L. detection and characterization of circulating tumor cells by the cellsearch approach. Methods Mol Biol. 2015, 1347, 263–278. [Google Scholar] [CrossRef] [PubMed]
  12. Budczies, J.; Klauschen, F.; Sinn, B.V.; Gyorffy, B.; Schmitt, W.D.; Darb-Esfahani, S.; Denkert, C. Cutoff finder: A comprehensive and straightforward Web application enabling rapid biomarker cutoff optimization. PLoS ONE. 2012, 7, e51862. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. DeLong, E.R.; DeLong, D.M.; Clarke-Pearson, D.L. Comparing the areas under two or more correlated receiver operating characteristic curves: A nonparametric approach. Biometrics 1988, 44, 837–845. [Google Scholar] [CrossRef] [PubMed]
  14. O’Connell, P.J.; McKenzie, A.; Fisicaro, N.; Rockman, S.P.; Pearse, M.J.; d’Apice, A.J. Endoglin: A 180-kD endothelial cell and macrophage restricted differentiation molecule. Clin Exp. Immunol. 1992, 90, 154–159. [Google Scholar] [CrossRef]
  15. Lastres, P.; Bellon, T.; Cabanas, C.; Sanchez-Madrid, F.; Acevedo, A.; Gougos, A.; Letarte, M.; Bernabeu, C. Regulated expression on human macrophages of endoglin, an Arg-Gly-Asp-containing surface antigen. Eur. J. Immunol. 1992, 22, 393–397. [Google Scholar] [CrossRef]
  16. Dignat-George, F.; Sabatier, F.; Blann, A.; Woywodt, A. Detection of circulating endothelial cells: CD146-based magnetic separation enrichment or flow cytometric assay? J. Clin. Oncol. 2007, 25, e1–e2. [Google Scholar] [CrossRef]
  17. Elshal, M.F.; Khan, S.S.; Takahashi, Y.; Solomon, M.A.; McCoy, J.P. CD146 (Mel-CAM), an adhesion marker of endothelial cells, is a novel marker of lymphocyte subset activation in normal peripheral blood. Blood 2005, 106, 2923–2924. [Google Scholar] [CrossRef]
  18. Pickl, W.F.; Majdic, O.; Fischer, G.F.; Petzelbauer, P.; Fae, I.; Waclavicek, M.; Stockl, J.; Scheinecker, C.; Vidicki, T.; Aschauer, H.; et al. MUC18/MCAM (CD146), an activation antigen of human T lymphocytes. J. Immunol. 1997, 158, 2107–2115. [Google Scholar]
  19. Cristofanilli, M.; Budd, G.T.; Ellis, M.J.; Stopeck, A.; Matera, J.; Miller, M.C.; Reuben, J.M.; Doyle, G.V.; Allard, W.J.; Terstappen, L.W.; et al. Circulating tumor cells, disease progression, and survival in metastatic breast cancer. N. Engl. J. Med. 2004, 351, 781–791. [Google Scholar] [CrossRef] [Green Version]
  20. Scher, H.I.; Heller, G.; Molina, A.; Attard, G.; Danila, D.C.; Jia, X.; Peng, W.; Sandhu, S.K.; Olmos, D.; Riisnaes, R.; et al. Circulating tumor cell biomarker panel as an individual-level surrogate for survival in metastatic castration-resistant prostate cancer. J. Clin. Oncol. 2015, 33, 1348–1355. [Google Scholar] [CrossRef]
  21. Bidard, F.C.; Fehm, T.; Ignatiadis, M.; Smerage, J.B.; Alix-Panabieres, C.; Janni, W.; Messina, C.; Paoletti, C.; Muller, V.; Hayes, D.F.; et al. Clinical application of circulating tumor cells in breast cancer: Overview of the current interventional trials. Cancer Metastasis Rev. 2013, 32, 179–188. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. El-Heliebi, A.; Hille, C.; Laxman, N.; Svedlund, J.; Haudum, C.; Ercan, E.; Kroneis, T.; Chen, S.; Smolle, M.; Rossmann, C. In situ detection and quantification of AR-V7, AR-FL, PSA, and KRAS point mutations in circulating tumor cells. Clin. Chem. 2018, 64, 536–546. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Swennenhuis, J.F.; Reumers, J.; Thys, K.; Aerssens, J.; Terstappen, L.W. Efficiency of whole genome amplification of single circulating tumor cells enriched by CellSearch and sorted by FACS. Genome Med. 2013, 5, 106. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Coumans, F.A.; Ligthart, S.T.; Uhr, J.W.; Terstappen, L.W. Challenges in the enumeration and phenotyping of CTC. Clin. Cancer Res. 2012, 18, 5711–5718. [Google Scholar] [CrossRef] [Green Version]
  25. Nanou, A.; Zeune, L.L.; Bidard, F.-C.; Pierga, J.-Y. Terstappen, L.W.M.M. HER2 expression on tumor-derived extracellular vesicles and circulating tumor cells in metastatic breast cancer. Breast Cancer Res. 2020, 22, 1–11. [Google Scholar] [CrossRef]
  26. Ciardiello, C.; Cavallini, L.; Spinelli, C.; Yang, J.; Reis-Sobreiro, M.; de Candia, P.; Minciacchi, V.R.; Di Vizio, D. Focus on extracellular vesicles: New frontiers of cell-to-cell communication in cancer. Int. J. Mol. Sci. 2016, 17, 175. [Google Scholar] [CrossRef] [Green Version]
  27. Koukourakis, M.I.; Giatromanolaki, A.; Simopoulos, C.; Polychronidis, A.; Sivridis, E. Lactate dehydrogenase 5 (LDH5) relates to up-regulated hypoxia inducible factor pathway and metastasis in colorectal cancer. Clinical Exp. Metastasis 2005, 22, 25–30. [Google Scholar] [CrossRef]
  28. Jurisic, V.; Radenkovic, S.; Konjevic, G. The actual role of LDH as tumor marker, biochemical and clinical aspects. Adv. Exp. Med. Biol. 2015, 867, 115–124. [Google Scholar]
  29. Koukourakis, M.I.; Giatromanolaki, A.; Sivridis, E.; Bougioukas, G.; Didilis, V.; Gatter, K.; Harris, A. Lactate dehydrogenase-5 (LDH-5) overexpression in non-small-cell lung cancer tissues is linked to tumour hypoxia, angiogenic factor production and poor prognosis. Br. J. Cancer 2003, 89, 877–885. [Google Scholar] [CrossRef]
  30. Kolev, Y.; Uetake, H.; Takagi, Y.; Sugihara, K. Lactate dehydrogenase-5 (LDH-5) expression in human gastric cancer: Association with hypoxia-inducible factor (HIF-1α) pathway, angiogenic factors production and poor prognosis. Ann. Surg. Oncol. 2008, 15, 2336–2344. [Google Scholar] [CrossRef]
  31. Zweitzig, D.R.; Smirnov, D.A.; Connelly, M.C.; Terstappen, L.W.; O’Hara, S.M.; Moran, E. Physiological stress induces the metastasis marker AGR2 in breast cancer cells. Mol. Cell. Biochem. 2007, 306, 255–260. [Google Scholar] [CrossRef] [PubMed]
  32. Kallergi, G.; Markomanolaki, H.; Giannoukaraki, V.; Papadaki, M.A.; Strati, A.; Lianidou, E.S.; Georgoulias, V.; Mavroudis, D.; Agelaki, S. Hypoxia-inducible factor-1α and vascular endothelial growth factor expression in circulating tumor cells of breast cancer patients. Breast Cancer Res. 2009, 11, R84. [Google Scholar] [CrossRef] [PubMed]
  33. Strijbos, M.H.; Gratama, J.W.; Kraan, J.; Lamers, C.H.; den Bakker, M.A.; Sleijfer, S. Circulating endothelial cells in oncology: Pitfalls and promises. Br. J. Cancer 2008, 98, 1731–1735. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Strijbos, M.H.; Gratama, J.W.; Schmitz, P.I.; Rao, C.; Onstenk, W.; Doyle, G.V.; Miller, M.C.; de Wit, R.; Terstappen, L.W.; Sleijfer, S. Circulating endothelial cells, circulating tumour cells, tissue factor, endothelin-1 and overall survival in prostate cancer patients treated with docetaxel. Eur. J. Cancer 2010, 46, 2027–2035. [Google Scholar] [CrossRef]
  35. Smirnov, D.A.; Foulk, B.W.; Doyle, G.V.; Connelly, M.C.; Terstappen, L.W.M.M.; O’Hara, S.M. Global gene expression profiling of circulating endothelial cells in patients with metastatic carcinomas. Cancer Res. 2006, 66, 2918–2922. [Google Scholar] [CrossRef] [Green Version]
  36. Gootjes, E.C.; Kraan, J.; Buffart, T.; Verhoef, C.; Verheul, H.M.; Sleijfer, S.; Group, O.S. CD276-positive circulating endothelial cells in advanced colorectal cancer. Am. Soc. Clin. Oncol. 2019, 37, 572. [Google Scholar] [CrossRef]
  37. Duval, A.; Helley, D.; Capron, L.; Youinou, P.; Renaudineau, Y.; Dubucquoi, S.; Fischer, A.M.; Hachulla, E. Endothelial dysfunction in systemic lupus patients with low disease activity: Evaluation by quantification and characterization of circulating endothelial microparticles, role of anti-endothelial cell antibodies. Rheumatology 2010, 49, 1049–1055. [Google Scholar] [CrossRef] [Green Version]
  38. Dignat-George, F.; Boulanger, C.M. The many faces of endothelial microparticles. Arter. Thromb. Vasc. Biol. 2011, 31, 27–33. [Google Scholar] [CrossRef] [Green Version]
  39. Hromada, C.; Muhleder, S.; Grillari, J.; Redl, H.; Holnthoner, W. Endothelial extracellular vesicles-promises and challenges. Front Physiol. 2017, 8, 275. [Google Scholar] [CrossRef]
  40. Jy, W.; Horstman, L.L.; Jimenez, J.J.; Ahn, Y.S.; Biro, E.; Nieuwland, R.; Sturk, A.; Dignat-George, F.; Sabatier, F.; Camoin-Jau, L.; et al. Measuring circulating cell-derived microparticles. J. Thromb Haemost. 2004, 2, 1842–1851. [Google Scholar] [CrossRef]
  41. Lacroix, R.; Robert, S.; Poncelet, P.; Dignat-George, F. Overcoming limitations of microparticle measurement by flow cytometry. Semin. Thromb Hemost. 2010, 36, 807–818. [Google Scholar] [CrossRef] [PubMed]
  42. Caruso, S.; Poon, I.K. Apoptotic cell-derived extracellular vesicles: More than just debris. Front. Immunol. 2018, 9, 1486. [Google Scholar] [CrossRef] [Green Version]
  43. De Wit, S.; Zeune, L.L.; Hiltermann, T.J.N.; Groen, H.J.M.; Dalum, G.V.; Terstappen, L. Classification of cells in CTC-enriched samples by advanced image analysis. Cancers 2018, 10, 377. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Bhindi, B.; Hermanns, T.; Wei, Y.; Yu, J.; Richard, P.O.; Wettstein, M.S.; Templeton, A.; Li, K.; Sridhar, S.S.; Jewett, M.A.; et al. Identification of the best complete blood count-based predictors for bladder cancer outcomes in patients undergoing radical cystectomy. Br. J. Cancer. 2016, 114, 207–212. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Taussky, D.; Soulieres, D.; Azoulay, L.; Yin, H.; Bahig, H.; Bahary, J.P.; Delouya, G. A combination of testosterone and white blood cell count as a predictive factor of overall survival in localized prostate cancer. Target Oncol. 2017, 12, 695–701. [Google Scholar] [CrossRef] [PubMed]
  46. Tibaldi, C.; Vasile, E.; Bernardini, I.; Orlandini, C.; Andreuccetti, M.; Falcone, A. Baseline elevated leukocyte count in peripheral blood is associated with poor survival in patients with advanced non-small cell lung cancer: A prognostic model. J. Cancer Res. Clin. Oncol. 2008, 134, 1143–1149. [Google Scholar] [CrossRef]
  47. Mandrekar, S.J.; Schild, S.E.; Hillman, S.L.; Allen, K.L.; Marks, R.S.; Mailliard, J.A.; Krook, J.E.; Maksymiuk, A.W.; Chansky, K.; Kelly, K.; et al. A prognostic model for advanced stage nonsmall cell lung cancer. Pooled analysis of north central cancer treatment group trials. Cancer 2006, 107, 781–792. [Google Scholar] [CrossRef]
  48. Bournazos, S.; Wang, T.T.; Ravetch, J.V. The role and function of fcgamma receptors on myeloid cells. Microbiol. Spectr. 2016, 4. [Google Scholar] [CrossRef] [Green Version]
  49. Pincetic, A.; Bournazos, S.; DiLillo, D.J.; Maamary, J.; Wang, T.T.; Dahan, R.; Fiebiger, B.M.; Ravetch, J.V. Type I and type II Fc receptors regulate innate and adaptive immunity. Nat. Immunol. 2014, 15, 707–716. [Google Scholar] [CrossRef]
  50. Uracz, W.; Pituch-Noworolska, A.; Popiela, T.; Zembala, M. The Fc receptors of normal and cancer patients monocytes. Immunobiology 1982, 163, 450–459. [Google Scholar] [CrossRef]
  51. Trzpis, M.; McLaughlin, P.M.; de Leij, L.M.; Harmsen, M.C. Epithelial cell adhesion molecule: More than a carcinoma marker and adhesion molecule. Am. J. Pathol. 2007, 171, 386–395. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Breuhahn, K.; Baeuerle, P.A.; Peters, M.; Prang, N.; Tox, U.; Kohne-Volland, R.; Dries, V.; Schirmacher, P.; Leo, E. Expression of epithelial cellular adhesion molecule (Ep-CAM) in chronic (necro-)inflammatory liver diseases and hepatocellular carcinoma. Hepatol Res. 2006, 34, 50–56. [Google Scholar] [CrossRef] [PubMed]
  53. De Boer, C.J.; van Krieken, J.H.; Janssen-van Rhijn, C.M.; Litvinov, S.V. Expression of Ep-CAM in normal, regenerating, metaplastic, and neoplastic liver. J. Pathol. 1999, 188, 201–206. [Google Scholar] [CrossRef]
  54. Amann, M.; Friedrich, M.; Lutterbuese, P.; Vieser, E.; Lorenczewski, G.; Petersen, L.; Brischwein, K.; Kufer, P.; Kischel, R.; Baeuerle, P.A.; et al. Therapeutic window of an EpCAM/CD3-specific BiTE antibody in mice is determined by a subpopulation of EpCAM-expressing lymphocytes that is absent in humans. Cancer Immunol. Immunother 2009, 58, 95–109. [Google Scholar] [CrossRef]
  55. Nelson, A.J.; Dunn, R.J.; Peach, R.; Aruffo, A.; Farr, A.G. The murine homolog of human Ep-CAM, a homotypic adhesion molecule, is expressed by thymocytes and thymic epithelial cells. Eur. J. Immunol. 1996, 26, 401–408. [Google Scholar] [CrossRef]
  56. Dailey, M.O. Expression of T iymphocyte adhesion molecules: Regulation during antigen-induced T cell activation and differentiation. Crit. Rev. Immunol. 1998, 18, 153–184. [Google Scholar] [CrossRef]
  57. McKelvey, K.J.; Powell, K.L.; Ashton, A.W.; Morris, J.M.; McCracken, S.A. Exosomes: mechanisms of uptake. J. Circ. Biomark. 2015, 4, 7. [Google Scholar] [CrossRef] [Green Version]
  58. Valadi, H.; Ekström, K.; Bossios, A.; Sjöstrand, M.; Lee, J.J.; Lötvall, J.O. Exosome-mediated transfer of mRNAs and microRNAs is a novel mechanism of genetic exchange between cells. Nat. Cell Biol. 2007, 9, 654. [Google Scholar] [CrossRef] [Green Version]
  59. Chennakrishnaiah, S.; Meehan, B.; D’Asti, E.; Montermini, L.; Lee, T.H.; Karatzas, N.; Buchanan, M.; Tawil, N.; Choi, D.; Divangahi, M. Leukocytes as a reservoir of circulating oncogenic DNA and regulatory targets of tumor-derived extracellular vesicles. J. Thromb. Haemostasis. 2018, 16, 1800–1813. [Google Scholar] [CrossRef] [Green Version]
  60. Zech, D.; Rana, S.; Buchler, M.W.; Zoller, M. Tumor-exosomes and leukocyte activation: An ambivalent crosstalk. Cell Commun. Signal. 2012, 10, 37. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Examples of tumor-derived extracellular vesicles (tdEVs) (a), endothelium-derived extracellular vesicles (edEVs) (b), leukocyte-derived extracellular vesicles (ldEVs) (c), and leukocytes (d), isolated with the CellSearch system and identified by the ACCEPT (Automated CTC (circulating tumor cell) Classification, Enumeration, and Phenotyping) software. The dots corresponding to the objects that fall within the defined gates are depicted in blue in the respective scatter plots, whereas all other objects are depicted in grey. For each class, examples of four thumbnails with different fluorescence intensities (1–4) are shown. Cytokeratin (CK) (a) and endoglin (CD105) (b) are depicted in green, CD45 in red, and 4′,6-diamidino-2-phenylindole (DAPI) in blue (c and d). Scale bars indicate 6.4 μm.
Figure 1. Examples of tumor-derived extracellular vesicles (tdEVs) (a), endothelium-derived extracellular vesicles (edEVs) (b), leukocyte-derived extracellular vesicles (ldEVs) (c), and leukocytes (d), isolated with the CellSearch system and identified by the ACCEPT (Automated CTC (circulating tumor cell) Classification, Enumeration, and Phenotyping) software. The dots corresponding to the objects that fall within the defined gates are depicted in blue in the respective scatter plots, whereas all other objects are depicted in grey. For each class, examples of four thumbnails with different fluorescence intensities (1–4) are shown. Cytokeratin (CK) (a) and endoglin (CD105) (b) are depicted in green, CD45 in red, and 4′,6-diamidino-2-phenylindole (DAPI) in blue (c and d). Scale bars indicate 6.4 μm.
Cells 09 02688 g001
Figure 2. Size distribution of edEVs (a) and tdEVs (c) detected in 395 patient samples, derived from the respective CellSearch fluorescence images using the ACCEPT software. Sampling distributions of the mean edEV (b) and tdEV (d) size of 395 patient samples. The black lines correspond to the respective cumulative distributions.
Figure 2. Size distribution of edEVs (a) and tdEVs (c) detected in 395 patient samples, derived from the respective CellSearch fluorescence images using the ACCEPT software. Sampling distributions of the mean edEV (b) and tdEV (d) size of 395 patient samples. The black lines correspond to the respective cumulative distributions.
Cells 09 02688 g002
Figure 3. Frequencies of CellSearch isolated CTCs and tdEVs (in black) from 7.5 mL of blood, and CECs and edEVs (in grey) from 4.0 mL of blood of 395 mCRC patients before the initiation of a new therapy. The horizontal lines correspond to median values.
Figure 3. Frequencies of CellSearch isolated CTCs and tdEVs (in black) from 7.5 mL of blood, and CECs and edEVs (in grey) from 4.0 mL of blood of 395 mCRC patients before the initiation of a new therapy. The horizontal lines correspond to median values.
Cells 09 02688 g003
Figure 4. Association of CTCs (a), tdEVs (b), CECs (c), edEVs (d), combined CTCs and edEVs (e), and combined tdEVs and edEVs (f) with progression-free survival (PFS) of 323 metastatic colorectal cancer (mCRC) patients before the initiation of a new therapy.
Figure 4. Association of CTCs (a), tdEVs (b), CECs (c), edEVs (d), combined CTCs and edEVs (e), and combined tdEVs and edEVs (f) with progression-free survival (PFS) of 323 metastatic colorectal cancer (mCRC) patients before the initiation of a new therapy.
Cells 09 02688 g004
Figure 5. Association of CTCs (a), tdEVs (b), CECs (c), edEVs (d), combined CTCs and edEVs (e), and combined tdEVs and edEVs (f) with overall survival of 395 mCRC patients before the initiation of a new therapy.
Figure 5. Association of CTCs (a), tdEVs (b), CECs (c), edEVs (d), combined CTCs and edEVs (e), and combined tdEVs and edEVs (f) with overall survival of 395 mCRC patients before the initiation of a new therapy.
Cells 09 02688 g005
Table 1. Univariable Cox proportional hazards regression analysis.
Table 1. Univariable Cox proportional hazards regression analysis.
VariableCut-OffPFS RiskOS Risk
HR (95% CI)p-ValueHR (95% CI)p-Value
Prior adjuvant therapy 0.7 (0.4–1.0)0.0430.7 (0.4–1.0)0.046
Number of metastatic sites 1.3 (1.2–1.5)<0.0011.3 (1.2–1.5)< 0.001
Primary tumor in situ 1.5 (1.1–2.0)0.0131.8 (1.3–2.4)< 0.001
Presence of KRAS mutation 1.0 (0.7–1.3)0.8391.1 (0.8–1.5)0.666
Presence of BRAF mutation 3.0 (2.0–4.6)<0.0012.7 (1.8–4.0)< 0.001
Presence of NRAS mutation 1.19 (0.60–2.40)0.6291.0 (0.5–2.0)1.000
Right-sidedness of primary tumor (transverse colon, hepatic flexure, ascending colon, cecum) 1.3 (1.0–1.8)0.0421.6 (1.2–2.1)0.001
Treatment arm 1.1 (0.9–1.4)0.4651.2 (0.9–1.5)0.313
ECOG performance status 1.4 (1.1–1.8)0.0121.4 (1.1–1.9)0.006
age 1.0 (1.0–1.0)0.0781.0 (1.0–1.0)0.290
gender 1.0 (0.7–1.3)0.9581.0 (0.7–1.3)0.773
CellSearch Populations
CTCs (log transformed)-1.2 (1.1–1.4)0.0041.6 (1.4–1.9)<0.001
≥ 3 CTCs31.4 (1.1–1.9)0.0122.2 (1.7–3.0)<0.001
tdEVs (log-transformed)-1.8 (1.4–2.3)<0.0012.5 (2.0–3.2)<0.001
≥ 40 tdEVs402.0 (1.5–2.6)<0.0012.7 (2.0–3.5)<0.001
CECs (log-transformed)-1.1 (0.9–1.5)0.3381.0 (0.8–1.2)0.767
≥ 66 CECs660.9 (0.7–1.3)0.7030.7 (0.5–1.0)0.051
edEVs (log-transformed)-1.7 (1.2–2.5)0.0041.9 (1.3–2.8)0.002
≥ 287 edEVs2871.8 (1.3–2.5)<0.0012.1 (1.5–2.8)<0.001
Other CellSearch Populations
Leukocytes a
(log-transformed)
-1.3 (1.1–1.5)0.0091.3 (1.1–1.6)0.001
ldEVs a
(log-transformed)
-1.1 (0.9–1.4)0.4371.2 (0.9–1.5)0.219
Nucleated a
(log-transformed)
-1.1 (0.9–1.4)0.2391.3 (1.0–1.6)0.048
Leukocytes b
(log-transformed)
-1.1 (0.8–1.5)0.7371.1 (0.8–1.5)0.751
ldEVs b
(log-transformed)
-1.0 (0.7–1.4)0.9761.1 (0.8–1.6)0.608
Nucleated b
(log-transformed)
-1.1 (0.9–1.5)0.4361.3 (1.0–1.7)0.091
CD105+ leukocytes b
(log-transformed)
-1.0 (0.9–1.2)0.7891.0 (0.8–1.2)0.927
CD105+ ldEVs b
(log-transformed)
-1.0 (0.8–1.2)0.8421.1 (0.9–1.4)0.203
CBC-Based Parameters
Leukocytes
(log-transformed)
-6.0 (2.1–16.9)0.0019.5 (3.5–25.9)<0.001
Platelets
(log-transformed)
-4.0 (1.6–10.3)0.0045.3 (2.1–13.4)<0.001
Hemoglobin, per 1 g/dL-0.8 (0.7–0.9)<0.0010.7 (0.6–0.8)<0.001
Alkaline phosphatase (ALP)
(log-transformed)
2.7 (1.6–4.4)<0.0013.8 (2.4–6.2)<0.001
Baseline serum lactate Dehydrogenase (LDH)
(log-transformed)
1.9 (1.3–2.8)0.0022.4 (1.6–3.5)<0.001
a Refers to populations co-isolated with the CTC kit; b refers to populations co-isolated with the circulating endothelial cell (CEC) kit.
Table 2. Multivariable Cox proportional hazards regression analysis (using backward and forward stepwise elimination).
Table 2. Multivariable Cox proportional hazards regression analysis (using backward and forward stepwise elimination).
Variable in EquationPFS RiskOS Risk
HR (95% CI)p-ValueHR (95% CI)p-Value
Presence of BRAF mutation 2.8 (1.8–4.4)< 0.0012.3 (1.4–3.6)<0.001
tdEVs (log-transformed) 1.8 (1.3–2.4)0.001
Number of metastatic sites1.4 (1.2–1.7)< 0.0011.3 (1.1–1.5)0.003
Right-sidedness of the primary tumor 1.6 (1.1–2.3)0.010
Platelets (log transformed) 4.6 (1.4–15.1)0.013
Hemoglobin, per 1 g/dL0.8 (0.7–1.0)0.048
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Nanou, A.; Mol, L.; Coumans, F.A.W.; Koopman, M.; Punt, C.J.A.; Terstappen, L.W.M.M. Endothelium-Derived Extracellular Vesicles Associate with Poor Prognosis in Metastatic Colorectal Cancer. Cells 2020, 9, 2688. https://doi.org/10.3390/cells9122688

AMA Style

Nanou A, Mol L, Coumans FAW, Koopman M, Punt CJA, Terstappen LWMM. Endothelium-Derived Extracellular Vesicles Associate with Poor Prognosis in Metastatic Colorectal Cancer. Cells. 2020; 9(12):2688. https://doi.org/10.3390/cells9122688

Chicago/Turabian Style

Nanou, Afroditi, Linda Mol, Frank A. W. Coumans, Miriam Koopman, Cornelis J. A. Punt, and Leon W. M. M. Terstappen. 2020. "Endothelium-Derived Extracellular Vesicles Associate with Poor Prognosis in Metastatic Colorectal Cancer" Cells 9, no. 12: 2688. https://doi.org/10.3390/cells9122688

APA Style

Nanou, A., Mol, L., Coumans, F. A. W., Koopman, M., Punt, C. J. A., & Terstappen, L. W. M. M. (2020). Endothelium-Derived Extracellular Vesicles Associate with Poor Prognosis in Metastatic Colorectal Cancer. Cells, 9(12), 2688. https://doi.org/10.3390/cells9122688

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop