Identification of a tumor microenvironment-associated prognostic gene signature in bladder cancer by integrated bioinformatic analysis.
Ontology highlight
ABSTRACT:
SUBMITTER: Liu Z
PROVIDER: S-EPMC8167492 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:BackgroundThe tumor microenvironment (TME) regulates the proliferation and metastasis of solid tumors and the effectiveness of immunotherapy against them. We investigated the prognostic role of TME-related genes based on transcriptomic data of bladder urothelial carcinoma (BLCA) and formulated a prediction model of TME-related signatures.MethodsMolecular subtypes were identified using the non-negative matrix factorization (NMF) algorithm based on TME-related genes from the TCGA database. TME-related genes with prognostic significance were screened with univariate Cox regression analysis and lasso regression. Nomogram was developed based on risk genes. Receiver operating characteristic (ROC) curve and decision curve analysis (DCA) were used for inner and outer validation of the model. Risk scores (RS) of patients were calculated and divided into high-risk group (HRG) and low-risk group (LRG) to compare the differences in clinical characteristics and PD-L1 treatment responsiveness between HRG and LRG.ResultsWe identified two molecular subtypes (C1 and C2) according to the NMF algorithm. There were significant differences in overall survival (OS) (p<0.05), progression-free survival (PFS) (p<0.05), and immune cell infiltration between the two subtypes. A total of eight TME-associated genes (CABP4, ZNF432, BLOC1S3, CXCL11, ANO9, OAS1, FBN2, CEMIP) with independent prognostic significance were screened to build prognostic risk models. Age (p<0.001), grade (p<0.001), and RS (p<0.001) were independent predictors of survival in BLCA patients. The developed RS nomogram was able to predict the prognosis of BLCA patients at 1, 3, and 5 years more potentially than the models of other investigators according to ROC and DCA. RS showed significantly higher values (p = 0.047) in patients with stable disease (SD)/progressive disease (PD) compared to patients with complete response (CR)/partial response (PR).ConclusionsWe successfully clustered and constructed predictive models for TME-associated genes and helped guide immunotherapy strategies.
| S-EPMC8921452 | biostudies-literature
Project description:In the worldwide, uterine corpus endometrial carcinoma (UCEC) is the sixth most common malignancy in women, and the number of women diagnosed is increasing. Kinase plays an important role in the occurrence and development of malignant tumors. However, the research about kinase in endometrial cancer is still unclear. Here, we first downloaded the gene expression data of 552 UCEC patients and 23 healthy endometrial tissues from The Cancer Genome Atlas (TCGA), obtained 538 kinase-related genes from the previous literature, and calculated 67 differentially expressed kinases. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were referenced to identify multiple important biological functions and signaling pathways related to 67 differentially expressed kinases. Using univariate Cox regression and Least absolute shrinkage and selection operator (LASSO), seven kinases (ALPK2, CAMKV, TTK, PTK6, MAST1, CIT, and FAM198B) were identified to establish a prognostic model of endometrial cancer. Then, patients were divided into high- and low-risk groups based on risk scores. Receiver operating characteristic (ROC) curves were plotted to evaluate that the model had a favorable predictive ability. Kaplan-Meier survival analysis suggested that high-risk groups experienced worse overall survival than low-risk groups. qRT-PCR and ISH assays confirmed the consistency between predicted candidate genes and real sample contents. CIBERSORT algorithm and ssGSEA were adopted to investigate the relationship between this signature and tumor immune microenvironment, and revealed that in low- and high-risk groups, the types of tumor-infiltrating immune cells and the immune cell-related functions were significantly different. In summary, a seven-gene signature risk model has been constructed, and could accurately predict the prognosis of UCEC, which may offer ideas and breakthrough points to the kinase-associated development of UCEC.
| S-EPMC9491090 | 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: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:Programmed cell death (PCD) refers to a molecularly regulated form of cell death that functions as an essential anticancer defense mechanism and serves as a target of anticancer therapies. Multiple types of PCD comprehensively regulate tumorigenesis and tumor progression and metastasis. However, a systemic exploration of the multiple types of PCD in cancers, especially bladder cancer, is lacking. In this study, we evaluated the expression pattern of genes associated with multiple types of PCD in bladder cancer using the "ssGSEA" method and conceptualized the multiple types of PCD as being collectively involved in "Pan-PCD". Based on the differentially expressed genes related to Pan-PCD, we developed a Pan-PCD-related prognostic signature (PPRPS) to predict patient prognosis via univariate and multivariate Cox regression analysis. The PPRPS is an independent prognostic factor, and the AUC (Area Under Curve) for 3-year overall survival was 0.748. Combined with age and stage, PPRPS displayed excellent predictive ability. Based on the PPRPS, higher levels of immune cell infiltration, tumor microenvironment, and immune checkpoint molecules were observed in the high-PPRPS group. Furthermore, PPRPS enabled accurate risk prediction for metastatic urothelial carcinoma after anti-PD-L1 monoclonal antibody treatment. Patients in the high-PPRPS group had poor prognoses. Docetaxel, staurosporine, and luminespib were identified as potentially effective drugs for high-PPRPS bladder cancer patients. In summary, we developed the Pan-PCD signature to improve the accuracy of bladder cancer prognostic predictions and to provide a novel classification method to guide treatment selection.
| S-EPMC9728529 | biostudies-literature
Project description:As major components of the tumor microenvironment (TME), tumor-associated macrophages (TAMs) play an exceedingly complicated role in tumor progression and tumorigenesis. However, few studies have reported the specific TAM gene signature in bladder cancer. Herein, this study focused on developing a TAM-related prognostic model in bladder cancer patients based on The Cancer Genome Atlas (TCGA) data. Weighted Gene Co-Expression Network Analysis (WGCNA) was used to identify key genes related to TAM (M2 macrophage). Gene ontology (GO) enrichment and the Kyoto Encyclopedia of Genes and Genomes (KEGG) signaling pathway analysis showed the functional categories of the key genes. Simultaneously, we used the Least Absolute Shrinkage and Selection Operator (LASSO) and univariate and multivariate Cox regressions to establish a TMA-related prognostic model containing six key genes: TBXAS1, GYPC, HPGDS, GAB3, ADORA3, and FOLR2. Subsequently, single-cell sequencing data downloaded from Gene Expression Omnibus (GEO) suggested that the six genes in the prognostic model were expressed in TAM specifically and may be involved in TAM polarization. In summary, our research uncovered six-TAM related genes that may have an effect on risk stratification in bladder cancer patients and could be regarded as potential TAM-related biomarkers.
| S-EPMC9582252 | biostudies-literature
Project description:The microenvironment plays a vital role in the tumor recurrence of neuroblastoma. This research aimed at exploring prognostic genes that are involved in neuroblastoma microenvironment. We used "estimate" R package to calculate the immune/stromal/ESTIMATE scores of each sample of ArrayExpress dataset E-MTAB-8248 based on the ESTIMATE algorithm. Then we found that immune/stromal/ESTIMATE scores were not correlated with age/chromosome 11q, but tumor stage, MYCN gene amplifications, and chromosome 1p. Samples were then divided into high- and low-score groups, and 280 common differentially expressed genes (DEGs) were identified. 64 potential prognostic genes were harvested through overall survival analysis from the common DEGs. 14 prognostic genes (ABCA6, SEPP1, SLAMF8, GPR171, ABCA9, ARHGAP15, IL7R, HLA-DPB1, GZMA, GPR183, CCL19, ITK, FGL2, and CD1C) were obtained after screening in two independent cohorts. GO and KEGG analysis discovered that common DEGs and 64 potential prognostic genes are mainly involved in T-cell activation, lymphocyte activation regulation, leukocyte migration, and the interaction of cytokines and cytokine receptors. Correlation analysis showed that all prognostic genes were negatively correlated with MYCN amplification. Cox analysis identified 5 independent prognostic genes (ARHGAP15, ABCA9, CCL19, SLAMF8, and CD1C).
| S-EPMC7501561 | biostudies-literature
Project description:BackgroundAccumulating evidences demonstrated tumor microenvironment (TME) of bladder cancer (BLCA) may play a pivotal role in modulating tumorigenesis, progression, and alteration of biological features. Currently we aimed to establish a prognostic model based on TME-related gene expression for guiding clinical management of BLCA.MethodsWe employed ESTIMATE algorithm to evaluate TME cell infiltration in BLCA. The RNA-Seq data from The Cancer Genome Atlas (TCGA) database was used to screen out differentially expressed genes (DEGs). Underlying relationship between co-expression modules and TME was investigated via Weighted gene co-expression network analysis (WGCNA). COX regression and the least absolute shrinkage and selection operator (LASSO) analysis were applied for screening prognostic hub gene and establishing a risk predictive model. BLCA specimens and adjacent tissues from patients were obtained from patients. Bladder cancer (T24, EJ-m3) and bladder uroepithelial cell line (SVHUC1) were used for genes validation. qRT-PCR was employed to validate genes mRNA level in tissues and cell lines.Results365 BLCA samples and 19 adjacent normal samples were selected for identifying DEGs. 2141 DEGs were identified and used to construct co-expression network. Four modules (magenta, brown, yellow, purple) were regarded as TME regulatory modules through WGCNA and GO analysis. Furthermore, seven hub genes (ACAP1, ADAMTS9, TAP1, IFIT3, FBN1, FSTL1, COL6A2) were screened out to establish a risk predictive model via COX and LASSO regression. Survival analysis and ROC curve analysis indicated our predictive model had good performance on evaluating patients prognosis in different subgroup of BLCA. qRT-PCR result showed upregulation of ACAP1, IFIT3, TAP1 and downregulation of ADAMTS9, COL6A2, FSTL1,FBN1 in BLCA specimens and cell lines.ConclusionsOur study firstly integrated multiple TME-related genes to set up a risk predictive model. This model could accurately predict BLCA progression and prognosis, which offers clinical implication for risk stratification, immunotherapy drug screen and therapeutic decision.
| S-EPMC8194149 | biostudies-literature
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.
| S-EPMC10839309 | biostudies-literature
Project description:Background: Bladder cancer (BLCA) is a common malignant tumor of the genitourinary system, and there is a lack of specific, reliable, and non-invasive tumor biomarker tests for diagnosis and prognosis evaluation. Homeobox genes play a vital role in BLCA tumorigenesis and development, but few studies have focused on the prognostic value of homeobox genes in BLCA. In this study, we aim to develop a prognostic signature associated with the homeobox gene family for BLCA. Methods: The RNA sequencing data, clinical data, and probe annotation files of BLCA patients were downloaded from the Gene Expression Omnibus database and the University of California, Santa Cruz (UCSC), Xena Browser. First, differentially expressed homeobox gene screening between tumor and normal samples was performed using the "limma" and robust rank aggregation (RRA) methods. The mutation data were obtained with the "TCGAmutation" package and visualized with the "maftools" package. Kaplan-Meier curves were plotted with the "survminer" package. Then, a signature was constructed by logistic regression analysis. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed using "clusterProfiler." Furthermore, the infiltration level of each immune cell type was estimated using the single-sample gene set enrichment analysis (ssGSEA) algorithm. Finally, the performance of the signature was evaluated by receiver-operating characteristic (ROC) curve and calibration curve analyses. Results: Six genes were selected to construct this prognostic model: TSHZ3, ZFHX4, ZEB2, MEIS1, ISL1, and HOXC4. We divided the BLCA cohort into high- and low-risk groups based on the median risk score calculated with the novel signature. The overall survival (OS) rate of the high-risk group was significantly lower than that of the low-risk group. The infiltration levels of almost all immune cells were significantly higher in the high-risk group than in the low-risk group. The average risk score for the group that responded to immunotherapy was significantly lower than that of the group that did not. Conclusion: We constructed a risk prediction signature with six homeobox genes, which showed good accuracy and consistency in predicting the patient's prognosis and response to immunotherapy. Therefore, this signature can be a potential biomarker and treatment target for BLCA patients.
| S-EPMC8334560 | biostudies-literature