Purpose: The m5C RNA methylation regulators are closely related to tumor proliferation, occurrence, and metastasis. This study aimed to investigate the gene expression, clinicopathological characteristics, and prognostic value of m5C regulators in triple-negative breast cancer (TNBC) and their correlation with the tumor immune microenvironment (TIM).Methods: The TNBC data, Luminal BC data and HER2 positive BC data set were obtained from The Cancer Genome Atlas and Gene Expression Omnibus, and 11 m5C RNA methylation regulators were analyzed. Univariate Cox regression and the least absolute shrinkage and selection operator regression models were used to develop a prognostic risk signature. The UALCAN and cBioportal databases were used to analyze the gene characteristics and gene alteration frequency of prognosis-related m5C RNA methylation regulators. Gene set enrichment analysis was used to analyze cellular pathways enriched by prognostic factors. The Tumor Immune Single Cell Hub (TISCH) and Timer online databases were used to explore the relationship between prognosis-related genes and the TIM.Results: Most of the 11 m5C RNA methylation regulators were differentially expressed in TNBC and normal samples. The prognostic risk signature showed good reliability and an independent prognostic value. Prognosis-related gene mutations were mainly amplified. Concurrently, the NOP2/Sun domain family member 2 (NSUN2) upregulation was closely related to spliceosome, RNA degradation, cell cycle signaling pathways, and RNA polymerase. Meanwhile, NSUN6 downregulation was related to extracellular matrix receptor interaction, metabolism, and cell adhesion. Analysis of the TISCH and Timer databases showed that prognosis-related genes affected the TIM, and the subtypes of immune-infiltrating cells differed between NSUN2 and NSUN6.Conclusion: Regulatory factors of m5C RNA methylation can predict the clinical prognostic risk of TNBC patients and affect tumor development and the TIM. Thus, they have the potential to be a novel prognostic marker of TNBC, providing clues for understanding the RNA epigenetic modification of TNBC.

The microbiome plays diverse roles in many diseases and can potentially contribute to cancer development. Breast cancer is the most commonly diagnosed cancer in women worldwide. Thus, we investigated whether the gut microbiota differs between patients with breast carcinoma and those with benign tumors. The DNA of the fecal microbiota community was detected by Illumina sequencing and the taxonomy of 16S rRNA genes. The α-diversity and β-diversity analyses were used to determine richness and evenness of the gut microbiota. Gene function prediction of the microbiota in patients with benign and malignant carcinoma was performed using PICRUSt. There was no significant difference in the α-diversity between patients with benign and malignant tumors ( P = 3.15e⁻¹ for the Chao index and P = 3.1e⁻¹ for the ACE index). The microbiota composition was different between the 2 groups, although no statistical difference was observed in β-diversity. Of the 31 different genera compared between the 2 groups, level of only Citrobacter was significantly higher in the malignant tumor group than that in benign tumor group. The metabolic pathways of the gut microbiome in the malignant tumor group were significantly different from those in benign tumor group. Furthermore, the study establishes the distinct richness of the gut microbiome in patients with breast cancer with different clinicopathological factors, including ER, PR, Ki-67 level, Her2 status, and tumor grade. These findings suggest that the gut microbiome may be useful for the diagnosis and treatment of malignant breast carcinoma.

Purpose: The role of 5-methylcytosine-related long non-coding RNAs (m5C-lncRNAs) in breast cancer (BC) remains unclear. Here, we aimed to investigate the prognostic value, gene expression characteristics, and correlation between m5C-lncRNA risk model and tumor immune cell infiltration in BC.Methods: The expression matrix of m5C-lncRNAs in BC was obtained from The Cancer Genome Atlas database, and the lncRNAs were analyzed using differential expression analysis as well as univariate and multivariate Cox regression analysis to eventually obtain BC-specific m5C-lncRNAs. A risk model was developed based on three lncRNAs using multivariate Cox regression and the prognostic value, accuracy, as well as reliability were verified. Gene set enrichment analysis (GSEA) was used to analyze the Kyoto Encyclopedia of Genes and Genomes signaling pathway enrichment of the risk model. CIBERSORT algorithm and correlation analysis were used to explore the characteristics of the BC tumor-infiltrating immune cells. Finally, reverse transcription-quantitative polymerase chain reaction was performed to detect the expression level of three lncRNA in clinical samples.Results: A total of 334 differential m5C-lncRNAs were identified, and three BC-specific m5C-lncRNAs were selected, namely AP005131.2, AL121832.2, and LINC01152. Based on these three lncRNAs, a highly reliable and specific risk model was constructed, which was proven to be closely related to the prognosis of patients with BC. Therefore, a nomogram based on the risk score was built to assist clinical decisions. GSEA revealed that the risk model was significantly enriched in metabolism-related pathways and was associated with tumor immune cell infiltration based on the analysis with the CIBERSORT algorithm.Conclusion: The efficient risk model based on m5C-lncRNAs associated with cancer metabolism and tumor immune cell infiltration could predict the survival prognosis of patients, and AP005131.2, AL121832.2, and LINC01152 could be novel biomarkers and therapeutic targets for BC.

Introduction: Solute carrier family 31 member 1(SLC31A1) has been reported as the copper importer, and was identified to be involved in the process of “cuproptosis”. However, the mechanism of SLC31A1 in breast cancer remains unclear.Methods: We examined the expression of SLC31A1 mRNA in breast cancer tissues and cell lines using Real-time PCR. The data for this study were obtained from The Cancer Genome Atlas (TCGA) database and analyzed via R 3.6.3. TIMER, UALCAN, GEPIA2, STRING, Metascape, Kaplan–Meier Plotter, starBase and miRNet websites were used for a comprehensive analysis of SLC31A1.Results: Our study suggested that SLC31A1 mRNA was over-expressed in breast tumor tissue and breast cancer cell lines, and which was closely related to poor relapse-free survival (RFS) and distant metastasis-free survival (DMFS). In addition, we constructed a co-expression network of SLC31A1. Functional enrichment analysis indicated that they were mainly involved in copper ion transport. Interestingly, SLC31A1 expression was positively associated with all m6A-related genes, especially with YTHDF3 (r = 0.479). Importantly, the LINC00511/miR-29c-3p/SLC31A1 axis was identified as the most potential pathway promoting breast cancer progress by affecting copper transport. Furthermore, the expression level of SLC31A1 in breast cancer was positively correlated with tumor immune cell infiltration, immune cell biomarkers and cancer-associated fibroblast (CAF).Conclusion: Up-regulation of SLC31A1 expression and regulation of copper ion transport mediated by LINC00511-miR-29-3p axis is related to poor prognosis and positively correlated with tumor immune infiltration in breast cancer.