Dataset Information

Bioinformatics identification of CCL8/21 as potential prognostic biomarkers in breast cancer microenvironment.

ABSTRACT:

Background

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.

Methods

In our study, we analyzed the expression and survival data of CC chemokines in patients with BC using several bioinformatics analyses tools.

Results

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.

Conclusion

Integrative bioinformatics analysis demonstrated CCL8/21 as potential prognostic biomarkers in BC microenvironment.

SUBMITTER: Chen B

PROVIDER: S-EPMC7687043 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

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.

Project description:PurposeThe prognosis of breast cancer (BC) patients who develop into brain metastases (BMs) is very poor. Thus, it is of great significance to explore the etiology of BMs in BC and identify the key genes involved in this process to improve the survival of BC patients with BMs.Patients and methodsThe gene expression data and the clinical information of BC patients were downloaded from TCGA and GEO database. Differentially expressed genes (DEGs) in TCGA-BRCA and GSE12276 were overlapped to find differentially expressed metastatic genes (DEMGs). The protein-protein interaction (PPI) network of DEMGs was constructed via STRING database. ClusterProfiler R package was applied to perform the gene ontology (GO) enrichment analysis of DEMGs. The univariate Cox regression analysis and the Kaplan-Meier (K-M) curves were plotted to screen DEMGs associated with the overall survival and the metastatic recurrence survival, which were identified as the key genes associated with the BMs in BC. The immune infiltration and the expressions of immune checkpoints for BC patients with brain relapses and BC patients with other relapses were analyzed respectively. The correlations among the expressions of key genes and the differently infiltrated immune cells or the differentially expressed immune checkpoints were calculated. The gene set enrichment analysis (GSEA) of each key gene was conducted to investigate the potential mechanisms of key genes involved in BC patients with BMs. Moreover, CTD database was used to predict the drug-gene interaction network of key genes.ResultsA total of 154 DEGs were identified in BC patients at M0 and M1 in TCGA database. A total of 667 DEGs were identified in BC patients with brain relapses and with other relapses. By overlapping these DEGs, 17 DEMGs were identified, which were enriched in the cell proliferation related biological processes and the immune related molecular functions. The univariate Cox regression analysis and the Kaplan-Meier curves revealed that CXCL9 and GPR171 were closely associated with the overall survival and the metastatic recurrence survival and were identified as key genes associated with BMs in BC. The analyses of immune infiltration and immune checkpoint expressions showed that there was a significant difference of the immune microenvironment between brain relapses and other relapses in BC. GSEA indicated that CXCL9 and GPR171 may regulate BMs in BC via the immune-related pathways.ConclusionOur study identified the key genes associated with BMs in BC patients and explore the underlying mechanisms involved in the etiology of BMs in BC. These findings may provide a promising approach for the treatments of BC patients with BMs.