Dataset Information

In silico analysis of differentially expressed genesets in metastatic breast cancer identifies potential prognostic biomarkers.

ABSTRACT:

Background

Identification of specific biological functions, pathways, and appropriate prognostic biomarkers is essential to accurately predict the clinical outcomes of and apply efficient treatment for breast cancer patients.

Methods

To search for metastatic breast cancer-specific biological functions, pathways, and novel biomarkers in breast cancer, gene expression datasets of metastatic breast cancer were obtained from Oncomine, an online data mining platform. Over- and under-expressed genesets were collected and the differentially expressed genes were screened from four datasets with large sample sizes (N > 200). They were analyzed for gene ontology (GO), KEGG pathway, protein-protein interaction, and hub gene analyses using online bioinformatic tools (Enrichr, STRING, and Cytoscape) to find enriched functions and pathways in metastatic breast cancer. To identify novel prognostic biomarkers in breast cancer, differentially expressed genes were screened from the entire twelve datasets with any sample sizes and tested for expression correlation and survival analyses using online tools such as KM plotter and bc-GenExMiner.

Results

Compared to non-metastatic breast cancer, 193 and 144 genes were differentially over- and under-expressed in metastatic breast cancer, respectively, and they were significantly enriched in regulating cell death, epidermal growth factor receptor signaling, and membrane and cytoskeletal structures according to the GO analyses. In addition, genes involved in progesterone- and estrogen-related signalings were enriched according to KEGG pathway analyses. Hub genes were identified via protein-protein interaction network analysis. Moreover, four differentially over-expressed (CCNA2, CENPN, DEPDC1, and TTK) and three differentially under-expressed genes (ABAT, LRIG1, and PGR) were further identified as novel biomarker candidate genes from the entire twelve datasets. Over- and under-expressed biomarker candidate genes were positively and negatively correlated with the aggressive and metastatic nature of breast cancer and were associated with poor and good prognosis of breast cancer patients, respectively.

Conclusions

Transcriptome datasets of metastatic breast cancer obtained from Oncomine allow the identification of metastatic breast cancer-specific biological functions, pathways, and novel biomarkers to predict clinical outcomes of breast cancer patients. Further functional studies are needed to warrant validation of their roles as functional tumor-promoting or tumor-suppressing genes.

SUBMITTER: Kim J

PROVIDER: S-EPMC8235641 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundExosomes are involved in cell-to-cell communication, neovascularization, cancer metastasis, and drug resistance, which all play an important role in the occurrence and progression of hepatocellular carcinoma (HCC). Because there are few mechanistic studies about the function of exosomes in HCC, the goals of this study were to identify exosome-related genes in HCC, to establish a reliable prognostic model for HCC, and to explore underlying mechanisms.MethodsThe exoRBase and The Cancer Genome Atlas (TCGA) databases were used to analyze differentially expressed genes (DEGs). Cox regression and least absolute shrinkage and selection operator analyses were used to identify DEGs closely related to the overall survival of patients with HCC. An exosome-related prognostic model was then constructed in TCGA and validated in the International Cancer Genome Consortium database. A nomogram was developed to predict survival. CIBERSORT was used to estimate the abundance of different types of immune cells. Immunotherapy-related DEGs were used to predict the effect of immunotherapy.ResultsForty-eight exosome-related DEGs were obtained; of them, five [exportin 1 (XPO1), lysosomal thiol reductase (IFI30), F-box protein 16 (FBXO16), calmodulin 1 (CALM1), MORC family CW-type zinc finger 3 (MORC3)] were selected to construct a predictive model. Patients with HCC were then divided into low- and high-risk groups using the best cut-off value, as determined by the X-tile software. Prognosis was significantly poorer in the high-risk than in the low-risk group (P=0.009; hazard ratio =2.65). Features related to exosomes were found to positively regulate immune response. Further analysis showed a higher risk score was associated with higher expression of immune checkpoint molecules, including programmed death ligand 1 (PD-L1), programmed death ligand 2 (PD-L2), T cell Ig and ITIM domain (TIGIT), and indoleamine-2,3-dioxygenase 1 (IDO1).ConclusionsThis study has identified a novel signature based on exosome-related genes that has potential as a prognostic biomarker for HCC. Our research provides an immunological perspective for the development of precision treatment for HCC.

Project description:Gastric cancer (GC) is one of the most common malignancies worldwide. To the best of our knowledge, no biomarkers have been widely accepted for the early diagnosis and prognostic prediction of GC. This study aimed to identify potential novel prognostic biomarkers for GC. The dataset GSE29272, which originates from the public database Gene Expression Omnibus, was employed in the present study. The online tool GEO2R was used to calculate the differentially expressed genes (DEGs) in GSE29272 between tumour tissues and adjacent tissues. CytoHubba and MCODE plugins of Cytoscape software were used to obtain hub genes and modules of DEGs. The online tools Database for Annotation, Visualisation and Integrated Discovery and Search Tool for the Retrieval of Interacting Genes were employed to conduct Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes pathway analysis, and to construct protein-protein interaction networks. A total of 117 DEGs were extracted from GSE29272. In addition, 15 hub genes and seven modules were identified in the 117 DEGs. The enrichment analysis revealed that they were mainly enriched in GO biological process and cellular component domains, and the 'ECM-receptor interaction', 'focal adhesion', 'metabolism of xenobiotics by cytochrome P450' and 'drug metabolism' pathways. The hub genes asporin (ASPN), collagen type I α1 chain (COL1A1), fibronectin 1 (FN1), versican (VCAN) and mucin 5AC (MUC5AC) were demonstrated to have prognostic value for patients with GC. The ASPN and VCAN genes were significantly associated with overall survival and disease-free survival (log-rank P=0.025, 0.038, 0.0014 and 0.015, respectively). COL1A1 and FN1 were significantly associated with overall survival (log-rank P=0.013 and 0.05, respectively), and MUC5AC was significantly associated with disease-free survival (log-rank P=0.027). Results from the present study suggested that ASPN, COL1A1, FN1, VCAN and MUC5AC may represent novel prognostic biomarkers for GC.

Project description:INTRODUCTION: MicroRNAs (miRNAs) are a group of small noncoding RNAs involved in the regulation of gene expression. As such, they regulate a large number of cellular pathways, and deregulation or altered expression of miRNAs is associated with tumorigenesis. In the current study, we evaluated the feasibility and clinical utility of circulating miRNAs as biomarkers for the detection and staging of breast cancer. METHODS: miRNAs were extracted from a set of 84 tissue samples from patients with breast cancer and eight normal tissue samples obtained after breast-reductive surgery. After reverse transcription and preamplification, 768 miRNAs were profiled by using the TaqMan low-density arrays. After data normalization, unsupervised hierarchical cluster analysis (UHCA) was used to investigate global differences in miRNA expression between cancerous and normal samples. With fold-change analysis, the most discriminating miRNAs between both tissue types were selected, and their expression was analyzed on serum samples from 20 healthy volunteers and 75 patients with breast cancer, including 16 patients with untreated metastatic breast cancer. miRNAs were extracted from 200 ?l of serum, reverse transcribed, and analyzed in duplicate by using polymerase chain reaction (qRT-PCR). RESULTS: UHCA showed major differences in miRNA expression between tissue samples from patients with breast cancer and tissue samples from breast-reductive surgery (P < 0.0001). Generally, miRNA expression in cancerous samples tends to be repressed when compared with miRNA expression in healthy controls (P = 0.0685). The four most discriminating miRNAs by fold-change (miR-215, miR-299-5p, miR-411, and miR-452) were selected for further analysis on serum samples. All miRNAs at least tended to be differentially expressed between serum samples from patients with cancer and serum samples from healthy controls (miR-215, P = 0.094; miR-299-5P, P = 0.019; miR-411, P = 0.002; and miR-452, P = 0.092). For all these miRNAs, except for miR-452, the greatest difference in expression was observed between serum samples from healthy volunteers and serum samples from untreated patients with metastatic breast cancer. CONCLUSIONS: Our study provides a basis for the establishment of miRNAs as biomarkers for the detection and eventually staging of breast cancer through blood-borne testing. We identified and tested a set of putative biomarkers of breast cancer and demonstrated that altered levels of these miRNAs in serum from patients with breast cancer are particularly associated with the presence of metastatic disease.

Project description:Introduction: Abnormal status of gene expression plays an important role in tumorigenesis, progression and metastasis of breast cancer. Mechanisms of gene silence or activation were varied. Methylation of genes may contribute to alteration of gene expression. This study aimed to identify differentially expressed hub genes which may be regulated by DNA methylation and evaluate their prognostic value in breast cancer by bioinformatic analysis. Methods: GEO2R was used to obtain expression microarray data from GSE54002, GSE65194 and methylation microarray data from GSE20713, GSE32393. Differentially expressed-aberrantly methylated genes were identified by FunRich. Biological function and pathway enrichment analysis were conducted by DAVID. PPI network was constructed by STRING and hub genes was sorted by Cytoscape. Expression and DNA methylation of hub genes was validated by UALCAN and MethHC. Clinical outcome analysis of hub genes was performed by Kaplan Meier-plotter database for breast cancer. IHC was performed to analyze protein levels of EXO1 and Kaplan-Meier was used for survival analysis. Results: 677 upregulated-hypomethylated and 361 downregulated-hypermethylated genes were obtained from GSE54002, GSE65194, GSE20713 and GSE32393 by GEO2R and FunRich. The most significant biological process, cellular component, molecular function enriched and pathway for upregulated-hypomethylated genes were viral process, cytoplasm, protein binding and cell cycle respectively. For downregulated-hypermethylated genes, the result was peptidyl-tyrosine phosphorylation, plasma membrane, transmembrane receptor protein tyrosine kinase activity and Rap1 signaling pathway (All p< 0.05). 12 hub genes (TOP2A, MAD2L1, FEN1, EPRS, EXO1, MCM4, PTTG1, RRM2, PSMD14, CDKN3, H2AFZ, CCNE2) were sorted from 677 upregulated-hypomethylated genes. 4 hub genes (EGFR, FGF2, BCL2, PIK3R1) were sorted from 361 downregulated-hypermethylated genes. Differential expression of 16 hub genes was validated in UALCAN database (p<0.05). 7 in 12 upregulated-hypomethylated and 2 in 4 downregulated-hypermethylated hub genes were confirmed to be significantly hypomethylated or hypermethylated in breast cancer using MethHC database (p<0.05). Finally, 12 upregulated hub genes (TOP2A, MAD2L1, FEN1, EPRS, EXO1, MCM4, PTTG1, RRM2, PSMD14, CDKN3, H2AFZ, CCNE2) and 3 downregulated genes (FGF2, BCL2, PIK3R1) contributed to significant unfavorable clinical outcome in breast cancer (p<0.05). High expression level of EXO1 protein was significantly associated with poor OS in breast cancer patients (p=0.03). Conclusion: Overexpression of TOP2A, MAD2L1, FEN1, EPRS, EXO1, MCM4, PTTG1, RRM2, PSMD14, CDKN3, H2AFZ, CCNE2 and downregulation of FGF2, BCL2, PIK3R1 might serve as diagnosis and poor prognosis biomarkers in breast cancer by more research validation. EXO1 was identified as an individual unfavorable prognostic factor. Methylation might be one of the major causes leading to abnormal expression of those genes. Functional analysis and pathway enrichment analysis of those genes would provide novel ideas for breast cancer research.