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.

Project description:BackgroundCervical cancer belongs to one of the most common female cancers; yet, the exact underlying mechanisms are still elusive. Recently, microarray and sequencing technologies have been widely used for screening biomarkers and molecular mechanism discovery in cancer studies. In this study, we aimed to analyse the microarray datasets using comprehensive bioinformatics tools and identified novel biomarkers associated with the prognosis of patients with cervical cancer.MethodsThe differentially expressed genes (DEGs) from Gene Expression Omnibus (GEO) datasets including GSE138080, GSE113942 and GSE63514 were analysed using GEO2R tool. The functional enrichment analysis was performed using g:Profiler tool. The protein-protein interaction (PPI) network construction and hub genes identification were performed using the STRING database and Cytoscape software, respectively. The hub genes were subjected to expression and survival analysis in the cervical cancer. The EdU incorporation and Cell Counting Kit-8 assays were performed to evaluate the effects of hub gene knockdown on the proliferation of cervical cancer cells.ResultsA total of 89 overlapping DEGs (63 up-regulated and 26 down-regulated genes) were identified in the microarray datasets. The functional enrichment analysis indicated that the overlapping DEGs were mainly associated with "DNA replication" and "cell cycle". Furthermore, the PPI network analysis revealed that the network contains 87 nodes and 309 edges. Sub-module analysis using the Molecular Complex Detection tool identified 21 hub genes from the PPI network. The expression levels of the 21 hub genes were all up-regulated in the cervical cancer tissues when compared to normal cervical tissues as analysed by GEPIA tool. The survival analysis showed that the low expression of cell division cycle 45 (CDC45), GINS complex subunit 2 (GINS2), minichromosome maintenance complex component 2 (MCM2) and proliferating cell nuclear antigen (PCNA) was significantly correlated with the shorter overall survival of patients with cervical cancer. Moreover, the protein expression levels of GINS2, MCM2 and PCNA, but not CDC45, were significantly up-regulated in the cervical cancer tissues when compared to normal cervical tissues. Finally, knockdown of MCM2 significantly suppressed the proliferation of HeLa and SiHa cells.ConclusionIn conclusion, we screened a total of 89 overlapping DEGs from the GEO datasets, and further analysis identified four hub genes (CDC45, GINS2, MCM2 and PCNA) that were likely associated with the prognosis of patients with cervical cancer. MCM2 knockdown repressed the cervical cancer cell proliferation. The current findings may provide novel insights into understanding the pathophysiology of cervical cancer and develop therapeutic targets for patients with cervical cancer.