Angiogenesis-related gene signatures reveal the prognosis of cervical cancer based on single cell sequencing and co-expression network analysis.

Abstract

Cervical cancer ranks first in female reproductive tract tumors in terms of morbidity and mortality. Yet the curative effect of patients with persistent, recurrent or metastatic cervical cancer remains unsatisfactory. Although antitumor angiogenic drugs have been recommended as the first-line treatment options for cervical cancer, there are no comprehensive prognostic indicators for cervical cancer based on angiogenic signature genes. In this study, we aimed to develop a model to assess the prognosis of cervical cancer based on angiogenesis-related (AG) signature genes, and to provide some reference for the comprehensive treatment of cervical cancer in the clinical setting. First we screened the AG gene set from GeneCard website, and then performed angiogenesis-related scores (AGS) per cell from single cell sequencing dataset GSE168652, followed by performing weighted gene co-expression network analysis (WGCNA) for cervical cancer patients according to angiogenesis phenotype. Thus, we established a prognostic model based on AGS by taking the intersection of WGCNA angiogenic module gene and differential gene (DEGs) of GSE168652. The GSE44001 was selected as an external validation set, followed by performing ROC curve analysis to assess its accuracy. The results showed that we successfully constructed a prognostic model related to the AG genes. Patients in the high-AGS group in both the train, test and the validation sets had a worse prognosis than those in the low-AGS group, had lower expression of most immune checkpoint-associated genes and lower tumor mutational burden as well. Patients in the low-AGS group were more sensitive to AMG.706, Bosutinib, and Lenalidomide while Imatinib, Pazopanib, and Sorafenib were more recommended to patients in the high-AGS group. Finally, TXNDC12 and ZC3H13, which have high hazard ratio and poor prognosis in the model, were highly expressed in cervical cancer cell lines and tissue. Meanwhile, the results showed that TXNDC12 promoted the migration of cervical cancer cells and the tubule-forming ability of endothelial cells. In conclusion, our model based on genes with AG features can effectively assess the prognosis of cervical cancer, and can also provide reference for clinicians to choose immune-related treatments.