Dataset Information

Comprehensive analysis of prognostic immune-related genes and drug sensitivity in cervical cancer.

ABSTRACT:

Background

Cervical cancer (CC) is the leading cause of cancer-related death in women. A limited number of studies have investigated whether immune-prognostic features can be used to predict the prognosis of CC. This study aimed to develop an improved prognostic risk scoring model (PRSM) for CC based on immune-related genes (IRGs) to predict survival and determine the key prognostic IRGs.

Methods

We downloaded the gene expression profiles and clinical data of CC patients from the TCGA and GEO databases. The ESTIMATE algorithm was used to calculate the score for both immune and stromal cells. Differentially expressed genes (DEGs) in different subpopulations were analyzed by "Limma". A weighted gene co-expression network analysis (WGCNA) was used to establish a DEG co-expression module related to the immune score. Immune-related gene pairs (IRGPs) were constructed, and univariate- and Lasso-Cox regression analyses were used to analyze prognosis and establish a PRSM. A log-rank test was used to verify the accuracy and consistency of the scoring model. Identification of the predicted key IRG was ensured by the application of functional enrichment, DisNor, protein-protein interactions (PPIs) and heatmap. Finally, we extracted the key prognostic immune-related genes from the gene expression data, validated the key genes by immunohistochemistry and analyzed the correlation between their expression and drug sensitivity.

Results

A new PRSM was developed based on 22 IRGPs. The prognosis of the low-risk group in the model group (P < 0.001) and validation group (P = 0.039) was significantly better than that in the high-risk group. Furthermore, M1 and M2 macrophages were highly expressed in the low-risk group. Retinoic acid-inducible gene-I (RIG-I)-like receptors (RLRs) and the Janus kinase-signal transducer and activator of transcription (JAK-STAT) signaling pathway were significantly enriched in the low-risk group. Three representative genes (CD80, CD28, and LCP2) were markers of CC prognosis. CD80 and CD28 may more prominent represent important indicators to improve patient prognosis. These key genes was positively correlated with drug sensitivity. Finally, we found that differences in the sensitivity to JNK inhibitors could be distinguished based on the use and risk grouping of this PRSM.

Conclusions

The prognostic model based on the IRGs and key genes have potential clinical significance for predicting the prognosis of CC patients, providing a foundation for clinical prognosis judgment and individualized treatment.

SUBMITTER: Pi YN

PROVIDER: S-EPMC8638517 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Cervical cancer and endometrial cancer remain serious threats to women's health. Even though some patients can be treated with surgery plus chemoradiotherapy as a conventional option, the overall efficacy is deemed unsatisfactory. As such, the development for new treatment approaches is truly necessary. In recent years, immunotherapy has been widely used in clinical practice and it is an area of great interest that researchers are keeping attention on. However, a thorough immune-related genes (IRGs) study for cervical cancer and endometrial cancer is still lacking. We therefore aim to make a comprehensive evaluation of IRGs through bioinformatics and large databases, and also investigate the relationship between the two types of cancer. We reviewed the transcriptome RNAs of IRGs and clinical data based on the TCGA database. Survival-associated IRGs in cervical/endometrial cancer were identified using univariable and multivariable Cox proportional-hazard regression analysis for developing an IRG signature model to evaluate the risk of patients. In the end, this model was validated based on the enrichment analyses through GO, KEGG, and GSEA pathways, Kaplan-Meier survival curve, ROC curves, and immune cell infiltration. Our results showed that out of 25/23 survival-associated IRGs for cervical/endometrial cancer, 13/12 warranted further examination by multivariate Cox proportional-hazard regression analysis and were selected to develop an IRGs signature model. As a result, enrichment analyses for high-risk groups indicated main enriched pathways were associated with tumor development and progression, and statistical differences were found between high-risk and low-risk groups as shown by Kaplan-Meier survival curve. This model could be used as an independent measure for risk assessment and was considered relevant to immune cell infiltration, but it had nothing to do with clinicopathological characteristics. In summary, based on comprehensive analysis, we obtained the IRGs signature model in cervical cancer (LTA, TFRC, TYK2, DLL4, CSK, JUND, NFATC4, SBDS, FLT1, IL17RD, IL3RA, SDC1, PLAU) and endometrial cancer (LTA, PSMC4, KAL1, TNF, SBDS, HDGF, LTB, HTR3E, NR2F1, NR3C1, PGR, CBLC), which can effectively evaluate the prognosis and risk of patients and provide justification in immunology for further researches.

Project description:AbstractGrowing evidence supports that the tumor microenvironment plays a key role in the development and progression of tumors. But immune microenvironment of hepatocellular carcinoma (HCC) has not yet been fully explored. In the present investigation, the clinical value and prognostic significance of immune-related genes in HCC were investigated.The immune and stromal scores of HCC were calculated through the application of Estimation of Stromal and Immune cells in Malignant Tumor tissues using Expression data Algorithm based on the Cancer Genome Atlas database. Differentially expressed genes were identified using the "edgeR" package of the R software. Functional annotation and pathway enrichment were performed using "ggplots2" and "clusterProfiler" packages in R software. Protein-protein interaction network was constructed using STRING, and the hub genes were identified through the Cytoscape. Survival analysis was performed using Kaplan-Meier methods. Tumor Immune Estimation Resource algorithm was used to view the immune landscape of the microenvironment in HCC.Firstly, the immune and stromal scores of HCC were calculated and we found that the immune and stromal scores of HCC were closely related to the patients' prognosis. Then the differentially expressed genes were identified respectively stratified by the median value of the immune and stromal scores, and the immune-related genes that related to the prognosis in HCC patients were further identified. Functional enrichment analysis and protein-protein interaction networks further showed that these genes mainly participated in immune-related biological process. In addition, dendritic cells were found to be the most abundant in the microenvironment of HCC through Tumor Immune Estimation Resource algorithm and were significantly associated with the patients' prognosis. To robust the results, the immune-related genes were validated in an independent dataset from the Gene Expression Omnibus database.We arrived at a more comprehensive understanding of the microenvironment of HCC and extracted 7 immune-related genes that were significantly associated with the recurrence survival of HCC.

Project description:Pancreatic cancer remains to have a high mortality, which is partly due to the lack of effective treatment strategies. In this study, genes with potential associations with immunophenotyping of pancreatic cancer were screened through bioinformatics analysis and the correlation between immune-related genes and the prognosis of pancreatic cancer patients was assessed. Firstly, differentially expressed immune genes were extracted from the pancreatic cancer-related datasets obtained for purposes of this study. The samples were processed by the "Consensus Cluster Plus" R package to determine the number of immune subtypes. Then, the pancreatic cancer immunophenotyping-related gene modules were determined. Differential analysis of immune gene modules was performed, and the function of genes related to pancreatic cancer immune subtypes was identified. The number of immune cells in the samples was calculated, followed by the differential expression analysis of immune cell numbers in each immune subtype of pancreatic cancer. The immune infiltration score was also estimated, and the correlation between the immune infiltration score and the patient prognosis with different immune subtypes was determined. Gene differences between each immune subtype were identified by differential expression analysis, and key immune genes affecting immunophenotyping were obtained. Following the analysis, 426 immune-related genes were identified to have potential involvement in the occurrence and development of pancreatic cancer, of which CD19 may be the most critical gene affecting the immunophenotyping of pancreatic cancer. CD19 played a significant role in the occurrence and development of IS2 and IS3 immune subtypes of pancreatic cancer through its action on B cells and T cells. Moreover, the expression of CD19 was increased in the collected pancreatic cancer tissues. Overall, our findings uncovered the critical role of CD19 in the prognosis of pancreatic cancer patients.

Project description:It is reported that ferroptosis has close relation with tumorigenesis and drug resistance. However, the clinical significance of ferroptosis in lung adenocarcinoma (LUAD) remains elusive, and the potential targets for ferroptosis-based treatment are limited. In this study, we constructed a 15-gene prognostic signature predicting overall survival based on ferroptosis-related genes (ferroptosis driver genes VDAC2, GLS2, FLT3, TLR4, PHKG2, phosphogluconate dehydrogenase (PGD), PANX1, KRAS, PEBP1, ALOX15, and ALOX12B, and suppressor genes ACSL3, CISD1, FANCD2, and SLC3A2) in The Cancer Genome Atlas (TCGA)-LUAD cohort. The signature's predictive ability was validated in the GSE68465 and GSE72094 cohorts by survival analysis and independent prognostic analysis with clinical features. Nomograms were provided for clinical reference. Functional analysis revealed that ferroptosis was closely related to cell cycle, cell metabolism, and immune pathways. Pan-cancer analysis comprehensively analyzed these 15 genes in 33 cancer types, indicating that the heterogeneity of 15 genes was evident across different cancer types. Besides, these genes were critical regulators modulating drug resistance, tumor microenvironment infiltration, and cancer stemness. Then, we screened 10 genes (TLR4, PHKG2, PEBP1, GLS2, FLT3, ALOX15, ACSL3, CISD1, FANCD2, and SLC3A2) as potential targets for further research because their biological functions in ferroptosis were consistent with their prognostic significance. Somatic mutation and copy number variation analysis revealed that the alteration rates of KRAS, PGD, and ALOX15 were more than 1% and significantly associated with overall survival in LUAD. Moreover, the expression of KRAS and PGD was positively related to tumor mutation burden, indicating that KRAS and PGD could serve as novel biomarkers for predicting immunotherapy response rate. Our study identified and validated a ferroptosis-related gene signature for LUAD, provided a 10-gene set for future research, and screened KRAS and PGD as potential novel immunotherapy biomarkers.

Project description:BackgroundAccumulating evidence demonstrated that tumor microenvironmental cells played important roles in predicting clinical outcomes and therapeutic efficacy. We aimed to develop a reliable immune-related gene signature for predicting the prognosis of ovarian cancer (OC).MethodsSingle sample gene-set enrichment analysis (ssGSEA) of immune gene-sets was used to quantify the relative abundance of immune cell infiltration and develop high- and low-abundance immune subtypes of 308 OC samples. The presence of infiltrating stromal/immune cells in OC tissues was calculated as an estimate score. We estimated the correlation coefficients among the immune subtype, clinicopathological feature, immune score, distribution of immune cells, and tumor mutation burden (TMB). The differentially expressed immune-related genes between high- and low-abundance immune subtypes were further used to construct a gene signature of a prognostic model in OC with lasso regression analysis.ResultsThe ssGSEA analysis divided OC samples into high- and low-abundance immune subtypes based on the abundance of immune cell infiltration, which was significantly related to the estimate score and clinical characteristics. The distribution of immune cells was also significantly different between high- and low-abundance immune subtypes. The correlation analysis showed the close relationship between TMB and the estimate score. The differentially expressed immune-related genes between high- and low-abundance immune subtypes were enriched in multiple immune-related pathways. Some immune checkpoints (PDL1, PD1, and CTLA-4) were overexpressed in the high-abundance immune subtype. Furthermore, the five-immune-related-gene-signature prognostic model (CCL18, CXCL13, HLA-DOB, HLA-DPB2, and TNFRSF17)-based high-risk and low-risk groups were significantly related to OC overall survival.ConclusionImmune-related genes were the promising predictors of prognosis and survival, and the comprehensive landscape of tumor microenvironmental cells of OC has potential for therapeutic schedule monitoring.