- Open Access
NHERF1 and tumor microenvironment: a new scene in invasive breast carcinoma
Journal of Experimental & Clinical Cancer Researchvolume 37, Article number: 96 (2018)
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.
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.
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).
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.
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 . 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 . 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 . 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 . TWIST1 is a basic helix-loop-helix transcription factor that contributes to carcinogenesis by triggering EMT, thereby influencing tumor invasion behavior . 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 . 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 . Furthermore, at membranous and cytoplasmic level, its overexpression and co-localization with VEGFR1, a marker correlated with high metastasis risk and relapse has been evaluated . 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 . In a recent study, NHERF1 expression was positively correlated with VEGFR2 expression in CRC, suggesting NHERF1 involvement in disease progression by VEGFR2 signaling pathway . 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 .
TMAs were constructed manually from formalin-fixed and paraffin-embedded tissues, according to standard procedures as previously described . 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 . 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% H2O2. 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 double-blinded 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/mm2). For cVEGF, the HIS was calculated by combining the quantity score (percentage of positive stained cells) with the staining intensity score . 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 . 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].
Immunofluorescent analysis was performed as described previously . 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 charge-coupled 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, hazard-ratio (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 .
The tumor characteristics of the 183 patients included in the study are listed in Table 1. The 48.6% of patients were younger than 49 years. The majority had an invasive ductal carcinoma (IDC) (90.7%), a moderate histological grade (G2) (45.9%), and tumors larger than 2 cm (55.2%) with axillary lymph node involvement (59.2%). Furthermore, most patients had ER (71%) and PgR (58.5%) positive status, high MIB1 (57.9%) and HER2/neu negative (50.6%); 10.4% were triple negative breast cancers. The treatment of only 106 patients was known. Of these, 25.5% received adjuvant chemotherapy, 15.1% of ER-positive patients received hormonotherapy, and 39.6% received adjuvant tamoxifen for 5 years after the end of chemotherapy.
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/mm2. 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 PgR-positive 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 over-expression 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.
Association between protein expressions analyzed
As regards the dichotomized variables, Chi-squared test for protein associations is summarized in Table 3. For the continuous variables, Kendall’s τ for rank-based correlations between markers and relative heatmap are shown in Fig. 2a and b.
Of 183 total tumors, 138 could be assessed for both cVEGF and cNHERF1. In particular, among these 79 tumors over-expressed cNHERF1 and 65.8% had high cVEGF expression as well; additionally, of 59 cases with negative cNHERF1, 66.1% had negative cVEGF expression (p < 0.001). Cytoplasmic VEGF negative expression was also related to nuclear (p = 0.001) and membranous NHERF1 expression (p = 0.001) in 91.1% and 94.4% of cases, respectively. Kendall’s τ for rank-based correlation underlined direct relationship between cVEGF and cNHERF1 (τ = 0.323, p < 0.001) and inverse between nNHERF1 (τ = − 0.247, p < 0.001) and mNHERF1 (τ = − 0.178, p = 0.012). One hundred and forty-seven tumors resulted evaluable for both cVEGFR1 and cNHERF1 expression. In particular, 53.7% the tumors showed cNHERF1 over-expression and 59.5% had also high cVEGFR1 expression (p = 0.027). No statistically significant correlation was observed between cVEGFR1 and nuclear and membranous NHERF1 expression. Kendall’s τ test for continuous data showed an inverse association between cVEGFR1 and mNHERF1 (τ = − 0.161, p = 0.024), while confirming the direct association between cVEGFR1 and cNHERF1 expression (τ = 0.189, p = 0.003). Nuclear HIF-1α and cNHERF1 expression was observed in 142 tumors. Low nHIF-1α expression was related to 60.4% of the cNHERF1 positive tumors (p = 0.021) while low nuclear and membranous NHERF1 was related to negative nHIF-1α expression (p = 0.001 and p = 0.046) for 87.6% and 91% of cases, respectively. Kendall’s τ test confirmed the presence of an inverse correlation between nHIF-1α and cNHERF1 (τ = − 0.191, p = 0.004) and a positive correlation between nHIF-1α and nNHERF1 (τ = 0.260, p < 0.001). Nuclear TWIST1 and cNHERF1 expression were detected in 130 tumors. Nuclear TWIST1 over-expression linked to a low cNHERF1 detection (p = 0.001) appeared in 71% of the cases while low expression of nNHERF1 and mNHERF1 was related to low nTWIST1 expression (p = 0.008 and p = 0.035) in 90.3% and 93.5% of tumors respectively. The data were also supported by Kendall’s τ test, which emphasized a reverse behavior of nTWIST1 and cNHERF1 (τ = − 0.200, p = 0.003) and a positive relation between nTWIST1 and nuclear (τ = 0.205, p = 0.005) and mNHERF1 (τ = 0.265, p < 0.001). 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).
Moreover, of the 183 total tumors, 145 could be assessed for both cVEGF and cVEGFR1, and 64.9% showed a positive relation between cVEGF and cVEGFR1 expression (p = 0.019). Of 141 tumor samples valuable for both cVEGF and nHIF-1α, 62.2% of the tumor samples presented cVEGF negative expression and nHIF-1α positive expression (p = 0.049). In addition, of the 128 tumors that could be assessed for both cVEGF and nTWIST1, 67.2% of the samples presented cVEGF negative expression and nTWIST1 positive expression (p = 0.013). High MVD was associated to cVEGFR1 positive expression (p = 0.013) and nHIF-1α/nTWIST1 negative expressions were together related (p = 0.001) in 61.7% and 78.5% of the tumors, respectively (Table 3). Furthermore, Kendall’s τ test also revealed a positive correlation between cVEGF and cVEGFR1 expression (τ = 0.236, p < 0.001), between cVEGFR1 expression and MVD (τ = 0.137, p = 0.024), and nHIF-1α and nTWIST1 expression (τ = 0.236, p = 0.001). While, an inverse relation between cVEGF and nHIF-1α (τ = − 0.212, p = 0.001), cVEGF and nTWIST1 expression (τ = − 0.213, p = 0.001),was observed (Fig. 2b).
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 overall-survival (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, p = 0.003, respectively) as significant for DFS. Moreover, for continuous data (HR = 1.09, 95%CI 1.02–1.17, p = 0.011), mNHERF1 expression was a significant indicator of DFS (Table 5) as well. HR for OS was not computed due to the low number of events.
The relationship among the different analyzed biomarkers expression and BCs survival was then investigated. Kaplan-Meier curves revealed that patients with nTWIST1- expression had a worse DFS than of patients with nTWIST1+ expression (5-years 77% vs 92%, p < 0.001) (Fig. 4a).
When nTWIST1 and mNHERF1 expressions were considered, Kaplan-Meier curves showed that patients with nTWIST1-/mNHERF1+ immunophenotype had a better DFS respect to patients with nTWIST1+/mNHERF1- (p < 0.001) (Fig. 4b). Analyzing patients with nTWIST1-/mNHERF1- vs nTWIST1+/mNHERF1+ expression, no statistically significant results were found (p = 0.510), (Fig. 4c).
In addition, the subgroup of patients nTWIST1-/cNHERF1- showed a higher DFS with respect to the subgroup nTWIST1+/cNHERF1+ (p = 0.004). When nTWIST1-/cNHERF1+ vs nTWIST1+/cNHERF1- phenotype were compared, the first group showed a worse DFS (p = 0.037) (Fig. 4d, e). Kaplan-Meier curves, computed for the other proteins analyzed, revealed that there were no statistically significant results comparing the DFS of the group of patients with cVEGF-/cNHERF1- vs cVEGF+/cNHERF1+ expression (p = 0.430) (Fig. 4f), and with cVEGFR1-/cNHERF1- vs cVEGFR1+/cNHERF1+ expression (p = 0.180) (Fig. 4g) and nHIF-1α-/cNHERF1+ vs nHIF-1α+/cNHERF1- expression (p = 0.950) (Fig. 4h). There were no significant correlations for protein levels and OS.
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 . In addition, the loss of nNHERF1 expression had also been associated with reduced survival, indicating it as an independent prognostic marker in BC . An increase of TME analyzed markers had been reported in these tumors, to confirm their aggressiveness. Moreover, it was shown that cNHERF1 over-expression 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 . Despite this is an innovative approach, therapeutical failures projected for novel solutions to improve and overcome actual limits . 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 . 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 . 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 . 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 . 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.
Thus, the analyses of nTWIST1+/cNHERF1+ patients selected a phenotype with a worse DFS compared to nTWIST1-/cNHERF1-, as well as the nTWIST1-/cNHERF1+ with respect to nTWIST1+/cNHERF1-, corroborating the pivotal role of cNHERF1 in cancer progression, repeatedly stated by previous researches [18, 23,24,25].
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.
epithelial to mesenchimal transition
human epidermal growth factor receptor 2
hypoxia inducible factor-1 alpha
Na+/H+ exchanger regulatory factor 1
vascular endothelial growth factor
vascular endothelial growth factor receptor 1
Quail DF, Joyce JA. Microenvironmental regulation of tumor progression and metastasis. Nat Med. 2013;19:1423–37.
Spill F, Reynolds DS, Kamm RD, Zaman MH. Impact of the physical microenvironment on tumor progression and metastasis. Curr Opin Biotechnol. 2016;40:41–8.
Ghoncheh M, Pournamdar Z, Salehiniya H. Incidence and Mortality and Epidemiology of Breast Cancer in the World. Asian Pac J Cancer Prev. 2016;17:43–6.
Alizadeh AA, Aranda V, Bardelli A, Blanpain C, Bock C, Borowski C, et al. Toward understanding and exploiting tumor heterogeneity. Nat Med. 2015;21:846–53.
Andreuzzi E, Colladel R, Pellicani R, Tarticchio G, Cannizzaro R, Spessotto P, et al. The angiostatic molecule Multimerin 2 is processed by MMP-9 to allow sprouting angiogenesis. Matrix Biol. 2017;64:40–53.
Saponaro C, Malfettone A, Ranieri G, Danza K, Simone G, Paradiso A, et al. VEGF, HIF-1α expression and MVD as an angiogenic network in familial breast cancer. PLoS One. 2013;8:e53070.
De Francesco EM, Pellegrino M, Santolla MF, Lappano R, Ricchio E, Abonante S, et al. GPER mediates activation of HIF1α/VEGF signaling by estrogens. Cancer Res. 2014;74:4053–64.
Kaufhold S, Bonavida B. Central role of Snail1 in the regulation of EMT and resistance in cancer: a target for therapeutic intervention. J Exp Clin Cancer Res. 2014;33:62.
Zhu QQ, Ma C, Wang Q, Song Y, Lv T. The role of TWIST1 in epithelial-mesenchymal transition and cancers. Tumour Biol. 2016;37:185–97.
Sun X, Cui M, Zhang A, Tong L, Wang K, Li K, et al. MiR-548c impairs migration and invasion of endometrial and ovarian cancer cells via downregulation of Twist. J Exp Clin Cancer Res. 2016;35:10.
Vesuna F, Bergman Y, Raman V. Genomic pathways modulated by Twist in breast cancer. BMC Cancer. 2017;17:52.
Schirosi L, De Summa S, Tommasi S, Paradiso A, Gasparini G, Popescu O, et al. VEGF and TWIST1 in a 16-biomarker immunoprofile useful for prognosis of breast cancer patients. Int J Cancer. 2017;141:1901–11.
Sossey-Alaoui K, Pluskota E, Davuluri G, Bialkowska K, Das M, Szpak D, et al. Kindlin-3 enhances breast cancer progression and metastasis by activating Twist-mediated angiogenesis. FASEB J. 2014;28:2260–71.
Dunn HA, Ferguson SS. PDZ Protein Regulation of G Protein-Coupled Receptor Trafficking and Signaling Pathways. Mol Pharmacol. 2015;88:624–39.
Mangia A, Scarpi E, Partipilo G, Schirosi L, Opinto G, Giotta F, et al. NHERF1 together with PARP1 and BRCA1 expression as a new potential biomarker to stratify breast cancer patients. Oncotarget. 2017;8:65730–42.
Georgescu MM, Mobley BC, Orr BA, Shang P, Lehman NL, Zhu X, et al. NHERF1/EBP50 and NF2 as diagnostic markers for choroid plexus tumors. Acta Neuropathol Commun. 2016;4:55.
Mangia A, Caldarola L, Dell'Endice S, Scarpi E, Saragoni L, Monti M, et al. The potential predictive role of nuclear NHERF1 expression in advanced gastric cancer patients treated with epirubicin/oxaliplatin/capecitabine first line chemotherapy. Cancer Biol Ther. 2015;16:1140–7.
Cardone RA, Bellizzi A, Busco G, Weinman EJ, Dell'Aquila ME, Casavola V, et al. The NHERF1 PDZ2 domain regulates PKA-RhoA-p38-mediated NHE1 activation and invasion in breast tumor cells. Mol Biol Cell. 2007;18:1768–80.
Malfettone A, Saponaro C, Paradiso A, Simone G, Mangia A. Peritumoral vascular invasion and NHERF1 expression define an immunophenotype of grade 2 invasive breast cancer associated with poor prognosis. BMC Cancer. 2012;12:106.
Malfettone A, Silvestris N, Paradiso A, Mattioli E, Simone G, Mangia A. Overexpression of nuclear NHERF1 in advanced colorectal cancer: association with hypoxic microenvironment and tumor invasive phenotype. Exp Mol Pathol. 2012;92:296–303.
Gu Y, Yu H, Hao C, Martin TA, Hargest R, He J, et al. NHERF1 regulates the progression of colorectal cancer through the interplay with VEGFR2 pathway. Oncotarget. 2017;8:7753–65.
Wolff AC, Hammond ME, Schwartz JN, Hagerty KL, Allred DC, Cote RJ, et al. American Society of Clinical Oncology/College of American Pathologists guideline recommendations for human epidermal growth factor receptor 2 testing in breast cancer. J Clin Oncol. 2007;25:118–45.
Mangia A, Chiriatti A, Bellizzi A, Malfettone A, Stea B, Zito FA, et al. Biological role of NHERF1 protein expression in breast cancer. Histopathology. 2009;55:600–8.
Schirosi L, De Summa S, Tommasi S, Paradiso A, Sambiasi D, Popescu O, et al. Immunoprofile from tissue microarrays to stratify familial breast cancer patients. Oncotarget. 2015;6:27865–79.
Paradiso A, Scarpi E, Malfettone A, Addati T, Giotta F, Simone G, et al. Nuclear NHERF1 expression as a prognostic marker in breast cancer. Cell Death Dis. 2013;4:e904.
Mangia A, Malfettone A, Saponaro C, Tommasi S, Simone G, Paradiso A. Human epidermal growth factor receptor 2, Na+/H+ exchanger regulatory factor 1, and breast cancer susceptibility gene-1 as new biomarkers for familial breast cancers. Hum Pathol. 2011;42:1589–95.
R Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2017. https://www.R-project.org/.
Zhang L, Wang H, Li C, Zhao Y, Wu L, Du X, et al. VEGF-A/Neuropilin 1 pathway confers Cancer Stemness via activating Wnt/β-catenin Axis in breast Cancer cells. Cell Physiol Biochem. 2017;44:1251–62.
Weddell JC, Chen S, Imoukhuede PI. EGFR1 promotes cell migration and proliferation through PLCγ and PI3K pathways. NPJ Syst Biol Appl. 2017;4:1.
Liang G, Liu Z, Tan L, Su AN, Jiang WG, Gong C. HIF1α-associated circDENND4C promotes proliferation of breast Cancer cells in hypoxic environment. Anticancer Res. 2017;37:4337–43.
Lambrechts D, Lenz HJ, de Haas S, Carmeliet P, Scherer SJ. Markers of response for the antiangiogenic agent bevacizumab. J Clin Oncol. 2013;31:1219–30.
Xu Y, Li Q, Li XY, Yang QY, Xu WW, Liu GL. Short-term anti-vascular endothelial growth factor treatment elicits vasculogenic mimicry formation of tumors to accelerate metastasis. J Exp Clin Cancer Res. 2012;31:16.
Sa-Nguanraksa D, Chuangsuwanich T, Pongpruttipan T, O-Charoenrat P. High vascular endothelial growth factor gene expression predicts poor outcome in patients with non-luminal A breast cancer. Mol Clin Oncol. 2015;3:1103–8.
Jin MS, Lee H, Park IA, Chung YR, Im SA, Lee KH, et al. Overexpression of HIF1α and CAXI predicts poor outcome in early-stage triple negative breast cancer. Virchows Arch. 2016;469:183–90.
Semenza GL. Targeting HIF-1 for cancer therapy. Nat Rev Cancer. 2003;3:721–32.
Grzegrzolka J, Wojtyra P, Biala M, Piotrowska A, Gomulkiewicz A, Rys J, et al. Correlation between expression of Twist and Podoplanin in ductal breast carcinoma. Anticancer Res. 2017;37:5485–93.
Rankin EB, Giaccia AJ. Hypoxic control of metastasis. Science. 2016;352:175–80.
Ye X, Weinberg RA. Epithelial-mesenchymal plasticity: a central regulator of cancer progression. Trends Cell Biol. 2015;25:675–86.
Tiwari N, Gheldof A, Tatari M, Christofori G. EMT as the ultimate survival mechanism of cancer cells. Semin Cancer Biol. 2012;22:194–207.
Beca F, Polyak K. Intratumor Heterogeneity in Breast Cancer. Adv Exp Med Biol. 2016;882:169–89.
Turashvili G, Brogi E. Tumor Heterogeneity in Breast Cancer. Front Med (Lausanne). 2017;4:227.
Yoshida GJ. Metabolic reprogramming: the emerging concept and associated therapeutic strategies. J Exp Clin Cancer Res. 2015;34:111.
Zhang S, Zhang D, Sun B. Vasculogenic mimicry: current status and future prospects. Cancer Lett. 2007;254:157–64.
Zang M, Zhang Y, Zhang B, Hu L, Li J, Fan Z, et al. CEACAM6 promotes tumor angiogenesis and vasculogenic mimicry in gastric cancer via FAK signaling. Biochim Biophys Acta. 2015;1852:1020–8.
Kaessmeyer S, Bhoola K, Baltic S, Thompson P, Plendl J. Lung cancer neovascularisation: Cellular and molecular interaction between endothelial and lung cancer cells. Immunobiology. 2014;219:308–14.
Hou Y, Wu Y, Farooq SM, Guan X, Wang S, Liu Y, et al. A critical role of CXCR2 PDZ-mediated interactions in endothelial progenitor cell homing and angiogenesis. Stem Cell Res. 2015;14:133–43.
Chen P, Wang Y, Yang L, Li C, Wang Y, Xie L, et al. Novel bioactivity of NHERF1 in corneal neovascularization. Graefes Arch Clin Exp Ophthalmol. 2012;250:1615–25.
Montserrat N, Gallardo A, Escuin D, Catasus L, Prat J, Gutiérrez-Avignó FJ, et al. Repression of E-cadherin by SNAIL, ZEB1, and TWIST in invasive ductal carcinomas of the breast: a cooperative effort? Hum Pathol. 2011;42:103–10.
Cates JM, Byrd RH, Fohn LE, Tatsas AD, Washington MK, Black CC. Epithelial-mesenchymal transition markers in pancreatic ductal adenocarcinoma. Pancreas. 2009;38:e1–6.
Molina JR, Agarwal NK, Morales FC, Hayashi Y, Aldape KD, Cote G, et al. PTEN, NHERF1 and PHLPP form a tumor suppressor network that is disabled in glioblastoma. Oncogene. 2012;31:1264–74.
Pan Y, Wang L, Dai JL. Suppression of breast cancer cell growth by Na+/H+ exchanger regulatory factor 1 (NHERF1). Breast Cancer Res. 2006;8:R63.
Kreimann EL, Morales FC, de Orbeta-Cruz J, Takahashi Y, Adams H, Liu TJ, et al. Cortical stabilization of beta-catenin contributes to NHERF1/ EBP50 tumor suppressor function. Oncogene. 2007;26:5290–9.
The authors would like to thank Francesco Fanelli and Valeria Albenzio for technical assistance. and Bianca Tino for language revision and editing.
This work was supported by funding from the Italian Ministry of Health “Ricerca Corrente 2017”.
Availability of data and materials
All data generated or analyzed during this study were included in this published article.
Ethics approval and consent to participate
This study was approved by the Ethic Committee of Istituto Tumori “Giovanni Paolo II” with the reference number 234/CE 13–11-2017. Before undergoing routine surgery, all patients signed an informed consent form authorizing the Institute to utilize their removed biological tissue for research purposes according to ethical standards.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.