Background Breast cancer is a highly heterogeneous disease, this poses challenges for classification and management. Long non-coding RNAs play acrucial role in the breast cancersdevelopment and progression, especially in tumor-related immune processes which have become the most rapidly investigated area. Therefore, we aimed at developing an immune-related lncRNA signature to improve the prognosis prediction of breast cancer. Methods We obtained breast cancer patient samples and corresponding clinical data from The Cancer Genome Atlas (TCGA) database. Immune-related lncRNAs were screened by co-expression analysis of immune-related genes which were downloaded from the Immunology Database and Analysis Portal (ImmPort). Clinical patient samples were randomly separated into training and testing sets. In the training set, univariate Cox regression analysis and LASSO regression were utilized to build a prognostic immune-related lncRNA signature. The signature was validated in the training set, testing set, and whole cohorts by the Kaplan–Meier log-rank test, time-dependent ROC curve analysis, principal component analysis, univariate andmultivariate Cox regression analyses. Results A total of 937 immune- related lncRNAs were identified, 15 candidate immune-related lncRNAs were significantly associated with overall survival (OS). Eight of these lncRNAs (OTUD6B-AS1, AL122010.1, AC136475.2, AL161646.1, AC245297.3, LINC00578, LINC01871, AP000442.2) were selected for establishment of the risk prediction model. The OS of patients in the low-risk group was higher than that of patients in the high-risk group (p = 1.215e − 06 in the training set; p = 0.0069 in the validation set; p = 1.233e − 07 in whole cohort). The time-dependent ROC curve analysis revealed that the AUCs for OS in the first, eighth, and tenth year were 0.812, 0.81, and 0.857, respectively, in the training set, 0.615, 0.68, 0.655 in the validation set, and 0.725, 0.742, 0.741 in the total cohort. Multivariate Cox regression analysis indicated the model was a reliable and independent indicator for the prognosis of breast cancer in the training set (HR = 1.432; 95% CI 1.204–1.702, p < 0.001), validation set (HR = 1.162; 95% CI 1.004–1.345, p = 0.044), and whole set (HR = 1.240; 95% CI 1.128–1.362, p < 0.001). GSEA analysis revealed a strong connection between the signature and immune-related biological processes and pathways. Conclusions We constructed and verified a robust signature of 8 immune-related lncRNAs for the prediction of breast cancer patient survival.

Cancer stem cells (CSCs) are subsets of cells with the ability of self-renewal and differentiation in neoplasm, which are considered to be related to tumor heterogeneity. It has been reported that CSCs act on tumorigenesis and tumor biology of triple-negative breast cancer (TNBC). However, the key genes that cause TNBC showing stem cell characteristics are still unclear. We combined the RNA sequencing (RNA-seq) data from The Cancer Genome Atlas (TCGA) database and mRNA expression-based stemness index (mRNAsi) to further analyze mRNAsi with regard to molecular subtypes, tumor depth, and pathological staging characteristics of breast cancer (BC). Secondly, we extract the differential gene expression of tumor vs. normal group and TNBC vs. other subtypes of BC group, respectively, and intersect them to achieve precise results. We used a weighted gene coexpression network analysis (WGCNA) to screen significant gene modules and the functions of selected genes including BIRC5, CDC25A, KIF18B, KIF2C, ORC1, RAD54L, and TPX2 were carried out through gene ontology (GO) functional annotation. The Oncomine, bc-GenExMiner v4.4, GeneMANIA, Kaplan-Meier Plotter (KM-plotter), and GEPIA were used to verify the expression level and functions of key genes. In this study, we found that TNBC had the highest stem cell characteristics in BC compared with other subtypes. The lower the mRNAsi score, the better the overall survival and treatment outcome. Seven key genes of TNBC were screened and functional annotation indicated that there were strong correlations between them, relating to nuclear division, organelle fission, mitotic nuclear division, and other events that determine cell fate. Among these genes, we found four genes that were highly associated with adverse survival events. Seven key genes identified in this study were found to be closely related to the maintenance of TNBC stemness, and the overexpression of four showed earlier recurrence. The overall survival (OS) curves of all key genes between differential expression level crossed at around nine-year follow-up, which was consistent with the trend of the OS curve related to mRNAsi. These findings may provide new ideas for screening therapeutic targets in order to depress TNBC stemness.

To study the expression of MRPL13 in breast cancer tissues using TCGA database, analyze the correlation between the expression and clinicopathological characteristics of patients, and explore the role of MRPL13 in the development of breast cancer (BC). Methods: The BC mRNA data and clinical information were downloaded from TCGA database. The correlation between MRPL13 expression and clinicopathological parameters was analyzed. Cox regression multivariate analysis was used to explore the factors affecting the prognosis of BC patients. The UALCAN database was used to analyze the expression level of MRPL13 in BC and its relationship with clinical pathological factors. The GSEA method was used to predict the possible regulatory pathways of MRPL13. Immune responses of MRPL13 expression were analyzed using TISIDB and CIBERSORT. Additionally, GEPIA, K-M survival analysis and data from the HPA were used to validate the outcomes. Results: The expression of MRPL13 in BC tissues was significantly higher than normal counterparts, patients with low MRPL13 expression had a better survival prognosis, also indicated an independent prognostic factor. GSEA analysis showed that the regulation of cell migration, positive regulation of endothelial cell migration, and Notch signaling pathway were enriched in tissues with low expression of MRPL13. Additionally, depleting MRPL13 expression inhibited invasion in MCF-10A and MCF-7 cells. Furthermore, PCR showed that MRPL13 affected VEGFA and MMP gene expression. CIBERSORT analysis revealed that the amount of NK cells decreased when MRPL13 expression was high. Conclusion: The expression of MRPL13 mRNA is upregulated in BC tissues, and the expression level of MRPL13 is significantly related to the clinicopathological factors of patients. High MRPL13 expression is a poor prognostic factor for BC, and it can be used as a molecular marker for prognosis judgment and as a potential therapeutic target.