Identification of an Immune Gene Signature Based on Tumor Microenvironment Characteristics in Colon Adenocarcinoma.
Ontology highlight
ABSTRACT:
SUBMITTER: Chen Y
PROVIDER: S-EPMC8020110 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:BackgroundColon adenocarcinoma (COAD) is a common gastrointestinal tumor and often occurs in the left colon with a poor prognosis. The progression of COAD is closely related to the tumor microenvironment, especially the hypoxia. Currently, few studies have reported the correlation between hypoxia-related genes and the prognosis of COAD patients. Furthermore, we constructed a prognostic model using four hypoxia-related genes to predict the prognosis of COAD patients.MethodsThe mRNA expression profiles and corresponding clinicopathological data were downloaded from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO). The string online analysis tool was used to construct a protein-protein interaction network (PPI) of hypoxia-related genes. Kaplan-Meier curve was used to analyze the relationship of hypoxia risk score and the overall survival of COAD patients, and the receiver operating characteristic (ROC) curve was used to assess the reliability.ResultsWe screened out four hypoxia genes, including TKTL1 (transketolase like 1), SLC2A3 (solute carrier family 2 member 3), ALDOB (aldolase, fructose-bisphosphate B) and ENO3 (enolase 3), which were used to construct a hypoxia risk model to predict the overall survival of COAD patients. Besides, we also found that the hypoxia risk score was correlated with the immunosuppression of tumor microenvironment.ConclusionThe model we constructed with four survival-related hypoxia genes, including TKTL1, SLC2A3, ALDOB and ENO3, could be used to predict the overall survival of COAD patients with high stability.
| S-EPMC8687688 | biostudies-literature
Project description:The composition and relative abundances of immune cells in the tumor microenvironment are key factors affecting the progression of lung adenocarcinomas (LUADs) and the efficacy of immunotherapy. Using the cancer gene expression dataset from The Cancer Genome Atlas (TCGA) program, we scored stromal and immune cells for tumor purity prediction by CIBERSORT and ESTMATE. Differential expression analysis was employed to identify 374 genes between the high-score group and the low-score group, which were utilized to conduct Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. Protein-protein interaction (PPI) and Cox regression analysis were performed on the differentially expressed genes (DEGs) to identify four key tumor microenvironment (TME) -related genes (CCR2, CCR4, P2RY12, and P2RY13). The expression levels of the four DEGs differed significantly among LUAD patients of different ages, genders, and TNM stages. We found that the infiltration of resting memory CD4+ T cells, memory B cells, and M0 macrophages into the TME was co-regulated by these four DEGs. These four genes were closely related to the prognosis of LUAD and affected the infiltration of immune cells into the TME, which had predictive prognostic value in LUAD.
| S-EPMC7880407 | biostudies-literature
Project description:Background Gastrointestinal malignant cancers affect many sites in the intestinal tract, including the colon. In this study, we purposed to improve prognostic predictions for colon cancer (CC) patients by establishing a novel biosignature of immune-related genes (IRGs) based on the tumor microenvironment (TME). Methods Using the estimation of stromal and immune cells in malignant tumor tissues using expression data (ESTIMATE) algorithm, we calculated the stromal and immune scores of every CC patient extracted from The Cancer Genome Atlas (TCGA). We then identified 4 immune-related messenger RNA (mRNA) biosignatures through a Cox and least absolute shrinkage and selection operator (LASSO) univariate analysis, and a Cox multivariate analysis. Relationships between tumor immune infiltration and the risk score were evaluated through the CIBERSORT algorithm and Tumor Immune Estimation Resource (TIMER) database. Results Our studies showed that individuals who had a high immune score (P=0.017) and low stromal score (P=0.041) had a favorable overall survival (OS) rate. By comparing high/low scores cohort, 220 differentially expressed genes (DEGs) were determined. Then an immune-related four-mRNA biosignature, including PDIA2, NAFTC1, VEGFC, and CD1B was identified. Kaplan-Meier, calibration, and receiver operating characteristic (ROC) curves verified the model’s performance. By using univariate and multivariate Cox analyses, we found each biosignature was an independent risk factor for assessing a CC patient’s survival. Three external GEO cohorts validated its good efficiency in estimating OS among individuals with CC. Moreover, the signature was also related to infiltration of several cells of the immune system in the tumor microenvironment. Conclusions The resultant model in our study included 4 IRGs associated with the TME. These IRGs can be utilized as an auxiliary variable to estimate and help improve the prognosis of individuals with CC.
| S-EPMC8576251 | biostudies-literature
Project description:Background: Lung cancer is the leading cause of cancer-related death globally. Hypoxia can suppress the activation of the tumor microenvironment (TME), which contributes to distant metastasis. However, the role of hypoxia-mediated TME in predicting the diagnosis and prognosis of lung adenocarcinoma (LUAD) patients remains unclear. Methods: Both RNA and clinical data from the LUAD cohort were downloaded from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. Both univariate and multivariate Cox regression analyses were used to further screen prognosis-related hypoxia gene clusters. Time-dependent receiver operation characteristic (ROC) curves were established to evaluate the predictive sensitivity and specificity of the hypoxia-related risk signature. The characterization of gene set enrichment analysis (GSEA) and TME immune cell infiltration were further explored to identify hypoxia-related immune infiltration. Results: Eight hypoxia-related genes (LDHA, DCN, PGK1, PFKP, FBP1, LOX, ENO3, and CXCR4) were identified and established to construct a hypoxia-related risk signature. The high-risk group showed a poor overall survival compared to that of the low-risk group in the TCGA and GSE68465 cohorts (p < 0.0001). The AUCs for 1-, 3-, and 5-year overall survival were 0.736 vs. 0.741, 0.656 vs. 0.737, and 0.628 vs. 0.649, respectively. The high-risk group was associated with immunosuppression in the TME. Conclusion: The hypoxia-related risk signature may represent an independent biomarker that can differentiate the characteristics of TME immune cell infiltration and predict the prognosis of LUAD.
| S-EPMC8634728 | biostudies-literature
Project description:The tumor microenvironment (TME) has been shown to be involved in angiogenesis, tumor metastasis, and immune response, thereby affecting the treatment and prognosis of patients. This study aims to identify genes that are dysregulated in the TME of patients with colon adenocarcinoma (COAD) and to evaluate their prognostic value based on RNA omics data. We obtained 512 COAD samples from the Cancer Genome Atlas (TCGA) database and 579 COAD patients from the independent dataset (GSE39582) in the Gene Expression Omnibus (GEO) database. The immune/stromal/ESTIMATE score of each patient based on their gene expression was calculated using the ESTIMATE algorithm. Kaplan-Meier survival analysis, Cox regression analysis, gene functional enrichment analysis, and protein-protein interaction (PPI) network analysis were performed. We found that immune and stromal scores were significantly correlated with COAD patients' overall survival (log rank p < 0.05). By comparing the high immune/stromal score group with the low score group, we identified 688 intersection differentially expressed genes (DEGs) from the TCGA dataset (663 upregulated and 25 downregulated). The functional enrichment analysis of intersection DEGs showed that they were mainly enriched in the immune process, cell migration, cell motility, Toll-like receptor signaling pathway, and PI3K-Akt signaling pathway. The hub genes were revealed by PPI network analysis. Through Kaplan-Meier and Cox analysis, four TME-related genes that were significantly related to the prognosis of COAD patients were verified in GSE39582. In addition, we uncovered the relationship between the four prognostic genes and immune cells in COAD. In conclusion, based on the RNA expression profiles of 1091 COAD patients, we screened four genes that can predict prognosis from the TME, which may serve as candidate prognostic biomarkers for COAD.
| S-EPMC8762242 | biostudies-literature
Project description:BackgroundThe 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).MethodsTranscriptome 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.ResultsGene 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.ConclusionsIn 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.
| S-EPMC8349485 | biostudies-literature
Project description:BackgroundLung cancer has become the most common cancer type and caused the most cancer deaths. Lung adenocarcinoma (LUAD) is one of the major types of lung cancer. Accumulating evidence suggests the tumor microenvironment is correlated with the tumor progress and the patient's outcome. This study aimed to establish a gene signature based on tumor microenvironment that can predict patients' outcomes for LUAD.MethodsDataset TCGA-LUAD, downloaded from the TCGA portal, were taken as training cohort, and dataset GSE72094, obtained from the GEO database, was set as validation cohort. In the training cohort, ESTIMATE algorithm was applied to find intersection differentially expressed genes (DEGs) among tumor microenvironment. Kaplan-Meier analysis and univariate Cox regression model were performed on intersection DEGs to preliminarily screen prognostic genes. Besides, the LASSO Cox regression model was implemented to build a multi-gene signature, which was then validated in the validation cohorts through Kaplan-Meier, Cox, and receiver operating characteristic curve (ROC) analyses. In addition, the correlation between tumor mutational burden (TMB) and risk score was evaluated by Spearman test. GSEA and immune infiltrating analyses were conducted for understanding function annotation and the role of the signature in the tumor microenvironment.ResultsAn eight-gene signature was built, and it was examined by Kaplan-Meier analysis, revealing that a significant overall survival difference was seen. The eight-gene signature was further proven to be independent of other clinico-pathologic parameters via the Cox regression analyses. Moreover, the ROC analysis demonstrated that this signature owned a better predictive power of LUAD prognosis. The eight-gene signature was correlated with TMB. Furthermore, GSEA and immune infiltrating analyses showed that the exact pathways related to the characteristics of eight-genes signature, and identified the vital roles of Mast cells resting and B cells naive in the prognosis of the eight-gene signature.ConclusionIdentifying the eight-gene signature (INSL4, SCN7A, STAP1, P2RX1, IKZF3, MS4A1, KLRB1, and ACSM5) could accurately identify patients' prognosis and had close interactions with Mast cells resting and B cells naive, which may provide insight into personalized prognosis prediction and new therapies for LUAD patients.
| S-EPMC7546815 | biostudies-literature
Project description:Background:Lung cancer has the highest morbidity and mortality worldwide, and lung adenocarcinoma (LADC) is the most common pathological subtype. Accumulating evidence suggests the tumor microenvironment (TME) is correlated with the tumor progress and the patient's outcome. As the major components of TME, the tumor-infiltrated immune cells and stromal cells have attracted more and more attention. In this study, differentially expressed immune and stromal signature genes were used to construct a TME-related prognostic model for predicting the outcomes of LADC patients. Methods:The expression profiles of LADC samples with clinical information were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO). The differentially expressed genes (DEGs) related to the TME of LADC were identified using TCGA dataset by Wilcoxon rank sum test. The prognostic effects of TME-related DEGs were analyzed using univariate Cox regression. Then, the least absolute shrinkage and selection operator (LASSO) regression was performed to reduce the overfit and the number of genes for further analysis. Next, the prognostic model was constructed by step multivariate Cox regression and risk score of each sample was calculated. Then, survival and Receiver Operating Characteristic (ROC) analyses were conducted to validate the model using TCGA and GEO datasets, respectively. The Kyoto Encyclopedia of Genes and Genomes analysis of gene signature was performed using Gene Set Enrichment Analysis (GSEA). Finally, the overall immune status, tumor purity and the expression profiles of HLA genes of high- and low-risk samples was further analyzed to reveal the potential mechanisms of prognostic effects of the model. Results:A total of 93 TME-related DEGs were identified, of which 23 DEGs were up-regulated and 70 DEGs were down-regulated. The univariate cox analysis indicated that 23 DEGs has the prognostic effects, the hazard ratio ranged from 0.65 to 1.25 (p < 0.05). Then, seven genes were screened out from the 23 DEGs by LASSO regression method and were further analyzed by step multivariate Cox regression. Finally, a three-gene (ADAM12, Bruton Tyrosine Kinase (BTK), ERG) signature was constructed, and ADAM12, BTK can be used as independent prognostic factors. The three-gene signature well stratified the LADC patients in both training (TCGA) and testing (GEO) datasets as high-risk and low-risk groups, the 3-year area under curve (AUC) of ROC curves of three GEO sets were 0.718 (GSE3141), 0.646 (GSE30219) and 0.643 (GSE50081). The GSEA analysis indicated that highly expressed ADAM12, BTK, ERG mainly correlated with the activation of pathways involving in focal adhesion, immune regulation. The immune analysis indicated that the low-risk group has more immune activities and higher expression of HLA genes than that of the high-risk group. In sum, we identified and constructed a three TME-related DEGs signature, which could be used to predict the prognosis of LADC patients.
| S-EPMC6886493 | biostudies-literature
Project description:Lung adenocarcinoma (LUAD) is a common subtype of lung cancer with a depressing survival rate. The reprogramming of tumor metabolism was identified as a new hallmark of cancer in tumor microenvironment (TME), and we made a comprehensive exploration to reveal the prognostic role of the metabolic-related genes. Transcriptome profiling data of LUAD were, respectively, downloaded from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database. Based on the extracted metabolic-related genes, a novel 5-gene metabolic prognostic signature (including GNPNAT1, LPGAT1, TYMS, LDHA, and PTGES) was constructed by univariate Cox regression and least absolute shrinkage and selection operator (LASSO) regression. This signature confirmed its robustness and accuracy by external validation in multiple databases. It could be an independent risk factor for LUAD, and the nomograms possessed moderately accurate performance with the C-index of 0.755 (95% confidence interval: 0.706-0.804) and 0.691 (95% confidence interval: 0.636-0.746) in training set and testing set. This signature could reveal the metabolic features according to the results of gene set enrichment analysis (GSEA) and meanwhile monitor the status of TME through ESTIMATE scores and the infiltration levels of immune cells. In conclusion, this gene signature is a cost-effective tool which could indicate the status of TME to provide more clues in the exploration of new diagnostic and therapeutic strategy.
| S-EPMC7335383 | biostudies-literature
Project description:BackgroundThe tumor microenvironment (TME) is associated with disease outcomes and treatment response in colon cancer. Here, we constructed a TME-related gene signature that is prognosis of disease survival and may predict response to immunotherapy in colon cancer.MethodsWe calculated immune and stromal scores for 385 colon cancer samples from The Cancer Genome Atlas (TCGA) database using the ESTIMATE algorithm. We identified nine TME-related prognostic genes using Cox regression analysis. We evaluated associations between protein expression, extent of immune cell infiltrate, and patient survival. We calculated risk scores and built a clinical predictive model for the TME-related gene signature. Receiver operating characteristic (ROC) curves were generated to assess the predictive power of the signature. We estimated the half-maximal inhibitory concentration (IC50) of chemotherapeutic drugs in patients using the pRRophetic algorithm. The expression of immune checkpoint genes was evaluated.ResultsHigh immune and stromal scores are significantly associated with poor overall survival (p < 0.05). We identified 773 differential TME-related prognostic genes associated with survival; these genes were enriched in immune-related pathways. Nine key prognostic genes were identified and were used to construct a TME-related prognostic signature: CADM3, LEP, CD1B, PDE1B, CCL22, ABI3BP, IGLON5, SELE, and TGFB1. This signature identified a high-risk group with worse survival outcomes, based on Kaplan-Meier analysis. A nomogram composed of clinicopathological factors and risk score exhibited good accuracy. Drug sensitivity analysis identified no difference in sensitivity between the high-risk and low-risk groups. High-risk patients had higher expression of PD-1, PDL-1, and CTLA-4 and lower expression of LAG-3 and VSIR. Infiltration of dendritic cells was higher in the high-risk group.ConclusionsWe identified a novel prognostic TME-related gene expression signature in colon cancer. Stratification of patients based on this gene signature could be used to improve outcomes and guide better therapy for colon cancer patients.
| S-EPMC8420973 | biostudies-literature