==== Front
BMC Cancer
BMC Cancer
BMC Cancer
1471-2407
BioMed Central London

11100
10.1186/s12885-023-11100-0
Research
Elevated levels of MMP12 sourced from macrophages are associated with poor prognosis in urothelial bladder cancer
Kerzeli Iliana K. 1
Kostakis Alexandros 1
Türker Polat 2
Malmström Per-Uno 2
Hemdan Tammer 2
Mezheyeuski Artur 34
Ward Douglas G. 5
Bryan Richard T. 5
Segersten Ulrika 2
Lord Martin 1
Mangsbo Sara M. sara.mangsbo@farmaci.uu.se

1
1 grid.8993.b 0000 0004 1936 9457 Department of Pharmacy, Science for Life Laboratory, Uppsala University, Uppsala, Sweden
2 grid.8993.b 0000 0004 1936 9457 Department of Surgical Sciences, Uppsala University, Uppsala, Sweden
3 grid.8993.b 0000 0004 1936 9457 Department of Immunology, Genetics and Pathology, Uppsala University, Uppsala, Sweden
4 grid.411083.f 0000 0001 0675 8654 Vall d’Hebron Institute of Oncology, Barcelona, Spain
5 grid.6572.6 0000 0004 1936 7486 Bladder Cancer Research Centre, Institute of Cancer & Genomic Sciences, College of Medical & Dental Sciences, University of Birmingham, Birmingham, UK
30 6 2023
30 6 2023
2023
23 60514 11 2022
21 6 2023
© The Author(s) 2023
https://creativecommons.org/licenses/by/4.0/ 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://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
Background

Urothelial bladder cancer is most frequently diagnosed at the non-muscle-invasive stage (NMIBC). However, recurrences and interventions for intermediate and high-risk NMIBC patients impact the quality of life. Biomarkers for patient stratification could help to avoid unnecessary interventions whilst indicating aggressive measures when required.

Methods

In this study, immuno-oncology focused, multiplexed proximity extension assays were utilised to analyse plasma (n = 90) and urine (n = 40) samples from 90 newly-diagnosed and treatment-naïve bladder cancer patients. Public single-cell RNA-sequencing and microarray data from patient tumour tissues and murine OH-BBN-induced urothelial carcinomas were also explored to further corroborate the proteomic findings.

Results

Plasma from muscle-invasive, urothelial bladder cancer patients displayed higher levels of MMP7 (p = 0.028) and CCL23 (p = 0.03) compared to NMIBC patients, whereas urine displayed higher levels of CD27 (p = 0.044) and CD40 (p = 0.04) in the NMIBC group by two-sided Wilcoxon rank-sum tests. Random forest survival and multivariable regression analyses identified increased MMP12 plasma levels as an independent marker (p < 0.001) associated with shorter overall survival (HR = 1.8, p < 0.001, 95% CI:1.3–2.5); this finding was validated in an independent patient OLINK cohort, but could not be established using a transcriptomic microarray dataset. Single-cell transcriptomics analyses indicated tumour-infiltrating macrophages as a putative source of MMP12.

Conclusions

The measurable levels of tumour-localised, immune-cell-derived MMP12 in blood suggest MMP12 as an important biomarker that could complement histopathology-based risk stratification.

As MMP12 stems from infiltrating immune cells rather than the tumor cells themselves, analyses performed on tissue biopsy material risk a biased selection of biomarkers produced by the tumour, while ignoring the surrounding microenvironment.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-023-11100-0.

Keywords

Biomarkers
Matrix metalloproteinase 12 (MMP12)
Urothelial bladder cancer
Single-cell transcriptomics
Proteomics
Proximity Extension Assay (PEA)
Macrophages
Uppsala UniversityOpen access funding provided by Uppsala University.

issue-copyright-statement© BioMed Central Ltd., part of Springer Nature 2023
==== Body
pmcBackground

Urothelial bladder cancer (UBC) is the 10th most prevalent cancer worldwide and the 4th among men [1]. Long-term surveillance of patients at risk of recurrence or progression is costly for healthcare providers and burdensome for patients, necessitating repeated cystoscopies and imaging. However, existing classifications, which depend solely on clinicopathological variables [2, 3], are insufficient for the precise prediction of treatment responses or disease trajectory. Consequently, there is an urgent and unmet need to identify accurate local and systemic predictive and/or prognostic biomarkers. Proximity Extension Assays (PEAs) are sensitive tools for biomarker discovery, but are yet to be explored in-depth in liquid biopsies from UBC patients.

UBC is associated with carcinogen exposure, often resulting in tumours characterised by high mutational load [4, 5]; such features may lead to immunogenicity and an immunotherapy-responsive disease by increasing the likelihood of HLA-presentable tumour neoantigens that are able to initiate an immune response [6–8]. High numbers of macrophages [9–11], neutrophils [12], and myeloid-derived suppressor cells [13, 14] in UBC are further associated with poor survival and limited treatment responses.

This study aimed to investigate whether a sensitive protein detection method such as PEA technology could aid in prognostic biomarker discovery in liquid biopsies by exploring the systemic and local proteomic contexture of urothelial cancer in relation to immune- and tumour-related processes. To achieve this, we utilised an immuno-oncology-focused multiplexed PEA panel to analyse plasma and urine samples from a Swedish UBC patient cohort. Findings were validated using PEA data from an independent UK UBC cohort [15], and publicly-available RNA-sequencing data were interrogated to identify gene expression changes and the cellular sources of the prognostic biomarkers.

Methods

A graphical illustration of the overall analysis workflow and aims was created with BioRender.com and is presented in Fig. 1.Fig. 1 A graphical overview of the methods, cohorts and analysis workflows presented in this study. An initial 1) screen of immune-oncology-related proteins in plasma/serum was performed using OLINK, which was followed by 2) the use of transcriptomics to further investigate and validate the proteomic findings. Patient numbers are referred to as n with public datasets presented next to each transcriptomic assay

Study populations

Patients included in the discovery cohort were diagnosed with UBC between 1986 and 1996 at the University Academic Hospital of Uppsala, Sweden. The study was performed in accordance with the Declaration of Helsinki and approved by the regional ethical review board of Uppsala (ID 2015/143/1). Plasma samples (collected in heparin tubes) from 90 of these patients and paired urine samples (morning urine) from 40 (33/40 passed QC) were stored at -80 °C and were not thawed until their analysis with PEA as described below. Reference pathology was performed by an experienced pathologist according to the 2004 WHO Classification of Bladder Tumours [16]. Clinical and histopathological characteristics are presented in Table 1. The clinical follow-up of the patients was conducted in accordance with the European and National Guidelines with regular cystoscopies and transurethral resection of the tumours when recurrences had been detected. For validation, we utilised PEA data derived from the pre-treatment serum samples of 79 UBC patients recruited to the West Midlands’ Bladder Cancer Prognosis Programme (BCPP), UK (ethical approval: 06/MRE04/65) [15, 17].Table 1 Clinical characteristics and demographic data of the main screening urothelial bladder cancer cohort were used in this study. Within the Recurrences group, “Few” is defined as less than three recurrent tumours within 18 months, whereas “Frequent” is defined as three or more recurrences within the same period [18]

	Plasma n = 90	Urine n = 33	
Sex	
 Female	23 (26%)	11 (33%)	
 Male	67 (74%)	22 (67%)	
 Age	72 (64, 77)	72 (64, 78)	
Stage Grade	
 Ta Low grade	21 (23%)	7 (21%)	
 Ta High grade	10 (11%)	3 (9%)	
 T1 Low grade	3 (3%)	1 (3%)	
 T1 High grade	27 (30%)	10 (30%)	
 T2, T3, T4 High grade	29 (32%)	12 (36%)	
Invasiveness status	
 NMIBC	61 (68%)	21 (64%)	
 MIBC	29 (32%)	12 (36%)	
Recurrences	
 No	26 (29%)	11 (33%)	
 Few	37 (41%)	10 (30%)	
 Frequent	16 (18%)	6 (18%)	
 Unknown	11 (12%)	6 (18%)	
Progressed to MIBC	
 No	65 (72%)	24 (73%)	
 Yes	18 (20%)	7 (21%)	
 Unknown	7 (8%)	2 (6%)	
n (%); Median (IQR)

Multiplex proximity extension assay

Levels of 92 immuno-oncology-related proteins (OLINK Immune Oncology panel I) in urine and plasma were quantified by PEA [19] (Proseek platform, OLINK) at the Clinical Biomarker facility (Science for Life Laboratory, Uppsala, Sweden). Briefly, 92 proteins were detected in 1 μl of sample in a multiplex assay utilising concomitant binding of matched oligonucleotide-conjugated antibody pairs. When in close proximity, the corresponding antibody oligonucleotide pairs hybridize, resulting in a unique DNA template. Following a polymerase amplification reaction, the sequences are quantified by microfluidic qPCR and presented as Normalised Protein eXpression (NPX), an arbitrary unit based on log2-transformed Ct values (http://www.olink.com).

Microarray gene expression analysis

To complement the small cohort used for the OLINK assay and to examine whether the observed differences could be noted in a larger cohort, made up of tissue biopsies, microarray data for the tumours from 308 UBC patients (in the form of non-normalised, raw intensity values along with their detection p-values) were retrieved from the gene expression omnibus (GEO) repository (GSE32894); clinical information regarding the patients and their samples was obtained from the associated publication [20]. Two of the samples were excluded from further analysis because the stage was not reported. Quantile normalisation, probe filtering and batch effect correction were performed according to the methods supplied by the original authors. Differential expression analysis was conducted in an R 4.1 environment (limma, v.3.48.0 [21], R). A linear model comprising a design matrix without intercept for stage of each sample, as well as a contrast matrix stating the comparisons to be made were fitted to the data to calculate fold changes, with positive values indicating higher expression in muscle invasive bladder cancer (MIBC) and negative values reflecting higher expression in non-muscle invasive bladder cancer (NMIBC). Patients with available survival data (n = 224) were selected for subsequent survival analyses (survival 3.2.7, R), and were split into highly- or lowly-expressing for each queried protein with cutoffs determined by maximally selected rank statistics (maxstat v.0.7.25, R).

Single-cell RNA-sequencing data analysis

Analysis of single-cell RNA-sequencing data from a murine UBC model and human UBC patients were explored with the goal of identifying the putative source of the mRNA transcripts coding for selected proteins from the OLINK assay. Two sets of single-cell RNA-sequencing data were downloaded and processed: unique molecular identifier (UMI) count matrices originating from bladders of a urothelial cancer mouse model [22] (GSE174182) and raw sequence reads derived from patient tumour material from the European Nucleotide Archive under the accession code PRJNA662018 [23]. To generate count matrices for patient samples, the raw sequences were aligned to the GRCh38 reference genome through CellRanger 2.2.0. Downstream processing followed the Seurat v4.0 pipeline [24]. Murine cells were filtered to retain higher quality cells (> 200 & < 8000 uniquely identified genes, < 25% of reads mapped to mitochondrial genes), whereas the filtering thresholds for the human samples remained identical to the original publication (< 6000 unique genes, > 1000 UMI counts per cell and < 10% reads stemming from mitochondrial DNA). SCTransform [25] was employed to normalise and cluster the murine-derived cells, with clusters annotated based on their top differentially expressed genes. As for the human-derived cells, an additional integration step was performed prior to clustering due to the highly patient-specific clusters observed.

Statistical analyses

For group comparisons, a two-sided Wilcoxon rank-sum test was used with adjustment for multiple testing by the Benjamini–Hochberg procedure (PEA data), or Bonferroni (mRNA data). Due to the internal normalisation procedures applied for the quantification of each protein, NPX values cannot be compared between proteins. Instead, the strength of association between protein markers was interrogated by Spearman rank correlation within or between respective sample types (serum and urine). Receiver operating characteristic (ROC) curves were constructed from logistic regression. Proteins important for overall survival (OS) were retrieved from a random survival forest model (RSF) [26]. The ranking of proteins was based on their importance in classifying the data in the model using their variable importance (VIMP) measurements. The assembly model used herein was built on the time-to-event data of the cohort and used the intersection of the 63 markers shared between the OLINK Immuno-Oncology panel I and the 422 OLINK markers investigated in the serum of UBC patients by Bryan et al. [15] (Supplementary Table 1). NPX values were scaled and trained on the Swedish cohort of 90 plasma samples, followed by validation against 79 serum samples from the UK cohort. The machine learning model was developed in the R milieu (RStudio Team 2019, Boston, USA) using the packages caret (v.6.0.86) and randomForestSRC (v.2.9.3). Maximally selected rank statistics (maxstat v.0.7.25, R) were used to dichotomize protein concentrations (high/low) in Kaplan–Meier log-rank test analyses. Multivariate overall survival estimates were calculated from Cox proportional hazard regression models comprising the covariates stage, grade, sex and age (survival 3.2.7, R).

Results

Large heterogeneity of proteomic profiles in urine and plasma of urothelial bladder cancer patients

Plasma and urine samples from a treatment-naïve UBC patient cohort were subjected to a highly sensitive targeted proteomic analysis (Immuno-oncology panel I, OLINK). Clinical and demographic data are presented in Table 1. From the initial screening of 92 protein panel markers, a subset was excluded from further analyses because NPX levels were below the limit of detection (LOD) in a large proportion of samples. Markers below LOD in plasma: IL-1α, IFN-γ, IL-33, IL-21, IL-35, CXCL12, FGF2, IL-2, IFN-β, IL-4, IL-13, CD28; and in urine: IL-5, IL-13, CXCL12, TNF, CD40L, IL-10, IL-4, IFN-β, NOS3, IL-21, IFN-γ, IL-2, IL-33, IL-35. Unsupervised hierarchical clustering and uniform manifold approximation and projection (UMAP) dimensional reduction analyses of the remaining proteins in plasma and urine revealed immuno-oncology profiles predominantly driven by differences among individuals as opposed to clinical variables such as stage, grade, sex or age (Fig. 2A-D). The strength of the association between protein markers was investigated by Spearman rank correlation. The analysis revealed several proteins with significantly correlated urine and plasma concentrations from the 33 patients with paired samples. Following adjustment for multiple testing, 12 proteins involved in vital immune cell functions remained, with the strongest correlations observed for the surface-bound receptor/ligands LAMP3, MICA/B, PD-1 and CD27 (Fig. 2E).Fig. 2 Immune-oncology protein marker profiles in plasma and urine of UBC patients. Unsupervised hierarchical clustering of patients based on (A) plasma and (B) urine protein NPX levels with corresponding clinical information presented above. Uniform manifold approximation and projection analysis of the overall protein marker levels in (C) plasma and (D) urine with Ta-T1 (NMIBC) and T2-T4 (MIBC) indicated in purple and orange, respectively. Circles represent low-grade and triangles high-grade tumours. E The correlation matrix for the 12 proteins significantly correlated between plasma and urine in paired samples (Spearman's rank-order correlation p < 0.05, p-adj < 0.05). Yellow represents a positive correlation and blue negative correlation with the size of the circle proportional to its R-value

Stage-specific protein elevation in plasma and urine

NPX values of proteins were compared between NMIBC and MIBC patients. Systemic levels of matrix metalloproteinase 7 (MMP7) and C–C motif chemokine ligand 23 (CCL23) were significantly elevated in the plasma of MIBC patients compared to patients diagnosed with NMIBC (p = 0.028 and p = 0.03, respectively) (Fig. 3A); however, these levels were not statistically different in urine (p = 0.87 and p = 0.59, respectively). Instead, protein levels of the two TNF-receptor superfamily members, CD27 and CD40, were found to be present at significantly lower levels in the urine of MIBC patients compared to urine from NMIBC patients (p = 0.044 and p = 0.040, respectively) (Fig. 3B). To be noted, following correction for multiple testing, neither of the reported proteins reached significance (q < 0.05, FDR). When assessed as stage-specific classifiers in terms of sensitivity and specificity by calculating the area under the curve (AUC), the two urine markers CD27 (AUC = 0.71) and CD40 (AUC = 0.72) performed relatively well in combination (AUC = 0.75), whereas the plasma markers CCL23 (AUC = 0.64) and MMP7 (AUC = 0.64) alone, or in combination (AUC = 0.65) performed poorly (Supplementary Fig. 1).Fig. 3 Differential protein expression between NMIBC and MIBC patients. A MMP7 and CCL23 NPX levels in plasma and (B) CD27 and CD40 levels in urine

High MMP12 plasma levels are associated with poor overall survival outcome

Compared to univariable or multivariate regression analysis with an individual protein, the RSF model included all proteins to predict survival status. To evaluate the accuracy of the model, data from a UK cohort [15] were used. The model included the 63 proteins shared between the Swedish and UK cohorts (Supplementary Table 1) and the respective patients’ overall survival status. Data from the Swedish cohort (90 plasma samples) were used as a discovery cohort to train the model with data provided by Bryan et al. [15] used as a validation cohort. Although not a perfect validation cohort (Swedish cohort = plasma, UK cohort = serum), the machine learning model presented an error rate of 32% given the different sample types. Here, MMP12 was ranked as the most important contributor to OS, as shown in the variable importance (VIMP) plot (Fig. 4A). The model’s main purpose was to identify prognostic biomarkers; thus, it did not include clinical variables. Instead, the top-ranked proteins were validated in a subsequent multivariate Cox-regression analysis comprising the covariates stage, grade, sex, and age (Supplementary Fig. 2, Supplementary Table 2A-B). When assessed as a continuous variable, MMP12 remained significantly associated with a poor outcome (HR = 1.8, CI:1.3–2.5, p < 0.001) in the training cohort (Fig. 4B) and in the validation cohort (HR = 2.20, CI:1.42–3.41, p < 0.001) (Supplementary Fig. 3, Supplementary Table 2C-D). The association between high levels of MMP12 and poor OS was maintained in the independent validation cohort (Fig. 4C). An interesting possibility with RSF models is the ability to study how MMP12 levels influence the predicted survival over time. Here, as a proof-of-concept, survival probability is presented at 6 months, 1 year and 5 years using the aforementioned model (63 markers, clinical information excluded) (Fig. 4D).Fig. 4 Proteins in plasma and serum are associated with overall survival. A The 10 most important plasma proteins for survival outcome were extracted from the plasma cohort RSF model (VIMP: Variable Importance). B Forest plot from the Swedish training data for MMP12, stage, grade, age and sex. High MMP12 levels are significantly associated with poor overall survival (hazard ratio > 1, p < 0.05). C Kaplan–Meier Survival Estimates curves from two urothelial cancer patient cohorts based on high or low MMP12 protein levels. Cut-off levels were calculated by maximally selected rank statistics. D RSF model predicted overall survival probability at different MMP12 NPX levels after 6 months, 1 year and 5 years

Transcriptomic validation of MMP7 and MMP12 expression variation between NMIBC and MIBC

Due to the clinical importance of systemic levels of MMPs in this cohort, illustrated by stage specificity and survival differences, we sought to investigate their expression at the transcript level in datasets generated with complementary methods. We utilised a microarray-based dataset from a UBC patient cohort [20], as well as two single-cell RNA-sequencing datasets: a clinically relevant murine urothelial cancer model developed by our group [22] and a patient cohort comprising four non-invasive and four invasive tumours [23].

Analysis of the microarray gene expression data from 213 NMIBC and 93 MIBC patient tumour biopsies revealed that MMP7 and MMP12 mRNA expression levels were higher in MIBC. An approximately 1.85-fold increase for MMP7 (p-adj = 0.0018) (Fig. 5A) and a 1.89-fold increase for MMP12 (p-adj = 0.0001) (Fig. 5B) were observed for MIBC. Survival analyses in this cohort showed that patients with high expression of MMP7 had statistically significant, poorer OS (p = 0.018) (Fig. 5C). The same trend was noticeable for patients with high expression of MMP12, however this comparison marginally failed to reach statistical significance (p = 0.055) (Fig. 5D). To be noted, the cohort examined was heavily skewed towards NMIBC patients.Fig. 5 MMP7 and MMP12 expression in UBC patient biopsy tissues and an autochthonous murine urothelial cancer model. Normalised gene expression levels in tumour samples (n = 306) were analysed with microarray probes for MMP7 (A) and MMP12 (B) in relation to the patient’s stage. Kaplan–Meier plots for selected patients (n = 224) of the aforementioned cohort based on high or low mRNA levels of MMP7 (C) and MMP12 (D). Mmp12 mRNA expression at the single-cell level from a murine urothelial cancer model is significantly higher at the MIBC stage, irrespective of sex (E). Bar plots illustrate the percentage and the total number of macrophages exhibiting non-zero expression of Mmp12 (F) for the different stages and sexes. The average log2 fold changes and the BH-adjusted p-values were computed using limma for the microarray data, while the respective fold changes and Bonferroni-corrected p-values were calculated through a Wilcoxon rank-sum test provided by Seurat for the single-cell RNA-sequencing data

Thus, we proceeded to study single-cell RNA data to fully capture the cellular source of the MMP findings. Subsequent examination of a public human urothelial cancer dataset (PRJNA662018) at the single-cell RNA level determined the source of these transcripts to be epithelial cells for MMP7 and macrophages for MMP12 transcripts (Supplementary Figs. 4 and 5). Notably, weak MMP12 mRNA expression from epithelial cells was present in less than 10% of the total cells.

Single-cell transcriptomic analysis of murine tumours further corroborated the source of human tumour-derived MMP12 expression to macrophages, as also murine Mmp12 expression was mapped to macrophages (Supplementary Fig. 5). However, the expression of Mmp7 was not adequately identified in the murine dataset. A total of 2869 cells were annotated as macrophages, 1355 of which originated from NMIBC and 1514 from MIBC murine samples (Fig. 5E). Mmp12 expression was increased > sevenfold in male-derived macrophages (p-adj = 1.63–40) in MIBC compared to NMIBC and almost sixfold in female-derived macrophages (p-adj = 1.47–16) in MIBC compared to NMIBC (Fig. 5E). In addition to the elevated levels of Mmp12 in MIBC, the percentages of Mmp12-expressing macrophages also increased as the disease progressed (Fig. 5F).

Discussion

The majority of UBC patients require intensive, long-term surveillance which increases morbidity and healthcare-related costs. With the current standard-of-care treatments for NMIBC, 5-year recurrence-free survival ranges from 23 to 43%, while the 5-year progression-free survival ranges from 54 to 93% [27]. However, based solely on clinical and pathological characteristics, patient tumour recurrence and survival outcomes are difficult to predict. Although NMIBC and MIBC are molecularly complex and heterogeneous [28], recent developments in transcriptomic classification offer promise for better treatment stratification and personalised therapy [29, 30]. These approaches rely upon laboratory-based deep sequencing of tumour tissue, necessitating invasive biopsies, specialist public or commercial genomics laboratories, and a turnaround time of 5–10 days. Accurate protein-based prognostic biomarkers assessed in liquid biopsies could overcome some of these limitations and easily permit repeated assessments during the patient’s disease course.

In this study, we used an assay which combines high technical sensitivity and specificity with a high-throughput multiplex format, permitting the assessment of 92 protein markers of important immune- and tumour-related processes; the clinical samples were collected at the first diagnosis of UBC before therapeutic interventions. Nevertheless, stage-specific patient classification or clustering based on protein levels was not achieved – a result of the heterogeneity in proteomic patient profiles which did not follow stage- or grade-related patterns.

Released extracellular vesicles or shed immune-related receptors/ligands such as MICA/B, which was one of the 12 identified immune-related receptor/ligand proteins positively correlating between plasma and urine, may have a systemic impact on the immune response with an increased risk of developing immune evasion later in the disease trajectory in several cancer types [31].

Interestingly, we found that receptor-bound proteins were elevated in our liquid based analysis. The origin of elevated soluble CD40 in the urine of NMIBC compared to MIBC patients is likely to be antigen-presenting cells and its decreased levels at MIBC might be a biomarker of immune suppression in the local microenvironment, there are discrepancies as to whether increased expression of CD40 has beneficial or adverse effects on cancer prognosis and therapy responses [32, 33]. On the other hand, soluble CD27 is reduced in the serum of cancer patients compared to healthy individuals, while in the serum of cancer patients post immunotherapy, soluble CD27 correlates with improved survival probability [34].

MMP7 and CCL23, involved in tissue degradation and T cell recruitment signalling, were found to be significantly elevated in the plasma of MIBC patients compared to NMIBC patients. Significantly higher MMP7 levels were also observed at MIBC compared to NMIBC in the UK cohort [15], in gene expression data of an independent patient cohort [20], and in other studies [35]. In our study, analysis of tumour tissue gene expression indicated that MMP7 was expressed at higher levels in tumours of MIBC patients compared to NMIBC.

Despite stage-specific differences, none of the aforementioned proteins reached a high UBC stage classification accuracy (AUC > 0.85). High levels of MMP7 have been reported to be associated with shorter overall survival in UBC [15, 36, 37]. However, neither MMP7 nor CCL23, CD27 or CD40 were significantly associated with overall survival when included as continuous variables in multivariable Cox-regression analyses (data not shown) in our dataset.

Instead, analysis of the data by application of an RSF model ranked increased plasma levels of MMP12 as the main driver for overall survival in the cohort, and as an independent prognostic biomarker associated with decreased overall survival of UBC patients following multivariable analysis comprising clinically relevant confounders, such as stage and grade. Due to the high inter-patient heterogeneity and collinearity between markers, a non-parametric and non-linear method, such as random forests, offers an alternative to classical approaches that are affected by many noise variables and high-order interactions between investigated variables [26]. Despite the differences in liquid biopsies (plasma and serum), sampling procedures, and the use of different multiplex panels, the association between high systemic MMP12 levels and poor survival were recapitulated in the UK cohort.

To evaluate if MMP7 and MMP12 can be identified as biomarkers when using clinical biopsies, we investigated their gene expression in 306 UBC patients [20]. Here, we found that MIBC had higher MMP7 mRNA levels overall. Interestingly, compared to MMP12 protein levels, MMP12 mRNA expression was significantly higher at the MIBC stage. However, when studying the impact of MMP12 on survival in tissue biopsies analysed by transcriptomics, we could not corroborate the findings of the liquid biopsy data. This could be due to the bias introduced in the selection of the malignant cells retrieved for pathological scoring as well as the difference between transcriptomic versus proteomic analyses or the high number of NMIBC patients compared to late stage patients in the transcriptomics cohort. A biopsy for clinical diagnosis also focuses on sampling the most malignant tumor area, rather than the complete tumour microenvironment.

To explore the source of MMP7 and MMP12 in UBC, we analysed UBC patient tumour samples that underwent single-cell sequencing and identified epithelial cells as the source of MMP7 and macrophages as the source of MMP12 mRNA transcripts. Previous studies have shown that increased numbers of tumour-infiltrating macrophages correlate with advanced tumour stage and Bacillus Calmette–Guérin (BCG) treatment failure in UBC [38, 39]. The importance of sex in UBC in the context of immunomodulatory therapies is still not fully explored [40], additionally, macrophages are usually expressed at a higher number and in a more active state in females compared to males [41], thus we included that variable in the single-cell RNA analysis for MMP12. However, we did not see any sex dependent difference in that aspect.

The family of MMPs is largely orthologous between humans and mice, sharing structure and function [42]. Utilising single-cell RNA-sequencing data of an autochthonous murine model of UBC, Mmp12 mRNA expression was attributed to tumour-infiltrating macrophages. As observed in the human single-cell transcriptomic dataset, expression was also elevated in MIBC compared to NMIBC in male and female mice. However, Mmp7 mRNA transcripts in the murine data set were only detected in a minor subset of urothelial cells and its expression was inadequate to warrant further inspection.

MMPs degrade the tissue matrix by extracellular proteolysis, thus facilitating tumour cell migration, invasion and metastasis. Moreover, MMPs are involved in several physiological and tumour-supporting cellular processes, including loss of cell adhesion, tumour angiogenesis, cell proliferation, epithelial-to-mesenchymal transition and apoptosis [43]. They have been extensively analysed in UBC, which has revealed potential roles for certain MMPs as diagnostic markers and prognostic factors at different stages of the disease course [44]. MMP12 is linked to reduced overall survival rates in numerous cancers [45–47]; in UBC, certain genetic polymorphisms of the gene have been shown to increase invasiveness [48], while higher mRNA expression of MMP12 in tumour tissue has been reported to correlate with higher tumour grade [49].

Conclusion

The data presented here, and validated across two independent OLINK patient cohorts, represent the first report of MMP12 as a potential independent prognostic biomarker for UBC survival. Further, the information about the lymph node metastasis status does not exist in the Swedish training cohort, which is a limitation of our study. Another limitation is that the Swedish cohort does not include sex and age-matched control samples collected and stored at the same timepoint.

Future studies should focus on validating these findings across defined patient subgroups using liquid versus tissue based analyses, investigating the detailed mechanistic role of MMP12 on tumour progression and immune response regulation, and assessing the possibilities for its inhibition as a route to improving long-term survival in UBC patients. Our proteomic finding and the non-tumour origin of MMP12 also highlights the necessity that further studies should implement a holistic approach (meaning that profiling must consider all cells present in the tumour microenvironment) to avoid the risk of a selection bias for biomarkers that are tumour shed. This could complete the picture on how immune and stroma cells play a role in tumour invasiveness.

Supplementary Information

Additional file 1: Supplementary Table 1. Intersect of multiplex panel protein markers between the Swedish and the UK cohort.

Additional file 2: Supplementary Table 2. Univariate and multivariate Cox regression analyses of the VIMP proteins in the Swedish cohort and MMP12 in the UK cohort.

Additional file 3: Supplementary Fig. 1.

Additional file 4: Supplementary Fig. 2.

Additional file 5: Supplementary Fig. 3.

Additional file 6: Supplementary Fig. 4.

Additional file 7: Supplementary Fig. 5.

Abbreviations

AUC Area under curve

BCG Bacillus Calmette–Guérin

BCPP Bladder cancer prognosis programme

DAPI 4′,6-Diamidino-2-phenylindole

GEO Gene expression omnibus

LOD Limit of detection

MIBC Muscle invasive bladder cancer

MMP Matrix metalloprotease

NMIBC Non-muscle invasive bladder cancer

NPX Normalised protein expression

OH-BBN N-butyl-N-(4-hydroxybutyl)nitrosamine

OS Overall survival

PANCK Pan-cytokeratin

PEA Proximity extension assay

ROC Receiver operating characteristic

RSF Random survival forest

TMA Tissue microarray

UBC Urothelial bladder cancer

UMAP Uniform manifold and projection

UMI Unique molecular identifier

VIMP Variable importance

WHO World health organisation

Acknowledgements

We are very grateful for the reference pathology performed by pathologist Prof. Christer Busch.

Authors’ contributions

Concept and design: Lord, Kerzeli, Segersten, Mangsbo; Clinical data and sample acquisition: Segersten, Türker, Ward, Bryan; Analysis, or interpretation of data: All authors; Drafting of the manuscript: Lord, Kerzeli, Türker, Kostakis, Mangsbo; Critical revision of the manuscript for important intellectual content: All authors; Statistical analysis, data processing, illustrations: Lord, Kostakis, Kerzeli, Mezheyeuski, Segersten, Mangsbo; Obtained funding: Segersten, Mangsbo; Administrative, technical, or material support: All authors; Supervision: Lord, Segersten, Mangsbo;

Funding

Open access funding provided by Uppsala University. This work was supported by the Swedish Society for Medical Research (S15-0065) to S.M. and the Swedish Cancer Foundation (CAN 2017/199) to P-U. M. The funding bodies did not influence the research performed.

Availability of data and materials

OLINK data files are available from the corresponding author on reasonable request. Microarray data for the tumours from 308 UBC patients (in the form of non-normalised, raw intensity values along with their detection p-values) were retrieved from the GEO repository (GSE32894, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE32894); clinical information regarding the patients and their samples was obtained from the associated publication [20]. Two sets of single-cell RNA-sequencing data were downloaded and processed: unique molecular identifier (UMI) count matrices originating from bladders of a urothelial cancer mouse model [22] (GSE174182, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE174182) and raw sequence reads derived from patient tumour material from the European Nucleotide Archive under the accession code PRJNA662018 (https://www.ebi.ac.uk/ena/browser/view/PRJNA662018) [23].

All data generated or analysed during this study are included in this published article and its supplementary information files.

Declarations

Ethics approval and consent to participate

Regarding patient samples acquired in Sweden, this study was approved by Uppsala’s regional ethical review board (ID 2015/143/1). The need for informed consent was waived by the Uppsala’s regional ethical review board. Individual patient data was handled anonymised and the clinical dataset has not been made available. Regarding patient samples acquired in the United Kingdom, the study was approved by National Research Ethics Service Committee East Midlands—Derby (UK ethics ref: 06/MRE04/65) and written informed consent was obtained from all patients. The study was performed in accordance with the Declaration of Helsinki.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Iliana K. Kerzeli and Alexandros Kostakis shared first authors.

Martin Lord and Sara M. Mangsbo shared last authors.
==== Refs
References

1. Sung H Ferlay J Siegel RL Laversanne M Soerjomataram I Jemal A Global Cancer Statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 Countries CA Cancer J Clin 2021 71 3 209 249 10.3322/caac.21660 33538338
2. Lammers RJ Hendriks JC Rodriguez Faba OR Witjes WP Palou J Witjes JA Prediction model for recurrence probabilities after intravesical chemotherapy in patients with intermediate-risk non-muscle-invasive bladder cancer, including external validation World J Urol 2016 34 2 173 180 10.1007/s00345-015-1598-0 26025189
3. Sylvester RJ Rodríguez O Hernández V Turturica D Bauerová L Bruins HM European Association of Urology (EAU) Prognostic Factor Risk Groups for Non-muscle-invasive Bladder Cancer (NMIBC) Incorporating the WHO 2004/2016 and WHO 1973 classification systems for grade: an update from the EAU NMIBC guidelines panel Eur Urol 2021 79 4 480 488 10.1016/j.eururo.2020.12.033 33419683
4. Marsit CJ Karagas MR Danaee H Liu M Andrew A Schned A Carcinogen exposure and gene promoter hypermethylation in bladder cancer Carcinogenesis 2006 27 1 112 116 10.1093/carcin/bgi172 15987713
5. Jin F Thaiparambil J Donepudi SR Vantaku V Piyarathna DWB Maity S Tobacco-specific carcinogens induce hypermethylation, DNA adducts, and DNA damage in bladder cancer Cancer Prev Res (Phila) 2017 10 10 588 597 10.1158/1940-6207.CAPR-17-0198 28851690
6. Lopez-Beltran A Cimadamore A Blanca A Massari F Vau N Scarpelli M Immune checkpoint inhibitors for the treatment of bladder cancer Cancers (Basel) 2021 13 1 131 10.3390/cancers13010131 33401585
7. Song D Powles T Shi L Zhang L Ingersoll MA Lu YJ Bladder cancer, a unique model to understand cancer immunity and develop immunotherapy approaches J Pathol 2019 249 2 151 165 10.1002/path.5306 31102277
8. Samstein RM Lee CH Shoushtari AN Hellmann MD Shen R Janjigian YY Tumor mutational load predicts survival after immunotherapy across multiple cancer types Nat Genet 2019 51 2 202 206 10.1038/s41588-018-0312-8 30643254
9. Fu H Zhu Y Wang Y Liu Z Zhang J Xie H Identification and validation of stromal immunotype predict survival and benefit from adjuvant chemotherapy in patients with muscle-invasive bladder cancer Clin Cancer Res 2018 24 13 3069 3078 10.1158/1078-0432.CCR-17-2687 29514839
10. Hanada T Nakagawa M Emoto A Nomura T Nasu N Nomura Y Prognostic value of tumor-associated macrophage count in human bladder cancer Int J Urol 2000 7 7 263 269 10.1046/j.1442-2042.2000.00190.x 10910229
11. Takayama H Nishimura K Tsujimura A Nakai Y Nakayama M Aozasa K Increased infiltration of tumor associated macrophages is associated with poor prognosis of bladder carcinoma in situ after intravesical bacillus Calmette-Guerin instillation J Urol 2009 181 4 1894 1900 10.1016/j.juro.2008.11.090 19237175
12. Marchioni M Primiceri G Ingrosso M Filograna R Castellan P De Francesco P The clinical use of the Neutrophil to Lymphocyte Ratio (NLR) in urothelial cancer: a systematic review Clin Genitourin Cancer 2016 14 6 473 484 10.1016/j.clgc.2016.04.008 27209348
13. Yang G Shen W Zhang Y Liu M Zhang L Liu Q Accumulation of myeloid-derived suppressor cells (MDSCs) induced by low levels of IL-6 correlates with poor prognosis in bladder cancer Oncotarget 2017 8 24 38378 38388 10.18632/oncotarget.16386 28418913
14. Zhang H Ye YL Li MX Ye SB Huang WR Cai TT CXCL2/MIF-CXCR2 signaling promotes the recruitment of myeloid-derived suppressor cells and is correlated with prognosis in bladder cancer Oncogene 2017 36 15 2095 2104 10.1038/onc.2016.367 27721403
15. Bryan RT, Gordon NS, Abbotts B, Zeegers MP, Cheng KK, James ND, Ward DG. Multiplex screening of 422 candidate serum biomarkers in bladder cancer patients identifies syndecan-1 and macrophage colony-stimulating factor 1 as prognostic indicators. Transl Cancer Res. 2017;6(Suppl4):S657–65. 10.21037/tcr.2017.
16. Sauter G Eble J Epstein J Sesterhenn I World Health Organization classification of tumours. Pathology and genetics of tumours of the urinary system and male genital organs 2004 Lyon IARC 1 359
17. Zeegers MP Bryan RT Langford C Billingham L Murray P Deshmukh NS The West Midlands bladder cancer prognosis programme: rationale and design BJU Int 2010 105 6 784 788 10.1111/j.1464-410X.2009.08849.x 19751260
18. Hemdan T Lindén M Lind SB Namuduri AV Sjöstedt E de Ståhl TD The prognostic value and therapeutic target role of stathmin-1 in urinary bladder cancer Br J Cancer 2014 111 6 1180 1187 10.1038/bjc.2014.427 25072257
19. Lundberg M Eriksson A Tran B Assarsson E Fredriksson S Homogeneous antibody-based proximity extension assays provide sensitive and specific detection of low-abundant proteins in human blood Nucleic Acids Res 2011 39 15 e102 10.1093/nar/gkr424 21646338
20. Sjödahl G Lauss M Lövgren K Chebil G Gudjonsson S Veerla S A molecular taxonomy for urothelial carcinoma Clin Cancer Res 2012 18 12 3377 3386 10.1158/1078-0432.CCR-12-0077-T 22553347
21. Ritchie ME Phipson B Wu D Hu Y Law CW Shi W limma powers differential expression analyses for RNA-sequencing and microarray studies Nucleic Acids Res 2015 43 7 e47 10.1093/nar/gkv007 25605792
22. Kerzeli IK Lord M Doroszko M Elgendy R Chourlia A Stepanek I Single-cell RNAseq and longitudinal proteomic analysis of a novel semi-spontaneous urothelial cancer model reveals tumor cell heterogeneity and pretumoral urine protein alterations PLoS One 2021 16 7 e0253178 10.1371/journal.pone.0253178 34232958
23. Chen Z Zhou L Liu L Hou Y Xiong M Yang Y Single-cell RNA sequencing highlights the role of inflammatory cancer-associated fibroblasts in bladder urothelial carcinoma Nat Commun 2020 11 1 5077 10.1038/s41467-020-18916-5 33033240
24. Hao Y Hao S Andersen-Nissen E Mauck WM Zheng S Butler A Integrated analysis of multimodal single-cell data Cell 2021 184 13 3573 87.e29 10.1016/j.cell.2021.04.048 34062119
25. Hafemeister C Satija R Normalization and variance stabilization of single-cell RNA-seq data using regularized negative binomial regression Genome Biol 2019 20 1 296 10.1186/s13059-019-1874-1 31870423
26. Ishwaran H Kogalur UB Blackstone EH Lauer MS Random survival forests Ann Appl Stat 2008 2 3 841 60, 20 10.1214/08-AOAS169
27. Ritch CR Velasquez MC Kwon D Becerra MF Soodana-Prakash N Atluri VS Use and validation of the AUA/SUO risk grouping for nonmuscle invasive bladder cancer in a contemporary cohort J Urol 2020 203 3 505 511 10.1097/JU.0000000000000593 31609178
28. Ward DG Arnold R Bryan RT Molecular subtypes of T1 bladder cancer: biomolecular characteristics versus clinical utility Eur Urol 2020 78 4 538 539 10.1016/j.eururo.2020.07.021 32773351
29. Lindskrog SV Prip F Lamy P Taber A Groeneveld CS Birkenkamp-Demtröder K An integrated multi-omics analysis identifies prognostic molecular subtypes of non-muscle-invasive bladder cancer Nat Commun 2021 12 1 2301 10.1038/s41467-021-22465-w 33863885
30. Kamoun A de Reyniès A Allory Y Sjödahl G Robertson AG Seiler R A consensus molecular classification of muscle-invasive bladder cancer Eur Urol 2020 77 4 420 433 10.1016/j.eururo.2019.09.006 31563503
31. Xing S de Ferrari Andrade L NKG2D and MICA/B shedding: a ‘tag game’ between NK cells and malignant cells Clin Transl Immunol 2020 9 12 e1230 10.1002/cti2.1230
32. Lee GH Askari A Malietzis G Bernardo D Clark SK Knight SC The role of CD40 expression in dendritic cells in cancer biology; a systematic review Curr Cancer Drug Targets 2014 14 7 610 620 10.2174/1568009614666140828103253 25163469
33. Daneshmandi S Pourfathollah AA Forouzandeh-Moghaddam M Enhanced CD40 and ICOSL expression on dendritic cells surface improve anti-tumor immune responses; effectiveness of mRNA/chitosan nanoparticles Immunopharmacol Immunotoxicol 2018 40 5 375 386 10.1080/08923973.2018.1510959 30265161
34. Huang J Jochems C Anderson AM Talaie T Jales A Madan RA Soluble CD27-pool in humans may contribute to T cell activation and tumor immunity J Immunol 2013 190 12 6250 6258 10.4049/jimmunol.1300022 23677477
35. Bolenz C Knauf D John A Erben P Steidler A Schneider SW Decreased Invasion of Urothelial Carcinoma of the Bladder by Inhibition of Matrix-Metalloproteinase 7 Bladder Cancer 2018 4 1 67 75 10.3233/BLC-170124 29430508
36. Svatek RS Shah JB Xing J Chang D Lin J McConkey DJ A multiplexed, particle-based flow cytometric assay identified plasma matrix metalloproteinase-7 to be associated with cancer-related death among patients with bladder cancer Cancer 2010 116 19 4513 4519 10.1002/cncr.25401 20572047
37. Szarvas T Hoffmann MJ Olah C Szekely E Kiss A Hess J MMP-7 Serum and tissue levels are associated with poor survival in platinum-treated bladder cancer patients Diagnostics (Basel) 2020 11 1 48 10.3390/diagnostics11010048 33396213
38. Sharifi L Nowroozi MR Amini E Arami MK Ayati M Mohsenzadegan M A review on the role of M2 macrophages in bladder cancer; pathophysiology and targeting Int Immunopharmacol 2019 76 105880 10.1016/j.intimp.2019.105880 31522016
39. Lima L Oliveira D Tavares A Amaro T Cruz R Oliveira MJ The predominance of M2-polarized macrophages in the stroma of low-hypoxic bladder tumors is associated with BCG immunotherapy failure Urol Oncol 2014 32 4 449 457 10.1016/j.urolonc.2013.10.012 24252833
40. Koti M Ingersoll MA Gupta S Lam CM Li X Kamat AM Sex differences in bladder cancer immunobiology and outcomes: a collaborative review with implications for treatment Eur Urol Oncol 2020 3 5 622 630 10.1016/j.euo.2020.08.013 32967818
41. Rubin JB Lagas JS Broestl L Sponagel J Rockwell N Rhee G Sex differences in cancer mechanisms Biol Sex Differ 2020 11 1 17 10.1186/s13293-020-00291-x 32295632
42. Jackson BC Nebert DW Vasiliou V Update of human and mouse matrix metalloproteinase families Hum Genomics 2010 4 3 194 201 10.1186/1479-7364-4-3-194 20368140
43. Egeblad M Werb Z New functions for the matrix metalloproteinases in cancer progression Nat Rev Cancer 2002 2 3 161 174 10.1038/nrc745 11990853
44. Szarvas T vom Dorp F Ergün S Rübben H Matrix metalloproteinases and their clinical relevance in urinary bladder cancer Nat Rev Urol 2011 8 5 241 254 10.1038/nrurol.2011.44 21487384
45. Klupp F Neumann L Kahlert C Diers J Halama N Franz C Serum MMP7, MMP10 and MMP12 level as negative prognostic markers in colon cancer patients BMC Cancer 2016 16 494 10.1186/s12885-016-2515-7 27431388
46. Ng KT Qi X Kong KL Cheung BY Lo CM Poon RT Overexpression of matrix metalloproteinase-12 (MMP-12) correlates with poor prognosis of hepatocellular carcinoma Eur J Cancer 2011 47 15 2299 2305 10.1016/j.ejca.2011.05.032 21683576
47. McGowan PM Duffy MJ Matrix metalloproteinase expression and outcome in patients with breast cancer: analysis of a published database Ann Oncol 2008 19 9 1566 1572 10.1093/annonc/mdn180 18503039
48. Kader AK Liu J Shao L Dinney CP Lin J Wang Y Matrix metalloproteinase polymorphisms are associated with bladder cancer invasiveness Clin Cancer Res 2007 13 9 2614 2620 10.1158/1078-0432.CCR-06-1187 17473191
49. Wallard MJ Pennington CJ Veerakumarasivam A Burtt G Mills IG Warren A Comprehensive profiling and localisation of the matrix metalloproteinases in urothelial carcinoma Br J Cancer 2006 94 4 569 577 10.1038/sj.bjc.6602931 16465195