AUTHOR=Wang Nana , Nanding Abiyasi , Jia Xiaocan , Wang Yuping , Yang Chaojun , Fan Jingwen , Dong Ani , Zheng Guowei , Ma Jiaxin , Shi Xuezhong , Yang Yongli TITLE=Mining of immunological and prognostic-related biomarker for cervical cancer based on immune cell signatures JOURNAL=Frontiers in Immunology VOLUME=13 YEAR=2022 URL=https://www.frontiersin.org/journals/immunology/articles/10.3389/fimmu.2022.993118 DOI=10.3389/fimmu.2022.993118 ISSN=1664-3224 ABSTRACT=Background

Immunotherapy has changed the therapeutic landscape of cervical cancer (CC), but has durable anti-tumor activity only in a subset of patients. This study aims to comprehensively analyze the tumor immune microenvironment (TIME) of CC and to mine biomarkers related to immunotherapy and prognosis.

Methods

The Cancer Genome Atlas (TCGA) data was utilized to identify heterogeneous immune subtypes based on survival-related immune cell signatures (ICSs). ICSs prognostic model was constructed by Cox regression analyses, and immunohistochemistry was conducted to verify the gene with the largest weight coefficient in the model. Meanwhile, the tumor immune infiltration landscape was comprehensively characterized by ESTIMATE, CIBERSORT and MCPcounter algorithms. In addition, we also analyzed the differences in immunotherapy-related biomarkers between high and low-risk groups. IMvigor210 and two gynecologic tumor cohorts were used to validate the reliability and scalability of the Risk score.

Results

A total of 291 TCGA-CC samples were divided into two ICSs clusters with significant differences in immune infiltration landscape and prognosis. ICSs prognostic model was constructed based on eight immune-related genes (IRGs), which showed higher overall survival (OS) rate in the low-risk group (P< 0.001). In the total population, time-dependent receiver operating characteristic (ROC) curves displayed area under the curve (AUC) of 0.870, 0.785 and 0.774 at 1-, 3- and 5-years. Immunohistochemical results showed that the expression of the oncogene (FKBP10) was negatively correlated with the degree of differentiation and positively correlated with tumor stage, while the expression of tumor suppressor genes (S1PR4) was the opposite. In addition, the low-risk group had more favorable immune activation phenotype and higher enrichment of immunotherapy-related biomarkers. The Imvigor210 and two gynecologic tumor cohorts validated a better survival advantage and immune efficacy in the low-risk group.

Conclusion

This study comprehensively assessed the TIME of CC and constructed an ICSs prognostic model, which provides an effective tool for predicting patient’s prognosis and accurate immunotherapy.