Identification and validation of tumor microenvironment-related prognostic biomarkers in breast cancer
Ontology highlight
ABSTRACT: Background
Stromal cells and immune cells in tumor microenvironment (TME) have been reported to have significant value in the diagnosis and prognosis of cancers. We aimed to identify key biomarkers predicting survival in the TME of breast cancer. Methods
Cell type enrichment analysis was performed to estimate cell scores using the xCell method with gene expression data from public database. Least absolute shrinkage and selection operator (LASSO) regression was used to identify key signature from the cell scores. Results
Totally, 50 cells in TME had different scores between 1,078 breast cancer tissues and 112 adjacent normal tissues. We identified a 4-cell signature predicting breast cancer survival, including myocytes, natural killer T cell (NKT), conventional dendritic cell (cDC) and sebocytes, which was validated in the test set. Further analysis showed that cDC score was a key signature predicting prognosis of breast cancer. cDC score was significantly associated with molecular classification and stage of breast cancer, as well as expression level of Ki67. Spearman’s correlation analysis found that cDC score was inversely correlated with the expression level of HER2. High cDC score may predicate better pathological complete response rate. Mechanism analysis indicated high cDC score was associated with elevated immune activity; IL-2 was a key gene associated with high cDC score; and Breast cancer patients with high IL-2 expression had a longer survival time. Conclusions
In conclusion, cDC score was a key signature predicting prognosis for breast cancer. cDCs may exert antitumor effects by upregulating IL-2.
SUBMITTER: Zhong S
PROVIDER: S-EPMC8798301 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Background: The tumor microenvironment (TME) has been reported to have significant value in the diagnosis and prognosis of cancers. This study aimed to identify key biomarkers in the TME of luminal breast cancer (BC). Methods: We obtained immune scores (ISs) and stromal scores (SSs) for The Cancer Genome Atlas (TCGA) luminal BC cohort from the online ESTIMATE (Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data) portal. The relationships between ISs and SSs and the overall survival of luminal BC patients were assessed by the Kaplan-Meier method. The differentially expressed messenger RNAs (DEmRNAs) related to the ISs and SSs were subjected to functional enrichment analysis. Additionally, a competing endogenous RNA (ceRNA) network was constructed with differentially expressed microRNAs (DEmiRNAs) and long noncoding RNAs (DElncRNAs). Furthermore, a protein-protein interaction (PPI) network was established to analyze the DEmRNAs in the ceRNA network. Then, survival analysis of biomarkers involved in the ceRNA network was carried out to explore their prognostic value. Finally, these biomarkers were validated using the luminal BC dataset from the Gene Expression Omnibus (GEO) database. Results: The results showed that ISs were significantly associated with longer survival times of luminal BC patients. Functional enrichment analysis showed that the DEmRNAs were mainly associated with immune response, antigen binding, and the extracellular region. In the PPI network, the top 10 DEmRNAs were identified as hub genes that affected the TME of luminal BC. Finally, two DEmiRNAs, two DElncRNAs, and 17 DEmRNAs of the ceRNA network associated with the TME were shown to have prognostic value. Subsequently, the expression of 15 prognostic biomarkers was validated in one additional dataset (GSE81002). In particular, one lncRNA (GVINP1) and five mRNAs (CCDC69, DOCK2, IKZF1, JCHAIN, and NCKAP1L) were novel biomarkers. Conclusions: Our studies demonstrated that ISs were associated with the survival of luminal BC patients, and a set of novel biomarkers that might play a prognostic role in the TME of luminal BC was identified.
| S-EPMC7735391 | biostudies-literature
Project description:BackgroundHypoxia is a common microenvironment condition in most malignant tumors and has been shown to be associated with adverse outcomes of cervical cancer patients. In this study, we investigated the effects of hypoxia-related genes on tumor progress to characterize the tumor hypoxic microenvironment.MethodsWe retrieved a set of hypoxia-related genes from the Molecular Signatures Database and evaluated their prognostic value for cervical cancer. A hypoxia-based prognostic signature for cervical cancer was then developed and validated using tumor samples from two independent cohorts (TCGA-CESC and CGCI-HTMCP-CC cohorts). Finally, we validated the hypoxia prediction of ccHPS score in eight human cervical cancer cell lines treated with the hypoxic and normoxic conditions, and 286 tumor samples with hypoxic category (more or less) from Gene Expression Omnibus (GEO) database with accession GSE72723.ResultsA risk signature model containing nine hypoxia-related genes was developed and validated in cervical cancer. Further analysis showed that this risk model could be an independent prognosis factor of cervical cancer, which reflects the condition of the hypoxic tumor microenvironment and its remodeling of cell metabolism and tumor immunity. Furthermore, a nomogram integrating the novel risk model and lymphovascular invasion status was developed, accurately predicting the 1-, 3- and 5-year prognosis with AUC values of 0.928, 0.916 and 0.831, respectively. These findings provided a better understanding of the hypoxic tumor microenvironment in cervical cancer and insights into potential new therapeutic strategies in improving cancer therapy.
| S-EPMC9165826 | biostudies-literature
Project description:Thyroid Carcinoma (THCA) is the most common endocrine tumor that is mainly treated using surgery and radiotherapy. In addition, immunotherapy is a recently developed treatment option that has played an essential role in the management of several types of tumors. However, few reports exist on the use of immunotherapy to treat THCA. The study downloaded the miRNA, mRNA and lncRNA data for THCA patients from the TCGA database ( https://portal.gdc.cancer.gov/ ). Thereafter, the tumor samples were divided into cold and hot tumors, based on the immune score of the tumor microenvironment. Moreover, the differentially expressed lncRNAs and miRNAs were obtained. Finally, the study jointly constructed a ceRNA network through differential analysis of the mRNA data for cold and hot tumors. The study first assessed the level of immune infiltration in the THCA tumor microenvironment then divided the samples into cold and hot tumors, based on the immune score. Additionally, a total of 568 up-regulated and 412 down-regulated DEGs were screened by analyzing the differences between hot and cold tumors. Thereafter, the study examined the differentially expressed genes for lncRNA and miRNA. The results revealed 629 differentially expressed genes related to lncRNA and 114 associated with miRNA. Finally, a ceRNA network of the differentially expressed genes was constructed. The results showed a five-miRNA hubnet, i.e., hsa-mir-204, hsa-mir-128, hsa-mir-214, hsa-mir-150 and hsa-mir-338. The present study identified the immune-related mRNA, lncRNA and miRNA in THCA then constructed a ceRNA network. These results are therefore important as they provide more insights on the immune mechanisms in THCA. The findings also provides additional information for possible THCA immunotherapy.
| S-EPMC8355328 | biostudies-literature
Project description:BackgroundBreast cancer is an invasive disease with complex molecular mechanisms. Prognosis-related biomarkers are still urgently needed to predict outcomes of breast cancer patients.MethodsOriginal data were download from The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO). The analyses were performed using perl-5.32 and R-x64-4.1.1.ResultsIn this study, 1086 differentially expressed genes (DEGs) were identified in the TCGA cohort; 523 shared DEGs were identified in the TCGA and GSE10886 cohorts. Eight subtypes were estimated using non-negative matrix factorization clustering with significant differences seen in overall survival (OS) and progression-free survival (PFS) (P < 0.01). Univariate Cox analysis and least absolute shrinkage and selection operator (LASSO) regression analysis were performed to develop a related risk score related to the 17 DEGs; this score separated breast cancer into low- and high-risk groups with significant differences in survival (P < 0.01) and showed powerful effectiveness (TCGA all group: 1-year area under the curve [AUC] = 0.729, 3-year AUC = 0.778, 5-year AUC = 0.781). A nomogram prediction model was constructed using non-negative matrix factorization clustering, the risk score, and clinical characteristics. Our model was confirmed to be related with tumor microenvironment. Furthermore, DEGs in high-risk breast cancer were enriched in histidine metabolism (normalized enrichment score [NES] = 1.49, P < 0.05), protein export (NES = 1.58, P < 0.05), and steroid hormone biosynthesis signaling pathways (NES = 1.56, P < 0.05).ConclusionsWe established a comprehensive model that can predict prognosis and guide treatment.
| S-EPMC8833129 | 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:Breast cancer (BC) is the most common malignancy among females worldwide. The tumor microenvironment usually prevents effective lymphocyte activation and infiltration, and suppresses infiltrating effector cells, leading to a failure of the host to reject the tumor. CC chemokines play a significant role in inflammation and infection. In our study, we analyzed the expression and survival data of CC chemokines in patients with BC using several bioinformatics analyses tools. The mRNA expression of CCL2/3/4/5/7/8/11/17/19/20/22 was remarkably increased while CCL14/21/23/28 was significantly down-regulated in BC tissues compared with normal tissues. Methylation could down-regulate expression of CCL2/5/15/17/19/20/22/23/24/25/26/27 in BC. Low expression of CCL3/4/23 was found to be associated with drug resistance in BC. Results from Kaplan-Meier plotter and BC Gene-Expression Miner v4.2 (bcGenExMiner) v4.2 demonstrated that BC patients with high CCL8 and low CCL19/21/22 expression were more likely to have a worse prognosis. CCL8 expression was significantly up-regulated in BC tissues compared with normal tissues. High CCL8 expression was significantly correlated with negative PR, negative ER, positive nodal status, triple-negative BC subtype, basal-like BC subtype, triple-negative and basal-like BC subtype and high grades. CCL21 was down-regulated in BC, while high levels of CCL21 was associated with negative PR, triple-negative subtype, basal-like subtype and low tumor grade. Functional analysis demonstrated that CCL8 and CCL21 were involved in carcinogenesis, tumor immune escape and chemoresistance in BC. Integrative bioinformatics analysis demonstrated CCL8/21 as potential prognostic biomarkers in BC microenvironment.
| S-EPMC7687043 | 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: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:BackgroundLncRNA dysregulation and the tumor microenvironment (TME) have been shown to play a vital role in the progression and prognosis of colon cancer (CC). We aim to reveal the potential molecular mechanism from the perspective of lncRNA in the TME and provide the candidate biomarkers for CC prognosis.MethodsESTIMATE analysis was used to divide the CC patients into high and low immune or stromal score groups. The expression array of lncRNA was re-annotated by Seqmap. Microenvironment-associated lncRNAs were filtered through differential analysis. The m6A-associated lncRNAs were screened by Pearson correlation analysis. Lasso Cox regression analyses were performed to construct the m6A- and tumor microenvironment-related lncRNA prognostic model (m6A-TME-LM). Survival analysis was used to assess the prognostic efficacy of candidate lncRNAs. Enrichment analyses annotated the candidate genes' functions.ResultsWe obtained 25 common differentially expressed lncRNAs (DELs) associated with immune microenvironment and m6A-related genes for subsequent lasso analysis. Four out of these DELs were selected for the m6A-TME-LM. All the four lncRNAs were related to overall survival, and a test set testified the result. Further stratification analysis of the m6A-TME-LM retained its ability to predict OS for male and chemotherapy adjuvant patients and performed an excellent prognostic efficacy in the TNM stage III and IV subgroups. Network analysis also found the four lncRNAs mediated co-expression network was associated with tumor development.ConclusionWe constructed the m6A-TME-LM, which could provide a better prognostic prediction of CC.
| S-EPMC8231021 | biostudies-literature
Project description:BackgroundEsophageal squamous cell carcinoma (ESCC) is one of the most fatal cancers in the world. The 5-year survival rate of ESCC is <30%. However, few biomarkers can accurately predict the prognosis of patients with ESCC. We aimed to identify potential survival-associated biomarkers for ESCC to improve its poor prognosis.MethodsImmuneAI analysis was first used to access the immune cell abundance of ESCC. Then, ESTIMATE analysis was performed to explore the tumor microenvironment (TME), and differential analysis was used for the selection of immune-related differentially expressed genes (DEGs). Weighted gene coexpression network analysis (WGCNA) was used for selecting the candidate DEGs. Least absolute shrinkage and selection operator (LASSO) Cox regression was used to build the immune-cell-associated prognostic model (ICPM). Kaplan-Meier curve of survival analysis was performed to evaluate the efficacy of the ICPM.ResultsBased on the ESTIMATE and ImmuneAI analysis, we obtained 24 immune cells' abundance. Next, we identified six coexpression module that was associated with the abundance. Then, LASSO regression models were constructed by selecting the genes in the module that is most relevant to immune cells. Two test dataset was used to testify the model, and we finally, obtained a seven-genes survival model that performed an excellent prognostic efficacy.ConclusionIn the current study, we filtered seven key genes that may be potential prognostic biomarkers of ESCC, and they may be used as new factors to improve the prognosis of cancer.
| S-EPMC8573319 | biostudies-literature