Elevated MMP9 expression in breast cancer is a predictor of shorter patient survival

Purpose MMP9 is a matricellular protein associated with extracellular matrix (ECM) remodelling, that promotes tumour progression, and modulates the activity of cell adhesion molecules and cytokines. This study aims to assess the prognostic value of MMP9 and its association with cytoskeletal modulators in early-stage invasive breast cancer (BC). Methods MMP9 expression was evaluated by immunohistochemistry using a well-characterised series of primary BC patients with long-term clinical follow-up. Association with clinicopathological factors, patient outcome and ECM remodelling BC-biomarkers were investigated. METABRIC dataset, BC-GenExMiner v4.0 and TCGA were used for the external validation of MMP9 expression. GSEA gene enrichment analyses were used to evaluate MMP9 associated pathways. Results MMP9 immunopositivity was observed in the stroma and cytoplasm of BC cells. Elevated MMP9 protein levels were associated with high tumour grade, high Nottingham Prognostic Index, and hormonal receptor negativity. Elevated MMP9 protein expression correlated significantly with cytokeratin 17 (Ck17), Epidermal Growth Factor Receptor (EGFR), proliferation (Ki67) biomarkers, cell surface adhesion receptor (CD44) and cell division control protein 42 (CDC42). Cytoplasmic MMP9 expression was an independent prognostic factor associated with shorter BC-specific survival. In the external validation cohorts, MMP9 expression was also associated with poor patients’ outcome. Transcriptomic analysis confirmed a positive association between MMP9 and ECM remodelling biomarkers. GSEA analysis supports MMP9 association with ECM and cytoskeletal pathways. Conclusion This study provides evidence for the prognostic value of MMP9 in BC. Further functional studies to decipher the role of MMP9 and its association with cytoskeletal modulators in BC progression are warranted. Electronic supplementary material The online version of this article (10.1007/s10549-020-05670-x) contains supplementary material, which is available to authorized users.


TCGA
The cancer genome atlas TNBC Triple negative breast cancer

Introduction
Matrix metalloproteinases (MMPs) are a family of proteases that have multiple biological functions in cancer development and progression and are abundantly up regulated in breast cancer (BC). MMP9, also known as gelatinase B, plays an important role in extracellular matrix (ECM) remodelling, protein cleavage, and is associated with tumour invasion, metastasis and modulation of tumour microenvironment [1,2]. MMP9 has the capability to degrade collagens, including Type IV collagen [3], which plays a role in basement membrane degradation promoting migration, invasion and metastases. MMP9 is secreted as an inactive pro-enzyme and activation of latent MMP9 is the critical step in its regulation [4]. In vitro and in vivo experiments in human and experimental models of cancers reveal that the increased MMP9 expression is related to tumour progression [5,6]. MMP9 expression is regulated by several molecular pathways such as extracellular signal-regulated kinase (ERK), mitogenactivated protein kinase (MAPK), and phosphoinositide-3-kinase-protein kinase (PI3K); pathways recognised to be altered in BC [7]. In BC, increased epidermal growth factor receptor (EGFR) expression, which is a poor prognostic marker, is implicated with up-regulation of MMP9 [8]. The regulation of the MMPs, particularly MMP9, by p53 has also been documented [9].
MMP9 is activated by the Cell Division Cycle 42 (CDC42), a Rho GTPase family member, which plays a major role in ECM adhesion and cytoskeletal organisation. Elevated CDC42 expression in BC dysregulates the epithelial architecture, which may initiate oncogenes. CDC42 gene silencing studies in BC xenografts showed that CDC42 knockdown decreased the tumour cell invasion and metastasis in vivo [10]. In addition, CDC42 activation stimulated trafficking of membrane-type 1 matrix metalloproteinase (MT1-MMP) in BC cells [11].
CD44 is a transmembrane glycoprotein cell surface adhesion receptor that promotes the secretion of active MMP9. MMP9 gene silencing is shown to change the expression of CD44 and significantly decreases migration and invasion of tumour cells [12]. Increased MMP9 mRNA expression was also observed in CD44 + BC cells compared to CD44 − cells.
In vitro experiments showed that, inhibition of the CD44-MMP axis may provide therapeutic targets for reducing the tumour spread which further establishes a positive association between MMP9 and CD44 expression [10]. Thus, these studies support a role for CD44 in regulating MMP9 and is strongly associated with aggressively behaving tumours. In addition, MMP9 is part of the Rosetta poor-prognosis signature for BC [13] and in silico analysis of BC DNA microarray datasets also showed a positive association of MMP9 with poor outcomes [14]. For these reasons in this study we investigated the association between MMP9, cytoskeletal modulators, and clinicopathological factors of BC at the protein and mRNA levels using multiple well-characterised early-stage BC cohorts.

Study cohort characteristics
This study obtained ethics approval by the North West-Greater Manchester Central Research Ethics Committee under the title; Nottingham Health Science Biobank (NHSB), reference number 15/NW/0685. All samples from Nottingham used in this study were pseudo-anonymised and collected prior to 2006 and stored in compliance with the UK Human Tissue Act. MMP9 protein expression was evaluated using a well-characterised cohort of early-stage (operable) primary invasive BC (n = 675) treated in Nottingham University Hospital NHS Trust as previously described [15]. Clinical and pathological data of patients (including hormone receptor status, histological tumour type, tumour grade, tumour size, lymph node status, Nottingham Prognostic Index and lymphovascular invasion (LVI)) were available. Tumour types were revisited and coded according to the recent WHO blue book [16]. BC-specific survival (BCSS) was maintained on a prospective basis. The expression of a large panel of BC progression/metastasis-related biomarkers, including the Ki67 [17], EGFR [18], CDC42 [19], CD44 [20], PIK3CA [21] and basal marker (cytokeratin 5/6 and 17) [22], was also studied. BC molecular subtypes based on the IHC profile of Oestrogen Receptor (ER), Progesterone Receptor (PR) and Human Epidermal Growth Factor 2 (HER2) were defined as previously described [23]: Luminal A: ER+/HER2− Low Proliferation (Ki67 < 10%), Luminal B: ER+/HER2− High Proliferation (Ki67 ≥ 10%) or ER +/ HER2 + , HER2-positive class: HER2 + regardless of ER status, Triple Negative Breast Cancer (TNBC): exhibiting negative expression of ER, PR, and HER2.

Immunohistochemistry
Specificity of MMP9 antibody was validated by western blotting prior to immunohistochemistry. Cell lysates blots of HEK293 and MCF7 cell lines which were used as positive and negative controls respectively (obtained from the American Type Culture Collection; Rockville, MD, USA) were incubated with anti-MMP9 antibody (Rabbit monoclonal [EP1254], Abcam) at 1:800 dilution for overnight (4 °C) and bands were detected using fluorescent secondary antibodies at (1:15,000) (IR Dye 800CW donkey antirabbit and 680RD donkey anti-mouse, LI-COR Biosciences, UK). Mouse β-Actin (A5441, Sigma-Aldrich; Clone AC-15; Sigma, UK) at 1:5000 was used as a house-keeping protein.
Blocking and visualisation were done as previously documented [27]. The specificity of the antibody was validated with a single specific band at the predicted molecular weight (70 kDa, Supplementary Fig. 1a).
Tumour samples were arrayed onto tissue microarrays (TMAs) as previously described [28]. Full-face BC tissue sections and TMAs were immunoassayed using Novolink Max Polymer Detection system (Leica, Newcastle, UK). In brief, 4 µm BC tissue sections were deparaffinised with xylene and rehydrated through 100% ethanol. Heat-induced (pH6) citrate antigen retrieval was performed and MMP9 antibody (1:100) was incubated for overnight at 4 °C. 3-3′ Diaminobenzidine tetrahydrochloride (Novolink DAB substrate buffer plus) was used as the chromogen. Slides were counterstained with Novolink haematoxylin for 6 min, dehydrated and cover slipped. Normal kidney tissue was used as a positive tissue control, whereas no primary antibody was used as a negative control.
TMA slides were digitally scanned at 20X magnification and viewed using NanoZoomer NDP viewer (Hamamatsu Photonics, Welwyn Garden City, UK). Both the percentage of staining and staining intensity of MMP9 cytoplasmic expression in invasive tumour cells and stromal cells were individually assessed to calculate the final histochemical score (H-score). Staining was double scored blindly by two researchers (NO and IMM) for 25% cores to assess interobserver concordance. Inter-observer agreement was determined, and the interclass correlation coefficient was 0.86, indicating an excellent concordance between scorers. Discordant cases were re-scored by both observers and a final score was agreed.

Statistical analysis
IBM SPSS 24.0 (SPSS IBM Corp, Chicago, IL, USA) software was used for statistical analysis and reported in line with REMARK guidelines [29]. The MMP9 H-score was dichotomised into high and negative/low expression using the median cut-point value. Chi-squared test was used to evaluate the association between MMP9 expression and the clinicopathological parameters.
Kaplan-Meier analysis with log-rank test for significance was performed to assess BCSS. Cox multivariate analysis was performed to test independence from standard prognostic factors in BC (nodal stage, tumour grade, tumour size, ER level of expression, and Ki67). A p value of < 0.05 was considered significant.
Gene expression was analysed in the subset of META-BRIC patients for which MMP9 expression was available. DGE between A) low Vs high MMP9 cytoplasmic expression, and B) low Vs high MMP9 stromal expression were calculated using the Robina implementation of Edge-R statistical tool [30] and DGE with > twofold-change, and adjusted p values < 0.05 were considered significant. The DGE were examined using the online WebGestalt platform and adjusted p < 0.01 considered statistically significant [31,32] and significantly enriched gene ontologies common for both cytoplasmic and stromal MMP9 protein expression. Furthermore, Venny 2.0 [33] was used to identify the overlapping DGEs common to both the cytoplasmic and stromal MMP9 protein expression.

MMP9 protein expression in BC
Full-face BC tissue (n = 10) sections were used to evaluate the pattern of MMP9 protein expression prior to staining of TMAs. This showed uniformly weak MMP9 expression in normal glandular epithelium (Fig. 1a) and ductal carcinoma in situ (DCIS; Fig. 1b). There was slightly increased immunoreactivity of MMP9 observed in the co-existing invasive BC cells (Fig. 1c), in the intravascular tumour cell emboli (Fig. 1d), stromal expression ( Fig. 1e) and Fig. 1f showing No Primary Antibody Control. MMP9 cytoplasmic and stromal protein expression (C+/S+) showed a positive correlation with MMP9 mRNA (Spearman's coefficient 0.218; p = 0.027), this association was confirmed using TCGA data [34,35] (Supplementary Fig. 1b).
On BC TMAs, a variable degree of MMP9 protein expression was observed ( Fig. 1g-k). The H-scores of both cytoplasmic and stromal expressions did not follow a normal distribution. Therefore, for dichotomisation into negative/low and high expression, the median H-scores (0 and 50, respectively) were used. Out of 675 informative TMA cores, 71% showed negative/low expression in the cytoplasm (Fig. 1g) while 29% showed high expression (Fig. 1h). In the stroma 53% showed negative/low expression (Fig. 1i), while 47% showed high expression (Fig. 1j). Positive immunoreactivity was observed in human kidney sections (Fig. 1k).

MMP9 genomic profiling
Consistent with the results obtained for MMP9 protein expression, in the METABRIC and TCGA datasets, MMP9 copy number gain (14.0%) and high mRNA expression (50.1%) was associated with negative ER and PR status, high histological grade, and poor NPI (all; p < 0.01; Table 4). MMP9 copy number gain and high MMP9 mRNA expression was associated with poor prognostic METABRIC Integrative clusters [24] such as 1, 5 and 9 (p < 0.0001). Associations between MMP9 copy number alterations, mRNA expression and clinicopathological variables are summarised in Table 4. External validation of the pooled data using BC-GenExMiner v4.0 was in agreement where high MMP9 mRNA expression was associated with ER, and PR negativity, high histological grade (all; p < 0.001) and poor NPI (p < 0.01). In PAM50 subtypes, high MMP9 expression was associated with basal-like and HER2 + classes (p < 0.0001).

Discussion
Several studies have examined the roles of MMP9 in cancer development, progression, and its impact on patients' survival and prognosis [2,[36][37][38]. Limited research however has been done to demonstrate the mRNA, copy number alterations and protein expression of MMP9 in BC and correlate the findings with clinicopathologic variables and cytoskeletal modulators in an extended cohort of BC patients. The cytoskeleton connects the cytoplasm and the plasma membrane and responds to external stimuli and signals. Cells that grow abnormally and acquire the ability to migrate and invade are the hallmark of    [39]. The degradation of ECM by MMP family members including MMP9 is believed to favour tumour growth, metastasis, invasion and cytoskeletal re-organisation [40]. Thus, the link between MMP9 and cytoskeletal modulators may have important clinical implications. The aim of this study was to determine whether elevated MMP9 expression at mRNA and protein level correlated with tumour grade, BC morphology, cytoskeletal modulators and patient outcome using a large clinical data set with long-term follow-up. High MMP9 expression was associated with an increase in cell proliferation activity indicated by high expression of Ki67, which is associated with poor patient outcome [41]. Also, elevated levels of CK17 was associated with increased expression of MMP9. Breast tumours expressing CK17 and CK5/6 show high mortality rate which clearly implicate a role in tailoring treatment decisions [22]. These findings strengthen the putative role of MMP9 in tumour progression.  The EGFR/PIK3 signalling pathway plays important roles in tumour progression, and these pathways are reported to be frequently altered in BC [42]. We observed that high MMP9 expression was positively associated with EGFR and PIK3 signalling pathways. Elevated EGFR/PIK3 expression was associated with poor patient outcome in TNBC and Luminal B subtypes. Moreover, the PIK3/Akt pathway triggers MMP9 secretion and promotes cell invasion [43]. This implies that MMP9 could have a role in tumourigenic pathways.
MMP9 expression levels were positively associated with the expression of CDC42. CDC42 is a member of the Rho family of GTPases, which plays a role in many of the cellular processes that are associated with tumour progression, such as cell migration, proliferation, cytoskeletal control and vesicular trafficking [44]. A study conducted by Sipes and colleagues; revealed that, CDC42 deletion significantly reduced MMP9 activity. CDC42 is associated with formation of invadopodia, which can act as 'guiding' structures to pave the way for further cytoskeletal protrusions [45,46]. Moreover, in vitro silencing of MMP9 decreased the migratory activity in Adenoid Cystic Carcinoma cells [47]. PIK3 and CDC42 mediate a positive feedback loop to regulate the tumour progression role [48]. Investigating the functional role of MMP9 in regulating the PIK3/CDC42 positive feedback loop in BC, might reveal a new role for MMP9 in the promotion of migration and invasion.
CD44 plays a major role in modulating migration/invasion processes during tumour advancement. Formation of CD44-MMP9 complex in prostate cancer promotes cellular  [49]. High MMP9 expression was positively associated with CD44 in our study. MMP9 acts as a processing enzyme for CD44 which promotes cell motility and, transcriptional knockdown of MMP9 inhibits this interaction [50]. Further studies will be needed to unravel the mechanisms by which MMP9 drives CD44 mediated invasion and tumour progression.
LVI is the presence of cancer cells in lymph vessels and is linked with a poor outcome in BC [51]. Daniele et al.; [52] showed that high levels of MMP9 expression was found in BC tumours with positive sentinel lymph nodes. The sentinel node is the proximal lymph node affected by metastatic cells since it is the first barrier receiving lymphatic drainage from the tumour. In the event of LVI, interactions between ECM and stromal non-tumoural cells induce migration, invasion and metastasis. Moreover, stromal fibroblast reported to secrete MMP9, which in turn may activate tumour cells [53]. Lymphatic networks within lymph nodes spread out before the onset of metastasis [54]. Although in our study we observed a weak association with stromal MMP9 and LVI, it warrants further investigations. Evaluation of stromal MMP9 expression may provide valuable information regarding the early LVI events.
The oncogenic expression of HER2 induces BC disease progression and invasiveness, which is hypothesised to increase MMP9 activity [53]. In this study, high MMP9 mRNA expression was associated with poor prognostic parameters including higher tumour grade, ER-/PR-, HER 2 + tumours and TNBC tumours. MMP9 was also highly expressed in basal type tumours over luminal A and luminal B tumours. This is consistent with another study which found MMP9 as differentially expressed between molecular subsets of tumours and as a feature of TNBC and HER2 + BCs [55]. Increased MMP9 expression was associated with the poor prognostic category of the NPI. Hence, as observed in our study, MMP9 is a marker indicative of unfavourable prognosis in BC. Stromal invasion requires degradation of the basement membrane. MMP9 cleave the basement membrane type IV collagen and promote tumour invasion and metastasis [56,57]. The correlation between high expression of MMP9 and collagen type XXVIIα1, XXIIIα1, XXIα2 may also induce basement degradation. Collagens, matricellular proteins and CDC42 at mRNA level showed strong positive association with MMP9. Results on pathway analysis confirmed the significant association with collagens, extracellular matrix and cytoskeletal part gene ontologies. This implies that MMP9 plays a role in tumourigenic pathways and could be a marker of poor prognosis in BC.
This study reveals that MMP9 at both proteomic and transcriptomic levels is associated with poor prognostic characteristics and short-term survival outcomes in BC. Cytoplasmic MMP9 expression on its own and combined Pathway analysis for gene ontology significantly associated with cytoplasmic and stromal MMP9 protein expression. The gene panel were significantly associated with extracellular matrix related gene ontologies cytoplasmic and stromal expression was predicative for shorter BCSS in the whole cohort, MMP9 did not show any association with patient outcome. Stromal cell genetic stability plays significant role in modulating the tumour microenvironment contributing to drug resistance and tumour relapse [58]. Moreover, expression of tumour markers in the stroma is found to be closely associated with tumour progression and patient outcome [59]. These findings suggest that MMP9 in both stromal cells and tumour cells might play an important role in the BC progression. In agreement with our study high expression of MMP9 was associated with poor patient survival [60]. Thus, the current study provides definitive evidence that MMP9 is an independent prognostic marker of poor short-term clinical outcome in primary BC and supports further mechanistic and translational studies to target MMP9.
Funding This research received no specific grant from any funding agency in the public, commercial, or not-for-profit sectors.

Compliance with ethical standards
Conflict of interest All authors declare that they have no conflicts of interest. Availability of data and materials The authors confirm the data that has been used in this work are available on reasonable request.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.