ORIGINAL RESEARCH article
Sec. Computational Genomics
Volume 11 - 2020 | https://doi.org/10.3389/fgene.2020.00547
Analysis of Expression and Its Clinical Significance of the Secreted Phosphoprotein 1 in Lung Adenocarcinoma
- 1Department of Thoracic Surgery, Zhongnan Hospital of Wuhan University, Wuhan, China
- 2Department of Biological Repositories, Zhongnan Hospital of Wuhan University, Wuhan, China
- 3Hubei Key Laboratory of Tumor Biological Behaviors, Hubei Cancer Clinical Study Center, Zhongnan Hospital of Wuhan University, Wuhan, China
- 4Department of Urology, Zhongnan Hospital of Wuhan University, Wuhan, China
- 5Department of Oncology, The First People’s Hospital of Tianmen, Tianmen, China
- 6Human Genetics Resource Preservation Center of Hubei Province, Wuhan, China
Objective: To explore the expression of secreted phosphoprotein 1 (SPP1) in lung adenocarcinoma (LUAD), and evaluate its relationship with clinicopathological characteristics and prognosis of LUAD, and analyze the advantages of SPP1 as a potential prognostic marker in LUAD.
Methods: The expression of SPP1 in normal lung tissue and LUAD was analyzed from the Cancer Cell Line Encyclopedia (CCLE), Gene Expression Omnibus (GEO), and Human Protein Atlas (HPA) databases. GSE68465 was used to explore the relationship between the SPP1 expression and clinicopathological characteristics and the prognosis of LUAD patients. The relationship between SPP1 and immune infiltration in LUAD was analyzed by the Tumor Immune Estimation Resource (TIMER) database. Gene enrichment analysis was performed in GSEA. The Cancer Genome Atlas (TCGA)-LUAD data was used to verify the results.
Results: In the cell line level, non-small cell lung cancer ranked ninth among cancer cell lines based on SPP1 expression. In the messenger RNA (mRNA) and protein levels, SPP1 expression was higher in LUAD tissues than that in normal control. SPP1 expression was related to gender, N stage, histological grade, and progression or relapse. In men, SPP1 expression were higher compared to that in women. The higher the N stage, the higher the SPP1 expression level. As LUAD progresses or relapses, SPP1 expression could increase. In the pathological grade, the SPP1 expression was higher in LUAD samples with moderate differentiation. In addition, the overall 5-year survival rates of the SPP1 high and low expression groups were 50.574 and 59.181% [P = 0.008; hazard ratio (HR) = 0.7057; 95% CI, 0.5467–0.9109], indicating that SPP1 had an impact on overall survival for LUAD patients. The relationship between SPP1 expression and CD4+ T cell, macrophage, neutrophil, and dendritic cell infiltration was weak in LUAD. SPP1 could be considered as an independent prognostic marker in LUAD (P = 0.003; HR = 1.150; 95% CI, 1.048–1.261) by multivariate Cox regression analysis. The results of GSEA indicated that samples with high SPP1 expression were enriched in protein secretion, mTORC1 signaling, angiogenesis, and glycolysis pathway. The analysis results obtained by TCGA-LUAD data were basically consistent with the results obtained by GSE68465.
Conclusions: SPP1 can not only affect the occurrence and development of LUAD but also may be an independent prognostic marker of LUAD. SPP1 is expected to be a new target for molecular targeted therapy.
Lung cancer is one of the malignant tumors that seriously endanger human health worldwide. The morbidity of lung cancer accounts for 11.6%, ranking first of all cancers, and 18.4% of the total number of cancer deaths are lung cancer (Bray et al., 2018). Non-small cell lung cancer (NSCLC) accounts for 80% of total lung cancer and is the most common histological type of lung cancer (Johnson et al., 2012). Lung adenocarcinoma (LUAD) accounts for ∼40% of all lung cancers, which is the most aggressive pathological type of NSCLC, and is an important cause of death from respiratory cancer (Noguchi et al., 1995). Compared with different types of lung cancers, LUAD has slower growth and no specific clinical manifestations in its early stage. It is easy to misdiagnosis because of the similar symptoms of common respiratory diseases (Lam et al., 2013; Berardi et al., 2016). Approximately 75% of patients are diagnosed with metastatic or advanced cancer and often lose the opportunity for radical surgery. With the improvements of modern clinical diagnostic methods, the survival rate is still unsatisfactory, only ∼15%, although the overall survival rate and the quality of life have been improved (Siegel et al., 2015; Franceschini et al., 2017). In recent years, in addition to improving the treatment of lung cancer patients, searching for effective tumor markers to guide clinical diagnosis has become a hot spot in LUAD research.
Secreted phosphoprotein 1 (SPP1), also known as osteopontin, is a secretory phosphorylated glycoprotein rich in arginine–glycine–aspartic acid sequence (RGD domain) (Wei et al., 2017). It is a widely expressed mucoprotein, which can be detected in normal tissues, body fluids, and cells (Liaw et al., 1998). SPP1 also participates in the regulation of physiological processes such as development, differentiation, inflammation, and the wound healing (Fedarko et al., 2001). In addition, more and more evidence showed that SPP1 was highly expressed in various kinds of tumors and plays a considerable role in the occurrence and metastasis of various tumors (Kijewska et al., 2017; Li Y. et al., 2017; Zeng et al., 2018; Zhang et al., 2018). The aim of this research was to explore the relevance between SPP1 and clinicopathological features of patients with LUAD and to determine whether it can be used as an independent prognostic marker of LUAD.
Materials and Methods
Cancer Cell Line Encyclopedia Analysis
Cancer Cell Line Encyclopedia (CCLE1) is a database of cancer cell lines maintained by the Broad Institute of MIT & Harvard, which is an open-access database including the RNA-Seq data for the current 1,457 cancer cell lines (Barretina et al., 2012). We used CCLE database to analyze the messenger RNA (mRNA) expression levels of the SPP1 in various kinds of cancer cell lines.
SPP1 mRNA Expressions in Normal Lung Tissues and LUAD Tissues
To compare the mRNA expression levels of SPP1 in LUAD with normal lung tissues, Gene Expression Omnibus (GEO2) database was employed for analyses. The mRNA expression levels were shown as the intensity of log2 median centered, and the SPP1 mRNA expression levels in LUAD tissues were compared with that in control normal lung tissues by Mann–Whitney U test. There were seven datasets that were used to compare the SPP1 expression in normal lung tissue and LUAD, including GSE2514 (20 LUAD samples and 19 normal lung tissue samples) (Stearman et al., 2005), GSE7670 (20 LUAD samples and 30 normal lung tissue samples) (Su et al., 2007; Chen et al., 2009), GSE10072 (58 LUAD samples and 49 normal lung tissue samples) (Landi et al., 2008), GSE19188 (45 LUAD samples and 65 normal lung tissue samples) (Hou et al., 2010), GSE31210 (226 LUAD samples and 20 normal lung tissue samples) (Okayama et al., 2012; Yamauchi et al., 2012), GSE32863 (58 LUAD samples and 58 normal lung tissue samples) (Selamat et al., 2012), and GSE83227 (132 LUAD samples and 17 normal lung tissue samples) (Bhattacharjee et al., 2001) from the GEO database. Each dataset was downloaded in Series Matrix File (s) format. The preprocess of each dataset adopted the original author’s preprocessing results.
The Human Protein Atlas Analysis
The protein levels of SPP1 expression in LUAD samples and normal lung samples from the Human Protein Atlas (HPA)3 database were analyzed by the immunohistochemistry (IHC) staining data (Uhlen et al., 2005, 2010, 2015). These data were then used to analyze the correlation between the mRNA levels of SPP1 expression and prognosis. The expression levels could be divided into four categories: high, medium, low, and not detected by the scoring system, which consists of the proportion of stained cells (>75, 25–75, or <25%) and the intensity of staining (strong, moderate, weak, or negative).
Data Mining From the GEO Database
GSE68465, including 19 normal lung tissues and 443 LUAD samples, was downloaded from the GEO database. GSE68465 was performed on GPL96 (Shedden et al., 2008). A large amount of clinicopathological information, such as gender, age, and tumor–node–metastasis (TNM) stage, was collected from the database as well. The associated probe for SPP1 was 209875_s_at. We then analyzed the survival information of samples with GSE68465 using the Kaplan–Meier method. GSE68465 were normalized by the MAS5 algorithm and log2 transformation were conducted.
Tumor Immune Estimation Resource (TIMER)4 systematically analyzed the comprehensive resource of immune infiltrates in multiple tumor types, including six types of immune infiltrates (CD4+ T cells, B cells, CD8+ T cells, macrophages, neutrophils, and dendritic cells) (Li et al., 2016; Li T. et al., 2017). We selected the parameters: “Gene,” gene symbol: “SPP1,” cancer types: “LUAD (lung adenocarcinoma),” immune infiltrates: “Default” to explore the relationship between the SPP1 expression and the proportion of immune infiltration. We selected the parameters “survival,” cancer types: “LUAD (Lung Adenocarcinoma),” clinical: (optional): “all,” Immune infiltrates: “Default,” gene symbol: “SPP1” to conduct survival analysis and build a multivariable Cox proportional hazards model. This model was fitted by the function coxph() using R package “survival.” Baseline variables considered to be univariate with outcome or clinically relevant to the results were included into multivariable Cox proportional hazards model (Stone et al., 2011).
Gene Set Enrichment Analysis
Gene Set Enrichment Analysis (GSEA) was used for gene enrichment analysis (Mootha et al., 2003; Subramanian et al., 2005). GSE68465 was included in GSEA, and lung cancer tissue samples were divided into low- and high-expression groups according to the median expression level of SPP1. The hallmark gene set in MSigDB database was used as reference gene set. The default weighted enrichment statistical method was applied to enrichment analysis, and the number of the gene set permutation for each analysis was set as 1,000 times. Pathways with false discovery rate (FDR) q < 0.25 and nominal P < 0.05 were considered as significantly enriched pathways.
Validation in the Cancer Genome Atlas Dataset
FPKM-UQ normalized mRNA expression of The Cancer Genome Atlas (TCGA)-LUAD and clinical data were downloaded from the GDC Data Portal5, and then, we screened the samples based on the complete survival information and gene expression data. The TCGA-LUAD gene expression data were processed by log2 transformation to compare the expression levels of SPP1 in normal lung tissue samples (N = 59) and LUAD samples (N = 504). According to the expression level of SPP1, the samples above the median expression level were set as high-expression group (N = 252), and the samples below the median expression level were set as low-expression group (N = 252) (Tang et al., 2018). Then, we explored the relationship between the SPP1 expression level and the prognosis of LUAD.
Mann–Whitney U test was used to compare the SPP1 expression in LUAD tissues and normal lung tissues from the GEO database. The relevance between clinicopathological characteristics and the SPP1 mRNA expression level was analyzed by χ2 test with the SPSS software (version 20). For the LUAD samples in GSE68465, according to the expression level of SPP1, the samples above the median expression level were set as high-expression group (N = 222), and the samples below the median expression level were set as low-expression group (N = 221) (Tian et al., 2018; Guo et al., 2020). Using the GraphPad Prism software (version 8) and Kaplan–Meier method, the relationship between the overall survival rate and the SPP1 expression level was analyzed. Univariate Cox regression model was conducted using SPSS software (version 20). Statistical significance was set as P < 0.05.
SPP1 mRNA Expressions in Normal Lung Tissues and Tumor
NSCLC ranked ninth among cancer cell lines based on the SPP1 expression according to the data from CCLE database (Figure 1A). In seven studies, the SPP1 mRNA expression levels in LUAD tissues were significantly higher than that in normal lung tissues (P < 0.05). Each study supported the above conclusion (Figures 1B–H).
Figure 1. (A) The messenger RNA (mRNA) expression levels of SPP1 in non-small cell lung cancer (NSCLC) and other cancer cell lines from the Cancer Cell Line Encyclopedia (CCLE) database. NSCLC ranked ninth among cancer cell lines based on SPP1 expression. (B–H) Comparison of SPP1 mRNA levels in normal lung tissues and lung adenocarcinoma (LUAD) tissues across seven analyses of LUAD, and each study showed that the SPP1 mRNA expression levels in LUAD tissues were significantly higher than that in normal lung tissues.
SPP1 Expression at the Protein Level
Based on the protein expression data from the HPA, the protein expression of SPP1 in LUAD tissues and control normal lung tissues were compared by antibody CAB002212. The SPP1 protein expression in normal lung tissues was “not detected” in two cases. However, the protein expression of SPP1 in LUAD was “medium” in two cases, “low” in three cases, and “not detected” in one case (Figures 2A,B). It was further confirmed that SPP1 protein was highly expressed in LUAD than in normal lung tissues. The survival curves of 262 samples with high SPP1 expression and 238 samples with low SPP1 expression were obtained by analyzing the relevance between the SPP1 mRNA expression and prognosis in the HPA database (Figure 2C). It was found that compared with the SPP1 low-expression group, the overall survival rate of SPP1 high-expression group was significantly lower (log-rank P = 0.011).
Figure 2. (A) The protein expression of secreted phosphoprotein 1 (SPP1) was detected in normal lung tissues from the Human Protein Atlas (HPA) database. (B) The protein expression of SPP1 was medium in lung adenocarcinoma (LUAD) from the HPA database. (C) The survival curves of 262 samples with high SPP1 mRNA expression and 238 samples with low SPP1 mRNA expression in HPA database.
Association Between Clinicopathological Characteristics and SPP1 mRNA Expression
The SPP1 mRNA expression level in LUAD tissues was significantly higher than that in normal lung tissues in GSE68465 (P = 0.039; Figure 3A). Based on the clinical and pathological information from GSE68465, the SPP1 expression in LUAD was not related to age, T stage, and smoking history (P > 0.05) but related to the gender, N stage, histological grade, and progression and relapse (P < 0.05; Table 1). In men, the SPP1 expression levels were higher compared to that in women. In the progression of lymph node metastasis of LUAD, the higher the N stage, the higher the SPP1 expression level. As LUAD progresses or relapses, the SPP1 expression could increase. In the pathological grade, the expression level of SPP1 was higher in lung adenocarcinoma samples whit moderately differentiated. Based on the clinical follow-up data of GSE68465, the overall 5-year survival rates of the SPP1 high-expression group and the SPP1 low-expression group were 50.574 and 59.181%, respectively [hazard ratio (HR) = 0.7057; 95% CI, 0.5467–0.9109; P = 0.008; Figure 3B]. By univariate Cox regression model, we concluded that age, gender, T stage, N stage, and SPP1 had an impact on overall survival rate of LUAD patients (P < 0.05; Supplementary Table S1).
Figure 3. (A) The messenger RNA (mRNA) expression of secreted phosphoprotein 1 (SPP1) in lung adenocarcinoma (LUAD) and normal lung tissues based on GSE68465. (B) The overall survival rates of the SPP1 high-expression group and the SPP1 low-expression group in LUAD based on GSE68465.
Table 1. Associations between SPP1 expression and clinicopathological factors of patients with LUAD (based on GSE68465).
The Relationship Between SPP1 and Immune Infiltrates
Based on the TIMER database, the relationship between SPP1 expression and CD4+ T cell, neutrophil, macrophage, and dendritic cell infiltration was weak in LUAD (P < 0.05; Figure 4A). According to the Kaplan–Meier plots for immune infiltration and SPP1 expression, the survival rate of high level of B cell and dendritic cell infiltration group was significantly higher than that of low level of infiltration group (P < 0.05). At the same time, the group with low SPP1 expression had significantly higher survival rates than that group with high SPP1 expression (P < 0.05; Figure 4B).
Figure 4. (A) The relationship between secreted phosphoprotein 1 (SPP1) expression and six immune infiltrates in lung adenocarcinoma (LUAD) from the Tumor Immune Estimation Resource (TIMER). (B) The survival rate of SPP1 and the LUAD immune subsets from the TIMER.
SPP1 Is an Independent Prognostic Marker for LUAD Patients
The clinical relevance of LUAD immune subsets was explored by a multivariable Cox proportional hazards model, and SPP1 could be considered as an independent prognostic marker for LUAD patients (HR = 1.150; 95% CI, 1.048–1.261; P = 0.003; Table 2).
Table 2. Multivariate Cox proportional hazard regression analyses of the relationship between clinicopathological characteristics and overall survival in LUAD from the Tumor Immune Estimation Resource (TIMER).
Gene Set Enrichment Analysis
The effect of SPP1 expression on biological pathway was analyzed by GSEA based on GSE68465. The results of GSEA indicated that samples with high SPP1 expression enriched in protein secretion, mTORC1 signaling, angiogenesis, and glycolysis pathway (Table 3 and Figure 5).
Figure 5. Gene set enrichment analysis (GSEA) showed that the samples with high secreted phosphoprotein 1 (SPP1) expression enriched in (A) protein secretion, (B) mTORC1 signaling, (C) angiogenesis, and (D) glycolysis pathway.
Validation in TCGA Dataset
From the TCGA-LUAD data, the SPP1 mRNA expression in normal lung tissue was significantly lower than that in LUAD samples (P < 0.001; Figure 6A). Based on the clinical follow-up data of GSE68465, the overall 5-year survival rates of the SPP1 high-expression group and the SPP1 low-expression group were 33.755 and 45.329%, respectively (HR = 0.7088; 95% CI, 0.5298–0.9483; log-rank P = 0.020; Figure 6B). According to the clinical and pathological information from TCGA-LUAD, the SPP1 expression in LUAD was related to the N stage. As the N stage progresses, the expression level of SPP1 is higher than before (P < 0.05; Supplementary Table S2). The analysis results obtained using TCGA-LUAD data were basically consistent with the results obtained using GSE68465.
Figure 6. (A) The messenger RNA (mRNA) expression of secreted phosphoprotein 1 (SPP1) in lung adenocarcinoma (LUAD) and normal lung tissues based on The Cancer Genome Atlas (TCGA)-LUAD. (B) The overall survival rates of the SPP1 high-expression group and the SPP1 low-expression group in LUAD based on TCGA-LUAD.
SPP1 is a secreted calcium binding phosphorylation protein, which can activate many kinds of signaling pathways, promote tumor growth, metastasis, and angiogenesis and regulate the expression of various carcinogenic and angiogenic molecules. More and more studies have shown that SPP1 is closely related to the migration and metastasis of the malignant tumors. In multiple malignant tumors such as gastric, esophageal, and colorectal cancers, SPP1 could be detected to have a significantly high expression. Moreover, the SPP1 expression level was significantly higher in tumors with higher malignancy than that in tumors with lower malignancy (Lin et al., 2015; Xu et al., 2017; Chen et al., 2018; Assidi et al., 2019). Xu et al. (2017) found that SPP1 could activate the epithelial–mesenchymal transition (EMT) pathway to promote the metastasis of colorectal cancer. Chen et al. (2018) detected the potential functional marker single nucleotide polymorphism (tagSNP) of SPP1 by PCR and found that the frequency of SPP1 rs4754 genotype was significantly different from that of the control group, and the rs4754 TT genotype increased the risk of gastric cancer. Qin et al. (2018) found that SPP1 was upregulated in tumor tissues and plasma of patients with head and neck cancer, and the overexpression of SPP1 or the elevated level of plasma SPP1 in head and neck cancer tissues was associated with high malignancy and poor prognosis. These studies showed that SPP1 played a prominent role in the development of multiple tumors.
Giopanou et al. (2019) found that loss of SPP1 expression could be protective for mice harboring KRASG12D-driven LUAD; however, overexpression of SPP1 could promote the occurrence of early tumors and cause tumor-related inflammation. Zhang et al. (2017) detected that SPP1 was highly expressed in tumor-associated macrophages (TAMs) and LUAD tissues, and SPP1 played a role in lung cancer escape and mediating macrophage polarization. Wang et al. (2019) detected that inhibiting the SPP1 expression in lung cancer may overcome the resistance of second-generation EGFR TKI, and the expression of SPP1 may lead to increased resistance of lung cancer. These studies indicated that there was a close relationship with SPP1 and LUAD; SPP1 affecting the mechanism of lung adenocarcinoma needs to be further explored.
In this study, bioinformatics analysis revealed that SPP1 was highly expressed in LUAD tissues, which was related to poor prognosis. At the cell line level, NSCLC ranked ninth among cancer cell lines based on SPP1 expression. At the mRNA and protein level, we confirmed that SPP1 expression in LUAD tissues was higher than that in normal control based on the data from CCLE, GEO and HPA databases. We further verified the difference in SPP1 between LUAD and normal lung tissues, and the prognosis of LUAD samples by GSE68465. We found that the expression of SPP1 was related to gender, N stage, histological grade, and progression and relapse. Moreover, age, gender, T stage, N stage, and SPP1 had an impact on overall survival for LUAD patients. Then, we found that, compared to the SPP1 low-expression group, the overall survival rate of SPP1 high-expression group was significantly lower. The relationship between SPP1 expression and CD4+ T cell, macrophage, neutrophil, and dendritic cell infiltration was weak in LUAD from the TIMER database. Meanwhile, SPP1 could be considered as an independent prognostic marker for LUAD by a multivariable Cox proportional hazards model. The results indicated that samples with high SPP1 expression enriched in protein secretion, mTORC1 signaling, angiogenesis, and glycolysis pathway by GSEA based on GSE68465. The analysis results obtained using TCGA-LUAD data were basically consistent with the results obtained using GSE68465.
In the lung cancer microenvironment, tumor-infiltrating immune cells can not only attack and kill lung cancer cells to inhibit tumor progression but also screen tumor cells that were more suitable for survival in immune-active hosts or change the tumor microenvironment and promote tumor progression (Merlo et al., 2016). CD4+ T cells in the tumor microenvironment were divided into two subgroups of Th1 cells and Th2 cells according to the phenotype. The interaction between lung cancer cells, Th1 cells, and Th2 cells led to an imbalance of Th1/Th2 cell proportion, increasing the risk of lung cancer progression (Bremnes et al., 2016). Activated Th1 cells released cytokines such as tumor necrosis factor alpha (TNF-α), interferon gamma (IFN-γ), and interleukin (IL)-2, which can inhibit tumor progression by mediating cellular immunity, inducing apoptosis of tumor cells, and inhibiting the formation of neovascularization in tumors, while activated Th2 cells released IL-4, IL-5, IL-10, and IL-13 and other cytokines, through the mediator fluid immunity, and played a role in promoting tumor growth (Huang et al., 2017; Mateu-Jimenez et al., 2017). The M2 phenotype TAMs participated in complex autocrine and paracrine pathways and promoted tumor progression (Li et al., 2018). Dendritic cells, as powerful antigen-presenting cells, were the central link for initiating, regulating, and maintaining immune responses, and played a key role in inducing antilung cancer immune responses (Sabado et al., 2017). Zhang et al. (2017) detected that SPP1 was highly expressed in TAMs and LUAD tissues. Based on the above researches, the mechanism by which SPP1 promotes the development of LUAD may be closely related to changes in the immune microenvironment.
SPP1 may be involved in the secretion of exosomes of tumor cells. Tumor cells exchanged materials and information with the surrounding microenvironment by releasing a large number of exosomes, microbubbles or microparticles, and apoptotic bodies, which participated in promoting tumor growth and progression (Martins et al., 2013; McAllister and Weinberg, 2014). Excessive activation of mTOR signaling pathway was ubiquitous in tumor cells. mTOR activation played an important role in tumor cell growth, protein synthesis, cell metabolism, and other physiological and pathological processes, and its over activation could lead to the proliferation of a variety of cancer cells (Gweon and Kim, 2013; Martini et al., 2014; Safdari et al., 2015). Park et al. (2016) confirmed that resveratrol could induce autophagy of breast cancer cell MCF-7 by inhibiting the expression of mTOR. Because of the rapid proliferation of tumor cells, the microenvironment was always in a relatively hypoxic state, which was the basis of pathological angiogenesis. Neovascularization provided oxygen and nutrition for tumor growth and transported and removed its metabolic waste, thus supporting tumor survival (Hattori et al., 2011). Tumor angiogenesis was the key to tumor growth, invasion, and metastasis (Folkman, 1971; Jain, 2013). Warburg (1928) found that normal cells chose oxidative phosphorylation to provide ATP, while tumor cells preferred glycolysis; even when oxygen was sufficient, tumor cells still provided energy through glycolysis. The change in glycometabolism pathway was one of the differences between tumor cells and normal cells. Tumor cells could produce more nucleotides, fatty acids, proteins, and ATP through enhanced aerobic glycolysis, which provided a material basis for the rapid proliferation of tumor cells (Lunt and Vander, 2011). At the same time, it could reduce the production of reactive oxygen species, improve the antioxidant capacity of cells, and reduce the apoptosis of cells. In addition, aerobic glycolysis could produce a lot of lactic acid, creating an acidic microenvironment for tumor cells, which was conducive to the invasion and metastasis of tumor cells (Lu et al., 2015). Based on the above researches, SPP1 may be involved in protein secretion, mTORC1 signaling, angiogenesis, and glycolysis pathways, which in turn affect the occurrence and development of LUAD.
Although there is no research confirmed at present, according to the results of bioinformatics analysis, it can be concluded that SPP1 has the potential to be a new target for molecular targeted therapy of LUAD. However, due to the lack of necessary experimental conditions, the specific mechanism of SPP1 in the development of LUAD has not been thoroughly explored. The research group will continue to explore the specific role of SPP1 in the development of LUAD after the experimental conditions are perfect.
To sum up, we found that SPP1 was highly expressed in LUAD, which was closely related to the clinical characteristics of LUAD patients. Compared with the SPP1 low-expression samples, the overall survival rate of the SPP1 high-expression samples was significantly lower. At the same time, we explored the relationship between SPP1 and immune infiltration and elucidated that SPP1 enriched in protein secretion, mTORC1 signaling, angiogenesis, and glycolysis pathway. Therefore, SPP1 can not only affect the occurrence and development of LUAD but also be an independent prognostic marker of LUAD. SPP1 is expected to become a new target for molecular targeted therapy.
Data Availability Statement
All datasets generated/analyzed for this study are included/have their accession numbers included in the article/Supplementary Material.
WH and SL conceived and designed the study. ZG, JH, and YW performed the analysis procedures. JH, ZG, WL, and X-PL analyzed the results. JY, WH, and SL contributed to the analysis tools. ZG and JH contributed to the writing of the manuscript. All authors have reviewed the manuscript.
The present study was supported by the 351 Talent Project of Wuhan University (Luojia Young Scholars: SL), the Zhongnan Hospital of Wuhan University Science, Technology and Innovation Cultivating Fund (grant no. cxpy2017041), and Key Projects of Hubei provincial Health and Family Planning Commission (grant no. WJ2017Z006).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2020.00547/full#supplementary-material
- ^ https://portals.broadinstitute.org/ccle/
- ^ http://www.ncbi.nlm.nih.gov/geo/
- ^ http://www.proteinatlas.org/
- ^ https://cistrome.shinyapps.io/timer/
- ^ https://portal.gdc.cancer.gov/
Assidi, M., Gomaa, W., Jafri, M., Hanbazazh, M., Al-Ahwal, M., Pushparaj, P., et al. (2019). Prognostic value of Osteopontin (SPP1) in colorectal carcinoma requires a personalized molecular approach. Tumour. Biol. 41:1391226037. doi: 10.1177/1010428319863627
Barretina, J., Caponigro, G., Stransky, N., Venkatesan, K., Margolin, A. A., Kim, S., et al. (2012). The cancer cell line encyclopedia enables predictive modelling of anticancer drug sensitivity. Nature 483, 603–607. doi: 10.1038/nature11003
Berardi, R., Rinaldi, S., Santoni, M., Newsom-Davis, T., Tiberi, M., Morgese, F., et al. (2016). Prognostic models to predict survival in patients with advanced non-small cell lung cancer treated with first-line chemo- or targeted therapy. Oncotarget 7, 26916–26924. doi: 10.18632/oncotarget.8309
Bhattacharjee, A., Richards, W. G., Staunton, J., Li, C., Monti, S., Vasa, P., et al. (2001). Classification of human lung carcinomas by mRNA expression profiling reveals distinct adenocarcinoma subclasses. Proc. Natl. Acad. Sci. U.S.A. 98, 13790–13795. doi: 10.1073/pnas.191502998
Bray, F., Ferlay, J., Soerjomataram, I., Siegel, R. L., Torre, L. A., and Jemal, A. (2018). Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 68, 394–424. doi: 10.3322/caac.21492
Bremnes, R. M., Busund, L. T., Kilvaer, T. L., Andersen, S., Richardsen, E., and Paulsen, E. E. (2016). The role of tumor-infiltrating lymphocytes in development, progression, and prognosis of non-small cell lung cancer. J. Thorac. Oncol. 11, 789–800. doi: 10.1016/j.jtho.2016.01.015
Chen, C. H., Lai, J. M., Chou, T. Y., Chen, C. Y., Su, L. J., Lee, Y. C., et al. (2009). VEGFA upregulates FLJ10540 and modulates migration and invasion of lung cancer via PI3K/AKT pathway. PLoS One 4:e5052. doi: 10.1371/journal.pone.0005052
Chen, L. Z., He, C. Y., Su, X., Peng, J. L., Chen, D. L., Ye, Z., et al. (2018). SPP1 rs4754 and its epistatic interactions with SPARC polymorphisms in gastric cancer susceptibility. Gene 640, 43–50. doi: 10.1016/j.gene.2017.09.053
Fedarko, N. S., Jain, A., Karadag, A., Van Eman, M. R., and Fisher, L. W. (2001). Elevated serum bone sialoprotein and osteopontin in colon, breast, prostate, and lung cancer. Clin. Cancer Res. 7, 4060–4066.
Franceschini, J. P., Jamnik, S., and Santoro, I. L. (2017). Survival in a cohort of patients with lung cancer: the role of age and gender in prognosis. J. Bras. Pneumol. 43, 431–436. doi: 10.1590/S1806-37562016000000298
Giopanou, I., Kanellakis, N. I., Giannou, A. D., Lilis, I., Marazioti, A., Spella, M., et al. (2019). Osteopontin drives KRAS-mutant lung adenocarcinoma. Carcinogenesis. doi: 10.1093/carcin/bgz190 [Epub ahead of print].
Guo, Q., Ke, X. X., Liu, Z., Gao, W. L., Fang, S. X., Chen, C., et al. (2020). Evaluation of the prognostic value of STEAP1 in lung adenocarcinoma and insights into its potential molecular pathways via bioinformatic analysis. Front Genet 11:242. doi: 10.3389/fgene.2020.00242
Gweon, E. J., and Kim, S. J. (2013). Resveratrol induces MMP-9 and cell migration via the p38 kinase and PI-3K pathways in HT1080 human fibrosarcoma cells. Oncol. Rep. 29, 826–834. doi: 10.3892/or.2012.2151
Hattori, N., Okochi-Takada, E., Kikuyama, M., Wakabayashi, M., Yamashita, S., and Ushijima, T. (2011). Methylation silencing of angiopoietin-like 4 in rat and human mammary carcinomas. Cancer Sci. 102, 1337–1343. doi: 10.1111/j.1349-7006.2011.01955.x
Hou, J., Aerts, J., den Hamer, B., van Ijcken, W., den Bakker, M., Riegman, P., et al. (2010). Gene expression-based classification of non-small cell lung carcinomas and survival prediction. PLoS One 5:e10312. doi: 10.1371/journal.pone.0010312
Huang, J., Shen, F., Huang, H., Ling, C., and Zhang, G. (2017). Th1high in tumor microenvironment is an indicator of poor prognosis for patients with NSCLC. Oncotarget 8, 13116–13125. doi: 10.18632/oncotarget.14471
Kijewska, M., Kocyk, M., Kloss, M., Stepniak, K., Korwek, Z., Polakowska, R., et al. (2017). The embryonic type of SPP1 transcriptional regulation is re-activated in glioblastoma. Oncotarget 8, 16340–16355. doi: 10.18632/oncotarget.14092
Landi, M. T., Dracheva, T., Rotunno, M., Figueroa, J. D., Liu, H., Dasgupta, A., et al. (2008). Gene expression signature of cigarette smoking and its role in lung adenocarcinoma development and survival. PLoS One 3:e1651. doi: 10.1371/journal.pone.0001651
Li, B., Severson, E., Pignon, J. C., Zhao, H., Li, T., Novak, J., et al. (2016). Comprehensive analyses of tumor immunity: implications for cancer immunotherapy. Genome. Biol. 17:174. doi: 10.1186/s13059-016-1028-7
Li, H., Huang, N., Zhu, W., Wu, J., Yang, X., Teng, W., et al. (2018). Modulation the crosstalk between tumor-associated macrophages and non-small cell lung cancer to inhibit tumor migration and invasion by ginsenoside Rh2. BMC Cancer 18:579. doi: 10.1186/s12885-018-4299-4
Li, T., Fan, J., Wang, B., Traugh, N., Chen, Q., Liu, J. S., et al. (2017). TIMER: a web server for comprehensive analysis of Tumor-Infiltrating immune cells. Cancer Res. 77, e108–e110. doi: 10.1158/0008-5472.CAN-17-0307
Liaw, L., Birk, D. E., Ballas, C. B., Whitsitt, J. S., Davidson, J. M., and Hogan, B. L. (1998). Altered wound healing in mice lacking a functional osteopontin gene (spp1). J. Clin. Invest. 101, 1468–1478. doi: 10.1172/JCI2131
Lin, J., Myers, A. L., Wang, Z., Nancarrow, D. J., Ferrer-Torres, D., Handlogten, A., et al. (2015). Osteopontin (OPN/SPP1) isoforms collectively enhance tumor cell invasion and dissemination in esophageal adenocarcinoma. Oncotarget 6, 22239–22257. doi: 10.18632/oncotarget.4161
Lu, J., Tan, M., and Cai, Q. (2015). The warburg effect in tumor progression: mitochondrial oxidative metabolism as an anti-metastasis mechanism. Cancer Lett. 356, 156–164. doi: 10.1016/j.canlet.2014.04.001
Mateu-Jimenez, M., Curull, V., Pijuan, L., Sanchez-Font, A., Rivera-Ramos, H., Rodríguez-Fuster, A., et al. (2017). Systemic and tumor th1 and th2 inflammatory profile and macrophages in lung cancer: influence of underlying chronic respiratory disease. J. Thorac. Oncol. 12, 235–248. doi: 10.1016/j.jtho.2016.09.137
Mootha, V. K., Lindgren, C. M., Eriksson, K. F., Subramanian, A., Sihag, S., Lehar, J., et al. (2003). PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat. Genet 34, 267–273. doi: 10.1038/ng1180
Okayama, H., Kohno, T., Ishii, Y., Shimada, Y., Shiraishi, K., Iwakawa, R., et al. (2012). Identification of genes upregulated in ALK-positive and EGFR/KRAS/ALK-negative lung adenocarcinomas. Cancer Res. 72, 100–111. doi: 10.1158/0008-5472.CAN-11-1403
Park, D., Jeong, H., Lee, M. N., Koh, A., Kwon, O., Yang, Y. R., et al. (2016). Resveratrol induces autophagy by directly inhibiting mTOR through ATP competition. Sci. Rep. 6:21772. doi: 10.1038/srep21772
Qin, X., Yan, M., Wang, X., Xu, Q., Wang, X., Zhu, X., et al. (2018). Cancer-associated fibroblast-derived IL-6 promotes head and neck cancer progression via the Osteopontin-NF-kappa b signaling pathway. Theranostics 8, 921–940. doi: 10.7150/thno.22182
Safdari, Y., Khalili, M., Ebrahimzadeh, M. A., Yazdani, Y., and Farajnia, S. (2015). Natural inhibitors of PI3K/AKT signaling in breast cancer: emphasis on newly-discovered molecular mechanisms of action. Pharmacol. Res. 93, 1–10. doi: 10.1016/j.phrs.2014.12.004
Selamat, S. A., Chung, B. S., Girard, L., Zhang, W., Zhang, Y., Campan, M., et al. (2012). Genome-scale analysis of DNA methylation in lung adenocarcinoma and integration with mRNA expression. Genome. Res. 22, 1197–1211. doi: 10.1101/gr.132662.111
Shedden, K., Taylor, J. M., Enkemann, S. A., Tsao, M. S., Yeatman, T. J., Gerald, W. L., et al. (2008). Gene expression-based survival prediction in lung adenocarcinoma: a multi-site, blinded validation study. Nat. Med. 14, 822–827. doi: 10.1038/nm.1790
Stearman, R. S., Dwyer-Nield, L., Zerbe, L., Blaine, S. A., Chan, Z., and Bunn, P. A. Jr., et al. (2005). Analysis of orthologous gene expression between human pulmonary adenocarcinoma and a carcinogen-induced murine model. Am. J. Pathol 167, 1763–1775. doi: 10.1016/S0002-9440(10)61257-6
Stone, G. W., Maehara, A., Lansky, A. J., de Bruyne, B., Cristea, E., Mintz, G. S., et al. (2011). A prospective natural-history study of coronary atherosclerosis. N. Engl. J. Med. 364, 226–235. doi: 10.1056/NEJMoa1002358
Su, L. J., Chang, C. W., Wu, Y. C., Chen, K. C., Lin, C. J., Liang, S. C., et al. (2007). Selection of DDX5 as a novel internal control for Q-RT-PCR from microarray data using a block bootstrap re-sampling scheme. BMC Genomics 8:140. doi: 10.1186/1471-2164-8-140
Subramanian, A., Tamayo, P., Mootha, V. K., Mukherjee, S., Ebert, B. L., Gillette, M. A., et al. (2005). Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. U.S.A. 102, 15545–15550. doi: 10.1073/pnas.0506580102
Tang, J., Kong, D., Cui, Q., Wang, K., Zhang, D., Gong, Y., et al. (2018). Prognostic genes of breast cancer identified by gene co-expression network analysis. Front. Oncol. 8:374. doi: 10.3389/fonc.2018.00374
Tian, T., Wang, M., Zhu, Y., Zhu, W., Yang, T., Li, H., et al. (2018). Expression, clinical significance, and functional prediction of MNX1 in breast cancer. Mol. Ther. Nucleic Acids 13, 399–406. doi: 10.1016/j.omtn.2018.09.014
Uhlen, M., Bjorling, E., Agaton, C., Szigyarto, C. A., Amini, B., Andersen, E., et al. (2005). A human protein atlas for normal and cancer tissues based on antibody proteomics. Mol. Cell Proteomics 4, 1920–1932. doi: 10.1074/mcp.M500279-MCP200
Uhlen, M., Fagerberg, L., Hallstrom, B. M., Lindskog, C., Oksvold, P., Mardinoglu, A., et al. (2015). Proteomics. tissue-based map of the human proteome. Science 347:1260419. doi: 10.1126/science.1260419
Uhlen, M., Oksvold, P., Fagerberg, L., Lundberg, E., Jonasson, K., Forsberg, M., et al. (2010). Towards a knowledge-based human protein atlas. Nat. Biotechnol. 28, 1248–1250. doi: 10.1038/nbt1210-1248
Wang, X., Zhang, F., Yang, X., Xue, M., Li, X., Gao, Y., et al. (2019). Secreted phosphoprotein 1 (SPP1) contributes to Second-Generation EGFR tyrosine kinase inhibitor resistance in Non-Small cell lung cancer. Oncol. Res. 27, 871–877. doi: 10.3727/096504018X15426271404407
Xu, C., Sun, L., Jiang, C., Zhou, H., Gu, L., Liu, Y., et al. (2017). SPP1, analyzed by bioinformatics methods, promotes the metastasis in colorectal cancer by activating EMT pathway. Biomed. Pharmacother. 91, 1167–1177. doi: 10.1016/j.biopha.2017.05.056
Yamauchi, M., Yamaguchi, R., Nakata, A., Kohno, T., Nagasaki, M., Shimamura, T., et al. (2012). Epidermal growth factor receptor tyrosine kinase defines critical prognostic genes of stage I lung adenocarcinoma. PLoS One 7:e43923. doi: 10.1371/journal.pone.0043923
Zhang, W., Fan, J., Chen, Q., Lei, C., Qiao, B., Liu, Q., et al. (2018). SPP1 and AGER as potential prognostic biomarkers for lung adenocarcinoma. Oncol. Lett. 15, 7028–7036. doi: 10.3892/ol.2018.8235
Keywords: secreted phosphoprotein 1, lung adenocarcinoma, prognosis, prognostic marker, osteopontin, bioinformatics
Citation: Guo Z, Huang J, Wang Y, Liu X-P, Li W, Yao J, Li S and Hu W (2020) Analysis of Expression and Its Clinical Significance of the Secreted Phosphoprotein 1 in Lung Adenocarcinoma. Front. Genet. 11:547. doi: 10.3389/fgene.2020.00547
Received: 15 February 2020; Accepted: 07 May 2020;
Published: 12 June 2020.
Edited by:Pinyi Lu, Biotechnology HPC Software Applications Institute (BHSAI), United States
Reviewed by:Xu Zhang, Southwest University, China
Jianxiang Li, Medical College of Soochow University, China
Yang Xu, Zhejiang University, China
Copyright © 2020 Guo, Huang, Wang, Liu, Li, Yao, Li and Hu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
†These authors share first authorship