Exploration of a Robust and Prognostic Immune Related Gene Signature for Cervical Squamous Cell Carcinoma.

Background: Cervical squamous cell carcinoma (CESC) is one of the most frequent malignancies in women worldwide. The level of immune cell infiltration and immune-related genes (IRGs) can significantly affect the prognosis and immunotherapy of CESC patients. Thus, this study aimed to identify an immune-related prognostic signature for CESC. Methods: TCGA-CESC cohorts, obtained from TCGA database, were divided into the training group and testing group; while GSE44001 dataset from GEO database was viewed as external validation group. ESTIMATE algorithm was applied to evaluate the infiltration levels of immune cells of CESC patients. IRGs were screened out through weighted gene co-expression network analysis (WGCNA). A multi-gene prognostic signature based on IRGs was constructed using LASSO penalized Cox proportional hazards regression, which was validated through Kaplan-Meier, Cox, and receiver operating characteristic curve (ROC) analyses. The abundance of immune cells was calculated using ssGSEA algorithm in the ImmuCellAI database, and the response to immunotherapy was evaluated using immunophenoscore (IPS) analysis and the TIDE algorithm. Results: In TCGA-CESC cohorts, higher levels of immune cell infiltration were closely associated with better prognoses. Moreover, a prognostic signature was constructed using three IRGs. Based on this given signature, Kaplan-Meier analysis suggested the significant differences in overall survival (OS) and the ROC analysis demonstrated its robust predictive potential for CESC prognosis, further confirmed by internal and external validation. Additionally, multivariate Cox analysis revealed that the three IRGs signature served as an independent prognostic factor for CESC. In the three-IRGs signature low-risk group, the infiltrating immune cells (B cells, CD4/8 + T cells, cytotoxic T cells, macrophages and so on) were much more abundant than that in high-risk group. Ultimately, IPS and TIDE analyses showed that low-risk CESC patients appeared to present with a better response to immunotherapy and a better prognosis than high-risk patients. Conclusion: The present prognostic signature based on three IRGs (CD3E, CD3D, LCK) was not only reliable for survival prediction but efficient to predict the clinical response to immunotherapy for CESC patients, which might assist in guiding more precise individual treatment in the future.

Zuo Z ，Xiong J ，Zeng C ，Jiang Y ，Xiong K ，Tao H ，Guo Y ... -  《Frontiers in Molecular Biosciences》

A prognostic signature based on immune-related genes for cervical squamous cell carcinoma and endocervical adenocarcinoma.

Liu J ，Wu Z ，Wang Y ，Nie S ，Sun R ，Yang J ，Cheng W ... -  《-》

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》

Identification and validation of an immune prognostic signature in colorectal cancer.

Colorectal cancer (CRC) is one of the most common malignant neoplasms worldwide. Although the significant efficacy of immunotherapy has been shown, only limited CRC patients benefit from it. Therefore, we aimed to establish a prognostic signature based on immune-related genes (IRGs) to predict overall survival (OS) and the potential response to immunotherapy in CRC patients. Gene expression profiles and clinical information of CRC patients were collected from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. The prognostic signature composed of IRGs was established using univariate Cox regression and the least absolute shrinkage and selection operator (LASSO) regression analysis. CIBERSORT was used to estimate the immune cell infiltration. A total of 24 survival-related IRGs were identified from 247 differentially expressed IRGs. Then, 16 IRGs were selected to establish the prognostic signature that stratified the patients into the high-risk and low-risk groups with statistically different survival outcomes. The AUCs of the time-dependent ROC curves indicated that the signature had a strong predictive accuracy in internal and external validation sets. Multivariate cox regression analysis suggested that the signature could also act as an independent prognostic factor for OS. The low-risk group had a higher proportion of immune cell infiltration than the high-risk group, such as CD4 memory resting T cells, activated dendritic cells, and resting dendritic cells. In addition, patients in the high-risk group exhibited higher tumor mutation burden and BRAF mutation. We developed an immune-related prognostic signature to predict the OS and immune status in CRC patients. We believed that our signature is conducive to better stratification and more precise immunotherapy for CRC patients.

Li M ，Wang H ，Li W ，Peng Y ，Xu F ，Shang J ，Dong S ，Bu L ，Wang H ，Wei W ，Hu Q ，Liu L ，Zhao Q ... -  《-》

Comprehensive FGFR3 alteration-related transcriptomic characterization is involved in immune infiltration and correlated with prognosis and immunotherapy response of bladder cancer.

Bladder cancer (BC) threatens the health of human beings worldwide because of its high recurrence rate and mortality. As an actionable biomarker, fibroblast growth factor receptor 3 (FGFR3) alterations have been revealed as a vital biomarker and associated with favorable outcomes in BC. However, the comprehensive relationship between the FGFR3 alteration associated gene expression profile and the prognosis of BC remains ambiguous. Genomic alteration profile, gene expression data, and related clinical information of BC patients were downloaded from The Cancer Genomics database (TCGA), as a training cohort. Subsequently, the Weighted Gene Co-expression Network Analysis (WGCNA) was conducted to identify the hub modules correlated with FGFR3 alteration. The univariate, multivariate, and least absolute shrinkage and selection operator (LASSO) Cox regression analyses were used to obtain an FGFR3 alteration-related gene (FARG) prognostic signature and FARG-based nomogram. The receiver operating characteristic (ROC) curve analysis was used for evaluation of the ability of prognosis prediction. The FARG signature was validated in four independent datasets, namely, GSE13507, GSE31684, GSE32548, and GSE48075, from Gene Expression Omnibus (GEO). Then, clinical feature association analysis, functional enrichment, genomic alteration enrichment, and tumor environment analysis were conducted to reveal differential clinical and molecular characterizations in different risk groups. Lastly, the treatment response was evaluated in the immunotherapy-related dataset of the IMvigor210 cohort and the frontline chemotherapy dataset of GSE48276, and the chemo-drug sensitivity was estimated via Genomics of Drug Sensitivity in Cancer (GDSC). There were a total of eleven genes (CERCAM, TPST1, OSBPL10, EMP1, CYTH3, NCRNA00201, PCDH10, GAP43, COLQ, DGKB, and SETBP1) identified in the FARG signature, which divided BC patients from the TCGA cohort into high- and low-risk groups. The Kaplan-Meier curve analysis demonstrated that BC patients in the low-risk group have superior overall survival (OS) than those in the high-risk group (median OS: 27.06 months vs. 104.65 months, p < 0.0001). Moreover, the FARG signature not only showed a good performance in prognosis prediction, but also could distinguish patients with different neoplasm disease stages, notably whether patients presented with muscle invasive phenotype. Compared to clinicopathological features, the FARG signature was found to be the only independent prognostic factor, and subsequently, a FARG-based prognostic nomogram was constructed with better ability of prognosis prediction, indicated by area under ROC curve (AUC) values for 1-, 3-, and 5-year OS of 0.69, 0.71, and 0.79, respectively. Underlying the FARG signature, multiple kinds of metabolism- and immune-related signaling pathways were enriched. Genomic alteration enrichment further identified that FGFR3 alterations, especially c.746C>G (p.Ser249Cys), were more prevalent in the low-risk group. Additionally, FARG score was positively correlated with ESTIMATE and TIDE scores, and the low-risk group had abundant enrichment of plasma B cells, CD8+ T cells, CD4+ naive T cells, and helper follicular T cells, implying that patients in the low-risk group were likely to make significant responses to immunotherapy, which was further supported by the analysis in the IMvigor210 cohort as there was a significantly higher response rate among patients with lower FARG scores. The analysis of the GDSC database finally demonstrated that low-risk samples were more sensitive to methotrexate and tipifarnib, whereas those in the high-risk group had higher sensitivities in cisplatin, docetaxel, and paclitaxel, instead. The novel established FARG signature based on a comprehensive FGFR3 alteration-related transcriptomic profile performed well in prognosis prediction and was also correlated with immunotherapy and chemotherapy treatment responses, which had great potential in future clinical applications.

Xu T ，Xu W ，Zheng Y ，Li X ，Cai H ，Xu Z ，Zou Q ，Yu B ... -  《Frontiers in Immunology》