Skip to main content

SNX10 and PTGDS are associated with the progression and prognosis of cervical squamous cell carcinoma



Cervical cancer (CC) is the primary cause of death in women. This study sought to investigate the potential mechanism and prognostic genes of CC.


We downloaded four gene expression profiles from GEO. The RRA method was used to integrate and screen differentially expressed genes (DEGs) between CC and normal samples. Functional analysis was performed by clusterprofiler. We built PPI network by Search Tool for the Retrieval of Interacting Genes Database (STRING) and selected hub modules via Molecular COmplex Detection (MCODE). CMap database was used to find molecules with therapeutic potential for CC. The hub genes were validated in GEO datasets, Gene Expession Profiling Interactive Analysis (GEPIA), immunohistochemistry, Cox regression analysis, TCGA methylation analysis and ONCOMINE were carried out. ROC curve analysis and GSEA were also performed to describe the prognostic significance of hub genes.


Functional analysis revealed that 147 DEGs were significantly enriched in binding, cell proliferation, transcriptional activity and cell cycle regulation. PPI network screened 30 hub genes, with CDK1 having the strongest connectivity with CC. Cmap showed that apigenin, thioguanine and trichostatin A might be used to treat CC(P < 0.05). Eight genes (APOD, CXCL8, MMP1, MMP3, PLOD2, PTGDS, SNX10 and SPP1) were screened out through GEPIA. Of them, only PTGDS and SNX10 had not appeared in previous studies about CC. The validation in GEO showed that PTGDS showed low expression while SNX10 presented high expression in tumor tissues. Their expression profiles were consistent with the results in immunohistochemistry. ROC curve analysis indicated that the model had a good diagnostic efficiency (AUC = 0.738). GSEA analysis demonstrated that the two genes were correlated with the chemokine signaling pathway (P < 0.05). TCGA methylation analysis showed that patients with lowly-expressed and highly-methylated PTGDS had a worse prognosis than those with highly-expressed and lowly-methylated PTGDS (p = 0.037). Cox regression analysis showed that SNX10 and PTGDS were independent prognostic indicators for OS among CC patients (P = 0.007 and 0.003).


PTGDS and SNX10 showed abnormal expression and methylation in CC. Both genes might have high prognostic value of CC patients.

Peer Review reports


An annual death toll of 265,700 makes cervical cancer (CC) the second deadliest malignancy in females [1]. Despite pre-cancerous screening and emerging treatments, CC remains the primary cause of death in women in developing countries [2]. When CC metastasizes and recurs, the prognosis gets even worse. Therefore, it is of great significance to create new treatments of CC based on its to-be-clarified molecular mechanism.

Gene expression microarray, as an efficient means of acquiring large-scale genetic data, is being widely used to study gene expression profiling in many human cancers. Upon microarray and databases, effective analytic tools have been designed to explore tumor-associated genes, molecular mechanisms and target therapies. The integration of databases containing gene expression chips allows in-depth study of molecular mechanisms [3, 4].

Thousands of differentially expressed genes (DEGs) in CC have been discovered [5,6,7]. However, the results on some mRNAs are inconsistent. Here we use an unbiased approach to solve this problem.

In our study, we screened DEGs from four profiles downloaded from GEO. PPI network was built by STRING Database and hub modules selected via plug-in MCODE. CMap was used to find potential genes associated with CC. We also validated hub genes with GEO, GEPIA, immunohistochemistry and ONCOMINE. ROC curve analysis and GESA were also done to tease out the significance of hub genes. The flow chart of this research was displayed in Fig. 1.

Fig. 1
figure 1

Flow chart of the present study


Screening DEGs

Keywords “cervical cancer geo accession” were put in the GEO database ( and the mRNA expression profiles of GSE6791, GSE63514, GSE39001 and GSE9750 were downloaded. The dataset details were shown in Table 1. We processed unqualified data by R package. The data is calibrated, standardized and log2-transformed. Gene expression analysis was performed using the “limma” R package [8] in the Bioconductor package. Relevant codes were placed into R. we selected four microarray datasets and analyzed them with limma. The |log2fold change (FC)| > 2 and adjusted p < 0.05 were set as cutoffs. RRA package was download ( [9] and R software was used for running the instruction code.

Table 1 Details for GEO cervical cancer data

Functional analysis based on DEGs

The biological function of DEGs was analyzed with DAVID ( database and clusterprofiler [10] (a package visualizing the biological profiles of genes). P < 0.05 was considered to be statistically significant.

PPI network integration

Search Tool for the Retrieval of Interacting Genes Database (STRING) [11] ( was used to assess PPI complex between identified and predicted proteins. In addition, the plug-in MCODE [12] of Cytoscape was conducted to select and visualize hub modules in PPI complex.

Identification of potential drugs

CMap [13] is a computer simulation method for predicting the potential drug that may induce or reverse a biological state encoded by the gene expression signature. The different probe components commonly between CC and normal samples were screened out with CMap database and divided into the up- and down-regulated groups. An enrichment score representing similarity was calculated. The positive score illustrated that the drug could induce cancer in human; the negative score illustrated the drug function oppositely and had potential therapeutic value.

Construction of a prognostic signature

Univariate Cox regression analysis was performed based on DEGs. The genes associated with prognosis were defined using the cutoff value of p < 0.05. Next, a multivariate Cox regression model was constructed with genes of p < 0.01. Cox regression with p < 0.05 was constructed to estimate the risk score of each patient on the expression of the DEGs. To further screen out prognosis-related genes of CC, we constructed a linear risk model. The prognostic score = expRNA1 × βRNA1 + expRNA2 × βRNA2 + expRNA3 × βRNA3 + ...expRNAn×βRNAn (expRNA is the expression level of each methylation-driven gene, and βRNA is the regression coefficient calculated by the multivariate Cox regression analysis). The prognostic risk value of each sample was calculated according to the formula, and then the median of the index value was cut off. Patients were separated into a low- and high-risk group according to their mean scores of prognostic risks. Kaplan-Meier survival analysis was conducted based on the low- and high-risk group. We also performed ROC curve analysis using 5 years as the predicted time to assess the predictive value of the outcome. The areas under the ROC curve, sensitivity and specificity were used to describe predictive values.

Validation of key genes

We used GEPIA [14] (Gene Expression Profiling Interactive Analysis) to analyze the expression and prognostic significance of DEGs. After reviewing literature, we screened real hub genes from DEGs. Subsequently the real hub genes were validated into GEO datasets, including GSE7803 and GSE29576, and ONCOMINE database ( GSE7803 included 21 cervical cancer samples and 10 normal cervix tissue samples. GSE29576 included 45 CC samples and 17 normal cervix tissue samples. ONCOMINE [15] dataset is a public online cancer microarray database that enables online analysis on relations between certain gene and various tumors according to DNA and RNA sequence data. The Human Protein Atlas (HPA) ( was also used to measure the expressing level of the real hub genes. ROC curve analysis was performed to distinguish normal and cancer tissues.

Survival analysis and mapping of methylation level

Survival analysis on gene methylation and expression was conducted through R package to identify key prognosis-associated genes in CC. To explore the relation between aberrant methylation and expression of genes, we extracted key genes with methylated expression from the downloaded data on TCGA CESC methylation. Then we evaluated the association between the methylated expression and the gene expression.

Gene set enrichment analysis (GSEA)

According to the hub gene expression level, the samples were then separated into two different groups. To explore the potential function of the DEGs, GSEA [16] ( was used to research a series of biological pathways that might be enriched in the gene rank derived from hub gene among the two groups. Annotated gene set of c2.cp.kegg.v6.0.symbols.gmt in Molecular Signatures Database (MSigDB, was selected as the reference. Additionally, we used “Clusterprofiler” package in R to handle the datasets, and the “Enrichplot” package to tease out the enriched pathways of the key genes. The adjusted-P < 0.05 was set as significant.


Identification of DEGs in CC

The CC expression microarray datasets (GSE6791, GSE9750, GSE39001 and GSE63514) were firstly standardized (Figure S1). With limma package, 256 DEGs were filtered from GSE6791 (60 downregulated and 196 upregulated), 236 DEGs from GSE9750 (136 downregulated and 100 upregulated), 98 DEGs from GSE39001 (38 upregulated and 60 downregulated), 489 DEGs from GSE63514 (177 upregulated and 312 downregulated). DEGs from the 4 microarray datasets were exhibited in volcano maps (Figure S2A-D) and heatmaps (Figure S3A-D). We analyzed the four microarray datasets via the limma package and then with RRA method according to their log-folding variation values ((|log2fold change (FC)| > 1 and adj. p < 0.05). The RRA method was based on the theory that genes in each experiment were randomly ordered. For the genes ranking higher in the experiment, the possibility of differential expression is inversely proportional to the value of P. Through analytic hierarchy analysis, we sorted out 74 up-regulated and 73 down-regulated genes (Table 2). Finally, the R-heatmap software was performed to plot the top 40 up- and down-regulated genes (Fig. 2).

Table 2 Screening DEGs in cervical cancer by integrated microarray
Fig. 2
figure 2

LogFC heatmap of the image data of each expression microarray. Notes: The abscissa is the geo ID, and the ordinate is the gene name. Red represents logFC>0, green represents logFC<0, and the values in the box represent the logFC values

Functional analysis of DEGs

The biological annotations of DEGs in CC were obtained with an online analysis tool named DAVID, which had GO analysis of up- and down-regulated genes (P<0.05). The GO analysis of DEGs covered three aspects: biological processes, molecular function and cellular composition (Figure S4A). The upregulated genes were significantly enriched in microtubule binding, tubulin binding and ATPase activity (Fig. 3a), and the down-regulated genes in serine-type peptidase activity, serine-type endopeptidase and serine hydrolase activity (Fig. 3b). These results indicated that most DEGs were prominently enriched in structural molecule activity, midbody, kinesin complex and microtubule motor activity. (Figure S4 B-C and Figure S5A). Clusterprofile was performed to analyze the DEGs. The result showed that the upregulated genes were mostly enriched in DNA replication, Oocyte meiosis and Cell cycle. (Fig. 3c), and the down-regulated genes in Arachidonic acid metabolism, prostate cancer and signaling pathways regulating pluripotency of stem cells (Fig. 3d). The pathway-gene network (Figure S5B) suggested that the cell cycle was the most important term in the biological processes of CC.

Fig. 3
figure 3

GO and KEGG analysis of the DEGs. a GO analysis of upregulated genes associated with cervical cancer. b GO analysis of downregulated genes associated with cervical cancer. c KEGG analysis of upregulated genes associated with cervical cancer. d KEGG analysis of downregulated genes associated with cervical cancer

PPI network construction and modules selection

The PPI network of DEGs was constructed, including 147 nodes (74 up-regulated and 73 down-regulated genes) and 562 edges (Fig. 4a). Degrees ≥30 was set as the cutoff. A total of 16 genes, such as CDK1, TOP2A, NCAPG, and KIF11, were showed the most significant difference in expression (Fig. 4b). A significant module was selected with plug-in MCODE, including 27 nodes and 343 edges (Figure S6A). GO and KEGG analysis indicated that the genes in the module were related to microtubule binding, tubulin binding, cell cycle and oocyte mitosis (Figure S6B and C).

Fig. 4
figure 4

PPI network analysis. a Using the STRING online database, a total of 147 DEGs were filtered into the DEGs PPI network. b Degree of the top 30 nodes in the PPI network. All these nodes are upregulated genes

Small molecule drugs screening

CMap network was used to analyze 147 DEGs into two groups (74 in up-regulated group and 73 in down-regulated group). After the signature query, the three compounds with the highest negative enrichment score (apigenin, thioguanine, and trichostatin A) were identified as potential therapeutic agents for CC (Table 3). The three-dimensional chemical structure of these three compounds is shown in Fig. 5.

Table 3 Results of CMap analysis
Fig. 5
figure 5

Three-dimensional diagram of the three most significant drugs. a Apigenin b Thioguanosine c Trichostatin A

Validation of hub genes

We validated DEGs at GEPIA website, including survival analysis and tissue sample expression analysis (Figure S7 and Figure S8). Eight genes (APOD, CXCL8, MMP1, MMP3, PLOD2, PTGDS, SNX10 and SPP1) had the same trend in both the above analysis. We literature-reviewed these eight genes, finding that only PTGDS and SNX10 had not been reported to be associated with CC. Therefore, we used GSE7803 and GSE29576 to validate PTGDS and SNX10 (Figure S9). The results showed that PTGDS had high expression levels in normal tissues and low expression levels in tumor tissues, while SNX10 showed an opposite profile. We further validated the two genes using immunohistochemistry (Fig. 6a-b) and ONCOMINE, obtaining the results consistent with those from the GEO database (Fig. 6c-d). The area under the curve of PTGDS was 0.919 and that of PTGDS was 0.905, suggesting that both can distinguish CC and normal tissue and have a good diagnostic efficiency (Fig. 7a). GSEA was performed to search KEGG pathways enriched in the TCGA samples. The top ten most enriched pathways included “hematopoietic cell lineage”, “adhesion molecules cams”, “vascular smooth muscle contraction”, “systemic lupus erythematosus”, “chemokine signaling pathway”, “t cell receptor signaling pathway”, “cytokine cytokine receptor interaction”, “calcium signaling pathway”, “neuroactive ligand receptor interaction” and “leukocyte transendothelial migration” (Fig. 7b) (adj.p < 0.05). In addition, the univariate and multivariate Cox regression analyses showed that SNX10 and PTGDS were independent prognostic indicators for OS among CESC patients (P = 0.007 and 0.003) (Table 4).

Fig. 6
figure 6

Validation of PTGDS and SNX10. a Immunohistochemistry of PTGDS. b Immunohistochemistry of SNX10. c Expression of PTGDS in ONCOMINE. d Expression of SNX10 in ONCOMINE

Fig. 7
figure 7

Validation of PTGDS and SNX10. a ROC curve analysis of the two genes. b GSEA of PTGDS and SNX10. c The methylated expression and gene expression of PTGDS. d The methylated expression and gene expression of SNX10. e Survival analysis of patients with methylated PTGDS expression

Table 4 Univariate analysis and multivariate analysis of SNX10 and PTGDS expression among cervical cancer patients

To find out the mechanism of abnormal gene expression, we analyzed the gene expression level and methylation level from the Illumina Human Methylation 450 platform based on TCGA data. The association between the methylated expression and the gene expression of the two key driving genes were shown in Fig. 7c-d. The survival analysis showed that the patients with low-expressed and hyper-methylated PTGDS had a worse prognosis than those with high-expressed and hypo-methylated PTGDS (P = 0.037) (Fig. 7e). However, SNX10 methylation has no statistical significance in survival analysis.

Establishment of cox regression model

Univariate cox regression analysis screened out seven genes, including APOD, CXCL8, MMP1, MMP3, PLOD2, PTGDS and SPP1 (Figure S10). Multivariate Cox regression analysis screened five genes, including SPP1, CXCL8, PTGDS, PLOD2 and MMP3 (Figure S11). The score for overall survival risk prediction was calculated as followed: Risk score = 0.143* SPP1 + 0.136* CXCL8–0.093* PTGDS+ 0.206* PLOD2 + 0.067* MMP3. Based on the risk score, CC patients were divided into low- and high-risk groups. Kaplan-Meier analysis suggested that low-risk patients had better outcomes than high-risk patients in the TCGA cohort (Fig. 8a). ROC curve analysis was also completed according to the 5-year survival of the area under the receiver operating characteristic curve (AUC) value. The specificity and sensitivity were both highest when the risk score was 0.738 (Fig. 8b). The distribution of risk score, survival status, and the expression levels of five genes was also presented (Fig. 8c-f). The expression levels of five genes in low- and high-risk groups were shown in Figure S12A. The univariate and multivariate Cox regression analyses showed that only the risk score based on five genes was independent prognostic indictor of CC (Figure S12B-C).

Fig. 8
figure 8

Survival prognosis model of the 5 hub genes. a Survival analysis showed that the patients in the high risk group had statistically significantly worse overall survival than those in low risk group in TCGA cohort. b ROC analysis was performed to find out the most optimal cutoff value to divide the CC patients into high risk and low risk group. c-d The risk scores for all patients in TCGA cohort are plotted in ascending order and marked as low risk (blue) or high risk (red), as divided by the threshold (vertical black line). e Eight expression and risk score distribution in TCGA cohort by z-score, with red indicating higher expression and light blue indicating lower expression

The Heatmap showed the expression levels of the five genes in two subgroup patients in the TCGA dataset. We observed significant difference in survival state (P < 0.001) and stage (P < 0.05) (Figure S12D).


CC brings on more than 265,700 deaths per year, making it the second deadlist malignancy in women. Nowadays, microarray and high-throughput sequencing technology are used to identify the potential targets in CC treatment. Previous studies often establish a single group or have a small-size, which restricts their reliability.

This study analyzed the expression profiles of four genes using R software and bioinformatics tools. A total of 147 differential genes were identified using RRA analysis, including 73 downregulated and 74 upregulated. GO analysis indicated that upregulated DEGs were associated with microtubule binding and downregulated DEGs with serine-type peptidase activity. KEGG analysis showed that these DEGs were primarily enriched in the cell cycle pathway.

Our findings echo with the previous studies. It has been reported that microtubule binding and cell cycle have an effect on CC [17]. Other studies showed that microtubule binding played a role in the biology of acute myeloid leukemia cells and colorectal cancer cells [18, 19]. Cell cycle also decides the abnormal proliferation of many tumor cells [20, 21].

PPI network displayed 30 hub genes associated with CC associated proteins. Next, we found CDK1 was in the center of PPI network of CESC. CDK1 was harbored by top module 1, suggesting that the top module 1 plays a crucial role in CC pathogenesis. Functional analysis indicated that the genes in this module were mainly enriched in microtubule binding, tubulin binding and cell cycle. Cyclin-dependent kinase 1 (CDK1) is a serine/threonine kinases that interacted with specific cyclins to regulate the cell cycle [22]. It has been reported that CDK1 regulated the development of CC and many other tumors [23]. Y. Zeng et al. found that the mitotic phosphorylation level of the transcriptional co-repressor Vgll4 was mediated by CDK1 to its tumor-suppressing activity. K. Bednarek et al. found that CDK1 was involved in the processes of laryngeal squamous cell carcinoma [24, 25].

Cmap showed that apigenin, thioguanine and trichostatin A could be used to treat CC. Apigenin and trichostatin A can inhibit breast cancer growth [26, 27]. 6-thioguanine has also potential therapeutic effects on tumors [28]. Our findings may help create the appropriate drugs for CC treatment.

Eight genes (APOD, CXCL8, MMP1, MMP3, PLOD2, PTGDS, SNX10 and SPP1) were screened out of DEGs through GEPIA. Of them, only PTGDS and SNX10 had not been reported in CC research. According to GEO validation results, PTGDS was lowly expressed and SNX10 highly expressed in tumor tissues, which was consistent with the results from immunohistochemistry. TCGA methylation analysis showed that the patients with lowly-expressed and highly-methylated PTGDS had a worse prognosis than those with highly-expressed and lowly-methylated PTGDS. Cox regression analysis showed that SNX10 and PTGDS were independent prognostic indicators for OS among CC patients. GSEA showed that the two genes were associated with the chemokine signaling pathway. Zhong G et al. suggested that chemokine signaling is involved in the invasion and migration of lung cancer cells [29]. Chemokine signaling has also been reported to maneuver the progression of breast and hepatobiliary cancer [30, 31]. In addition, the prognostic signature was constructed based on the eight hub genes. Of them, five genes (SPP1, CXCL8, PTGDS, PLOD2, and MMP3) exhibited significant prognostic value. The Cox regression analysis showed that only the risk score of the five genes was an independent prognostic indicator of CC.

Interestingly, all of the above genes are associated with cervical cancer. Yan R et al. found that CXCL8 had prognostic value in cervical carcinoma patients [32]. Tian R et al. identified the function of MMP1 in cervical cancer [33]. Xie B et al. defined that genetic polymorphisms in MMP 3 connected with the clinical outcome of cervical cancer in a Chinese Han population [34]. Xu F et al. found that increased PLOD2 expression facilitated epithelial-to-mesenchymal transition (EMT) and focal adhesion formation, thus promoting the migration and invasion of cervical cancer cells [35]. Chen X et al. found SPP1 inhibition enhanced the chemosensitivity of cervical cancer cell lines to cisplatin [36]. Using microarray analysis, Song JY et al. found that APOD played in the invasion of cervical cancer [37]. These genes play different roles in other tumors. For example, Allina DO et al. demonstrated the diagnostic significance of APOD for prostatic neoplasms [38]. Shen T et al. held that CXCL8 induced EMT in colon cancer [39]. Wang Y et al. found that CXCL8 regulated the development of breast cancer [40]. Ha H et al. found that CXCL8 was also involved in inflammatory diseases in addition to tumors [41]. Liu M et al. argued that MMP1 promoted the growth and metastasis of esophageal squamous cell carcinoma [42]. MMP1 also participated in breast and ovarian cancer [43, 44]. Banik D et al. demonstrated that MMP3 regulated tumor progression [45]. Ji Y et al. demonstrated that C/EBPβ promoted tumor cell invasion and metastasis of colorectal cancer [46]. PLOD2 is implicated in cervical cancer [35] and renal cell carcinoma [47]. SPP1 linked with lung adenocarcinoma, gastric cancer and colorectal cancer [48,49,50].

Sorting nexin 10 (SNX10) can suppress the progression of colorectal cancer, liver cancer and stomach cancer [51, 52]. Cervantes-Anaya N et al. demonstrated SNX10/V-ATPase pathway regulated ciliogenesis in vitro and in vivo [52]. Prostaglandin D synthase (PTGDS) has also been intensely studied [53]. Omori K et al. demonstrated that PTGDS attenuated the malignance of tumor endothelial cells and regulated the processes in non-small cell lung cancer and gastric cancer [54, 55]. The present study is the first to report the expression and prognostic calue of these two genes in CC. Their methylation is associated with CC prognosis, a finding that has never been reported before.

This study has some limitations. First, the analysis is entirely based on open databases, so its results should be validated with functional experiments. Second, the five genetic profiles are based on a single cohort with relatively small sample size. Further studies should involve larger independent cohorts.


Our study indicated that two novel genes PTGDS and SNX10 showed abnormal expression and methylation associated with CC development and explored their prognostic value. However, biofunctions of two genes remained to be unveiled with more in-depth research. The two genes might serve as potential prognostic biomarkers and therapeutic targets in the treatment of CC.

Availability of data and materials

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



Cervical cancer


Differentially expressed gene


Gene Expression Omnibus




Gene Ontology


Kyoto Encyclopedia of Genes and Genomes


Connectivity Map


Protein–protein interaction


Search Tool for the Retrieval of Interacting Genes Database


Molecular Complex Detection


Gene Expression Profiling Interactive Analysis


Receiver operating characteristic


Gene set enrichment analysis


  1. Torre LA, Bray F, Siegel RL, Ferlay J, Lortet-Tieulent J, Jemal A. Global cancer statistics, 2012. CA Cancer J Clin. 2015;65(2):87–108.

    Article  Google Scholar 

  2. Maguire R, Kotronoulas G, Simpson M, Paterson C. A systematic review of the supportive care needs of women living with and beyond cervical cancer. Gynecol Oncol. 2015;136(3):478–90.

    Article  PubMed  Google Scholar 

  3. Petryszak R, Burdett T, Fiorelli B, Fonseca NA, Gonzalez-Porta M, Hastings E, et al. Expression atlas update--a database of gene and transcript expression from microarray- and sequencing-based functional genomics experiments. Nucleic Acids Res. 2014;42(Database issue):D926–32.

    CAS  Article  PubMed  Google Scholar 

  4. Nannini M, Pantaleo MA, Maleddu A, Astolfi A, Formica S, Biasco G. Gene expression profiling in colorectal cancer using microarray technologies: results and perspectives. Cancer Treat Rev. 2009;35(3):201–9.

    CAS  Article  PubMed  Google Scholar 

  5. Zhai Y, Kuick R, Nan B, Ota I, Weiss SJ, Trimble CL, et al. Gene expression analysis of preinvasive and invasive cervical squamous cell carcinomas identifies HOXC10 as a key mediator of invasion. Cancer Res. 2007;67(21):10163–72.

    CAS  Article  PubMed  Google Scholar 

  6. den Boon JA, Pyeon D, Wang SS, Horswill M, Schiffman M, Sherman M, et al. Molecular transitions from papillomavirus infection to cervical precancer and cancer: role of stromal estrogen receptor signaling. Proc Natl Acad Sci U S A. 2015;112(25):E3255–64.

    CAS  Article  Google Scholar 

  7. Pappa KI, Polyzos A, Jacob-Hirsch J, Amariglio N, Vlachos GD, Loutradis D, et al. Profiling of discrete gynecological cancers reveals novel transcriptional modules and common features shared by other cancer types and embryonic stem cells. PLoS One. 2015;10(11):e0142229.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  8. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  9. Piwecka M, Rolle K, Belter A, Barciszewska AM, Zywicki M, Michalak M, et al. Comprehensive analysis of microRNA expression profile in malignant glioma tissues. Mol Oncol. 2015;9(7):1324–40.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  10. Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. Omics. 2012;16(5):284–7.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  11. von Mering C, Huynen M, Jaeggi D, Schmidt S, Bork P, Snel B. STRING: a database of predicted functional associations between proteins. Nucleic Acids Res. 2003;31(1):258–61.

    CAS  Article  Google Scholar 

  12. Bandettini WP, Kellman P, Mancini C, Booker OJ, Vasu S, Leung SW, et al. MultiContrast Delayed Enhancement (MCODE) improves detection of subendocardial myocardial infarction by late gadolinium enhancement cardiovascular magnetic resonance: a clinical validation study. J Cardiovasc Magn Reson. 2012;14(1):83.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Lamb J. The connectivity map: a new tool for biomedical research. Nat Rev Cancer. 2007;7(1):54–60.

    Article  PubMed  Google Scholar 

  14. Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45(W1):W98–w102.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  15. Rhodes DR, Kalyana-Sundaram S, Mahavisno V, Varambally R, Yu J, Briggs BB, et al. Oncomine 3.0: genes, pathways, and networks in a collection of 18,000 cancer gene expression profiles. Neoplasia (New York, NY). 2007;9(2):166–80.

    CAS  Article  Google Scholar 

  16. Subramanian A, Kuehn H, Gould J, Tamayo P, Mesirov JP. GSEA-P: a desktop application for gene set enrichment analysis. Bioinformatics (Oxford, England). 2007;23(23):3251–3.

    CAS  Article  Google Scholar 

  17. Su Z, Yang H, Zhao M, Wang Y, Deng G, Chen R. MicroRNA-92a promotes cell proliferation in cervical cancer via inhibiting p21 expression and promoting cell cycle progression. Oncol Res. 2017;25(1):137–45.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Pera B, Calvo-Vidal MN, Ambati S, Jordi M, Kahn A, Diaz JF, et al. High affinity and covalent-binding microtubule stabilizing agents show activity in chemotherapy-resistant acute myeloid leukemia cells. Cancer Lett. 2015;368(1):97–104.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  19. Huda MN, Erdene-Ochir E, Pan CH. Assay for phosphorylation and microtubule binding along with localization of tau protein in colorectal cancer cells. J Vis Exp. 2017;(128).

  20. Roy D, Sheng GY, Herve S, Carvalho E, Mahanty A, Yuan S, et al. Interplay between cancer cell cycle and metabolism: Challenges, targets and therapeutic opportunities. Biomed Pharmacother. 2017;89:288–96.

    CAS  Article  Google Scholar 

  21. Newell M, Baker K, Postovit LM, Field CJ. A critical review on the effect of docosahexaenoic acid (DHA) on cancer cell cycle progression. Int J Mol Sci. 2017;18(8).

  22. Zhang P, Kawakami H, Liu W, Zeng X, Strebhardt K, Tao K, et al. Targeting CDK1 and MEK/ERK overcomes apoptotic resistance in BRAF-mutant human colorectal cancer. Mol Cancer Res. 2018;16(3):378–89.

    CAS  Article  PubMed  Google Scholar 

  23. Luo Y, Wu Y, Peng Y, Liu X, Bie J, Li S. Systematic analysis to identify a key role of CDK1 in mediating gene interaction networks in cervical cancer development. Ir J Med Sci. 2016;185(1):231–9.

    CAS  Article  PubMed  Google Scholar 

  24. Zeng Y, Stauffer S, Zhou J, Chen X, Chen Y, Dong J. Cyclin-dependent kinase 1 (CDK1)-mediated mitotic phosphorylation of the transcriptional co-repressor Vgll4 inhibits its tumor-suppressing activity. J Biol Chem. 2017;292(36):15028–38.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  25. Bednarek K, Kiwerska K, Szaumkessel M, Bodnar M, Kostrzewska-Poczekaj M, Marszalek A, et al. Recurrent CDK1 overexpression in laryngeal squamous cell carcinoma. Tumour Biol. 2016;37(8):11115–26.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  26. Tseng TH, Chien MH, Lin WL, Wen YC, Chow JM, Chen CK, et al. Inhibition of MDA-MB-231 breast cancer cell proliferation and tumor growth by apigenin through induction of G2/M arrest and histone H3 acetylation-mediated p21(WAF1/CIP1) expression. Environ Toxicol. 2017;32(2):434–44.

    CAS  Article  PubMed  Google Scholar 

  27. Chen L, Jin T, Zhu K, Piao Y, Quan T, Quan C, et al. PI3K/mTOR dual inhibitor BEZ235 and histone deacetylase inhibitor Trichostatin A synergistically exert anti-tumor activity in breast cancer. Oncotarget. 2017;8(7):11937–49.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Munshi PN, Lubin M, Bertino JR. 6-thioguanine: a drug with unrealized potential for cancer therapy. Oncologist. 2014;19(7):760–5.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  29. Zhong G, Chen L, Yin R, Qu Y, Bao Y, Xiao Q, et al. Chemokine (CC motif) ligand 21/CC chemokine receptor type 7 triggers migration and invasion of human lung cancer cells by epithelialmesenchymal transition via the extracellular signalregulated kinase signaling pathway. Mol Med Rep. 2017;15(6):4100–8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  30. Ehling J, Tacke F. Role of chemokine pathways in hepatobiliary cancer. Cancer Lett. 2016;379(2):173–83.

    CAS  Article  PubMed  Google Scholar 

  31. Xiao G, Wang X, Wang J, Zu L, Cheng G, Hao M, et al. CXCL16/CXCR6 chemokine signaling mediates breast cancer progression by pERK1/2-dependent mechanisms. Oncotarget. 2015;6(16):14165–78.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Yan R, Shuai H, Luo X, Wang X, Guan B. The clinical and prognostic value of CXCL8 in cervical carcinoma patients: immunohistochemical analysis. Biosci Rep. 2017;37(5).

  33. Tian R, Li X, Gao Y, Li Y, Yang P, Wang K. Identification and validation of the role of matrix metalloproteinase-1 in cervical cancer. Int J Oncol. 2018;52(4):1198–208.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  34. Xie B, Zhang Z, Wang H, Chen Z, Wang Y, Liang H, et al. Genetic polymorphisms in MMP 2, 3, 7, and 9 genes and the susceptibility and clinical outcome of cervical cancer in a Chinese Han population. Tumour Biol. 2016;37(4):4883–8.

    CAS  Article  PubMed  Google Scholar 

  35. Xu F, Zhang J, Hu G, Liu L, Liang W. Hypoxia and TGF-beta1 induced PLOD2 expression improve the migration and invasion of cervical cancer cells by promoting epithelial-to-mesenchymal transition (EMT) and focal adhesion formation. Cancer Cell Int. 2017;17(1):54.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  36. Chen X, Xiong D, Ye L, Yang H, Mei S, Wu J, et al. SPP1 inhibition improves the cisplatin chemo-sensitivity of cervical cancer cell lines. Cancer Chemother Pharmacol. 2019;83(4):603–13.

    CAS  Article  PubMed  Google Scholar 

  37. Song JY, Lee JK, Lee NW, Jung HH, Kim SH, Lee KW. Microarray analysis of normal cervix, carcinoma in situ, and invasive cervical cancer: identification of candidate genes in pathogenesis of invasion in cervical cancer. Int J Gynecol. 2008;18(5):1051–9.

    CAS  Article  Google Scholar 

  38. Allina DO, Andreeva YY, Zavalishina LE, Moskvina LV, Frank GA. Estimation of the diagnostic potential of APOD, PTOV1, and EPHA4 for prostatic neoplasms. Arkh Patol. 2016;78(5):9–14.

    CAS  Article  PubMed  Google Scholar 

  39. Shen T, Yang Z, Cheng X, Xiao Y, Yu K, Cai X, et al. CXCL8 induces epithelial-mesenchymal transition in colon cancer cells via the PI3K/Akt/NF-kappaB signaling pathway. Oncol Rep. 2017;37(4):2095–100.

    CAS  Article  PubMed  Google Scholar 

  40. Wang Y, Liu J, Jiang Q, Deng J, Xu F, Chen X, et al. Human adipose-derived mesenchymal stem cell-secreted CXCL1 and CXCL8 facilitate breast tumor growth by promoting angiogenesis. Stem Cells (Dayton, Ohio). 2017;35(9):2060–70.

    CAS  Article  Google Scholar 

  41. Ha H, Debnath B, Neamati N. Role of the CXCL8-CXCR1/2 axis in cancer and inflammatory diseases. Theranostics. 2017;7(6):1543–88.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  42. Liu M, Hu Y, Zhang MF, Luo KJ, Xie XY, Wen J, et al. MMP1 promotes tumor growth and metastasis in esophageal squamous cell carcinoma. Cancer Lett. 2016;377(1):97–104.

    CAS  Article  PubMed  Google Scholar 

  43. Yokoi A, Yoshioka Y, Yamamoto Y, Ishikawa M, Ikeda SI, Kato T, et al. Malignant extracellular vesicles carrying MMP1 mRNA facilitate peritoneal dissemination in ovarian cancer. Nat Commun. 2017;8(1):14470.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Shen CJ, Kuo YL, Chen CC, Chen MJ, Cheng YM. MMP1 expression is activated by slug and enhances multi-drug resistance (MDR) in breast cancer. PLoS One. 2017;12(3):e0174487.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  45. Banik D, Netherby CS, Bogner PN, Abrams SI. MMP3-mediated tumor progression is controlled transcriptionally by a novel IRF8-MMP3 interaction. Oncotarget. 2015;6(17):15164–79.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Ji Y, Li J, Li P, Wang L, Yang H, Jiang G. C/EBPbeta promotion of MMP3-dependent tumor cell invasion and association with metastasis in colorectal cancer. Genet Test Mol Biomarkers. 2018;22(1):5–10.

    CAS  Article  PubMed  Google Scholar 

  47. Kurozumi A, Kato M, Goto Y, Matsushita R, Nishikawa R, Okato A, et al. Regulation of the collagen cross-linking enzymes LOXL2 and PLOD2 by tumor-suppressive microRNA-26a/b in renal cell carcinoma. Int J Oncol. 2016;48(5):1837–46.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  48. Zhang Y, Du W, Chen Z, Xiang C. Upregulation of PD-L1 by SPP1 mediates macrophage polarization and facilitates immune escape in lung adenocarcinoma. Exp Cell Res. 2017;359(2):449–57.

    CAS  Article  PubMed  Google Scholar 

  49. Xu C, Sun L, Jiang C, Zhou H, Gu L, Liu Y, et al. SPP1, analyzed by bioinformatics methods, promotes the metastasis in colorectal cancer by activating EMT pathway. Biomed Pharmacother. 2017;91:1167–77.

    CAS  Article  Google Scholar 

  50. Chen LZ, He CY, Su X, Peng JL, Chen DL, Ye Z, et al. SPP1 rs4754 and its epistatic interactions with SPARC polymorphisms in gastric cancer susceptibility. Gene. 2018;640:43–50.

    CAS  Article  PubMed  Google Scholar 

  51. Zhang J, Wu Y, Jin HY, Guo S, Dong Z, Zheng ZC, et al. Prognostic value of sorting nexin 10 weak expression in stomach adenocarcinoma revealed by weighted gene co-expression network analysis. World J Gastroenterol. 2018;24(43):4906–19.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  52. Cervantes-Anaya N, Ponciano-Gomez A, Lopez-Alvarez GS, Gonzalez-Reyes C, Hernandez-Garcia S, Cabanas-Cortes MA, et al. Downregulation of sorting nexin 10 is associated with overexpression of miR-30d during liver cancer progression in rats. Tumour Biol. 2017;39(4):1010428317695932.

    CAS  Article  PubMed  Google Scholar 

  53. Munkholm K, Peijs L, Kessing LV, Vinberg M. Reduced mRNA expression of PTGDS in peripheral blood mononuclear cells of rapid-cycling bipolar disorder patients compared with healthy control subjects. Int J Neuropsychopharmacol. 2014;18(5).

  54. Zhang B, Bie Q, Wu P, Zhang J, You B, Shi H, et al. PGD2/PTGDR2 signaling restricts the self-renewal and tumorigenesis of gastric cancer. Stem Cells (Dayton, Ohio). 2018;36(7):990–1003.

    CAS  Article  Google Scholar 

  55. He LP, Chen YF, Yang J. Investigation on the role and mechanism of prostagland in D2 synthase in non-small cell lung cancer. Zhonghua Yi Xue Za Zhi. 2017;97(38):3022–7.

    CAS  PubMed  Google Scholar 

  56. Pyeon D, Newton MA, Lambert PF, den Boon JA, Sengupta S, Marsit CJ, Woodworth CD, Connor JP, Haugen TH, Elaine M. Smith, Karl T. Kelsey, Lubomir P. Turek, Paul Ahlquist. Fundamental Differences in Cell Cycle Deregulation in Human Papillomavirus–Positive and Human Papillomavirus–Negative Head/Neck and Cervical Cancers. Can Res 2007;67 (10):4605–19.

  57. Scotto L, Narayan G, Nandula SV, Arias-Pulido H, Subramaniyam S, Schneider A, Kaufmann AM, Wright JD, Pothuri B, Mansukhani M, Murty VV. Identification of copy number gain and overexpressed genes on chromosome arm 20q by an integrative genomic approach in cervical cancer: potential role in progression. Genes Chromosomes Cancer. 2008;47(9):755–65.

  58. Espinosa AM, Alfaro A, Roman-Basaure E, Guardado-Estrada M, Palma Í, Serralde C, Medina I, Juárez E, Bermúdez M, Márquez E, Borges-Ibáñez M, Muñoz-Cortez S, Alcántara-Vázquez A, Alonso P, Curiel-Valdez J, Kofman S, Villegas N, Berumen J, Scheurer M. Mitosis Is a Source of Potential Markers for Screening and Survival and Therapeutic Targets in Cervical Cancer. PLoS ONE. 2013;8(2):e55975.

  59. den Boon JA, Pyeon D, Wang SS, Horswill M, Schiffman M, Sherman M, Zuna RE, Wang Z, Hewitt SM, Pearson R, Schott M, Chung L, He Q, Lambert P, Walker J,  Newton ML, Wentzensen N, Ahlquist P. Molecular transitions from papillomavirus infection to cervical precancer and cancer: Role of stromal estrogen receptor signaling. Proc Natl Acad Sci. 2015;112 (25):E3255-E3264.

Download references


Not applicable.


This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

Author information

Authors and Affiliations



Authors SLF, MXX and PPJ designed the project. Authors YC and FG contributed on data analysis, PPJ, JHL and WS prepared the main manuscript. ZYM made revisions in language and finalization of the manuscript. All authors have reviewed and approved the manuscript.

Corresponding authors

Correspondence to Manxin Xie or Shilong Fu.

Ethics declarations

Ethics approval and consent to participate

This study does not involve the use of any animal or human data or tissue. The approval or the consent were not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1

. Standardization of gene expression. The blue bar represents the data before normalization, and the red bar represents the normalized data.(A) The standardization of GSE6791 data, (B) the standardization of GSE9750 data, (C) the standardization of GSE39001 data (D) the standardization of GSE63514 data.

Additional file 2: Figure S2

. Volcanic maps of DEGs from all samples. The red dots represent upregulated genes screened according to |fold change| ≥ 2.0 and a corrected P-value ≤0.05. The green points represent downregulated genes and the cutoff value is the same as that of upregulated genes. The black points represent genes with no significant difference. (A)GSE6791 data, (B) GSE9750 data, (C) GSE39001 data and (D) GSE63514 data.

Additional file 3: Figure S3

. Heatmap of DEGs from 4 datasets. Red cubes represent upregulated genes, green cubes represent downregulated genes, and black represent genes with no significance. Genes underexpressed are painted with gray. DEGs are screened by criterion:|fold change| ≥ 2.0 and a corrected P-value ≤0.05. (A) GSE6791 data, (B) GSE9750 data, (C) GSE39001 data and (D) GSE63514 data.

Additional file 4: Figure S4

. GO enrichment analysis of DEGs in cervical cancer. (A) GO analysis divided DEGs into three functional groups: molecular function, biological processes, and cell composition. (B) The top 10 GO terms of DEGs in CC. The outer circle shows a scatter plot for each term of the logFC of the assigned genes. Red circles display upregulation and blue ones downregulation.(C) Distribution of DEGs in cervical cancer for different GO-enriched functions.

Additional file 5: Figure S5.

(A) Collections between significant GO groups and DEGs. Blue rounds represent the GO groups, green rounds represent downregulated genes and red rounds represent upregulated genes. (B) Collections between significant KEGG groups and DEGs. Blue rounds represent signaling pathway, green rounds represent downregulated genes and red rounds represent upregulated genes.

Additional file 6: Figure S6

. (A) A significant module selected from protein–protein interaction network. (B) GO analysis of these significant molecule.(C) KEGG analysis of these significant molecule.

Additional file 7: Figure S7.

Validation of genes expression in GEPIA. (A) APOD, (B) CXCL8, (C) MMP1, (D) MMP3, (E) PLOD2, (F) PTGDS, (G) SNX10, (H) SPP1.

Additional file 8: Figure S8.

Survival analysis of genes in GEPIA. (A) APOD, (B) CXCL8, (C) MMP1, (D) MMP3, (E) PLOD2, (F) PTGDS, (G) SNX10, (H) SPP1.

Additional file 9: Figure S9.

Gene expression in GSE7803 and GSE29570. (A) PTGDS in GSE7803. (B) SNX10 in GSE7803. (C) PTGDS in GSE29570. (D) SNX10 in GSE29570.

Additional file 10: Figure S10.

Seven genes were screened by Univariate cox regression analysis.

Additional file 11: Figure S11.

Five genes were screened by Multivariate cox regression analysis.

Additional file 12: Figure S12.

Regression analysis of the 5 genes identified. (A) The expression level of the 5 genes in low- and high-risk groups. (B) The univariate Cox proportional hazards regression. (C) The multivariate Cox proportional hazards regression. (D) The heatmap of the 5 genes in high- and low-risk patients in TCGA dataset.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Jiang, P., Cao, Y., Gao, F. et al. SNX10 and PTGDS are associated with the progression and prognosis of cervical squamous cell carcinoma. BMC Cancer 21, 694 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Cervical squamous cell carcinoma
  • Methylation
  • Prognosis
  • Risk model
  • SNX10