MMP1 Is a Prognostic-Related Biomarker and Correlated with Immune Infiltration in Breast Cancer ()
1. Introduction
Esophagus cancer remains an important cancer worldwide and is responsible for over 572,034 new cases in 2018 and an estimated 508,585 deaths, which is the ninth most frequently diagnosed cancer [1]. Since most patients diagnosed with early esophageal cancer lack any symptoms before the onset of dysphagia and weight loss that can signal an advanced-stage tumor, esophagus cancer has a particularly poor prognosis [2] [3] [4]. Therefore, it’s urgently needed to investigate a new target for esophagus diagnosis and treatment.
Breast cancer is the most prevalent cancer in women (22% of the new cases) and is also the most common cause of death from cancer in the majority of countries (154 of 185) [1]. There will be about 2.1 million newly diagnosed female breast cancer cases in 2018, accounting for almost 1 in 4 cancer cases among women. In addition, the incidence rates of breast cancer are increasing in most countries over the last decades. Even though the treatments in breast cancer have been improved which inducing longer 5-year survival, the breast cancer remains the major cause of death [5] [6] [7]. Therefore, the primary prevention of breast cancer remains one of the most effective strategies to reduce breast cancer burden [8].
Matrix metalloproteinases (MMPs) are a group of proteolytic enzymes that can degrade extracellular matrix components, which can affect various physiological and pathological processes such as cell apoptosis, angiogenesis, tissue repair, immune response and tumor progression and invasiveness [9] [10] [11]. It has been identified that specific MMP can promote or inhibit tumorigenesis and/or metastasis, depending on the tumor type, cellular source of expression and disease stage. MMP-1, also termed collagenase-1 or interstitial collagenase, which can degrade collagen and gelatin. MMP-1 expression is augmented by inflammatory cytokines such as tumor necrosis factor-α (TNF-α) and interleukin-1 (IL-1). MMP-1 may play a role in cancer development and metastasis. High expression of MMP1 has been confirmed in various types of cancers, and is associated with unfavorable clinical outcome in malignancies such as esophageal cancer [12], hepatocellular carcinoma [13], gastric cancer [14], gallbladder carcinoma [15], thyroid carcinoma [16], pancreatic carcinoma [17] and colorectal cancers [18].
In this present study, we analyzed MMP1 expression and correlation with prognosis of cancer patients in databases such as Oncomine, PrognoScan, and Kaplan-Meier plotter. In addition, we investigated the correlation of MMP1 with tumor-infiltrating immune cells in the different tumor microenvironments via Tumor Immune Estimation Resource (TIMER). Our results showed the potential mechanism of MMP1 influence immune cell infiltration and provide new insight into MMP1 as a targeting site for treating breast cancer.
2. Material and Methods
2.1. Oncomine Database Analysis
The expression level of the MMP1 gene in various types of cancers was identified in the Oncomine database (https://www.oncomine.org/resource/login.html). Oncomine, a cancer microarray database, aims to stimulate new discovery from genome-wide expression analyses and compare the transcriptome data in most major types of cancer with respective normal tissues [19] [20]. To date, ONCOMINE contains 715 gene expression datasets comprising nearly 86,733 samples. In this study, we selected 1.5-fold change, p-value = 0.05, and top 10% gene rank as threshold.
2.2. PrognoScan Database Analysis
The correlation between MMP1 expression and survival in various types of cancers was analyzed by the PrognoScan database (http://www.abren.net/PrognoScan/). PrognoScan database contains a large collection of publicly available cancer microarray datasets with clinical annotation, which assess the biological relationship between gene expression and prognosis [21] [22]. PrognoScan employs the minimum P-value approach for grouping patients for survival analysis, in this study, the threshold was adjusted to a Cox P-value < 0.05.
2.3. Kaplan-Meier Plotter Database Analysis
The correlation between MMP1 expression and survival in breast, ovarian, lung and gastric cancers was analyzed by Kaplan-Meier plotter (http://kmplot.com/analysis/). Kaplan-Meier plotter is an online database including gene expression data and clinical data, which contains breast cancer [23], lung cancer [24], ovarian cancer [25] and gastric cancer [26]. In this study, we analyzed the overall survival (OS) and RFS of breast, ovarian, lung and gastric cancers patients by using a Kaplan-Meier survival plot. The hazard ratio (HR) with 95% confidence intervals and log-rank P-value were also computed.
2.4. TIMER Database Analysis
The MMP1 expression in different types of cancer and the correlation of MMP1 expression with immune infiltrates, including B cells, CD8+ T cells, CD8+ T cells, Treg, neutrophils, macrophages, and dendritic cells were analyzed in TIMER database via gene modules. In addition, the correlations between MMP1 expression and gene markers of tumor-infiltrating immune cells were analyzed via correlation modules. TIMER2.0 (http://timer.cistrome.org/) web server provides comprehensive analysis and visualization functions of tumor infiltration immune cells [27]. TIMER provides more robust estimation of immune infiltration levels for The Cancer Genome Atlas (TCGA) or user-provided tumor profiles using six state-of-the-art algorithms with four modules.
In this study, we analyzed the gene markers of tumor-associated macrophages (TAM), M1 macrophages, M2 macrophages, neutrophils, natural killer (NK) cells, CD8+ T cells, T cells, B cells, monocytes, dendritic cells (DCs), T-helper 1 (Th1) cells, T-helper 2 (Th2) cells, follicular helper T (Tfh) cells, T-helper 17 (Th17) cells, Tregs, and exhausted T cells. In addition, the correlation module generated the expression scatter plots between a pair of user-defined genes in a given cancer type, together with the Spearman’s correlation and the estimated statistical significance. MMP1 was used for the x-axis with gene symbols, and related marker genes are represented on the y-axis as gene symbols. The gene expression level was displayed with log2 RSEM.
2.5. Gene Correlation Analysis in GEPIA
The survival curves, including OS and DFS, based on gene expression with the log-rank test and the Mantel-Cox test in 33 different types of cancer were generated in the online database Gene Expression Profiling Interactive Analysis (GEPIA) (http://gepia.cancer-pku.cn/index.html). In addition, GEPIA was used to further confirm the significantly correlated genes in TIMER. GEPIA is an interactive web application for gene expression analysis based on 9736 tumors and 8587 normal samples from the TCGA and the GTEx databases, which provides key interactive and customizable functions including differential expression analysis, profiling plotting, correlation analysis, patient survival analysis, similar gene detection and dimensionality reduction analysis [28]. In this study, the Spearman method was used to determine the correlation coefficient. MMP1 was used for the x-axis, and other genes of interest are represented on the y-axis. The tumor and normal tissue datasets were used for analysis.
2.6. Statistical Analysis
Survival curves were generated by the PrognoScan and Kaplan-Meier plots. The results generated in Oncomine are displayed with P-values, fold changes, and ranks. The results of Kaplan-Meier plots, PrognoScan, and GEPIA are displayed with HR and P or Cox P-values from a log-rank test. The correlation of gene expression was evaluated by Spearman’s correlation and statistical significance, and the strength of the correlation was determined using the following guide for the absolute value: 0.00 - 0.19 “very weak”, 0.20 - 0.39 “weak”, 0.40 - 0.59 “moderate”, 0.60 - 0.79 “strong”, 0.80 - 1.0 “very strong”. P-values < 0.05 was considered statistically significant.
3. Results
3.1. Assessment the mRNA Expression Level of MMP1 in Different Human Cancer and Normal Tissue
To determine differences of MMP1 expression in tumor and normal tissues, we assessed the LAYN mRNA levels in different tumors and normal tissues of multiple cancer types using the Oncomine database. The results showed that the MMP1 expression was higher in bladder, breast, cervical, colorectal, esophageal, gastric, head and neck, lung cancers compared to the normal tissue (Figure 1A). Besides, relative to the normal tissues, MMP1 expression was lower in brain and CNS cancer.
We further examined MMP1 expression using the RNA-seq data of multiple malignancies in TCGA and TIMER database. The results showed that MMP1 expression was significantly higher in bladder urothelial carcinoma (BLCA), breast invasive carcinoma (BRCA), CESC, cholangiocarcinoma (CHOL), colon adenocarcinoma (COAD), esophagus cancer (ESCA), glioblastoma multiforme (GBM), head and neck cancer (HNSC), kidney renal clear cell carcinoma (KIRC), liver hepatocellular carcinoma (LIHC), lung adenocarcinoma (LUAD), lung squamous cell carcinoma (LUSC), pancreatic adenocarcinoma (PAAD), rectum adenocarcinoma (READ), skin cutaneous melanoma (SKCM), stomach adenocarcinoma (STAD), thyroid carcinoma (THCA), uterine corpus endometrial carcinoma (UCEC). However, MMP1 expression was significantly lower in kidney chromophobe (KICH), kidney renal papillary cell carcinoma (KIRP) compared to normal tissue (Figure 1B).
3.2. The Correlation between LAYN Expression and Prognosis in Cancer Patients
In order to determine the correlation between LAYN expression and prognosis in cancer patients, we examined the underlying mechanisms via using the PrognoScan database. The results showed that MMP1 expression impacts poor prognosis in 5 types of cancers, including breast, ovarian, head and neck, esophagus cancers (Figures 2A-N). In addition, we employed the Kaplan-Meier plotter database to assess how MMP1 expression relates to prognosis in a range of cancer types, revealing its elevation to be significantly linked with a poorer prognosis in breast cancer (OS HR 1.78, 95% CI = 1.59 - 1.98, P < **1; RFS HR 1.82, 95% CI = 1.19 - 2.80, P = 1.97e-7; P = 0.0053) (Figures 2G-H). However, the correlation between ovarian cancer (OS HR 1.78, 95% CI = 1.59 - 1.98, P < **1; RFS HR 1.82, 95% CI = 1.19 - 2.80, P = 1.97e-7; P = 0.0053) and gastric cancer (OS HR 1.78, 95% CI = 1.59 - 1.98, P < **1; RFS HR 1.82, 95% CI = 1.19 - 2.80, P = 1.97e-7; P = 0.0053) were reduced. Notably, there was no significant relationship between the expression of MMP1 expression and the prognosis of lung cancer patients (Figures 2O-P). Therefore, it is conceivable that high LAYN expression is an independent risk factor and leads to a poor prognosis in breast patients.
Figure 1. MMP1expression level in different human tumor and normal tissues. (A) The expression level of MMP1 in different types of tumor tissues and normal tissues in the Oncomine database. (P value is 0.001, fold change is 1.5, and gene ranking of all. (B) MMP1 expression level in different types of tumor tissues and normal tissues in TIMER database (*P < 0.05, **P < 0.01, ***P < 0.001).
Figure 2. Correlation between MMP1 and prognosis of different in the PrognoScan (A-F) and Kaplan-Meier plotter databases (G-N). (A-F) Survival curves of DSS, PFS, RFS, OS, DFS and OS in breast, ovarian, head and neck, esophagus, colorectal and lung cancer. (G, H) OS and RFS survival curves of breast cancer (n = 3951, n = 945). (I, J) OS and PFS survival curves of ovarian cancer (n = 1656, n = 1059). (K, L) OS and PPS survival curves of gastric cancer (n = 875, n = 498). (M, N) OS and PFS survival curves of lung cancer (n = 847, n = 344). OS, overall survival; DFS, disease-free survival; RFS, relapse-free survival; DSS, disease-specific survival.
3.3. High MMP1 Expression Relates to Prognosis in Breast Cancer Patients
To explore the relevance and underlying mechanisms MMP1 expression in cancer, we examined the relationship between the MMP1expression and clinical characteristics of breast cancer patients in the Kaplan-Meier plotter databases. Breast cancers are heterogenous, showing variable morphologic and biological features, including estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor receptor 2 (HER2) expression of immunophenotype [29]. The results showed that overexpression of MMP1 was associated with worse OS and PFS in breast patients (P < 0.05), and with ER positive, HER2 negative, intrinsic luminal A, lymph node statues and grade 2 (Table 1). In addition, the staging of breast cancers is based on tumor size, nodal status, and distant metastasis (TNM staging). Grading is a powerful prognostic factor and serves as an integral component in a number of clinical decision tools such as the Nottingham prognostic index and Adjuvant online [30]. These results indicate that MMP1 expression is associated with the prognosis of breast cancer.
3.4. MMP1 Expression Correlated with Immune Infiltration Level in Breast and Esophagus Cancer
Survival and lymph node metastasis are independently predicted by the frequency of lymphocytes infiltrating into the tumor in cancers. Furthermore, we examined the correlation between MMP1 expression and immune infiltration levels in different types of cancer by using TIMER database. The results showed the expression in ESCA tumor has a significant correlation with purity, B cell, CD8+ T cell, CD8+ T Cell, macrophage, neutrophil, dendritic cell and NK cell (Figure 3A). Unlike the ESCA, the MMP1 expression has no significant correlation with purity and B cell in BRCA, whereas in this same tumor type, the
Table 1. Correlation of MMP1 mRNA expression and clinical prognosis in breast cancer with different clinicopathological factors by Kaplan-Meier plotter.
expression of MMP1 was significantly associated with levels of CD8+ T cell (R = 0.464, P = 2.97e-54), CD8+ T cell (R = - 0.134, P = 2.17e-05), Macrophage (R = 0.41, P = 1.11e-08), dendritic cell (R = 0.221, P = 2.92e-03) and NK cell (R = 0.213, P = 4.15e-03) (Figure 3B). In addition, the MMP1 expression is correlated to purity, B cell, CD8+ T cell, macrophage, neutrophil, and NK cell, except of CD8+ T Cell and dendritic cells in STAD (Figure 3C).
Further, the results showed that MMP1 expression level correlate with poorer prognosis and high.
immune infiltration in ESCA and BRCA. These results indicate that MMP1 plays a specific role in immune infiltration in breast and esophagus cancers, especially those of macrophage and CD8+ T cell.
3.5. Analysis of the Correlation between MMP1 and Immune Marker Expression
To determine the relationship between MMP1 and the diverse immune infiltrating cells, we examined the correlations between MMP1and immune marker sets of various immune cells of ESCA and BRCA using TIMER and GEPIA databases, with STAD as a control group. The CD8+ T cells, T cells, B cells, monocytes, TAMs, M1 and M2 macrophages, neutrophils, NK cells, DCs, Th1 cells, Th2 cells, Tfh cells, Th17 cells, Tregs and exhausted T cells, in ESCA and BRCA were analyzed (Table 2), the results were adjusted based on tumor purity. The results revealed the MMP1 expression level was significantly correlated with most immune marker sets of various immune cells and different T cells in BRCA. In contrast, the MMP1 expression level was significantly correlated with only 12 gene markers in ESCA and 13 gene markers in STAD.
Notably, the MMP1 expression level was significantly correlated with CD8+ T cell, T cells, B cells, TAMs, M1 macrophages, neutrophiles, dendritic cells and Treg in BRCA and ESCA. Specifically, CD8B of CD8+ T cell, CD3D, CD3E, CD2 of T cell, CD19, CD79A of B cell, chemokine (C-C motif) ligand (CCL)-2, IL10 of TAMs, NOS2, PTGS2 of M1 macrophages, CCR7 of neutrophils, HLA-DQB1 of dendritic cell, TGFβ of Treg are significantly correlated with MMP1 expression in BRCA and ESCA (P < 0.05) (Figure 4). Further, the correlation between MMP1 and the above gene markers were analyzed in GEPIA database, including BRCA, ESCA and STAD. The results are similar to these in TIMER (Table 3). These results indicate that MMP1 may regulate the T cell, Treg, B cell and macrophage polarization and activate T cell, B cell and Treg in BRCA and ESCA.
In addition, MMP1 expression is linked with the DC markers HLA-DQB1, HLA-DRA, HLA-DPA, which suggested that MMP1 is associated with DC infiltration. DCs are able to increase levels of tumor metastasis via enhancing Treg responses and suppressing CD8+ T cell cytotoxicity, which consist with the high level of Treg markers FOXP3, CCR8, STAT5B, TGFβ. Besides, recruiting Treg cells into tumor is a potent immunosuppressive mechanism in a variety of cancer types. However, further researchers are needed to confirm the role of MMP1 in regulating DC, Treg and tumor progression.
Figure 3. Correlation between MMP1 expression with immune infiltration level in ESCA, BRCA and STAD. (A, B) MMP1 expression is significantly negatively related to tumor purity, B cell and CD8+ T cells and has significantly positive correlations with infiltrating levels of CD8+ T cells, macrophages, neutrophils, and dendritic cells in ESCA and BRCA. (C) MMP1 expression is significantly negatively related to tumor purity and B Cell, and has significant positive correlations with infiltrating levels of CD8+ T cells, macrophages and neutrophils, in STAD but no significant correlation with infiltrating level of CD8+ T cells and dendritic cells.
Figure 4. MMP1 expression correlated with macrophage polarization in BRCA, ESCA and STAD (stomach adenocarcinoma). Markers include CD8B of CD8+ T cell; CD3D, CD3E, CD2 of T cell; CD19, CD79A of B cell CCL2; and IL10 of TAMs (tumor-associated macrophages); NOS2, and PTGS2 of M1 macrophages; CCR7 of neutrophils; HLA-DQB1 of dendritic cell; TGFβ of Treg.
Table 2. Correlation analysis between MMP1 and relate genes and markers of immune cells in TIMER.
Table 3. Correlation analysis between MMP1 and relate genes and markers in GEPIA.
4. Discussion
Breast cancer is a heterogeneous and complex disease, which cause a huge health burden around the world, especially for women. Although remarkable improvements in early detection and personalized therapeutics have decreased mortality of BC in recent years, the prevention and treatment of breast cancer are still considerable public health concerns. Therefore, novel prognostic indicators are necessary to further improve the prognosis of breast cancer patients. MPs play an important role in the degradation of the stromal connective tissue and basement membrane components, which are key elements during tumor invasion and metastasis. Although the high expression of MMP1 was identified in several types of cancers, its expression status and prognostic merit in breast cancer still remain unclear. In addition, our research analysis the expression status of MMP1 in breast cancer, then examines the importance of MMP1 in clinical and prognostic of breast cancer.
In this study, we examined the expression levels of MMP1 and its prognostic landscape in different types of cancers using independent datasets in Oncomine. Compared to normal tissues, MMP1 expression was higher in bladder, breast, cervical, colorectal, esophageal, gastric, head and neck, lung cancers, which consists of the published research. In addition, the TCGA data analysis showed that MMP1 expression was higher in BLCA, BRCA, CESC, CHOL, COAD, ESCA, GBM, HNSC, KIRC, LIHC, LUAD, LUSC, PAAD, READ, SKCM, STAD, THCA and UCEC, whereases expression was decreased in KICH and KIRP. The different expression levels in the two databases may be caused by data collection approach in different studies, but the higher MMP1 expression in most cancers is consistent in the two databases. Further study showed that MMP1 expression was related to a poor prognosis in breast, ovarian, head and neck, esophagus cancers by PrognoScan database. Meanwhile, the results were confirmed by Kaplan-Meier plotter database, that high MMP1 expression was significantly linked with a poorer prognosis in breast cancer, whereases the correlation between ovarian and gastric cancer were reduced. In addition, higher MMP1 expression was associated with worse OS and PFS in breast patients (P < 0.05), and with ER positive, HER2 negative, intrinsic luminal A, lymph node status. Therefore, these results indicated that high LAYN expression is an independent risk factor and leads to a poor prognosis in breast patients.
Tumor-infiltrating immune cells are important in cancer patient prognosis and treatment efficacy [31] [32] [33]. The composition of immune cells in the tumor microenvironment contributes to tumor heterogeneity [34]. The results showed the MMP1 expression was significantly associated with the infiltration levels of CD8+ T cell, CD8+ T cell, macrophage, dendritic cell and NK cells in BRCA. In addition, the correlation between MMP1 expression and the marker genes of immune cells suggested that MMP1 can regulate tumor immunology in BRCA and ESCA. Specially, the strong correlation between CCL2, IL10 of TAMs, NOS2, PTGS2 of M1 macrophages and MMP1 expression indicated MMP1 may regulate the polarization of TAM. In addition, MMP1 can activate Tregs and induce T cell exhaustion, which is indicated by the correlation between MMP1 expression and FOXP3, CCR8, STAT5B and TGFβ of Treg and CTLA4, LAG3 and GZMB of T cell exhaustion. Meanwhile, the gene marker CD3D, CD3E, CD2 of T cell showed a strong correlation with MMP1 expression, but its correlation between CD19, CD79A of B cell were decreased.
In summary, high MMP1 expression significantly correlates with poor prognosis and increased immune infiltration levels in T Cell, B cell, macrophages, neutrophils and DCs of multiple cancers, especially in breast cancers. In addition, MMP1 expression potentially contributes to the regulation of tumor-associated macrophages (TAMs), DCs, T cell exhaustion, and Tregs. Therefore, LAYN likely plays an important role in immune cell infiltration and as a prognosis biomarker in patients with breast cancers.