Identification of 10 differentially expressed genes involved in the tumorigenesis of cervical cancer via next-generation sequencing

Background: The incidence and mortality of cervical cancer remain high in female malignant tumors worldwide. There is still a lack of diagnostic and prognostic markers for cervical carcinoma. This study aimed to screen differentially expressed genes (DEGs) between normal and cervical cancer tissues to identify candidate genes for further research. Methods: Uterine cervical specimens were resected from our clinical patients after radical hysterectomy. Three patients' transcriptomic datasets were built by the next generation sequencing (NGS) results. DEGs were selected through the edgeR and DESeq2 packages in the R environment. Functional enrichment analysis, including GO/DisGeNET/KEGG/Reactome enrichment analysis, was performed. Normal and cervical cancer tissue data from the public databases TCGA and GTEx were collected to compare the expression levels of 10 selected DEGs in tumor and normal tissues. ROC curve and survival analysis were performed to compare the diagnostic and prognostic values of each gene. The expression levels of candidate genes were verified in 15 paired clinical specimens via quantitative real-time polymerase chain reaction. Results: There were 875 up-regulated and 1,482 down-regulated genes in cervical cancer samples compared with the paired adjacent normal cervical tissues according to the NGS analysis. The top 10 DEGs included APOD, MASP1, ACKR1, C1QTNF7, SFRP4, HSPB6, GSTM5, IGFBP6, F10 and DCN. GO, DisGeNET and Reactome analyses revealed that the DEGs were related to extracellular matrix and angiogenesis which might influence tumorigenesis. KEGG enrichment showed that PI3K-Akt signaling pathway might be involved in cervical cancer tumorigenesis and progression. The expression levels of selected genes were decreased in tumors in both the public database and our experimental clinical specimens. All the candidate genes showed excellent diagnostic value, and the AUC values exceeded 0.90. Additionally, APOD, ACKR1 and SFRP4 expression levels could help predict the prognosis of patients with cervical cancer. Conclusions: In this study, we selected the top 10 DEGs which were down-regulated in cervical cancer tissues. All of them had dramatically diagnostic value. APOD, ACKR1 and SFRP4 were associated with the survivals of cervical cancer. C1QTNF7, HSPB6, GSTM5, IGFBP6 and F10 were first reported to be candidate genes of cervical carcinoma.