NHERF1 and tumor microenvironment: a new scene in invasive breast carcinoma

Background Tumor microenvironment (TME) includes many factors such as tumor associated inflammatory cells, vessels, and lymphocytes, as well as different signaling molecules and extracellular matrix components. These aspects can be de-regulated and consequently lead to a worsening of cancer progression. In recent years an association between the scaffolding protein Na+/H+ exchanger regulatory factor 1 (NHERF1) and tumor microenvironment changes in breast cancer (BC) has been reported. Methods Subcellular NHERF1 localization, vascular endothelial growth factor (VEGF), its receptor VEGFR1, hypoxia inducible factor 1 alpha (HIF-1α), TWIST1 expression and microvessel density (MVD) in 183 invasive BCs were evaluated, using immunohistochemistry on tissue microarrays (TMA). Immunofluorescence was employed to explore protein interactions. Results Cytoplasmic NHERF1(cNHERF1) expression was directly related to cytoplasmic VEGF and VEGFR1 expression (p = 0.001 and p = 0.027 respectively), and inversely to nuclear HIF-1α (p = 0.021) and TWIST1 (p = 0.001). Further, immunofluorescence revealed an involvement of tumor cells with NHERF1 positive staining in neo-vascular formation, suggesting a “mosaic” structure development of these neo-vessels. Survival analyses showed that loss of nuclear TWIST1 (nTWIST1) expression was related to a decrease of disease free survival (DFS) (p < 0.001), while nTWIST1-/mNHERF1+ presented an increased DFS with respect to nTWIST1+/mNHERF1- phenotype (p < 0.001). Subsequently, the analyses of nTWIST1+/cNHERF1+ phenotype selected a subgroup of patients with a worse DFS compared to nTWIST1-/cNHERF1- patients (p = 0.004). Conclusion Resulting data suggested a dynamic relation between NHERF1 and TME markers, and confirmed both the oncosuppressor role of membranous NHERF1 expression and the oncogene activity of cytoplasmic NHERF1. Electronic supplementary material The online version of this article (10.1186/s13046-018-0766-7) contains supplementary material, which is available to authorized users.


Background
Tumor microenvironment (TME) is a complex system in which different cell types, playing an important role in carcinogenesis and tumor-host relationship, co-exist. A bidirectional communication between cells and microenvironment has been highlighted in normal tissue homeostasis and also in the tissues involved in different diseases [1]. In this context the tumor generates and proliferates, surrounded by blood vessels, fibroblasts, inflammatory and immune cells, lymphocytes, different signaling molecules and extracellular matrix. Thus, the tumor can directly affect the microenvironment by the release of extracellular signals and immune escape mechanisms [2]. The interactions are often intricate and their characterization could permit to identify prognostic and predictive biomarkers useful for the development of new pharmacological targets and to provide novel therapeutical opportunities for cancer patients. This is also true for breast cancer (BC), one of the most diffused cancer in women, with a high incidence rate in all countries [3]. Recently, the heterogeneity of cellular and non-cellular components of TME has been underlined, including growth factors, cytokines, RNA, DNA, metabolites and matricellular proteins. They are involved in cancer onset and sustenance by supplying metabolites, energy, drug resistance environment, metastatic, growth and angiogenetic signals as well as evading immune surveillance [4,5]. Alterations and local changes in the concentration of angiogenic factors such as vascular endothelial growth factor (VEGF), its receptor (VEGFR1) and transcriptional targets of hypoxia inducible factor 1 alpha (HIF-1α), may trigger angiogenesis, by increasing microvessel density (MVD) in hypoxic breast TME [6,7]. Among the factors distinguishing the TME, the epithelial-mesenchymal transition (EMT) plays an important role as well, and to date the involvement of some EMT markers has been extensively investigated [8]. TWIST1 is a basic helix-loop-helix transcription factor that contributes to carcinogenesis by triggering EMT, thereby influencing tumor invasion behavior [9]. TWIST1 over-expression was confirmed as having an essential role in tumor initiation, stemness, angiogenesis, invasion, metastasis, and chemo-resistance in a variety of carcinomas [10][11][12][13].
During recent years, our focus has been on the adaptor protein Na + /H + exchanger regulating factor 1 (NHERF1), a scaffold protein, with two PDZ and one EB (Ezrin-binding) domains, able to link different molecules and to drive several signal transduction pathways [14]. Its implication in different tumors such as predictive/ prognostic biomarkers has been extensively investigated [15][16][17], but little is known about NHERF1-TME relation. Our previous studies underlined its possible involvement in the orchestrate TME signaling pathway. In fact, its overexpression is implicated in the increase of invasive phenotype in breast cancer cells, in synergy with exposure to the hypoxic tumor microenvironment and with HIF-1α expression [18]. Furthermore, at membranous and cytoplasmic level, its overexpression and colocalization with VEGFR1, a marker correlated with high metastasis risk and relapse has been evaluated [19]. A similar behavior has been also observed in advanced colorectal cancer (CRC), where overexpression of nuclear NHERF1 in association with hypoxic microenvironment and tumor invasive phenotype has been reported [20]. In a recent study, NHERF1 expression was positively correlated with VEGFR2 expression in CRC, suggesting NHERF1 involvement in disease progression by VEGFR2 signaling pathway [21]. The present study aimed to evaluate the possible involvement of NHERF1 and its association with the major TME markers as VEGF, VEGFR1, HIF-1α, MVD and TWIST1, in 183 invasive BC tumors. Moreover, we investigated possible specific immunophenotypes useful in identifying a subgroup of BC patients with a different clinical outcome, related to these TME markers. A better characterization of TME could allow the identification of new combinations of prognostic biomarkers valuable for clinical management and the development of novel therapeutic strategies.

Patients and tissues
The tissue microarray (TMA) utilized in this retrospective study included 549 samples of invasive BCs from 183 patients. The patients were subjected to primary surgery with nodal dissection at the Institute IRCCS Istituto Tumori "Giovanni. Paolo II" of Bari, Italy, in the years 2002-2003. For each tumor, 3 punches of cancerous tissues were available. The median age of patients was 49 years (range 24-83). For the subgroup of patients with follow-up (64,5%), the median follow-up time was 82 months (range 1-162 months). TNM classification, tumor size, grade, perineoplastic invasion, estrogen receptor (ER), progesterone receptor (PgR), proliferative activity (MIB1) and HER2 status were provided by the Pathology Department of our Institute. The clinicopathological characteristics of the patients are listed in Table 1. Tumors with ER or PgR expression were classified as positive when nuclear staining was found in > 10%. MIB1 nuclear staining was used to assess the proliferative activity, with a cut off value of 20% positive cells to indicate the tumors with MIB1 > 20% as highly proliferating. This cut off is the median value of scores relative to all breast tumors samples analysed during these years within our Institute. HER2 protein expression was studied using a monoclonal antibody (MoAb clone CB11; Novocastra Laboratories, Ltd., Newcastle, UK) and scored in accordance with the HercepΤest scoring system (Food and Drug Administration): 0 indicated no membranous immunoreactivity or < 10% of cells reactive; 1+, an incomplete membranous reactivity in > 10% of cells; 2+, ≥10% of cells with weak to moderate complete membranous reactivity; and 3+, strong and complete membranous reactivity in > 10% of cells. Cytoplasmic immunoreactivity was ignored. Cases scoring 0 and 1+ were classified as negative. HER2 was considered to be positive if immunostaining was 3+ or if a score 2+ showed gene amplification by fluorescence in situ hybridization (FISH). In FISH analyses, each copy of the HER2 gene and its centromere 17 (CEP17) reference was counted. The interpretation was in accordance with the criteria of 2007 ASCO/ CAP guidelines for HER2 testing in BC, therefore positive if the HER2/CEP17 ratio was higher than 2,2 [22].
TMAs were constructed manually from formalin-fixed and paraffin-embedded tissues, according to standard procedures as previously described [15]. In brief, tumor target areas were selected from one donor block per patient, three 0.5 mm cores were punched out and transferred to the recipient TMA blocks. Each sample was arrayed in triplicate to minimize tissue loss and to overcome tumor heterogeneity, thus the three cores were representative of the whole specimen.

Immunohistochemistry (IHC) and evaluation
Four μm-thick slices were cut from the TMA blocks and transferred to slides. The TMA slides were processed and stained by the indirect immunoperoxidase method [23]. In brief, TMA slides were deparaffinised and partially rehydrated through absolute ethanol and 95% ethanol series. Antigen retrieval was performed by the 0. 01 M citrate buffer (pH 6.0) at 98°C in a water bath from a minimum of 20 to a maximum of 45 min. The slides were then allowed to cool for 30 min and the endogenous peroxidase activity was blocked for 10 min with 3% H 2 O 2 . The primary antibodies, diluted in PBS/ BSA 1%, were incubated on the slides at 4°C overnight in a moist chamber. For anti-VEGFR1, 1 h incubation at room temperature was required. A polymer-based IHC detection system was used as the amplification system (EnVision + System-HRP Labelled Polymer Anti-Rabbit or Anti-Mouse secondary antibody, Dako, Carpinteria, CA, USA) according to the manufacture's instruction. The bound antibody was visualized by incubating the sections in 3-amino-9-ethylcarbazole (AEC + Substrate Chromogen, Dako, Carpinteria, CA, USA) for 15 min, except for anti-CD34, which requires the use of 3,3′-diaminobenzidine (Liquid DAB + Substrate Chromogen System, Dako, Carpinteria, CA, USA) for 8-10 min. Cell nuclei were counterstained with Mayer's Haematoxylin (Bio-Optica, MI, Italy) and the slides were mounted with aqueous mounting medium (Faramount Aqueous Mounting Medium, Dako, Carpinteria, CA, USA). The different analysed biomarkers, dilution, source/clone, the staining localization of antibody and the cut off [median value or immunohistochemical score (HIS)] used to classify positive versus negative cases are shown in Additional file 1: Table S1.
All immunostained samples were scored by doubleblinded independent observers who had no patient outcome and clinicopathological data information, and the mean of the three readings for each patient was calculated. If one core was uninformative, lost or contained no tumor tissue, the overall score applied was that of the remaining cores. The results from two observers were identical in most cases, and discrepancies were resolved by re-examination and consensus.
NHERF1 immunostaining was predominantly cytoplasmic (cNHERF1), however in some cases an intense nuclear (nNHERF1) and membranous (mNHERF1) staining were also demonstrated. These were scored separately and their significance was evaluated. VEGF and VEGFR1 were mainly observed in the cytoplasm of breast cancer tissues (cVEGF and cVEGFR1), and the staining was evaluated as HIS or percentage of immunoreactive cells, respectively. Nuclear localization of HIF-1α (nHIF-1α) and TWIST1 (nTWIST1) was mainly observed. Cytoplasmic staining of HIF-1α and TWIST1 was occasionally reported even if it was not evaluated for the aims of the study. The median value of immunoreactive cells was used as cut off for cNHERF1 (≥40%), nNHERF1 (> 0%), mNHERF1 (> 0%), nTWIST1 (≥4%), cVEGFR1 (> 0%), nHIF-1α (> 0%) and MVD (≥15 microvessels/mm 2 ). For cVEGF, the HIS was calculated by combining the quantity score (percentage of positive stained cells) with the staining intensity score [6]. The quantity score ranged from 0 to 4: 0 = no immunoreactivity; 1 ≤ 25% cells stained; 2 = 26-50% cells stained; 3 = 51-75% cells stained; and 4 = ≥76% cells stained. The staining intensity was scored as: 0 (negative), 1 (weak), 2 (moderate) and 3 (strong). Raw data were converted to IHS by adding the quantity score (0-4) to the staining intensity score (0-3). Theoretically, the scores can range from 0 to 7. An IHS of 6-7 was considered a strong immunoreactivity; 3-5, moderate; 1-2, weak; and 0, negative. For our analyses, tumors presenting a moderate or strong score were cVEGF positive (IHS:3-7). Finally, microvessel counting was performed by identifying the areas which represented the highest vascular density, the so called "hot spots". The MVD measurements were made in the fields with a higher density of CD-34 positive cells and cell clusters at 200× magnification, as previously described [6]. Positive and negative controls were included in each staining run as indicated in the data sheet of each antibody. The accuracy, reliability, scoring strategy and reproducibility assessments of these antibodies have been validated in previous studies already published [6,[24][25][26].
Immunofluorescence Immunofluorescent analysis was performed as described previously [18]. Briefly, formalin-fixed and paraffin embedded tissue serial sections of 3 μm in thickness were deparaffinized with xylene, and rehydrated in an ethanol series. Antigen retrieval was carried out immersing slide in a 0.01 M saline citrated buffer (pH 6.0) at 95°C for 40 min, then tissues were permeabilized with 0.1% Triton X100-Phosphate Buffered Saline for 15 min, blocked 30 min with 1% Bovine Serum Albumin-Phosphate Buffered Saline and incubated overnight at 4°C in a humidified chamber with a mouse monoclonal anti-CD34 (clone QBEND-10, dilution 1:20; Novocastra Lab. Ltd., UK) together with a rabbit polyclonal anti-NHERF1 (PA1-090, 1 μg/100 μL dilution; Affinity Bio-Reagents). The slides were then incubated at room temperature for 1 h with the Alexa Fluor 488 and Alexa Fluor 568 immunoglobulin G secondary conjugated antibodies (1: 2000 dilution; Molecular Probes Inc., Eugene, OR, USA) and mounted with DAPI (ProLong® Gold antifade reagent; Molecular Probes Inc.). Positive control slides that were run simultaneously were used for assessing the quality of immunoreactivity. For negative controls, slide sections that were immunopositive were treated with 1% Bovine Serum Albumin instead of the primary antibody, and no reactivity was observed in any of these controls. Images were obtained on a BX40 microscope (Olympus, Tokyo, Japan) with a SenSys 1401E-Photometrics chargecoupled device camera. To verify protein colocalization, each acquired stack was merged by transforming the three channels corresponding to red (tetramethylrhodamine B isothiocyanate), green (fluorescein isothiocyanate) and blue (4′,6-diamidino-2-phenylindole) into a single three-color stack by using the "RGB merge" command of ImageJ software (National Institutes of Health Bethesda, MD).

Follow up and statistical analysis
Associations of tumor markers expression with various clinicopathological features and the among markers themselves were determined by the Chi-squared test for the independence of categorical variables. The presence of non-linear relationship between two continuous variables was assessed by the Kendall rank test.
The results from the immunohistochemical analyses of our biomarkers were assessed in relation to disease-free survival (DFS) and overall survival (OS). DFS (in months) was defined as the time from diagnosis to the date of loco-regional or distant recurrence, second invasive breast carcinoma, second primary cancer and/or death without evidence of breast cancer or to the date of last contact. OS (in months) was defined as the time from diagnosis to the date of last contact or of death from any cause. Univariate analyses of DFS and OS were carried out considering both tumor markers expression and clinicopathological variables: five-year DFS and OS and their 95% confidence interval were based on Kaplan-Meier method, hazardratio (HR) were computed using the Cox proportional hazard regression model for DFS only (for OS it could not be computed due to the low number of deaths), and p-values of the log-rank test were used to infer on the equality of probability of an event (relapse for DFS or death for OS) in the two groups. Multivariate analyses on DFS (OS could not be analyzed due to the low number of deaths) were performed for tumor markers expressed as categorical and continuous variables using the Cox proportional hazard regression model; HRs were computed along with their p-values.
Data analysis was carried out using the statistical packages survival and survminer of the statistical language R, version 3.4 [27].

Protein expression profiling
The immunohistochemical expression of NHERF1, TWIST1, VEGFR1, VEGF, HIF-1α and MVD, was evaluated on TMA sections. The frequency of expression of these biomarkers, in the whole cohort of tumor samples is reported in Additional file 2: Table S2. NHERF1 expression was detected in the apical membrane, cytoplasm, and nucleus of tumor samples. These different localizations were scored separately and their significance was evaluated. All breast cancers showed NHERF1 protein localized in the cytoplasm of tumor cells, of these 49% (76/154) overexpressed cNHERF1, according to the median value ≥40%. Membranous NHERF1 and nNHERF1expression were detected in 13% (20/151) and 20% (30/149) of cases, respectively (median value > 0 for both compartments). Cytoplasmic VEGF was expressed in 96% (144/150) of tumor cells, of these 54% (81/151) overexpressed it, according to the IHS median value ≥3. Cytoplasmic VEGFR1 and nuclear HIF1α were expressed in about 50% (80/161) and 37% (56/153) of tumor cells and, according to the median value > 0 for both markers. Nuclear TWIST1 was expressed in 60% (84/140) of tumor cells, of these 51% (72/140) overexpressed it, according to the median value ≥4. Tumors showed microvessel presence and 54% (85/157) of them were characterized by a higher MVD, according to the median value ≥15 microvessels/mm 2 . Figure 1 shows an example of immunohistochemical staining pattern of NHERF1 and VEGF, VEGFR1, HIF-1α, TWIST1 and MVD expression in TMA tumor cores. NHERF1 was present in cytoplasmic compartment and less frequently at nuclear level, according with the median values of expression.
VEGF, VEGFR1, HIF-1α, TWIST1 and MVD expression was heterogeneous for intensity and expression patterns in the different tumor samples examined.

Relationship between tumor markers and clinicopathological features
A summary of significant association between tumor marker expressions and clinicopathological features are reported in Table 2. Among the BCs with cNHERF1 overexpression, 67.1% exhibited a significant association with tumor size > 2 cm (p = 0.023) and 68.4% with MIB1 positive expression (p = 0.027). Further, negative cNHERF1 expression was related to PgR-positive expression (p = 0.023) and HER2/neu negative (p = 0.002) in 67.1% and 65.2% of cases respectively. The lack of nuclear NHERF1 was noticeably associated to tumor size > 2 cm (p = 0.014), while its presence was linked to MIB1 negative expression (p < 0.001). Analyzing NHERF1 reactivity at membrane level, its absence was observed in IDCs (p = 0.002) and positive lymph node status (p = 0.013); conversely mNHERF1 expression was significantly associated with tumor size ≤2 cm (p = 0. 039) and MIB1 negative (p = 0.019). A statistically significant relation was observed between mNHERF1 and PgRpositive expression (p = 0.034).
High cVEGF expression was present in 61.7% of younger patients (p = 0.006), and it showed a significant association with positive lymph node status (p = 0.010), higher tumor histological grade (p = 0.020), MIB1 positive expression (p < 0.001). CytoplasmicVEGF negative detection was related to HER2/neu negative expression (p = 0.006). Furthermore, a statistically significant relation was  observed between cVEGF negative and ER-positive expression (p = 0.006). Expression of cVEGFR1 was statistically prevalent in invasive tumor phenotype (p = 0.041) and PgR-positive expression (p = 0.009). G2 tumors presented high nHIF-1α expression (p = 0.001), and this overexpression was directly associated with PgR-positive expression (p = 0.037) while nHIF-1α negative expression was associated to MIB1 positive expression (p = 0.016). Low nTWIST1 expression was statistically prevalent in high MIB1 tumors (p = 0.015). There were no statistically significant associations between MVD and the clinicopathological features.
No statistically significant association was found between MVD and NHERF1 for both dichotomized and continuous variables.
The analysis of the localization of VEGF and its receptor VEGFR1, by means of immunofluorescence studies, indicated that they co-localized with NHERF1, when both proteins were over-expressed within cytoplasmic and cytoplasmic and/or membranous compartments respectively, in invasive cellular clusters. Furthermore, our studies indicated that nHIF-1α and nTWIST1 co-localized with NHERF1 too (Fig. 3a). It was interesting to observe that while an immunofluorescence assay showed a high MVD, a neo-vascular formation with a mixed labeling CD34/ NHERF1 was also noted (Fig. 3b).

Expression of the proteins and patient outcome
Univariate analyses were carried out for all the clinicopathological characteristics and the expression of cNHERF1, nNHERF1, mNHERF1, cVEGF, cVEGFR1, MVD, nHIF-1α and nTWIST1 proteins, as dichotomized variables. These were correlated to disease-free survival (DFS) and overallsurvival (OS) ( Table 4). The subgroup of patients with negative expression of nTWIST1 had a shorter 5-year DFS (p < 0.001); the same held true for continuous data (p = 0.022, data not shown). Moreover, univariate analysis of clinicopathological characteristics in the entire cohort revealed that only MIB1 was significantly associated with worse DFS (p = 0.003).
According to the Cox proportional hazard regression model, multivariate analysis of the entire cohort (Table 5), identified nTWIST1 expression, in both categorical and continuous data (HR = 0,14, 95% confidence interval: CI 0.05-0.43, p = 0.001; and HR = 0.92, 95% CI 0.87-0.97, Fig. 3 Immunofluorescence analysis of NHERF1 and TME biomarkers expressions in invasive breast cancer. a A representative tissue sample stained with NHERF1 and VEGF, VEGFR1, HIF-1α and TWIST1 antibodies and detected with Alexa Fluor 568 (red) and Alexa Fluor 488 (green) secondary antibodies and nucleus in blue (DAPI stained), respectively, prior to fluorescence microscopy analysis (200× magnification). Arrowheads indicate co-localization sites of NHERF1 with the others markers. b A representative tissue sample stained with NHERF1 and CD34 antibodies and detected with Alexa Fluor 568 (red) and Alexa Fluor 488 (green) secondary antibodies and nucleus in blue (DAPI stained), respectively, prior to fluorescence microscopy analysis. The square includes microvessel formation with a double staining NHERF1-CD34 The relationship among the different analyzed biomarkers expression and BCs survival was then investigated. Kaplan-Meier curves revealed that patients with nTWIST1expression had a worse DFS than of patients with nTWIST1+ expression (5-years 77% vs 92%, p < 0.001) (Fig. 4a).

Discussion
The study regarding cancer context has become of great significance because of the therapeutic implications linked to it. In recent years, the complexity of TME has been evidenced enough for it to be considered as a real system in an ever-changing state.
In this study, a panel of factors involved in TME changes such as VEGF, VEGFR1, HIF-1a, TWIST1 and MVD, that have been related to BC progression, aggressiveness and invasion was analyzed [6,12,[28][29][30]. Further, so as to better characterize this microenvironment, their association with the scaffolding protein NHERF1 and with survival was observed, allowing to identify a possible subgroup of patients that could benefit from a specific therapy.
Mainly, in this BC cohort there was evidence of an increase of cNHERF1 expression with a loss of mNHERF1 detection and a small nNHERF1 presence, in conformity with a more aggressive behavior. Our group had already shown that the switch from apical membranous to cytoplasmic expression was linked to breast carcinogenesis, with a high cNHERF1 expression, present in more aggressive tumors [23]. In addition, the loss of nNHERF1 expression had also been associated with reduced survival, indicating it as an independent prognostic marker in BC [25]. An increase of TME analyzed markers had  been reported in these tumors, to confirm their aggressiveness. Moreover, it was shown that cNHERF1 overexpression was related to some aggressive clinical parameters, as demonstrated by association with extensive tumor size and high proliferative activity. Membranous and nuclear NHERF1 expression reaffirmed their protective role, being related to low proliferative activity and smaller tumors. Moreover, our results supported that cVEGF might have a potential additional predictive value, as demonstrated in this population, including younger patients with larger tumors, extensive lymph-node metastases and high grade and proliferative activity. Its receptor cVEGFR1 over-expression was associated to an invasive phenotype as well. VEGF and VEGFR1 have been indicated as strong biomarker indicators to select patients who could benefit from anti-angiogenic therapy [31]. Despite this is an innovative approach, therapeutical failures projected for novel solutions to improve and overcome actual limits [32]. It was e also demonstrated, in this study, that cNHERF1 overexpression was positively related with cVEGF expression. Interestingly, a significant direct association between increased levels of cNHERF1 and cVEGFR1 was also noticed. At the same time, membranous and nuclear NHERF1 reconfirmed their defensive role, being negatively related to cVEGF expression. VEGF is well known as the principal modulator in the process of cancer growth, invasion and metastasis, through new microvessel formation and lymphovascular invasion [33]. Previous studies carried out by our équipe had already reported a NHERF1/VEGFR1 relationship with poor clinical outcome in invasive BCs, indicating a possible interaction between the two proteins [19]. VEGF and VEGFR1 high expression is part of highly hypoxic and aggressive environments, involving high HIF-1α and TWIST1 levels of expression too. Several studies have identified a relation to poorer prognosis/therapy response in the altered expression of HIF-1α and TWIST1 [34][35][36]. Moreover, the hypoxic environment in BCs is one of the key regulators of the network of the biological elements of matrix, involved in metastasis process and it affects both the early and late stages of metastasis [37]. The dissemination process induced by hypoxia is supplied by paracrine and autocrine stimuli and it is interconnected with EMT program and characterized by high expression of different proteins [38,39]. The results in this study indicated a direct relation between nHIF-1α and nTWIST1 expression, but a controversial action was present, being both nHIF-1α and nTWIST1 loss interrelated to more aggressive phenotype. Furthermore, an inverse relation was found between cNHERF1 and nHIF-1α and cNHERF1 and nTWIST1 expression, even if a trend consistent with cNHERF1 positive expression was expected, being it correlated to BCs more aggressive phenotype [18,23,25]. Additionally, atypically, high nHIF-1α and nTWIST1 expression was present in tumors with low cVEGF expression. Different factors could have contributed to these results, especially knowing that TME is characterized by continuous changing. A possible reason for these discrepancies could be the heterogeneity of BC. Cancer is inherently a heterogeneous disease and so heterogeneity has become its main feature. Clinical and histo-pathological parameters, biomarkers and genetic heterogeneity can be affected by the different characteristics of a single patient (age, menopausal state, general health status etc), but also by epigenetic changes (histone modification, DNA methylation etc) [40,41] and metabolic reprogramming [42]. Expression of biomarkers can be highly variable within an individual tumor, causing interpretation problems and discordant results. Regarding nTWIST1 negative expression, it was observed that most of these tumors presented more aggressive clinic-pathological characteristics (high proliferative activity, positive lymph node status, large tumor size, high histological grade and HER2-postive status), despite the fact that there was not a significant statistical evidence for any of these.
Evidence of all this scenario was the presence of new vessels formation, detected by immunohistochemistry and immunofluorescence assay. MVD showed only a direct association with cVEGFR1 over-expression, while no other significant relation has been found. Remarkably, in the tumor area, a neo-vascular formation with a "mosaic" structure was observed with a double staining of CD34 and NHERF1, which could indicate a probable tumor cell recruitment in neo blood vessels formation with endothelial cells. This is a well-known mechanism, which in addition to sprouting angiogenesis, has been recognized to contribute to tumor vascularization. In fact, new blood vessels can be due to vasculogenic mimicry formation, that comprise not only vessels formed by endothelial cells, but also tubules formed by malignant cancer cells [43][44][45]. NHERF1 involvement in neovascularization phenomena has already been reported [46,47], but this is the first evidence of a possible direct participation of tumor cells NHERF1 positive to new-vessels formation.
In a subgroup of patients with long-term follow-up, multivariate analysis showed that loss nTWIST1 expression was related to a decrease of DFS. Other authors had reported a worse DFS linked to TWIST1 negative expression, suggesting that its deregulation was involved in poor patient outcome [48,49]. The Kaplan Meyer curves were used to evaluate the relationship between the biomarkers. Interestingly, when nTWIST1-/mNHERF1+ phenotype were considered, it was discovered that these patients showed a better DFS with respect to nTWIST1 +/mNHERF1-phenotype. In addition, data showed a worse outcome in patients nTWIST1-/mNHERF1-vs nTWIST1+/mNHERF1+, despite it being not statistically significant. These findings confirmed the oncosuppressor role of NHERF1 expression, when it was localized at the plasma membrane [50][51][52], suggesting that mNHERF1 expression might interfere and influence clinical outcome.

Conclusions
In conclusion, the study reported a closer interaction of cNHERF1 with cVEGF, cVEGFR1 and nTWIST1, that denoted an active participation of tumor microenvironment to cancer aggressiveness. From this point of view, it is very interesting to consider the mimicry involving NHERF1 + tumor cells in neo-microvessels formation, and this opens a new stimulating scenario for further NHERF1 studies. Moreover, the identification of nTWIST1-/mNHERF1+ phenotype as a subgroup with an increased DFS and nTWIST1+/cNHERF1+, nTWIST1-/cNHERF1+ phenotypes as patient subgroups with poor outcome, reinforced both the tumor suppressor role of mNHERF1 and oncogenic activity of cNHERF1. Furthermore, these date suggest the potential prognostic role of NHERF1. Finally, we believe that our intriguingly results warrant additional studies aimed to evaluate tumor microenvironment regulation in combination with NHERF1 as a possible targeted-oriented therapeutic approach.