HER2 expression on tumor-derived extracellular vesicles and circulating tumor cells in metastatic breast cancer.

BACKGROUND: Tumor-derived extracellular vesicles (tdEVs) and circulating tumor cells (CTCs) in the blood of metastatic cancer patients associate with poor outcomes. In this study, we explored the human epidermal growth factor receptor 2 (HER2) expression on CTCs and tdEVs of metastatic breast cancer patients.
METHODS: Blood samples from 98 patients (CLCC-IC-2006-04 study) were originally processed with the CellSearch® system using the CTC kit and anti-HER2 as an additional marker in the staining cocktail. CTCs and tdEVs were automatically enumerated from the generated CellSearch images using the open-source ACCEPT software.
RESULTS: CTCs and tdEVs were subdivided based on their cytokeratin (CK) and HER2 phenotype into CK+HER2-, CK-HER2+, and CK+HER2+. The inclusion of anti-HER2 increased the percentage of informative samples with ≥ 1 detectable CTC from 89 to 95%. CK- CTCs and tdEVs correlated equally well with the clinical outcome as CK+ CTCs and tdEVs. Inter- and intra-patient heterogeneity was found for the CTC/tdEV phenotypes, and the presence of 2 or 3 classes of CTCs/tdEVs was associated with worse prognosis compared to a uniform CTC/tdEV phenotype present (1 class). The use of ≥ 7% HER2+CK+ tdEVs can predict HER2 expression of the tissue with 74% sensitivity and specificity using the HER2 amplification status of the primary tumor as a classification variable.
CONCLUSIONS: HER2 can be detected on CTCs and tdEVs not expressing CK, and these CK- CTCs/tdEVs have similar clinical relevance to CTCs and tdEVs expressing CK. tdEVs perform better than CTCs in predicting the HER2 status of the primary tissue. CTC and tdEV heterogeneity in the blood of patients is inversely associated with overall survival.

Background

The availability of more and more targeted treatments necessitates the development of techniques to screen cancer patients for the presence of the respective therapeutic targets in the primary and metastatic lesions. However, the invasiveness of sampling of biopsies from all different metastatic lesions and the associated patient morbidity demands the need for a less invasive test. EpCAM+ circulating tumor cells (CTCs) isolated from the peripheral blood of cancer patients constitute clinically relevant and easily accessible tumor material from a single tube of blood [1, 2]. The possibility of longitudinal measurements and the screening of CTCs for therapeutic targets can provide clinicians with the mutational status of the tumor and the presence of treatment targets in real-time facilitating their decisions on treatment monitoring [3, 4]. The reliability of such an assessment increases with the number of CTCs available. In metastatic breast cancer, ~ 50% of patients have ≥ 5 CTCs detected with the CellSearch system and the percentage of patients with ≥ 10 and ≥ 100 CTCs decreases rapidly [1].

The development of the open-source ACCEPT software allowed the automated enumeration of all objects in the fluorescence images and the more objective assessment of treatment targets on CTCs [5] eliminating inter- and intra-operator variations [6].

Large tumor-derived extracellular vesicles (tdEVs) are co-isolated with CTCs and are present in an order of magnitude higher frequencies than CTCs [7, 8]; thus, they can increase the number of “readable” patient samples for the assessment of therapeutic targets.

The original study used in the present analysis is the CLCC-IC-2006-04 (NCT00898014) with 267 stage IV breast cancer patients receiving first-line chemotherapy [9]. The primary objective was to predict overall survival (OS) and progression-free survival (PFS) by counting all CK+ CTCs before the initiation of the second course of chemotherapy [9]. The secondary objective was to evaluate the HER2 status on CTCs over the course of a treatment. Ligthart et al. [3] and Zeune et al. [5] quantified the HER2 expression on the pre-scored by the operator CTCs in an automated manner demonstrating a more objective and less biased assessment among the different operators.

In the present study, we wanted to address two important questions:

Is HER2 expressed on tdEVs?

Is HER2 expressed on CTCs lacking CK expression? And if so, do these populations associate with the clinical outcome of the patients?

To address the aforementioned research questions, the open-source ACCEPT software was used and the enumeration of all CTCs and tdEVs falling in the different immunophenotypic classes was performed in an automated manner.

Methods

CellSearch images of breast cancer patients

Digitally stored CellSearch® (Menarini Silicon Biosystems, Huntingdon Valley, PA, USA) image files from 98 metastatic breast cancer patients from a previously reported study (CLCC-IC-2006-04 study, ClinicalTrials.gov Identifier: NCT00898014) were re-analyzed [3, 9]. The original study included 267 patients with either HER2+ or HER2− primary tumor, whose characteristics at the baseline have been previously described [9]. Of these samples, 114 were further screened for HER2 expression on CTCs. These samples had been processed with the CellSearch system using the CTC kit, as previously described [7] including an additional fluorescein (FITC)-conjugated monoclonal antibody recognizing HER2 (clone Her81) in the staining mixture [11]. The specific clone recognizes an epitope in the extracellular domain of HER2, and it does not cross block binding of the FDA-cleared trastuzumab used in the clinics [11]. Ninety-eight out of 114 image datasets (where the anti-HER2 had been included in the staining mix) were available for the present analysis.

Classification and automated enumeration of CTCs and tdEVs using ACCEPT

For the automated CTC and tdEV count enumeration, the digitally stored fluorescence image files were processed with the Automated CTC Classification, Enumeration and PhenoTyping (ACCEPT) software v1.1 (http://github.com/LeonieZ/ACCEPT) using the “Full Detection” function. The ACCEPT software detects all objects, present in the fluorescence images with a size larger than 4 pixels, and it extracts for each of them measurements of 10 morphological and fluorescence intensity features per fluorescence channel [12]. The operator can define the classes of their interest by designing linear gates using the desired features. Once optimized, the gates can be applied in all samples, and the counts of objects falling within them can be automatically extracted. The application of the same gates for all samples allows fast enumeration and elimination of inter- and intra-operator variability and bias leading to a more objective consensus [6].

The herein used gates took into account the roundness (eccentricity and perimeter to the area), min and max sizes (area and perimeter), mean fluorescence intensity of all different channels, and the overlap of CK/HER2 with DNA (in case of CTCs). The exact gates that were applied can be found in Supplementary Table S1, Additional file 3.

Statistical analysis

Statistical analysis was performed in IBM SPSS Statistics v24.0 (SPSS Inc., Chicago, IL, USA) and MedCalc v19 (MedCalc Software, Ostend, Belgium). For the matched CTC and tdEV counts of the same subclass, a two-tailed Spearman’s rho test was performed to evaluate their relation through a monotonic function and the non-parametric Wilcoxon signed ranks test was used to test the equality of their distributions. The non-parametric Mann-Whitney U test was used to compare the distributions for the absolute and relative frequencies of CTCs and tdEVs with different immunophenotypes. The open-source web application Cutoff Finder (http://molpath.charite.de/cutoff/) was used to calculate the hazard ratios (HRs) for overall survival (OS) including 95% confidence intervals (CIs), over a wide range of cutoff values for CK+/− CTCs and tdEVs. Cutoff Finder uses the R code to provide optimization and visualization tools for cutoff determination [13]. Receiver operating characteristic (ROC) curves were used to evaluate the performance of %HER2+ CTCs (or tdEVs) in predicting HER2 amplification status of the primary tumor as assessed by FISH (available data for 92/98 patients). The areas under the curve (AUCs) were compared using the non-parametric DeLong approach [14]. We determined a minimum threshold for %HER2+CK+ CTCs or tdEVs to predict a HER2+ tissue as the value that led to equal sensitivity and specificity (value that led to minimum |sensitivity − specificity|). OS was defined as the elapsed time between blood draw and death, and it was available for 94/98 patients. The patients who were still alive at the last follow-up were censored. Kaplan-Meier (KM) survival curves of OS were used to compare (using the log-rank test) patients with favorable and unfavorable manual and automated CTC and tdEV counts. Cox proportional hazards regression analysis was used to determine the univariable HRs with 95% CIs of the categorical and continuous CTCs and tdEVs for OS. Continuous CTC and tdEV counts were log-transformed to achieve a better model fit.

Results

Patient subpopulation included in the present analysis

The patients included in the present study were not chosen based on their CTC counts; that would lead to a selection bias. There was no significant difference in regard to the probability distribution for OS (p = 0.863, log-rank test) and the distributions for the manual CTC counts (p = 0.600, Mann-Whitney U test) between the complete patient cohort and the patient subpopulation included in the present analysis. The Kaplan-Meier of OS for the complete patient cohort of the study and the patient subpopulation chosen as well as the box plots of the manual CTC counts of the corresponding patients are shown in panels A and B of Additional file 1, respectively.

ACCEPT display and automated enumeration of CTCs and tdEVs

After the inclusion of the HER2-FITC antibody, 3 different subclasses of CTCs and tdEVs were detected in the CellSearch images (Additional file 2). Linear ACCEPT gates were designed (Supplementary Table S1, Additional file 3) and applied in the image datasets for the automated enumeration of CTCs/tdEVs of each subclass.

The automated CTC and tdEV counts of the 98 breast cancer patients are shown in dot plots in Fig. 1. Comparison of the applied ACCEPT gates for CTC automated enumeration and the manual CK+ CTC counts scored by the operator in the original study showed a high correlation (rS = 0.83, p < 0.01) and similar association to the clinical outcome (Additional file 4).

Fig. 1

Frequencies of CTC and tdEV subpopulations in metastatic breast cancer patients. Dot plots of different subclasses (indicated in the x-axis) of CTCs (in gray) and tdEVs (in black) in 98 metastatic breast cancer patients. Each class/immunophenotype is indicated by a different color on the x-axis. Each dot corresponds to the counts of one patient. Horizontal black lines indicate the median values

Ninety-three of 98 (95%) patients had ≥ 1 detectable CTCs, and 98/98 (100%) of patients had ≥ 1 tdEV (Fig. 1). Whereas CK+HER− tdEVs are found in a ~ 15-fold higher frequency compared to CK+HER− CTCs (p < 0.01, Wilcoxon signed ranks test), the distributions of CK+HER2+ and CK−HER2+ tdEVs are not significantly different compared to the respective distributions of CTCs (Wilcoxon signed ranks test). Furthermore, the correlation between CTC and tdEV subclasses was very strong only in case of CK+HER2− populations (rS = 0.84, p < 0.01), whereas in case of CK−HER2+ and CK+HER2+ populations, a moderate correlation was found (0.5 < rS < 0.6, p < 0.01).

Association of CK+ and CK− CTCs and tdEVs with clinical outcome

The original CellSearch CTC definition requires the expression of CK within the nucleated cells of a size larger than 4 μm, after EpCAM enrichment. Similarly, our previously reported tdEV definition enclosed all events isolated with the CellSearch system that were positive for CK and negative for CD45, without a nucleus and with a size range between 1 and 14 μm [8]. The CK- and CD45-negative cell and EV populations that are expressing a treatment target, such as HER2, could indicate tumorigenic origin. To determine whether these lacking CK and expressing HER2 objects are actually cancer-associated, we named them CK− CTCs and tdEVs and compared their prognostic power with CTCs and tdEVs expressing CK. For that comparison, KM plots for OS were generated (Fig. 2) stratifying patients based on their CK+ CTCs (panel a), CK+ tdEVs (panel b), CK− CTCs (panel c), or CK− tdEVs (panel d). The selected cutoffs were 5 for CK+ CTCs (established cutoff value), 1 for CK− CTCs (since we wanted to evaluate whereas the presence itself is associated with worse prognosis), 20 for CK+ tdEVs (as the upper bound of the normal range that we reported previously [8]), and 10 for CK− tdEVs (assuming that each class contributes 1/3 in the total tdEV counts). However, the screening of all possible cutoff values to dichotomize patients into two risk groups and their associated HR of OS are shown in the respective plots of HRs over the CK+/− CTC and tdEV distributions (Additional file 5).

Fig. 2

Association of CK+ and CK− CTCs and tdEVs with the clinical outcome of metastatic breast cancer patients. KM curves for OS dichotomizing patients based on their automated CK+ CTCs (a) and tdEVs (b) and their CK− CTCs (c) and tdEVs (d)

CK− CTCs (or tdEVs) perform similarly to CK+ CTCs (or tdEVs) in predicting the event of the patients as shown by the patient stratification and the resulting HR with patients with increased CK−/CK+ CTCs/tdEVs having the tendency to have worse clinical outcome compared to the patients with lower counts in their blood.

The univariable Cox regression of the continuous log-transformed CK+ and CK− CTC and tdEV counts (Supplementary Table S2, Additional file 3) showed that only CK+ tdEVs are significant predictors of OS.

Heterogeneity of CTC/tdEV immunophenotypes and association with clinical outcome

Inter- and intra-patient heterogeneity was observed regarding the CK and HER2 phenotype of CTCs and tdEVs. The Venn diagrams in Fig. 3a show the number of patients with one or more subclasses (CK+HER2−, CK+HER2+, and CK−HER2+) of CTCs or tdEVs detected. Interestingly, the majority of patients had all three subclasses of CTCs present, whereas very few patients had solely one class of CTCs. More specifically, 56/93 (60%) of patients had CTCs of all three immunophenotypes, followed by 22/93 (24%) of patients with CTCs of 2 different immunophenotypes and only 15/93 (16%) of patients with only 1 class of CTCs present. In case of tdEVs, a threshold of 10 was considered to be the “noise,” and only when a class was present in a frequency above that number, it was considered to be present. That led to 18/98 (18%) of patients having all 3 subclasses of tdEVs present, 22/98 (22%) of them having 2 subclasses present, and 38/98 (39%) of them having one class of tdEVs present.

Fig. 3

Inter- and intra-patient CTC and tdEV heterogeneity (a) and association with clinical outcome (b). The numbers of patients with the respective CTC (left) and tdEV (right) classes present in their blood are shown as numbers in the corresponding parts of the Venn diagrams (a). For a tdEV class to be considered present, a threshold of 10 tdEV counts was used. Patients with uniform CTC and tdEV immunophenotypes (1 class) have a better prognosis compared to patients with more heterogeneous CTC and tdEV immunophenotypes (2 or 3 classes) present (b)

Patients were grouped based on the CTC/tdEV classes present in their blood, and the respective Kaplan-Meier curves of overall survival were plotted to demonstrate whereas the presence of the heterogeneous immunophenotypes is related to worse prognosis (Fig. 3b). The presence of heterogeneous immunophenotypes of CTCs/tdEVs in the blood of metastatic cancer patients (2 or 3 classes present) was associated with a worse clinical outcome compared to the patient population with a uniform CTC/tdEV immunophenotype (1 class present).

HER2 amplification of tissue versus HER2 expression on CTCs and tdEVs

For 92 out of 98 patients, the HER2 status of the tissue was assessed by fluorescence in situ hybridization (FISH) with 39 patients (42%) having HER2+ and 53 (58%) HER2− tissue.

When we compare these 2 groups of patients in terms of their automated CTC and tdEV counts of each subclass, only CK+HER+ tdEV distribution is significantly increased (p < 0.05, Mann-Whitney U test) in HER2+ patients (Panel B, Additional file 6). The rest of the CTC and tdEV distributions (Panels A and B, Additional file 6) are not significantly different. However, the differences become more profound when comparing the 2 groups in regard to %HER2+CK+/−, %HER2+CK+, and %HER2+CK− CTCs and tdEVs estimated over the total CTC and tdEV counts, respectively. Patients with HER2+ tissue have significantly lower percentage of CK+HER2− CTCs and tdEVs and significantly higher percentage of CK+HER2+ CTCs and tdEVs compared to patients with HER2− tissues (Panels C and D, Additional file 6). No statistically significant difference could be found concerning the percentage of CK−HER2+ CTCs and tdEVs.

To evaluate which HER2+ population is in a better concordance with the HER2+ primary tumor, ROC curves of the different HER2+ CTC (Fig. 4a) and tdEV (Fig. 4b) proportions were constructed treating HER2+ tissue as the classification variable. %HER2+CK+ CTCs and tdEVs performed better as indicated by the larger AUCs (0.69 for CTCs and 0.79 for tdEVs, which were not significantly different, DeLong test). The asterisks in panels a and b indicate the selected threshold of %HER2+CK+ CTCs and tdEVs, for which the test has sensitivity ≈ specificity. In case of CTCs, more than 23% HER2+CK+ CTCs (Fig. 4c) could predict a HER2+ tissue with a sensitivity of 65% and specificity of 66%. Likewise for tdEVs, more than 7% HER2+CK+ tdEVs (Fig. 4d) could predict HER2+ tissue with a sensitivity and specificity of 74%. The influence of the increasing number of CTCs and tdEVs on the accuracy of the respective test is summarized in Supplementary Table S3, Additional file 3. The accuracy of ≥ 23% HER2+CK+ CTC test improved and reached 80% when ≥ 50 CTCs were detected; however, only 27/98 (28%) patients accounted for that CTC load. In case of %HER2+CK+ tdEVs, the accuracy of the test was constantly above 70% and gradually improved from 74 up to 91% with increasing tdEV counts.

Fig. 4

Prediction of HER2 status of tissue from CTCs and tdEVs. ROC curves of %HER2+CK− (light gray lines), HER2+CK+ (dark gray lines), and total HER2+CK+/− (black lines) CTCs (a) and tdEVs (b) treating HER2+ tissue (as assessed by FISH) as the classification variable. %HER2+CK+ populations performed the best as shown by the largest AUCs. The asterisks indicate the points, where sensitivity ≈ specificity for CTCs and tdEVs (23% HER2+CK+ CTCs leading to 65% sensitivity and 66% specificity; 7% HER2+CK+ tdEVs leading to 74% sensitivity and specificity). Scatter plot of total CTCs (c) and tdEVs (d) for each patient (x-axis). Samples were sorted on the percentage of HER2+CK+ CTCs or tdEVs respectively indicated by the dark gray lines. On the top of the panels and along the x-axis, the HER2 status of the tissue is indicated as positive (+) or negative (−), as evaluated by FISH. The vertical black dashed lines indicate the 23% HER2+ CTCs (and 7% HER2+ tdEVs) threshold, right of which the tissue of the patient could be considered as HER2+

Discussion

The HER2 is found overexpressed in around 15–30% of breast cancer cases mainly because of the HER2/neu oncogene amplification [15, 16]. Randomized clinical trials have demonstrated its predictive and prognostic value with HER2+ breast cancer patients having a remarkably improved clinical outcome when treated with anti-HER2 therapy, such as trastuzumab or lapatinib next to chemotherapy, as compared to HER2+ patients treated with solely chemotherapy [17-20]. Furthermore, HER2− breast cancer patients randomized in two treatment arms with and without anti-HER2 therapy did not show any clinical benefits further supporting the predictive value of HER2 [21]. These observations necessitate the use of an accurate test to assess the HER2 tumor status and facilitate the clinician’s treatment decision-making. The current state of the art evaluates HER2 status on solid biopsies (either tumor needle biopsy or whole tumor after resection) by immunohistochemistry or/and FISH. Several guidelines have improved the accuracy of HER2 evaluation [22]. Nevertheless, there are cases of discorcondant immunohistochemistry and FISH results or high HER2 heterogeneity of the tumor [23, 24] preventing an objective consensus on HER2 status among different operators.

Non-invasive liquid biopsies have emerged to be promising alternatives to solid biopsies providing clinically relevant tumor material in real-time. The increasing load of CK+ CTCs and tdEVs as identified by the CellSearch system is strongly associated with worsening progression-free and overall survival [1, 2, 8, 10, 25]. The CellSearch system immunomagnetically enriches CTCs targeting the epithelial cell adhesion molecule (EpCAM); consequently, all EpCAM-negative or low expressing CTCs are missed by the system. Many approaches have been introduced to overcome this limitation by enrichment of CTCs based on their physical properties such as size, density, and charge [26], but no obvious advantages have been demonstrated in clinical studies so far. Importantly, EpCAMlow, CK+ CTCs identified after size-based separation of EpCAM-depleted blood samples did not correlate with clinical outcomes in metastatic non-small cell lung, prostate, and breast cancer [27, 28]. In the present study, we investigated whether the expression of HER2 can be assessed through tdEVs and whether CTCs and tdEVs are missed by the system due to the lack of CK expression. The whole workflow from EpCAM enrichment to CTC and tdEV scoring can be done in a fully automated manner, once a blood sample of 7.5 mL is available. Consequently, the test is not subjected to the quality of the sample, the handling, the staining methodology, and the judgment of the technician. It can be performed in a timely fashion and can facilitate the treatment monitoring of the patient.

The present analysis validates the association of CK+ tdEVs with the overall survival of metastatic breast cancer patients since it is based on an independent patient cohort than the one that we previously reported (IMMC01 study) [8]. HER2 expression could be detected on CK+ tdEVs that were co-isolated with CTCs using the CellSearch system.

The ACCEPT software allows the visualization and automated enumeration of CK− objects, which are not presented to the operator by the FDA-cleared CellSearch image analysis algorithm (CellTracks Analyzer II). Interestingly, the inclusion of the HER2-FITC antibody allowed the detection of CD45−, CK−, and HER2+ CTCs and tdEVs in the EpCAM-enriched blood samples of metastatic breast cancer patients increasing the percentage of patient samples with detectable CTCs from 89 to 95%. CK+ tdEVs were already present in 99% of patient samples, which increased to 100% after the inclusion of anti-HER2.

Although we have not genetically proven that these CK− CTCs and tdEVs are indeed cancerous, their similar correlation with OS (Fig. 2) tends to confirm that hypothesis. The aneuploidy of these CK− cells should be addressed by FISH in a prospective clinical study, since the cartridges, whose image datasets were used in the present analysis, are no longer available. Fehm et al. [29] first demonstrated the malignancy of EpCAM-enriched CK+, CD45−, and DAPI+ cells (CellSearch CTC definition) by FISH from the blood of patients with carcinomas. A protocol to perform FISH directly in the CellSearch cartridges was developed later by Swennenhuis et al. [30], who evaluated the technique in blood samples of castration-resistant prostate cancer patients.

The observation of CK−HER2+ CTCs and tdEVs raises questions about other CK−HER2− CTC and tdEV populations present in the EpCAM-enriched samples. The findings of Crespo et al. that aneuploid CD45−, CK−, and AR+ CTCs in EpCAM-enriched blood samples of castration-resistant prostate cancer patients are associated with worse OS further support our hypothesis [31]. The detection of CK− aneuploid cells in blood samples of ovarian, breast, and colorectal cancer patients has been also described by Pecot et al. [32]. All aforementioned studies further encourage the inclusion of additional antibodies to detect nucleated events of unknown cell lineage in the EpCAM-enriched samples and decrease false-negative rates of CTCs.

Importantly, all three immunophenotypes of CTCs (HER2+CK−, HER2−CK+, and HER2+CK+) were present in the majority of EpCAM-enriched patient samples, and the presence of heterogeneous CTC and tdEV immunophenotypes was associated with worse clinical outcome of the patients (Fig. 3).

Our next question was whereas a minimum threshold of a HER2+ population could predict the HER2 status of the tissue as a real-time liquid biopsy. The %HER2+CK+ tdEVs performed the best achieving higher sensitivity and specificity than CTCs (Fig. 4), most likely because of the higher frequencies of tdEVs better reflecting the tumor heterogeneity. Based on that test, a patient can be affirmed to have a HER2+ tumor with at least 74% accuracy if ≥ 7% of total EpCAM+ tdEVs are CK+ and HER2+. However, it is important to mention that the determination of the HER2 amplification status of the primary tissue was done at the time of diagnosis with most of the patients being as stages I–III. The CTC evaluation reported in the original study (and the present analysis) occurred after the progression of the patients to stage IV (one of the inclusion criteria of the patients enrolled in the study), which can be years after the initial diagnosis and the tumor resection. It is well known that the phenotype and genotype of the primary tumor can substantially differ from the respective genotypic and phenotypic features of the metastatic lesions as the disease progresses [33-36] with CTCs resembling better the mutational status of metastatic lesions than of primary tumor [37]. Therefore, the HER2 expression on the primary tumor can differ substantially from the status of the metastatic lesions as well the CTCs and tdEVs isolated having an impact on the concordance found. Whether the HER2 assessment via a liquid biopsy (CTCs or/and tdEVs) can better predict response to anti-HER2-targeted therapies compared to the current assessment via a solid biopsy remains to be addressed.

Molecular characterization of the EpCAM+ tdEVs of patients undergoing HER2-targeted therapies at follow-up time points can contribute to better comprehend the underlying mechanism of tumor resistance to anti-HER2 treatment [38], which is observed in more than 70% of HER+ breast cancer patients within a year from the initiation of the treatment [39]. Ciravolo et al. have already suggested a mechanism of anti-HER2 resistance by the increased binding efficiency of HER2+ exosomes to trastuzumab in progressive HER2+ breast cancers as compared to earlier stages of breast cancer [40]. Another mechanism has been suggested by Al-Nedawi et al. with tdEVs transferring the oncogenic form of epidermal growth factor receptor EGFRvIII to cells without that immunophenotype [41].

Conclusions

The inclusion of anti-HER2 in the CellSearch assay allowed the detection of EpCAM+CK− CTCs and tdEVs with similar prognostic power as EpCAM+CK+ CTCs and tdEVs in metastatic breast cancer emerging the importance of including more detection antibodies next to CK for the identification of CK− CTCs and tdEVs. Importantly, the presence of heterogeneous CTC and tdEVs (2 or 3 different immunophenotypes) in the blood of cancer patients was associated with worse overall survival compared to patients with uniform CTCs and tdEVs (1 immunophenotype). These findings should be validated in a prospective study, where also the malignancy of these populations should be addressed by FISH. Our results enrich the already available data of real-time liquid biopsy and encourage the screening of tdEVs for treatment targets in clinics, since they seem to better reflect the HER2 phenotype of the primary tumor than CTCs opening the path towards a more rational and objective choice of the patients who will or will not be subjected to HER2 targeting therapies.

Additional file 1: Figure S1. Comparison between the patient subpopulation used in the present analysis and the complete patient cohort of the CLCC-IC-2006-04 study. Probability distribution for Overall Survival of complete patient cohort (in black) and patient subpopulation used (in pink) as shown by the respective Kaplan Meier plots (Panel A). Box plot with overlapping data of manual CTC counts of complete patient cohort (in black) and patient subpopulation used (in pink).

Additional file 2: Figure S2. ACCEPT display of different subclasses of CTCs (left) and tdEVs (right) isolated by the CellSearch system from metastatic breast cancer patients. Examples of CK+HER2- (Panel A), CK-HER2+ (Panel B) and CK+HER2+ (Panel C) CTCs and tdEVs. The red contours around the objects in the respective channel indicate the detected and segmented signal by the ACCEPT image analysis algorithm. All objects were isolated from 7.5 mL of blood using the CellSearch system. Scale bars indicate 6.4 μm.

Additional file 3: Supplementary Table S1. Linear ACCEPT gates used for the automated enumeration of CTCs and tdEVs. Supplementary Table S2. Univariable Cox regression analyses of CK+ and CK- CTCs and tdEVs after log transformation. Supplementary Table S3. Sensitivity, specificity and accuracy of ≥ 23% CTCs and ≥ 7% tdEVs, double positive for CK and HER2, as tests to predict the HER2 status of the tissue. The accuracy increases with the total CTCs and tdEVs detected (≥ 1, 5, 10, 20, 50, 100) at the cost of number of eligible patients to be assessed.

Additional file 4: Figure S3. Correlation of manual with automated CK+ CTC counts and association with clinical outcome of patients. Scatter plot of CK+ manual CTCs (mCTCs) versus CK+ automated CTCs (aCTCs) showing strong correlation (Panel A). KM plots of OS (Panel B) for patients with < and ≥ 5 CTCs. The dichotomization of patients was done based on either manual (black and grey lines) or automated (red and green) CTC counts showing equivalent association to OS.

Additional file 5: Figure S4. Overview plots of HRs (with 95% CI) for all possible cut-off values for CK+ CTCs (Panel A), CK+ tdEVs (Panel B), CK- CTCs (Panel C) and CK- tdEVs (Panel D). The rug plots at the bottom of Panels A-D correspond to the value distributions of CK+ CTCs, CK+ tdEVs, CK- CTCs and CK- tdEVs respectively. For CK+ CTCs, a larger percentage of cut-off values (31%, Panel A) could significantly dichotomize patients with a higher and lower risk as compared to CK- CTCs (13%, Panel C). The opposite was observed for tdEVs with a larger percentage of cut-off values for CK- tdEVs (30%, Panel D) leading to a significant dichotomization of patients as compared to CK+ tdEVs (14%, Panel B).

Additional file 6: Figure S5. Comparison of patients with HER2+ and HER2- tissues in regards to their relative and absolute frequencies of CTCs and tdEVs of different phenotypes. Box plots with data overlap depicting the automated CTC counts (Panel A), automated tdEV counts (Panel B), % of CTCs (Panel C) and % of tdEVs (Panel D) of the 3 different immunophenotypes (indicated in the x-axis) in metastatic breast cancer patients split based on the HER2 status of their tissue. Each dot corresponds to the counts of one patient (green dots: patients with HER2+ tissue (N=39), red dots: patients with HER2- tissue (N= 53), there was no available HER2 status for N=6 patients). Lower and upper bounds of box plots correspond to the 1st (Q1) and 3rd (Q3) quartile of data, horizontal black lines indicate median vaues. Whiskers indicate 1.5*(Q3-Q1). * indicates significant (p < 0.05) and ** highly significant (p < 0.001) statistical difference (Mann-Whitney U test).