- Primary research
- Open Access
Identification of a tumor microenvironment-related gene signature to improve the prediction of cervical cancer prognosis
Cancer Cell International volume 21, Article number: 182 (2021)
Previous studies have found that the microenvironment of cervical cancer (CESC) affects the progression and treatment of this disease. Thus, we constructed a multigene model to assess the survival of patients with cervical cancer.
We scored 307 CESC samples from The Cancer Genome Atlas (TCGA) and divided them into high and low matrix and immune scores using the ESTIMATE algorithm for differential gene analysis. Cervical cancer patients were randomly divided into a training group, testing group and combined group. The multigene signature prognostic model was constructed by Cox analyses. Multivariate Cox analysis was applied to evaluate the significance of the multigene signature for cervical cancer prognosis. Prognosis was assessed by Kaplan–Meier curves comparing the different groups, and the accuracy of the prognostic model was analyzed by receiver operating characteristic-area under the curve (ROC-AUC) analysis and calibration curve. The Tumor Immune Estimation Resource (TIMER) database was used to analyze the relationship between the multigene signature and immune cell infiltration.
We obtained 420 differentially expressed genes in the tumor microenvironment from 307 patients with cervical cancer. A three-gene signature (SLAMF1, CD27, SELL) model related to the tumor microenvironment was constructed to assess patient survival. Kaplan–Meier analysis showed that patients with high risk scores had a poor prognosis. The ROC-AUC value indicated that the model was an accurate predictor of cervical cancer prognosis. Multivariate cox analysis showed the three-gene signature to be an independent risk factor for the prognosis of cervical cancer. A nomogram combining the three-gene signature and clinical features was constructed, and calibration plots showed that the nomogram resulted in an accurate prognosis for patients. The three-gene signature was associated with T stage, M stage and degree of immune infiltration in patients with cervical cancer.
This research suggests that the developed three-gene signature may be applied as a biomarker to predict the prognosis of and personalized therapy for CESC.
Despite expanding the screening scale of cervical cancer (CESC) and promoting the use of vaccines, there were 569,847 new cases of CESC and 311,365 deaths from CESC worldwide in 2018. The incidence and mortality rate of CESC were the fourth highest among all cancers in women, suggesting a serious societal burden . High-risk human papillomavirus (HR-HPV) persistent infection is the main cause of cervical cancer development . Cervical squamous cell carcinoma and cervical adenocarcinoma are the main types of cervical cancer. At present, the conventional treatment of cervical cancer includes radiotherapy, chemotherapy and surgery, but patients at advanced stages are prone to developing radiotherapy and chemotherapy resistance. Many patients with cervical cancer are already in the advanced stage at the time of diagnosis, which is often accompanied by high invasion rates, and the 3-year or 5-year mortality rate is between 52 and 79% [3,4,5]. Although some diagnostic and prognostic markers for cervical cancer have been discovered, the overall efficacy of these markers remains a challenge because of the heterogeneity of cervical cancer. Therefore, it is necessary to find prognostic markers for cervical cancer at the molecular level, which will offer new insight into the treatment and prognosis of cervical cancer.
Previous articles have revealed that the balance of the host tumor microenvironment (TME) plays a vital role in the process of cancers. The TME consists of fibroblasts, vascular endothelial cells, immune cells, cytokines, growth factors, hormones, extracellular matrix components, and so on . The TME affects invasion, migration, proliferation, apoptosis, and response to drugs in CESC cells [7,8,9,10,11]. The Estimation of Stromal and Immune cells in MAlignant Tumor tissues using Expression data (ESTIMATE) algorithm uses gene expression data to analyze the proportion of immune cells and stromal cells in the TME . The algorithm is based on single-sample gene set enrichment analysis (ssGSEA), and the immune score, matrix score and tumor purity score are obtained. TME scores have been reported to aid in the prognosis of glioblastoma and colorectal cancer patients [13, 14].
As the largest cancer gene database at present, The Cancer Genome Atlas (TCGA) database has rich sequencing and microarray data, including gene expression, miRNA expression, copy number variation, DNA methylation, SNP and standardized clinical data. Previous studies have found that multigene models can predict tumor prognosis base on TCGA. In breast cancer, Qi et al. found that a prognostic model of a four-gene signature predicted the survival of breast cancer patients . A recent article involving papillary thyroid cancer showed that a prognostic model of an seven-gene signature predicted patient outcomes and was associated with stage and metastasis . In addition, He et al. indicated a prognostic model of an eight-gene signature to be an independent risk factor for the prognosis of lung cancer, providing new ideas for accurate treatment of patients . Multigene models are also important for individualized systemic treatment, prolonged survival and early detection of cervical cancer. A recent study showed that a 70-gene signature can be used as a predictor of the therapeutic effect of individualized systemic treatment of advanced-stage cervical cancer, which can help prolong the survival time of patients . Xie et al. found that an 8-gene signature can predict the prognosis of patients with cervical cancer after radiotherapy . Cai et al. reported that a glycolysis-related six-gene signature can be used as a prognostic marker to accurately predict the prognosis of cervical cancer patients . Ma et al. found that a 4-gene signature in peripheral plasma is a new marker for the diagnosis of cervical cancer and has good diagnostic performance . However, a multigene model associated with the tumor microenvironment (TME) has not yet been reported for assessing survival in cervical cancer.
In this study, we used the ESTIMATE algorithm to determine the immune and matrix scores of CESC samples to identify hub genes in the TME. We constructed a prognostic model of multigene signatures by performing a multivariate Cox regression analysis of prognosis-associated hub genes in the TME using the ‘‘survival’’ package in R software. According to our results, the three-gene signature (SLAMF1, CD27, SELL) may be applied as a biomarker to predict prognosis and to provide new insight for therapy.
We obtained mRNA expression data and patient clinical information from the TCGA database (https://cancergenome.nih.gov/). The clinical data included pathological tumor type, patient age, tumor stage, tumor grade, follow-up time and status. The exclusion criteria were set as follows: (1) histologic diagnosis is not CESC; (2) samples lack gene expression profile for analysis; and (3) clinical information is not available for analysis. Overall, a total of 307 CESC patients from the TCGA database were included in our study. The patients were randomly divided into a training group and a testing group at a ratio of 1:1.
Screening of differentially expressed genes in the TME
The ratio of immune to stromal cells in the TME was calculated using the ESTIMATE algorithm, and the samples were separated into high- and low-score groups based on the median value for differential analysis. Gene expression values were analyzed using the limma package in R 3.5.3. The expression value of each gene was log2-converted for further study, and the differential genes with a log2-fold change (log2FC) > 1 and a false discovery rate (FDR) < 0.05 were considered significant.
Biological functionality analyses
Gene Ontology (GO) functional annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis are commonly used to assess the functions of genes. We used the clusterProfiler package in R to perform GO and KEGG analyses on the differential genes in the microenvironment. The STRING 11.0 database was utilized to discover protein–protein interactions (PPIs), and the Molecular Complex Detection (MCODE) 1.5.1 plugin was used to find core genes . The PPI network was visualized using Cytoscape 3.7.1 .
Analysis of immune cell infiltration
The Tumor Immune Estimation Resource (TIMER) system can be found at https://cistrome.shinyapps.io/timer/. The TIMER database system was used to analyze the immune infiltration of the samples from the TCGA database, including 10,897 tissues and 32 cancer types . We used the TIMER database to study the relationship between the risk score of three-gene signature and immune cell infiltration in CESC tissues.
Construction of a multigene signature prognostic model by Cox model analysis was carried out using the ‘‘survival’’ package in R software. The gene signature can be predicted as follows: risk score = ∑βmRNAi × ExpmRNAi (where β is the coefficient and Exp is expression of the gene). Survival analysis was used to compare the relationship between different risk scores (high score vs low score) and prognosis. The accuracy of the prognostic evaluation of the risk score of the three-gene signature was evaluated based on the ROC analysis. The same formula was used to calculate the risk score in the testing group and the combined group (training group + testing group) to verify the accuracy of the model. Multivariate cox analysis of the multigene signature was performed to assess whether it is an independent prognostic risk factor. Nomogram analysis and calibration plot were conducted using the “rms” package. The association between the risk score and immune cell infiltration was analyzed by Spearman correlation. The data were analyzed using R software. Statistical significance was set as a P-value < 0.05.
Matrix score and immune score of the CESC microenvironment
We downloaded the clinical information and gene expression matrices of 307 CESC patients from TCGA. The mean age of these patients was 48 years. A total of 254 patients were diagnosed with squamous cell carcinoma (SCC), 46 patients were diagnosed with adenocarcinoma (AC), and 6 patients were diagnosed with adenosquamous carcinoma (ASC). To study the ratio of stromal cells and immune cells in the TME, we used the ESTIMATE algorithm to analyze gene expression data. The matrix score ranged from − 2437.39 to 804.22, and the average value was -868.78; the immune score ranged from − 1209.74 to 3419.33, and the average value was 828.83. The results showed that SCC had the highest matrix score, followed by ASC, and AC had the lowest matrix score (Fig. 1a, P = 0.002). The order of the immune score from high to low was SCC > AC > ASC (Fig. 1b, P = 2.72E−09), suggesting that the pathological type of CESC affects the matrix score and the immune score.
Expression of genes in different matrix and immune score groups
We included 307 patients in the cohort to study the expression of genes in different stromal and immune score groups. Genes with a fold change > 2 and an FDR < 0.05 were included in the subsequent analysis. As a result, 1029 genes were highly expressed in the high matrix-score group, and 10 genes were expressed at low levels (Fig. 2a, c). There were 625 genes with high expression in the high immune-score group and 260 genes with low expression (Fig. 2b, d). The Venn diagram shows that a total of 416 genes were upregulated in the stromal/immune high-score group, and 4 genes were downregulated in the stromal/immune high-score group. These genes were considered to be differentially expressed genes in the microenvironment (Fig. 2e, f).
GO, KEGG, and PPI analyses of differentially expressed genes
To study the role of these genes in the CESC microenvironment, we used GO and KEGG analyses. The biological process (BP) category contained T cell activation and leukocyte cell–cell adhesion. T cell receptor complex and immunological synapse were included in the GO cellular component (CC) category. The molecular function (MF) category included cytokine receptor activity and CCR chemokine receptor binding (Table 1; Fig. 3). The results of the KEGG pathway enrichment are displayed in Table 2; the genes were mainly enriched in immune-related pathways, including the T cell receptor signaling pathway, intestinal immune network for IgA production and B cell receptor signaling pathway (Fig. 4). The PPI network analysis was based on the STRING database and included 392 nodes and 7228 edges. Then, we used the MCODE plugin in Cytoscape to find the core genes. A total of 59 genes were considered to be core genes and were used for subsequent analysis (Fig. 5).
Identification of the three-gene signature characteristics
In the training group, we performed univariate Cox regression analysis to identify the core genes associated with overall survival, which showed that six core genes were associated with overall survival (P < 0.05). Then, we constructed a prognostic model of multigene signatures by performing a multivariate Cox regression analysis of prognosis-associated core genes using the ‘‘survival’’ package in R software. The results show that the three-gene signature (SLAMF1, CD27, SELL) can be predicted as follows: Risk Score = [SLAMF1 × (0.7380)] + [CD27 × (− 0.1894)] + [SELL × (− 0.1008)] (Table 3). The patient’s risk score was evaluate based on expression of the three genes, and the patients were divided into groups of high and low scores according to the median value (Fig. 6). The results showed that the overall survival time of the high-risk score patients was shorter than that of the low-risk score patients; the 5-year survival rate of the high-risk score patients was 55.3%, and that of the low-risk score patients was 76.1% (Fig. 6d). ROC analysis showed that the three-gene signature can accurately predict the survival of patients with cervical cancer in the training group (AUC = 0.785, Fig. 6e). In addition, we performed verification analysis of the three-gene signature in the testing group and combined group. Survival analysis showed that the 5-year survival rate in the high-risk group was lower than that in the low-risk group (Fig. 6d). The AUCs for the testing group and combined group were 0.714 and 0.736, respectively (Fig. 6e). Prognostic analysis also showed that the three-gene signature can predict prognosis among patients at different ages and stages (Fig. 7a–f).
The three-gene signature is an independent risk factor of prognosis
Univariate analysis showed that TNM stage, T stage and risk score were associated with prognosis in patients with cervical cancer (Fig. 8a). After adjusting for age, grade, TNM stage, T stage, M stage, N stage, and radiation therapy of cervical cancer patients, multivariate analysis showed that the risk score remained an independent prognostic risk factor (Fig. 8b, HR = 6.267, 95% CI = 2.060–19.067, P = 0.001). Finally, we constructed a nomogram based on clinical features and risk score to assess the prognosis of cervical cancer patients (Fig. 9a), whereby a high score indicates that the patient has a poor prognosis. Calibration plots showed that the nomogram resulted in an accurate prognosis for patients with cervical cancer (Fig. 9b–d).
Relationship between the three-gene signature and clinicopathological features
We analyzed the risk scores and clinical features of patients with cervical cancer, including age, grade, TNM stage, T stage, M stage, and N stage. The results showed a higher risk score in the high-grade T stage (P = 0.027) and metastatic tissues (P = 0.018) (Fig. 10). We also examined the relationship between risk score and immune cell infiltration levels and found that the risk score was significantly negatively correlated with the degree of infiltration of B cells, CD8 + T cells, CD4 + T cells, macrophages, neutrophils, and dendritic cells in CESC (Fig. 11).
CESC is one of the main causes of cancer death in females, especially in sub-Saharan Africa and Southeast Asia, with a mortality rate of more than 50%, which is 18 times that of developed countries [1, 3]. Because of the poor prognosis of cervical cancer, it is important to find molecular markers to predict prognosis. With the development of bioinformatics technology, an increasing number of researchers are evaluating the significance of multigene signatures for the prognosis of different tumors through open databases. For example, Yang et al. found that a seven-gene model can accurately predict the prognosis of head and neck squamous cell carcinoma . Using lung cancer data from TCGA, a seven-gene model was able to predict prognosis . An eight-gene prognostic model was reported to have predictive power for cervical cancer prognosis . Ding et al. identified a potential three-gene signature (MSRB3, CENPM, ZIC2) that can predict the prognosis of cervical cancer through the TCGA dataset, which is expected to provide a basis for the prognostic evaluation of patients in the future . The tumor microenvironment affects the prognosis and treatment of tumors, especially with regard to immunotherapy, though a multigene model associated with the tumor microenvironment (TME) has not been reported for cervical cancer prognosis.
In this study, we downloaded mRNA expression data and clinical follow-up information for 307 patients with cervical cancer from TCGA. We divided the cervical cancer samples into high and low immunity and matrix scores using the ESTIMATE algorithm and obtained 420 differential genes in the cervical cancer microenvironment. GO and KEGG analyses revealed differential gene enrichment in immune-related functions and pathways. In addition, we constructed a three-gene signature (SLAMF1, CD27, SELL) to predict prognosis based on Cox analysis in the training group. The overall survival time of the high-risk score group was shorter than that of the low-risk score group, including patients in different subgroups (age and stage). The three-gene signature was able to accurately predict the prognosis of cervical cancer. In addition, the results were further verified in the testing and combined groups. Multivariate Cox analysis showed that the three-gene signature was an independent risk factor for cervical cancer prognosis. Based on age, grade, TNM stage, T stage, M stage, N stage and risk score, a nomogram was constructed to score prognosis, and calibration plots showed that the nomogram could accurately predict patient survival. Furthermore, the results showed that the three-gene signature scored high for the high-grade T stage and metastatic tissue and was associated with the degree of immune infiltration in the TME.
SLAMF1, also known as CD150, belongs to the SLAM receptor family. SLAMF1 is commonly expressed in the immune system, and it was recently found to be expressed in tumors of the central nervous system . A previous study showed that a two-gene prognostic model can accurately predict the prognosis of patients with chronic lymphocytic leukemia, in which SLAMF1 is a risk gene, and patients with high risk scores had a poorer prognosis than did those with low risk scores . In breast cancer, Lin et al. suggested associations between allelic expression imbalance polymorphisms and breast cancer risk (SLAMF1 rs1061217) . CD27, a member of the tumor necrosis factor receptor family, is present in CD4, CD8 T lymphocytes and NK cells, plays a significant role in tumor immunotherapy [31, 32]. In multiple myeloma, Kaplan–Meier analysis showed that patients with high CD27 expression had a longer overall survival time than patients with low CD27 expression . SELL, also called CD62L, is located on the surface of T cells and promotes T cell homing to lymph nodes. Vlad et al. revealed that CD62L is downregulated in chronic lymphocytic leukemia and associated with disease progression . Additionally, CD62L can enhance the efficacy of tumor immunotherapy for cancer treatment . In bladder cancer, the serum level of SELL has been found to be useful as a tool for diagnostic staging and grading . The results of these studies were consistent with our findings in CESC.
To investigate the mechanisms in the CESC microenvironment, we performed GO and KEGG analysis for the differentially expressed genes, as well as PPI network analyses and immune infiltration correlation analysis. The results revealed that genes were concentrated on immune-related pathways, such as the T cell receptor signaling pathway, the intestinal immune network for IgA production, and the B cell receptor signaling pathway. T cell receptor signaling has been reported to play a role in the progression of extranodal NK/T cell lymphoma and to activate the ITK/NF-κB/GATA-3 axis to promote chemoresistance in T cell lymphomas [37, 38]. Yang et al. found that the intestinal immune network for IgA production affects the cell proliferation and migration of liver cancer cells . Krysiak et al. found that the B cell receptor signaling pathway affects the biology, prognosis, and treatment of follicular lymphoma . Another study revealed that B cell receptor signaling affects the prognosis of patients with chronic lymphocytic leukemia . Based on these results, the hub genes are speculated to play an important role in the CESC microenvironment.
However, there are still some limitations in our research. Although the prognostic significance of the three-gene signature in patients with cervical cancer was analyzed by bioinformatics, additional samples and external datasets will need to be further studied and verified. In our future work, we will externally verify the prognostic accuracy of the three-gene signature in other databases and local patient data via PCR. Functional testing is needed, and determining the mechanisms of the TME requires further research. Our research team is currently working towards some of these goals.
For the first time, we analyzed CESC sample data using the ESTIMATE algorithm and found that a three-gene signature consisting of core genes in the cervical cancer microenvironment can be used to predict patient prognosis. In addition, the three-gene signature is an independent risk factor for the prognosis of cervical cancer. Gene functional analysis revealed associations with immune responses. Indeed, this research indicates that three-gene signature may be applied as biomarkers to predict the prognosis and therapeutic targets of CESC.
Availability of data and materials
Research Data are available at TCGA.
High-risk human papillomavirus
The cancer genome atlas
Tumor immune estimation resource
Receiver operating characteristic
Area under the curve
Kyoto encyclopedia of genes and genomes
Molecular complex detection
False discovery rate
Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.
Wang X, Huang X, Zhang Y. Involvement of human papillomaviruses in cervical cancer. Front Microbiol. 2018;9:2896.
Sankaranarayanan R, Swaminathan R, Brenner H, Chen K, Chia KS, Chen JG, Law SCK, Ahn Y-O, Xiang YB, Yeole BB, et al. Cancer survival in Africa, Asia, and Central America: a population-based study. Lancet Oncol. 2010;11(2):165–73.
Small W, Bacon MA, Bajaj A, Chuang LT, Fisher BJ, Harkenrider MM, Jhingran A, Kitchener HC, Mileshkin LR, Viswanathan AN, et al. Cervical cancer: a global health crisis. Cancer. 2017;123(13):2404–12.
Wakatsuki M, Kato S, Kiyohara H, Ohno T, Karasawa K, Tamaki T, Ando K, Shiba S, Kamada T, Nakano T. The prognostic value of rectal invasion for stage IVA uterine cervical cancer treated with radiation therapy. BMC Cancer. 2016;16:244.
Wu T, Dai Y. Tumor microenvironment and therapeutic response. Cancer Lett. 2017;387:61–8.
Mao Y, Dong L, Zheng Y, Dong J, Li X. Prediction of recurrence in cervical cancer using a nine-lncRNA signature. Front Genet. 2019;10:284.
Yadav SS, Prasad SB, Prasad CB, Pandey LK, Pradhan S, Singh S, Narayan G. CXCL12 is a key regulator in tumor microenvironment of cervical cancer: an in vitro study. Clin Exp Metas. 2016;33(5):431–9.
Han S, Wang Y, Shi X, Zong L, Liu L, Zhang J, Qian Q, Jin J, Ma Y, Cui B, et al. Negative roles of B7–H3 and B7–H4 in the microenvironment of cervical cancer. Exp Cell Res. 2018;371(1):222–30.
Yadav SS, Prasad CB, Prasad SB, Pandey LK, Singh S, Pradhan S, Narayan G. Anti-tumor activity of staurosporine in the tumor microenvironment of cervical cancer: an in vitro study. Life Sci. 2015;133:21–8.
Dou YN, Chaudary N, Chang MC, Dunne M, Huang H, Jaffray DA, Milosevic M, Allen C. Tumor microenvironment determines response to a heat-activated thermosensitive liposome formulation of cisplatin in cervical carcinoma. J Contr Rel Off J Contr Rel Soc. 2017;262:182–91.
Yoshihara K, Shahmoradgoli M, Martinez E, Vegesna R, Kim H, Torres-Garcia W, Trevino V, Shen H, Laird PW, Levine DA, et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun. 2013;4:2612.
Jia D, Li S, Li D, Xue H, Yang D, Liu Y. Mining TCGA database for genes of prognostic value in glioblastoma microenvironment. Aging. 2018;10(4):592–605.
Liu W, Ye H, Liu YF, Xu CQ, Zhong YX, Tian T, Ma SW, Tao H, Li L, Xue LC, et al. Transcriptome-derived stromal and immune scores infer clinical outcomes of patients with cancer. Oncology Lett. 2018;15(4):4351–7.
Qi L, Yao Y, Zhang T, Feng F, Zhou C, Xu X, Sun C. A four-mRNA model to improve the prediction of breast cancer prognosis. Gene. 2019;721:144100.
Lin P, Guo YN, Shi L, Li XJ, Yang H, He Y, Li Q, Dang YW, Wei KL, Chen G. Development of a prognostic index based on an immunogenomic landscape analysis of papillary thyroid cancer. Aging. 2019;11(2):480–500.
He R, Zuo S. A robust 8-gene prognostic signature for early-stage non-small cell lung cancer. Front Oncol. 2019;9:693.
Nguyen NNY, Choi TG, Kim J, Jung MH, Ko SH, Shin Y, Kang I, Ha J, Kim SS, Jo YH. A 70-gene signature for predicting treatment outcome in advanced-stage cervical cancer. Mol Ther Oncol. 2020;19:47–56.
Xie F, Dong D, Du N, Guo L, Ni W, Yuan H, Zhang N, Jie J, Liu G, Tai G. An 8-gene signature predicts the prognosis of cervical cancer following radiotherapy. Mol Med Rep. 2019;20(4):2990–3002.
Cai L, Hu C, Yu S, Liu L, Yu X, Chen J, Liu X, Lin F, Zhang C, Li W, et al. Identification and validation of a six-gene signature associated with glycolysis to predict the prognosis of patients with cervical cancer. BMC Cancer. 2020;20(1):1133.
Ma G, Song G, Zou X, Shan X, Liu Q, Xia T, Zhou X, Zhu W. Circulating plasma microRNA signature for the diagnosis of cervical cancer. Cancer Biomarkers Sect A Dis Markers. 2019;26(4):491–500.
Bader GD, Hogue CW. An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinformatics. 2003;4:2.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.
Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, Li B, Liu XS. TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells. Can Res. 2017;77(21):e108–10.
Yang B, Fu L, Xu S, Xiao J, Li Z, Liu Y. A nomogram based on a gene signature for predicting the prognosis of patients with head and neck squamous cell carcinoma. Int J Biol Markers. 2019;34(3):309–17.
Wang Z, Wang Z, Niu X, Liu J, Wang Z, Chen L, Qin B. Identification of seven-gene signature for prediction of lung squamous cell carcinoma. OncoTargets Ther. 2019;12:5979–88.
Ding TT, Ma H, Feng JH. A three-gene novel predictor for improving the prognosis of cervical cancer. Oncol Lett. 2019;18(5):4907–15.
Romanets-Korbut O, Najakshin AM, Yurchenko M, Malysheva TA, Kovalevska L, Shlapatska LM, Zozulya YA, Taranin AV, Horvat B, Sidorenko SP. Expression of CD150 in tumors of the central nervous system: identification of a novel isoform. PLoS ONE. 2015;10(2):e0118302.
Schweighofer CD, Coombes KR, Barron LL, Diao L, Newman RJ, Ferrajoli A, O’Brien S, Wierda WG, Luthra R, Medeiros LJ, et al. A two-gene signature, SKI and SLAMF1, predicts time-to-treatment in previously untreated patients with chronic lymphocytic leukemia. PLoS ONE. 2011;6(12):e28277.
Lin W, Lin HD, Guo XY, Lin Y, Su FX, Jia WH, Tang LY, Zheng W, Long JR, Ren ZF. Allelic expression imbalance polymorphisms in susceptibility chromosome regions and the risk and survival of breast cancer. Mol Carcinog. 2017;56(1):300–11.
Han BK, Olsen NJ, Bottaro A. The CD27-CD70 pathway and pathogenesis of autoimmune disease. Semin Arthritis Rheum. 2016;45(4):496–501.
Bullock TN. Stimulating CD27 to quantitatively and qualitatively shape adaptive immunity to cancer. Curr Opin Immunol. 2017;45:82–8.
Alaterre E, Raimbault S, Goldschmidt H, Bouhya S, Requirand G, Robert N, Boireau S, Seckinger A, Hose D, Klein B, et al. CD24, CD27, CD36 and CD302 gene expression for outcome prediction in patients with multiple myeloma. Oncotarget. 2017;8(58):98931–44.
Vlad A, Deglesne PA, Letestu R, Saint-Georges S, Chevallier N, Baran-Marszak F, Varin-Blank N, Ajchenbaum-Cymbalista F, Ledoux D. Down-regulation of CXCR4 and CD62L in chronic lymphocytic leukemia cells is triggered by B-cell receptor ligation and associated with progressive disease. Can Res. 2009;69(16):6387–95.
Watson HA, Durairaj RRP, Ohme J, Alatsatianos M, Almutairi H, Mohammed RN, Vigar M, Reed SG, Paisey SJ, Marshall C, et al. L-selectin enhanced T cells improve the efficacy of cancer immunotherapy. Front Immunol. 2019;10:1321.
Phadke GS, Satterwhite-Warden JE, Choudhary D, Taylor JA, Rusling JF. A novel and accurate microfluidic assay of CD62L in bladder cancer serum samples. Analyst. 2018;143(22):5505–11.
Miyata-Takata T, Chuang SS, Takata K, Toji T, Maeda Y, Sato Y, Yoshino T. Expression of T-cell receptor signalling pathway components in extranodal NK/T-cell lymphoma. Histopathology. 2018;73(6):1030–8.
Wang T, Lu Y, Polk A, Chowdhury P, Murga-Zamalloa C, Fujiwara H, Suemori K, Beyersdorf N, Hristov AC, Lim MS, et al. T-cell receptor signaling activates an ITK/NF-kappaB/GATA-3 axis in T-cell lymphomas facilitating resistance to chemotherapy. Clin Cancer Res Off J Am Assoc Cancer Res. 2017;23(10):2506–15.
Yang Z, Tao Y, Xu X, Cai F, Yu Y, Ma L. Bufalin inhibits cell proliferation and migration of hepatocellular carcinoma cells via APOBEC3F induced intestinal immune network for IgA production signaling pathway. Biochem Biophys Res Commun. 2018;503(3):2124–31.
Krysiak K, Gomez F, White BS, Matlock M, Miller CA, Trani L, Fronick CC, Fulton RS, Kreisel F, Cashen AF, et al. Recurrent somatic mutations affecting B-cell receptor signaling pathway genes in follicular lymphoma. Blood. 2017;129(4):473–83.
Dezorella N, Katz BZ, Shapiro M, Polliack A, Perry C, Herishanu Y. SLP76 integrates into the B-cell receptor signaling cascade in chronic lymphocytic leukemia cells and is associated with an aggressive disease course. Haematologica. 2016;101(12):1553–62.
The researchers acknowledge the contribution of Zhiyue Kang from Guangxi Medical University Cancer Hospital for helpful discussions and support.
National Natural Science Foundation of China (Grant No. 81660563); National Key R & D program (Grant No. 2017YFC0907103).
Ethics approval and consent to participate
Consent for publication
The authors declare no conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Chen, Q., Qiu, B., Zeng, X. et al. Identification of a tumor microenvironment-related gene signature to improve the prediction of cervical cancer prognosis. Cancer Cell Int 21, 182 (2021). https://doi.org/10.1186/s12935-021-01867-2
- Cervical cancer
- Tumor microenvironment
- Prognostic signature