- Research article
- Open access
- Published:
Prognostic value of S1PR1 and its correlation with immune infiltrates in breast and lung cancers
BMC Cancer volume 20, Article number: 766 (2020)
Abstract
Background
Sphingosine-1-phosphate receptor (S1PR1) is involved in vascular development, a key process in tumorigenesis. This study aimed to evaluate its roles in tumor development and prognosis.
Methods
S1PR1 expression levels were analyzed using TIMER and Oncomine database, and the prognostic significance of S1PR1 was assessed using PrognoScan and Kaplan-Meier plotter databases. The relationship between S1PR1 and tumor-infiltrated immune cells was analyzed using TIMER.
Results
S1PR1 expression was remarkably lower in breast and lung cancer tissues than in the corresponding normal tissues. Lower expression was related to poor overall survival and disease-free survival in breast invasive carcinoma (BRCA), lung adenocarcinoma (LUAD), and lung squamous cell carcinoma (LUSC). A functional network analysis confirmed the function of S1PR1 in regulating vasculogenesis. In addition, S1PR1 levels were significantly negative with regard to the tumor purity of BRCA (r = − 0.508, P = 1.76e-66), LUAD (r = − 0.353, P = 6.05e-16), and LUSC (r = − 0.402, P = − 5.20e-20). Furthermore, S1PR1 levels were significantly positive with regard to infiltrating CD8+ (r = 0.38, P = 5.91e-35) and CD4+ T cells (r = 0.335, P = 1.03e-26), macrophages (r = 0.219, P = 3.67e-12), neutrophils (r = 0.168, P = 2.03e-7), and dendritic cells (DCs) (r = 0.208, P = 9.14e-11) in BRCA; S1PR1 levels were significantly positive with regard to CD8+ T cells (r = 0.308, P = 3.61e-12), macrophages (r = 0.376, P = 1.01e-17), neutrophils (r = 0.246, P = 4.15e-8), and DCs (r = 0.207, P = 4.16e-6) in LUAD; and positive with regard to B cells (r = 0.356, P = 1.57e-15), CD8+ (r = 0.459, P = 3.83e-26) and CD4+ T cells (r = 0.338, P = 3.98e-14), macrophages (r = 0.566, P = 2.61e-45), neutrophils (r = 0.453, P = 1.79e-25), and DCs (r = 0.56, P = 2.12e-40) in LUSC.
Conclusions
S1PR1 levels are positively correlated with multiple immune markers in breast and lung cancer. These observed correlations between S1PR1 and the prognosis and immune cell infiltration provide a foundation for further research on its immunomodulatory role in cancer.
Background
Sphingosine-1-phosphate (S1P), produced by sphingosine kinase (Sphk), is a biologically active signaling lipid [1]. S1P regulates vascular development and function, including vascular maturation [2, 3]. S1P receptor (S1PR1) is a biologically active sphingolipid metabolite that mediates S1P activity and promotes cell proliferation and survival [4, 5]. S1PR1 is widely expressed in vascular endothelial cells and is required for embryonic vascular development and maturation [6]. Estrogen (the growth-stimulating hormone in breast cancer cells) was shown to stimulate endothelial cell growth via S1PR1 [7, 8]. In the tumor microenvironment, S1P exhibits multiple functions: (a) it increases the survival of macrophages; (b) it serves as the “come-and-get-me” signal of dead cells, attracting and enhancing macrophage migration by combining with S1PR1; (c) it stimulates the polarization of TAM/M2 macrophages by activating S1PR1/2/4 [9,10,11]. Accumulating evidence demonstrated that tumor progression requires new blood vessel growth, which is achieved by producing angiogenic factors that can activate vascular endothelial cells [12]. Tumor cells release angiogenic stimuli, such as vascular endothelial growth factor (VEGF)-a, which leads to angiogenesis and tumor growth [13]. Studies have shown that S1PR1 inhibits VEGF signaling by promoting the interaction between VE-cadherin and VEGFR2, thereby inhibiting VEGF-induced vascular sprouting [14, 15].
However, the role of S1PR1 in tumorigenesis and its prognostic value are unclear. A preclinical study on human breast cancer cells found that S1PR1 antibody can enhance the cytotoxic and anti-proliferative effect of carboplatin on MDA-MB-231 and SK-BR-3 (HER2 subtype) cells, respectively [16]. Lei et al. found that S1PR1 signaling has tumor-suppressive effects and survival benefits in breast cancer [17]. Therefore, it is necessary to clarify the role of S1PR1 in tumor development and progression. Transcriptome analysis can be used to predict important issues, such as the intrinsic subtype of the primary tumor, tumor grade, drug reactivity, and recurrence risk [18,19,20].
Herein, we used Oncomine, Kaplan-Meier plotter, PrognoScan, UALCAN and GEPIA datasets to analyze S1PR1 expression and its relationship with the prognosis of cancer patients. Furthermore, we studied the correlation between S1PR1 and tumor-infiltrated immune cells in the tumor microenvironment using TIMER. Our results shed light on the important role of S1PR1 in breast and lung cancer, and determined that it is closely related to tumor immunity.
Methods
Oncomine database analysis
The Oncomine database (https://www.oncomine.org/resource/login.html) was used to evaluate the expression level of S1PR1 in various types of cancers [21]. The thresholds were a P-value of 0.0001, fold change of 2.0 and data type was mRNA.
PrognoScan database analysis
The PrognoScan database (www.prognoscan.org/) was used to test S1PR1 expression and survival in various types of cancers [22]. The threshold was an adjusted Cox P-value of < 0.05.
C-BioPortal database analysis
c-BioPortal (http://cbioportal.org) contains multidimensional cancer genomics data sets [23]. S1PR1 mutations and copy number variation (CNV) in breast and lung cancers were analyzed using c-BioPortal. The OncoPrint tab was used to obtain an overview of the genetic alterations for each sample.
Kaplan-Meier plotter
Kaplan-Meier Plotter (https://kmplot.com/) was applied to assess the prognostic value of S1PR1. Grouped according to the median expression of S1PR1 (high vs low expression), all patients were analyzed for overall survival (OS) and progression-free survival (PFS), and Kaplan-Meier was used to draw a survival chart [24].
Immune infiltrates analysis using the TIMER
TIMER 2.0 (https://cistrome.shinyapps.io/timer/) was used to analyze immune infiltrates across different types of cancer [25]. Especially, the expression of S1PR1 in different cancer types, and the correlation between the expression of S1PR1 and the abundance of immune invasion was determined. In addition, the correlation between S1PR1 expression and tumor infiltrating immune cell gene markers was also explored through related modules.
Gene correlation analysis using GEPIA
GEPIA (http://gepia.cancer-pku.cn/index.html) was used to confirm the genes with significantly correlated expression levels in TIMER [26]. The Spearman method was used to determine the correlation coefficients. The tumor tissue datasets were used for analysis.
LinkedOmics database analysis
The LinkedOmics database (http://www.linkedomics.org/login.php) was used to analyze S1PR1 co-expression based on Pearson’s correlation coefficients. The results were visually evaluated using volcano plots and heat maps. The function module of LinkedOmics was used to analyze gene ontology (GO) biological processes (BP) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways by a gene set enrichment analysis (GSEA). The rank criterion was FDR < 0.05 and 500 simulations were performed [27].
UALCAN database analysis
UALCAN (http://ualcan.path.uab.edu) included the Cancer Genome Atlas (TCGA) level RNA sequences. Clinical data from 31 cancer types were used to analyze the relative expression of genes in tumor and normal samples according to tumor stage, tumor grade or other clinicopathological characteristics [28].
S1PR1 mRNA expression level analysis
Gene expression data of breast invasive carcinoma (BRCA), lung adenocarcinoma (LUAD), and lung squamous cell carcinoma (LUSC) in TCGA were downloaded in UCSC Xena (https://xenabrowser.net). S1PR1 mRNA expression level was compared between cancerous and normal tissue using Mann-Whitney test with P < 0.05 setting as a cut-off.
Statistical analysis
Gene expression data in the Oncomine database was analyzed using p-value, fold change, and mRNA data type. The survival curves were generated via Kaplan-Meier plots and PrognoScan database are displayed with HR and P or Cox P-values from a log-rank test. Spearman correlation analysis was used to evaluate the correlation of gene expression in TIMER and LinkedOmics databases. P < 0.05 was considered statistically significant.
Results
S1PR1 mRNA expression levels in different types of human cancers
The Oncomine database was used to analyze S1PR1 mRNA levels in tumor tissues and normal tissues of various cancer types. S1PR1 expression was lower in most tumor tissues, including sarcoma, bladder, brain, central nervous system, breast, colorectal, leukemia, lung, myeloma, and ovarian cancer tissues, than in normal tissues (Fig. 1a). The mRNA-seq data from TCGA were analyzed using TIMER to verify these findings. Data from TCGA shown that the differential expression of S1PR1 between the tumor and adjacent normal tissues is shown in Fig. 1b. Compared with adjacent normal tissues, S1PR1 expression was significantly reduced in bladder urothelial carcinoma (BLCA), BRCA, cholangiocarcinoma (CHOL), colon adenocarcinoma (COAD), esophageal carcinoma (ESCA), head and neck squamous cell carcinoma (HNSC), kidney chromophobe (KICH), kidney renal papillary cell carcinoma (KIRP), liver hepatocellular carcinoma (LIHC), LUAD, LUSC, prostate adenocarcinoma (PRAD), rectum adenocarcinoma (READ), skin cutaneous melanoma (SKCM), stomach adenocarcinoma (STAD), and uterine corpus endometrial carcinoma (UCEC). However, S1PR1 expression was significantly higher in kidney renal clear cell carcinoma (KIRC) and thyroid carcinoma (THCA) than in adjacent normal tissues (Fig. 1b). These data showed that alterations in S1PR1 expression depend on the tumor type, suggesting that this gene exerts diverse functions in various tumors.
Prognostic evaluation of S1PR1 in cancers
We investigated whether S1PR1 expression is related to prognosis. The effect of S1PR1 expression on survival was evaluated by PrognoScan. Two probes (204642_at and 239401_at) matching S1PR1 were detected. Notably, in three breast cancer cohorts (GSE1456-GPL96, GSE7378, and GSE12276), low S1PR1 expression was significantly associated with a poorer prognosis breast cancer (Fig. 2a–f). We used the Kaplan-Meier plotter database to further examine the prognostic value of S1PR1 in breast cancer. Poor prognosis based on recurrence-free survival (RFS) in breast cancer was significantly correlated with low S1PR1 expression (HR = 0.67, P = 7.1e-13), but a significant correlation was not observed for overall survival (OS) (HR = 0.86, P = 0.17) and post-progression survival PPS (HR = 1.03, P = 0.82) (Fig. 2g–i). Its determined that the low expression of S1PR1 is an independent risk factor for poor prognosis of breast cancer. In addition, low S1PR1 expression was also related to poor prognosis in two cohorts of patients with lung cancer (GSE31210 and GSE8894), as determined using two probes (204642_at and 239401_at) (Fig. 2j–l). Kaplan-Meier plotter database also showed that low expression of S1PR1 was an independent risk factor for poor prognosis of lung cancer (overall survival, HR = 0.7, P = 6.9e-08; recurrence-free survival, HR = 0.71, P = 0.00035), but not related to post-progression survival in lung cancer (HR = 0.82, P = 0.14) (Fig. 2m–o).
Furthermore, we found that low S1PR1 expression was associated with a poor prognosis in patients with soft tissue, blood, and brain cancers (Fig. S1a–c). In contrast, low S1PR1 expression was an independent risk factor for a good prognosis in gastric cancer (Fig. S1d–g). These results confirmed the prognostic value of S1PR1 in specific types of cancer; both high and low S1PR1 expression was associated with prognosis depending on the type of cancer. Based on the consistent results for the associations between S1PR1 expression and survival in breast and lung cancers, we focused on the precise effects of S1PR1 in these two cancer types, as well as the underlying mechanisms.
Correlations between clinical characteristics and S1PR1 expression in breast cancer and lung cancer
We used the Kaplan-Meier plotter to study the relationship between S1PR1 expression and clinical characteristics in patients with breast cancer and lung cancer. Low expression of S1PR1 was associated with worse overall survival (OS) in male and female patients with lung adenocarcinoma (P < 0.05) (Table 1). In particular, low S1PR1 mRNA expression was correlated with worse OS in stage 1 (P = 9.20E-13) and early-stage (AJCC stage M) (P = 0.013) lung cancer (Table 1). Low S1PR1 mRNA expression was related to poor OS in patients with (P = 0.023) or without (P = 0.00075) a smoking history (Table 1). In addition, low S1PR1 mRNA expression was related to worse OS in patients who did not receive chemotherapy or radiotherapy. These findings strongly suggest that low S1PR1 mRNA expression is correlated with poor OS in lung cancer (Table 1). In BRCA, low S1PR1 mRNA expression was related to poor OS in ER-positive or HER2-negative patients and in the luminal androgen receptor subtype (Table 2). Taken together, high expression of S1PR1 could be considered a good prognostic indictor for breast and lung cancers depending on the clinical characteristics.
Decreased expression of S1PR1 in breast cancer and lung cancer patients
We further analyzed the expression of S1PR1 in breast and lung cancers. Gene expression data of breast invasive carcinoma (BRCA), lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC) in TCGA were downloaded and S1PR1 mRNA expression level was compared between tumor and normal tissue. As shown in Fig. 3a, the expression of S1PR1 was significantly decreased in tumor tissues of BRCA, LUAD and LUSC (Fig. 3a). In comparison with normal control tissues, breast cancer and lung cancer tissues presented lower expression of S1PR1, which was also observed by GEPIA analysis (Fig. 3b). Furthermore, we analyzed TCGA data using the UALCAN database. Compared to normal tissues, S1PR1 mRNA expression was significantly decreased in primary tumors and tumor stages (stage 1, stage 2, stage 3, and stage 4) of BRCA, LUAD, and LUSC (Fig. 3c–e). Taken together, these data confirmed the down-regulation of S1PR1 expression in breast cancer and lung cancer patients.
Regulators of S1PR1 in breast cancer and lung cancer
We used the LinkedOmics function module to detect the S1PR1 regulatory network to further understand the biological role of S1PR1 in breast cancer and lung cancer. Figure 4a–c shows genes with significantly positive (dark red dots) and negative (dark green dots) correlations with S1PR1 (false discovery rate, FDR < 0.01). The top 50 positively and negatively related genes are shown in a heat map in Fig. 4d–f. A Gene Ontology (GO)-based gene set enrichment analysis (GSEA) showed that genes that are co-expressed with S1PR1 are enriched for vasculogenesis and the purinergic receptor signaling pathway, while genes related to mitochondria and RNA transcript processing were inhibited in breast cancer (Fig. 4g). Similarly, GO annotation results showed that genes co-expressed with S1PR1 are primarily associated with vasculogenesis, the purinergic receptor signaling pathway, and the phospholipase C-activating G protein coupled receptor signaling pathway, while tRNA metabolic process, RNA modification, and RNA transcript processing were inhibited in lung cancer (Fig. 4h–i). A KEGG pathway analysis showed enrichment for hematopoietic cell lineage, Staphylococcus aureus infection, and renin secretion pathways in both breast cancer and lung cancer. Spliceosome, DNA replication, and proteasome pathways were inhibited in both tumor types (Fig. 4j-l). These results suggest that S1PR1 contributes to various processes in tumor development at least partially through regulate vasculogenesis.
Genomic alterations in S1PR1 in breast cancer and lung cancer
cBioPortal database was used to determine the types and frequencies of S1PR1 alterations in BRCA, LUAD, and LUSC. S1PR1 was altered in 4% of patients with BRCA. These alterations included mRNA missense mutations, amplifications, and deletions (Fig. 5a). S1PR1 was altered in 6% of patients with LUAD and 2.3% of patients with LUSC, including mRNA missense mutations, truncating mutations, amplifications, and deletions (Fig. 5a). Moreover, S1PR1 CNV was associated with OS in LUAD but not with OS or DFS in BRCA and LUSC (Fig. 5b–d). These results suggest that mutations in S1PR1 are associated with prognosis in LUAD.
Relationship between immune and S1PR1 expression in breast cancer and lung cancer
Tumor infiltrating lymphocytes (TIL) are lymphocytes that leave the blood circulation and migrate to the vicinity of the tumor. The amount of TIL in the tumor is an important indicator to predict the prognosis of cancer patients and the response to immunotherapy [29, 30]. Tumor purity is a key factor in analyses of immune infiltration by genomic approaches [31]. Therefore, we use TIMER to investigate whether the expression of S1PR1 in breast cancer and lung cancer is related to immune infiltration. We found a significant negative correlation between the S1PR1 expression level and tumor purity in both breast cancer and lung cancer (Fig. 6a–f, Left). S1PR1 was a determinant of immune infiltration in BRCA (tumor purity; r = − 0.508, P = 1.76e-66), including subtypes of BRCA (BRCA-Basal: r = − 0.5411, P = 1.28e-06; BRCA-Her2: r = − 0.505, P = 4.44e-06 and BRCA-Luminal: r = − 0.557, P = 9.15e-46). S1PR1 was related to immune infiltration in lung cancer, including LUAD (tumor purity; r = − 0.353, P = 6.05e-16) and LUSC (tumor purity; r = − 0.402, P = 5.20e-20).
Furthermore, the relationship between S1PR1 and specific immune infiltrates in breast cancer and lung cancer were analyzed. The S1PR1 expression level was significantly positively correlated with levels of infiltrating CD8+ T cells (r = 0.38, P = 5.97e-35), CD4+ T cells (r = 0.335, P = 1.03e-26), macrophages (r = 0.219, P = 3.67e-12), neutrophils (r = 0.168 P = 2.03e-07), and DCs (r = 0.208, P = 9.14e-11) in BRCA (Fig. 6a). In BRCA-Basal, there were slight positive correlations between S1PR1 expression levels and levels of infiltrating CD8+ T cells (r = 0.279, P = 1.76e-03) and CD4+ T cells (r = 0.237, P = 8.52e-03). Similarly, there were positive correlations with infiltrating levels of CD8+ T cells (r = 0.546, P = 1.13e-05), CD4+ T cells (r = 0.529, P = 2.00e-05), neutrophils (r = 0.342, P = 8.57e-03), and DCs (r = 0.488, P = 1.35e-04) in BRCA-Her2. S1PR1 expression levels were positively correlated with levels of infiltrating CD8+ T cells (r = 0.147, P = 3.43e-21), CD4+ T cells (r = 0.316, P = 6.26e-14), macrophages (r = 0.151, P = 4.14e-04), neutrophils (r = 0.147, P = 6.67e-04), and DCs (r = 0.213, P = 6.44e-07) in BRCA-Luminal tumors (Fig. 6a). We also found that S1PR1 expression levels were positively correlated with levels of infiltrating CD8+ T cells (r = 0.308, P = 3.61e-12), macrophages (r = 0.376, P = 1.01e-17), neutrophils (r = 0.246, P = 4.15e-08), and DCs (r = 0.207, P = 4.16e-06) in LUAD. In addition, there were positive correlations with levels of infiltrating B cells (r = 0.358, P = 1.27e-15), CD8+ T cells (r = 0.459, P = 3.83e-26), CD4+ T cells (r = 0.338, P = 3.98e-14), macrophages (r = 0.586, P = 2.61e-45), neutrophils (r = 0.453, P = 1.79e-25), and DCs (r = 0.56, P = 2.12e-40) in LUSC. These results strongly suggest that S1PR1 plays a special role in the immune infiltration of breast and lung cancers, and has a particularly strong effect on T cells, macrophages, neutrophils and DCs. These observed correlations between S1PR1 and various types of immune cells in breast and lung cancers indicated that S1PR1 may have high prognostic value.
Correlations between S1PR1 expression and immune markers
We further evaluated the correlations between S1PR1 and markers of various immune cells in breast cancer and lung cancer using TIMER (Table 3) and GEPIA databases (Table S1). The correlations between S1PR1 expression and immune marker genes for different immune cell populations, including CD8+ T cells, T cells (general), B cells, monocytes, TAMs, M1, and M2 macrophages, neutrophils, NK cells, DCs, and various functional T cells, such as Th1 cells, Th2 cells, Tfh cells, Th17 cells, and Tregs, as well as exhausted T cells were analyzed by TIMER. After adjusting for tumor purity, S1PR1 expression levels were significantly positively correlated with marker sets for various immune cells, except for NK cells, Th17, and T cell exhaustion in BRCA (Table 3 and Fig. 7). However, S1PR1 expression levels were highly positively correlated with most immune marker sets and both T cell populations and exhausted T cells in LUAD and LUSC (Table 3 and Fig. 7). We further analyzed the correlation between S1PR1 expression and the markers using the GEPIA database, including data for BRCA, LUAD, and LUSC. The results for correlations between S1PR1 and markers of immune infiltrating cells were similar to those of the TIMER analysis (Table S1). This further confirms that S1PR1 is significantly related to immune infiltrating cells in lung and breast cancer, suggesting that high levels of S1PR1 could induce immune activity in the lung and breast cancer microenvironment.
Discussion
We systematically analyzed the expression levels of S1PR1 and the prognostic value in different types of cancers. Compared with levels in normal tissues, S1PR1 expression was significantly lower in BLCA, BRCA, CHOL, COAD, ESCA, HNSC, KICH, KIRP, LIHC, LUAD, LUSC, PRAD, READ, SKCM, STAD, and UCEC and was significantly higher in KIRC and THCA. Accordingly, S1PR1 expression patterns depend on the type of cancer. Prognostic data from Kaplan-Meier plotter showed that low levels of S1PR1 are significantly related to poor prognosis in breast cancer and lung cancer.
The down-regulation of S1PR1 was associated with worse prognosis in breast cancer and lung cancer and was significantly related to clinical characteristics, such as gender, population, smoking status, and stage. These results suggested that S1PR1 is a prognostic biomarker in breast cancer and lung cancer. However, some literatures have reported the oncogenic role of S1PR1 in breast cancer. Lee H demonstrated that Stat3-induced S1PR1 expression, as well as S1P/S1PR1 pathway, is important for persistent Stat3 activation in cancer cells and the tumor microenvironment and for malignant progression [32]. This may be one of the molecular mechanisms by which S1PR1 mediates such a complex biological response. We considered that the main reason for this inconsistency is that our study analyzed the expression of S1PR1 at the overall level. We further verified the significant reduction of S1PR1 expression in breast cancer and lung cancer patients through TCGA analysis. Another study has also claimed a survival function benefit of S1P/S1PR signaling in BRCA patients, which might explain the obstacle to relative antagonist therapy in clinics [17]. A recent study determined that attenuated endothelial S1PR1 function led to increased tumor growth and metastasis, whereas S1PR1 overexpression led to smaller tumors, and strategies to enhance S1PR1 function in the tumor vasculature may potentiate the efficacy of cytotoxic and targeted anticancer therapies [33]. These studies support our findings that high expression of S1PR1 is beneficial for tumor survival.
The tumor microenvironment refers to non-cancer cells in and around tumors; infiltrated of immune cells in the tumor microenvironment plays a vital function in the occurrence and development of tumors [34, 35]. Lymphocyte infiltration in the tumor microenvironment is an independent predictor of cancer patient survival and lymph node metastasis [29, 30]. Studies have shown that S1PR1 can affect the proliferation and differentiation of lymphocytes in the tumor microenvironment [36]. The evaluation of immune cell infiltration in breast and lung cancers using the TIMER database revealed strong negative correlations between S1PR1 and tumor purity in BRCA, LUAD, and LUSC. Furthermore, the S1PR1 expression level was positively correlated with levels of CD8+ T, CD4+ T, neutrophils, macrophages, and DCs in BRCA. The correlation between S1PR1 expression and immune cell marker genes suggests that S1PR1 regulates lung cancer tumor immunity through multiple immune cell populations. These results indicate that high levels of S1PR1 could increase the cytotoxicity of the immune system and immune activation in BRCA, LUAD and LUSC by increasing the infiltration of CTLs, CD4 + T cells, and DCs. On the contrary, low expression of S1PR1 could lead to reduced infiltrated effector cells in the tumor microenvironment. As shown in recently reports, endothelial loss of S1PR1 led to a reduction in CD45+ cells, macrophages, and DCs, which influences tumor growth and metastasis [33]. In addition, S1P is involved in enhancing endocytosis and migration of mature dendritic cells through S1PR3, an event that may increase the immune response to cancer cells. Our findings are consistent with such reports, and these discoveries imply that S1PR1 plays an important role in recruiting and governing immune infiltration in BRCA, LUAD and LUSC.
To further elucidate the molecular mechanisms underlying the role of S1PR1 in breast and lung cancers, we used GSEA to identify pathways that are enriched in genes co-expressed with S1PR1. We found that S1PR1 was significantly associated with vasculogenesis, the purinergic receptor signaling pathway, and metabolism of nucleic acids in tumor conditions. This conclusion is consistent with previous research reports that showed that S1PR1 regulates vasculogenesis [7]. Recent studies have provided potential explanations for the associations between S1PR1 expression, immune infiltration, and poor prognosis. Angiogenesis mimicry (VM) system is a blood vessel-like network in which tumor cells are co-expressed with endothelial cells and tumor markers [37]. VM is closely related to a variety of human malignancies, including breast cancer [38]. Angiogenesis mimicry leads to worse prognosis, increased tumor metastasis, low 5-year overall survival, and increased mortality [39]. This shows that S1PR1 defects promote the occurrence of VM, and the knockout of S1PR1 in breast cancer cells increases the number of VMs. More importantly, tumor cells with low S1PR1 expression receive nutrition through VM, and accelerate tumor growth in animal models [40]. Recent research has shown that S1PR1 signaling is an important vascular factor affecting tumor progression, metastasis, and responses to chemotherapy and immunotherapy [33]. Strategies to enhance S1PR1 function in the tumor vasculature may enhance the cytotoxic killing effect and chemotherapy effect of targeted anti-cancer therapy.
A limitation of our study was the lack of in vitro and animal experiments to confirm the role of S1PR1 in the growth and progression of breast cancer and lung cancer and its relationship with the infiltration of immune cells in the tumor microenvironment. Therefore, further research is needed to verify the role of S1PR1 in breast cancer and lung cancer using these models.
Conclusions
In conclusion, decreased S1PR1 expression was related to poor prognosis together with reduction of effect immune cell infiltration in breast and lung cancers. In addition, the down-regulation of S1PR1 may have profound effects on vasculogenic mimicry in tumor microenvironment. Accordingly, S1PR1 presumably plays a critical part in immune infiltration and acts as a prognostic marker in breast cancer and lung cancer.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Abbreviations
- S1P:
-
Sphingosine-1-phosphate
- S1PR1:
-
Sphingosine-1-phosphate receptor
- CNV:
-
Copy number variation
- GO:
-
Gene Ontology
- BP:
-
Biological processes
- KEGG:
-
Kyoto Encyclopedia of Genes and Genomes
- GSEA:
-
Gene set enrichment analysis
- TCGA:
-
The Cancer Genome Atlas
- TIL:
-
Tumor infiltrating lymphocytes
- BLCA:
-
Bladder urothelial carcinoma
- BRCA :
-
Breast invasive carcinoma
- CHOL:
-
Cholangiocarcinoma
- COAD:
-
Colon adenocarcinoma
- ESCA:
-
Esophageal carcinoma
- HNSC:
-
Head and neck squamous cell carcinoma
- KICH:
-
Kidney chromophobe
- KIRP:
-
Kidney renal papillary cell carcinoma
- LIHC:
-
Liver hepatocellular carcinoma
- LUAD:
-
Lung adenocarcinoma
- LUSC:
-
Lung squamous cell carcinoma
- PRAD:
-
Prostate adenocarcinoma
- READ:
-
Rectum adenocarcinoma
- SKCM:
-
Skin cutaneous melanoma
- STAD:
-
Stomach adenocarcinoma
- UCEC:
-
Uterine corpus endometrial carcinoma
- KIRC:
-
Kidney renal clear cell carcinoma
- THCA:
-
Thyroid carcinoma
- OS:
-
Overall survival
- RFS:
-
Relapse-Free Survival
- PPS:
-
Post-Progression survival
- DSS:
-
Disease-specific survival
- DFS:
-
Disease-free survival
References
Sun K, Zhang Y, D'Alessandro A, Nemkov T, Song A, Wu H, et al. Sphingosine-1-phosphate promotes erythrocyte glycolysis and oxygen release for adaptation to high-altitude hypoxia. Nat Commun. 2016;15(7):12086.
Asano Y, Stawski L, Hant F, Highland K, Silver R, Szalai G, et al. Endothelial Fli1 deficiency impairs vascular homeostasis: a role in scleroderma vasculopathy. Am J Pathol. 2010;176(4):1983–98.
Allende ML, Yamashita T, Proia RL. G-protein-coupled receptor S1P1 acts within endothelial cells to regulate vascular maturation. Blood. 2003;102(10):3665–7.
Spiegel S, Milstien S. The outs and the ins of sphingosine-1-phosphate in immunity. Nat Rev Immunol. 2011;11(6):403–15.
Spiegel S, Milstien S. Sphingosine 1-phosphate, a key cell signaling molecule. J Biol Chem. 2002;277(29):25851–4.
Proia RL, Hla T. Emerging biology of sphingosine-1-phosphate: its role in pathogenesis and therapy. J Clin Invest. 2015;125(4):1379–87.
Sukocheva O, Wadham C, Gamble J, Xia P. Sphingosine-1-phosphate receptor 1 transmits estrogens' effects in endothelial cells. Steroids. 2015;104:237–45.
Sukocheva OA. Expansion of Sphingosine kinase and Sphingosine-1-phosphate receptor function in Normal and Cancer cells: from membrane restructuring to mediation of estrogen signaling and stem cell programming. Int J Mol Sci. 2018;19(2):420.
Rodriguez YI, Campos LE, Castro MG, Aladhami A, Oskeritzian CA, Alvarez SE. Sphingosine-1 phosphate: a new modulator of immune plasticity in the tumor microenvironment. Front Oncol. 2016;6:218.
Nakajima M, Nagahashi M, Rashid OM, Takabe K, Wakai T. The role of sphingosine-1-phosphate in the tumor microenvironment and its clinical implications. Tumour Biol. 2017;39(4):1010428317699133.
Sukocheva OA, Furuya H, Ng ML, Friedemann M, Menschikowski M, Tarasov VV, et al. Sphingosine kinase and sphingosine-1-phosphate receptor signaling pathway in inflammatory gastrointestinal disease and cancers: a novel therapeutic target. Pharmacol Ther. 2020;207:107464.
Folkman J. Angiogenesis: an organizing principle for drug discovery? Nat Rev Drug Discov. 2007;6(4):273–86.
Ferrara N. Pathways mediating VEGF-independent tumor angiogenesis. Cytokine Growth Factor Rev. 2010;21(1):21–6.
Jung B, Obinata H, Galvani S, Mendelson K, Ding B, Skoura A, et al. Flow-regulated endothelial S1P receptor-1 signaling sustains vascular development. Dev Cell. 2012;23(3):600–10.
Gaengel K, Niaudet C, Hagikura K, Hagikura K, Laviña B, Muhl L, et al. The sphingosine-1-phosphate receptor S1PR1 restricts sprouting angiogenesis by regulating the interplay between VE-cadherin and VEGFR2. Dev Cell. 2012;23(3):587–99.
Xiao S, Yang J. Preclinical study of the antitumor effect of sphingosine-1-phosphate receptor 1 antibody (S1PR1-antibody) against human breast cancer cells. Investig New Drugs. 2019;37(1):57–64.
Lei FJ, Cheng BH, Liao PY, et al. Survival benefit of sphingosin-1-phosphate and receptors expressions in breast cancer patients. Cancer Med. 2018;7(8):3743–54.
Glinsky GV, Glinskii AB, Stephenson AJ, Hoffman RM, Gerald WL. Gene expression profiling predicts clinical outcome of prostate cancer. J Clin Invest. 2004;113(6):913–23.
Chen B, Tang H, Chen X, Zhang G, Wang Y, Xie X, et al. Transcriptomic analyses identify key differentially expressed genes and clinical outcomes between triple-negative and non-triple-negative breast cancer. Cancer Manag Res. 2018;11:179–90 Published 2018 Dec 21.
Tang H, Huang X, Wang J, Yang L, Kong Y, Gao G, et al. circKIF4A acts as a prognostic factor and mediator to regulate the progression of triple-negative breast cancer. Mol Cancer. 2019;18(1):23.
Rhodes DR, Kalyana-Sundaram S, Mahavisno V, Varambally R, Yu J, Briggs BB, et al. Oncomine 3.0: genes, pathways, and networks in a collection of 18,000 cancer gene expression profiles. Neoplasia. 2007;9(2):166–80.
Lánczky A, Nagy Á, Bottai G, Munkácsy G, Szabó A, Santarpia L, et al. miRpower: a web-tool to validate survival-associated miRNAs utilizing expression data from 2178 breast cancer patients. Breast Cancer Res Treat. 2016;160(3):439–46.
Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO, et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal. 2013;6(269):pl1.
Szász AM, Lánczky A, Nagy Á, Förster S, Hark K, Green JE, et al. Cross-validation of survival associated biomarkers in gastric cancer using transcriptomic data of 1,065 patients. Oncotarget. 2016;7(31):49322–33.
Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, et al. TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. 2017;77(21):e108–10.
Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45(W1):W98–W102.
Vasaikar SV, Straub P, Wang J, Zhang B. LinkedOmics: analyzing multi-omics data within and across 32 cancer types. Nucleic Acids Res. 2018;46(D1):D956–63.
Chandrashekar DS, Bashel B, Balasubramanya SAH, Creighton CJ, Ponce-Rodriguez I, Chakravarthi BVSK, et al. UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses. Neoplasia. 2017;19(8):649–58.
Ohtani H. Focus on TILs: prognostic significance of tumor infiltrating lymphocytes in human colorectal cancer. Cancer Immun. 2007;7:4.
Azimi F, Scolyer RA, Rumcheva P, Moncrieff M, Murali R, McCarthy SW, et al. Tumor-infiltrating lymphocyte grade is an independent predictor of sentinel lymph node status and survival in patients with cutaneous melanoma. J Clin Oncol. 2012;30(21):2678–83.
Yoshihara K, Shahmoradgoli M, Martínez E, Vegesna R, Kim H, Torres-Garcia W, et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun. 2013;4:2612.
Lee H, Deng J, Kujawski M, Yang C, Liu Y, Herrmann A, et al. STAT3-induced S1PR1 expression is crucial for persistent STAT3 activation in tumors. Nat Med. 2010;16(12):1421–8.
Cartier A, Leigh T, Liu CH, Hla T. Endothelial sphingosine 1-phosphate receptors promote vascular normalization and antitumor therapy. Proc Natl Acad Sci U S A. 2020;117(6):3157–66.
Zhu J, Petit PF, Van den Eynde BJ. Apoptosis of tumor-infiltrating T lymphocytes: a new immune checkpoint mechanism. Cancer Immunol Immunother. 2019;68(5):835–47.
Sharma P, Hu-Lieskovan S, Wargo JA, Ribas A. Primary, adaptive, and acquired resistance to Cancer immunotherapy. Cell. 2017;168(4):707–23.
Aoki M, Aoki H, Ramanathan R, Hait NC, Takabe K. Corrigendum to "Sphingosine-1-phosphate signaling in immune cells and inflammation: roles and therapeutic potential". Mediat Inflamm. 2016;2016:2856829.
Valdivia A, Mingo G, Aldana V, Pinto MP, Ramirez M, Retamal C, et al. Fact or fiction, it is time for a verdict on Vasculogenic mimicry? Front Oncol. 2019;9:680.
El Hallani S, Boisselier B, Peglion F, Rousseau A, Colin C, Idbaih A, et al. A new alternative mechanism in glioblastoma vascularization: tubular vasculogenic mimicry. Brain. 2010;133(Pt 4):973–82. https://doi.org/10.1093/brain/awq044.
Cao Z, Bao M, Miele L, Sarkar FH, Wang Z, Zhou Q. Tumour vasculogenic mimicry is associated with poor prognosis of human cancer patients: a systemic review and meta-analysis. Eur J Cancer. 2013;49(18):3914–23.
Liu S, Ni C, Zhang D, Sun H, Dong X, Che N, et al. S1PR1 regulates the switch of two angiogenic modes by VE-cadherin phosphorylation in breast cancer. Cell Death Dis. 2019;10(3):200.
Acknowledgements
We thank Xiaotao Jiang for his help in data analysis.
Funding
YFL was supported by National Natural Science Foundation of China (No. 81700512), and Natural Science Foundation of Guangdong Province of China (No. 2016A030310252).
Author information
Authors and Affiliations
Contributions
LMZ, DLC, and YFL designed the study. LMZ, LLX and ZYY collected and analyzed the data. LMZ and YFL organized the manuscript. LLX, LJL and SHS reviewed the papers and revised the manuscript. All the authors have read and approved the final manuscript. All authors contributed toward data analysis, drafting, and revising the paper, and agree to be accountable for all aspects of the work.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Additional file 1: Figure S1.
Prognostic potential of S1PR1 in different cancers. (a–c) Kaplan-Meier survival curves comparing the high and low expression of S1PR1 in soft tissue cancer (a), blood cancer (b), and brain cancer (c) in the PrognoScan. (d–f) Survival curves of gastric cancer analyzed with mRNA-seq data of TCGA of Kaplan-Meier plotter databases. OS = Overall survival; RFS = Relapse-Free Survival; PPS = Post-progression survival.
Additional file 2: Table S1.
Correlations between S1PR1 and related genes and markers, as evaluated using GEPIA.
Rights and permissions
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.
About this article
Cite this article
Zhong, L., Xie, L., Yang, Z. et al. Prognostic value of S1PR1 and its correlation with immune infiltrates in breast and lung cancers. BMC Cancer 20, 766 (2020). https://doi.org/10.1186/s12885-020-07278-2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12885-020-07278-2