Jie Han1, Yihui Rong2, Xudong Gao2. 1. Department of Hepatology, Qilu Hospital, Shandong University, Shandong, China. 2. Infection Disease Center of Peking University International Hospital, Beijing, China.
Abstract
OBJECTIVE: To investigate SPARC (osteonectin), cwcv and kazal like domains proteoglycan 1 (SPOCK1) gene expression across The Cancer Genome Atlas (TCGA) cancers, both in cancer versus normal tissues and in different stages across the cancer types. METHODS: This integrated bioinformatics study used data from several bioinformatics databases (Cancer Cell Line Encyclopedia, Genotype-Tissue Expression, TCGA, Tumor Immune Estimation Resource [TIMER]) to define the expression pattern of the SPOCK1 gene. A survival analysis was undertaken across the cancers. The search tool for retrieval of interacting genes (STRING) database was used to identify proteins that interacted with SPOCK1. Gene Set Enrichment Analysis was conducted to determine pathway enrichment. The TIMER database was used to explore the correlation between SPOCK1 and immune cell infiltration. RESULTS: This multiomic analysis showed that the SPOCK1 gene was expressed differently between normal tissues and tumours in several cancers and that it was involved in cancer progression. The overexpression of the SPOCK1 gene was associated with poor clinical outcomes. Analysis of gene expression and tumour-infiltrating immune cells showed that SPOCK1 correlated with several immune cells across cancers. CONCLUSIONS: This research showed that SPOCK1 might serve as a new target for several cancer therapies in the future.
OBJECTIVE: To investigate SPARC (osteonectin), cwcv and kazal like domains proteoglycan 1 (SPOCK1) gene expression across The Cancer Genome Atlas (TCGA) cancers, both in cancer versus normal tissues and in different stages across the cancer types. METHODS: This integrated bioinformatics study used data from several bioinformatics databases (Cancer Cell Line Encyclopedia, Genotype-Tissue Expression, TCGA, Tumor Immune Estimation Resource [TIMER]) to define the expression pattern of the SPOCK1 gene. A survival analysis was undertaken across the cancers. The search tool for retrieval of interacting genes (STRING) database was used to identify proteins that interacted with SPOCK1. Gene Set Enrichment Analysis was conducted to determine pathway enrichment. The TIMER database was used to explore the correlation between SPOCK1 and immune cell infiltration. RESULTS: This multiomic analysis showed that the SPOCK1 gene was expressed differently between normal tissues and tumours in several cancers and that it was involved in cancer progression. The overexpression of the SPOCK1 gene was associated with poor clinical outcomes. Analysis of gene expression and tumour-infiltrating immune cells showed that SPOCK1 correlated with several immune cells across cancers. CONCLUSIONS: This research showed that SPOCK1 might serve as a new target for several cancer therapies in the future.
Cancer is among the most lethal diseases worldwide and it is an increasing threat to human health.
The number of cancer patients continues to expand due to the growth and aging of the global population.
Over the last few decades, many efforts have been invested in cancer prevention, early diagnosis and treatments to reduce the disease burden.
The process of carcinogenesis is extremely complex, with normal cells becoming pathological with disrupted apoptosis and dysregulation of metabolism.
The accumulation of gene alterations is key to oncogenesis and closely related to the prognosis of cancer patients.
Hence identification of genes that are involved can be utilized as new markers.
An understanding of the mechanism of the altered expression of these genes will enable them to be exploited as novel therapeutic targets and shed new light on cancer treatments.The SPARC (osteonectin), cwcv and kazal like domains proteoglycan 1 (SPOCK1) gene encodes a protein that is a member of the secreted protein acidic and rich in cysteine (SPARC) family.
Other members of the gene family include SPOCK2, SPOCK3, SPARC like 1 (SPARCL1), SPARC related modular calcium binding 1 (SMOC1) and SMOC2.
The proteins encoded by the members of the SPARC gene family are similar in structure, consisting of an N-terminus, follistatin-like domain and a C-terminus.
They have similar functions, for example they play important roles in the process of cancer cell adhesion and cell–matrix interactions, cell proliferation, migration and apoptosis.
Several studies have explored the function of SPOCK1 in light of its similar structure to SPARC and found that it participates in cancer cell invasion in oesophageal squamous cell carcinoma,
colorectal cancer[11,12] and gallbladder cancer.
Research has also shown that SPOCK1 plays an important role in prostate cancer recurrence.
The over-expression of the SPOCK1 gene promoted glioma cell proliferation, migration and invasion via the phosphatidylinositol-4,5-bisphosphate 3-kinase (PI3K)/protein kinase B (AKT) and Wnt/β-catenin signalling pathways; and SPOCK1 silencing reversed this process.
Further research showed that SPOCK1 was upregulated in osimertinib-resistant lung cancer cells and knockdown of SPOCK1 inhibited osimertinib-resistant cell growth and overcame resistance.
Recent research demonstrated that SPOCK1 could induce mesenchymal epithelial transition factor receptor-dependent epithelial–mesenchymal transition (EMT) signalling in lapatinib-resistant gastric cancer, which is responsible for gastric cancer drug resistance.
The over-expression of SPOCK1 also promoted proliferation and invasion and blocked apoptosis of HCC cells.
SPOCK1 can function on its own, as well as by working together with other elements like fibronectin.
Since the expression of SPOCK1 is low among normal tissue specimens,[10-13] it could be used as a potential biomarker for the early detection and precise prognosis in a variety of cancers. The underlying oncogenic mechanisms of the SPOCK1 gene during cancer initiation and progression needs further investigation and evaluation.All of the previous studies suggested that SPOCK1 plays a critical role in human tumorigenesis,[10-19] but details about SPOCK1-mediated cancer progression are lacking. This current bioinformatics analysis investigated SPOCK1 gene expression across The Cancer Genome Atlas (TCGA) cancers, both in cancer versus normal tissues and in different stages across cancers. It also explored important genes and proteins that have close interactions with SPOCK1.
Materials and methods
Expression of the SPOCK1 gene across diverse cancer types
Clinicopathological data, RNA-Seq data (HTSeq–FPKM) and immune subtypes data acquisition was achieved by downloading from TCGA
and The Genotype-Tissue Expression (GTEx).
The GTEx database was used as a supplementation to the noncancerous control group for cancer types that had less than 30 adjacent normal tissues in TCGA. An online tool, the Broad Institute Cancer Cell Line Encyclopedia (CCLE),
was used for visualization of the pan-cancer expression of SPOCK1 across diverse cancer types.The cancer types that were investigated were bladder urothelial carcinoma (BLCA), cholangiocarcinoma (CHOL), colon adenocarcinoma (COAD), head and neck squamous cell carcinoma (HNSC), kidney renal clear cell carcinoma (KIRC), lung adenocarcinoma (LUAD), pancreatic adenocarcinoma (PAAD), prostate adenocarcinoma (PRAD), sarcoma (SARC), stomach adenocarcinoma (STAD) and uveal melanoma (UVM).
Role of SPOCK1 gene over-expression in cancer prognosis
The expression levels of SPOCK1 across several types of cancers and their subtypes and its differential gene expression analysis at different pathological stages was analysed. The current study also explored the correlation between the expression of SPOCK1 and cancer prognosis in cancers. Overall survival (OS) and the disease-free survival (DFS) analysis was undertaken on the basis of SPOCK1 expression using Gene Expression Profiling Interactive Analysis.
Specifically, paired Student’s t-test was applied to compare the mRNA expression from the SPOCK1 gene in normal and tumour samples. Correlation between SPOCK1 expression and different individual cancer stages was analysed using analysis of variance. The OS and DFS of two groups (cut-off set as the median expression levels of SPOCK1) were compared using the Kaplan–Meier method. Survival maps were plotted with a significance level of P < 0.05 and false discovery rate (FDR) adjustment.
Protein–protein interaction network analysis
The search tool for retrieval of interacting genes (STRING; https://string-db.org/) database was utilized to generate protein–protein interaction (PPI) networks.
Functional interactions between SPOCK1 and other genes were analysed using Pearson's correlation coefficient. In the PPI network, nodes represent proteins and the edges represent interactions between the proteins. The median confidence score was set at 0.4. R software was used to perform Gene Set Enrichment Analysis (GSEA; https://www.gsea-msigdb.org/gsea/index.jsp) enrichment analysis where an FDR <0.05 and an enrichment score >0.65 were set as the cut-off criteria.
Correlation between SPOCK1 and tumour-infiltrating immune cells
In order to explore the correlation between SPOCK1 and tumour-infiltrating immune cells, the Tumor Immune Estimation Resource (TIMER) was utilized.
Spearman's rank correlation coefficient was used and a P-value < 0.05 was considered statistically significant.
Results
The CCLE database was used to visualize the expression of the SPOCK1 gene across diverse cancer types (Figure 1). Figure 2 illustrates the expression of the SPOCK1 gene (transcripts per million) in GTEx normal tissues, TCGA normal tissues and TCGA cancer tissues (Tables 1 and 2). In the dot plot in Figure 2, the dots represent the expression levels of the samples.
Figure 1.
The Broad Institute Cancer Cell Line Encyclopedia was used to visualize the pan-cancer expression of the SPARC (osteonectin), cwcv and kazal like domains proteoglycan 1 (SPOCK1) gene across diverse cancer types. Central bold horizontal lines represent the median values. Extremities of the box represent the 25th and 75th percentiles. Error bars represent minimum and maximum outliers. The colour version of this figure is available at: http://imr.sagepub.com
Figure 2.
SPARC (osteonectin), cwcv and kazal like domains proteoglycan 1 (SPOCK1) gene expression in Genotype-Tissue Expression normal tissues, The Cancer Genome Atlas (TCGA) normal tissues and TCGA cancer tissues. ACC, adrenocortical carcinoma; BLCA, bladder urothelial carcinoma; BRCA, breast invasive carcinoma; CESC, cervical squamous cell carcinoma and endocervical adenocarcinoma; CHOL, cholangiocarcinoma; COAD, colon adenocarcinoma; DLBC, lymphoid neoplasm diffuse large B-cell lymphoma; ESCA, oesophageal carcinoma; GBM, glioblastoma multiforme; HNSC, head and neck squamous cell carcinoma; KICH, kidney chromophobe; KIRC, kidney renal clear cell carcinoma; KIRP, kidney renal papillary cell carcinoma; LAML, acute myeloid leukaemia; LGG, brain lower grade glioma; LIHC, liver hepatocellular carcinoma; LUAD, lung adenocarcinoma; LUSC, lung squamous cell carcinoma; MESO, mesothelioma; OV, ovarian serous cystadenocarcinoma; PAAD, pancreatic adenocarcinoma; PCPG, pheochromocytoma and paraganglioma; PRAD, prostate adenocarcinoma; READ, rectum adenocarcinoma; SARC, sarcoma; SKCM, skin cutaneous melanoma; STAD, stomach adenocarcinoma; TGCT, testicular germ cell tumours; THCA, thyroid carcinoma; THYM, thymoma; UCEC, uterine corpus endometrial carcinoma; UCS, uterine carcinosarcoma; UVM, uveal melanoma. The colour version of this figure is available at: http://imr.sagepub.com
Table 1.
SPARC (osteonectin), cwcv and kazal like domains proteoglycan 1 (SPOCK1) gene expression in Genotype-Tissue Expression normal tissues and The Cancer Genome Atlas normal tissues across various cancers.
The Broad Institute Cancer Cell Line Encyclopedia was used to visualize the pan-cancer expression of the SPARC (osteonectin), cwcv and kazal like domains proteoglycan 1 (SPOCK1) gene across diverse cancer types. Central bold horizontal lines represent the median values. Extremities of the box represent the 25th and 75th percentiles. Error bars represent minimum and maximum outliers. The colour version of this figure is available at: http://imr.sagepub.comSPARC (osteonectin), cwcv and kazal like domains proteoglycan 1 (SPOCK1) gene expression in Genotype-Tissue Expression normal tissues, The Cancer Genome Atlas (TCGA) normal tissues and TCGA cancer tissues. ACC, adrenocortical carcinoma; BLCA, bladder urothelial carcinoma; BRCA, breast invasive carcinoma; CESC, cervical squamous cell carcinoma and endocervical adenocarcinoma; CHOL, cholangiocarcinoma; COAD, colon adenocarcinoma; DLBC, lymphoid neoplasm diffuse large B-cell lymphoma; ESCA, oesophageal carcinoma; GBM, glioblastoma multiforme; HNSC, head and neck squamous cell carcinoma; KICH, kidney chromophobe; KIRC, kidney renal clear cell carcinoma; KIRP, kidney renal papillary cell carcinoma; LAML, acute myeloid leukaemia; LGG, brain lower grade glioma; LIHC, liver hepatocellular carcinoma; LUAD, lung adenocarcinoma; LUSC, lung squamous cell carcinoma; MESO, mesothelioma; OV, ovarian serous cystadenocarcinoma; PAAD, pancreatic adenocarcinoma; PCPG, pheochromocytoma and paraganglioma; PRAD, prostate adenocarcinoma; READ, rectum adenocarcinoma; SARC, sarcoma; SKCM, skin cutaneous melanoma; STAD, stomach adenocarcinoma; TGCT, testicular germ cell tumours; THCA, thyroid carcinoma; THYM, thymoma; UCEC, uterine corpus endometrial carcinoma; UCS, uterine carcinosarcoma; UVM, uveal melanoma. The colour version of this figure is available at: http://imr.sagepub.comSPARC (osteonectin), cwcv and kazal like domains proteoglycan 1 (SPOCK1) gene expression in Genotype-Tissue Expression normal tissues and The Cancer Genome Atlas normal tissues across various cancers.BLCA, bladder urothelial carcinoma; BRCA, breast invasive carcinoma; CESC, cervical squamous cell carcinoma and endocervical adenocarcinoma; CHOL, cholangiocarcinoma; COAD, colon adenocarcinoma; ESCA, oesophageal carcinoma; HNSC, head and neck squamous cell carcinoma; KIRC, kidney renal clear cell carcinoma; LIHC, liver hepatocellular carcinoma; LUAD, lung adenocarcinoma; LUSC, lung squamous cell carcinoma; PAAD, pancreatic adenocarcinoma; PRAD, prostate adenocarcinoma; STAD, stomach adenocarcinoma; READ, rectum adenocarcinoma.SPARC (osteonectin), cwcv and kazal like domains proteoglycan 1 (SPOCK1) gene expression in The Cancer Genome Atlas cancer tissues.BLCA, bladder urothelial carcinoma; BRCA, breast invasive carcinoma; CESC, cervical squamous cell carcinoma and endocervical adenocarcinoma; CHOL, cholangiocarcinoma; COAD, colon adenocarcinoma; ESCA, oesophageal carcinoma; HNSC, head and neck squamous cell carcinoma; KIRC, kidney renal clear cell carcinoma; LIHC, liver hepatocellular carcinoma; LUAD, lung adenocarcinoma; LUSC, lung squamous cell carcinoma; PAAD, pancreatic adenocarcinoma; PRAD, prostate adenocarcinoma; STAD, stomach adenocarcinoma; READ, rectum adenocarcinoma.The SPOCK1 gene has ten splice variants: SPOCKl-201, SPOCK1-001, SPOCK1-002, SPOCK1-003, SPOCK1-004, SPOCK1-005, SPOCK1-006, SPOCK1-007, SPOCK1-008, SPOCK1-009. The predominant SPOCK1 isoforms expressed in all contexts were SPOCK1-201, SPOCK1-001 and SPOCK1-003 (Figure 3).
Figure 3.
Different isoforms of the SPARC (osteonectin), cwcv and kazal like domains proteoglycan 1 (SPOCK1) gene.
Different isoforms of the SPARC (osteonectin), cwcv and kazal like domains proteoglycan 1 (SPOCK1) gene.The levels of SPOCK1 gene expression between TCGA tumour tissues and normal tissues were significantly different in cancers such as CHOL, PAAD, PRAD and STAD (P < 0.05) (Figure 4 and Table 3). The expression of the SPOCK1 gene was significantly different among pathological stages in cancers such as BLCA, KIRC, LUAD and STAD (P < 0.05) (Figure 5 and Table 4).
Figure 4.
Comparison of SPARC (osteonectin), cwcv and kazal like domains proteoglycan 1 (SPOCK1) gene expression between tumour tissues (red boxes) and normal tissues (grey boxes). CHOL, cholangiocarcinoma; PAAD, pancreatic adenocarcinoma; PRAD, prostate adenocarcinoma; STAD, stomach adenocarcinoma. *P < 0.05, paired Student’s t-test. Central black horizontal lines represent the median values. Extremities of the box represent the 25th and 75th percentiles. Error bars represent minimum and maximum outliers. The colour version of this figure is available at: http://imr.sagepub.com
Table 3.
The levels of SPARC (osteonectin), cwcv and kazal like domains proteoglycan 1 (SPOCK1) gene expression in cancer tissues compared with normal tissues.
Violin plots showing the comparison of SPARC (osteonectin), cwcv and kazal like domains proteoglycan 1 (SPOCK1) gene expression between tumour tissues at different pathological stages. (a) BLCA, bladder urothelial carcinoma; (b) KIRC, kidney renal clear cell carcinoma; (c) LUAD, lung adenocarcinoma; (d) STAD, stomach adenocarcinoma. Cancer stages were compared using analysis of variance. The colour version of this figure is available at: http://imr.sagepub.com
Table 4.
The levels of SPARC (osteonectin), cwcv and kazal like domains proteoglycan 1 (SPOCK1) gene expression in tumour tissues at different pathological stages.
Comparison of SPARC (osteonectin), cwcv and kazal like domains proteoglycan 1 (SPOCK1) gene expression between tumour tissues (red boxes) and normal tissues (grey boxes). CHOL, cholangiocarcinoma; PAAD, pancreatic adenocarcinoma; PRAD, prostate adenocarcinoma; STAD, stomach adenocarcinoma. *P < 0.05, paired Student’s t-test. Central black horizontal lines represent the median values. Extremities of the box represent the 25th and 75th percentiles. Error bars represent minimum and maximum outliers. The colour version of this figure is available at: http://imr.sagepub.comThe levels of SPARC (osteonectin), cwcv and kazal like domains proteoglycan 1 (SPOCK1) gene expression in cancer tissues compared with normal tissues.aStudent’s t-test.TPM, transcripts per million; CHOL, cholangiocarcinoma; COAD, colon adenocarcinoma; HNSC, head and neck squamous cell carcinoma; KIRC, kidney renal clear cell carcinoma; LIHC, liver hepatocellular carcinoma; LUSC, lung squamous cell carcinoma; PRAD, prostate adenocarcinoma.Violin plots showing the comparison of SPARC (osteonectin), cwcv and kazal like domains proteoglycan 1 (SPOCK1) gene expression between tumour tissues at different pathological stages. (a) BLCA, bladder urothelial carcinoma; (b) KIRC, kidney renal clear cell carcinoma; (c) LUAD, lung adenocarcinoma; (d) STAD, stomach adenocarcinoma. Cancer stages were compared using analysis of variance. The colour version of this figure is available at: http://imr.sagepub.comThe levels of SPARC (osteonectin), cwcv and kazal like domains proteoglycan 1 (SPOCK1) gene expression in tumour tissues at different pathological stages.aAnalysis of variance.BLCA, bladder urothelial carcinoma; CESC, cervical squamous cell carcinoma and endocervical adenocarcinoma; COAD, colon adenocarcinoma; ESCA, oesophageal carcinoma; HNSC, head and neck squamous cell carcinoma; KIRC, kidney renal clear cell carcinoma; LUAD, lung adenocarcinoma; READ, rectum adenocarcinoma; UCEC, uterine corpus endometrial carcinoma; UCS, uterine carcinosarcoma; UVM, uveal melanoma.The OS time between cancers with higher and lower SPOCK1 gene expression were compared in TCGA cancer types (Figure 6). Data demonstrated a shorter OS with a worse prognosis in patients with cancers with higher SPOCK1 gene expression levels compared with cancers with lower SPOCK1 gene expression levels for the following cancers: COAD, HNSC, KIRC, LUAD and UVM (P < 0.05). For DFS time in TCGA tumour types (Figure 7), data demonstrated that higher SPOCK1 gene expression levels resulted in worse DFS prognosis compared with lower levels of SPOCK1 gene expression in the following tumours: BLCA, COAD, KIRC, SARC and UVM (P < 0.05).
Figure 6.
Kaplan–Meier analysis of the overall survival (OS) time between higher SPARC (osteonectin), cwcv and kazal like domains proteoglycan 1 (SPOCK1) gene expression and lower SPOCK1 gene expression in The Cancer Genome Atlas tumour types. The other lines represent the 25% and 75% cut-offs. (a) COAD, colon adenocarcinoma; (b) HNSC, head and neck squamous cell carcinoma; (c) KIRC, kidney renal clear cell carcinoma; (d) LGG, brain lower grade glioma; (e) LUAD, lung adenocarcinoma; (f) UVM, uveal melanoma. The colour version of this figure is available at: http://imr.sagepub.com
Figure 7.
Kaplan–Meier analysis of the disease-free survival (DFS) time between higher SPARC (osteonectin), cwcv and kazal like domains proteoglycan 1 (SPOCK1) gene expression and lower SPOCK1 gene expression in The Cancer Genome Atlas tumour types. (a) BLCA, bladder urothelial carcinoma; (b) COAD, colon adenocarcinoma; (c) KIRC, kidney renal clear cell carcinoma; (d) LGG, brain lower grade glioma; (e) SARC, sarcoma; UVM, (f) uveal melanoma. The colour version of this figure is available at: http://imr.sagepub.com
Kaplan–Meier analysis of the overall survival (OS) time between higher SPARC (osteonectin), cwcv and kazal like domains proteoglycan 1 (SPOCK1) gene expression and lower SPOCK1 gene expression in The Cancer Genome Atlas tumour types. The other lines represent the 25% and 75% cut-offs. (a) COAD, colon adenocarcinoma; (b) HNSC, head and neck squamous cell carcinoma; (c) KIRC, kidney renal clear cell carcinoma; (d) LGG, brain lower grade glioma; (e) LUAD, lung adenocarcinoma; (f) UVM, uveal melanoma. The colour version of this figure is available at: http://imr.sagepub.comKaplan–Meier analysis of the disease-free survival (DFS) time between higher SPARC (osteonectin), cwcv and kazal like domains proteoglycan 1 (SPOCK1) gene expression and lower SPOCK1 gene expression in The Cancer Genome Atlas tumour types. (a) BLCA, bladder urothelial carcinoma; (b) COAD, colon adenocarcinoma; (c) KIRC, kidney renal clear cell carcinoma; (d) LGG, brain lower grade glioma; (e) SARC, sarcoma; UVM, (f) uveal melanoma. The colour version of this figure is available at: http://imr.sagepub.comThe gene interaction network showed that the SPOCK1 gene was co-expressed with the TNFAIP6, MLLT10, MEGF11, SGCB, DCTN1, SLC24A2, CNTN1 and LSAMP genes (Figure 8). The SPOCK1 gene had physical interactions with the GTF2E2 and TNF genes. The SPOCK1 gene had shared protein domains with the TACSTD2, SPOCK3 and SPOCK2 genes; and co-localized with the SLC25A33, SFXN3, MEGF11, SGCB, GRIA4, SLC24A2 and LSAMP genes. Data from the STRING database analysis revealed that all of the proteins from these genes (CHD1L, CGREF1, SERTM1, MMP14, MMP16, SPARC, SMOC1, TG, STS, and SPARCL1) were in the same protein network with SPOCK1 (Figure 9).
Figure 8.
Genes in correlations with the SPARC (osteonectin), cwcv and kazal like domains proteoglycan 1 (SPOCK1) gene.
Figure 9.
The protein–protein interactions with the SPARC (osteonectin), cwcv and kazal like domains proteoglycan 1 (SPOCK1) gene.
Genes in correlations with the SPARC (osteonectin), cwcv and kazal like domains proteoglycan 1 (SPOCK1) gene.The protein–protein interactions with the SPARC (osteonectin), cwcv and kazal like domains proteoglycan 1 (SPOCK1) gene.The GSEA enrichment analysis revealed that the enrichment of the SPOCK1 gene was found to mediate upstream glycolysis, G2M checkpoint, mammalian target of rapamycin complex 1 signalling, mitotic spindle activities, EMT and myc targets; and to mediate downstream E2 factor targets, tumour necrosis factor-α signalling via nuclear factor kappa-light-chain-enhancer of activated B cells, oestrogen response, hypoxia and the interleukin-2-signal transducer and activator of transcription 5 signalling process (Figure 10).
Figure 10.
The Gene Set Enrichment Analysis enrichment analysis of the SPARC (osteonectin), cwcv and kazal like domains proteoglycan 1 (SPOCK1) gene. The colour version of this figure is available at: http://imr.sagepub.com
The Gene Set Enrichment Analysis enrichment analysis of the SPARC (osteonectin), cwcv and kazal like domains proteoglycan 1 (SPOCK1) gene. The colour version of this figure is available at: http://imr.sagepub.comFigure 11 presents the correlation between the expression of the SPOCK1 gene and tumour-infiltrating immune cells in BLCA, KIRC, LUAD and STAD based on data from the TIMER database. The results showed that in BLCA, the expression of the SPOCK1 gene had significant inverse correlations with tumour purity (partial correlation = –0.352, P < 0.0001) and the infiltration of B cells (partial correlation = –0.152, P < 0.0001; and significant positive correlations with the infiltration of CD4+ T cells (partial correlation = 0.17, P = 0.001096) and macrophages (partial correlation = 0.352, P < 0.0001).
Figure 11.
Correlation analysis between tumour-infiltrating immune cell signatures and the SPARC (osteonectin), cwcv and kazal like domains proteoglycan 1 (SPOCK1) gene. Spearman's rank correlation coefficient was used and a P-value <0.05 was considered statistically significant. BLCA, bladder urothelial carcinoma; KIRC, kidney renal clear cell carcinoma; LUAD, lung adenocarcinoma; STAD, stomach adenocarcinoma. The colour version of this figure is available at: http://imr.sagepub.com
Correlation analysis between tumour-infiltrating immune cell signatures and the SPARC (osteonectin), cwcv and kazal like domains proteoglycan 1 (SPOCK1) gene. Spearman's rank correlation coefficient was used and a P-value <0.05 was considered statistically significant. BLCA, bladder urothelial carcinoma; KIRC, kidney renal clear cell carcinoma; LUAD, lung adenocarcinoma; STAD, stomach adenocarcinoma. The colour version of this figure is available at: http://imr.sagepub.comIn KIRC, the expression of the SPOCK1 gene had a significant inverse correlation with tumour purity (partial correlation = −0.171, P = 0.000224); and significant positive correlations with the infiltration of B cells (partial correlation =0.175, P = 0.000164), neutrophils (partial correlation = 0.134, P = 0.003934) and dendritic cells (partial correlation = 0.123, P = 0.00889).In LUAD, the expression of the SPOCK1 gene had significant inverse correlations with tumour purity (partial correlation = –0.169, P = 0.000166) and the infiltration of B cells (partial correlation =–0.151, P = 0.000881); and significant positive correlations with the infiltration of neutrophils (partial correlation = 0.179, P < 0.0001) and dendritic cells (partial correlation = 0.114, P = 0.011486).In STAD, the expression of the SPOCK1 gene had significant inverse correlations with tumour purity (partial correlation = −0.129, P = 0.011966) and the infiltration of B cells (partial correlation = –0.111, P = 0.033804); and significant positive correlations with the infiltration of CD8+ T cells (partial correlation = 0.116, P = 0.026166), CD4+ T cells (partial correlation = 0.265, P < 0.0001), macrophages (partial correlation = 0.583, P < 0.0001), neutrophils (partial correlation = 0.212, P < 0.0001) and dendritic cells (partial correlation = 0.368, P < 0.0001).
Discussions
The SPOCK1 gene is overexpressed in prostate cancer, lung cancer, ovarian cancer, gastric cancer, colorectal cancer and breast cancer and plays a key role in cancer progression,[9,17,27-32] which suggests that SPOCK1 could be a novel gene of interest in the search of new therapeutic targets. The SPOCK1 gene participates in regulating cancer cell proliferation, cell cycle regulation, apoptosis, adhesion, cell–matrix interactions, metastasis and drug resistance in these cancers.[9,17,21-32] However, the underlying mechanism of action of the SPOCK1 gene is still not completely understood. To address this lack of data, this current study implemented a comprehensive bioinformatics approach to discover the function of the SPOCK1 gene in a pan-cancer setting. Analysis of SPOCK1 gene expression was performed across all TCGA cancer types compared with TCGA normal tissues in addition to GTEx normal tissues. The results of this current study demonstrated that the SPOCK1 gene was upregulated in several types of cancer, showing that the SPOCK1 gene is functionally active in these tumours. The current analysis showed that in CHOL, PAAD, PRAD and STAD, the expression of the SPOCK1 gene was notably higher in cancer tissues than in normal tissues. In BLCA, KIRC, LUAD and STAD, the differences between the pathological stages were statistically significant. The stages of the disease and the degree of deterioration are related, indicating that the SPOCK1 gene might govern the initiation and progression processes of these cancers. In the survival analysis, cancers (COAD, HNSC, KIRC, LUAD and UVM), which have highly elevated the SPOCK1 gene expression, presented with worse overall survival than those that showed lower SPOCK1 gene expression. As for the DFS time in the TCGA tumour types, this current research showed that higher expression levels of the SPOCK1 gene resulted in worse DFS prognosis in BLCA, COAD, KIRC, SARC and UVM, emphasizing the importance of the SPOCK1 gene for prognosis. Interestingly, in lower grade glioma of the brain, lower expression of the SPOCK1 gene leads to worse OS and DFS, which might have correlations with the relatively higher expression of the SPOCK1 gene in the brain;
and the underlying reason for this unusual phenomena deserves further investigation.The CHD1L, CGREF1, SERTM1, MMP14, MMP16, SPARC, SMOC1, TG, STS, and SPARCL1 genes correlated with the SPOCK1 gene in the protein–protein network. Among the SPOCK1 gene protein partners, matrix metalloproteinase 14 (MMP14) and MMP16 are involved in the breakdown of the extracellular matrix.
The PPI network showed interactions between the SPOCK1 gene and general transcription factor IIE subunit 2 (GTF2E2), tumour necrosis factor (TNF), lysyl oxidase like 2 (LOXL2) and multiple EGF like domrix in cancer metastasis (MEGF11); and with CHD1L, which is involved in DNA replication, repair and transcription processes.
The GSEA enrichment analysis showed that the SPOCK1 gene might have positive interactions with the EMT process. Previous research also indicated that the SPOCK1 gene mediates EMT signalling to regulate cancer cell progression and drug resistance.
Checkpoint immunotherapy has been studied comprehensively in recent years.
However, the underlying mechanisms of how tumour cells and immune cells interact remains elusive. This current study showed that the correlation between the SPOCK1 gene and tumour-infiltrating immune cells might be of novel therapeutic value to the immunotherapy of cancers.In conclusion, the SPOCK1 gene has been demonstrated to be an oncogene that is involved in major oncogenic processes including cell-cycle control, DNA repair, apoptosis and metastasis. It may affect the migration and invasion of cancer cells in several cancers through the EMT process. It has been confirmed that the overexpression of the SPOCK1 gene can promote the occurrence and development of tumours, suggesting that it may become a new anti-tumour therapeutic target. This current systematic pan-cancer analysis of SPOCK1 gene expression in several cancer databases has provided evidence of the relationship between the altered expression of the SPOCK1 gene and clinical outcomes. This current study uncovered the importance of SPOCK1 gene expression and possible SPOCK1-related proteins and pathways in cancer progression. Therefore, this current analysis may provide valuable insights into the SPOCK1 gene as a potential biomarker and therapeutic target for various human cancers.
Authors: John N Weinstein; Eric A Collisson; Gordon B Mills; Kenna R Mills Shaw; Brad A Ozenberger; Kyle Ellrott; Ilya Shmulevich; Chris Sander; Joshua M Stuart Journal: Nat Genet Date: 2013-10 Impact factor: 38.330
Authors: Taiwen Li; Jingyu Fan; Binbin Wang; Nicole Traugh; Qianming Chen; Jun S Liu; Bo Li; X Shirley Liu Journal: Cancer Res Date: 2017-11-01 Impact factor: 12.701
Authors: Mahmoud Ghandi; Franklin W Huang; Judit Jané-Valbuena; Gregory V Kryukov; Christopher C Lo; E Robert McDonald; Jordi Barretina; Ellen T Gelfand; Craig M Bielski; Haoxin Li; Kevin Hu; Alexander Y Andreev-Drakhlin; Jaegil Kim; Julian M Hess; Brian J Haas; François Aguet; Barbara A Weir; Michael V Rothberg; Brenton R Paolella; Michael S Lawrence; Rehan Akbani; Yiling Lu; Hong L Tiv; Prafulla C Gokhale; Antoine de Weck; Ali Amin Mansour; Coyin Oh; Juliann Shih; Kevin Hadi; Yanay Rosen; Jonathan Bistline; Kavitha Venkatesan; Anupama Reddy; Dmitriy Sonkin; Manway Liu; Joseph Lehar; Joshua M Korn; Dale A Porter; Michael D Jones; Javad Golji; Giordano Caponigro; Jordan E Taylor; Caitlin M Dunning; Amanda L Creech; Allison C Warren; James M McFarland; Mahdi Zamanighomi; Audrey Kauffmann; Nicolas Stransky; Marcin Imielinski; Yosef E Maruvka; Andrew D Cherniack; Aviad Tsherniak; Francisca Vazquez; Jacob D Jaffe; Andrew A Lane; David M Weinstock; Cory M Johannessen; Michael P Morrissey; Frank Stegmeier; Robert Schlegel; William C Hahn; Gad Getz; Gordon B Mills; Jesse S Boehm; Todd R Golub; Levi A Garraway; William R Sellers Journal: Nature Date: 2019-05-08 Impact factor: 49.962
Authors: Lóránd Váncza; Katalin Karászi; Bálint Péterfia; Lilla Turiák; Katalin Dezső; Anna Sebestyén; Andrea Reszegi; Gábor Petővári; András Kiss; Zsuzsanna Schaff; Kornélia Baghy; Ilona Kovalszky Journal: Front Oncol Date: 2022-02-03 Impact factor: 6.244
Figure 1.
Figure 3.
Figure 8.
Figure 9.