Next Article in Journal
Regulation of Secondary Metabolism in the Penicillium Genus
Next Article in Special Issue
Sevoflurane and Desflurane Exposure Enhanced Cell Proliferation and Migration in Ovarian Cancer Cells via miR-210 and miR-138 Downregulation
Previous Article in Journal
Differential Co-Expression Analyses Allow the Identification of Critical Signalling Pathways Altered during Tumour Transformation and Progression
Previous Article in Special Issue
Anti-Stem Cell Property of Pterostilbene in Gastrointestinal Cancer Cells
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 21
Issue 24
10.3390/ijms21249460
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Increased Stromal Infiltrating Lymphocytes Are Associated with the Risk of Disease Progression in Mesenchymal Circulating Tumor Cell-Positive Primary Breast Cancer Patients

by
Bozena Smolkova
1,†,
Zuzana Cierna
2,3,†,
Katarina Kalavska
4,5,
Svetlana Miklikova
1,
Jana Plava
1,
Gabriel Minarik
6,
Tatiana Sedlackova
7,8,
Dana Cholujova
1,
Paulina Gronesova
1,
Marina Cihova
1,
Karolina Majerova
1,
Marian Karaba
9,
Juraj Benca
9,10,
Daniel Pindak
9,11,
Jozef Mardiak
4 and
Michal Mego
4,5,*
1
Department of Molecular Oncology, Cancer Research Institute, Biomedical Research Center of the Slovak Academy of Sciences, Dubravska Cesta 9, 845 05 Bratislava, Slovakia
2
Department of Pathology, Faculty of Medicine, Comenius University, Sasinkova 4, 811 08 Bratislava, Slovakia
3
Department of Pathology, Faculty Hospital, A. Zarnova 11, 917 75 Trnava, Slovakia
4
2nd Department of Oncology, Faculty of Medicine, Comenius University and National Cancer Institute, Klenova 1, 833 10 Bratislava, Slovakia
5
Translational Research Unit, Faculty of Medicine, Comenius University, Klenova 1, 833 10 Bratislava, Slovakia
6
Department of Molecular Biology, Faculty of Natural Sciences, Comenius University in Bratislava, Ilkovicova 6, 842 15 Bratislava, Slovakia
7
Comenius University Science Park, Ilkovicova 8, 841 04 Bratislava, Slovakia
8
Geneton Ltd., Ilkovicova 8, 841 04 Bratislava, Slovakia
9
Department of Oncosurgery, National Cancer Institute, Klenova 1, 83310 Bratislava, Slovakia
10
Department of Medicine, St. Elizabeth University, Namestie 1. maja 1, 811 02 Bratislava, Slovakia
11
Department of Oncosurgery, Slovak Medical University, Limbova 12, 83103 Bratislava, Slovakia
 
add Show full affiliation list
 
remove Hide full affiliation list
*
Author to whom correspondence should be addressed.
†
Bozena Smolkova and Zuzana Cierna share first authorship.
Int. J. Mol. Sci. 2020, 21(24), 9460; https://doi.org/10.3390/ijms21249460
Submission received: 1 October 2020 / Revised: 7 December 2020 / Accepted: 10 December 2020 / Published: 12 December 2020
(This article belongs to the Special Issue Attacking Cancer Progression and Metastasis 2.0)
Browse Figures
Versions Notes

Abstract

:
Circulating tumor cells (CTCs) and the immune infiltration of tumors are closely related to clinical outcomes. This study aimed to verify the influence of stromal lymphocyte infiltration and the immune context of tumor microenvironment on the hematogenous spread and prognosis of 282 chemotherapy naïve primary BC patients. To detect the presence of mesenchymal CTCs, RNA extracted from CD45-depleted peripheral blood was interrogated for the expression of mesenchymal gene transcripts. Tumor-infiltrating lymphocytes (TILs) were detected in the stromal areas by immunohistochemistry, using CD3, CD8, and CD45RO antibodies. The concentrations of 51 plasma cytokines were measured by multiplex bead arrays. TILs infiltration in mesenchymal CTC-positive patients significantly decreased their progression-free survival (HR = 4.88, 95% CI 2.30–10.37, p < 0.001 for CD3high; HR = 6.17, 95% CI 2.75–13.80, p < 0.001 for CD8high; HR = 6.93, 95% CI 2.86–16.81, p < 0.001 for CD45ROhigh). Moreover, the combination of elevated plasma concentrations of transforming growth factor beta-3 (cut-off 662 pg/mL), decreased monocyte chemotactic protein-3 (cut-off 52.5 pg/mL) and interleukin-15 (cut-off 17.1 pg/mL) significantly increased the risk of disease recurrence (HR = 4.838, 95% CI 2.048–11.427, p < 0.001). Our results suggest a strong impact of the immune tumor microenvironment on BC progression, especially through influencing the dissemination and survival of more aggressive, mesenchymal CTC subtypes.

1. Introduction

Breast cancer (BC) represents a heterogeneous disease that develops in a very complex microenvironment composed of several types of benign cells, some of which are involved in the immunogenicity of BC [1]. Although BC has not been considered as an “immunogenic“ malignancy, the occurrence of tumor-infiltrating lymphocytes (TILs) has been consistently documented, with an impact on prognosis [2]. TILs subpopulations are differently expressed in BC tumors, but their main component seems to be represented by CD3+ T-cells [3]. Higher levels of tumor-infiltrating effector T-cells are associated with better clinical outcomes in selected BC subtypes [4,5,6,7]. TILs are also a valuable predictive biomarker. BC patients with elevated TIL amounts had a significantly increased pathological complete response rate upon neoadjuvant chemotherapy compared to patients with poorly infiltrated tumors [7,8].
On the other hand, Liu et al. recently showed that increased stromal TILs are associated with circulating tumor cells (CTCs) and metastatic relapse in BC patients after neoadjuvant chemotherapy [9]. The metastatic process is responsible for around 90% of cancer-related deaths and remains the biggest challenge of the BC treatment. The presence of CTCs in peripheral blood is a negative prognostic marker for primary as well as metastatic BC [10,11,12,13,14]. The role of CTCs as an independent prognostic factor in progression-free survival (PFS) and overall survival (OS) was confirmed by several retrospective studies [15,16]. A critical step in the metastatic process is an activation of epithelial-to-mesenchymal transition (EMT) in tumor cells, leading to the expression of mesenchymal traits in epithelial cancer cells [17]. EMT results in several hybrid phenotypes, possessing both epithelial and mesenchymal features. Therefore, the epithelial cell adhesion molecule (EpCAM)-based methods, which have consistently confirmed the prognostic value of CTCs, can underestimate highly heterogeneous non-epithelial CTC sub-populations, frequently associated with poor prognosis [14,18,19,20,21]. CTCs can enter the blood circulation as single cells or as multicellular clusters composed of homotypic or heterotypic cells. Heterotypic CTC clusters incorporate stromal or immune cells together with cancer cells, which increases their likelihood to efficiently metastasize [22]. The role of neutrophils in enhancing the metastatic potential of CTC-neutrophil clusters as well as their role in the establishment of the metastatic niche was reviewed recently [23]. Due to their dual phenotype and the role in cancer biology, it was hypothesized that their function is dictated in a context-dependent fashion [24].
Among other features, the tumor microenvironment is responsible for immune cell recruitment [25]. However, the contribution of different TIL subpopulations to the biological and clinical tumor behavior remains unclear [26]. Immune cells move into tissues under the influence of specific cytokines, chemokines, and adhesion molecules that exert a high degree of complexity. They can directly stimulate immune effector and stromal cells and enhance anti-tumor immunity. On the contrary, cancer cells can also directly inhibit the immune cell function by decreasing cellular receptors expression or by suppressing the adaptive immune system response. The anti- or pro-tumor action of immune cells is amplified by cytokines released into circulation playing a role in increased invasiveness, progression, and prognosis of cancer. As shown recently, factors secreted by the tumor, stromal, or immune cells may affect the composition of the patient’s serum [27]. Therefore, circulating cytokine levels can become a non-invasive marker of immune derangement.
Based on these and our previous findings, we hypothesize that the TILs location, in addition to their quantity, may contribute differently to BC outcomes with important consequences for their function and prognostic value. To confirm this hypothesis, we analyzed the influence of stromal lymphocyte infiltration and the immune context of tumor microenvironment on the hematogenous spread and prognosis of chemotherapy naïve primary BC patients.

2. Results

Patients’ clinical characteristics are shown in Table 1. Most of the patients were older than 50 years (77%, n = 217), T-stage T1 (68.4%, n = 193), low or intermediate grade (62.8%, n = 174). The majority of them were diagnosed with invasive ductal carcinomas (88.6%, n = 247), out of which most were hormone receptor-positive (86.5%, n = 244) and HER-2/neu negative (85.8%, n = 242). The most frequent molecular subtype was luminal A (53.2%, n = 150), followed by luminal B (21.3%, n = 60) and HER2 positive (14.2%, n = 40) subtype. The triple-negative subtype represented the smallest group (11.3%, n = 32). More than half of the patients were LN negative (64.3%, n = 180) without lymphovascular invasion (76.4%, n = 175). The association of individual clinical characteristics with PFS is shown in Table S1.

2.1. Association of TILs with Clinical Characteristics and Progression-Free Survival

Stromal TILs, namely CD3, CD8, and CD45RO, were evaluated in the formalin-fixed paraffin-embedded (FFPE) tumor tissues by immunohistochemistry (IHC) as a percentage of area occupied by CD3, CD8, and CD45RO positive mononuclear inflammatory cells over the total intratumoral stromal area. Intratumoral TILs defined as lymphocytes in tumor nests directly interacting with tumor cells were not evaluated (Figure 1).
Lymphocyte predominant tumors, defined using cut-off values of 60% and 50%, respectively, were rare. Based on CD3, it was 2.9% (n = 7) and 5.4% (n = 13) of tumors, 2.3% (n = 6) and 6.1% (n = 16) based on CD8, and 21.4% (n = 55) and 25.3% (n = 65) based on CD45RO TILs.
The mean percent ± standard deviations (SD) for individual TILs were 13.0% ± 15.3 for CD3; 12.3% ±14.9 for CD8, and 30.5% ± 27.2 for CD45RO, while medians were 7% (ranging from 1–75%), 5% (ranging from 1–60%), and 20% (ranging from 1–90%), respectively. Individual TILs positively correlated with each other, CD3 and CD8, r = 0.773, p < 0.001; CD3 and CD45RO, r = 0.801, p < 0.001; CD8 and CD45RO, r = 0.681, p < 0.001.
Based on the ROC curve analysis (Figure S1a–c), the values of CD3 above 6% were defined as CD3high, values below 6% as CD3low. Similarly, the values of CD8 above 6% were defined as CD8high, values below 6% as CD8low. Finally, the values of CD45RO above 12.5% were defined as CD45ROhigh, values below 12.5% as CD45ROlow. Those values represent cut-offs for further analysis.
The TILs mean values differed significantly across clinical categories (Table 1). A higher stromal infiltration was uniformly identified in patients with adverse outcomes, e.g., high grade, negative HR status, high Ki-67 proliferation index, and triple-negative tumor subtype. Less pronounced, but still significant differences were identified in IDC, HER2-, and bcl2-positive patients as well as those with LVI.
High CD8 and CD45RO stromal infiltration were associated with shorter PFS (HR 1.83, 95% CI 1.03–3.27, p = 0.040, and HR 2.16, 1.12–4.17, p = 0.022, respectively) (Table S2, Figure 2). The subtype-specific prognostic significance of individual TILs is shown in Table S3.

2.2. Association between TILs and CTCs and Their Prognostic Significance

To determine the mRNA expression of EMT-inducing transcription factors (TF) in CD45 depleted fraction of peripheral blood, we compared the expression levels in patient samples with those of 60 healthy donors (Table S4). Among the patient samples, Twist1 and Slug transcripts were overexpressed in four (1.8%) and 44 (16.5%) samples, respectively (Table S5). Relative to the highest Snail and Zeb1 transcript levels detected in healthy donor samples, none of the patients overexpressed these gene transcripts. In total, CTC EMT positivity was detected in 17.6% (n = 47) of patients. The presence of CTC EMT was associated with decreased PFS (Figure 3) (HR 2.75, 95% CI 1.53–4.94, p = 0.001).
We did not find differences in stromal TILs infiltration between CTC EMT positive and negative patients (Table 1). CTC EMT-positive CD3high (HR 4.88, 95% CI 2.30–10.37, p < 0.001), CD8high (HR 6.17, 95% CI 2.75–13.80, p < 0.001), and CD45ROhigh (HR 6.93, 95% CI 2.86–16.82, p < 0.001) patients had significantly shorter PFS in comparison with all other combinations (Figure 4).

2.3. Association between TIL Infiltration and Plasma Cytokine Levels

Cytokine measurement was done in the subgroup of 147 patients only. Although TILs strongly positively correlated with each other, correlations between plasma cytokine levels and TILs were rather weak, although several of them were significant (Figure 5). The most significant was the correlation between CD45RO and G-CFS (r = 0.242, p = 0.006). For all other identified correlations, the significance levels were close to the borderline value p < 0.05. All CD3, CD8, and CD45RO TILs correlated positively with IL-17, MIP-1α, IL-5, and G-CSF.
All cytokines were dichotomized by the median values. The univariate Cox-proportional hazard model analysis was used for the selection of the most significant ones for multivariate regression analysis. Variables that achieved the p-value threshold p < 0.1 were included in the forward stepwise multivariate models. Three plasma cytokines, namely TGF-β3 (HR = 1.985, 95% CI 0.916–4.304, p = 0.082), MCP-3 (HR = 2.599, 95% CI 1.130–5.977, p = 0.025), and IL-15 (HR = 2.172, 95% CI 0.929–5.077, p = 0.073) fulfilled these criteria. Their combination (elevated plasma concentration of TGF-β3, cut-off 662 pg/mL marked as TGF-β3high, decreased MCP-3, cut-off 52.5 pg/mL marked as MCP-3low, and IL-15, cut-off 17.1 pg/mL marked as IL-15low) was significantly associated with an elevated risk of disease progression (HR = 4.838, 95% CI 2.048–11.427, p < 0.001) (Figure 6).

2.4. Variables with the Most Significant Impact on Progression-Free Survival

The multivariate Cox logistic regression was employed to predict the recurrence probability in primary BC patients controlled for clinical predictor variables (listed in Table 1), stromal TILs infiltration, and cytokine levels.
In addition to the known clinicopathological parameters, namely LN positivity and high Ki-67 proliferation index, plasma cytokines (TGF-β3high, MCP-3low, and IL-15low), and the interaction between stromal CD3, CD8, and CD45RO TILs infiltration and CTC EMT positivity were factors significantly associated with PFS (Table 2).

3. Discussion

The prognostic role of lymphocytic infiltrates in BC was proposed in 1992 by Aaltomma et al. [28]. Since then, retrospective and prospective studies have shown that the presence of TILs is a predictive marker for higher responses to neoadjuvant chemotherapy and better survival, particularly in triple-negative and HER2-positive early BC [7,29]. Moreover, the presence of TILs in BC was identified as an independent predictor of the response to anthracycline/taxane neoadjuvant chemotherapy [30]. TILs were also shown to be an independent prognostic factor for disease-free survival, distant recurrence-free interval, and overall survival in triple-negative BC patients [4]. In the latter and ER+/HER2+ tumors, TILs have been associated with a significant reduction in the relative risk of death from the disease [5]. As shown above, the prognostic utility of TILs is related to the intrinsic subtypes and clinicopathological characteristics. Kurozumi et al. have shown high TILs expression to be a poor prognostic marker in ER-positive patients but a good prognostic marker in ER-negative patients [31]. Therefore, additional research should address the exact immune subsets of TILs and their prognostic utility in specific patient subgroups.
Spatial TILs heterogeneity within the tumor can also influence their prognostic significance. Intratumoral TILs are defined as lymphocytes in direct cell-to-cell contact with carcinoma cells, while stromal TILs are located in the stroma between the carcinoma cells, but do not directly interact with them. TILs located intratumorally and those identified in stromal areas do not differ only in their contact with the cancer cells, but also in the degree of heterogeneity and the reaction to the signals from the microenvironment [32,33]. Importantly, the microenvironment consisting of a plethora of various non-cancerous cell types, extracellular matrix proteins as well as soluble molecules, can direct the fate of disseminated cancer cells, including cancer stem cells [33]. Based on the recommendations of an International TILs working group, the stromal TILs should be taken into account in the TILs evaluation [34]. The pooled data analysis of Loi et al. has shown the prognostic value of high stromal TILs in an early-stage, node-negative subgroup of triple-negative patients, who have low rates of recurrence and death [35]. Similarly, Rathore et al. have shown CD3+, CD4+, and CD8+ intratumoral and stromal TILs to predict favorable survival outcomes in infiltrating ductal carcinoma, where patients with intratumoral CD4+ and stromal CD8+ cells showed the highest survival [36]. Moreover, analyzing the patients treated with adjuvant anthracycline-based chemotherapy, Koletsa et al. have shown the correlation of higher stromal TILs density with a lower risk of relapse [37]. Oppositely, Liu et al. identified an increased stromal TILs and CD4+ T cell infiltration as an unfavorable prognostic factor measured by the rate of metastatic relapse [9]. To sum up, stromal TILs have an important prognostic and predictive value, especially in high-risk clinical BC subtypes. Recently, experts at the 16th St. Gallen Conference recommended their routine reporting in triple-negative patients [38]. Conflicting results published so far have been influenced by several factors, such as heterogeneity in lymphocyte distribution, technical slide-related issues, minimal assessable stroma, and other factors excellently discussed recently [39].
Our present data show the association of increased stromal TILs infiltration with shorter PFS, particularly in CTC EMT positive patients. A significant association between TILs and the presence of CTCs was demonstrated also in primary ovarian cancer patients [40]. In primary invasive BC, the prevalence of CTC was associated with an elevated number of intratumoral/peritumoral Tregs [41]. Increased infiltration of stromal TILs, CD4+, and CD8+ T cells after neoadjuvant chemotherapy was significantly correlated with the CTC presence [9]. There are no available published data assessing the prognostic value of stromal TIL infiltration in the context of hematogenous dissemination in chemotherapy naïve primary BC. To our best knowledge, our study is the first one concerning this association.
Previously, we have shown the association between some plasma cytokines and their receptors with CTCs in peripheral blood of early BC patients [42]. In the present study, we assessed the same extensive number of cytokines out of which a combination of three, namely TGF-β3, MCP-3, and IL-15, was independently associated with an elevated risk of disease progression. The metastatic process is a complex phenomenon in which cytokines are crucial players. TGF-β, one of the EMT triggers, is well documented in several cancers, including BC [43,44,45]. Moreover, platelet-derived TGF-β was shown to activate TGF-β/Smad and NF-κ B pathways in tumor cells, thereby initiating and/or stabilizing their transition into an invasive mesenchymal-like phenotype [46]. MCP-3 (also known as CCL7) is a chemotactic factor and lymphocyte attractant, playing a pivotal role in tumorigenesis. It promotes EMT progression via the TGF-β pathway and facilitates tumor invasion and metastasis [47]. However, it was shown previously, that downregulated serum MCP3 levels in BC could reflect a lowered tumor immune surveillance by eosinophils, impaired maintenance of T-cell memory, and a reduced attraction of leucocyte subsets, which potentially recognize and destroy tumor cells [48]. As a pleiotropic cytokine, IL-15 plays an important role in innate and adaptive immunity and is able to activate the antimetastatic activities of NK cells by mediating the cross-talk with patrolling monocytes [49]. Decreased serum IL-15 levels can be a consequence of the attenuated antitumor response proposed by Jabri and Abadie [50]. Given the high prognostic significance of plasma concentrations of these three cytokines, we hypothesize that they can represent a surrogate marker of tumor immune derangement, driven by the influence of the tumor microenvironment and its immune infiltration. However, the limitation of our study is a relatively low number of peripheral blood samples for cytokine assessment. This is even more pronounced in the subgroup analysis. Therefore, we consider our findings for a combination of these plasma cytokine levels preliminary and propose they should be confirmed by other research groups.

4. Materials and Methods

4.1. Patients

In this translational study (Protocol TRU-SK 002; Chair: M.Mego), 282 primary BC patients with stages I–III, who were undergoing definitive surgery, were enrolled between March 2012 and February 2015. From each patient, peripheral blood for CTCs detection and cytokine assessment was obtained. The corresponding paraffin-embedded tumor tissues were collected for the TILs expression examination. Each patient was given a complete diagnostic evaluation to exclude the presence of distant metastasis. Patients with concurrent malignancy other than non-melanoma skin cancer in the previous 5 years were excluded as well. Clinicopathological data including information on age, tumor stage, histology, regional lymph node involvement, hormone receptor status, and HER2 status were also recorded.
The study was approved by the Institutional Review Board (IRB) of the National Cancer Institute of Slovakia. Healthy donors (N = 60) were age-matched women without BC who were recruited and consented according to the IRB-approved protocol.

4.2. Detection of CTCs

CTCs were detected in peripheral blood by a quantitative real-time polymerase chain reaction (qRT-PCR) based assay utilizing CD45 positive (CD45+) cell depletion for CTCs enrichment, as described previously [51]. Peripheral blood was subjected to CD45 depletion using the RossetteSepTM kit (Stem Cell Technologies, Vancouver, Canada) according to the manufacturer’s instructions. CD45-depleted cells were mixed with the TRIzolVR LS Reagent (Invitrogen Corporation, Carlsbad, CA, USA) and stored at −80 °C until RNA extraction. RNA from CD45-depleted cells was extracted with the TRIzolVR LS Reagent (Invitrogen Corporation, Carlsbad, CA, USA). The RNA concentration was determined by absorbance readings at 260 nm. RNA extracted from HeLa (ATCC® CCL-2™), HCT 116 (ATCC® CCL-247™), MCF7 (ATCC® HTB-22™), and MDA-MB-231 (ATCC® CRM-HTB-26™) cell lines were used as positive controls.
Direct quantitative reverse transcription was used to detect EMT-inducing TF gene transcripts (TWIST, SNAI1, SLUG, and ZEB1). In brief, 1 μL of RNA was used in 20 μL of reaction volume containing 7.75 μL of water, 0.25 μL of One-Step RT-PCR enzyme mix, a combination of Omniscript and Sensiscript (both QIAGEN, Hilden, Germany) in the 1:1 ratio, 10 μL of Maxima Probe/ROX qPCR Master Mix (2X) (Thermo Fisher Scientific), and 1 μL of the assay. The following TaqMan assays were purchased from Life Technologies, USA: TWIST1: Hs00361186_m1; SNAI1: Hs00195591_m1; SNAI2: Hs00161904_m1; ZEB1: Hs01566408_m1; and GAPDH: Hs99999905_m1. Amplicons or probes spanned intron-exon boundaries. Amplification was performed on the Eppendorf Realplex Real-Time PCR system (Eppendorf, Germany) using the following cycling program: 50 °C for 30 min of reverse transcription, 95 °C for 10 min; 40 cycles of 95 °C for 15 s; and 60 °C for 60 s. Target cDNA was quantified using the delta-Ct method with the formula: 2 ^ (Ct target—Ct GAPDH).
Patient samples with higher EMT-associated TF gene transcripts than those of healthy donors were considered as CTC EMT positive. The highest expression levels of the EMT-inducing TF gene transcripts relative to that of GAPDH were 7.5 × 10−4, 3.8 × 10−2, and 1.7 × 10−1 for TWIST1, SNAI1, and ZEB1, while SLUG transcripts were not detected in any of the healthy donor samples. These values were used as a “cut-off” to determine CTC positivity (See Tables S4 and S5) [42].

4.3. Tumor Pathology

A pathology review was conducted at the Department of Pathology, Faculty of Medicine, Comenius University in Bratislava. The TILs assessment included tumor specimens from 271 patients. All specimens were classified according to the WHO Classification of 2012 [52].
According to the tumor histology, one or two representative tumor areas were identified in the hematoxylin and eosin sections. Sections were matched to their corresponding wax blocks (the donor blocks), and 3 mm diameter cores of the tumor were removed from these donor blocks with the multipurpose sampling tool Harris Uni-Core (Sigma-Aldrich, Steinheim, Germany) and inserted into the recipient master block. The recipient block was cut into 5-μm sections, and the sections were transferred to coated slides.
Slides were deparaffinized and rehydrated in a phosphate-buffered saline solution (10 mM, pH 7.2). The tissue epitopes were demasked using the automated water bath heating process in Dako PT Link (Dako, Glostrup, Denmark); the slides were incubated in a pH 6.0 citrate retrieval buffer at 98 °C for 20 min. The slides were subsequently incubated for 60 min at room temperature with the primary mouse monoclonal antibody against CD3 (Dako, M7254), CD8 (Dako, M7103), CD45RO (Dako, M0742) diluted 1:100 (CD3) or 1:200 (CD8, CD45RO). Immunostaining using anti-mouse/anti-rabbit immuno-peroxidase polymer (EnVision FLEX/HRP, Dako, Glostrup, Denmark) was performed for 30 min at room temperature. For visualization, the diaminobenzidine substrate-chromogen solution was used (DAB, Dako, Glostrup, Denmark) for 5 min. Finally, the slides were counterstained with hematoxylin (Sigma-Aldrich, Steinheim, Germany). For the negative control, breast tissue was subjected to the same procedure without staining with the primary antibody.
Tumor cores were independently assessed by a pathologist who was blinded to clinicopathological data. In cases of disagreement, the result was reached by consensus. The evaluation of TILs was performed according to the recommendations published by Salgado et al. [34]. In brief, we scored only stromal TILs as a percentage of the area occupied by CD3, CD8, and CD45RO positive mononuclear inflammatory cells over the total intratumoral stromal area. Intratumoral TILs defined as lymphocytes in tumor nests directly interacting with tumor cells were not evaluated. Altogether 248 samples were successfully scored for CD3, 270 for CD8, and 264 for CD45 expression. For better visualization, we used immunohistochemistry, similar to other published studies [5,30,53,54]. The percentage of positive stromal cells was assessed for all three markers.

4.4. Plasma Cytokines Assessment

In plasma samples of 147 patients, an analysis of 51 plasma cytokines and angiogenic factors was performed as published previously [42]. Human Group I and II cytokines and TGF-β panels were assessed using multiplex bead arrays (Bio-Plex 200 system, Bio-Rad Laboratories, Hercules, CA, USA). Human Group I 27-plex panel included the following targets: IL-1beta, IL-1r alpha, IL-2, IL-4, IL-5, IL-6, IL-7, IL-8, IL-9, IL-10, IL-12 (p70), IL-13, IL-15, IL-17, Basic FGF, Eotaxin, G-CSF, GM-CSF, IFN-gamma, IP-10, MCP-1 (MCAF), MIP-1alpha, MIP-1beta, PDGF-BB, RANTES, TNF-alpha, and VEGF (Bio-Plex Pro Human Cytokine 21-plex Immunoassay, Bio-Rad Laboratories, Hercules, CA, USA). Group II 21-plex, ICAM-1, and VCAM-1 panel contained targets: IL-1alpha, IL-2Ralpha, IL-3, IL-12 (p40), IL-16, IL-18, CTACK, GRO-alpha, HGF, IFN-alpha2, LIF, MCP-3, M-CSF, MIF, MIG, beta-NGF, SCF, SCGF-beta, SDF-1alpha, TNF-beta, and TRAIL, ICAM-1, and VCAM-1 (Bio-Plex Pro Human Cytokine 21-plex Immunoassay, Bio-Rad Laboratories, Hercules, CA, USA). TGF-β1, TGF-β2, and TGF-β3 were analyzed using the Bio-Plex Pro TGF-beta 3-plex Immunoassay (Bio-Rad Laboratories, Hercules, CA, USA). Premixed cytokine standards and samples were diluted following the manufacturer’s instructions and incubated with agitation (300 rpm, RT) with color-coded magnetic beads conjugated with monoclonal antibodies in the 96-well filter plate for 30 min (2 h for the TGF-β assay). As all three TGF-β isoforms are secreted as inactive complexes, samples were first activated with 1 N HCl for 10 min, then neutralized with 1.2 N NaOH/0.5M HEPES (Applichem, Darmstadt, Germany), and assayed immediately after the neutralization step. Following three washes, samples were incubated with a biotinylated detection antibody on a plate shaker (300 rpm agitation, RT) for 30 min in the dark (1 h for TGF-β). Each captured analyte was detected by the addition of streptavidin-phycoerythrin and quantified using a BioPlex suspension array reader (Bio-Rad Laboratories, Hercules, CA, USA) equipped with a 532 nm reporter laser and 635 nm classification laser diode. Cytokine concentrations (pg/mL) were calculated with the Bio-Plex Manager 4.0 software (Bio-Rad Lab, CA USA) using 5-parameter logistic (5PL) curve fitting. Medians of all measured cytokines are listed in Table S6.

4.5. Statistical Analysis

The patients’ characteristics were summarized using the median (range) for continuous variables and frequency (percentage) for categorical variables. The values of cytokines were dichotomized with the cut-off level of the median into two categories: Low (values below median) or high (values above the median). Normality of distribution was tested by the Kolmogorov-Smirnoff test. If normally distributed, sample means were tested by the Student t-test or the analysis of variance (ANOVA) with Bonferroni’s or Tamhane’s corrections, depending on the homogeneity of variance in the univariate analysis. For non-normally-distributed data, nonparametric Mann-Whitney U or the Kruskal-Wallis H test was used. Pearson’s or Spearman’s correlations were applied according to the normality of data. The univariate analyses were performed for categorical variables using χ2 or Fisher’s exact test.
The receiver operator characteristic (ROC) analyses were applied to calculate the cut-off value for individual TILs, providing the highest sensitivity and specificity and to evaluate their prognostic accuracy. The cut-off values for clinical variables were chosen according to clinically significant values.
The median follow-up period was calculated as a median observation time of all patients, as well as of those being still alive at the time of the last follow-up. PFS was calculated as the interval from the date of sampling (mostly date of surgery) to the date of progression, death, or date of the last adequate follow-up. PFS rates were estimated using the Kaplan-Meier product-limit method and the differences between survival curves were evaluated by the Log-rank test.
Estimates of hazard ratios were calculated using the univariate Cox proportional hazard regression analysis. Factors affecting PFS were determined by the multivariate Cox proportional hazard model, applied to estimate the hazard ratio of each covariate and to adjust for potential confounders. Each model included age, clinicopathological characteristics (Table 1), a combination of cytokines preselected in the univariate analysis, and studied TILs expression. A backward model selection was conducted, and the final fitted model is presented.
All presented p-values were two-sided and p < 0.05 was considered significant. Statistical analyses were performed using the IBM SPSS statistics version 23.0 software for Windows (IBM Corp. Released 2015. IBM SPSS Statistics for Windows, Version 23.0. Armonk, NY, USA: IBM Corp.)

5. Conclusions

In the present study, we evaluated the expression of stromal TILs, specifically T lymphocytes (CD3), cytotoxic T lymphocytes (CD8), and memory T lymphocytes (CD45RO) in breast tumor tissues and we correlated their expression in the stromal areas of tumors with the presence of CTC EMT in peripheral blood of patients. Previously, we have shown that abnormalities in T-cell-mediated immunity found in inflammatory CTC positive BC patients could potentially initiate and impact the dissemination of tumor cells [55]. Here, we have demonstrated the possible influence of stromal TILs infiltration on the hematogenous spread in primary BC patients. Moreover, we also interrogated changes in the plasma cytokine profile, which can serve as a surrogate marker of tumor-induced immune derangement. A combination of TGF- βhigh, MCP-3low, and IL-15low at a given cut-off has the potential, after further validation on a bigger sample size, to serve as a non-invasive prognostic biomarker.

Supplementary Materials

The following are available online at https://www.mdpi.com/1422-0067/21/24/9460/s1, Figure S1: ROC curves for CD3, CD8, and CD45RO; Table S1: Univariate Cox logistic regression analysis for individual clinical variables; Table S2: Univariate Cox logistic regression analysis for individual studied TILs and PFS; Table S3: Univariate Cox logistic regression analysis for studied TILs, stratified by clinicopathological categories and PFS; Table S4: qRT-PCR results for control samples; Table S5: qRT-PCR results for patient samples; Table S6: Measured cytokines median levels (pg/mL).

Author Contributions

Conceptualization, M.M. and B.S.; methodology, Z.C., G.M., T.S., D.C., and P.G.; validation, B.S. and M.M.; formal analysis, K.K., G.M., T.S., D.C., and P.G.; investigation, M.M., B.S., and S.M.; resources, M.M., J.M., and S.M.; data curation, B.S., G.M., and T.S.; writing—original draft preparation, S.M., B.S., M.M., M.C., and J.P.; writing—review and editing, S.M., B.S., M.M., M.K., J.B., and D.P.; visualization, B.S., J.P., S.M.; and K.M., supervision, M.M. and J.M.; project administration, M.M.; funding acquisition, M.M., S.M., and B.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Slovak Research and Development Agency (APVV), grant number APVV-16-0010, APVV-16-0178, by ERA-NET EuroNanoMed II INNOCENT and by the Scientific Grant Agency (VEGA), contract numbers 1/0724/11, 1/0044/15, 1/0271/17, and 2/0052/18.

Acknowledgments

We would like to acknowledge Denisa Maňasová for her excellent technical help. We are grateful to all patients for their participation in the study.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

BCBreast cancer
CAFCancer-associated fibroblast
CTCCirculating tumor cell
ECMExtracellular matrix
EMTEpithelial-mesenchymal transition
EpCAMEpithelial cell adhesion molecule
Basic FGFbasic fibroblast growth factor
G-CSFGranulocyte colony-stimulating factor
IDCInvasive ductal carcinoma
IFN-gammaInterferon gamma
ILInterleukin
LVILymphovascular invasion
MCP-3Monocyte chemotactic protein 3
MDSCMyeloid-derived suppressor cell
METMesenchymal-epithelial transition
MIP-1alphaMacrophage inflammatory protein-1α
OSOverall survival
PFSProgression-free survival
TAMTumor-associated macrophage
TGF-betaTransforming growth factor-beta
TILTumor-infiltrating lymphocyte
TNF-alphaTumor necrosis factor-alpha
TNF-betaTumor necrosis factor-beta

References

  1. Harbeck, N.; Penault-Llorca, F.; Cortes, J.; Gnant, M.; Houssami, N.; Poortmans, P.; Ruddy, K.; Tsang, J.; Cardoso, F. Breast cancer. Nat. Rev. Dis. Primers 2019, 5, 66. [Google Scholar] [CrossRef]
  2. Loi, S.; Sirtaine, N.; Piette, F.; Salgado, R.; Viale, G.; Van Eenoo, F.; Rouas, G.; Francis, P.; Crown, J.P.; Hitre, E.; et al. Prognostic and predictive value of tumor-infiltrating lymphocytes in a phase III randomized adjuvant breast cancer trial in node-positive breast cancer comparing the addition of docetaxel to doxorubicin with doxorubicin-based chemotherapy: BIG 02-98. J. Clin. Oncol. 2013, 31, 860–867. [Google Scholar] [CrossRef]
  3. Althobiti, M.; Aleskandarany, M.A.; Joseph, C.; Toss, M.; Mongan, N.; Diez-Rodriguez, M.; Nolan, C.C.; Ashankyty, I.; Ellis, I.O.; Green, A.R.; et al. Heterogeneity of tumour-infiltrating lymphocytes in breast cancer and its prognostic significance. Histopathology 2018, 73, 887–896. [Google Scholar] [CrossRef]
  4. Adams, S.; Gray, R.J.; Demaria, S.; Goldstein, L.; Perez, E.A.; Shulman, L.N.; Martino, S.; Wang, M.; Jones, V.E.; Saphner, T.J.; et al. Prognostic value of tumor-infiltrating lymphocytes in triple-negative breast cancers from two phase III randomized adjuvant breast cancer trials: ECOG 2197 and ECOG 1199. J. Clin. Oncol. 2014, 32, 2959–2966. [Google Scholar] [CrossRef]
  5. Ali, H.R.; Provenzano, E.; Dawson, S.J.; Blows, F.M.; Liu, B.; Shah, M.; Earl, H.M.; Poole, C.J.; Hiller, L.; Dunn, J.A.; et al. Association between CD8+ T-cell infiltration and breast cancer survival in 12,439 patients. Ann. Oncol. 2014, 25, 1536–1543. [Google Scholar] [CrossRef]
  6. Loi, S.; Michiels, S.; Salgado, R.; Sirtaine, N.; Jose, V.; Fumagalli, D.; Kellokumpu-Lehtinen, P.L.; Bono, P.; Kataja, V.; Desmedt, C.; et al. Tumor infiltrating lymphocytes are prognostic in triple negative breast cancer and predictive for trastuzumab benefit in early breast cancer: Results from the FinHER trial. Ann. Oncol. 2014, 25, 1544–1550. [Google Scholar] [CrossRef]
  7. Ingold Heppner, B.; Untch, M.; Denkert, C.; Pfitzner, B.M.; Lederer, B.; Schmitt, W.; Eidtmann, H.; Fasching, P.A.; Tesch, H.; Solbach, C.; et al. Tumor-Infiltrating Lymphocytes: A Predictive and Prognostic Biomarker in Neoadjuvant-Treated HER2-Positive Breast Cancer. Clin. Cancer Res. 2016, 22, 5747–5754. [Google Scholar] [CrossRef] [Green Version]
  8. Issa-Nummer, Y.; Loibl, S.; von Minckwitz, G.; Denkert, C. Tumor-infiltrating lymphocytes in breast cancer: A new predictor for responses to therapy. Oncoimmunology 2014, 3, e27926. [Google Scholar] [CrossRef]
  9. Liu, J.; Xu, Y.; Yu, M.; Liu, Z.; Xu, Y.; Ma, G.; Zhou, W.; Kong, P.; Ling, L.; Wang, S.; et al. Increased Stromal Infiltrating Lymphocytes are Associated with Circulating Tumor Cells and Metastatic Relapse in Breast Cancer Patients After Neoadjuvant Chemotherapy. Cancer Manag. Res. 2019, 11, 10791–10800. [Google Scholar] [CrossRef]
  10. Weigelt, B.; Bosma, A.J.; Hart, A.A.; Rodenhuis, S.; van’t Veer, L.J. Marker genes for circulating tumour cells predict survival in metastasized breast cancer patients. Br. J. Cancer 2003, 88, 1091–1094. [Google Scholar] [CrossRef] [Green Version]
  11. Cristofanilli, M.; Budd, G.T.; Ellis, M.J.; Stopeck, A.; Matera, J.; Miller, M.C.; Reuben, J.M.; Doyle, G.V.; Allard, W.J.; Terstappen, L.W.; et al. Circulating tumor cells, disease progression, and survival in metastatic breast cancer. N. Engl. J. Med. 2004, 351, 781–791. [Google Scholar] [CrossRef] [Green Version]
  12. Molloy, T.J.; Bosma, A.J.; Baumbusch, L.O.; Synnestvedt, M.; Borgen, E.; Russnes, H.G.; Schlichting, E.; van’t Veer, L.J.; Naume, B. The prognostic significance of tumour cell detection in the peripheral blood versus the bone marrow in 733 early-stage breast cancer patients. Breast Cancer Res. 2011, 13, R61. [Google Scholar] [CrossRef] [Green Version]
  13. Krawczyk, N.; Meier-Stiegen, F.; Banys, M.; Neubauer, H.; Ruckhaeberle, E.; Fehm, T. Expression of stem cell and epithelial-mesenchymal transition markers in circulating tumor cells of breast cancer patients. BioMed Res. Int. 2014, 2014, 415721. [Google Scholar] [CrossRef]
  14. Mego, M.; Karaba, M.; Minarik, G.; Benca, J.; Silvia, J.; Sedlackova, T.; Manasova, D.; Kalavska, K.; Pindak, D.; Cristofanilli, M.; et al. Circulating Tumor Cells With Epithelial-to-mesenchymal Transition Phenotypes Associated With Inferior Outcomes in Primary Breast Cancer. Anticancer Res. 2019, 39, 1829–1837. [Google Scholar] [CrossRef]
  15. Liu, M.C.; Mego, M.; Nakamura, S.; Nole, F.; Pierga, J.; Toi, M. Clinical validity of circulating tumor cell (CTC) enumeration in 841 subjects with metastatic breast cancer (MBC). J. Clin. Oncol. 1059, 29 (Suppl. 15), 10592. [Google Scholar] [CrossRef]
  16. Bidard, F.C.; Peeters, D.J.; Fehm, T.; Nole, F.; Gisbert-Criado, R.; Mavroudis, D.; Grisanti, S.; Generali, D.; Garcia-Saenz, J.A.; Stebbing, J.; et al. Clinical validity of circulating tumour cells in patients with metastatic breast cancer: A pooled analysis of individual patient data. Lancet. Oncol. 2014, 15, 406–414. [Google Scholar] [CrossRef]
  17. Beerling, E.; Seinstra, D.; de Wit, E.; Kester, L.; van der Velden, D.; Maynard, C.; Schäfer, R.; van Diest, P.; Voest, E.; van Oudenaarden, A.; et al. Plasticity between Epithelial and Mesenchymal States Unlinks EMT from Metastasis-Enhancing Stem Cell Capacity. Cell Rep. 2016, 14, 2281–2288. [Google Scholar] [CrossRef] [Green Version]
  18. Papadaki, M.A.; Stoupis, G.; Theodoropoulos, P.A.; Mavroudis, D.; Georgoulias, V.; Agelaki, S. Circulating tumor cells with stemness and epithelial-to-mesenchymal transition features are chemoresistant and predictive of poor outcome in metastatic breast cancer. Mol. Cancer Ther. 2019, 18, 437–447. [Google Scholar] [CrossRef] [Green Version]
  19. Guan, X.; Ma, F.; Li, C.; Wu, S.; Hu, S.; Huang, J.; Sun, X.; Wang, J.; Luo, Y.; Cai, R. The prognostic and therapeutic implications of circulating tumor cell phenotype detection based on epithelial–mesenchymal transition markers in the first-line chemotherapy of HER2-negative metastatic breast cancer. Cancer Commun. 2019, 39, 1–10. [Google Scholar] [CrossRef]
  20. Zhang, S.; Wu, T.; Peng, X.; Liu, J.; Liu, F.; Wu, S.; Liu, S.; Dong, Y.; Xie, S.; Ma, S. Mesenchymal phenotype of circulating tumor cells is associated with distant metastasis in breast cancer patients. Cancer Manag. Res. 2017, 9, 691. [Google Scholar] [CrossRef] [Green Version]
  21. Miklikova, S.; Minarik, G.; Sedlackova, T.; Plava, J.; Cihova, M.; Jurisova, S.; Kalavska, K.; Karaba, M.; Benca, J.; Smolkova, B.; et al. Inflammation-Based Scores Increase the Prognostic Value of Circulating Tumor Cells in Primary Breast Cancer. Cancers 2020, 12, 1134. [Google Scholar] [CrossRef]
  22. Aceto, N. Bring along your friends: Homotypic and heterotypic circulating tumor cell clustering to accelerate metastasis. Biomed. J. 2020, 43, 18–23. [Google Scholar] [CrossRef]
  23. Saini, M.; Szczerba, B.M.; Aceto, N. Circulating Tumor Cell-Neutrophil Tango along the Metastatic Process. Cancer Res. 2019, 79, 6067–6073. [Google Scholar] [CrossRef] [Green Version]
  24. Fridlender, Z.G.; Albelda, S.M. Tumor-associated neutrophils: Friend or foe? Carcinogenesis 2012, 33, 949–955. [Google Scholar] [CrossRef] [Green Version]
  25. Fares, J.; Fares, M.Y.; Khachfe, H.H.; Salhab, H.A.; Fares, Y. Molecular principles of metastasis: A hallmark of cancer revisited. Signal Transduct. Target. Ther. 2020, 5, 28. [Google Scholar] [CrossRef]
  26. Garner, H.; de Visser, K.E. Immune crosstalk in cancer progression and metastatic spread: A complex conversation. Nat. Rev. Immunol. 2020, 20, 483–497. [Google Scholar] [CrossRef]
  27. Jabeen, S.; Espinoza, J.A.; Torland, L.A.; Zucknick, M.; Kumar, S.; Haakensen, V.D.; Lüders, T.; Engebraaten, O.; Børresen-Dale, A.-L.; Kyte, J.A.; et al. Noninvasive profiling of serum cytokines in breast cancer patients and clinicopathological characteristics. Oncoimmunology 2018, 8, e1537691. [Google Scholar] [CrossRef] [Green Version]
  28. Aaltomaa, S.; Lipponen, P.; Eskelinen, M.; Kosma, V.M.; Marin, S.; Alhava, E.; Syrjanen, K. Lymphocyte infiltrates as a prognostic variable in female breast cancer. Eur. J. Cancer 1992, 28A, 859–864. [Google Scholar] [CrossRef]
  29. Dushyanthen, S.; Beavis, P.A.; Savas, P.; Teo, Z.L.; Zhou, C.; Mansour, M.; Darcy, P.K.; Loi, S. Relevance of tumor-infiltrating lymphocytes in breast cancer. BMC Med. 2015, 13, 202. [Google Scholar] [CrossRef] [Green Version]
  30. Denkert, C.; Loibl, S.; Noske, A.; Roller, M.; Muller, B.M.; Komor, M.; Budczies, J.; Darb-Esfahani, S.; Kronenwett, R.; Hanusch, C.; et al. Tumor-associated lymphocytes as an independent predictor of response to neoadjuvant chemotherapy in breast cancer. J. Clin. Oncol. 2010, 28, 105–113. [Google Scholar] [CrossRef]
  31. Kurozumi, S.; Matsumoto, H.; Kurosumi, M.; Inoue, K.; Fujii, T.; Horiguchi, J.; Shirabe, K.; Oyama, T.; Kuwano, H. Prognostic significance of tumour-infiltrating lymphocytes for oestrogen receptor-negative breast cancer without lymph node metastasis. Oncol. Lett. 2019, 17, 2647–2656. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Ravelli, A.; Roviello, G.; Cretella, D.; Cavazzoni, A.; Biondi, A.; Cappelletti, M.R.; Zanotti, L.; Ferrero, G.; Ungari, M.; Zanconati, F. Tumor-infiltrating lymphocytes and breast cancer: Beyond the prognostic and predictive utility. Tumor Biol. 2017, 39, 1010428317695023. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Ingangi, V.; Minopoli, M.; Ragone, C.; Motti, M.L.; Carriero, M.V. Role of Microenvironment on the Fate of Disseminating Cancer Stem Cells. Front. Oncol. 2019, 9, 82. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Salgado, R.; Denkert, C.; Demaria, S.; Sirtaine, N.; Klauschen, F.; Pruneri, G.; Wienert, S.; Van den Eynden, G.; Baehne, F.L.; Penault-Llorca, F.; et al. Harmonization of the evaluation of tumor infiltrating lymphocytes (TILs) in breast cancer: Recommendations by an international TILs-working group 2014. Ann. Oncol. 2015, 26, 259–271. [Google Scholar] [CrossRef] [PubMed]
  35. Loi, S.; Drubay, D.; Adams, S.; Pruneri, G.; Francis, P.A.; Lacroix-Triki, M.; Joensuu, H.; Dieci, M.V.; Badve, S.; Demaria, S.; et al. Tumor-Infiltrating Lymphocytes and Prognosis: A Pooled Individual Patient Analysis of Early-Stage Triple-Negative Breast Cancers. J. Clin. Oncol. 2019, 37, 559–569. [Google Scholar] [CrossRef]
  36. Rathore, A.S.; Kumar, S.; Konwar, R.; Makker, A.; Negi, M.P.S.; Goel, M.M. CD3+, CD4+ & CD8+ tumour infiltrating lymphocytes (TILs) are predictors of favourable survival outcome in infiltrating ductal carcinoma of breast. Indian J. Med. Res. 2014, 140, 361–369. [Google Scholar]
  37. Koletsa, T.; Kotoula, V.; Koliou, G.A.; Manousou, K.; Chrisafi, S.; Zagouri, F.; Sotiropoulou, M.; Pentheroudakis, G.; Papoudou-Bai, A.; Christodoulou, C.; et al. Prognostic impact of stromal and intratumoral CD3, CD8 and FOXP3 in adjuvantly treated breast cancer: Do they add information over stromal tumor-infiltrating lymphocyte density? Cancer Immunol. Immunother. 2020, 69, 1549–1564. [Google Scholar] [CrossRef]
  38. Burstein, H.J.; Curigliano, G.; Loibl, S.; Dubsky, P.; Gnant, M.; Poortmans, P.; Colleoni, M.; Denkert, C.; Piccart-Gebhart, M.; Regan, M. Estimating the benefits of therapy for early-stage breast cancer: The St. Gallen International Consensus Guidelines for the primary therapy of early breast cancer 2019. Ann. Oncol. 2019, 30, 1541–1557. [Google Scholar] [CrossRef] [Green Version]
  39. Kos, Z.; Roblin, E.; Kim, R.S.; Michiels, S.; Gallas, B.D.; Chen, W.; van de Vijver, K.K.; Goel, S.; Adams, S.; Demaria, S.; et al. Pitfalls in assessing stromal tumor infiltrating lymphocytes (sTILs) in breast cancer. NPJ Breast Cancer 2020, 6, 17. [Google Scholar] [CrossRef]
  40. Chebouti, I.; Bankfalvi, A.; Friedrich, C.; Kimmig, R.; Kasimir-Bauer, S. Association between tumor infiltrating immune cells, circulating tumor cells in blood and disseminated tumor cells in the bone marrow in patients with primary ovarian cancer. Cancer Res. 2015, 75, 369. [Google Scholar]
  41. Xue, D.; Xia, T.; Wang, J.; Chong, M.; Wang, S.; Zhang, C. Role of regulatory T cells and CD8(+) T lymphocytes in the dissemination of circulating tumor cells in primary invasive breast cancer. Oncol. Lett. 2018, 16, 3045–3053. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Mego, M.; Cholujova, D.; Minarik, G.; Sedlackova, T.; Gronesova, P.; Karaba, M.; Benca, J.; Cingelova, S.; Cierna, Z.; Manasova, D.; et al. CXCR4-SDF-1 interaction potentially mediates trafficking of circulating tumor cells in primary breast cancer. BMC Cancer 2016, 16, 127. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Zarzynska, J.M. Two faces of TGF-beta1 in breast cancer. Mediat. Inflamm. 2014, 2014, 141747. [Google Scholar] [CrossRef] [PubMed]
  44. Hao, Y.; Baker, D.; Ten Dijke, P. TGF-β-Mediated Epithelial-Mesenchymal Transition and Cancer Metastasis. Int. J. Mol. Sci. 2019, 20, 2767. [Google Scholar] [CrossRef] [Green Version]
  45. Panis, C.; Herrera, A.C.; Victorino, V.J.; Aranome, A.M.; Cecchini, R. Screening of circulating TGF-β levels and its clinicopathological significance in human breast cancer. Anticancer Res. 2013, 33, 737–742. [Google Scholar]
  46. Labelle, M.; Begum, S.; Hynes, R.O. Direct signaling between platelets and cancer cells induces an epithelial-mesenchymal-like transition and promotes metastasis. Cancer Cell 2011, 20, 576–590. [Google Scholar] [CrossRef] [Green Version]
  47. Liu, Y.; Cai, Y.; Liu, L.; Wu, Y.; Xiong, X. Crucial biological functions of CCL7 in cancer. PeerJ 2018, 6, e4928. [Google Scholar] [CrossRef]
  48. Carlsson, A.; Wingren, C.; Ingvarsson, J.; Ellmark, P.; Baldertorp, B.; Fernö, M.; Olsson, H.; Borrebaeck, C.A. Serum proteome profiling of metastatic breast cancer using recombinant antibody microarrays. Eur. J. Cancer 2008, 44, 472–480. [Google Scholar] [CrossRef]
  49. Kubo, H.; Mensurado, S.; Gonçalves-Sousa, N.; Serre, K.; Silva-Santos, B. Primary Tumors Limit Metastasis Formation through Induction of IL15-Mediated Cross-Talk between Patrolling Monocytes and NK Cells. Cancer Immunol. Res. 2017, 5, 812–820. [Google Scholar] [CrossRef] [Green Version]
  50. Jabri, B.; Abadie, V. IL-15 functions as a danger signal to regulate tissue-resident T cells and tissue destruction. Nat. Rev. Immunol. 2015, 15, 771–783. [Google Scholar] [CrossRef]
  51. Mego, M.; Mani, S.A.; Lee, B.N.; Li, C.; Evans, K.W.; Cohen, E.N.; Gao, H.; Jackson, S.A.; Giordano, A.; Hortobagyi, G.N.; et al. Expression of epithelial-mesenchymal transition-inducing transcription factors in primary breast cancer: The effect of neoadjuvant therapy. Int. J. Cancer 2012, 130, 808–816. [Google Scholar] [CrossRef] [PubMed]
  52. Lakhani, S.R.; Ellis, I.O.; Schnitt, S.J.; Tan, P.H.; van de Vijver, M.J. WHO Classification of Tumours of the Breast; IARC Press: Lyon, France, 2012. [Google Scholar]
  53. Baker, K.; Lachapelle, J.; Zlobec, I.; Bismar, T.A.; Terracciano, L.; Foulkes, W.D. Prognostic significance of CD8+ T lymphocytes in breast cancer depends upon both oestrogen receptor status and histological grade. Histopathology 2011, 58, 1107–1116. [Google Scholar] [CrossRef] [PubMed]
  54. Chen, Z.; Chen, X.; Zhou, E.; Chen, G.; Qian, K.; Wu, X.; Miao, X.; Tang, Z. Intratumoral CD8(+) cytotoxic lymphocyte is a favorable prognostic marker in node-negative breast cancer. PLoS ONE 2014, 9, e95475. [Google Scholar] [CrossRef] [PubMed]
  55. Mego, M.; Gao, H.; Cohen, E.N.; Anfossi, S.; Giordano, A.; Sanda, T.; Fouad, T.M.; De Giorgi, U.; Giuliano, M.; Woodward, W.A.; et al. Circulating Tumor Cells (CTC) Are Associated with Defects in Adaptive Immunity in Patients with Inflammatory Breast Cancer. J. Cancer 2016, 7, 1095–1104. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Ijms 21 09460 g001 550
Figure 1. Immunohistochemical expression of CD3 (a,d), CD8 (b,e) and CD45RO (c,f) in stromal tumor-infiltrating lymphocytes (TILs). Breast cancer tissue areas with low (a–c), and high stromal TIL (d–f) infiltration are marked by arrows. Original magnification ×400, visualization with 3,3′-diaminobenzidine.
Figure 1. Immunohistochemical expression of CD3 (a,d), CD8 (b,e) and CD45RO (c,f) in stromal tumor-infiltrating lymphocytes (TILs). Breast cancer tissue areas with low (a–c), and high stromal TIL (d–f) infiltration are marked by arrows. Original magnification ×400, visualization with 3,3′-diaminobenzidine.
Ijms 21 09460 g001
Ijms 21 09460 g002 550
Figure 2. The Kaplan-Meier progression-free survival (PFS) estimates for studied stromal TILs expression. In contrast to patients with CD3 infiltration (a) those with CD8high (b) and CD45ROhigh (c) had significantly shorter PFS than those with CD8low or CD45low, respectively (p = 0.037 and p = 0.019 by the Log-rank test).
Figure 2. The Kaplan-Meier progression-free survival (PFS) estimates for studied stromal TILs expression. In contrast to patients with CD3 infiltration (a) those with CD8high (b) and CD45ROhigh (c) had significantly shorter PFS than those with CD8low or CD45low, respectively (p = 0.037 and p = 0.019 by the Log-rank test).
Ijms 21 09460 g002
Ijms 21 09460 g003 550
Figure 3. The Kaplan-Meier PFS estimates for circulating tumor cells (CTC) epithelial-to-mesenchymal transition (EMT). CTC EMT positive patients had significantly shorter PFS (p < 0.001 by the Log-rank test).
Figure 3. The Kaplan-Meier PFS estimates for circulating tumor cells (CTC) epithelial-to-mesenchymal transition (EMT). CTC EMT positive patients had significantly shorter PFS (p < 0.001 by the Log-rank test).
Ijms 21 09460 g003
Ijms 21 09460 g004 550
Figure 4. The Kaplan-Meier PFS estimates for CTC EMT stratified by individual stromal TILs categories. CTC EMT positive patients with (a) CD3high, (b) CD8high, and (c) CD45ROhigh stromal infiltration had significantly shorter PFS than all other combinations (p ≤ 0.001 for CD3, CD8, and CR45RO, respectively, by the Log-rank test).
Figure 4. The Kaplan-Meier PFS estimates for CTC EMT stratified by individual stromal TILs categories. CTC EMT positive patients with (a) CD3high, (b) CD8high, and (c) CD45ROhigh stromal infiltration had significantly shorter PFS than all other combinations (p ≤ 0.001 for CD3, CD8, and CR45RO, respectively, by the Log-rank test).
Ijms 21 09460 g004
Ijms 21 09460 g005 550
Figure 5. Heatmap showing correlation coefficients for significant correlations (highlighted by asterisks) found between CD3, CD8, and CD45RO TILs infiltration (%) and plasma concentrations of selected cytokines (pg/mL).
Figure 5. Heatmap showing correlation coefficients for significant correlations (highlighted by asterisks) found between CD3, CD8, and CD45RO TILs infiltration (%) and plasma concentrations of selected cytokines (pg/mL).
Ijms 21 09460 g005
Ijms 21 09460 g006 550
Figure 6. The Kaplan-Meier PFS estimates for the combination of TGF-β3high, MCP-3low, and IL-15low plasma levels. Patients with this combination of plasma cytokine levels had shorter PFS, p < 0.001 by the Log-rank test.
Figure 6. The Kaplan-Meier PFS estimates for the combination of TGF-β3high, MCP-3low, and IL-15low plasma levels. Patients with this combination of plasma cytokine levels had shorter PFS, p < 0.001 by the Log-rank test.
Ijms 21 09460 g006
Table
Table 1. Mean values of studied TILs stratified by clinicopathological characteristics.
Table 1. Mean values of studied TILs stratified by clinicopathological characteristics.
VariablesCategoriesN%CD3
Mean% ± SD		pCD8
Mean% ± SD		pCD45RO
Mean% ± SD		p
All patients 28210013.00 ± 15.30 12.29 ± 14.89 30.48 ± 27.18
Age (years)≤506523.015.38 ± 17.950.47814.97 ± 17.710.40534.28 ± 30.440.369
>5021777.012.33 ± 14.46 11.45 ± 13.85 29.32 ± 26.08
T-stageT119368.411.14 ± 13.160.10111.81 ± 14.610.36827.60 ± 24.600.086
T2 and more8931.616.97 ± 18.58 13.29 ± 15.52 36.51 ± 31.22
HistologyIDC24787.613.86 ± 15.870.02013.18 ± 15.530.03531.99 ± 27.660.016
Others3512.46.46 ± 7.35 5.68 ± 5.39 18.62 ± 19.63
GradeLow and intermediate17462.89.08 ± 11.38<0.00110.13 ± 13.020.00323.79 ± 22.66<0.001
High10337.219.88 ± 18.48 16.15 ± 17.10 42.05 ± 30.31
N stageN018064.312.27 ± 15.460.12411.80 ± 15.050.24728.61 ± 27.330.046
N+10035.714.31 ± 15.13 12.91 ± 14.68 33.97 ± 27.00
LVI Absent17576.411.91 ± 14.880.06811.95 ± 14.890.02428.05 ± 26.190.034
Present5423.616.72 ± 17.27 14.94 ± 16.13 37.48 ± 29.43
HR status $Negative 3813.524.74 ± 19.34<0.00118.19 ± 16.660.00554.54 ± 29.85<0.001
Positive24486.511.251 ± 13.82 11.31 ± 14.39 26.43 ± 24.55
HER2 statusNegative24285.812.00 ± 14.800.01811.69 ± 14.700.06128.63 ± 26.310.013
Amplified4014.218.43 ± 16.96 15.92 ± 15.72 41.49 ± 29.92
p53Negative17763.011.27 ± 13.690.12411.60 ± 14.960.22728.92 ± 26.390.281
Positive10437.015.65 ± 17.35 13.52 ± 14.84 33.03 ± 28.52
bcl2Negative8028.417.85 ± 18.900.03314.08 ± 15.460.13339.67 ± 31.880.008
Positive20271.611.10 ± 13.23 11.58 ± 14.65 26.90 ± 24.28
Ki-67 &Low18063.88.94 ± 11.35<0.00110.36 ± 14.100.00122.81 ± 21.67<0.001
High10236.219.83 ± 18.43 15.59 ± 15.70 43.56 ± 30.54
Tumor subtypesLuminal A15053.28.87 ± 10.90<0.00110.53 ± 14.130.00622.08 ± 20.80<0.001
Luminal B6021.313.35 ± 16.48 10.51 ± 13.54 30.39 ± 26.80
HER2 positive4014.218.43 ± 16.96 15.92 ± 15.72 41.49 ± 29.92
Triple-negative3211.324.38 ± 20.43 18.66 ± 17.38 53.06 ± 31.64
CTC EMTNegative 22082.412.98 ± 15.330.80612.47 ± 15.200.77429.69 ± 27.230.414
Positive4717.613.07 ± 15.31 11.74 ± 14.60 31.55 ± 26.76
Cytokines *Low risk13391.712.54 ± 14.360.54811.40 ± 13.580.05627.77 ± 24.840.655
High risk128.313.60 ± 19.64 9.40 ± 18.73 27.30 ± 31.98
The total number of samples analyzed in the study was n = 282; solely cases with valid information on individual variables were included in the table; $ negative for both or positive for either with cut-off 1%; & cut-off 20%; * cut-off values: Transforming growth factor beta-3 (TGF-β3) > 662 pg/mL; monocyte chemotactic protein-3 (MCP-3) < 52.5 pg/mL and interleukin-15 (IL-15) < 17.1 pg/mL, data from two of 147 patients were not included, as they did not have successfully measured at least one cytokine used to set the cut-off value. Abbreviations: IDC: Invasive ductal carcinoma; LVI: Lymphovascular invasion; HR: Hormonal receptor.
Table
Table 2. Cox proportional hazard regression models for the association between clinical predictor variables, plasma cytokines, the interaction between individual TILs and CTC EMT, and PFS.
Table 2. Cox proportional hazard regression models for the association between clinical predictor variables, plasma cytokines, the interaction between individual TILs and CTC EMT, and PFS.
ModelVariableHR95% CIp
1LN+6.4462.56–16.24<0.001
Ki-67 > 20% 12.004.22–34.11<0.001
Cytokines *5.9282.12–16.580.001
CD3high and CTC EMT positivity5.2772.09–13.30<0.001
2LN+4.2561.79–10.100.001
Ki-67 > 20% 8.4823.25–22.13<0.001
Cytokines *6.3872.35–17.38<0.001
CD8high and CTC EMT positivity3.6551.49–9.000.005
3LN+3.7961.50–9.580.005
Ki-67 > 20% 7.2512.59–20.30<0.001
Cytokines *5.1721.80–14.870.002
CD45ROhigh and CTC EMT positivity4.9221.83–13.230.002
Categorical variables entered in step 1: Age categories; T-stage; Histological grade; Ki-67, cut-off 20%; HER2 status; HR status, cut-off 1%; Grade; N-stage; * TGF-β3high, MCP-3low, and IL-15low; interaction of CD3, CD8, and CD45RO values with CTC EMT positivity.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Smolkova, B.; Cierna, Z.; Kalavska, K.; Miklikova, S.; Plava, J.; Minarik, G.; Sedlackova, T.; Cholujova, D.; Gronesova, P.; Cihova, M.; et al. Increased Stromal Infiltrating Lymphocytes Are Associated with the Risk of Disease Progression in Mesenchymal Circulating Tumor Cell-Positive Primary Breast Cancer Patients. Int. J. Mol. Sci. 2020, 21, 9460. https://doi.org/10.3390/ijms21249460

AMA Style

Smolkova B, Cierna Z, Kalavska K, Miklikova S, Plava J, Minarik G, Sedlackova T, Cholujova D, Gronesova P, Cihova M, et al. Increased Stromal Infiltrating Lymphocytes Are Associated with the Risk of Disease Progression in Mesenchymal Circulating Tumor Cell-Positive Primary Breast Cancer Patients. International Journal of Molecular Sciences. 2020; 21(24):9460. https://doi.org/10.3390/ijms21249460

Chicago/Turabian Style

Smolkova, Bozena, Zuzana Cierna, Katarina Kalavska, Svetlana Miklikova, Jana Plava, Gabriel Minarik, Tatiana Sedlackova, Dana Cholujova, Paulina Gronesova, Marina Cihova, and et al. 2020. "Increased Stromal Infiltrating Lymphocytes Are Associated with the Risk of Disease Progression in Mesenchymal Circulating Tumor Cell-Positive Primary Breast Cancer Patients" International Journal of Molecular Sciences 21, no. 24: 9460. https://doi.org/10.3390/ijms21249460

APA Style

Smolkova, B., Cierna, Z., Kalavska, K., Miklikova, S., Plava, J., Minarik, G., Sedlackova, T., Cholujova, D., Gronesova, P., Cihova, M., Majerova, K., Karaba, M., Benca, J., Pindak, D., Mardiak, J., & Mego, M. (2020). Increased Stromal Infiltrating Lymphocytes Are Associated with the Risk of Disease Progression in Mesenchymal Circulating Tumor Cell-Positive Primary Breast Cancer Patients. International Journal of Molecular Sciences, 21(24), 9460. https://doi.org/10.3390/ijms21249460

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

Article Metrics

Back to TopTop