Background: Tumor microenvironment is essential for breast cancer progression and metastasis. Our study sets out to examine the genes affecting stromal and immune infiltration in breast cancer progression and prognosis.Materials and Methods: This work provides an approach for quantifying stromal and immune scores by using ESTIMATE algorithm based on gene expression matrix of breast cancer patients in TCGA database. We found differentially expressed genes (DEGs) through limma R package. Functional enrichments were accessed through Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. Besides, we constructed a protein-protein network, identified several hub genes in Cytoscape, and discovered functionally similar genes in GeneMANIA. Hub genes were validated with prognostic data by Kaplan-Meier analysis both in The Cancer Genome Atlas (TCGA) database and Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) database and a meta-analysis of hub genes prognosis data was utilized in multiple databases. Furthermore, their relationship with infiltrating immune cells was evaluated by Tumor IMmune Estimation Resource (TIMER) web tool. Cox regression was utilized for overall survival (OS) and recurrence-free survival (RFS) in TCGA database and OS in METABRIC database in order to evaluate the impact of stromal and immune scores on patients prognosis.Results: One thousand and eighty-five breast cancer patients were investigated and 480 differentiated expressed genes (DEGs) were found based on the analysis of mRNA expression profiles. Functional analysis of DEGs revealed their potential functions in immune response and extracellular interaction. Protein-protein interaction network gave evidence of 10 hub genes. Some of the hub genes could be used as predictive markers for patients prognosis. In this study, we found that tumor purity and specific immune cells infiltration varied in response to hub genes expression. The multivariate cox regression highlighted the fact that immune score played a detrimental role in overall survival (HR = 0.45, 95% CI: 0.27–0.74, p = 0.002) and recurrence-free survival (HR = 0.41, 95% CI: 0.22–0.77, p = 0.006) in TCGA database. These result was confirmed in METABRIC database that immune score was a protector of OS (HR = 0.88, 95% CI: 0.77–0.99, p = 0.039).Conclusions: Our findings promote a better understanding of the potential genes behind the regulation of tumor microenvironment and cells infiltration. Immune score should be considered as a prognostic factor for patients' survival.

Row 6 shows the estimate of the probability of not metastasizing given the subjects made the decision recommended by the model, and Row 7 shows 95% confidence interval for that estimate. Rows 8 and 9 show the same values of individuals who did not make the decision recommended by the model.

The variables in the LSDS-5YDM that are included in this study.

Row 6 shows the estimate of the probability of not metastasizing given the subjects made the decision recommended by the model, and Row 7 shows 95% confidence interval for that estimate. Rows 8 and 9 show the same values of individuals who did not make the decision recommended by the model.

Breast cancer (BC) is the most common cancer among women with high morbidity and mortality. Therefore, new research is still needed for biomarker detection. GSE101124 and GSE182471 datasets were obtained from the Gene Expression Omnibus (GEO) database to evaluate differentially expressed circular RNAs (circRNAs). The Cancer Genome Atlas (TCGA) and Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) databases were used to identify the significantly dysregulated microRNAs (miRNAs) and genes considering the Prediction Analysis of Microarray classification (PAM50). The circRNA-miRNA-mRNA relationship was investigated using the Cancer-Specific CircRNA, miRDB, miRTarBase, and miRWalk databases. The circRNA–miRNA–mRNA regulatory network was annotated using Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database. The protein-protein interaction network was constructed by the STRING database and visualized by the Cytoscape tool. Then, raw miRNA data and genes were filtered using some selection criteria according to a specific expression level in PAM50 subgroups. A bottleneck method was utilized to obtain highly interacted hub genes using cytoHubba Cytoscape plugin. The Disease-Free Survival and Overall Survival analysis were performed for these hub genes, which are detected within the miRNA and circRNA axis in our study. We identified three circRNAs, three miRNAs, and eighteen candidate target genes that may play an important role in BC. In addition, it has been determined that these molecules can be useful in the classification of BC, especially in determining the basal-like breast cancer (BLBC) subtype. We conclude that hsa_circ_0000515/miR-486-5p/SDC1 axis may be an important biomarker candidate in distinguishing patients in the BLBC subgroup of BC.