Next Article in Journal
Predictive Factors for Chemoradiation-Induced Oral Mucositis and Dysphagia in Head and Neck Cancer: A Scoping Review
Next Article in Special Issue
Cystoscopy Accuracy in Detecting Bladder Tumors: A Prospective Video-Confirmed Study
Previous Article in Journal
Analysis of Health-Related Quality of Life Reporting in Phase III RCTs of Advanced Genitourinary Tumors
Previous Article in Special Issue
MicroRNA-371a-3p—The Novel Serum Biomarker in Testicular Germ Cell Tumors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of SPP1 as a Prognostic Biomarker and Immune Cells Modulator in Urothelial Bladder Cancer: A Bioinformatics Analysis

by
Taoufik Nedjadi
1,*,
Mohamed Eldigire Ahmed
2,
Hifzur R. Ansari
1,
Sihem Aouabdi
1 and
Jaudah Al-Maghrabi
3,4
1
King Abdullah International Medical Research Centre, King Saud Bin Abdulaziz University for Health Sciences, Ministry of National Guard-Health Affairs, Jeddah 21423, Saudi Arabia
2
Faculty of Basic Sciences, King Saud Bin Abdulaziz University for Health Sciences, National Guard-Health Affairs, Jeddah 21423, Saudi Arabia
3
Department of Pathology, Faculty of Medicine, King Abdulaziz University, Jeddah 21589, Saudi Arabia
4
Department of Pathology, King Faisal Specialist Hospital & Research Center, Jeddah 21499, Saudi Arabia
*
Author to whom correspondence should be addressed.
Cancers 2023, 15(23), 5704; https://doi.org/10.3390/cancers15235704
Submission received: 9 October 2023 / Revised: 11 November 2023 / Accepted: 19 November 2023 / Published: 4 December 2023
(This article belongs to the Special Issue Biomarker in Urologic Cancer)

Abstract

:

Simple Summary

Bladder cancer (BC) is one of the most common urological cancers affecting more men than women. The absence of biomarkers for early detection and prognostication add to the burden of the disease. The SPP1 protein is differentially expressed in many malignancies and could serve as a potential prognostic biomarker. However, its expression profile, clinical significance and prognostic value, as well as its relationship with immune infiltration, have not been comprehensively studied in bladder cancer. Herewith, we aimed to investigate the prognostic value of SPP1 expression and its association with other cancer-promoting factors in bladder cancer. Our data indicated that SPP1 is differentially expressed at an early stage of BC development and is associated with an unfavorable prognosis. Interaction between SPP1 and other proteins/genes and immune cell infiltrates may play an important role in the pathogenesis and could as a potential therapeutic target for bladder cancer.

Abstract

Secreted phosphoprotein-1 (SPP1) expression is differentially altered in many malignancies and could serve as a potential prognostic biomarker. Recent findings indicated that SPP1 possesses a broader role in bladder cancer (BC) pathogenesis than previously envisioned; however, the underlying mechanisms governing its expression, cellular localization, prognostic value and immune-related role in bladder cancer remain poorly understood. The expression and the prognosis value of SPP1 were assessed using immunohistochemistry (IHC) staining on a tissue microarray. SPP1 expression was correlated with the clinicopathological parameters, and survival analysis was calculated using a Kaplan–Meier plotter. Bioinformatics analysis of TCGA data was queried using UALCAN, CIBERSORT and TIMER datasets to decipher the biological processes enrichment pattern, protein–protein interactions and characterize tumor-infiltrating immune cells, respectively. IHC revealed that SPP1 expression is significantly associated with tumor type, stage, grade and smoking status. The Kaplan–Meier survival curve showed that low SPP1 expression is an unfavorable prognostic indicator in bladder cancer patients (p = 0.02, log-rank). The significant increased expression of the SPP1 level is associated with evident hypomethylation of the gene promoter in cancer compared to normal tissues in the TCGA-bladder dataset. Missense mutation is the most frequent genetic alteration of the SPP1 gene. Protein–protein interactions demonstrated that SPP1 shares the same network with many important genes and is involved in many signaling pathways and biological processes. TIMER reported a significant correlation between SPP1 expression and multiple immune cells infiltration. Furthermore, the expression of SPP1 was found to be positively correlated with a number of immune checkpoint genes such as PD-1 and CTLA4. The current investigation indicates that the SPP1 protein could serve as a prognostic biomarker and merit further investigation to validate its clinical usefulness in patients with bladder cancer.

1. Introduction

Bladder cancer is the second most common urological malignancy after prostate cancer and represents a major cause of high mortality and comorbidity in men [1,2]. Recent reports revealed that the incidence of bladder cancer is increasing at an alarming pace in the developing countries especially in Egypt and Middle Eastern countries [3]. Pathologically, the disease is classified into non-muscle-invasive bladder cancer (NMIBC) and muscle-invasive bladder cancer (MIBC). The overall survival (OS) remains very poor especially for MIBC and metastatic patients with a 5-year OS of less than 50% with radical radiotherapy [4]. Also, the high rates of recurrence and progression to a muscle-invasive phenotype are the main obstacles that compromise the clinical management of the disease [5]. Currently, the assessment of prognosis relies enormously on the conventional staging and grading systems which lack enough accuracy to predict patients’ outcome [6]. This is mainly due to the large molecular heterogeneity of the disease [7]. Altogether, this situation highlights the complexity of the disease and necessitates the urgent need for additional prognostic biomarkers to improve patients’ outcome. The discovery of novel prognostic biomarkers is of paramount importance for high-risk patient stratification in order to identify those who need a more radical therapy at the outset. The molecular stratification would help in setting personalized management approaches and enhance the concept of precision oncology. Several attempts have been made to develop reliable screening biomarkers using urine, tissue or blood-based samples. However, most of them lack sufficient sensitivity and specificity for detecting cancer at an early stage and predicting outcomes [8,9]. Hence, it is crucial to dedicate additional efforts in order to identify novel biomarkers of great clinical significance to improve prognosis, optimize targeted therapy and alleviate the burden of the disease.
Secreted phosphoprotein-1 (SPP1), also known as osteopontin (OPN), is a secreted glycophosphoprotein and a member of the small integrin-binding ligand N-linked glycoproteins (SIBLING) that is encoded by the SPP1 gene localized on the long arm of chromosome 4q13 [10]. Biologically, secreted SPP1 has been associated with bone regeneration, cell adhesion, migration, invasion, chemotaxis and apoptosis [11]. The cellular functions of SPP1 have also been investigated in cancer and other disorders including obesity, diabetes and cardiovascular disease [12,13,14]. Aberrant over-expression of the SPP1 protein was reported in several solid tumors including prostate cancer [15], lung cancer [16], gastric cancer [17], breast cancer [18] and hepatocellular [19] and it is down-regulated in colorectal and endometrial carcinomas [20,21]. It has been reported that SPP1 regulates the anchorage-independent growth of tumor cells, invasion and chemo-resistance [22]. Interestingly, the expression of SPP1 is highly regulated in invasive cancer cells, indicating its intriguing role in cancer progression and distant metastasis [23]. Furthermore, recent findings reported that SPP1 could mediate the immune response and influence immunotherapy response in lung cancer through the modulation of macrophage M2 polarization and PD-L1 expression [24,25]. In the context of bladder cancer, the involvement of SPP1 in the pathogenesis of the disease is still poorly understood. Zaravinos et al. (2011) reported differential overexpression of the SPP1 gene transcript using DNA microarray analysis, while Ke et al. (2011) demonstrated an overexpression of the SPP1 protein in urothelial malignancy using immunohistochemistry [26,27]. Taken together, these findings suggested an important role played by SPP1 cancer development and progression and prompted us to investigate the protein expression patterns and the prognostic value, as well as the molecular pathways and biological function of SPP1. Additionally, we used a combination of online-available datasets and experimental data to comprehensively analyze SPP1-interacting proteins and its potential involvement in immune infiltration in bladder cancer.

2. Materials and Methods

2.1. Study Population and Sample Collection

A consecutive series of 182 bladder cancer patients, including 152 male patients (83.5%) and 30 female patients (16.5%) with a mean age of 61 years, were investigated. A tissue microarray (TMA) was constructed using formalin fixed paraffin-embedded (FFPE) blocks from consented bladder cancer patients who underwent bladder resection at King Abdulaziz University Hospital, Jeddah (KSA), between 2005 and 2010. None of the patients had received chemotherapy or radiotherapy prior to the surgical resection. The histological variants of bladder cancer include muscle-invasive bladder cancer (MIBC, 43.9%) and non-muscle-invasive bladder cancer (NMIBC, 39%). For correlation analysis, the clinical and pathological data of cancer patients, including age, sex, tumor grade, TNM stage smoking status and other parameters, were collected from patients’ medical records. The histopathological characteristics of the study cohort are described in Table 1a. The median patient age and follow-up was 62 (IQR = 17) and 14.4 months (IQR = 38), respectively (Table 1b).

2.2. TCGA Gene Expression Analysis

UALCAN, a web-based tool was used to evaluate the expression of SPP1 across various cancer types, as well as the methylation status of the SPP1 gene in bladder cancer. These data were extrapolated using UALCAN, a publicly available interactive online portal (http://ualcan.path.uab.edu/index.html accessed on 23 June 2023). The portal enables easy access to the mRNA expression of SPP1 in both cancer tissues and the matching normal counterparts based on The Cancer Genome Atlas (TCGA) database.

2.3. Tissue Microarray (TMA) Construction

Cylindrical tissue cores (1 mm in diameter) were extracted from FFPE donor blocks and arranged in duplicate in new recipient blocks. Tissue specimens were sectioned at a thickness of 5 μM, and the slides were then subjected to immunohistochemical staining.

2.4. Immunohistochemistry (IHC)

Immunohistochemical staining was performed using the Bench-Mark XT automated system (Ventana Medical Systems, Inc., Tucson, AZ, USA) as previously described [28]. Briefly, after a de-waxing step in xylene and a rehydration step, the FFPE tissue sections were then subjected to antigen retrieval using CC1 conditioning buffer for 60 min. The slides were then incubated with the primary polyclonal human anti-SPP1 antibody (Spring Bioscience, Pleasanton, CA, USA. cat#E3284) at a dilution of 1:100 for 32 min. Subsequently, the slides were treated with the DAB chromogen detection kit, counterstained with hematoxylin then with bluing reagent for 4 min as per the manufacturer’s instructions. After dehydration in graded ethanol and xylene, the sections were permanently mounted with DPX mounting solution and overlaid with coverslips.

2.5. Staining Evaluation and Scoring

SPP1 protein expression scoring was performed in a blind fashion to the patients’ clinicopathological data using a 40× objective. The information collected during the scoring process included the sub-cellular localization, the intensity of the staining on scale 0 to 3 (0 = negative; 1 = weak; 2 = moderate; and 3 = strong) and the extent of the staining (percentage of tumor cells showing positive immunoreactivity: 0–100% of cells). A staining index score was calculated by multiplying the staining intensity by the percentage of positive tumor cells [28]. The staining index varied from 0 to 300.

2.6. STRING Protein–Protein Interaction (PPI) Analysis

The search for known and predicted protein–protein interaction networks was performed using the STRING database [29]. These data provide evidence on both the functional and physical protein interactions of SPP1. The minimum required interaction score was set at the highest confidence interval (>0.900).

2.7. SPP1 Gene Co-Expression Analysis

The molecular alterations of SPP1 in 32 cancer studies were investigated using the TCGA Pan-Cancer Atlas through the cBioPortal cancer genomics website (https://www.cbioportal.org/ accessed on 9 June 2023). We used the expression data from cBioPortal entitled “blca_tcga_pub_2017” which was based on 412 muscle-invasive bladder cancer patients (PMID: 28988769). We used raw RSEM values without normalization. The gene–gene interaction network for humans was downloaded from the BIOGRID (version 4.4.209) database. The expression matrix and interaction network was used in CEMiTool, a Bioconductor package for gene co-expression module identification (PMID: 29458351) with default parameters and filter = TRUE, Pearson correlation and apply_vst (variance stabilizing transformation) = TRUE. We obtained 8 co-expression modules. SPP1 belongs to Module 3 with 468 genes (Supplementary Figure S1). The CEMiTool run gave ~15 most significant hub genes for Module 3. The gene set enrichment analysis was performed by using the clusterProfiler package (PMID: 22455463). All cluster network and KEGG pathway enrichment plots were generated in R.

2.8. Analysis of Bladder Cancer Infiltrating Immune Cells

We used the TIMER 2.0 database (PMID: 32442275) for tumor (BLCA)-infiltrating immune cells using SPP1 as the target gene. We selected only significant Spearman’s correlations (p < 0.05) with “Purity Adjustment” for bladder cancer.

2.9. Relashionship between SPP1 and Immune Checkpoint Genes in BC

The correlation between the SPP1 and immune checkpoint genes was assessed using the cBioPortal database on the TCGA bladder cancer data. Thirty-three immune checkpoint genes were screened and analyzed using the R package.

2.10. Statistical Analysis

The correlation between the SPP1 protein expression patterns and patients’ clinicopathological parameters was analyzed using the Chi-square test (χ2). For the association between SPP1 levels and overall patients’ survival, univariate survival analysis based on the Kaplan–Meier survival curve was used and the p-value was calculated using the log-rank test. The statistical analysis of the data was processed using SPSS (version 21) and data were considered significant for values of p < 0.05.

3. Results

3.1. SPP1 Profiling and Molecular Alterations in TCGA Dataset

We initially looked at the gene expression profile of SPP1 mRNA levels across a panel of multiple primary cancers and their matching normal tissues (24 in total) using the TCGA database. The data indicated no clear pattern of SPP1 expression, while the SPP1 mRNA level is over-expressed in the common cancers including breast, colon and liver, and we noticed a significant down-regulation of SPP1 expression in renal, pancreatic, cardiac and cartilage cancers compared to the normal tissues (Figure 1A). Compared with normal bladder tissues, bladder cancer showed an increased expression level of SPP1 concomitantly associated with decreased promoter methylation of the SPP1 gene (Figure 1B,C). The molecular changes in SPP1 were investigated using the cBioPortal for cancer genomics. The genetic alterations range from mutations and amplification to deep deletions. The data indicated that the SPP1 gene was altered in 2.19% of 411 bladder cases which were found to harbor the highest proportion of SPP1 gene amplifications with 1.22% relative to other cancer types (Figure 2A). The mutation frequency was also considerable (0.73%), most of which were missense mutations. These included D95Y and R248W, as well as an X73 splice (Figure 2B).

3.2. Expression of SPP1 in Bladder Cancer Patients

The TCGA data indicated that SPP1 may play an important role in bladder cancer pathogenesis and prompted us to investigate the expression pattern and analyze the prognostic value of SPP1 in bladder cancer. In order to do so, tissue microarrays were constructed from 182 Saudi bladder cancer patients and immunohistochemistry staining was performed using an antibody against SPP1. The IHC staining showed that our target protein exhibits both cytoplasmic and nuclear localization. Tumor cells showed diffuse cytoplasmic patterns of staining. The cytoplasmic staining intensity (Table 1c) varied from negative expression (12.8%) to weak expression (60.5%) to moderate expression (25%) to strong expression (1.7%) in some cores (Figure 3). On the other hand, the nuclear SPP1 staining was less prominent compared to the cytoplasmic staining where the expression pattern was observed as follows: 18.9% of the cores were negative, 73.3% exhibited weak staining and 7.8% moderate staining (Figure 4). For correlation analysis, the median intensity was taken as the cut-off point for low and high staining. Therefore, our data indicated that the cytoplasmic positivity of SPP1 staining was detected in 13% of all analyzed cores. For nuclear staining, 92% of the analyzed cohort were scored as low whereas only 8.0% of the specimens were considered as high. Interestingly, a positive association was observed between the staining of SPP1 in the cytoplasmic and nuclear compartments, meaning that patients who expressed high nuclear SPP1 had significantly higher cytoplasmic SPP1 scores (p < 0.001) (Supplementary Figure S2).

3.3. Relationship between SPP1 Expression and Clinicopathological Characteristics

Next, we sought to analyze the association between SPP1 expression and patients’ clinicopathological features. Our data indicated no association between the cytoplasmic SPP1 score and age, gender and marital status (p > 0.05). However, significant correlations were observed between the cytoplasmic level of SPP1 and the tumor grade (p = 0.046), stage (p = 0.018) and lymph node (p = 0.047) (Table 2). Interestingly, a signification positive correlation was apparent between SPP1 expression and smoking status (p = 0.002). Furthermore, non-muscle-invasive bladder cancer patients (NMIBC) showed a significant high expression of SPP1 (p = 0.019), supporting that increased expression of SPP1 could be an early event in bladder cancer development. Moreover, patients with high nuclear SPP1 expression have a low tumor grade (p = 0.007).

3.4. Correlation of SPP1 Expression and Survival Outcomes

The Kaplan–Meier survival curve demonstrated that the reduced expression of cytoplasmic SPP1 is significantly associated with poor patient outcomes (log-rank, p = 0.022). The median survival of patients exhibiting high cytoplasmic SPP1 staining was 29 months (95% confidence interval, 16.5–41.5 months), whereas the median survival of patients showing low cytoplasmic SPP1 staining was 35.2 months (95% confidence interval, 24.5–46 months) (Figure 5A). No survival difference was observed for bladder cancer patients when nuclear SPP1 expression was taken in consideration (Figure 5B).

3.5. Gene Enrichment and Functional Interactions of SPP1

We next sought to investigate the mechanism by which SPP1 could promote tumor development in bladder cancer. The search for known and predicted protein–protein interactions was performed using the cBioPortal database on TGCA data. A comprehensive quantitative network of SPP1 interactors indicated that our protein of interest belongs to the M3 node which embraces several interacting and co-expressing proteins (Figure 6A). Furthermore, a protein–protein interaction analysis obtained from the STRING website showed a direct interaction between SPP1 and CD44, MMP3, MMP7, TIMP1 and fibronectin-1, all of which are directly involved in the extracellular matrix remodeling and promotion of tumor metastasis [30] (Supplementary Figure S3A). The results revealed that SPP1 could bind to the afore-mentioned network of proteins with the highest confidence value (0.900). Additionally, SPP1 can also recognize albumin (ALB), integrin-A (ITGAV) and integrin-B (ITGB1) proteins whose expression was found to be highly increased at an early stage of bladder cancer development, using the TCGA-BLCA search engine (Supplementary Figure S3B–E). For SPP1-interacting proteins, the analysis of molecular function demonstrated that it is linked to several groups of proteins that are involved in chemokine receptor binding, chemokine activity, cytokine activity and immune receptor activity (Figure 6B). Moreover, KEGG analysis was performed in order to obtain a comprehensive understanding on the molecular pathway mechanisms, as well as the functional annotations associated with SPP1. KEGG gene enrichment analysis demonstrated that SPP1 is mainly enriched in cytokine receptor interaction, chemokine signaling pathway, hematopoietic cell lineage and the IL-17 and TLR signaling pathways (Figure 6C). Using the TIMER, we analyzed the correlation between SPP1 expression and immune cells’ infiltration. As shown in Figure 7, our data indicated that SPP1 expression was positively correlated with the infiltration of CD4+ memory-activated T cells, CD8+ T cells, neutrophils, dendritic cells and both M1 and M2 macrophages. Additionally, a significant inverse correlation was observed between SPP1 expression and B cells, endothelial cells and monocytes.

3.6. Correlation between SPP1 and Immune Checkpoint Genes

Next, we assessed the relationship between SPP1 expression and a list of immune checkpoint genes. Our data indicated that SPP1 expression exhibited a positive association with most of the immune checkpoint genes including PD-1, CTLA4, CD48, BTLA, TIGIT, LAG3 and LAIR1. However, an inverse correlation between SPP1 expression and a few immune checkpoint genes including LGALS9, ICOS, TNFSF1, TNFRSF1 and CD86 was also reported (Figure 8).

4. Discussion

Bladder carcinoma is a deadly disease especially when it spreads to the muscle layers of the bladder (MIBC) and the surrounding organs. Recent statistical data indicated that over half a million persons have developed bladder cancer and more than 200,000 patients died of the disease in 2020 [31]. The disease is characterized by its high recurrence rates which necessitate an active long-term surveillance/follow-up regimen [1]. The absence of effective molecular markers for early detection, disease progression and risk stratification renders the therapeutic management of the disease a real challenge. Therefore, new effective early diagnostic and/or prognostic biomarkers that can improve the clinical management of bladder cancer are urgently needed. Thus far, several potential biomarkers have been thoroughly studied [32] and others are currently under investigation, such as SPP1 particularly in bladder cancer. In the present study, we assessed the expression pattern of SPP1 and investigated whether its expression impacts patients’ survival or not. The correlation with the clinical and pathological features in a cohort of Saudi bladder cancer patients was also evaluated. We showed that SPP1 is differentially expressed between healthy controls and disease conditions where the pattern of expression was tissue dependent. The gene expression and DNA methylation of SPP1 from TCGA data demonstrated that the upregulation of SPP1 expression is associated with the decreased methylation of the same gene (Figure 1B,C). In this regard, many studies have linked methylation of the CpG island sites of the promoter region to gene silencing and to alternative gene splicing in cancer [33] which might be the case for the SPP1 gene in bladder cancer. DNA methylation has been suggested as a biomarker for high-accuracy detection and for outcome prediction in cancer [34]. Interestingly, computational analysis of the TCGA pan-cancer data characterized both gene expression and DNA methylation as the most effective predictors of the somatic mutation state [35], as previously demonstrated with the oncogenic mutations of the IDH1 and IDH2 genes in cancer [36]. The molecular analysis demonstrated that SPP1 gene amplification was the highest genetic alteration of SPP1 across many cancer types, whereas the frequency of SPP1 mutation (missense mutation) was low in bladder cancer compared to uterine, skin, lung, cervical, gastric and colorectal cancers. The impact of the SPP1 gene mutations and/or amplifications on bladder cancer pathogenesis remains to be further verified. It is well documented that bladder cancer is among the tumor types with the highest mutational burden involving several deletions of tumor suppressor genes and the amplification of oncogenes [37]. The amplifications of chromosome 6p22 and 11q are the most common alterations in bladder cancer [38]. Similarly, abnormal amplifications of the ERBB2, mdm2 and FGFR3 genes were found at a high frequency in bladder cancer [39,40,41]. Alterations of these genes appear to contribute, either directly or indirectly, to the development of bladder cancer and also to have clinical implications in the management of cancer patients.
Our IHC data revealed that high SPP1 expression was observed in around 13.0% of the study cohort, whereas the majority of the study population (87.0%) showed low or no SPP1 expression. Furthermore, patients with lower cytoplasmic SPP1 expression had a significantly reduced survival time compared with those who have higher cytoplasmic SPP1 (Figure 5A). In concordance with our data, reduced expression of SPP1 was reported as an unfavorable prognosticator in both colorectal cancer and endometrial carcinoma [20,21]. These data are also consistent with previous observations by Collins et al. (2012) who demonstrated that higher SPP1 expression is beneficial since it was associated with an improved outcome in patients with pancreatic cancer [42]. By contrast, several other investigations highlighted that high SPP1 levels were correlated with poor outcomes in many solid malignancies such as prostate cancer, lung cancer, gastric cancer and breast cancer [43,44,45,46]. For instance, Wong et al. (2017) showed that the overexpression of SPP1 in bladder cancer patients with high T-stage and tumor grades was significantly associated with decreased patient survival time [47]. The differential expression of SPP1 in the tissue and plasma samples [48] and its association with poor prognosis in cancer patients with an advanced disease stage indicated that SPP1 may serve as a valuable prognostic biomarker in different solid tumors [49,50] and a potential target for cancer therapy [27,51].
Interestingly, the prognostic value of nuclear SPP1 expression in the analyzed cohort could not be established (Figure 5B). However, low cytoplasmic expression of SPP1 correlated with muscle-invasive bladder cancer patients and high cytoplasmic SPP1 expression was predominantly found in patients with non-muscle-invasive bladder cancer (NMIBC) which represents an early stage in bladder cancer development (Table 2). Our results suggested for the first time that high levels of SPP1 may play a role at an early stage of bladder cancer development (NMIBC). Moreover, higher cytoplasmic SPP1 expression was significantly correlated with low-grade tumors and the early tumor stage (Ta) (Table 2). Our data demonstrated that among patients expressing high cytoplasmic SPP1 levels, 78% of them are Stage Ta and T1 and only 22% of the patients have a high tumor stage (T2 and T3). Interestingly, the current data reported the association between the smoking status and the expression level of the cytoplasmic SPP1 protein (p = 0.002). These data are in line with a recent study by Jiang et al. (2022) who showed that cigarette smoking increased the expression of OPN through the Jak2/Stat3 pathway in lung cancer [52]. Epidemiological studies demonstrated the strong link between tobacco cigarette smoking and the risk of bladder cancer development [53]. In a meta-analysis conducted on 15 cohort studies, Hou et al. (2017) indicated the strong association between smoking status and the prognosis of patients with bladder cancer [54].
These findings indicate that increased SPP1 expression might be an early event in bladder cancer development. A study by Hussain et al. (2017) reported the up-regulation of cytoplasmic SPP1 expression at a late stage of bladder cancer development [55]. Published reports indicated that cytoplasmic SPP1 is associated with a poor outcome and resistance to radiation via activation of the Jak2-Stat3 pathway. Biologically, SPP1 was reported to play a cancer-promoting role by regulating cell proliferation, motility, invasion and angiogenesis [56,57]. These data indicated that SPP1 could have a significant clinical impact, and patients with higher cytoplasmic levels should be managed differently and monitored carefully. Strikingly, nuclear SPP1 is poorly investigated and not much data are available on the functional role and the clinical implication of the altered expression of nuclear SPP1. The available data indicated that the nuclear isoform of SPP1 (OPN-c) can bind to polo-like kinase-1 and regulate apoptosis through binding to p53 [58]. Furthermore, Zduniak et al. (2015) reported that OPN-c is predominantly expressed in the nucleus in breast cancer and correlated with the tumor grade and a poor outcome in early breast cancer [59]. In our bladder cancer cohort, high nuclear SPP1 expression was seen in only 8.0% of the patients and was significantly associated with a low tumor grade (p = 0.007). In colorectal cancer, Assidi et al. (2019) reported a weak nuclear SPP1 expression in 23% of the analyzed cohort. Zduniak et al. (2015) used IHC staining to analyze the expression of SPP1-c which is a splice variant of the SPP1 gene, in a cohort of 671 breast cancer patients. The investigators revealed a strong staining intensity of nuclear SPP1-c which was significantly associated with poor outcomes in patients with early breast cancer [58]. Altogether, these data suggested that the intracellular localization (cytoplasmic and nuclear) of SPP1 or one of its splice variants merit further investigation, as it might have a clinical impact for cancer patients.
Thus far, the molecular mechanisms that underpin the link between SPP1 expression and bladder cancer progression and outcome remain largely unknown. However, previous studies indicated that SPP1 regulates cancer cells’ proliferation and motility. Using SPP1-targeted siRNA to knock down the levels of SPP1 expression in the T24 cell line, Xu et al. (2015) revealed that SPP1 plays an important role in promoting bladder cancer growth and invasiveness [60]. SPP1 has also been involved in cell adhesion, survival, angiogenesis and metastasis [61]. Recent findings established the link between SPP1 and the epithelial-to-mesenchymal transition (EMT) through the NRP2 protein pathway in bladder cancer cell lines [62]. The involvement of SPP1 in these biological and molecular functions are mainly mediated through binding of multiple cell surface receptors and the activation of various signal transduction pathways that have an impact on all hallmarks of carcinogenesis [63]. For instance, activation of the mitogen-activated protein kinase (MAPK) pathway, mTOR pathway and phosphatidylinositol 3-kinase/Akt signaling pathway have been implicated in OPN-associated cancer development [63]. Moreover, increased evidence supporting the connection between SSP1 and immune cells’ infiltration in promoting cancer progression has been established [64]. Immune checkpoint molecules such as PD-1, PD-L1 and CTLA4, which are targets for cancer immunotherapy, play an important role in maintaining self-tolerance and modulating the cancer immune response [65]. In this study, SPP1 demonstrated a positive correlation with several immune checkpoint genes including PD-1 and CTLA4, suggesting the engagement of SPP1 in carcinogenesis through the recruitment of multiple immune cells and activation of immune-related genes. This process is mediated via binding with integrins and CD44 receptors and activation of the Jak1/Stat1 signaling pathway [66,67]. Recent findings by Zheng et al. (2021) revealed that SPP1 expression was associated with EGFR mutation which conferred resistance to immunotherapy and promoted cancer progression. This pathway was characterized by low CD8+ T cells and high M2-type macrophages creating an immunosuppressed microenvironment [24]. In liver cancer, SPP1-positive macrophages had a specific interaction with the fibroblast-associated TME. This distinct communication between these cells promotes the immunosuppressive atmosphere in the TME and facilitates metastatic tumor progression [68]. SPP1-expressing macrophages have been demonstrated to influence the invasive potential of colorectal cancer cells via HLA-G [69]. Emerging data from Bill et al. (2023) found that the expression ratio of SPP1 and CXCL9 is critically important in defining the polarity of macrophages. This ratio has a strong prognostic value as it coordinates a network of pro-tumor variables and determines the antitumor immunity in the TME [70].
Taken together, the current study complements the previous reports on the major role of SPP1 in cancer pathogenesis. SPP1 might be clinically beneficial since its expression is significantly associated with the tumor grade, stage and survival of bladder cancer and other malignancies. Deepti et al. (2022) reported that the SPP1 mRNA level was associated with cancer staging and was a good discriminator between malignant and non-malignant tissues in cervical cancer [71]. This result was consolidated by recent findings indicating that increased SPP1 expression was associated with a poor outcome for patients with lung and breast cancers and promoted drug resistance and EMT transition in prostate cancer [72,73]. Our study sheds light on the prognostic value of SPP1 expression in bladder cancer. Additionally, we also revealed the potential interacting proteins and enrichment pathways that may orchestrate the tumor-promoting role of SPP1. We finally highlighted the association between SPP1 and immune infiltrating cells and immune checkpoint genes in bladder cancer. Given the limitations associated with our study, related mainly to the small sample size, incomplete clinical and pathological patients’ data and the heterogeneity of the disease, the present study highlights the importance of SPP1 as a valuable prognosticator and potential molecular target for bladder cancer therapy.

5. Conclusions

The current study sheds light on the importance and the functional significance of SPP1 in bladder cancer. Systematic analysis of data emerging from our cohort and the TCGA dataset revealed that SPP1 may play an important role in cancer development through interaction with several proteins, pathways and TME-associated cells. This work also highlights the importance of allocating additional priority to characterize SPP1 isoforms in cancer concomitantly with clinical research to validate the prognostic value of SPP1 and its potential therapeutic target in bladder cancer in a larger patient cohort and multi-institutional networks.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers15235704/s1, Figure S1: Module set enrichment. Figure S2: Box plot analysis. Figure S3: Functional interaction of SPP1 in bladder cancer.

Author Contributions

T.N. performed the study design, IHC data acquisition and scoring and manuscript drafting. M.E.A. contributed in statistical data analysis. H.R.A. helped with the bioinformatics data analysis. S.A. participated in study design, data interpretation and critically corrected the manuscript. J.A.-M. contributed in TMA design and IHC staining and scoring. All authors contributed in manuscript drafting and approved the final version. All authors have read and agreed to the published version of the manuscript.

Funding

The authors would like to acknowledge King Abdullah International Medical Research Centre, Kingdom of Saudi Arabia, for the financial support (protocol number NRJ22J/244 07). The funding body had no role in study design, data collection and analysis, interpretation of data or in the writing of the manuscript.

Institutional Review Board Statement

This study complied with the principles of the Declaration of Helsinki. The study was approved by the Ethics Committee of King Abdulaziz University Hospital, Jeddah, Saudi Arabia (protocol Ref#77-21). This study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Institutional Research Ethics Committee of King Abdulaziz University Hospital, Jeddah, Saudi Arabia (Ref#149-14).

Informed Consent Statement

Written informed consent was obtained from all participants in this study and both clinical and follow-up data were retrieved according to the permission and guidelines of the Ethics Committee.

Data Availability Statement

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kamat, A.M.; Hahn, N.M.; Efstathiou, J.A.; Lerner, S.P.; Malmström, P.U.; Choi, W.; Guo, C.C.; Lotan, Y.; Kassouf, W. Bladder cancer. Lancet 2016, 388, 2796–2810. [Google Scholar] [CrossRef] [PubMed]
  2. Lenis, A.T.; Lec, P.M.; Chamie, K.; Mshs, M.D. Bladder Cancer: A Review. JAMA 2020, 324, 1980–1991. [Google Scholar] [CrossRef] [PubMed]
  3. Antoni, S.; Ferlay, J.; Soerjomataram, I.; Znaor, A.; Jemal, A.; Bray, F. Bladder Cancer Incidence and Mortality: A Global Overview and Recent Trends. Eur. Urol. 2017, 71, 96–108. [Google Scholar] [CrossRef] [PubMed]
  4. Hussain, S.A.; Porta, N.; Hall, E.; Salawu, A.; Lewis, R.; Sreenivasan, T.; Wallace, J.; Crundwell, M.; Jenkins, P.; Tremlett, J.; et al. Outcomes in Patients with Muscle-invasive Bladder Cancer Treated with Neoadjuvant Chemotherapy Followed by (Chemo)radiotherapy in the BC2001 Trial. Eur. Urol. 2020, 79, 307–315. [Google Scholar] [CrossRef] [PubMed]
  5. Witjes, J.A.; Bruins, H.M.; Cathomas, R.; Compérat, E.M.; Cowan, N.C.; Gakis, G.; Hernández, V.; Espinós, E.L.; Lorch, A.; Neuzillet, Y.; et al. European Association of Urology Guidelines on Muscle-invasive and Metastatic Bladder Cancer: Summary of the 2020 Guidelines. Eur. Urol. 2021, 79, 82–104. [Google Scholar] [CrossRef]
  6. Jordan, B.; Meeks, J.J. T1 bladder cancer: Current considerations for diagnosis and management. Nat. Rev. Urol. 2019, 16, 23–34. [Google Scholar] [CrossRef]
  7. Kamat, A.M.; Hegarty, P.K.; Gee, J.R.; Clark, P.E.; Svatek, R.S.; Hegarty, N.; Shariat, S.F.; Xylinas, E.; Schmitz-Dräger, B.J.; Lotan, Y.; et al. ICUD-EAU International Consultation on Bladder Cancer 2012: Screening, diagnosis, and molecular markers. Eur. Urol. 2013, 63, 4–15. [Google Scholar] [CrossRef]
  8. Soria, F.; Krabbe, L.-M.; Todenhöfer, T.; Dobruch, J.; Mitra, A.P.; Inman, B.A.; Gust, K.M.; Lotan, Y.; Shariat, S.F. Molecular markers in bladder cancer. World J. Urol. 2019, 37, 31–40. [Google Scholar] [CrossRef]
  9. Cheng, L.; Davison, D.D.; Adams, J.; Lopez-Beltran, A.; Wang, L.; Montironi, R.; Zhang, S. Biomarkers in bladder cancer: Translational and clinical implications. Crit. Rev. Oncol. 2014, 89, 73–111. [Google Scholar] [CrossRef]
  10. Gimba, E.R.; Tilli, T.M. Human osteopontin splicing isoforms: Known roles, potential clinical applications and activated signaling pathways. Cancer Lett. 2013, 331, 11–17. [Google Scholar] [CrossRef]
  11. Shen, H.; Weber, G.F. The osteopontin-c splice junction is important for anchorage-independent growth. Mol. Carcinog. 2014, 53, 480–487. [Google Scholar] [CrossRef]
  12. Kahles, F.; Findeisen, H.M.; Bruemmer, D. Osteopontin: A novel regulator at the cross roads of inflammation, obesity and diabetes. Mol. Metab. 2014, 3, 384–393. [Google Scholar] [CrossRef] [PubMed]
  13. Maser, R.E.; Lenhard, M.J.; Pohlig, R.T.; Balagopal, P.B. Osteopontin and osteoprotegerin levels in type 2 diabetes and their association with cardiovascular autonomic function. J. Diabetes Its Complicat. 2016, 30, 507–510. [Google Scholar] [CrossRef] [PubMed]
  14. Mirzaei, A.; Mohammadi, S.; Ghaffari, S.H.; Yaghmaie, M.; Vaezi, M.; Alimoghaddam, K.; Ghavamzadeh, A. Osteopontin b and c Splice isoforms in Leukemias and Solid Tumors: Angiogenesis Alongside Chemoresistance. Asian Pac. J. Cancer Prev. APJCP 2018, 19, 615–623. [Google Scholar] [PubMed]
  15. Wisniewski, T.; Zyromska, A.; Makarewicz, R.; Zekanowska, E. Osteopontin And Angiogenic Factors As New Biomarkers Of Prostate Cancer. Urol. J. 2019, 16, 134–140. [Google Scholar] [PubMed]
  16. Shi, L.; Wang, X. Role of osteopontin in lung cancer evolution and heterogeneity. Semin. Cell Dev. Biol. 2017, 64, 40–47. [Google Scholar] [CrossRef] [PubMed]
  17. Di Bartolomeo, M.; Pietrantonio, F.; Pellegrinelli, A.; Martinetti, A.; Mariani, L.; Daidone, M.G.; Bajetta, E.; Pelosi, G.; de Braud, F.; Floriani, I.; et al. Osteopontin, E-cadherin, and β-catenin expression as prognostic biomarkers in patients with radically resected gastric cancer. Gastric Cancer 2016, 19, 412–420. [Google Scholar] [CrossRef]
  18. Gao, Y.L.; Xing, L.Q.; Ren, T.J.; Hou, J.F.; Xue, Q.; Liu, C.; Han, Y.M. The expression of osteopontin in breast cancer tissue and its relationship with p21ras and CD44V6 expression. Eur. J. Gynaecol. Oncol. 2016, 37, 41–47. [Google Scholar]
  19. Sun, T.; Li, P.; Sun, D.; Bu, Q.; Li, G. Prognostic value of osteopontin in patients with hepatocellular carcinoma: A systematic review and meta-analysis. Medicine 2018, 97, e12954. [Google Scholar] [CrossRef]
  20. Assidi, M.; Gomaa, W.; Jafri, M.; Hanbazazh, M.; Al-Ahwal, M.; Pushparaj, P.; Al-Harbi, A.; Al-Qahtani, M.; Buhmeida, A.; Al-Maghrabi, J. Prognostic value of Osteopontin (SPP1) in colorectal carcinoma requires a personalized molecular approach. Tumor Biol. 2019, 41, 1010428319863627. [Google Scholar] [CrossRef]
  21. Al Maghrabi, H.; Gomaa, W.; Al-Maghrabi, J. Increased osteopontin expression in endometrial carcinoma is associated with better survival outcome. Ginekol. Pol. 2020, 91, 73–78. [Google Scholar] [CrossRef]
  22. Chen, X.; Xiong, D.; Ye, L.; Yang, H.; Mei, S.; Wu, J.; Chen, S.; Mi, R. SPP1 inhibition improve cisplatin chemo-resistance of servical cancer cell lines. Cancer Chemother. Pharmacol. 2019, 83, 603–613. [Google Scholar] [CrossRef] [PubMed]
  23. Kovacheva, M.; Zepp, M.; Schraad, M.; Berger, S.; Berger, M.R. Conditional Knockdown of Osteopontin Inhibits Breast Cancer Skeletal Metastasis. Int. J. Mol. Sci. 2019, 20, 4918. [Google Scholar] [CrossRef] [PubMed]
  24. Zheng, Y.; Hao, S.; Xiang, C.; Han, Y.; Shang, Y.; Zhen, Q.; Zhao, Y.; Zhang, M.; Zhang, Y. The Correlation between SPP1 and Immune Escape of EGFR Mutant Lung Adenocarcinoma Was Explored by Bioinformatics Analysis. Front. Oncol. 2021, 11, 592854. [Google Scholar] [CrossRef] [PubMed]
  25. Zhang, Y.; Du, W.; Chen, Z.; Xiang, C. Upregulation of PD-L1 by SPP1 mediates macrophage polarization and facilitates immune escape in lung adenocarcinoma. Exp. Cell Res. 2017, 359, 449–457. [Google Scholar] [CrossRef]
  26. Ke, H.-L.; Chang, L.-L.; Yang, S.-F.; Lin, H.-H.; Li, C.-C.; Wu, D.-C.; Wu, W.-J. Osteopontin overexpression predicts poor prognosis of upper urinary tract urothelial carcinoma. Urol. Oncol. 2011, 29, 703–709. [Google Scholar] [CrossRef]
  27. Zaravinos, A.; Lambrou, G.I.; Volanis, D.; Delakas, D.; Spandidos, D.A. Spotlight on differentially expressed genes in urinary bladder cancer. PLoS ONE 2011, 6, e18255. [Google Scholar] [CrossRef]
  28. Buhmeida, A.; Dallol, A.; Merdad, A.; Al-Maghrabi, J.; Gari, M.A.; Abu-Elmagd, M.M.; Chaudhary, A.G.; Abuzenadah, A.M.; Nedjadi, T.; Ermiah, E.; et al. High fibroblast growth factor 19 (FGF19) expression predicts worse prognosis in invasive ductal carcinoma of breast. Tumor Biol. 2014, 35, 2817–2824. [Google Scholar] [CrossRef]
  29. Szklarczyk, D.; Gable, A.L.; Lyon, D.; Junge, A.; Wyder, S.; Huerta-Cepas, J.; Simonovic, M.; Doncheva, N.T.; Morris, J.H.; Bork, P.; et al. STRING v11: Protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019, 47, D607–D613. [Google Scholar] [CrossRef]
  30. Theocharis, A.D.; Manou, D.; Karamanos, N.K. The extracellular matrix as a multitasking player in disease. FEBS J. 2019, 286, 2830–2869. [Google Scholar] [CrossRef]
  31. Sung, H.; Ferlay, J.; Siegel, R.L.; Laversanne, M.; Soerjomataram, I.; Jemal, A.; Bray, F. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 2021, 71, 209–249. [Google Scholar] [CrossRef] [PubMed]
  32. Hong, M.; He, G.; Goh, S.; Low, A.W.X.; Tay, K.J.; Lim, T.K.H.; Yeong, J.; Khor, L.Y.; Lim, T.S. Biomarkers for Precision Urothelial Carcinoma Diagnosis: Current Approaches and the Application of Single-Cell Technologies. Cancers 2021, 13, 260. [Google Scholar] [CrossRef] [PubMed]
  33. Lin, S.; Yi, S.; Qiu, P. Comprehensive analysis of TCGA data reveals correlation between DNA methylation and alternative splicing. BMC Genom. 2022, 23, 758. [Google Scholar] [CrossRef]
  34. Hao, X.; Luo, H.; Krawczyk, M.; Wei, W.; Wang, W.; Wang, J.; Flagg, K.; Hou, J.; Zhang, H.; Yi, S.; et al. DNA methylation markers for diagnosis and prognosis of common cancers. Proc. Natl. Acad. Sci. USA 2017, 114, 7414–7419. [Google Scholar] [CrossRef]
  35. Crawford, J.; Christensen, B.C.; Chikina, M.; Greene, C.S. Widespread redundancy in-omics profiles of cancer mutation states. bioRxiv 2021. [Google Scholar] [CrossRef] [PubMed]
  36. Waitkus, M.S.; Yan, H. Targeting Isocitrate Dehydrogenase Mutations in Cancer: Emerging Evidence and Diverging Strategies. Clin. Cancer Res. 2021, 27, 383–388. [Google Scholar] [CrossRef] [PubMed]
  37. Kandoth, C.; McLellan, M.D.; Vandin, F.; Ye, K.; Niu, B.; Lu, C.; Xie, M.; Zhang, Q.; McMichael, J.F.; Wyczalkowski, M.A.; et al. Mutational landscape and significance across 12 major cancer types. Nature 2013, 502, 333–339. [Google Scholar] [CrossRef]
  38. Hurst, C.D.; Knowles, M.A. Mutational landscape of non-muscle-invasive bladder cancer. Urol. Oncol. 2022, 40, 295–303. [Google Scholar] [CrossRef]
  39. Nedjadi, T.; Al-Maghrabi, J.; Assidi, M.; Dallol, A.; Al-Kattabi, H.; Chaudhary, A.; Al-Sayyad, A.; Al-Ammari, A.; Abuzenadah, A.; Buhmeida, A.; et al. Prognostic value of HER2 status in bladder transitional cell carcinoma revealed by both IHC and BDISH techniques. BMC Cancer 2016, 16, 653. [Google Scholar] [CrossRef]
  40. Brunelli, M.; Tafuri, A.; Cima, L.; Cerruto, M.A.; Milella, M.; Zivi, A.; Buti, S.; Bersanelli, M.; Fornarini, G.; Vellone, V.G.; et al. MDM2 gene amplification as selection tool for innovative targeted approaches in PD-L1 positive or negative muscle-invasive urothelial bladder carcinoma. J. Clin. Pathol. 2022, 75, 39–44. [Google Scholar] [CrossRef]
  41. Tan, T.Z.; Rouanne, M.; Tan, K.T.; Huang, R.Y.; Thiery, J.P. Molecular Subtypes of Urothelial Bladder Cancer: Results from a Meta-cohort Analysis of 2411 Tumors. Eur. Urol. 2019, 75, 423–432. [Google Scholar] [CrossRef]
  42. Collins, A.L.; Rock, J.; Malhotra, L.; Frankel, W.L.; Bloomston, M. Osteopontin expression is associated with improved survival in patients with pancreatic adenocarcinoma. Ann. Surg. Oncol. 2012, 19, 2673–2678. [Google Scholar] [CrossRef] [PubMed]
  43. Yu, A.; Guo, K.; Qin, Q.; Xing, C.; Zu, X. Clinicopathological and prognostic significance of osteopontin expression in patients with prostate cancer: A systematic review and meta-analysis. Biosci. Rep. 2021, 41, BSR20203531. [Google Scholar] [CrossRef] [PubMed]
  44. Wang, J.; Shi, S.-M.; Su, Z.-B.; Yu, D.-J.; Tu, J.-W.; Zhu, J.Q.; Zhao, J.-P.; Sheng, L.; Wang, S.-B.; Sheng, Y.-J.; et al. Increased osteopontin protein expression may be correlated with poor prognosis in non-small-cell lung cancer: A meta analysis. J. Cancer Res. Ther. 2016, 12, 277–282. [Google Scholar] [CrossRef]
  45. Zhang, X.; Tsukamoto, T.; Mizoshita, T.; Ban, H.; Suzuki, H.; Toyoda, T.; Tatematsu, M. Expression of osteopontin and CDX2: Indications of phenotypes and prognosis in advanced gastric cancer. Oncol. Rep. 2009, 21, 609–613. [Google Scholar]
  46. Göthlin Eremo, A.; Lagergren, K.; Othman, L.; Montgomery, S.; Andersson, G.; Tina, E. Evaluation of SPP1/osteopontin expression as predictor of recurrence in tamoxifen treated breast cancer. Sci. Rep. 2020, 10, 1451. [Google Scholar] [CrossRef] [PubMed]
  47. Wong, J.P.; Wei, R.; Lyu, P.; Tong, O.L.; Zhang, S.D.; Wen, Q.; Yuen, H.F.; El-Tanani, M.; Kwok, H.F. Clinical and in vitro analysis of Osteopontin as a prognostic indicator and unveil its potential downstream targets in bladder cancer. Int. J. Biol. Sci. 2017, 13, 1373–1386. [Google Scholar] [CrossRef]
  48. Zhao, L.; Wang, Y.; Qu, N.; Huang, C.; Chen, L. Significance of plasma osteopontin levels in patients with bladder urothelial carcinomas. Mol. Diagn. Ther. 2012, 16, 311–316. [Google Scholar] [CrossRef]
  49. Zhao, H.; Chen, Q.; Alam, A.; Cui, J.; Suen, K.C.; Soo, A.P.; Eguchi, S.; Gu, J.; Ma, D. The role of osteopontin in the progression of solid organ tumor. Cell Death Dis. 2018, 9, 356. [Google Scholar] [CrossRef]
  50. Cai, X.; Zhang, H.; Li, T. The role of SPP1 as a prognostic biomarker and therapeutic target in head and neck squamous cell carcinoma. Int. J. Oral Maxillofac. Surg. 2022, 51, 732–741. [Google Scholar] [CrossRef]
  51. Wei, R.; Wong, J.P.C.; Kwok, H.F. Osteopontin -a promising biomarker for cancer therapy. J. Cancer 2017, 8, 2173–2183. [Google Scholar] [CrossRef] [PubMed]
  52. Jiang, Y.J.; Chao, C.C.; Chang, A.C.; Chen, P.C.; Cheng, F.J.; Liu, J.F.; Liu, P.I.; Huang, C.L.; Guo, J.H.; Huang, W.C.; et al. Cigarette smoke-promoted increases in osteopontin expression attract mesenchymal stem cell recruitment and facilitate lung cancer metastasis. J. Adv. Res. 2022, 41, 77–87. [Google Scholar] [CrossRef] [PubMed]
  53. Zeegers, M.P.; Tan, F.E.; Dorant, E.; van Den Brandt, P.A. The impact of characteristics of cigarette smoking on urinary tract cancer risk: A meta-analysis of epidemiologic studies. Cancer 2000, 89, 630–639. [Google Scholar] [CrossRef] [PubMed]
  54. Hou, L.; Hong, X.; Dai, M.; Chen, P.; Zhao, H.; Wei, Q.; Li, F.; Tan, W. Association of smoking status with prognosis in bladder cancer: A meta-analysis. Oncotarget 2017, 8, 1278–1289. [Google Scholar] [CrossRef] [PubMed]
  55. Hussain, S.A.; Palmer, D.H.; Syn, W.-K.; Sacco, J.J.; Greensmith, R.M.; Elmetwali, T.; Aachi, V.; Lloyd, B.H.; Jithesh, P.V.; Arrand, J.; et al. Gene expression profiling in bladder cancer identifies potential therapeutic targets. Int. J. Oncol. 2017, 50, 1147–1159. [Google Scholar] [CrossRef]
  56. Wang, M.; Sun, X.; Xin, H.; Wen, Z.; Cheng, Y. SPP1 promotes radiation resistance through JAK2/STAT3 pathway in esophageal carcinoma. Cancer Med. 2022, 11, 4526–4543. [Google Scholar] [CrossRef]
  57. Amilca-Seba, K.; Sabbah, M.; Larsen, A.K.; Denis, J.A. Osteopontin as a Regulator of Colorectal Cancer Progression and Its Clinical Applications. Cancers 2021, 13, 3793. [Google Scholar] [CrossRef]
  58. Junaid, A.; Moon, M.C.; Harding, G.E.; Zahradka, P. Osteopontin localizes to the nucleus of 293 cells and associates with polo-like kinase-1. Am. J. Physiol. Cell Physiol. 2007, 292, C919–C926. [Google Scholar] [CrossRef]
  59. Zduniak, K.; Ziolkowski, P.; Ahlin, C.; Agrawal, A.; Agrawal, S.; Blomqvist, C.; Fjällskog, M.-L.; Weber, G.F. Nuclear osteopontin-c is a prognostic breast cancer marker. Br. J. Cancer 2015, 112, 729–738. [Google Scholar] [CrossRef]
  60. Xu, S.-T.; Guo, C.; Ding, X.; Fan, W.-J.; Zhang, F.-H.; Xu, W.-L.; Ma, Y.-C. Role of osteopontin in the regulation of human bladder cancer proliferation and migration in T24 cells. Mol. Med. Rep. 2015, 11, 3701–3707. [Google Scholar] [CrossRef]
  61. Ortiz-Martínez, F.; Sanmartín, E.; Pomares-Navarro, E.; Pérez-Balaguer, A.; Andrés, L.; Sánchez-Payá, J.; Aranda, F.I.; Lerma, E.; Peiro, G. Osteopontin Regulates VEGFA and ICAM-1 mRNA Expression in Breast Carcinoma. Am. J. Clin. Pathol. 2015, 143, 812–822. [Google Scholar] [CrossRef] [PubMed]
  62. Schulz, A.; Gorodetska, I.; Behrendt, R.; Fuessel, S.; Erdmann, K.; Foerster, S.; Datta, K.; Mayr, T.; Dubrovska, A.; Muders, M.H. Linking NRP2 with EMT and Chemoradioresistance in Bladder Cancer. Front. Oncol. 2019, 9, 1461. [Google Scholar] [CrossRef] [PubMed]
  63. Shevde, L.A.; Samant, R.S. Role of osteopontin in the pathophysiology of cancer. Matrix Biol. 2014, 37, 131–141. [Google Scholar] [CrossRef]
  64. Moorman, H.R.; Poschel, D.; Klement, J.D.; Lu, C.; Redd, P.S.; Liu, K. Osteopontin: A Key Regulator of Tumor Progression and Immunomodulation. Cancers 2020, 12, 3379. [Google Scholar] [CrossRef] [PubMed]
  65. Topalian, S.L.; Drake, C.G.; Pardoll, D.M. Immune checkpoint blockade: A common denominator approach to cancer therapy. Cancer Cell 2015, 27, 450–461. [Google Scholar] [CrossRef] [PubMed]
  66. Lin, Y.H.; Yang-Yen, H.F. The osteopontin-CD44 survival signal involves activation of the phosphatidylinositol 3-kinase/Akt signaling pathway. J. Biol. Chem. 2001, 276, 46024–46030. [Google Scholar] [CrossRef] [PubMed]
  67. Zhang, N.; Li, F.; Gao, J.; Zhang, S.; Wang, Q. Osteopontin accelerates the development and metastasis of bladder cancer via activating JAK1/STAT1 pathway. Genes Genom. 2020, 42, 467–475. [Google Scholar] [CrossRef]
  68. Sathe, A.; Mason, K.; Grimes, S.M.; Zhou, Z.; Lau, B.T.; Bai, X.; Su, A.; Tan, X.; Lee, H.; Suarez, C.J.; et al. Colorectal Cancer Metastases in the Liver Establish Immunosuppressive Spatial Networking between Tumor-Associated SPP1+ Macrophages and Fibroblasts. Clin. Cancer Res. 2023, 29, 244–260. [Google Scholar] [CrossRef]
  69. Ozato, Y.; Kojima, Y.; Kobayashi, Y.; Hisamatsu, Y.; Toshima, T.; Yonemura, Y.; Masuda, T.; Kagawa, K.; Goto, Y.; Utou, M.; et al. Spatial and single-cell transcriptomics decipher the cellular environment containing HLA-G+ cancer cells and SPP1+ macrophages in colorectal cancer. Cell Rep. 2023, 42, 111929. [Google Scholar] [CrossRef]
  70. Bill, R.; Wirapati, P.; Messemaker, M.; Roh, W.; Zitti, B.; Duval, F.; Kiss, M.; Park, J.C.; Saal, T.M.; Hoelzl, J.; et al. CXCL9:SPP1 macrophage polarity identifies a network of cellular programs that control human cancers. Science 2023, 381, 515–524. [Google Scholar] [CrossRef]
  71. Deepti, P.; Pasha, A.; Kumbhakar, D.V.; Doneti, R.; Heena, S.K.; Bhanoth, S.; Poleboyina, P.K.; Yadala, R.S.D.A.; Pawar, S.C. Overexpression of Secreted Phosphoprotein 1 (SPP1) predicts poor survival in HPV positive cervical cancer. Gene 2022, 824, 146381. [Google Scholar] [CrossRef] [PubMed]
  72. Ji, X.; Liu, Y.; Mei, F.; Li, X.; Zhang, M.; Yao, B.; Wu, R.; You, J.; Pei, F. SPP1 overexpression is associated with poor outcomes in ALK fusion lung cancer patients without receiving targeted therapy. Sci. Rep. 2021, 11, 14031. [Google Scholar] [CrossRef] [PubMed]
  73. Pang, X.; Zhang, J.; He, X.; Gu, Y.; Qian, B.Z.; Xie, R.; Yu, W.; Zhang, X.; Li, T.; Shi, X.; et al. SPP1 Promotes Enzalutamide Resistance and Epithelial-Mesenchymal-Transition Activation in Castration-Resistant Prostate Cancer via PI3K/AKT and ERK1/2 Pathways. Oxidative Med. Cell. Longev. 2021, 2021, 5806602. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Expression pattern of SPP1 in pan-cancer. (A) The expression of the SPP1 mRNA level in multiple TCGA cancers and matching normal tissues. p < 0.001, except Thym cancer (p = 0.7) and KICH cancer (p = 0.53). (B) Increased mRNA expression level of SPP1 in bladder cancer. (C) Promoter methylation status of SPP1 in bladder cancer and matching normal tissues. All data were analyzed using the UALCAN web tool.
Figure 1. Expression pattern of SPP1 in pan-cancer. (A) The expression of the SPP1 mRNA level in multiple TCGA cancers and matching normal tissues. p < 0.001, except Thym cancer (p = 0.7) and KICH cancer (p = 0.53). (B) Increased mRNA expression level of SPP1 in bladder cancer. (C) Promoter methylation status of SPP1 in bladder cancer and matching normal tissues. All data were analyzed using the UALCAN web tool.
Cancers 15 05704 g001
Figure 2. Molecular alterations of SPP1 in cancers. (A) High amplifications/mutations of the SPP1 gene in bladder cancer compared other cancer types. (B) Positions and mutation frequency in SPP1 in bladder cancer.
Figure 2. Molecular alterations of SPP1 in cancers. (A) High amplifications/mutations of the SPP1 gene in bladder cancer compared other cancer types. (B) Positions and mutation frequency in SPP1 in bladder cancer.
Cancers 15 05704 g002
Figure 3. Cytoplasmic expression of SPP1 in bladder carcinoma. Immunohistochemical staining of the bladder cancer tissue microarray using an SPP1 antibody. Figures showing: no expression (A,B), weak (C,D), moderate (E,F) and strong expression (G,H) of SPP1. Images were taken using 10× and 40× magnification objectives (scale bar equals 1 mm).
Figure 3. Cytoplasmic expression of SPP1 in bladder carcinoma. Immunohistochemical staining of the bladder cancer tissue microarray using an SPP1 antibody. Figures showing: no expression (A,B), weak (C,D), moderate (E,F) and strong expression (G,H) of SPP1. Images were taken using 10× and 40× magnification objectives (scale bar equals 1 mm).
Cancers 15 05704 g003
Figure 4. Nuclear SPP1 expression in bladder carcinoma. Immunohistochemical staining of the bladder cancer tissue microarray using an SPP1 antibody. Figures showing: no expression (C,D), and strong expression of SPP1 (A,B). Images were taken with 10× and 40× magnification objectives (scale bar equals 1 mm).
Figure 4. Nuclear SPP1 expression in bladder carcinoma. Immunohistochemical staining of the bladder cancer tissue microarray using an SPP1 antibody. Figures showing: no expression (C,D), and strong expression of SPP1 (A,B). Images were taken with 10× and 40× magnification objectives (scale bar equals 1 mm).
Cancers 15 05704 g004
Figure 5. SPP1 expression and patients’ survival. Kaplan–Meier survival curve for bladder cancer patients expressing cytoplasmic (A) and nuclear (B) patterns of SPP1 (low expression vs. high expression). Low SPP1 immunostaining is associated with poor overall survival (log-rank p = 0.022).
Figure 5. SPP1 expression and patients’ survival. Kaplan–Meier survival curve for bladder cancer patients expressing cytoplasmic (A) and nuclear (B) patterns of SPP1 (low expression vs. high expression). Low SPP1 immunostaining is associated with poor overall survival (log-rank p = 0.022).
Cancers 15 05704 g005
Figure 6. Enrichment analysis of SPP1 in bladder cancer. (A) Identification of SPP1-interacting genes. Protein–protein interaction map and hub genes of SPP1. The size of the hub is proportional to the expression level. (B) Identification of the SPP1 co-expression network. The figure was generated using the online cBioPortal database. (C) KEGG functional enrichment analysis of SPP1. The figure was generated using the online cBioPortal database.
Figure 6. Enrichment analysis of SPP1 in bladder cancer. (A) Identification of SPP1-interacting genes. Protein–protein interaction map and hub genes of SPP1. The size of the hub is proportional to the expression level. (B) Identification of the SPP1 co-expression network. The figure was generated using the online cBioPortal database. (C) KEGG functional enrichment analysis of SPP1. The figure was generated using the online cBioPortal database.
Cancers 15 05704 g006
Figure 7. Correlation between immune cells and SPP1 expression. TIMER analysis of the correlation between SPP1 expression and immune cells’ infiltration. Purity-adjusted Spearman’s rho across various cell types by different algorithms.
Figure 7. Correlation between immune cells and SPP1 expression. TIMER analysis of the correlation between SPP1 expression and immune cells’ infiltration. Purity-adjusted Spearman’s rho across various cell types by different algorithms.
Cancers 15 05704 g007
Figure 8. Relationship between SPP1 expression and immune checkpoint genes in bladder cancer. (A) Correlation analysis between SPP1 expression and immune checkpoint genes. (B) The expression of immune checkpoint genes in relation to SPP1 expression. Data were analyzed using the cBioPortal cancer genomics website on TCGA data. The p-value significance codes: *** ≤0.001, ** ≤0.01, * ≤0.05.
Figure 8. Relationship between SPP1 expression and immune checkpoint genes in bladder cancer. (A) Correlation analysis between SPP1 expression and immune checkpoint genes. (B) The expression of immune checkpoint genes in relation to SPP1 expression. Data were analyzed using the cBioPortal cancer genomics website on TCGA data. The p-value significance codes: *** ≤0.001, ** ≤0.01, * ≤0.05.
Cancers 15 05704 g008
Table 1. (a): Baseline clinical and pathological characteristics of bladder cancer patients. (b): Median and IQR of patients’ age and follow-up duration. (c): Expression intensity of SPP1.
Table 1. (a): Baseline clinical and pathological characteristics of bladder cancer patients. (b): Median and IQR of patients’ age and follow-up duration. (c): Expression intensity of SPP1.
(a)
Variablen = 182%
Gender
  Male15283.5
  Female3016.5
Age
  <607641.8
  60+10658.2
Marital status
  Married17495.6
  Single84.4
Type of disease
  MIBC8043.9
  NMIBC7139.0
  Unknown3117.1
Grade
  Low grade9251.1
  High grade8848.9
Stage
  Ta4122.5
  15831.8
  24625.3
  3126.3
  42513.7
Lymph node
  None15886.8
  Positive2413.2
Smoking
  Yes4826.4
  No13473.6
Status
  Alive12468.5
  Dead5731.5
(b)
VariableMedianIQR
Age (Years)62.017.0
Follow-up duration (Months)14.438.0
(c)
Staining intensity (%)
Negative Weak Moderate Strong
Cytoplasmic 12.860.525 1.7
Nuclear 18.973.37.80
Table 2. Correlation between cytoplasmic SPP1 expression and patients’ clinicopathological characteristics.
Table 2. Correlation between cytoplasmic SPP1 expression and patients’ clinicopathological characteristics.
VariableLow
159 (87.4)
High
23 (12.6)
p Value
Gender
  Male130 (81.7)22 (95.6)0.093
  Female29 (18.3)1 (4.4)
Age
  <6064 (40.3)12 (52.2)0.278
  60+95 (59.7)11 (47.8)
Marital status
  Married151 (94.9)23 (100)0.271
  Single8 (5.1)0 (0)
Type of disease
  MIBC75 (47.2)5 (21.7)0.019 *
  NMIBC56 (35.2)15 (65.2)
  Undecided28 (17.6)3 (13.1)
Grade
  Low grade77 (48.4)16 (69.9)0.046 *
  High grade82 (51.6)7 (30.4)
Stage
  035 (22.0)6 (26.1)0.018 *
  146 (28.9)12 (52.2)
  243 (27.1)3 (13.0)
  310 (6.3)2 (8.7)
  425 (15.7)0 (0)
Lymph node
  None135 (84.9)23 (100)0.047 *
  Positive24 (15.1) 0 (0)
Smoking
  Yes34 (21.4)14 (60.9)0.002 **
  No125 (78.6)9 (39.1)
Status
  Alive106 (67.1)18 (78.3)0.281
  Dead52 (32.9)5 (21.7)
* p < 0.05; ** p < 0.01.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Nedjadi, T.; Ahmed, M.E.; Ansari, H.R.; Aouabdi, S.; Al-Maghrabi, J. Identification of SPP1 as a Prognostic Biomarker and Immune Cells Modulator in Urothelial Bladder Cancer: A Bioinformatics Analysis. Cancers 2023, 15, 5704. https://doi.org/10.3390/cancers15235704

AMA Style

Nedjadi T, Ahmed ME, Ansari HR, Aouabdi S, Al-Maghrabi J. Identification of SPP1 as a Prognostic Biomarker and Immune Cells Modulator in Urothelial Bladder Cancer: A Bioinformatics Analysis. Cancers. 2023; 15(23):5704. https://doi.org/10.3390/cancers15235704

Chicago/Turabian Style

Nedjadi, Taoufik, Mohamed Eldigire Ahmed, Hifzur R. Ansari, Sihem Aouabdi, and Jaudah Al-Maghrabi. 2023. "Identification of SPP1 as a Prognostic Biomarker and Immune Cells Modulator in Urothelial Bladder Cancer: A Bioinformatics Analysis" Cancers 15, no. 23: 5704. https://doi.org/10.3390/cancers15235704

APA Style

Nedjadi, T., Ahmed, M. E., Ansari, H. R., Aouabdi, S., & Al-Maghrabi, J. (2023). Identification of SPP1 as a Prognostic Biomarker and Immune Cells Modulator in Urothelial Bladder Cancer: A Bioinformatics Analysis. Cancers, 15(23), 5704. https://doi.org/10.3390/cancers15235704

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop