Skip to main content
Download PDF
  • Research article
  • Open access
  • Published:

The prognostic value of JUNB-positive CTCs in metastatic breast cancer: from bioinformatics to phenotypic characterization

Breast Cancer Research volume 21, Article number: 86 (2019) Cite this article

Abstract

Background

Circulating tumor cells (CTCs) are important for metastatic dissemination of cancer. They can provide useful information, regarding biological features and tumor heterogeneity; however, their detection and characterization are difficult due to their limited number in the bloodstream and their mesenchymal characteristics. Therefore, new biomarkers are needed to address these questions.

Methods

Bioinformatics functional enrichment analysis revealed a subgroup of 24 genes, potentially overexpressed in CTCs. Among these genes, the chemokine receptor CXCR4 plays a central role. After prioritization according to the CXCR4 corresponding pathways, five molecules (JUNB, YWHAB, TYROBP, NFYA, and PRDX1) were selected for further analysis in biological samples. The SKBR3, MDA-MB231, and MCF7 cell lines, as well as PBMCs from normal (n = 10) blood donors, were used as controls to define the expression pattern of all the examined molecules. Consequently, 100 previously untreated metastatic breast cancer (mBC) patients (n = 100) were analyzed using the following combinations of antibodies: CK (cytokeratin)/CXCR4/JUNB, CK/NFYA/ΥWHΑΒ (14-3-3), and CK/TYROBP/PRDX1. A threshold value for every molecule was considered the mean expression in normal PBMCs.

Results

Quantification of CXCR4 revealed overexpression of the receptor in SKBR3 and in CTCs, following the subsequent scale (SKBR3>CTCs>Hela>MCF7>MDA-MB231). JUNB was also overexpressed in CTCs (SKBR3>CTCs>MCF7>MDA-MB231>Hela). According to the defined threshold for each molecule, CXCR4-positive CTCs were identified in 90% of the patients with detectable tumor cells in their blood. In addition, 65%, 75%, 14.3%, and 12.5% of the patients harbored JUNB-, TYROBP-, NFYA-, and PRDX-positive CTCs, respectively. Conversely, none of the patients revealed YWHAB-positive CTCs. Interestingly, JUNB expression in CTCs was phenotypically and statistically enhanced compared to patients’ blood cells (p = 0.002) providing a possible new biomarker for CTCs. Furthermore, the detection of JUNB-positive CTCs in patients was associated with poorer PFS (p = 0.015) and OS (p = 0.002). Moreover, JUNB staining of 11 primary and 4 metastatic tumors from the same cohort of patients revealed a dramatic increase of JUNB expression in metastasis.

Conclusions

CXCR4, JUNB, and TYROBP were overexpressed in CTCs, but only the expression of JUNB was associated with poor prognosis, providing a new biomarker and a potential therapeutic target for the elimination of CTCs.

Background

The malignant nature of circulating tumor cells (CTCs) has been recently demonstrated using xenograft models [1, 2], and it is widely accepted that these cells are responsible for tumor metastasis. CTCs have been detected in several tumor types, including breast, prostate, colorectal, non-small cell lung cancer (NSCLC), and small cell lung cancer (SCLC) [3,4,5,6,7,8]. Their detection has been associated with poor patients’ clinical outcome, irrespectively of the disease stage. In addition, the presence of CTCs before and after the completion of adjuvant chemotherapy in breast cancer patients is associated with poor clinical outcome [6, 9, 10].

Several studies have shown that CTCs present phenotypic and molecular differences from the corresponding primary tumor cells [11,12,13] and their changes in response to treatment could provide a useful alternate to solid biopsy [13, 14]. In addition, their detection and enumeration can give useful information regarding patient’s prognosis providing a surrogate marker for treatment efficacy [15,16,17]. Nevertheless, it is well known that chemotherapy cannot effectively eliminate all CTCs; therefore, the identification of new molecules on these cells could offer new therapeutic targets.

Furthermore, only a subgroup of patients with detectable CTCs will eventually relapse. This is attributed to the fact that some of the identified CTCs are destined to die [18]. Therefore, only a proportion of the detected tumor cells have metastatic potential; however, it is currently unknown which is the dangerous phenotype.

CTCs hold epithelial and/or mesenchymal (EMT) characteristics as well as stem cell features [19,20,21], and the acquisition of an EMT phenotype makes them invisible with common detecting platforms, based on the expression of epithelial markers. Consequently, there is an urgent need for new biomarkers to improve the identification and characterization of these cells.

Bioinformatics analysis is a new tool, based on Gene Expression databases, which can provide information about a number of genes-proteins that could potentially be upregulated in CTCs. Thus, the aim of the current study was the identification of novel biomarker/therapeutic targets in CTCs isolated from metastatic breast cancer patients based on bioinformatics analysis.

According to this analysis, CXCR4 is the most important molecule, potentially overexpressed in CTCs, while JUNB, YWHAB, TYROBP, NFYA, and PRDX1 were also found to be related to CXCR4 in different pathways and they could be potentially upregulated in CTCs [22].

Table 1 Patient characteristics
Full size table

It is important that CXCR4 has been found to be a prognostic marker in various tumor types, including BC [23]. It is upregulated in tumor tissue as compared to normal tissues. It has also been reported that CXCR4 is important for CTC-seeding and metastatic potential in early-stage BC of node-positive patients. Its expression was related to epithelial-mesenchymal transition (EMT), hence, associated with an acutely malignant phenotype [24].

JUNB has been also associated to invasion/metastasis in solid tumors including BC [25,26,27] and represented as an important target in diseases associated with EMT, including cancer and fibrosis [28, 29]. It has been implicated in the earliest events of the development of resistance to kinase inhibitors in BC [30]; however, the expression and the prognostic values of this transcription factor in CTCs have not been studied. Furthermore, the rest of the examined molecules (YWHAB, TYROBP, NFYA, and PRDX1) have been related to cancer progression and invasion [31,32,33,34,35], but to the best of our knowledge, they have not been studied in CTCs.

Methods

Dataset integration and bioinformatics analysis

We collated an ensemble of nine human genome microarrays (GSE22820; GSE19783; GSE31364; GSE9574; GSE18672; GSE27562; GSE16443; GSE15852; GSE12763), from the Gene Expression Omnibus (GEO) database [36]. The integration of the datasets followed the process described in the literature [22], resulting in 498 primary breast cancer tissues and 124 blood samples from breast cancer patients, as well as 104 normal breast tissue samples and 85 peripheral blood samples from healthy donors. Three comparative analyses were performed in the pooled set of samples, using the “Significance Analysis of Microarrays” (SAM) [37] as the feature selection approach of choice with the siggenes package of R/Bioconductor using the same parameters for all genes entries, in order to reveal comparison-specific differentially expressed genes. The false discovery rate (FDR) [38] was used as the criterion for determining the set of genes that exhibit differential expression with critical value set to 0.01 for all comparisons.

Patient samples and cytospin preparation

One hundred newly diagnosed and treatment-naïve patients with metastatic breast cancer (mBC) were enrolled in the study. All patients were enrolled in one center (Department of Medical Oncology, University General Hospital of Heraklion, Crete, Greece) and were treated with front-line chemotherapy, according to the national guidelines regardless of the CTCs’ results. Patients’ characteristics are shown in Table 1. In addition, 10 female normal blood donors were also included as controls. All blood samples were obtained at the middle of vein puncture, after the first 5 ml of blood was discarded. These precautions were undertaken in order to avoid contamination of the blood sample with epithelial cells from the skin during sample collection. All patients gave their informed consent to participate in the study, which has been approved by the Ethics and Scientific Committees of our Institution.

Table 2 Mean intensity of JUNB, CXCR4, TYROBP, PRDX1, NFYA, and YWHAB in breast cancer cell lines, PBMCs, and CTCs
Full size table

Twenty milliliters of blood in EDTA was obtained from all patients, and peripheral blood mononuclear cells (PBMCs) were isolated with Ficoll-Hypaque density gradient (d = 1077 g/mol) centrifugation at 1800 rpm for 30 min. PBMCs were washed three times with PBS and centrifuged at 1500 rpm for 10 min. Aliquots of 500,000 cells were cyto-centrifuged at 2000 rpm for 2 min on glass slides. Cytospins were dried up and stored at – 80 °C. Two slides with 106 cells in total were analyzed from each patient for every molecule and for CK/CD45 staining. This number of PBMCs according to previous studies of our group could be representative of the characterization of CTCs in patients’ blood [18, 19, 39].

Cell cultures

All cell lines were obtained from ATCC (American Type Culture Collection, USA) and used for spiking experiments. The MCF7 adenocarcinoma cells were cultured in 1:1 Dulbecco’s modified Eagle medium (DMEM Glutamax) (GIBCO-BRL Co, MD USA) supplemented with 10% fetal bovine serum (FBS) (GIBCO-BRL), 16 ng/ml insulin, and 50 mg/ml penicillin/streptomycin (GIBCO-BRL). MDA-MB-231 cells were cultured in DMEM, supplemented with 10% FBS and 50 mg/ml penicillin/streptomycin. SKBR3 breast cancer cells were cultured in RPMI supplemented with 10% FBS. Hela were cultured in 1:1 Dulbecco’s modified Eagle medium (DMEM Glutamax) supplemented with 10% FBS and 50 mg/ml penicillin/streptomycin. Cells were maintained in a humidified atmosphere of 5% CO2 in the air. Sub-cultivation was performed with 0.25% trypsin and 5 mM EDTA (GIBCO-BRL). All experiments were performed during the logarithmic growth phase.

Double staining experiments and confocal laser scanning microscopy

The presence of CK-positive cells in PBMCs’ cytospin was investigated using the mouse A45-B/B3 (detecting CK8, CK18, and CK19) (Micromet Munich, Germany) antibody and anti-CD45 (common leukocyte antigen) (Santa Cruz, Santa Cruz, CA, USA) in order to exclude possible ectopic expression of cytokeratin on hematopoietic cells. The cyto-morphological criteria proposed by Meng et al. [40] (i.e., high nuclear/cytoplasmic ratio, larger cells than white blood cells, etc.) were used in order to characterize a CK-positive cell as a CTC.

Fixed cells with acetone/methanol (9:1 v/v) were stained with A45-B/B3 antibody for 1 h, followed by Alexa488 anti-mouse for 45 min. Afterwards, CD45 anti-rabbit was used for 1 h. Consequently, the samples were stained with the corresponding Alexa555 anti-rabbit fluorochrome. Finally, cells were stained with DAPI (Invitrogen, Carlsbad, CA, USA) conjugated with antifade and analyzed with confocal laser scanning microscopy.

Triple immunofluorescence

Triple immunofluorescence was performed to cytospins’ preparations from all patients harboring CK (+)/CD45 (−) CTCs. For the CK/CXCR4/JUNB staining, cells were fixed with 3% paraformaldehyde (PFA) for 30 min at room temperature (RT). Permeabilization was achieved with 0.5% Triton X-100 for 10 min at RT. After blocking with PBS supplemented with 10% (v/v) FBS for 1 h, cells were incubated with A45-B/B3 for 1 h and then with Alexa 488 for 45 min. Consequently, slides were incubated with JUNB antibody conjugated with Alexa647 (Santa Cruz, CA, USA) for 1 h. CXCR4 anti-rabbit antibody (ABCAM, Cambridge, MA USA) was added to the samples for 1 h, followed by the corresponding Alexa555 anti-rabbit (Molecular Probes, Invitrogen, Carlsbad, CA, USA) for 45 min. Finally, cells were stained with DAPI conjugated with antifade.

For the triple staining with the CK/PRDX1 (Santa Cruz, USA)/TYROBP (Abcam) or the CK/NFYA (Abcam)/YWHAB (Santa Cruz), cells were initially incubated with acetone/methanol (9:1) solution for 10 min and, then, washed 3 times with PBS. Blocking solution with 10% (v/v) FBS in PBS for 1 h was used to eliminate non-specific binding. Consequently, slides were incubated with the corresponding primary and secondary antibodies. Zenon technology (FITC-conjugated IGg1 antibody) (Molecular Probes, Invitrogen) was used for CK detection with the A45-B/B3 antibody. Zenon antibodies were prepared within 30 min before use. Cells were also stained with DAPI conjugated with antifade.

Positive and negative controls were used in each experiment (Additional file 1: Figure S1, Additional file 2: Figure S2), using breast cancer cell lines’ cytospins by omitting one of the first antibodies. Therefore, each experiment included three different negative controls and one positive for all the antibodies.

Slides were, then, analyzed with confocal laser scanning microscopy, while the quantification of the protein expression was performed using the software platform of the ARIOL microscope. The platform automatically measured the intensity per pixel of each fluorochrome in every distinct isolated CTC.

Statistical analysis of the clinical data

All statistical tests were performed at the 5% level of significance. SPSS version 20 (SPSS Inc., Chicago, IL) statistical software was used for the analysis. Overall survival (OS) was defined as the time from treatment initiation, until death from any cause. Progression-free survival (PFS) was defined from the enrolment to the study until disease relapse or death whatever occurred first. Kaplan-Meier curves and Cox regression analysis for PFS and OS were compared using the log-rank test to provide a univariate and multivariate assessment of the prognostic value of selected clinical risk factors.

Immunohistochemistry: JUNB staining in tissue samples

Eleven primary and four metastatic tumors were available for JUNB staining from the same cohort of patients. One patient had available primary breast tumor and pleural metastasis. Three-micrometer-thick formalin-fixed paraffin-embedded (FFPE) tissue sections on charged glass slides were deparaffinized in xylene and rehydrated in a graded series of ethanols. Epitope retrieval was heat-induced in a vegetable steamer treated for 45 min in a solution buffered with citric acid (pH 6). Endogenous peroxidase was blocked by applying UltraVision Hydrogen Peroxide Block (Thermo Scientific, Waltham, MA) for 15 min. Nonspecific protein-binding sites were blocked with UltraVision Protein Block (Thermo Scientific, Waltham, MA) for 5 min. Sections were stained with Jun B mouse monoclonal antibody clone C-11 (1:50) (Santa Cruz, CA, USA) for 1 h at room temperature. Immunodetection was performed using UltraVision Quanto Detection System HRP Polymer DAB (Thermo Scientific, Waltham, MA, USA) according to the manufacturer’s instructions. 3,3′Diaminobenzidine Quanto Chromogen (Thermo Scientific, Waltham, MA) was used as chromogen. Slides were counterstained with hematoxylin.

Samples’ histoscore (H-score) was calculated by a semi-quantitative assessment of both the intensity of staining (graded as 0, no-staining; 1, weak; 2, medium; or 3, strong) and the percentage of positive cells. The percentage of cells at each staining intensity level was calculated, and an H-score was assigned summing the individual H-scores for each intensity level using the following formula: [1 × (% cells 1+) + 2 × (% cells 2+) + 3 × (% cells 3+)]. The range of possible scores was from 0 to 300.

Results

Bioinformatics results and functional enrichment evaluation

Based on our previous work [22], we performed a cross-platform comparative study of integrated public datasets, which revealed gene signatures exhibiting differentially overexpression patterns. In particular, the integrated datasets were of different “origin” (peripheral blood and tissue) and disease status (healthy and cancerous), and three different comparisons were performed (Fig. 1): cancer versus healthy tissue samples (Comparison 1, C1), cancer versus healthy blood samples (C2), and cancer tissue versus healthy blood samples (C3). These comparative analyses aimed to reveal persisting differences in the expression profile of certain genes in healthy and cancer status and possibly identify a limited set of candidate CTC biomarkers. For the C1 comparison, a total of 3725 genes were identified, exhibiting significant differential overexpression in cancer cells over control tissue samples. In the same way, 79 overexpressed genes were derived from the C2 comparison in cancerous peripheral blood (PB) samples and 245 overexpressed genes were extracted from the C3 comparison in cancer tissue samples over control PB samples. Based on our hypothesis, the peripheral blood carries information regarding both the primary and metastasis, as well as other cancer-induced alterations. We considered the intersection (C1∩C2∩C3) of the three gene signatures as a target of our biological assessment, since it includes only genes overexpressed in PB, related to cancer-associated factors, eliminating genes overexpressed in normal PB. This common set includes 24 genes, and CXCR4 plays a central role. A functional enrichment evaluation was further performed. Functional enrichment evaluation is generally recognized as a secondary analysis on large gene sets, resulted from high-throughput genomic methods. This step allowed the validation of the gene signature status, based on the biological significance of the 24 extracted genes. We further analyzed each of the 24 genes for their “intrinsic” properties and the 24-gene signature as whole according to their direct or indirect link with (i) CXCR4 pathway and (ii) CXCR4 biological features. Relevant biological information was extracted from four databases as follows: (a) pathways were acquired from the G2SBC (Genes-to-Systems Breast Cancer) [41], KEGG (Kyoto Encyclopedia of Genes and Genomes) [42], Gene Set Enrichment Analysis (GSEA) [43], and WebGestalt (WEB-based GEne SeT AnaLysis Toolkit) [44] in the category of “Pathway Commons” (PC); (b) molecular alterations in breast cancer, and the shortest pathways from CXCR4 to other genes, were delivered from G2SBC; (c) gene ontology (GO) terms were yielded from G2SBC and GSEA; and (d) oncogenic molecular signature overlaps were computed from GSEA. Only genes participating in enriched biological terms at the P ≤ 0.05 level, after multiple test correction [38], were considered as candidate genes and, if they shared multiple biological functions with CXCR4, were prioritized for further analysis in biological samples. Five genes, named JUNB, YWHAB, TYROBP, NFYA, and PRDX1, were evaluated as highly ranked genes, which could be tested in samples from mBC patients and cell lines (Fig. 1).

Fig. 1
figure 1

Bioinformatics and functional enrichment analysis for biomarker discovery. Twenty-four genes were obtained from a data-driven computational procedure (data integration, cross-site comparisons) and were further explored by using functional enrichment analysis. CXCR4 plays a central role in this analysis. Based on four databases (G2SBC, WebGestalt, GSEA, and KEGG), significant biological information was validated and five genes of the 24-signature were prioritized according to their direct or indirect association with the CXCR4 pathway. These six molecules (CXCR4, JUNB, YWHAB, TYROBP, NFYA, and PRDX1) were experimentally evaluated in biological samples

Full size image

Expression pattern of CXCR4 and JUNB in patients’ CTCs

Quantification of the mean intensity per pixel of each fluorochrome by the ARIOL system provided an expression pattern of every investigated protein in all the examined breast cancer cells lines from different subtypes (Luminal: MCF7, HER-positive: SKBR3, basal-like: MDA-MB 231). In addition, the expression of each protein in normal donors’ PBMCs was quantified and compared to patients’ monocytes. Hela cells were used in the analysis, because they were suggested by the manufacturers as positive controls for some of the antibodies (Fig. 2).

Fig. 2
figure 2

Expression of cytokeratin, CXCR4, and JUNB in CTCs isolated from breast cancer patients. a Quantification of CXCR4 expression (intensity per pixel) in breast cancer cell lines, Hela cells, PBMCs, and patients’ CTCs, using ARIOL system automate software. b Quantification of JUNB expression (intensity per pixel) in breast cancer cell lines, in Hela cells, PBMCs, and patients’ CTCs, using ARIOL system automate software. c Cytospins obtained from metastatic breast cancer patients were triple stained with CK (green), CXCR4 (orange), and JUNB (red) antibodies plus DAPI (blue) and analyzed with ARIOL system (magnification X40). d Representative images from confocal laser scanning microscopy. Patients’ samples were triple stained with CK (green), CXCR4 (red), and JUNB (blue) antibodies plus DAPI (not shown due to the absence of corresponding laser) magnification (× 60)

Full size image

In the breast cancer cell lines MCF7 and SKBR3, the intensity of CXCR4 was statistically increased (p = 0.017 and p = 0.0001 respectively) compared to normal PBMCs (Fig. 2a). In addition, the mean intensity of CXCR4 in normal donors’ PBMCs was 10.62 ± 1.01, while in patients’ blood, the mean intensity was enhanced (25.4 ± 1.59; p = 0.0001; Table 2) compared to normal PBMCs. Furthermore, the mean intensity of isolated CTCs in the patients was 33.57 ± 9.25 which was statistically higher compared to normal donors’ PBMCs (p = 0.05) and to patients’ PBMCs (p = 0.0001). Moreover, the expression of CXCR4 was higher in CTCs compared to MCF7 (p = 0.0027) and to MDA-MB 231 (p = 0.011) (Table 2, Fig. 2a). Setting as the threshold value the mean intensity of CXCR4 in normal donors’ PBMCs, 90% (18 out of 20) of the patients with detectable tumor cells in their blood harvested CTCs positive for CXCR4. In addition, 82.61% of the total isolated CTCs have an intensity higher than the threshold and considered positive for CXCR4.

JUNB expression in CTCs

JUNB intensity was statistically higher in all breast cancer cell lines (SKBR3: p = 0.0001, MCF7: p = 0.001, MDA-MB231: p = 0.0001) compared to normal donors’ PBMCs (Table 2, Fig. 2b). In addition, the intensity of JUNB in isolated patients’ CTCs was statistically higher compared to normal (p = 0.007) and patients’ (p = 0.002) PBMCs. As it is shown in Fig. 2c, d, JUNB expression in CTCs was enhanced compared to patients’ PBMCs, providing a potential biomarker to distinct cancer cells in patients’ blood.

Furthermore, the intensity of JUNB was increased compared to MCF7 (p = 0.009) and MDA-MB231 (p = 0.0008) cells. Using as threshold value the mean intensity of JUNB in normal PBMCs, 65% (13 out of 20) of the patients with detectable tumor cells in their blood harbored CTCs positive for JUNB. In addition, 78.57% of the total isolated CTCs were JUNB-positive.

Expression pattern of TYROBP and PRDX1 in breast cancer cell lines and in patients’ CTCs

Quantification of TYROBP expression in breast cancer cell lines and in PBMCs from normal donors and patients revealed that this protein was downregulated in MCF7 (p = 0.008) cells compared to normal PBMCs. In addition, a significant difference was observed between the mean intensity of TYROBP in patients’ PBMCs and breast cancer cell lines [SKBR3 (p = 0.022), MCF7 (p = 0.007)] (Table 2 and Fig. 3a).

Fig. 3
figure 3

Expression of cytokeratin, TYROBP, and PRDX1 in CTCs isolated from breast cancer patients. a Quantification of TYROBP expression (intensity per pixel) in breast cancer cell lines, Hela cells, PBMCs, and patients’ CTCs, using ARIOL system automate software. b Quantification of PRDX1 expression (intensity per pixel) in breast cancer cell lines, in Hela cells, in PBMCs, and in patients’ CTCs, using ARIOL system automate software. c Cytospins obtained from metastatic breast cancer patients were triple stained with CK (green), PRDX1(orange), and TYROBP (red) antibodies plus DAPI (blue) and analyzed with ARIOL system (magnification × 40). d Representative images from confocal laser scanning microscopy. Patients’ samples were triple stained with CK (green), TYROBP (blue), and PRDX1 (red) antibodies plus DAPI (not shown due to the absence of corresponding laser) and analyzed with confocal laser scanning microscopy (magnification × 60)

Full size image

The mean intensity of TYPORB expression in patients’ CTCs (8.89 ± 3.53) was higher than that of breast cancer cell lines [MCF 7 (3.61 ± 0.1, p = 0.05), SKBR3 (3.97 ± 0.13, p = 0.007), MDA-MB231 (5.08 ± 0.09, p = 0.11)]. It was also higher compared to normal (5.17 ± 0.67, p = 0.322) and to patients’ (4.82 ± 0.21, p = 0.03) PBMCs.

Using as threshold value the intensity of TYROBP in normal PBMCs, 64.3% of the isolated CTCs were positive for TYROBP. In addition, 75% (12 out of 16) of the patients with detectable tumor cells in their blood harvested TYROBP-positive CTCs.

PRDX1 expression in CTCs

PRDX1 was also quantified in the same cohort of patients (Table 2). The results revealed that PRDX1 expression was higher in MDA-MB231 cells (22.0 ± 1.64, p = 0.12) compared to normal PBMCs (15.9 ± 3.62). In addition, the mean PRDX1 expression in CTCs (29.34 ± 24.6) was numerically, but not statistically higher than that observed in breast cancer cell lines and in the normal PBMCs (15.9 ± 3.62) (Table 2, Fig. 3b).

Using as threshold value the mean intensity of PRDX1 in normal PBMCs, 20% of the isolated CTCs were positive for PRDX1. In addition, only 12.5% (2 out of 16) of the patients with detectable tumor cells in their blood had detectable PRDX1-positive CTCs.

Expression pattern of NFYA and YWHAB in breast cancer cell lines and in patients’ CTCs

Evaluation of mean intensity of NFYA revealed that the highest expression among the breast cancer cell lines was observed in MDA-MB 231 [(Table 2), 7.26 ± 0.15)] which was statistically different compared to normal donors’ PBMCs (5.9 ± 0.52, p = 0.016). The intensity of NFYA expression in patients’ CTCs (5.20 ± 1.36) was higher compared to MCF7 (4.93 ± 0.016) and SKBR3 (3.95 ± 0.021) cells, but it was lower than in normal PBMCs (5.9 ± 0.52, p = 0.005; Fig. 4a).

Fig. 4
figure 4

Expression of cytokeratin, NFYA, and YWHAB in CTCs isolated from breast cancer patients. a Quantification of NFYA expression (intensity per pixel) in breast cancer cell lines, Hela cells, PBMCs, and patients’ CTCs, using ARIOL system automate software. b Quantification of YWHAB expression (intensity per pixel) in breast cancer cell lines, in Hela cells, in PBMCs, and in patients’ CTCs, using ARIOL system automate software. c Cytospins obtained from metastatic breast cancer patients were triple stained with CK (green), YWHAB (orange), and NFYA (red) antibodies plus DAPI (blue) and analyzed with ARIOL system (magnification × 40). d Representative images from confocal laser scanning microscopy. Patients’ samples were triple stained with CK (green), NFYA (blue), and YWHAB (red) antibodies plus DAPI (not shown due to the absence of the corresponding laser) and analyzed with confocal laser scanning microscopy (magnification × 60)

Full size image

Using as a threshold the intensity of NFYA expression in normal donors’ PBMCs, only 6.45% of the whole number of isolated CTCs were positive for this protein. In addition, only 14.3% (2 out of 14) of the patients with detectable CK-positive cells were considered positive for NFYA.

YWHAB expression in CTCs

The intensity of YWHAB was also evaluated in breast cancer cell lines and in patients’ CTCs. The mean intensity of YWHAB expression in normal donors’ PBMCs (26.13 ± 10.10) was higher than that observed in MCF7 (7.60 ± 0.12, p = 0.043) and SKBR3 (4.04 ± 0.15, p = 0.058) cell lines. In contrast, MDA-MB231 (26.45 ± 3.52, p = 0.043) cells revealed significant higher expression compared to normal PBMCs.

In patients’ CTCs, the mean intensity of YWHAB expression was 5.70 ± 0.28, which was significantly higher than the intensity observed in SKBR3 (4.04 ± 0.15, p = 0.01) cells, whereas it was lower compared to MCF7 (7.6 ± 0.12, p = 0.0001) and MDA-MB231 (26.45 ± 3.52, p = 0.002) cells. The mean intensity of YWHAB expression in CTCs was numerically but not significantly lower than in normal PBMCs (p = 0.07; Fig. 4b). None of the patients could be considered positive for YWHAB, taking as threshold value the normal donors’ PBMCS mean intensity.

Clinical impact of the evaluated molecules

Analysis of the clinical data of all the examined patients (100) revealed that after a median follow-up period of 24 months (range 1–101), 60 patients died due to disease progression. The median overall survival (OS) was significantly lower, in the group of patients harboring JUNB-positive CTCs (HR = 2.308, p = 0.026, Cox regression and p = 0.02: Kaplan-Meier analysis), compared to patients without JUNB-expression in their CTCs [(17 months (range 0–38) vs 24.5 months (range 0–101); Fig. 5 (A)]. Particularly in the group of patients with JUNB-positive CTCs, 90% died during the follow-up period (9 out of 10 with available follow-up data) while in the cohort of patients without JUNB expression, 75% (45 out of 90) died during the follow-up period. Similarly, the median progression-free survival (PFS) was significantly lower in patients harboring JUNB-positive CTCs compared to patients without JUNB-positive CTCs [3.5 months (range 0–30) vs 9 months (range 0–72), p = 0.015, Kaplan-Meier analysis] (Fig. 5 (B)). Interestingly, the presence of CK/CD45 cells in these patients did not correlate with clinical outcome. Furthermore, statistical analysis of the rest of the examined proteins did not lead to clinical correlation with OS or PFS. Multivariate analysis revealed that JUNB expression in CTCs is an independent prognostic factor (p = 0.016, HR 2.2484) for OS in breast cancer patients.

Fig. 5
figure 5

I (A) Overall survival in patients with JUNB-positive CTCs compared to JUNB-negative patients (p = 0.02, Kaplan-Meier analysis). In the group of patients harboring JUNB-positive CTCs, 90% died during follow-up period vs 75% of deaths in the cohort of patients without JUNB expression. I (B) Kaplan-Meier analysis revealed that the progression-free survival in patients harboring JUNB-positive CTCs was statistically lower compared to patients without JUNB expression (p = 0.015). Patients without JUNB-positive CTCs experienced longer PFS vs JUNB-negative patients [3.5 (range 0–30) vs 9 (range 0–72) months]. II (A) JUNB staining in primary tumor obtained from BC patient. JUNB-positive tumor cells (brown) were very rare in the sample. Inner frame is shown control sample from a squamous cell carcinoma (magnification × 40). II (B) JUNB staining in metastatic tumor obtained from mBC patient. The majority of the tumor cells were positive for JUNB (magnification × 40)

Full size image

Further analysis of the other clinicopathological characteristics revealed that reduced OS was significantly associated with performance status (PS) (p = 0.001, HR = 2.391, Cox regression).

JUNB expression in primary tumors

Eleven primary tumors from the same cohort of patients were screened for JUNB expression. Eight of them were negative for JUNB. Mostly all of the blood samples from these patients were also negative for CTCs. Only one patient in this group harbored JUNB-positive CTCs.

Three primary tumor samples revealed low expression of JUNB (H-scores 5, 6, and 9). One of them had also JUNB-positive CTCs.

On the other hand, all the examined samples (four) from metastatic tumors were positive for JUNB and the H-scores were very high (20, 120, 105, and 140). Three of these patients harbored also JUNB-positive CTCs in their blood.

It is interesting that in a patient with available primary tumor and pleura metastasis, the sample from the primary tumor was completely negative for JUNB expression (H-score 0), while in metastasis the H-score was extreme (120). This patient harbored also JUNB-positive CTCs (Table 3).

Table 3 JUNB expression in primary and metastatic tumors
Full size table

Discussion

It is widely accepted that the detection of CTCs is a valuable parameter for the assessment of prognosis in cancer patients. CTCs can also be an important therapeutic target since in a randomized phase II trial it was observed that in patients with HER2-negative early breast cancer targeting HER2-positive CTCs with trastuzumab after the completion of adjuvant chemotherapy resulted in a significant improvement of patients’ DFS [45]. However, just the enumeration of CTCs is not sufficient for defining prognosis in cancer patients. Conversely, their characterization, which could be based on the expression of different proteins, is critical, since it allows the identification of CTCs undergoing EMT or expressing stem cell markers. In addition, it has been recently reported that certain CTCs’ phenotypes can be related to metastatic procedure [1, 2]. Therefore, it is obvious that there is an unmet need for the development of new biomarkers for their detection and characterization. Bioinformatics analysis, based on gene expression databases, could help to this end.

Current study investigated whether bioinformatics analysis using public databases could define putative biomarkers with clinical relevance in CTCs, isolated from metastatic breast cancer patients. The presented bioinformatics data revealed that 6 genes could be potentially overexpressed in CTCs derived from breast cancer patients. All these genes are related to different CXCR4 signal transduction pathways. CXCR4 is a chemokine receptor that binds CXCL12. It has been shown that CXCR4 is upregulated in a number of tumors, such as breast, melanoma, ovarian, and prostate cancer, whereas its expression is low in normal tissues [46, 47]. It is also interesting to note that CXCR4 is induced by angiogenic factors, such as VEGF, and previous studies from our group had shown that VEGF and VEGFR2 were upregulated in CTCs derived from metastatic breast cancer patients [39]. Interestingly in the present study, 90% of the patients with detectable CTCs in their blood harbored CTCs with CXCR4-positive phenotype. Furthermore, quantification of this expression (intensity per pixel) revealed that CXCR4 is significantly upregulated in patients’ CTCs compared to MCF7 (p = 0.027) and MDA-MB231 (p = 0.011) cells and in comparison to normal PBMCs (p = 0.05). This observation strongly suggests that CXCR4 receptor is implicated to the metastatic process. These results confirm the conclusions of our bioinformatics analysis, and they are in agreement with previous studies regarding the CXCR4 mRNA expression in CTCs from breast cancer patients [24, 48, 49]. CXCR4 expression in these studies was associated with the presence of > 3 positive lymph nodes and with EMT markers such as Vimentin and Snail [24] whereas Mego et al. [48] reported that CXCR4 expression was related to CTCs exhibiting epithelial markers. In addition, CXCR4 was upregulated in CTCs derived from different tumor types such as small cell lung cancer, and this expression was associated with poorer PFS [50]. However, in the current study although 82.6% of the whole number of detected CTCs were CXCR4-positive, this detection was not associated with the clinical outcome. Conversely, patients with JUNB expression in their CTCs revealed poorer OS (p = 0.002) and PFS (p = 0.015) compared to patients without this phenotype. Moreover, multivariate analysis showed that the detection of JUNB in CTCs could be emerged as an independent prognostic factor, associated with poor clinical outcome in terms of OS (p = 0.016). This clinical relevance of JUNB-positive CTCs raises the question of its significance as potential therapeutic target for the elimination of these cells. However, analysis of the different breast cancer patients’ subtypes (hormone receptor-positive, HER2-positive, triple negative) separately did not reveal statistical differences.

The data of the current study are in line with previous studies reporting the role of JUNB in cancer progression, such as Hodgkin’s disease and anaplastic large cell lymphomas [51, 52]. Moreover, it has been shown that JUNB is induced by ALK-NPM, participating the mTOR pathway [53] and is required for cell cycle re-entry, after quiescence, and it cooperates with c-jun for the development of fibrosarcoma [54]. In addition, there are some recent data indicating that JUNB is implicated in the earliest events of development of resistance to kinase inhibitors in breast cancer [30]. However in the current study, JUNB overexpression could not be attributed to the treatment, because the samples were obtained before the initiation of first-line treatment. Combining all these data with our findings strongly supports that JUNB plays a critical role in cancer metastasis. This assumption was reinforced by the analysis of 11 primary tumors from the same group of patients, showing that only three of them revealed low JUNB expression, while the rest were completely negative. However, all the metastatic tumors (four) were strongly positive for JUNB. Furthermore, in one patient with available primary tumor and metastasis, JUNB was dramatically increased (from H-score 0 to H-score 120) in metastatic tissue, denoting the critical role of this molecule in cancer progression.

Our data also indicate an increased intensity of PRDX1 expression in patients’ CTCs compared to normal PBMCs, but this difference was not statistically significant. This observation is in line with previous studies demonstrating that PRDX is increased in the most aggressive [triple negative (TN)] subtype of breast cancer patients [55]. It is interesting that previous studies have shown that PRDX provide a protective role in cancer cells regarding doxorubicin-induced toxicity and they are associated with overall mortality in breast cancer patients [33, 35].

The intensity of TYROBP expression was also increased in CTCs compared to normal PBMCs and MCF7 cells, implying an upregulation of this protein in aggressive metastatic cancer cells. The strong statistical difference of TYROBP intensity between MCF7 and MDA-MB 231 (p = 0.00001) reinforced this assumption.

Conversely, although bioinformatics analysis revealed that NFYA and YWHAB genes’ expression could be upregulated in CTCs, our results did not confirm this assumption. This could be attributed to post-transcriptional or post-translational modifications that prohibit the upregulation of protein levels of NFYA and YWHAB.

Conclusions

Bioinformatics analysis could be a useful tool for the identification of new biomarkers and therapeutic targets. Experimental approach to this analysis confirmed the overexpression of CXCR4, JUNB, and TYROBP. In addition, the presence of JUNB-positive CTCs emerged as an independent prognostic factor for OS in breast cancer patients.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

BC:

Breast cancer

CK:

Cytokeratin

CTCs:

Circulating tumor cells

CXCR4:

C-X-C chemokine receptor type 4

EMT:

Epithelial-mesenchymal transition

NFYA:

Nuclear transcription factor Y subunit alpha

OS:

Overall survival

PBMCs:

Peripheral blood mononuclear cells

PFS:

Progression-free survival

PRDX:

Peroxiredoxin

TYROBP (DAP12):

Tyrosine kinase-binding protein

References

  1. Baccelli I, Schneeweiss A, Riethdorf S, Stenzinger A, Schillert A, Vogel V, et al. Identification of a population of blood circulating tumor cells from breast cancer patients that initiates metastasis in a xenograft assay. Nat Biotechnol. 2013;31(6):539–44.

    Article  CAS  PubMed  Google Scholar 

  2. Hodgkinson CL, Morrow CJ, Li Y, Metcalf RL, Rothwell DG, Trapani F, et al. Tumorigenicity and genetic profiling of circulating tumor cells in small-cell lung cancer. Nat Med. 2014;20(8):897–903.

    Article  CAS  PubMed  Google Scholar 

  3. Cristofanilli M, Budd GT, Ellis MJ, Stopeck A, Matera J, Miller MC, et al. Circulating tumor cells, disease progression, and survival in metastatic breast cancer. N Engl J Med. 2004;351(8):781–91.

    Article  CAS  PubMed  Google Scholar 

  4. Huang Y, Fernandez SV, Goodwin S, Russo PA, Russo IH, Sutter TR, et al. Epithelial to mesenchymal transition in human breast epithelial cells transformed by 17beta-estradiol. Cancer Res. 2007;67(23):11147–57.

    Article  CAS  PubMed  Google Scholar 

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

    Article  PubMed  Google Scholar 

  6. Stathopoulou A, Vlachonikolis I, Mavroudis D, Perraki M, Kouroussis C, Apostolaki S, et al. Molecular detection of cytokeratin-19-positive cells in the peripheral blood of patients with operable breast cancer: evaluation of their prognostic significance. 17. J Clin Oncol. 2002;20(16):3404–12.

    Article  CAS  PubMed  Google Scholar 

  7. Kallergi G, Vetsika EK, Aggouraki D, Lagoudaki E, Koutsopoulos A, Koinis F, et al. Evaluation of PD-L1/PD-1 on circulating tumor cells in patients with advanced non-small cell lung cancer. Ther Adv Med Oncol. 2018;10:1758834017750121.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Messaritakis I, Politaki E, Kotsakis A, Dermitzaki EK, Koinis F, Lagoudaki E, et al. Phenotypic characterization of circulating tumor cells in the peripheral blood of patients with small cell lung cancer. PLoS One. 2017;12(7):e0181211.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Lobodasch K, Frohlich F, Rengsberger M, Schubert R, Dengler R, Pachmann U, et al. Quantification of circulating tumour cells for the monitoring of adjuvant therapy in breast cancer: an increase in cell number at completion of therapy is a predictor of early relapse. Breast. 2007;16(2):211–8.

    Article  PubMed  Google Scholar 

  10. Xenidis N, Perraki M, Kafousi M, Apostolaki S, Bolonaki I, Stathopoulou A, et al. Predictive and prognostic value of peripheral blood cytokeratin-19 mRNA-positive cells detected by real-time polymerase chain reaction in node-negative breast cancer patients. J Clin Oncol. 2006;24(23):3756–62.

    Article  CAS  PubMed  Google Scholar 

  11. Agelaki S, Dragolia M, Markonanolaki H, Alkahtani S, Stournaras C, Georgoulias V, et al. Phenotypic characterization of circulating tumor cells in triple negative breast cancer patients. Oncotarget. 2017;8(3):5309–22.

    Article  PubMed  Google Scholar 

  12. Ignatiadis M, Rothe F, Chaboteaux C, Durbecq V, Rouas G, Criscitiello C, et al. HER2-positive circulating tumor cells in breast cancer. PLoS One. 2011;6(1):e15624.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Kallergi G, Agelaki S, Papadaki MA, Nasias D, Matikas A, Mavroudis D, et al. Expression of truncated human epidermal growth factor receptor 2 on circulating tumor cells of breast cancer patients. Breast Cancer Res. 2015;17:113.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Yu M, Bardia A, Aceto N, Bersani F, Madden MW, Donaldson MC, et al. Cancer therapy. Ex vivo culture of circulating breast tumor cells for individualized testing of drug susceptibility. Science. 2014;345(6193):216–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Medicine KJ. Cancer's circulation problem. Science. 2010;327(5969):1072–4.

    Article  Google Scholar 

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Pantel K, Brakenhoff RH, Brandt B. Detection, clinical relevance and specific biological properties of disseminating tumour cells. Nat Rev Cancer. 2008;8(5):329–40.

    Article  CAS  PubMed  Google Scholar 

  18. Kallergi G, Konstantinidis G, Markomanolaki H, Papadaki MA, Mavroudis D, Stournaras C, et al. Apoptotic circulating tumor cells in early and metastatic breast cancer patients. 1. Mol Cancer Ther. 2013;12(9):1886–95.

    Article  CAS  PubMed  Google Scholar 

  19. Kallergi G, Papadaki MA, Politaki E, Mavroudis D, Georgoulias V, Agelaki S. Epithelial to mesenchymal transition markers expressed in circulating tumour cells of early and metastatic breast cancer patients. 5. Breast Cancer Res. 2011;13(3):R59.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Papadaki MA, Kallergi G, Zafeiriou Z, Manouras L, Theodoropoulos PA, Mavroudis D, et al. Co-expression of putative stemness and epithelial-to-mesenchymal transition markers on single circulating tumour cells from patients with early and metastatic breast cancer. BMC Cancer. 2014;14:651.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Theodoropoulos PA, Polioudaki H, Agelaki S, Kallergi G, Saridaki Z, Mavroudis D, et al. Circulating tumor cells with a putative stem cell phenotype in peripheral blood of patients with breast cancer. Cancer Lett. 2010;288(1):99–106.

    Article  CAS  PubMed  Google Scholar 

  22. Sfakianakis S, Bei ES, Zervakis M, Vassou D, Kafetzopoulos D. On the identification of circulating tumor cells in breast cancer. IEEE J Biomed Health Inform. 2014;18(3):773–82.

    Article  PubMed  Google Scholar 

  23. Teixido J, Martinez-Moreno M, Diaz-Martinez M, Sevilla-Movilla S. The good and bad faces of the CXCR4 chemokine receptor. Int J Biochem Cell Biol. 2018;95:121–31.

    Article  CAS  PubMed  Google Scholar 

  24. Markiewicz A, Ksiazkiewicz M, Welnicka-Jaskiewicz M, Seroczynska B, Skokowski J, Szade J, et al. Mesenchymal phenotype of CTC-enriched blood fraction and lymph node metastasis formation potential. PLoS One. 2014;9(4):e93901.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Gokulnath M, Swetha R, Thejaswini G, Shilpa P, Selvamurugan N. Transforming growth factor-beta1 regulation of ATF-3, c-Jun and JunB proteins for activation of matrix metalloproteinase-13 gene in human breast cancer cells. Int J Biol Macromol. 2017;94(Pt A:370–7.

    Article  CAS  PubMed  Google Scholar 

  26. Pei H, Guo Z, Wang Z, Dai Y, Zheng L, Zhu L, et al. RAC2 promotes abnormal proliferation of quiescent cells by enhanced JUNB expression via the MAL-SRF pathway. Cell Cycle. 2018;17(9):1115–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Sundqvist A, Morikawa M, Ren J, Vasilaki E, Kawasaki N, Kobayashi M, et al. JUNB governs a feed-forward network of TGFbeta signaling that aggravates breast cancer invasion. Nucleic Acids Res. 2018;46(3):1180–95.

    Article  CAS  PubMed  Google Scholar 

  28. Gervasi M, Bianchi-Smiraglia A, Cummings M, Zheng Q, Wang D, Liu S, et al. JunB contributes to Id2 repression and the epithelial-mesenchymal transition in response to transforming growth factor-beta. J Cell Biol. 2012;196(5):589–603.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Lian S, Shao Y, Liu H, He J, Lu W, Zhang Y, et al. PDK1 induces JunB, EMT, cell migration and invasion in human gallbladder cancer. Oncotarget. 2015;6(30):29076–86.

    PubMed  PubMed Central  Google Scholar 

  30. Hicks M, Hu Q, Macrae E, DeWille J. JUNB promotes the survival of Flavopiridol treated human breast cancer cells. Biochem Biophys Res Commun. 2014;450(1):19–24.

    Article  CAS  PubMed  Google Scholar 

  31. Gong Y, Cao Y, Song L, Zhou J, Wang C, Wu B. HMGB3 characterization in gastric cancer. Genet Mol Res. 2013;12(4):6032–9.

    Article  CAS  PubMed  Google Scholar 

  32. Liu TA, Jan YJ, Ko BS, Chen SC, Liang SM, Hung YL, et al. Increased expression of 14-3-3beta promotes tumor progression and predicts extrahepatic metastasis and worse survival in hepatocellular carcinoma. Am J Pathol. 2011;179(6):2698–708.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. McDonald C, Muhlbauer J, Perlmutter G, Taparra K, Phelan SA. Peroxiredoxin proteins protect MCF-7 breast cancer cells from doxorubicin-induced toxicity. Int J Oncol. 2014;45(1):219–26.

    Article  CAS  PubMed  Google Scholar 

  34. Pan X, Li D, Huo J, Kong F, Yang H, Ma X. LINC01016 promotes the malignant phenotype of endometrial cancer cells by regulating the miR-302a-3p/miR-3130-3p/NFYA/SATB1 axis. Cell Death Dis. 2018;9(3):303.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Seibold P, Hall P, Schoof N, Nevanlinna H, Heikkinen T, Benner A, et al. Polymorphisms in oxidative stress-related genes and mortality in breast cancer patients--potential differential effects by radiotherapy? Breast. 2013;22(5):817–23.

    Article  PubMed  Google Scholar 

  36. Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30(1):207–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Tusher VG, Tibshirani R, Chu G. Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci U S A. 2001;98(9):5116–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Benjamini Y, Drai D, Elmer G, Kafkafi N, Golani I. Controlling the false discovery rate in behavior genetics research. Behav Brain Res. 2001;125(1–2):279–84.

    Article  CAS  PubMed  Google Scholar 

  39. Kallergi G, Markomanolaki H, Giannoukaraki V, Papadaki MA, Strati A, Lianidou ES, et al. Hypoxia-inducible factor-1alpha and vascular endothelial growth factor expression in circulating tumor cells of breast cancer patients. Breast Cancer Res. 2009;11(6):R84.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Meng S, Tripathy D, Frenkel EP, Shete S, Naftalis EZ, Huth JF, et al. Circulating tumor cells in patients with breast cancer dormancy. Clin Cancer Res. 2004;10(24):8152–62.

    Article  PubMed  Google Scholar 

  41. Mosca E, Alfieri R, Merelli I, Viti F, Calabria A, Milanesi L. A multilevel data integration resource for breast cancer study. BMC Syst Biol. 2010;4:76.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Wang J, Vasaikar S, Shi Z, Greer M, Zhang B. WebGestalt 2017: a more comprehensive, powerful, flexible and interactive gene set enrichment analysis toolkit. Nucleic Acids Res. 2017;3;45(W1):W130–W137.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102(43):15545–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Ogata H, Goto S, Sato K, Fujibuchi W, Bono H, Kanehisa M. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 1999;27(1):29–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Georgoulias V, Bozionelou V, Agelaki S, Perraki M, Apostolaki S, Kallergi G, et al. Trastuzumab decreases the incidence of clinical relapses in patients with early breast cancer presenting chemotherapy-resistant CK-19mRNA-positive circulating tumor cells: results of a randomized phase II study. Ann Oncol. 2012;23(7):1744–50.

    Article  CAS  PubMed  Google Scholar 

  46. Balkwill F. Cancer and the chemokine network. Nat Rev Cancer. 2004;4(7):540–50.

    Article  CAS  PubMed  Google Scholar 

  47. Muller A, Homey B, Soto H, Ge N, Catron D, Buchanan ME, et al. Involvement of chemokine receptors in breast cancer metastasis. Nature. 2001;410(6824):50–6.

    Article  CAS  PubMed  Google Scholar 

  48. Mego M, Cholujova D, Minarik G, Sedlackova T, Gronesova P, Karaba M, et al. CXCR4-SDF-1 interaction potentially mediates trafficking of circulating tumor cells in primary breast cancer. BMC Cancer. 2016;16:127.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Powell AA, Talasaz AH, Zhang H, Coram MA, Reddy A, Deng G, et al. Single cell profiling of circulating tumor cells: transcriptional heterogeneity and diversity from breast cancer cell lines. PLoS One. 2012;7(5):e33788.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Salgia R, Weaver RW, McCleod M, Stille JR, Yan SB, Roberson S, et al. Prognostic and predictive value of circulating tumor cells and CXCR4 expression as biomarkers for a CXCR4 peptide antagonist in combination with carboplatin-etoposide in small cell lung cancer: exploratory analysis of a phase II study. Investig New Drugs. 2017;35(3):334–44.

    Article  CAS  Google Scholar 

  51. Devary Y, Gottlieb RA, Lau LF, Karin M. Rapid and preferential activation of the c-Jun gene during the mammalian UV response. Mol Cell Biol. 1991;11(5):2804–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Mathas S, Hinz M, Anagnostopoulos I, Krappmann D, Lietz A, Jundt F, et al. Aberrantly expressed c-Jun and JunB are a hallmark of Hodgkin lymphoma cells, stimulate proliferation and synergize with NF-kappa B. EMBO J. 2002;21(15):4104–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Staber PB, Vesely P, Haq N, Ott RG, Funato K, Bambach I, et al. The oncoprotein NPM-ALK of anaplastic large-cell lymphoma induces JUNB transcription via ERK1/2 and JunB translation via mTOR signaling. Blood. 2007;110(9):3374–83.

    Article  CAS  PubMed  Google Scholar 

  54. Kovary K, Bravo R. The jun and fos protein families are both required for cell cycle progression in fibroblasts. Mol Cell Biol. 1991;11(9):4466–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Karihtala P, Kauppila S, Soini Y, Arja J. Oxidative stress and counteracting mechanisms in hormone receptor positive, triple-negative and basal-like breast carcinomas. BMC Cancer. 2011;11:262.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Funding

This study was supported, regarding consumables and study design, by the OncoSeed Program PN 3175 funded by European Union and the Greek General Secretary of Research and Technology. Cretan Association for Biomedical Research (CABR) provided support in data interpretation. Authors also appreciate funding of this work by the International Scientific Partnership Program ISPP at the King Saud University through ISPP#009 which provide help in data interpretation.Financial support was also provided by Hellenic Oncology Research Group (HORG). This study was also co-financed by the European Regional Development Fund of the European Union and Greek national funds through the Operational Program Competitiveness, Entrepreneurship and Innovation under the call RESEARCH-CREATE-INNOVATE (project code: T1EDK-03505, e-MASS)

Author information

Authors and Affiliations

  1. Laboratory of Τumor Cell Βiology, Medical School, University of Crete, Heraklion, Greece

    Galatea Kallergi & Vassilis Georgoulias

  2. Department of Biochemistry, Medical School, University of Crete, Voutes, 70013, Heraklion, Crete, Greece

    Galatea Kallergi, Nefeli Zacharopoulou, Saad Alkahtani, Saud Alarifi & Christos Stournaras

  3. Department of Oncology, Hematology, Rheumatology, Immunology and Pulmology, University Hospital, Tübingen, Germany

    Vasileia Tsintari

  4. Computational BioMedicine Laboratory, Institute of Computer Science, Foundation for Research and Technology, Heraklion, Greece

    Stelios Sfakianakis

  5. Digital Image and Signal Processing Laboratory, School of Electrical and Computer Engineering, Technical University of Crete, Chania, Greece

    Ekaterini Bei & Michalis Zervakis

  6. Department of Pathology, University General Hospital of Heraklion, Heraklion, Crete, Greece

    Eleni Lagoudaki & Anastasios Koutsopoulos

  7. Department of Zoology, Science College, King Saud University, Riyadh, Saudi Arabia

    Saad Alkahtani & Saud Alarifi

  8. Hellenic Oncology Research Group (HORG), Athens, Greece

    Galatea Kallergi & Vassilis Georgoulias

Authors
  1. Galatea Kallergi
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Vasileia Tsintari
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Stelios Sfakianakis
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Ekaterini Bei
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Eleni Lagoudaki
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Anastasios Koutsopoulos
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Nefeli Zacharopoulou
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Saad Alkahtani
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Saud Alarifi
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Christos Stournaras
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Michalis Zervakis
    View author publications

    You can also search for this author in PubMed Google Scholar

  12. Vassilis Georgoulias
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

GK participated in the design and coordination of the study. She performed IF experiments. She performed the data interpretation, and she also drafted the manuscript. VT performed most of the immunofluorescence experiments. SA participated to biostatistics analysis and help drafted the manuscript. EB performed the functional enrichment analysis. EL and AK performed the experiments and the evaluation of the primary tumors. NZ and CS provided general support and help to draft the manuscript. SA and SA provided technical and economic support. MZ provided general support and participated in the study design and biostatistics analysis. VG provided general support and participated in the study design. He also participated in the manuscript preparation. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Galatea Kallergi.

Ethics declarations

Ethics approval and consent to participate

The protocol has been approved by the Ethics and Scientific Committees of our institution, and all patients and healthy blood donors gave their informed consent to participate in the study.

Consent for publication

Our study does not content any individual person’s data.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Figure S1. Positive and negative controls of CK/CD45 and CK/JUNB/CXCR4 stainings. (IA): Positive controls for CK/CD45 staining: Cytospins with SKBR3 cells spiked in normal donor’s PBMCs were stained with CK (green) anti-mouse, Alexa 488 anti-mouse, CD45 (blue) anti-rabbit, and Alexa 633 anti-rabbit antibodies. (IB): Negative controls for CK/CD45 staining: Cytospins were stained with all the above antibodies except the primary CK anti-mouse antibody. (IIA) Positive controls for CK/JUNB/CXCR4 staining: Cytospins with SKBR3 cells were stained with CK (green), JUNB (blue), CXCR4 (red) antibodies, and the corresponding fluorochromes. (IB) Negative controls for CK/JUNB/CXCR4 staining: Cytospins were stained with all the corresponding antibodies except the primary CXCR4 antibody. (IIC) Negative controls for CK/JUNB/CXCR4 staining: Cells were stained with all the corresponding antibodies except the primary JUNB antibody. (TIF 4960 kb)

Additional file 2:

Figure S2. (IA) Positive controls for CK/TYROBP/PRDX1 staining: MDA-MB 231 cells were stained with CK (green), TYROBP (blue), PRDX1 (red) antibodies, and the corresponding fluorochromes. (IB) Negative controls for CK/TYROBP/PRDX1 staining: Cells were stained with all the corresponding antibodies except the primary TYROBP antibody. (IC) Negative controls for CK/TYROBP/PRDX1 staining: Cells were stained with all the corresponding antibodies except the primary PRDX1 antibody. (IIA) Positive controls for CK/NFYA/YWHAB staining: MDA-MB 231 cells were stained with CK (green), NFYA (blue), YWHAB (red) antibodies, and the corresponding fluorochromes. (IIB) Negative controls for CK/NFYA/YWHAB staining: Cells were stained with all the corresponding antibodies except the primary NFYA antibody. (IIC) Negative controls for CK/NFYA/YWHAB staining: Cells were stained with all the corresponding antibodies except the primary YWHAB antibody. (TIF 4741 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Kallergi, G., Tsintari, V., Sfakianakis, S. et al. The prognostic value of JUNB-positive CTCs in metastatic breast cancer: from bioinformatics to phenotypic characterization. Breast Cancer Res 21, 86 (2019). https://doi.org/10.1186/s13058-019-1166-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13058-019-1166-4

Keywords

Breast Cancer Research

ISSN: 1465-542X

Contact us