Identification of tumor microenvironment-related prognostic genes in colorectal cancer based on bioinformatic methods.
Ontology highlight
ABSTRACT:
SUBMITTER: Liu Y
PROVIDER: S-EPMC8298640 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:BackgroundTo identify prognostic genes which were associated with adrenocortical carcinoma (ACC) tumor microenvironment (TME).Methods and materialsTranscriptome profiles and clinical data of ACC samples were collected from The Cancer Genome Atlas (TCGA) database. We use ESTIMATE (estimation of stromal and Immune cells in malignant tumor tissues using expression data) algorithm to calculate immune scores, stromal scores and estimate scores. Heatmap and volcano plots were applied for differential analysis. Venn plots were used for intersect genes selection. We used protein-protein interaction (PPI) networks and functional analysis to explore underlying pathways. After performing stepwise regression method and multivariate Cox analysis, we finally screened hub genes associated with ACC TME. We calculated risk scores (RS) for ACC cases based on multivariate Cox results and evaluated the prognostic value of RS shown by receiver operating characteristic curve (ROC). We investigated the association between hub genes with immune infiltrates supported by algorithm from online TIMER database.ResultsGene expression profiles and clinical data were downloaded from TCGA. Lower immune scores were observed in disease with distant metastasis (DM) and locoregional recurrence (LR) than other cases (P = .0204). Kaplan-Meier analysis revealed that lower immune scores were significantly associated with poor overall survival (OS) (P = .0495). We screened 1649 differentially expressed genes (DEGs) and 1521 DEGs based on immune scores and stromal scores, respectively. Venn plots helped us find 1122 intersect genes. After analysing by cytoHubba from Cytoscape software, 18 hub genes were found. We calculated RS and ROC showed significantly predictive accuracy (area under curve (AUC) = 0.887). ACC patients with higher RS had worse survival outcomes (P < .0001). Results from TIMER (tumor immune estimation resource) database revealed that HLA-DOA was significantly related with immune cells infiltration.ConclusionWe screened a list of TME-related genes which predict poor survival outcomes in ACC patients from TCGA database.
| S-EPMC6997077 | biostudies-literature
Project description:ObjectiveGrowing evidence has highlighted that the immune and stromal cells that infiltrate in pancreatic cancer microenvironment significantly influence tumor progression. However, reliable microenvironment-related prognostic gene signatures are yet to be established. The present study aimed to elucidate tumor microenvironment-related prognostic genes in pancreatic cancer.MethodsWe applied the ESTIMATE algorithm to categorize patients with pancreatic cancer from TCGA dataset into high and low immune/stromal score groups and determined their differentially expressed genes. Then, univariate and LASSO Cox regression was performed to identify overall survival-related differentially expressed genes (DEGs). And multivariate Cox regression analysis was used to screen independent prognostic genes and construct a risk score model. Finally, the performance of the risk score model was evaluated by Kaplan-Meier curve, time-dependent receiver operating characteristic and Harrell's concordance index.ResultsThe overall survival analysis demonstrated that high immune/stromal score groups were closely associated with poor prognosis. The multivariate Cox regression analysis indicated that the signatures of four genes, including TRPC7, CXCL10, CUX2, and COL2A1, were independent prognostic factors. Subsequently, the risk prediction model constructed by those genes was superior to AJCC staging as evaluated by time-dependent receiver operating characteristic and Harrell's concordance index, and both KRAS and TP53 mutations were closely associated with high risk scores. In addition, CXCL10 was predominantly expressed by tumor associated macrophages and its receptor CXCR3 was highly expressed in T cells at the single-cell level.ConclusionsThis study comprehensively investigated the tumor microenvironment and verified immune/stromal-related biomarkers for pancreatic cancer.
| S-EPMC8277941 | biostudies-literature
Project description:AimImmune cells that infiltrate the tumor microenvironment (TME) are associated with cancer prognosis. The aim of the current study was to identify TME related gene signatures related to the prognosis of sarcoma (SARC) by using the data from The Cancer Genome Atlas (TCGA).MethodsImmune and stromal scores were calculated by estimation of stromal and immune cells in malignant tumor tissues using expression data algorithms. The least absolute shrinkage and selection operator (lasso) based cox model was then used to select hub survival genes. A risk score model and nomogram were used to predict the overall survival of patients with SARC.ResultsWe selected 255 patients with SARC for our analysis. The Kaplan-Meier method found that higher immune (p = 0.0018) or stromal scores (p = 0.0022) were associated with better prognosis of SARC. The estimated levels of CD4+ (p = 0.0012) and CD8+ T cells (p = 0.017) via the tumor immune estimation resource were higher in patients with SARC with better overall survival. We identified 393 upregulated genes and 108 downregulated genes (p < 0.05, fold change >4) intersecting between the immune and stromal scores based on differentially expressed gene (DEG) analysis. The univariate Cox analysis of each intersecting DEG and subsequent lasso-based Cox model identified 11 hub survival genes (MYOC, NNAT, MEDAG, TNFSF14, MYH11, NRXN1, P2RY13, CXCR3, IGLV3-25, IGHV1-46, and IGLV2-8). Then, a hub survival gene-based risk score gene signature was constructed; higher risk scores predicted worse SARC prognosis (p < 0.0001). A nomogram including the risk scores, immune/stromal scores and clinical factors showed a good prediction value for SARC overall survival (C-index = 0.716). Finally, connectivity mapping analysis identified that the histone deacetylase inhibitors trichostatin A and vorinostat might have the potential to reverse the harmful TME for patients with SARC.ConclusionThe current study provided new indications for the association between the TME and SARC. Lists of TME related survival genes and potential therapeutic drugs were identified for SARC.
| S-EPMC7882740 | 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:Background and Objective: Esophageal cancer (ESCA) is a commonly occurring cancer worldwide with poor survival and limited therapeutic options. Due to the lack of biomarkers that facilitate early detection, its treatment remains a great challenge. This study aims at identifying the tumor microenvironment (TME)-related genes, which might affect prognosis and accelerate clinical treatment for ESCA patients. Methods: We integrated the expression profiles from ESCA patients in The Cancer Genome Atlas. Then, we determined the stromal and immune scores of each sample using the R package. The Gene Expression Omnibus database was used to validate the expression profile of the key genes. Results: Tumor mutational burden showed a significant difference between the groups of ESCA patients with high and low ESTIMATE scores. We identified 859 intersection genes among patients with different immune and stromal scores. Moreover, gene ontology analysis demonstrated that these 859 intersection genes were closely related to adaptive immune response and regulation of lymphocyte activation. Kyoto Encyclopedia of Genes and Genomes showed the enrichment of cytokine-cytokine receptor interaction and chemokine signaling pathway in the TME. Furthermore, the protein-protein interaction network consisted of 175 nodes. We selected 35 hub genes, including ITGAM, CXCL10, CCR2, CCR5, and CCR1. Of these, 23 intersection genes predicted the overall survival rate. C1QA and FCER1G correlated with overall survival of the ESCA patients in the two databases. Conclusion: We identified a set of stromal and immune score-related prognostic differentially expressed genes that could influence the complexity of the TME. C1QA and FCER1G were identified and validated with respect to their role in the progression of ESCA.
| S-EPMC8262216 | biostudies-literature
Project description:Endometrial cancer (EC) is one of the most common gynecologic malignancies. To identify potential prognostic biomarkers for EC, we analyzed the relationship between the EC tumor microenvironment and gene expression profiles. Using the ESTIMATE R tool, we found that immune and stromal scores correlated with clinical data and the prognosis of EC patients. Based on the immune and stromal scores, 387 intersection differentially expressed genes were identified. Eight immune-related genes were then identified using two machine learning algorithms. Functional enrichment analysis revealed that these genes were mainly associated with T cell activation and response. Kaplan-Meier survival analysis showed that expression of TMEM150B, CACNA2D2, TRPM5, NOL4, CTSW, and SIGLEC1 significantly correlated with overall survival times of EC patients. In addition, using the TIMER algorithm, we found that expression of TMEM150B, SIGLEC1, and CTSW correlated positively with the tumor infiltration levels of B cells, CD8+ T cells, CD4+ T cells, macrophages, and dendritic cells. These findings indicate that the composition of the tumor microenvironment affects the clinical outcomes of EC patients, and suggests that it may provide a basis for development of novel prognostic biomarkers and immunotherapies for EC patients.
| S-EPMC7066904 | biostudies-literature
Project description:BackgroundOvarian cancer is one of the lethal gynecological diseases in women. However, using tumor microenvironment related genes to identify prognostic signature of ovarian cancer has not been discussed in detail.MethodsThe mRNA profiles of 386 ovarian cancer patients were retrieved from The Cancer Genome Atlas. Univariate Cox regression and LASSO Cox regression analyses were performed and 14 optimized prognostic genes related to tumor microenvironment were identified.ResultsThe multivariate Cox hazards regression showed risk score was an independent prognostic signature for ovarian cancer. Nomogram model could reliably predict the patients' survival. Furthermore, M1 macrophages, M2 macrophages, and follicular helper T cells, differentially expressed between the high- and low-risk groups, were found to be associated with the risk score.ConclusionCTL-associated antigen 4 (CTLA4) and indoleamine 2,3-Dioxygenase 1 (IDO1), which were previously shown to be important immune checkpoints, probably contribute to the immunosuppressive microenvironment aberration. This study may shed light on the prognosis of ovarian cancer.
| S-EPMC8284754 | biostudies-literature
Project description:Background and aimColorectal cancer (CRC) ranks the second most common cause of cancer-related mortality worldwide. Ferroptosis, a recently discovered form of programmed cell death different from other, raises promising novel opportunities for therapeutic intervention of CRC. This study intended to systematically assess the prognosis value and multiple roles of the ferroptosis-related genes in the tumor immune microenvironment of CRC.Materials and methodsOf 1,192 CRC patients with complete information from the public datasets (TCGA CRC, GEO GSE39582 and GSE17538 cohorts) were selected for analysis. Firstly, K-means consensus clustering was performed to identify ferroptosis-associated subtypes in CRC patients. Subsequently, we constructed a risk signature based on ferroptosis-related genes in TCGA cohort and acquired its validation in two GEO cohorts. Additionally, we established a nomogram integrating the risk signature and clinical factors to improve risk assessment of CRC patients.ResultsFive molecular subtypes were identified by consensus clustering for ferroptosis-related genes. There were significant differences in the overall survival, immune cells infiltration status and PD1/PD-L1 mRNA among the five clusters. Then, a risk signature based on the ten-gene was constructed which could distinguish effectively high-risk group among CRC patients in both training and validation sets. The high-risk patients were more likely to have an inhibitory immune microenvironment and lower stemness features. A prognostic nomogram integrated risk signature and clinicopathological features could be used as a more accurate prognostic prediction visualization tool than TNM stage alone.ConclusionThis ferroptosis risk signature may accurately differentiate between different risk populations and predict the prognosis of CRC. Besides, this study elucidated the crucial role of ferroptosis in tumor immune microenvironment.
| S-EPMC8286063 | biostudies-literature
Project description:Lung adenocarcinoma (LUAD) is a major subtype of non-small cell lung cancer. Despite significant progress in its diagnosis and treatment, the mortality and morbidity rate of LUAD remains high worldwide. The aim of the present study was to perform a systematic investigation of the tumor microenvironment (TME) and identify TME-related genes of prognostic value in patients with LUAD. Firstly, the immune scores and stromal scores of patients with LUAD from The Cancer Genome Atlas were calculated using the Estimation of STromal and Immune cells in MAlignant Tumors using Expression data algorithm, and a total of 281 prognostic TME-related genes were identified. Subsequently, functional analysis and protein-protein interaction network analysis revealed that these genes were mainly related to immune response, inflammatory response and chemotaxis. Finally, two independent LUAD cohorts from the Gene Expression Omnibus database were used to validate these genes, and 4 genes (GTPase IMAP family member 1, T-cell surface glycoprotein CD1b, integrin alpha-L and leukocyte surface antigen CD53) were identified, and downregulation of these genes was indicated to be associated with poor overall survival rate in patients with LUAD. In conclusion, a comprehensive analysis of TME was performed and 4 prognostic TME-related genes in patients with LUAD were identified.
| S-EPMC7377199 | biostudies-literature
Project description:Background: Esophageal squamous cell carcinoma (ESCC) is the most prevalent histological type of esophageal cancer, but there is a lack of definite prognostic markers for this cancer. Methods: We used the ESTIMATE algorithm to access the tumor microenvironment (TME) of ESCC cases deposited in the TCGA database, and identified TME-related prognostic genes using Cox regression analysis. A least absolute shrinkage and selector operation or LASSO algorithm was used to identify key prognostic genes. Risk scores were calculated, and a clinical predictive model was constructed to evaluate the prognostic value of TME-related genes. Results: We found that high immune and stromal scores were significantly associated with poor overall survival (p < 0.05). We identified a total of 1,151 TME-related differently expression genes, among which 67 were prognosis-related genes. Through the LASSO method, 13 key prognostic genes were selected, namely, ADAMTS16, LOC51089, CH25H, CORO2B, DLGAP1, GYS2, HAL, MXRA8, NPTX1, OTX1, RET, SLC24A2, and SPI1, and a 13-gene risk score was constructed. A higher score was indicative of a poorer prognosis than a lower risk score (hazard ratio = 8.21, 95% confidence interval: 2.56-26.31; P < 0.001). The risk score was significantly correlated with immune/stromal scores and various types of infiltrating immune cells, including CD8 cells, regulatory T cells, and resting macrophages. Conclusion: We characterized the tumor microenvironment in ESCC, and identified the key prognosis genes. The risk score based on the expression profiles of these genes is proposed as an indicator of TME status and is instrumental in predicting patient prognosis.
| S-EPMC7767869 | biostudies-literature