Discussion
In summary, our findings suggested novel molecular subtypes of cervical cancer with distinct immunometabolic profiles and uncovered a novel therapeutic target.
Methods
Bulk sequencing and single-cell sequencing data were downloaded from the Cancer Genome Atlas (TCGA) database and Gene Expression Omnibus (GEO) database respectively. Immune and metabolism related genes were collected from Immport and Kyoto encyclopedia of genes and genomes (KEGG) database respectively. Unsupervised consensus clustering was performed to identify the molecular subtypes. Cibersort was applied to evaluate the immune cells infiltration status. Differential expression analysis and Gene set enrichment analysis (GSEA) were performed to characterize the molecular pattern of different subtypes. Multivariate Cox regression analysis was used for prognosis prediction model construction and receiver operating characteristic (ROC) curve was used for performance evaluation. The hub genes in the model were verified in single-cell sequencing dataset and clinical specimens. In vitro experiments were performed to validate the findings in our research.
Results
Three subtypes were identified with prognostic implications. C1 subgroup was in an immunosuppressive state with activation of mitochondrial cytochrome P450 metabolism, C2 had poor immune cells infiltration and was characterized by tRNA anabolism, and the C3 subgroup was in an inflammatory state with activation of aromatic amino acid synthesis. The area under the ROC curve of the constructed model was 0.8, which showed better performance than clinical features. IMPDH1 was found to be significantly upregulated in tumor tissue and it was demonstrated that IMPDH1 could be a novel therapeutic target in vitro.