Abstract:Background: Cervical cancer is the fourth common cancer among women. Its prognosis needs our more attention. Our purpose was to identity new prognostic gene sets to help other researchers develop more effective treatment for cervical cancer patients and improve the prognosis of patients. Methods: We used gene set variation analysis (GSVA) to calculate the enrichment scores of gene sets and identified three subtypes of cervical cancer through the Cox regression model, k-means clustering algorithm, and nonnegative matrix factorization method (NMF). Chi-square test was utilized to test whether a certain clinical characteristic is different among divided subtypes. We further screened the prognostic gene sets using differential analysis, univariate Cox regression analysis, and least absolute shrinkage and selection operator (LASSO) regression. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were used to analyze which pathways and function the genes from screened gene sets enriched. Search Tool for the Retrieval of Interacting Genes (STRING) was used to draw the protein-protein interaction network, and Cytoscape was used to visualize the hub genes of protein-protein interaction network. Results: We identified three novel subtypes of cervical cancer in The Cancer Genome Atlas (TCGA) samples and validated in Gene Expression Omnibus (GEO) samples. There were significant variations between the three subtypes in histological type, T stage, M stage, and N stage. T_GSE36888_UNTREATED_VS_IL2_TREATED_STAT5_AB_KNOCKIN_TCELL_2H_UP and N_HALLMARK_ANGIOGENESIS were screened prognostic gene sets. The prognostic model was as follows: riskScore = T_GSE36888_UNTREATED_VS_IL2_TREATED_STAT5_AB_KNOCKIN_TCELL_2H_UP∗ 2.617 + N_HALLMARK_ANGIOGENESIS∗ 4.860. Survival analysis presented that in these two gene sets, high enrichment scores were all significantly related to worse overall survival. The hub genes from T gene set included CXCL1, CXCL2, CXCL8, ALDOA, TALDO1, LDHA, CCL4, FCAR, FCER1G, SAMSN1, LILRB1, SH3PXD2B, PPM1N, PKM, and FKBP4. As for N gene sets, the hub genes included ITGAV, PTK2, SPP1, THBD, and APOH. Conclusions: Three novel subtypes and two prognostic gene sets were identified. 15 hub genes for T gene set and 5 hub genes for N gene set were discovered. Based on these findings, we can develop more and more effective treatments for cervical cancer patients. Based on the gene enriched pathways, we can development specific drugs targeting the pathways.

GS-TCGA: Gene Set-Based Analysis of The Cancer Genome Atlas

Identification of Gene Expression Pattern Related to Breast Cancer Survival Using Integrated TCGA Datasets and Genomic Tools.

TCGA4U: A Web-Based Genomic Analysis Platform to Explore and Mine TCGA Genomic Data for Translational Research

Gene Set linkage analysis: a tool for interpreting the overall functional impacts of observed transcriptomic changes

Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles

GSCA: an integrated platform for gene set cancer analysis at genomic, pharmacogenomic and immunogenomic levels

DrGaP: a powerful tool for identifying driver genes and pathways in cancer sequencing studies.

UALCAN: A Portal for Facilitating Tumor Subgroup Gene Expression and Survival Analyses

The Cancer Genome Atlas Clinical Explorer: a web and mobile interface for identifying clinical–genomic driver associations

TCGAnalyzeR: An Online Pan-Cancer Tool for Integrative Visualization of Molecular and Clinical Data of Cancer Patients for Cohort and Associated Gene Discovery

Inferring Pathway Crosstalk Networks Using Gene Set Co-Expression Signatures.

GSVA: gene set variation analysis for microarray and RNA-Seq data

Identification of Gene Signatures for COAD Using Feature Selection and Bayesian Network Approaches

An Integrated Pipeline For Tcga Data Analysis

Abstract 5275: an Integrated Pipeline for TCGA Data Analysis

Assessment of Gene Set Enrichment Analysis using curated RNA-seq-based benchmarks

CGPA: multi-context insights from the cancer gene prognosis atlas

Abstract SY12-03: Multi-tumor Analysis of TCGA Data Identifies Expression Commonalities Across Tumor Types.

IGSA: Individual Gene Sets Analysis, Including Enrichment and Clustering.

A Common Gene Expression Signature Analysis Method for Multiple Types of Cancer

Identification of Specific Cervical Cancer Subtypes and Prognostic Gene Sets in Tumor and Nontumor Tissues Based on GSVA Analysis