Analysis of Expression and Its Clinical Significance of the Secreted Phosphoprotein 1 in Lung Adenocarcinoma

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.

Introduction

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 (CCLE¹) 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 (GEO²) database was employed for analyses. The mRNA expression levels were shown as the intensity of log₂ 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)³ 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 log₂ transformation were conducted.

TIMER Analysis

Tumor Immune Estimation Resource (TIMER)⁴ 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 Portal⁵, 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 log₂ 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.

Statistical Analysis

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 χ² 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.

Results

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

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

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

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

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

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

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).

TABLE 3

Table 3. Gene set enrichment in SPP1 high-expression phenotype (based on GSE68465).

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

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.

Discussion

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.

Conclusion

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.

Author Contributions

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.

Funding

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.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2020.00547/full#supplementary-material

Footnotes

1. ^ https://portals.broadinstitute.org/ccle/
2. ^ http://www.ncbi.nlm.nih.gov/geo/
3. ^ http://www.proteinatlas.org/
4. ^ https://cistrome.shinyapps.io/timer/
5. ^ https://portal.gdc.cancer.gov/

References

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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.

Google Scholar

Folkman, J. (1971). Tumor angiogenesis: therapeutic implications. N. Engl. J. Med. 285, 1182–1186. doi: 10.1056/NEJM197111182852108

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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].

CrossRef Full Text | PubMed Abstract | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Jain, R. K. (2013). Normalizing tumor microenvironment to treat cancer: bench to bedside to biomarkers. J. Clin. Oncol. 31, 2205–2218. doi: 10.1200/JCO.2012.46.3653

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnson, J. L., Pillai, S., and Chellappan, S. P. (2012). Genetic and biochemical alterations in non-small cell lung cancer. Biochem. Res. Int. 2012:940405. doi: 10.1155/2012/940405

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Lam, E. W., Brosens, J. J., Gomes, A. R., and Koo, C. Y. (2013). Forkhead box proteins: tuning forks for transcriptional harmony. Nat Rev Cancer 13, 482–495. doi: 10.1038/nrc3539

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y., Du, W., Han, J., and Ge, J. (2017). LAMP3 promotes the invasion of osteosarcoma cells via SPP1 signaling. Mol Med Rep 16, 5947–5953. doi: 10.3892/mmr.2017.7349

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Lunt, S. Y., and Vander, H. M. (2011). Aerobic glycolysis: meeting the metabolic requirements of cell proliferation. Annu. Rev. Cell Dev. Biol. 27, 441–464. doi: 10.1146/annurev-cellbio-092910-154237

PubMed Abstract | CrossRef Full Text | Google Scholar

Martini, M., De Santis, M. C., Braccini, L., Gulluni, F., and Hirsch, E. (2014). PI3K/AKT signaling pathway and cancer: an updated review. Ann. Med. 46, 372–383. doi: 10.3109/07853890.2014.912836

PubMed Abstract | CrossRef Full Text | Google Scholar

Martins, V. R., Dias, M. S., and Hainaut, P. (2013). Tumor-cell-derived microvesicles as carriers of molecular information in cancer. Curr. Opin. Oncol. 25, 66–75. doi: 10.1097/CCO.0b013e32835b7c81

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

McAllister, S. S., and Weinberg, R. A. (2014). The tumour-induced systemic environment as a critical regulator of cancer progression and metastasis. Nat. Cell Biol. 16, 717–727. doi: 10.1038/ncb3015

PubMed Abstract | CrossRef Full Text | Google Scholar

Merlo, A., Dalla, S. S., Dolcetti, R., Zanovello, P., and Rosato, A. (2016). Reverse immunoediting: when immunity is edited by antigen. Immunol. Lett. 175, 16–20. doi: 10.1016/j.imlet.2016.04.015

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Noguchi, M., Morikawa, A., Kawasaki, M., Matsuno, Y., Yamada, T., and Hirohashi, S. (1995). Small adenocarcinoma of the lung. Histologic characteristics and prognosis. Cancer 75, 2844–2852.

Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Sabado, R. L., Balan, S., and Bhardwaj, N. (2017). Dendritic cell-based immunotherapy. Cell Res. 27, 74–95. doi: 10.1038/cr.2016.157

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Siegel, R. L., Miller, K. D., and Jemal, A. (2015). Cancer statistics, 2015. CA Cancer J. Clin. 65, 5–29. doi: 10.3322/caac.21254

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Warburg, O. (1928). The chemical constitution of respiration ferment. Science 68, 437–443. doi: 10.1126/science.68.1767.437

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, R., Wong, J., and Kwok, H. F. (2017). Osteopontin – a promising biomarker for cancer therapy. J. Cancer 8, 2173–2183. doi: 10.7150/jca.20480

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Zeng, B., Zhou, M., Wu, H., and Xiong, Z. (2018). SPP1 promotes ovarian cancer progression via Integrin beta1/FAK/AKT signaling pathway. Onco. Targets Ther. 11, 1333–1343. doi: 10.2147/OTT.S154215

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Du, W., Chen, Z., and Xiang, C. (2017). Upregulation of PD-L1 by SPP1 mediates macrophage polarization and facilitates immune escape in lung adenocarcinoma. Exp. Cell. Res. 359, 449–457. doi: 10.1016/j.yexcr.2017.08.028

PubMed Abstract | CrossRef Full Text | Google Scholar