BackgroundConsiderable evidence suggests that the heterogeneity of ovarian cancer (OC) is a major cause of treatment failure. Single-cell RNA sequencing (scRNA-seq) is a powerful tool to analyse the heterogeneity of the tumour at the single-cell level, leading to a better understanding of cell function at the genetic and cellular levels.MethodsOC scRNA-seq data were extracted from the Gene Expression Omnibus (GEO) database and the FindCluster () package used for cell cluster analysis. The GSVA package was used for single-sample gene set enrichment analysis (ssGSEA) analysis to obtain a Hallmark gene set score and bulk RNA-seq data were used to analyse the key genes of OC-associated immune cell subsets. CIBERSORT was used to identify immune scores of cells and the “WGCNA” package for the weighted correlation network analysis (WGCNA). KEGG (Kyoto Encyclopedia of Genes and Genomes) and GO (Gene Ontology) analyses of subtype groups were performed by GSEA. Then, univariate Cox and lasso regression were performed to further establish a signature. Finally, qPCR and immunohistochemistry staining were used to evaluate the expression of signature genes in OC.ResultsTwo scRNA-seq (GSE154600 and GES158937) datasets were integrated to obtain 20 cell clusters. T cells or NK cells (cluster 5, 6, 7, 11), B cells (cluster 16, 19, 20) and myeloid cells (cluster 4, 9, 10) were clustered according to immune cell markers. The ssGSEA revealed that M1- and M2-like myeloid cell-related genes were significantly upregulated in P3 and P4 patients in the GSE154600 data. Immune cell analysis in TCGA-OC showed that a high abundance of M1-like tumour-associated macrophages (TAMS) predicts better survival. WGCNA, univariate Cox and lasso Cox regression established a two-gene signature (RiskScore=-0.059*CXCL13-0.034*IL26). Next, the TCGA-test and TCGA-OC were used to test the risk prediction ability of the signature, showing a good effect in the datasets. Moreover, the qPCR and immunohistochemistry staining revealed that the expression of CXCL13 and IL26 was reduced in OC tissues.ConclusionA two-gene signature prognostic stratification system (CXCL13 and IL26) was developed based on the heterogeneity of OC immune cells to accurately evaluate the prognostic risk.

Background Ovarian cancer (OC) is the most lethal gynaecological tumor. Changes in glycolysis have been proven to play an important role in OC progression. We aimed to identify a novel glycolysis-related gene signature to better predict the prognosis of patients with OC. Methods mRNA and clinical data were obtained from The Cancer Genome Atlas (TCGA), International Cancer Genome Consortium (ICGC) and Genotype Tissue Expression (GTEx) database. The “limma” R package was used to identify glycolysis-related differentially expressed genes (DEGs). Then, a multivariate Cox proportional regression model and survival analysis were used to develop a glycolysis-related gene signature. Furthermore, the TCGA training set was divided into two internal test sets for validation, while the ICGC dataset was used as an external test set. A nomogram was constructed in the training set, and the relative proportions of 22 types of tumor-infiltrating immune cells were evaluated using the “CIBERSORT” R package. The enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were determined by single-sample gene set enrichment analysis (ssGSEA) with the “GSVA” R package. Finally, the expression and function of the unreported signature genes ISG20 and SEH1L were explored using immunohistochemistry, western blotting, qRT-PCR, proliferation, migration, invasion and xenograft tumor assays. Results A five-gene signature comprising ANGPTL4, PYGB, ISG20, SEH1L and IRS2 was constructed. This signature could predict prognosis independent of clinical factors. A nomogram incorporating the signature and three clinical features was constructed, and the calibration plot suggested that the nomogram could accurately predict the survival rate. According to ssGSEA, the signature was associated with KEGG pathways related to axon guidance, mTOR signalling, tight junctions, etc. The proportions of tumor-infiltrating immune cells differed significantly between the high-risk group and the low-risk group. The expression levels of ISG20 and SEH1L were lower in tumor tissues than in normal tissues. Overexpression of ISG20 or SEH1L suppressed the proliferation, migration and invasion of Caov3 cells in vitro and the growth of xenograft tumors in vivo. Conclusion Five glycolysis-related genes were identified and incorporated into a novel risk signature that can effectively assess the prognosis and guide the treatment of OC patients.

Background Ovarian cancer (OC) is an invasive gynaecologic cancer with a high cancer-related death rate. The purpose of this study was to establish an invasion-related multigene signature to predict the prognostic risk of OC. Methods We extracted 97 invasion-related genes from The Cancer Genome Atlas (TCGA) database. Then, the ConsensusClusterPlus and limma packages were used to calculate differentially expressed genes (DEGs). To calculate the immune scores of the molecular subtypes, we used ESTIMATE to evaluate the stromal score, immune score and ESTIMATE score. MCP-counter and the GSVA package ssgsea were used to evaluate the types of infiltrating immune cells. Survival and nomogram analyses were performed to explore the prognostic value of the signature. Finally, qPCR, immunohistochemistry staining and functional assays were used to evaluate the expression and biological abilities of the signature genes in OC. Results Based on the consistent clustering of invasion-related genes, cases in the OC datasets were divided into two subtypes. A significant difference was observed in prognosis between the two subtypes. Most genes were highly expressed in the C1 group. Based on the C1 group genes, we constructed an invasion-related 6-gene prognostic risk model. Furthermore, to verify the signature, we used the TCGA-test and GSE32062 and GSE17260 chip datasets for testing and finally obtained a good risk prediction effect in those datasets. Moreover, the results of the qPCR and immunohistochemistry staining assays revealed that KIF26B, VSIG4 and COL6A6 were upregulated and that FOXJ1, MXRA5 and CXCL9 were downregulated in OC tissues. The functional study showed that the expression of KIF26B, VSIG4, COL6A6, FOXJ1, MXRA5 and CXCL9 can regulate the migration and invasion abilities of OC cells. Conclusion We developed a 6-gene prognostic stratification system (FOXJ1, MXRA5, KIF26B, VSIG4, CXCL9 and COL6A6) that is independent of clinical features. These results suggest that the signature could potentially be used to evaluate the prognostic risk of OC patients.

PurposeTo explore the relationship between retinoic acid receptor gamma (RARG) and ovarian cancer (OC) cell proliferation and the prognosis of patients.MethodsThe transcriptome and clinical information of 379 OC and 88 normal ovarian samples were downloaded from the Cancer Genome Atlas (TCGA) database and the Genotype Tissue Expression (GTEx) database. We compared the mRNA level of RARG between ovrian normal and tumor tissues with the Wilcoxon rank sum test.The R package “limma” was used to analyze the differences in RARG expression between different clinical subgroups. Kaplan−Meier analysis was applied to evaluate the correlation between RARG and prognosis of patients. A nomogram was established to predict the effect of RARG on prognosis of OC patients. Immunohistochemistry and qRT−PCR experiments were conducted to determine the differential expression of RARG between ovarian normal and tumor tissues. Finally, we altered RARG expression using specific siRNA and lentiviral expression vectors to explore the function of RARG by CCK-8, cell cycle, colony formation, and xenograft assays in nude mice.ResultsRARG was highly expressed in ovarian tumors and was an independent predictor of poor overall survival outcomes. Subgroup analysis showed the high expression of RARG was related to FIGO stage III-IV (P=0.027), overall survival time <5 years (P=0.013) and dead status (P=0.041). The Kaplan-Meier curve indicated that patients with high RARG expression level had poor prognosis. The area under the curve (AUC) of RAGR expression for predicting patient survival rates at 1, 5 and 9 years were 0.659, 0.616 and 0.627, respectively. The GSEA enrichment analysis revealed that RARG was involved in ovarian cancer progression through multiple pathways. In cellular experiments in vitro, downregulation of RARG expression significantly suppressed the proliferation and colony formation capacity of OC cells. In cellular experiments in vivo, knockdown of RARG significantly reduced tumor growth in nude mice, decreased expression levels of Ki-67 and proliferation cell nuclear antigen (PCNA).ConclusionsHigh expression of RARG could promote OC cell proliferation and was an independent predictor of poor prognosis. RARG might work as a potential molecular target and biomarker for individualized diagnosis and treatment in OC patients.

Background: 5-Methylcytidine (m5C) is the most common RNA modification and plays an important role in multiple tumors including cervical cancer (CC). We aimed to develop a novel gene signature by identifying m5C modification subtypes of CC to better predict the prognosis of patients.Methods: We obtained the expression of 13 m5C regulatory factors from The Cancer Genome Atlas (TCGA all set, 257 patients) to determine m5C modification subtypes by the “nonnegative matrix factorization” (NMF). Then the “limma” package was used to identify differentially expressed genes (DEGs) between different subtypes. According to these DEGs, we performed Cox regression and Kaplan-Meier (KM) survival analysis to establish a novel gene signature in TCGA training set (128 patients). We also verified the risk prediction effect of gene signature in TCGA test set (129 patients), TCGA all set (257 patients) and GSE44001 (300 patients). Furthermore, a nomogram including this gene signature and clinicopathological parameters was established to predict the individual survival rate. Finally, the expression and function of these signature genes were explored by qRT-PCR, immunohistochemistry (IHC) and proliferation, colony formation, migration and invasion assays.Results: Based on consistent clustering of 13 m5C-modified genes, CC was divided into two subtypes (C1 and C2) and the C1 subtype had a worse prognosis. The 4-gene signature comprising FNDC3A, VEGFA, OPN3 and CPE was constructed. In TCGA training set and three validation sets, we found the prognosis of patients in the low-risk group was much better than that in the high-risk group. A nomogram incorporating the gene signature and T stage was constructed, and the calibration plot suggested that it could accurately predict the survival rate. The expression levels of FNDC3A, VEGFA, OPN3 and CPE were all high in cervical cancer tissues. Downregulation of FNDC3A, VEGFA or CPE expression suppressed the proliferation, migration and invasion of SiHa cells.Conclusions: Two m5C modification subtypes of CC were identified and then a 4-gene signature was established, which provide new feasible methods for clinical risk assessment and targeted therapies for CC.