A Novel Prognostic Risk Model for Cervical Cancer Based on Immune Checkpoint HLA-G-Driven Differentially Expressed Genes.

Human leukocyte antigen G (HLA-G) is a potential checkpoint molecule that plays a key role in cervical carcinogenesis. The purpose of this study was to construct and validate a prognostic risk model to predict the overall survival (OS) of cervical cancer patients, providing a reference for individualized clinical treatment that may lead to better clinical outcomes. HLA-G-driven differentially expressed genes (DEGs) were obtained from two cervical carcinoma cell lines, namely, SiHa and HeLa, with stable overexpression of HLA-G by RNA sequencing (RNA-seq). The biological functions of these HLA-G-driven DEGs were analysed by GO enrichment and KEGG pathway using the "clusterProfiler" package. The protein-protein interactions (PPIs) were assessed using the STRING database. The prognostic relevance of each DEG was evaluated by univariate Cox regression using the TCGA-CESC dataset. After the TCGA-CESC cohort was randomly divided into training set and testing set, and a prognostic risk model was constructed by LASSO and stepwise multivariate Cox regression analysis in training set and validated in testing set or in different types of cervical cancer set. The predictive ability of the prognostic risk model or nomogram was evaluated by a series of bioinformatics methods. A total of 1108 candidate HLA-G-driven DEGs, including 391 upregulated and 717 downregulated genes, were obtained and were enriched mostly in the ErbB pathway, steroid biosynthesis, and MAPK pathway. Then, an HLA-G-driven DEG signature consisting of the eight most important prognostic genes CD46, LGALS9, PGM1, SPRY4, CACNB3, PLIN2, MSMO1, and DAGLB was identified as a key predictor of cervical cancer. Multivariate Cox regression analysis showed that this signature is an independent risk factor for the overall survival of CESC patients. Kaplan-Meier survival analysis showed that the 5-year overall survival rate is 23.0% and 84.6% for the high-risk and low-risk patients, respectively (P<0.001). The receiver operating characteristic (ROC) curve of this prognostic model with an area under the curve (AUC) was 0.896 for 5 years, which was better than that of other clinical traits. This prognostic risk model was also successfully validated in different subtypes of cervical cancer, including the keratinizing squamous cell carcinoma, non-keratinizing squamous cell carcinoma, squamous cell neoplasms, non-squamous cell neoplasms set. Single-sample gene set enrichment (ssGSEA) algorithm and Tumor Immune Dysfunction and Exclusion (TIDE) analysis confirmed that this signature influence tumour microenvironment and immune checkpoint blockade. A nomogram that integrated risk score, age, clinical stage, histological grade, and pathological type was then built to predict the overall survival of CESC patients and evaluated by calibration curves, AUC, concordance index (C-index) and decision curve analysis (DCA). To summarize, we developed and validated a novel prognostic risk model for cervical cancer based on HLA-G-driven DEGs, and the prognostic signature showed great ability in predicting the overall survival of patients with cervical cancer.

Xu HH ，Wang HL ，Xing TJ ，Wang XQ ... -  《Frontiers in Immunology》

Construction of a novel mRNA-signature prediction model for prognosis of bladder cancer based on a statistical analysis.

Bladder cancer (BC) is a common malignancy neoplasm diagnosed in advanced stages in most cases. It is crucial to screen ideal biomarkers and construct a more accurate prognostic model than conventional clinical parameters. The aim of this research was to develop and validate an mRNA-based signature for predicting the prognosis of patients with bladder cancer. The RNA-seq data was downloaded from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO). Differentially expressed genes (DEGs) were screened in three datasets, and prognostic genes were identified from the training set of TCGA dataset. The common genes between DEGs and prognostic genes were narrowed down to six genes via Least Absolute Shrinkage and Selection Operator (LASSO) regression, and stepwise multivariate Cox regression. Then the gene-based risk score was calculated via Cox coefficient. Time-dependent receiver operating characteristic (ROC) and Kaplan-Meier (KM) survival analysis were used to assess the prognostic power of risk score. Multivariate Cox regression analysis was applied to construct a nomogram. Decision curve analysis (DCA), calibration curves, and time-dependent ROC were performed to assess the nomogram. Finally, functional enrichment of candidate genes was conducted to explore the potential biological pathways of candidate genes. SORBS2, GPC2, SETBP1, FGF11, APOL1, and H1-2 were screened to be correlated with the prognosis of BC patients. A nomogram was constructed based on the risk score, pathological stage, and age. Then, the calibration plots for the 1-, 3-, 5-year OS were predicted well in entire TCGA-BLCA patients. Decision curve analysis (DCA) indicated that the clinical value of the nomogram was higher than the stage model and TNM model in predicting overall survival analysis. The time-dependent ROC curves indicated that the nomogram had higher predictive accuracy than the stage model and risk score model. The AUC of nomogram time-dependent ROC was 0.763, 0.805, and 0.806 for 1-year, 3-year, and 5-year, respectively. Functional enrichment analysis of candidate genes suggested several pathways and mechanisms related to cancer. In this research, we developed an mRNA-based signature that incorporated clinical prognostic parameters to predict BC patient prognosis well, which may provide a novel prognosis assessment tool for clinical practice and explore several potential novel biomarkers related to the prognosis of patients with BC.

Li J ，Cao J ，Li P ，Yao Z ，Deng R ，Ying L ，Tian J ... -  《BMC CANCER》

Construction of mRNA prognosis signature associated with differentially expressed genes in early stage of stomach adenocarcinomas based on TCGA and GEO datasets.

Stomach adenocarcinomas (STAD) are the most common malignancy of the human digestive system and represent the fourth leading cause of cancer-related deaths. As early-stage STAD are generally mild or asymptomatic, patients with advanced STAD have short overall survival. Early diagnosis of STAD has a considerable influence on clinical outcomes. The mRNA expression data and clinical indicators of STAD and normal tissues were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database. The gene expression differences were analyzed by R packages, and gene function enrichment analysis was performed. Kaplan-Meier method and univariate Cox proportional risk regression analysis were used to screen differential expressed genes (DEGs) related to survival of STAD patients. Multivariate Cox proportional risk regression analysis was used to further screen and determine the prognostic DEGs in STAD patients, and to construct a multigene prognostic prediction signature. The accuracy of predictive signature was tested by receiver operating characteristic (ROC) curve software package, and the nomogram of patients with STAD was drawn. Cox regression was used to investigate the correlation between multigene prognostic signature and clinical factors. The predictive performance of this model was compared with two other models proposed in previous studies using KM survival analysis, ROC curve analysis, Harrell consistency index and decision curve analysis (DCA). qRT-PCR and Western blot were used to verify the expression levels of prognostic genes. The pathways and functions of possible involvement of features were predicted using the GSEA method. A total of 569 early-stage specific DEGs were retrieved from TCGA-STAD dataset, including 229 up-regulated genes and 340 down-regulated genes. Enrichment analysis showed that the early-stage specific DEGs were associated with cytokine-cytokine receptor interaction, neuroactive ligand-receptor interaction, and calcium signaling pathway. Multiple Cox regression algorithm was used to identify 10 early-stage specific DEGs associated with overall survival (P < 0.01) of STAD patients, and a multi-mRNA prognosis signature was established. The patients were divided into high-risk group and low-risk group according to the risk score. In the training set, the prognostic signature was positively correlated with tumor size and stage (P < 0.05), survival curve (P < 0.001) and time-dependent ROC (AUC = 0.625). In the training dataset and test dataset, the both signatures had good predictive efficiencies. Cox regression and DCA analysis revealed that the prognostic signature was an independent factor and had a better predict effect than the conventional TNM stage classification method and the earlier published biomarkers on the prognosis of STAD patients. In this study, based on the early-stage specifically expressed genes, the prognostic signature constructed through TCGA and GEO datasets may become an indicator for clinical prognosis assessment of STAD and a new strategy for targeted therapy in the future.

Jiang F ，Lin H ，Yan H ，Sun X ，Yang J ，Dong M ... -  《-》

Identification of a necroptosis-related prognostic gene signature associated with tumor immune microenvironment in cervical carcinoma and experimental verification.

Cervical carcinoma (CC) has been associated with high morbidity, poor prognosis, and high intratumor heterogeneity. Necroptosis is the significant cellular signal pathway in tumors which may overcome tumor cells' apoptosis resistance. To investigate the relationship between CC and necroptosis, we established a prognostic model based on necroptosis-related genes for predicting the overall survival (OS) of CC patients. The gene expression data and clinical information of cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC) patients were obtained from The Cancer Genome Atlas (TCGA). We identified 43 differentially expressed necroptosis-related genes (NRGs) in CESC by examining differential gene expression between CESC tumors and normal tissues, and 159 NRGs from the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Gene ontology (GO) and KEGG enrichment analysis illustrated that the genes identified were mainly related to cell necrosis, extrinsic apoptosis, Influenza A, I - kappaB kinase/NF - kappaB, NOD - like receptor, and other signaling pathways. Subsequently, least absolute shrinkage and selection operator (LASSO) regression and univariate and multivariate Cox regression analyses were used to screen for NRGs that were correlated with patient prognosis. A prognostic signature that includes CAMK2A, CYBB, IL1A, IL1B, SLC25A5, and TICAM2 was established. Based on the prognostic model, patients were stratified into either the high-risk or low-risk subgroups with distinct survival. Receiver operating characteristic (ROC) curve analysis was used to identify the predictive accuracy of the model. In relation to different clinical variables, stratification analyses were performed to demonstrate the associations between the expression levels of the six identified NRGs and the clinical variables in CESC. Immunohistochemical (IHC) validation experiments explored abnormal expressions of these six NRGs in CESC. We also explored the relationship between risk score of this necroptosis signature and expression levels of some driver genes in TCGA CESC database and Gene Expression Omnibus (GEO) datasets. Significant relationships between the six prognostic NRGs and immune-cell infiltration, chemokines, tumor mutation burden (TMB), microsatellite instability (MSI), and immune checkpoints in CESC were discovered. In conclusion, we successfully constructed and validated a novel NRG signature for predicting the prognosis of CC patients and might also play a crucial role in the progression and immune microenvironment in CC.

Sun K ，Huang C ，Li JZ ，Luo ZX ... -  《World Journal of Surgical Oncology》

A multi-omics-based investigation of the prognostic and immunological impact of necroptosis-related mRNA in patients with cervical squamous carcinoma and adenocarcinoma.

Necroptosis is a kind of programmed necrosis mode that plays a double-edged role in tumor progression. However, the role of necroptosis-related Messenger RNA (mRNA) in predicting the prognosis and immune response of cervical squamous carcinoma and adenocarcinoma (CESC) has not been fully studied. Firstly, the incidence of somatic mutation rate and copy number variation for 74 necroptosis-related mRNAs (NRmRNAs) were analyzed. Secondly, CESC patients were divided into four stable clusters based on the consensus clustering results and analyzed for correlations with a series of clinical factors. Subsequently, a total of 291 The Cancer Genome Atlas samples were randomly divided into either training or validation cohorts. A Cox proportional hazard model consisting of three NRmRNAs (CXCL8, CLEC9A, and TAB2) was constructed by univariate, least absolute shrinkage and selection operator and multivariate COX regression analysis to identify the prognosis and immune response. Its performance and stability were further validated in another testing dataset (GSE44001) from Gene Expression Omnibus database. The results of the receiver operating characteristic curve, principal component analysis, t-SNE, and nomogram indicated that the prognostic model we constructed can serve as an independent prognostic factor. The combination of the prognostic model and the classic TNM staging system could improve the performance in predicting the survival of CESC patients. In addition, differentially expressed genes from high and low-risk patients are screened by R software for functional analysis and pathway enrichment analysis. Besides, single-sample gene set enrichment analysis revealed that tumor-killing immune cells were reduced in the high-risk group. Moreover, patients in the low-risk group are more likely to benefit from immune checkpoint inhibitors. The analysis of tumor immune dysfunction and exclusion scores, M6A-related genes, stem cell correlation and Tumor mutational burden data with clinical information has quantified the expression levels of NRmRNAs between the two risk subgroups. According to tumor immune microenvironment scores, Spearman's correlation analysis, and drug sensitivity, immunotherapy may have a higher response rate and better efficacy in patients of the low-risk subgroup. In conclusion, we have reported the clinical significance of NRmRNAs for the prognosis and immune response in CESC patients for the first time. Screening of accurate and effective prognostic markers is important for designing a multi-combined targeted therapeutic strategy and the development of individualized precision medicine.

Zou J ，Lin Z ，Jiao W ，Chen J ，Lin L ，Zhang F ，Zhang X ，Zhao J ... -  《Scientific Reports》