Dataset Information

Identification of prognostic immune-related gene signature associated with tumor microenvironment of colorectal cancer.

ABSTRACT:

Background

The tumor microenvironment (TME) has significantly correlation with tumor occurrence and prognosis. Our study aimed to identify the prognostic immune-related genes (IRGs)in the tumor microenvironment of colorectal cancer (CRC).

Methods

Transcriptome and clinical data of CRC cases were downloaded from TCGA and GEO databases. Stromal score, immune score, and tumor purity were calculated by the ESTIMATE algorithm. Based on the scores, we divided CRC patients from the TCGA database into low and high groups, and the differentially expressed genes (DEGs) were identified. Immune-related genes (IRGs) were selected by venn plots. To explore underlying pathways, protein-protein interaction (PPI) networks and functional enrichment analysis were used. After utilizing LASSO Cox regression analysis, we finally established a multi-IRGs signature for predicting the prognosis of CRC patients. A nomogram consists of the thirteen-IRGs signature and clinical parameters was developed to predict the overall survival (OS). We investigated the association between prognostic validated IRGs and immune infiltrates by TIMER database.

Results

Gene expression profiles and clinical information of 1635 CRC patients were collected from the TCGA and GEO databases. Higher stromal score, immune score and lower tumor purity were observed positive correlation with tumor stage and poor OS. Based on stromal score, immune score and tumor purity, 1517 DEGs, 1296 DEGs, and 1892 DEGs were identified respectively. The 948 IRGs were screened by venn plots. A thirteen-IRGs signature was constructed for predicting survival of CRC patients. Nomogram with a C-index of 0.769 (95%CI, 0.717-0.821) was developed to predict survival of CRC patients by integrating clinical parameters and thirteen-IRGs signature. The AUC for 1-, 3-, and 5-year OS were 0.789, 0.783 and 0.790, respectively. Results from TIMER database revealed that CD1B, GPX3 and IDO1 were significantly related with immune infiltrates.

Conclusions

In this study, we established a novel thirteen immune-related genes signature that may serve as a validated prognostic predictor for CRC patients, thus will be conducive to individualized treatment decisions.

SUBMITTER: Wang Y

PROVIDER: S-EPMC8349485 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Neuroblastoma (NB) is the most prevalent malignant solid tumor in children. Tumor metabolism, including lipid, amino acid, and glucose metabolism, is intricately linked to the genesis and progression of tumors. This study aimed to establish a prognostic gene signature for NB patients, based on metabolism-related genes, and to investigate a treatment approach that could enhance the survival rate of high-risk NB patients. From the NB dataset GSE49710, we identified metabolism-related gene markers utilizing the "limma" R package and univariate Cox analysis combined with least absolute shrinkage and selection operator (LASSO) regression analysis. We explored the correlation between these gene markers and the overall survival of NB patients. Gene set enrichment analysis (GSEA) and single-sample GSEA algorithms were used to assess the differences in metabolism and immune status. Furthermore, we examined the association between metabolic subgroups and drug responsiveness. Concurrently, data downloaded from TARGET and MTAB were used for external verification. Using multicolor immunofluorescence and immunohistochemistry, we investigated the relationship between the lipid metabolism-related gene ELOVL6 with both the International Neuroblastoma Staging System classification of NB and survival rate. Finally, we explored the effect of high ELOVL6 expression on the immune microenvironment in NB using flow cytometry. We identified an eight-gene signature comprising metabolism-related genes in NB: ELOVL6, OSBPL9, RPL27A, HSD17B3, ACHE, AKR1C1, PIK3R1, and EPHX2. This panel effectively predicted disease-free survival, and was validated using an internal dataset from GSE49710 and two external datasets from the TARGET and MTAB databases. Moreover, our findings confirmed that ELOVL6 fosters an immunosuppressive microenvironment and contributes to the malignant progression in NB. The eight-gene signature is significant in predicting the prognosis of NB, effectively classifying patients into high- and low-risk groups. This classification may guide the development of innovative treatment strategies for these patients. Notably, the signature gene ELOVL6 markedly encourages an immunosuppressive microenvironment and malignant progression in NB.

Project description: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.

Project description:BackgroundThe current study probed prognosis-related potential for m6A-related lncRNAs signatures within colon tumor immune microenvironment (TIM).MethodsAfter downloading transcriptomic datasets for colon cancer (CC) patients from The Cancer Genome Atlas (TCGA), they were divided, in a 1:1 ratio, within training or test datasets. m6A-related lncRNAs were then scrutinized across such dataset using Pearson correlation assessment before generating a m6A-related lncRNAs prognosis-related model using the training dataset. The latter was then validated with the test and the whole dataset. In addition, we compared the differences of TIM and the estimated IC50 of drug response between the high- and low-risk groups.ResultsOverall survival (OS) resulted as linked with 11 m6A-related lncRNAs, while within the developed prognosis-related model, areas-under-curves were as follows: within training dataset, values at 3-, 4-, and 5-years were 0.777, 0.819, and 0.805, accordingly, and for test one, they were 0.697, 0.682, and 0.706, respectively. Finally, the values for the whole dataset were 0.675 (3-year), 0.682 (4-years), and 0.679 (5-years), accordingly. Moreover, CC cases categorized within low-risk cohort demonstrated enhanced OS (p < .0001), lower metastasis (p = 2e-06) and lower T stage (p = .0067), more instability for microsatellite status (p = .012), and downregulation for PD-L1, PD-1, CTLA-4, LAG3, and HAVCR2 (p < .05). In addition, risk scorings were significantly linked to the degree of infiltrative intensity for CD8 and CD4 (memory resting) T-cells, T-regulatory (Tregs), and Mast cells triggering (p < .05). Patients with low infiltrative propensity for CD4 T-cells also had better OS (p = .016). Moreover, six representative drugs were found to be sensitive for treating CC patients.ConclusionA robust m6A-related prognostic model with great performances was developed before exploring the TIM characteristics and its potential therapeutic drugs, which might improve the prognosis and therapeutic efficacy.

Project description:ObjectiveStromal cells and immune cells have important clinical significance in the microenvironment of colorectal cancer (CRC). This study is aimed at developing a CRC gene signature on the basis of stromal and immune scores.MethodsA cohort of CRC patients (n = 433) were adopted from The Cancer Genome Atlas (TCGA) database. Stromal/immune scores were calculated by the ESTIMATE algorithm. Correlation between prognosis/clinical characteristics and stromal/immune scores was assessed. Differentially expressed stromal and immune genes were identified. Their potential functions were annotated by functional enrichment analysis. Cox regression analysis was used to develop an eight-gene risk score model. Its predictive efficacies for 3 years, 5 years, overall survival (OS), and progression-free survival interval (PFI) were evaluated using time-dependent receiver operating characteristic (ROC) curves. The correlation between the risk score and the infiltering levels of six immune cells was analyzed using TIMER. The risk score was validated using an independent dataset.ResultsImmune score was in a significant association with prognosis and clinical characteristics of CRC. 736 upregulated and two downregulated stromal and immune genes were identified, which were mainly enriched into immune-related biological processes and pathways. An-eight gene prognostic risk score model was conducted, consisting of CCL22, CD36, CPA3, CPT1C, KCNE4, NFATC1, RASGRP2, and SLC2A3. High risk score indicated a poor prognosis of patients. The area under the ROC curves (AUC) s of the model for 3 years, 5 years, OS, and PFI were 0.71, 0.70, 0.73, and 0.66, respectively. Thus, the model possessed well performance for prediction of patients' prognosis, which was confirmed by an external dataset. Moreover, the risk score was significantly correlated with immune cell infiltration.ConclusionOur study conducted an immune-related prognostic risk score model, which could provide novel targets for immunotherapy of CRC.

Project description:Purpose: The purpose of this study was to construct a novel risk scoring model with prognostic value that could elucidate tumor immune microenvironment of hepatocellular carcinoma (HCC). Samples and methods: Data were obtained through The Cancer Genome Atlas (TCGA) database. Univariate Cox analysis, least absolute shrinkage and selection operator (LASSO) analysis, and multivariate Cox analysis were carried out to screen for glycolysis-related long noncoding RNAs (lncRNAs) that could provide prognostic value. Finally, we established a risk score model to describe the characteristics of the model and verify its prediction accuracy. The receiver operating characteristic (ROC) curves of 1, 3, and 5 years of overall survival (OS) were depicted with risk score and some clinical features. ESTIMATE algorithm, single-sample gene set enrichment analysis (ssGSEA), and CIBERSORT analysis were employed to reveal the characteristics of tumor immune microenvironment in HCC. The nomogram was drawn by screening indicators with high prognostic accuracy. The correlation of risk signature with immune infiltration and immune checkpoint blockade (ICB) therapy was analyzed. After enrichment of related genes, active behaviors and pathways in high-risk groups were identified and lncRNAs related to poor prognosis were validated in vitro. Finally, the impact of MIR4435-2HG upon ICB treatment was uncovered. Results: After screening through multiple steps, four glycolysis-related lncRNAs were obtained. The risk score constructed with the four lncRNAs was found to significantly correlate with prognosis of samples. From the ROC curve of samples with 1, 3, and 5 years of OS, two indicators were identified with high prognostic accuracy and were used to draw a nomogram. Besides, the risk score significantly correlated with immune score, immune-related signature, infiltrating immune cells (i.e. B cells, etc.), and ICB key molecules (i.e. CTLA4,etc.). Gene enrichment analysis indicated that multiple biological behaviors and pathways were active in the high-risk group. In vitro validation results showed that MIR4435-2HG was highly expressed in the two cell lines, which had a significant impact on the OS of samples. Finally, we corroborated that MIR4435-2HG had intimate relationship with ICB therapy in hepatocellular carcinoma. Conclusion: We elucidated the crucial role of risk signature in immune cell infiltration and immunotherapy, which might contribute to clinical strategies and clinical outcome prediction of HCC.

Project description:BackgroundEsophageal carcinoma (ESCA) is a common malignancy with a poor prognosis. Previous research has suggested that necroptosis is involved in anti-tumor immunity and promotes oncogenesis and cancer metastasis, which in turn affects tumor prognosis. However, the role of necroptosis in ESCA is unclear. This study aimed to investigate the relationships between necroptosis-related genes (NRGs) and ESCA.Methods and resultsThe clinical data and gene expression profiles of ESCA patients were extracted from The Cancer Genome Atlas (TCGA), and 159 NRGs were screened from the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. We then identified 52 differentially expressed NRGs associated with ESCA and used them for further analysis. Gene ontology (GO) and KEGG functional enrichment analyses showed that these NRGs were mostly associated with the regulation of necroptosis, Influenza A, apoptosis, NOD-like receptor, and NF-Kappa B signaling pathway. Next, univariate and multivariate Cox regression and LASSO analysis were used to identify the correlation between NRGs and the prognosis of ESCA. We constructed a prognostic model to predict the prognosis of ESCA based on SLC25A5, PPIA, and TNFRSF10B; the model classified patients into high- and low-risk subgroups based on the patient's risk score. Furthermore, the receiver operating characteristic (ROC) curve was plotted, and the model was affirmed to perform moderately well for prognostic predictions. In addition, Gene Expression Omnibus (GEO) datasets were selected to validate the applicability and prognostic value of our predictive model. Based on different clinical variables, we compared the risk scores between the subgroups of different clinical features. We also analyzed the predictive value of this model for drug sensitivity. Moreover, Immunohistochemical (IHC) validation experiments explored that these three NRGs were expressed significantly higher in ESCA tissues than in adjacent non-tumor tissues. In addition, a significant correlation was observed between the three NRGs and immune-cell infiltration and immune checkpoints in ESCA.ConclusionsIn summary, we successfully constructed and validated a novel necroptosis-related signature containing three genes (SLC25A5, PPIA, and TNFRSF10B) for predicting prognosis in patients with ESCA; these three genes might also play a crucial role in the progression and immune microenvironment of ESCA.

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.

Project description:Colorectal cancer (CRC) is a heterogeneous disease and one of the most prevalent malignancies worldwide. Previous research has demonstrated that mitophagy is crucial to developing colorectal cancer. This study aims to examine the association between mitophagy-related genes and the prognosis of CRC patients. Gene expression profiles and clinical information of CRC patients were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. Univariate Cox regression and the least absolute shrinkage and selection operator (LASSO) regression analysis were applied to establish a prognostic signature using mitophagy related genes. Kaplan-Meier and receiver operating characteristic (ROC) curves were used to analyze patient survival and predictive accuracy. Meanwhile, we also used the Genomics of Drug Sensitivity in Cancer (GDSC) database and Tumor Immune Dysfunction and Exclusion (TIDE) algorithm to estimate the sensitivity of chemotherapy, targeted therapy and immunotherapy. ATG14 overexpression plasmid was used to regulate the ATG14 expression level in HCT116 and SW480 cell lines, and cell counting kit-8, colony formation and transwell migration assay were performed to validate the function of ATG14 in CRC cells. A total of 22 mitophagy-driven genes connected with CRC survival were identified, and then a novel prognostic signature was established based on 10 of them (AMBRA1, ATG14, MAP1LC3A, MAP1LC3B, OPTN, VDAC1, ATG5, CSNK2A2, MFN1, TOMM22). Patients were divided into high-risk and low-risk groups based on the median risk score, and the survival of patients in the high-risk group was significantly shorter in both the training cohort and two independent cohorts. ROC curve showed that the area under the curves (AUC) of 1-, 3- and 5-year survival were 0.66, 0.66 and 0.64, respectively. Multivariate Cox regression analysis confirmed the independent prognostic value of the signature. Then we constructed a Nomogram combining the risk score, age and M stage, which had a concordance index of survival prediction of 0.77 (95% CI 0.71-0.83) and more robust predictive accuracy. Results showed that CD8+ T cells, regulatory T cells and activated NK cells were significantly more enriched in the high-risk group. Furthermore, patients in the high-risk group are more sensitive to targeted therapy or chemotherapy, including bosutinib, elesclomol, lenalidomide, midostaurin, pazopanib and sunitinib, while the low-risk group is more likely to benefit from immunotherapy. Finally, in vitro study confirmed the oncogenic significance of ATG14 in both HCT116 and SW480 cells, whose overexpression increased CRC cell proliferation, colony formation, and migration. In conclusion, we developed a novel mitophagy-related gene signature that can be utilized not only as an independent predictive biomarker but also as a tool for tailoring personalizing treatment for CRC patients, and we confirmed ATG14 as a novel oncogene in CRC.