MMP9 is a prognostic marker for bladder cancer and relates to the Immune Cell Infiltration in the Tumor Microenvironment

doi:10.21203/rs.3.rs-1550825/v1

PDF

Research Article

MMP9 is a prognostic marker for bladder cancer and relates to the Immune Cell Infiltration in the Tumor Microenvironment

https://doi.org/10.21203/rs.3.rs-1550825/v1

This work is licensed under a CC BY 4.0 License

Version 1

posted 15 Apr, 2022

You are reading this latest preprint version

Background:

Bladder cancer (BLCA) is one of the most common human cancers while its incidence is on the rise, especially in female. The tumor microenvironment (TME) acts a significant part in the development and progression of the tumor. We surveyed the expression profiles of BLCA in the Cancer Genome Atlas Program database and investigated genes that might affect the composition of tumor microenvironment (TME). We then analyzed the relationship between genes and tumor prognosis, looking for immune prognostic markers which might lead to new styles for the diagnosis and treatment of BLCA.

Methods

We got the RNA-seq data of BLCA from The Cancer Genome Atlas (TCGA); utilized the Estimation of Stromal and Immune Cells in the Malignant Tissues Expression (ESTIMATE) data algorithm to compute the immune and stromal score. The differentially expressed genes were performed to COX regression analysis and used to construct a protein-protein interaction (PPI) network.

Results

Intersection analysis demonstrated that Matrix Metalloproteinase-9 (MMP9) is the only crossover gene. We computed the percent of tumor infiltrating immune cells (TICs) in BLCA through the Cell-type Identification by Estimating Relative Subsets of RNA Transcripts (CIBERSORT) algorithm. The differences in the high and low expression group of MMP9 in TICs were established. The expression levels of MMP9 in each kind of TICs were then correlated. Via survival, clinical, and enrichment analyses, MMP9 was found to be an immune prognostic factor for BLCA and act a significant part in the tumor microenvironment. Ultimately, the correlation analysis between MMP9 and the immune checkpoint genes indicated that MMP9 can forecast the effect of immunotherapy on BLCA.

Conclusion

MMP9 is an immune prognostic marker gene. MMP9 is involved in regulating the immunological process of BLCA through multiple pathways and is closely related to PD-1 and CTLA4 immune targets. The expression of MMP9 may impact the prognosis of BLCA, indicating that it may be a prognostic-related marker of BLCA.

bladder cancer

MMP9

Bioinformatics analysis

Immune infiltration

Immune checkpoint

Bladder cancer (BLCA) is one of the most common cancers among men, and the incidence also has been on the rise in female recently.[1]. There are many risk factors associated with BLCA. Tobacco exposure is thought to be the most important one[2]. Early diagnosis and treatment can greatly decrease the risk of metastasis and disease-specific mortality. Checkpoint inhibitors, targeted therapies, and immunotherapy provide are now available for patients in various stages of the disease. Thus, identification of biomarkers and therapeutic targets that can accurately identify and predict the prognosis of BLCA has become paramount importance.

The tumor microenvironment is complicated and dynamic environment. Besides addition to basal cells, fibrocytes and endotheliocytes, TME also includes many innate and adaptive immune function cells [3]. Recent studies have been concentrating on the mechanism of tumor microenvironment to identify biomarkers that can serve as predictors for immunotherapeutic response, and to explore new venues for immunotherapeutic interventions [4]. In this research, we downloaded data from the Cancer Genome Atlas (TCGA); used ESTIMATE data algorithm to assign stromal score and immune score for BLCA samples to screen for differential genes. Concurrently, we used Cell-type Identification by Estimating Relative Subsets of RNA Transcripts (CIBERSORT) algorithm to find out tumor infiltrating immune cells (TICs) in the samples to screen the marker Matrix Metalloproteinase-9 (MMP9) gene[5].

The analytic process is shown as Figure 1. The data of TCGA database was download first, and then we used the ESTIMATE method to score the samples and use the CIBERSORT method to evaluate TICs. Then differential genes were screened out via immune and stromal score by using the ESTIMATE method, which applied to perform the PPI networks and univariate COX regression analysis. Then we dug out the MMP9 as the only overlapping gene. Finally, we complete survival analysis, clinical stage correlation analysis, enrichment analysis, TIC correlation analysis and correlation analysis of immune checkpoint genes of MMP9.

2.1 Data Sources

The transcriptome RNA-sequence and clinical data of 412 BLCA patients from TCGA were downloaded and used to follow-up analysis.

2.2 Immune Score and Stromal Score.

We used the R version 4.1.1 of the " ESTIMATE " packages to score each sample's immune score and stromal score. The summation of the above two score was then used for further analysis. The ratio of the scores to the corresponding components is positively correlated. A high score indicated that the components appeared frequently in TME.

2.3 Survival and clinical stage correlation analysis of high and low score group.

Based on above two scores and gene expression levels, we analyzed the clinical data of the 412 cases utilizing the R tool "survival" and "survminer" packages to select the optimal cutoff value and divide them into high and low expression groups. The Kaplan-Meier survival curve was then plotted, and significant differences were assessed through log rank test. p<0.05 was deemed significant. Additionally, we compared samples in the different clinical stages using either Wilcoxon or Kruskal-Wallis rank tests.

2.4 Immune score and stromal score to screen out differential expressed genes

The "limma" package of R tool was used to determine differentially expressed genes (DEGs) between the high and low score groups. A threshold was set to log |fold change|>2 and a false discovery rate (FDR) of <0.01. The "heat map" package of R tool was then chosen to draw a heat map of DEGs.

2.5 Enrichment analysis

We used the "clusterProfiler", “enrichplot”, and “ggplot2” packages of R tool to draw charts for gene ontology (GO) enrichment analysis and for the Kyoto Encyclopedia of Gene and Genomes (KEGG) pathway enrichment analysis of the 280 DEGs. P values and FDR<0.05 denoted significant enrichment.

2.6 Building protein-protein interaction networks

The PPI networks were built via STRING (https://string-db.org/) and visualized through Cytoscape version 3.9.0.

2.7 Construct COX regression analysis

The DEGs was used to perform univariate COX regression analysis through the " survival " packages of R tool.

2.8 Gene set enrichment analysis.

GSEA software version 4.1.0 was used to make the KEGG enrichment analysis for the transcriptome of all tumor samples. P<0.05 and FDR<0.05 represented acceptable results.

2.9 Analysis about relative proportion of TIC and correlation analysis of immune checkpoint genes

CIBERSORT algorithm was employed to compute the relative scale of TICs in the BLCA specimen. Samples with difference in the expression of candidate biomarkers p<0.05 were selected. Correlation analysis of checkpoint genes was then performed.

Finally, immune and stromal scorings were derived using the ESTIMATE method. The percentage of TICs was computed through CIBERSOR. The MMP9 as a TME-related gene was identified by utilizing the univariate Cox regression analysis and protein-protein interaction (PPI) networks.

Analysis process

Correlation between ESTIMATE score and prognosis

The immune and stromal score both have negative correlation with the survival rate of BLCA patients (Figure 2 (a, b)), but the estimate score has no correlation with the survival rate of BLCA (Figure 2 (c)). Next, we compare the clinicopathological staging and gender of the two high and low score groups. Shown as Figure 3, stromal score and estimate score is correlated with stage II and III, stage II and IV in tumor stage respectively (Figure 3 (a, b)), and immune score has no significant correlation with tumor stage (Figure 3 (c)); however, all the score was uncorrelated with the clinical stage (Figure 3 (d-f)), indicating that two types of scores are relevant to the progression of BLCA. Moreover, all scores are also related to gender. Our result shows that the prevalence of men is higher than that of women, which is consistent with previous conclusions.

Filter DEGs by score

There are 461 up-regulated DEGs and 35 down-regulated DEGs screened totally in immune score analysis (Figure 4 (a, c, d)). 496 up-regulated DEGs and 54 down-regulated DEGs screened in stromal score analysis (Figure 4 (b, c, d)). At the same time, we acted GO and KEGG pathway analysis of 280 co-expressed genes. The results of GO analysis show that the co-expressed genes are mainly enriched in the regulation of immune effect process, immune cell proliferation and cytokine activation, etc. and immune response (Figure 4 (e)). KEGG analysis results show that DEGs are mainly enriched in cytokines, the formation of extracellular bactericidal network of neutrophils and phagosomes (Figure 4 (f)), which are all related to immunity. Summing up the two above enrichment analysis results, we conclude that the function of differential genes may be involved in immune reaction, showing that immune reaction is essential to the TME of BLCA.

PPI networks and COX regression analysis of co-expressed genes

We employed the STRING online tool to build a PPI network in order to analyze the relationship between DEGs, and then visualized it via Cytoscape version 3.9.0 (Figure 5 (a)). And then We the top 30 genes were selected in the number of protein interaction network nodes (Figure 5 (b)). What’s over, we completed univariate COX regression analysis to find the relationship between the above genes and BLCA, and finally got 4 genes (Figure 5 (c)). Then, we used the Venn diagram to take the intersection that obtained from the top 30 genes in the number of protein interaction network nodes and the genes gotten via the univariate COX regression analysis (Figure 5 (d)). Finally, we screened out the MMP9.

Expression difference of MMP9 and clinical correlation analysis

After obtaining the gene, we analyzed the expression differences of MMP9 in the tissue of normal and tumor and studied its relationship with BLCA. MMP9 is significantly differently expressed in the tumor tissue and normal tissue, it expressed highly in the tumor tissue (P=0.003) (Figure 6 (a)). The results of the survival analysis showed that the survival rate of patients about MMP9 gene high expression is lower than that with low expression in BLCA (P=0.021)（Figure 6 (b). Besides, MMP9 has significant differences in expression with tumor stage, including stage II and stage III and stage IV (P<0.001) (Figure 6 (c, d)). The higher the staging, the higher the expression. In summary, it may be related to the prognosis of BLCA about the expression of MMP9.

Enrichment analysis of MMP9

Then we make enrichment analysis via GSEA version 4.1.0. The co-expressed genes are mostly enriched in chemokine signaling pathways, cytokine, and cell adhesion molecules (CAMs), all of which are related to immune response in the high expression group of MMP9(Figure 7). Interestingly, the enrichment analysis of MMP9 low expression is not significant. The enrichment analysis results show that there is a very close relationship between MMP9 and TME for bladder cancer.

Correlation between MMP9 and TICs

We analyzed the proportion of 22 types of TICs determined in BLCA case by the way of CIBERSORT algorithm in R tool (Figure 8(a-b)). Then we made a difference and correlation analysis between the high expression groups and low expression groups of MMP9 and TICs to further confirm the role of MMP9 in TME. Difference analysis showed that 7 types of TICs were related to the expression of MMP9 (Figure 8(c)). Correlation analysis externalized that the 10 types of TICs were correlated with the expression level of MMP9 (Figure 8(d)). Then we obtained 6 types of TICs related to the expression of MMP9 when taking the intersection of the two analysis results (Figure 8(e)). Among them, macrophages M0 and eosinophils are positively correlated with the expression of MMP9. On the other hand, activated dendritic cells, mast cells, CD8 T cells, and monocytes are negatively correlated with the expression of MMP9. The above results further prove that MMP9 may have great importance in the immune response of bladder cancer TME.

Correlation analysis of immune checkpoint genes

We performed the correlation analysis of immune checkpoint genes through the R tool, suggesting that it is positively correlated with the expression levels of MMP9 and 36 immune check sites (Figure 9).

Current research show that tumor microenvironment is essential to the development of cancer [3]. It is believed that the characteristics of cancer are determined the cancer cells alone, same as the stromal cells and immune cells surrounding the tumor, and the cytokines they secrete [6]. Even though the diagnostics and treatment methods advance rapidly, the prognosis of BLCA is not satisfactory. The research of BLCA on TME has received increasing attention to scholars, but there is still less research of immune response in BLCA related to immune cells and stromal cells, which are the components of TME. Therefore, it’s particularly significant to explore the function and mechanism of tumor microenvironment in BLCA, looking for effective immune prognostic markers and create new approaches for the diagnosis and therapies of BLCA. In this paper, we probed the expression profiles of BLCA case in the TCGA database, explored genes that affect the composition of tumor microenvironment in BLCA, and analyzed the relationship between genes and prognosis. We use ESTIMATE algorithm to perform stromal score and immune score on BLCA case. The stromal score was negatively correlated with lifetime of BLCA patients. The correlation analysis between the Stromal score and clinical stages indicated that the score shows an upward trend from stage II to stage III, then to stage IV, but it is no significant difference between stage I and other stages. We surmised that this result may arise from the deficiency of samples. There is no significant correlation between the immune score and the tumor stage, but the tumor stage of the two scores is significantly correlated when the two groups are compared. At the same time, we found that there is a significant correlation between the score and gender, which is consistent with the previous research that the occurrence of bladder cancer is related to gender [7]. Then we identified 262 up-regulated and 18 down-regulated DEGs based on the stromal score and immune score. The enrichment analysis of GO and KEGG indicated that DEGs are mainly relevant to the immune response process. Finally, we got the crossover gene MMP9 by constructing a protein interaction network and univariate COX regression analysis. Besides, we studied the differential expression of MMP9 between normal tissues and tumor tissues, and the results also showed significant differences. The survival analysis revealed that the expression of MMP9 was negatively correlated with the survival rate of BLCA patients. GSEA enrichment analysis revealed that high MMP9 expression was mostly enriched in CAMs and chemokine signaling pathways and so on, which are immune-related pathways. All above showed that MMP9 is essential to the immune response of TME of BLCA.

MMPs is a zinc-dependent endopeptidase family with 24 different members of humans, and MMP9 is one of its members [8]. The main role of MMPs is to regulate the extracellular stromal (ECM) by directly degrading ECM proteins (such as collagen, proteoglycan and fibronectin); and to release biologically active proteins such as cytokines, growth factors and chemokines [9]. Among them, the ability of MMP9 degrading extracellular stromal components and regulating the activity of a variety of soluble proteins is essential to the reproduction, growth, development, inflammation, vascular and proliferative diseases [10]. A number of studie have shown that MMP9 has a promoting effect on gastric cancer, colorectal cancer, lung cancer, breast cancer and pancreatic cancer [11–15]. Christine Mehner et al. found that MMP9 promotes tumor angiogenesis and accelerates the metastasis in of breast cancer [16]. Yousef, EM and others found that the expression of MMP9 is different in different subtypes of breast cancer [17]. Other studies have shown that MMP9 induces mesenchymal transformation of pancreatic cancer cells through protease-activated receptor 1 (PAR1) activation and accelerates the development of pancreatic cancer [18]. On the other hand, MMP9 plays the role of a tumor suppressor gene in colitis-related cancers (CAC) by activating the EGFR-Sp1 signaling pathway and the MMP9-Notch1-ARF-p53 axis to maintain the integrity of the epithelial mucosa. The research results are not consistent. Therefore, we speculate that MMP9 has different roles in different types of tumors. Studies have found that compared with the urine of healthy patients, the activity and total amount of MMP9 in the urine of patients with bladder cancer have increased and increased significantly in patients of highly aggressive bladder cancer [19]. It suggests that MMP9 plays a role in promoting cancer in BLCA.

In terms of immune infiltration, the expression of MMP9 has positive correlation with M0 macrophages and eosinophils, and inversely related with CD8 T cells, activated dendritic cells, mast cells, and monocytes. These cells participate in the immune response from various stages, further confirming the role of MMP9 in tumor TME. Correlation analysis between MMP9 and TICs shows that MMP9 can promote cancer by regulating TME of BLCA. Therefore, MMP9 inhibitors can exert anti-cancer effects by inhibiting the expression of MMP9. A variety of studies have proved that natural and synthetic MMP inhibitors (MMPIs) exert anti-cancer effects by inhibiting the cascade of MMPs [20, 21]. In inflammatory and neoplastic diseases such as ulcerative colitis and colorectal cancer, selective MMP9 inhibitors can inhibit the up-regulation of MMP9 to exert a tumor suppressor effect, and MMP9 is expected to become a new anti-cancer target [22]. In addition, we also studied the correlation between the expression of MMP9 and the expression of common immune checkpoint genes to predict the effect of immunotherapy. The results showed that MMP9 is related to most of the immune examination sites of bladder cancer, including CD274, CTLA-4, NFRSF9, LAG3, HAVCR2, etc. Among them, MMP9 has a significant positive correlation with PD-1 and CTLA-4, which indicates that MMP9 may be involved in immunotherapy. A recent study showed that the small molecule MMP2/MMP9 inhibitor SB-3CT enhanced the efficacy of PD-1 or CTLA-4 blockers in the treatment of primary and metastatic tumors [23]. Although anti-PD-1 immunotherapy can improve the survival rate of patients with sensitive tumors, most advanced patients will develop resistance, and we know little about the mechanism of resistance [24]. Studies have shown that MMP-mediated proteolytic lysis may limit the immunosuppressive ability of PD-1 ligands of tumor-related cells. MMP9 can promote the ability of PD-1 response by reducing the expression of tumor PD-L1[25]. Another study of mouse melanoma showed that the TGFβ pathway promotes the MMP9-dependent lysis of PD-L1 surface expression, leading to anti-PD-1 resistance[26]. From this we theorized that MMP9 has a great clinical significance in tumor immunotherapy. This study is based on the prediction model of MMP9 that can predict the status of immune cell infiltration in BLCA and can also be used as a novel biomarker to predict the prognosis of BLCA. But our research has limitations. We need more independent data sets and further experiments to verify our conclusions.

MMP9 plays a critical role in the tumor microenvironment of BLCA and can be used as a novel prognostic marker for BLCA. MMP9 may also be an immune prognostic gene of the TME in BLCA. Further research is needed to confirm the function of MMP9 in BLCA's TME.

BLCA: bladder cancer; TME: tumor microenvironment; TCGA: cancer genome map; DEGs: differentially expressed genes; MMP: stromal metalloproteinase; FDR: false discovery rate; PPI: protein-protein interaction; GO: gene ontology; KEGG: Kyoto Encyclopedia of Genes and Genomes; TICs: tumor infiltrating immune cells; ECM: extracellular stromal.

Acknowledgements

Not applicable.

Authors' contributions

HL and WL constructed the idea of this research and completed the processing and visualization of the data. HL and NF completed the writing of the first draft, SZ and SC participated in the writing and revision of the first draft, and JH and JL reviewed and revised the first draft. All authors have read and agreed to the final manuscript.

Funding

Not applicable.

Availability of data and materials

The datasets generated and/or analyzed during the current study are available in the Cancer Genome Atlas (TCGA) repository, https://portal.gdc.cancer.gov/

Ethics approval and consent to participate

The data of the article is obtained from the public database, and no ethical approval has been applied.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Sung H, Ferlay J, Siegel RL, et al.: Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA: a cancer journal for clinicians 2021; 71(3): 209–249.
Cumberbatch MG, Rota M, Catto JWF, La Vecchia C: The Role of Tobacco Smoke in Bladder and Kidney Carcinogenesis: A Comparison of Exposures and Meta-analysis of Incidence and Mortality Risks. European Urology 2016; 70(3): 458–466.
Wu T, Dai Y: Tumor microenvironment and therapeutic response. Cancer Letters 2017; 387: 61–68.
Gajewski TF, Schreiber H, Fu Y: Innate and adaptive immune cells in the tumor microenvironment. Nature Immunology 2013; 14(10): 1014–1022.
Yoshihara K, Shahmoradgoli M, Martínez E, et al.: Inferring tumour purity and stromal and immune cell admixture from expression data. Nature Communications 2013; 4(1).
Hanahan D, Coussens LM: Accessories to the Crime: Functions of Cells Recruited to the Tumor Microenvironment. Cancer cell 2012; 21(3): 309–322.
Lenis AT, Lec PM, Chamie K, Mshs MD: Bladder Cancer: A Review. JAMA 2020; 324(19): 1980–1991.
Stöcker W, Bode W: Structural features of a superfamily of zinc-endopeptidases: the metzincins. Current opinion in structural biology 1995; 5(3): 383–390.
Connective Tissue Remodeling: Cross-Talk between Endothelins and Ma- trix Metalloproteinases.
Biochemistry and Molecular Biology of Gelatinase B or Matrix Metalloproteinase-9 (MMP-9.
Vesna VDNS: Detection of gelatinase B activity in serum of gastric cancer patients. World journal of gastroenterology: WJG 2006; 12(1): 105–109.
Tekin C, Aberson HL, Waasdorp C, et al.: Macrophage-secreted MMP9 induces mesenchymal transition in pancreatic cancer cells via PAR1 activation. Cellular Oncology 2020; 43(6): 1161–1174.
Blanco-Prieto S, Barcia-Castro L, Páez De La Cadena M, et al.: Relevance of matrix metalloproteases in non-small cell lung cancer diagnosis. BMC Cancer 2017; 17(1).
Liang S, Chang L: Serum matrix metalloproteinase-9 level as a biomarker for colorectal cancer: a diagnostic meta-analysis. Biomarkers in Medicine 2018; 12(4): 393–402.
Li H, Qiu Z, Li F, Wang C: The relationship between MMP-2 and MMP-9 expression levels with breast cancer incidence and prognosis. Oncology letters 2017; 14(5): 5865–5870.
Mehner C, Hockla A, Miller E, Ran S, Radisky DC, Radisky ES: Tumor cell-produced matrix metalloproteinase 9 (MMP-9) drives malignant progression and metastasis of basal-like triple negative breast cancer. Oncotarget 2014; 5(9): 2736–2749.
Yousef EM, Tahir MR, St-Pierre Y, Gaboury LA: MMP-9 expression varies according to molecular subtypes of breast cancer. BMC cancer 2014; 14(1): 609–609.
Tekin C, Aberson HL, Waasdorp C, et al.: Macrophage-secreted MMP9 induces mesenchymal transition in pancreatic cancer cells via PAR1 activation. Cellular Oncology 2020; 43(6): 1161–1174.
Mondal S, Adhikari N, Banerjee S, Amin SA, Jha T: Matrix metalloproteinase-9 (MMP-9) and its inhibitors in cancer: A minireview. European Journal of Medicinal Chemistry 2020; 194: 112260.
Martens E, Leyssen A, Van Aelst I, et al.: A monoclonal antibody inhibits gelatinase B/MMP-9 by selective binding to part of the catalytic domain and not to the fibronectin or zinc binding domains. Biochimica et Biophysica Acta (BBA) - General Subjects 2007; 1770(2): 178–186.
Devy L, Huang L, Naa L, et al.: Selective Inhibition of Matrix Metalloproteinase-14 Blocks Tumor Growth, Invasion, and Angiogenesis. Cancer Research 2009; 69(4): 1517–1526.
Marshall DC, Lyman SK, McCauley S, et al.: Selective Allosteric Inhibition of MMP9 Is Efficacious in Preclinical Models of Ulcerative Colitis and Colorectal Cancer. PLOS ONE 2015; 10(5): e0127063.
Ye Y, Kuang X, Xie Z, et al.: Small-molecule MMP2/MMP9 inhibitor SB-3CT modulates tumor immune surveillance by regulating PD-L1. Genome Medicine 2020; 12(1).
Topalian SL, Hodi FS, Brahmer JR, et al.: Safety, Activity, and Immune Correlates of Anti–PD-1 Antibody in Cancer. New England Journal of Medicine 2012; 366(26): 2443–2454.
Dezutter-Dambuyant C, Durand I, Alberti L, et al.: A novel regulation of PD-1 ligands on mesenchymal stromal cells through MMP-mediated proteolytic cleavage. OncoImmunology 2015; 5(3): e1091146.
Zhao F, Evans K, Xiao C, et al.: Stromal Fibroblasts Mediate Anti–PD-1 Resistance via MMP-9 and Dictate TGFβ Inhibitor Sequencing in Melanoma. Cancer Immunology Research 2018; 6(12): 1459–1471.

No competing interests reported.

PDF

Version 1

posted 15 Apr, 2022

You are reading this latest preprint version

MMP9 is a prognostic marker for bladder cancer and relates to the Immune Cell Infiltration in the Tumor Microenvironment

Status:

Version 1

Abstract

Background:

Methods

Results

Conclusion

Figures

Introduction

Materials And Methods

Result

Discussion

Conclusion

Abbreviations

Declarations

References

Competing Interests

Status:

Version 1