Transcriptome of bone marrow-Derived stem cells reveals new inflammatory mediators related to increased survival in patients with multiple myeloma.

Although multiple myeloma (MM) is a neoplasm that leads affected individuals to death, little is known about why some patients survive much longer than others. In this context, we investigated the transcriptomic profile of bone marrow hematopoietic stem cells obtained from MM patients and compared the clinical outcomes of death and survival six months after bone marrow transplantation. The leukapheresis products of 39 patients with MM eligible for autologous transplantation were collected and analyzed. After extraction, the RNA was analyzed using the GeneChip Human Exon 1.0 Array method. The transcriptome profile was analyzed in silico, and the differentially expressed signaling pathways of interest were validated. The results showed a difference in the expression of inflammation-related genes, immune response processes, and the oxidative stress pathway. The in silico study also pointed out the involvement of the NFκB transcription factor in the possible modulation of these genes. We chose to validate molecules participating in these processes, including the cytokines TNF-α, IFN-γ, and TGF-ß1; in addition, we measured the levels of oxidative stress mediators (pro-oxidant profile and the total antioxidant capacity). TNF-α levels were significantly reduced in patients who died and were over 50 years old at diagnosis, as well as in patients with plasmacytoma. Increased TNF-α was detected in patients with very high levels of ß2-microglobulin. IFN-γ reduction was observed in patients with a complete response to treatment compared to those with a very good response. Patients with plasmacytoma who died also had an increased pro-oxidant profile. These data show the profile of inflammatory response markers that are altered in patients with MM who die quickly and serve as a basis for the development of future studies of markers to predict better survival in this disease.

Epigenetic Alterations in DCIS Progression: What Can lncRNAs Teach Us?

Some transcripts that are not translated into proteins can be encoded by the mammalian genome. Long noncoding RNAs (lncRNAs) are noncoding RNAs that can function as decoys, scaffolds, and enhancer RNAs and can regulate other molecules, including microRNAs. Therefore, it is essential that we obtain a better understanding of the regulatory mechanisms of lncRNAs. In cancer, lncRNAs function through several mechanisms, including important biological pathways, and the abnormal expression of lncRNAs contributes to breast cancer (BC) initiation and progression. BC is the most common type of cancer among women worldwide and has a high mortality rate. Genetic and epigenetic alterations that can be regulated by lncRNAs may be related to early events of BC progression. Ductal carcinoma in situ (DCIS) is a noninvasive BC that is considered an important preinvasive BC early event because it can progress to invasive BC. Therefore, the identification of predictive biomarkers of DCIS-invasive BC progression has become increasingly important in an attempt to optimize the treatment and quality of life of patients. In this context, this review will address the current knowledge about the role of lncRNAs in DCIS and their potential contribution to the progression of DCIS to invasive BC.

Genome-wide analysis reveals a rhamnolipid-dependent modulation of flagellar genes in Pseudomonas aeruginosa PAO1.

Pseudomonas aeruginosa is an opportunistic pathogen and an important model organism for the study of bacterial group behaviors, including cell motility and biofilm formation. Rhamnolipids play a pivotal role in biofilm formation and motility phenotypes in P. aeruginosa, possibly acting as wetting agents and mediating chemotactic stimuli. However, no biochemical mechanism or gene regulatory network has been investigated in regard to rhamnolipids' modulation of those group behaviors. Using DNA microarrays, we investigated the transcriptomic profiles in the stationary phase of growth of wild-type P. aeruginosa PAO1 and a rhlA-mutant strain, unable to produce rhamnolipids. A total of 134 genes were differentially expressed, comprising different functional categories, indicating a significant physiological difference between the rhamnolipid-producing and -non-producing strains. Interestingly, several flagellar genes are repressed in the mutant strain, which directly relates to the inability of the rhlA-minus strain to develop a swarming-motility phenotype. Supplementation with exogenous rhamnolipids has partially restored flagellar gene expression in the mutant strain. Our results show significant evidence that rhamnolipids, the major biosynthetic products of rhlABC pathway, seem to modulate gene expression in P. aeruginosa.

KMT2A-MLLT1 and the Novel SEC16A-KMT2A in a Cryptic 3-Way Translocation t(9;11;19) Present in an Infant With Acute Lymphoblastic Leukemia.

About 25% of the patients with the translocation t(11;19)(q23;p13.3)/KMT2A-MLLT1 present three-way or more complex fusions, associated with a worse prognosis, suggesting that a particular mechanism creates functional KMT2A fusions for this condition. In this work, we show a cryptic three-way translocation t(9;11;19). Interestingly, long-distance inverse polymerase chain reaction sequencing revealed a KMT2A-MLLT1 and the yet unreported out-of-frame SEC16A-KMT2A fusion, associated with low SEC16A expression and KMT2A overexpression, in an infant with B-acute lymphoblastic leukemia presenting a poor prognosis. Our case illustrates the importance of molecular cytogenetic tests in selecting cases for further investigations, which could open perspectives regarding novel therapeutic approaches for poor prognosis childhood leukemias.

Can torque teno virus be a predictor of SARS-CoV-2 disease progression in cancer patients?

INTRODUCTION: Cancer patients with SARS-CoV-2 infection can experience a broad range of clinical manifestations and outcomes. Previous studies have demonstrated an association between torque teno virus (TTV) load and deficiencies of the immune system. The impact of SARS-CoV-2 and TTV viral loads in cancer patients is unknown. METHODS: In this retrospective study, 157 cancer patients and 191 noncancer controls were analysed for SARS-CoV-2 RNA and TTV DNA presence. RESULTS: SARS-CoV-2 RNA was detected in 66.2% of cancer patients and in 68.6% of noncancer control subjects. In SARS-CoV-2-positive patients, TTV was detectable in 79.8% of cancer patients, while in controls, TTV was detected in 71.7% of subjects. No statistically significant correlation was found between TTV and SARS-CoV-2 loads in cancer patients. However, the 100-day survival rate in cancer patients who died from COVID-19 was significantly lower in the TTV-positive group than in the TTV-negative group (P = 0.0475). In the cancer TTV-positive group, those who died also had a higher load of TTV than those who did not die (P = 0.0097). CONCLUSIONS: Our findings indicated that the presence of TTV in nasopharyngeal swabs from cancer patients was related to a higher number of deaths from COVID-19 and to a higher TTV DNA load.

Translational Results of Zo-NAnTax: A Phase II Trial of Neoadjuvant Zoledronic Acid in HER2-Positive Breast Cancer.

Breast cancer is a heterogeneous disease with distinct clinical and molecular characteristics. Scientific advances in molecular subtype differentiation support the understanding of cellular signaling, crosstalk, proliferation, survival, migration, and invasion mechanisms, allowing the development of new molecular drug targets. The breast cancer subtype with super expression and/or amplification of human growth factor receptor 2 (HER2) is clinically aggressive, but prognosis significantly shifted with the advent of anti-HER2 targeted therapy. Zoledronic-acid (ZOL) combined with a neoadjuvant Trastuzumab-containing chemotherapy regimen (Doxorubicin, Cyclophosphamide followed by Docetaxel, Trastuzumab) increased the pCR rate in a RH-positive/ HER2-positive subgroup, according to the phase II Zo-NAnTax trial. To verify genes that could be related to this response, a microarray assay was performed finding 164 differentially expressed genes. Silico analysis of these genes showed signaling pathways related to growth factors, apoptosis, invasion, and metabolism, as well as differentially expressed genes related to estrogen response. In addition, the RAC3 gene was found to interact with the MVD gene, a member of the mevalonate pathway. Taken together, these results indicate that RH-positive/ HER2-positive patients present gene alterations before treatment, and these could be related to the improvement of pCR.

Twist1 Influences the Expression of Leading Members of the IL-17 Signaling Pathway in HER2-Positive Breast Cancer Cells.

Breast cancer (BC) is a heterogeneous disease composed of multiple subtypes with different molecular characteristics and clinical outcomes. The metastatic process in BC depends on the transcription factors (TFs) related to epithelial-mesenchymal transition (EMT), including the master regulator Twist1. However, its role beyond EMT in BC subtypes remains unclear. Our study aimed to investigate the role of Twist1, beyond EMT, in the molecular subtypes of BC. In patients, we observed the overexpression of TWIST1 in the HER2+ group. The silencing of TWIST1 in HER2+ BC cells resulted in the upregulation of 138 genes and the downregulation of 174 genes compared to control cells in a microarray assay. In silico analysis revealed correlations between Twist1 and important biological processes such as the Th17-mediated immune response, suggesting that Twist1 could be relevant for IL-17 signaling in HER2+ BC. IL-17 signaling was then examined, and it was shown that TWIST1 knockdown caused the downregulation of leading members of IL-17 signaling pathway. Taken together, our findings suggest that Twist1 plays a role on IL-17 signaling in HER2+ BC.

Label-free proteomic analysis of breast cancer molecular subtypes.

To better characterize the cellular pathways involved in breast cancer molecular subtypes, we performed a proteomic study using a label-free LC-MS strategy for determining the proteomic profile of Luminal A, Luminal-HER2, HER2-positive, and triple-negative (TN) breast tumors compared with healthy mammary tissue. This comparison aimed to identify the aberrant processes specific for each subtype and might help to refine our understanding regarding breast cancer biology. Our results address important molecular features (both specific and commonly shared) that explain the biological behavior of each subtype. Changes in proteins related to cytoskeletal organization were found in all tumor subtypes, indicating that breast tumors are under constant structural modifications to invade and metastasize. We also found changes in cell-adhesion processes in all molecular subtypes, corroborating that invasiveness is a common property of breast cancer cells. Luminal-HER2 and HER2 tumors also presented altered cell cycle regulation, as shown by the several DNA repair-related proteins. An altered immune response was also found as a common process in the Luminal A, Luminal-HER2, and TN subtypes, and complement was the most important pathway. Analysis of the TN subtype revealed blood coagulation as the most relevant biological process.

A novel HLA-C*02:10 variant, HLA-C*02:10:09, identified in a healthy individual from Brazil.

Identification of the novel HLA-C*02:10:09 allele that differs from HLA-C*02:10:01:01 at one position in exon 1.

Implementing Standard Diagnosis and Treatment for Locally Advanced Breast Cancer Through Global Research in Latin America: Results From a Multicountry Pragmatic Trial.

PURPOSE: Breast cancer mortality rates in Latin America (LA) are higher than those in the United States, possibly because of advanced disease presentation, health care disparities, or unfavorable molecular subtypes. The Latin American Cancer Research Network was established to address these challenges and to promote collaborative clinical research. The Molecular Profiling of Breast Cancer Study (MPBCS) aimed to evaluate the clinical characteristics and treatment outcomes of LA participants with locally advanced breast cancer (LABC). PATIENTS AND METHODS: The MPBCS enrolled 1,449 participants from Argentina, Brazil, Chile, Mexico, and Uruguay. Through harmonized procedures and quality assurance measures, this study evaluated clinicopathologic characteristics, neoadjuvant chemotherapy response, and survival outcomes according to residual cancer burden (RCB) and the type of surgery. RESULTS: Overall, 711 and 480 participants in the primary surgery and neoadjuvant arms, respectively, completed the 5-year follow-up period. Overall survival was independently associated with RCB (worse survival for RCBIII-adjusted hazard ratio, 8.19, P < .001, and RCBII [adjusted hazard ratio, 3.69, P < .008] compared with RCB0 [pathologic complete response or pCR]) and type of surgery (worse survival in mastectomy than in breast-conserving surgery [BCS], adjusted hazard ratio, 2.97, P = .001). The hormone receptor-negative-human epidermal growth factor receptor 2-positive group had the highest proportion of pCR (48.9%). The analysis of the ASCO Quality Oncology Practice Initiative breast module revealed high compliance with pathologic standards but lower adherence to treatment administration standards. Notably, compliance with trastuzumab administration varied widely among countries (33.3%-88.7%). CONCLUSION: In LABC, we demonstrated the survival benefit of BCS and the prognostic effect of the response to available neoadjuvant treatments despite an important variability in access to key treatments. The MPBCS represents a significant step forward in understanding the real-world implementation of oncologic procedures in LA.

Iroquois Family Genes in Gastric Carcinogenesis: A Comprehensive Review.

Gastric cancer (GC) is the fifth leading cause of cancer-associated death worldwide, accounting for 768,793 related deaths and 1,089,103 new cases in 2020. Despite diagnostic advances, GC is often detected in late stages. Through a systematic literature search, this study focuses on the associations between the Iroquois gene family and GC. Accumulating evidence indicates that Iroquois genes are involved in the regulation of various physiological and pathological processes, including cancer. To date, information about Iroquois genes in GC is very limited. In recent years, the expression and function of Iroquois genes examined in different models have suggested that they play important roles in cell and cancer biology, since they were identified to be related to important signaling pathways, such as wingless, hedgehog, mitogen-activated proteins, fibroblast growth factor, TGFß, and the PI3K/Akt and NF-kB pathways. In cancer, depending on the tumor, Iroquois genes can act as oncogenes or tumor suppressor genes. However, in GC, they seem to mostly act as tumor suppressor genes and can be regulated by several mechanisms, including methylation, microRNAs and important GC-related pathogens. In this review, we provide an up-to-date review of the current knowledge regarding Iroquois family genes in GC.

Integrated Network Analysis of microRNAs, mRNAs, and Proteins Reveals the Regulatory Interaction between hsa-mir-200b and CFL2 Associated with Advanced Stage and Poor Prognosis in Patients with Intestinal Gastric Cancer.

Intestinal gastric cancer (IGC) carcinogenesis results from a complex interplay between environmental and molecular factors, ultimately contributing to disease development. We used integrative bioinformatic analysis to investigate IGC high-throughput molecular data to uncover interactions among differentially expressed genes, microRNAs, and proteins and their roles in IGC. An integrated network was generated based on experimentally validated microRNA-gene/protein interaction data, with three regulatory circuits involved in a complex network contributing to IGC progression. Key regulators were determined, including 23 microRNA and 15 gene/protein hubs. The regulatory circuit networks were associated with hallmarks of cancer, e.g., cell death, apoptosis and the cell cycle, the immune response, and epithelial-to-mesenchymal transition, indicating that different mechanisms of gene regulation impact similar biological functions. Altered expression of hubs was related to the clinicopathological characteristics of IGC patients and showed good performance in discriminating tumors from adjacent nontumor tissues and in relation to T stage and overall survival (OS). Interestingly, expression of upregulated hub hsa-mir-200b and its downregulated target hub gene/protein CFL2 were related not only to pathological T staging and OS but also to changes during IGC carcinogenesis. Our study suggests that regulation of CFL2 by hsa-miR-200b is a dynamic process during tumor progression and that this control plays essential roles in IGC development. Overall, the results indicate that this regulatory interaction is an important component in IGC pathogenesis. Also, we identified a novel molecular interplay between microRNAs, proteins, and genes associated with IGC in a complex biological network and the hubs closely related to IGC carcinogenesis as potential biomarkers.

High Expression of THY1 in Intestinal Gastric Cancer as a Key Factor in Tumor Biology: A Poor Prognosis-Independent Marker Related to the Epithelial-Mesenchymal Transition Profile.

Gastric cancer (GC) is an important cancer-related death worldwide. Among its histological subtypes, intestinal gastric cancer (IGC) is the most common. A previous work showed that increased expression of the THY1 gene was associated with poor overall survival in IGC. Furthermore, it was shown that IGC tumor cells with high expression of THY1 have a greater capacity for tumorigenesis and metastasis in vitro. This study aimed to identify molecular differences between IGC with high and low expression of THY1. Using a feature selection method, a group of 35 genes were found to be the most informative gene set for THY1high IGC tumors. Through a classification model, these genes differentiate THY1high from THY1low tumors with 100% of accuracy both in the test subset and the independent test set. Additionally, this group of 35 genes correctly clustered 100% of the samples. An extensive validation of this potential molecular signature in multiple cohorts successfully segregated between THY1high and THY1low IGC tumors (>95%), proving to be independent of the gene expression quantification methodology. These genes are involved in central processes to tumor biology, such as the epithelial-mesenchymal transition (EMT) and remodeling of the tumor tissue composition. Moreover, patients with THY1high IGC demonstrated poor survival and a more advanced clinicopathological staging. Our findings revealed a molecular signature for IGC with high THY1 expression. This signature showed EMT and remodeling of the tumor tissue composition potentially related to the biology of IGC. Altogether, our results indicate that THY1high IGC tumors are a particular subset of tumors with a specific molecular and prognosis profile.

Wnt/ß-catenin pathway regulates ABCB1 transcription in chronic myeloid leukemia.

BACKGROUND: The advanced phases of chronic myeloid leukemia (CML) are known to be more resistant to therapy. This resistance has been associated with the overexpression of ABCB1, which gives rise to the multidrug resistance (MDR) phenomenon. MDR is characterized by resistance to nonrelated drugs, and P-glycoprotein (encoded by ABCB1) has been implicated as the major cause of its emergence. Wnt signaling has been demonstrated to be important in several aspects of CML. Recently, Wnt signaling was linked to ABCB1 regulation through its canonical pathway, which is mediated by ß-catenin, in other types of cancer. In this study, we investigated the involvement of the Wnt/ß-catenin pathway in the regulation of ABCB1 transcription in CML, as the basal promoter of ABCB1 has several ß-catenin binding sites. ß-catenin is the mediator of canonical Wnt signaling, which is important for CML progression. METHODS: In this work we used the K562 cell line and its derived MDR-resistant cell line Lucena (K562/VCR) as CML study models. Real time PCR (RT-qPCR), electrophoretic mobility shift assay (EMSA), chromatin immunoprecipitation (ChIP), flow cytometry (FACS), western blot, immunofluorescence, RNA knockdown (siRNA) and Luciferase reporter approaches were used. RESULTS: ß-catenin was present in the protein complex on the basal promoter of ABCB1 in both cell lines in vitro, but its binding was more pronounced in the resistant cell line in vivo. Lucena cells also exhibited higher ß-catenin levels compared to its parental cell line. Wnt1 and ß-catenin depletion and overexpression of nuclear ß-catenin, together with TCF binding sites activation demonstrated that ABCB1 is positively regulated by the canonical pathway of Wnt signaling. CONCLUSIONS: These results suggest, for the first time, that the Wnt/ß-catenin pathway regulates ABCB1 in CML.

The protein-protein interaction network of intestinal gastric cancer patients reveals hub proteins with potential prognostic value.

BACKGROUND: Gastric cancer (GC) is the third leading cause of cancer worldwide. According to the Lauren classification, gastric adenocarcinoma is divided into two subtypes: diffuse and intestinal. The development of intestinal gastric cancer (IGC) can take years and involves multiple factors. OBJECTIVE: To investigate the protein profile of tumor samples from patients with IGC in comparison with adjacent nontumor tissue samples. METHODS: We used label-free nano-LC-MS/MS to identify proteins from the tissues samples. The results were analyzed using MetaCore™ software to access functional enrichment information. Protein-protein interactions (PPI) were predicted using STRING analysis. Hub proteins were determined using the Cytoscape plugin, CytoHubba. Survival analysis was performed using KM plotter. We identified 429 differentially expressed proteins whose pathways and processes were related to protein folding, apoptosis, and immune response. RESULTS: The PPI network of these proteins showed enrichment modules related to the regulation of cell death, immune system, neutrophil degranulation, metabolism of RNA and chromatin DNA binding. From the PPI network, we identified 20 differentially expressed hub proteins, and assessed the prognostic value of the expression of genes that encode them. Among them, the expression of four hub genes was significantly associated with the overall survival of IGC patients. CONCLUSIONS: This study reveals important findings that affect IGC development based on specific biological alterations in IGC patients. Bioinformatics analysis showed that the pathogenesis of IGC patients is complex and involves different interconnected biological processes. These findings may be useful in research on new targets to develop novel therapies to improve the overall survival of patients with IGC.

Reduced Osteogenic Differentiation Potential In Vivo in Acute Myeloid Leukaemia Patients Correlates with Decreased BMP4 Expression in Mesenchymal Stromal Cells.

The osteogenic differentiation potential of mesenchymal stromal cells (hMSCs) is an essential process for the haematopoiesis and the maintenance of haematopoietic stem cells (HSCs). Therefore, the aim of this work was to evaluate this potential in hMSCs from AML patients (hMSCs-AML) and whether it is associated with BMP4 expression. The results showed that bone formation potential in vivo was reduced in hMSCs-AML compared to hMSCs from healthy donors (hMSCs-HD). Moreover, the fact that hMSCs-AML were not able to develop supportive haematopoietic cells or to differentiate into osteocytes suggests possible changes in the bone marrow microenvironment. Furthermore, the expression of BMP4 was decreased, indicating a lack of gene expression committed to the osteogenic lineage. Overall, these alterations could be associated with changes in the maintenance of HSCs, the leukaemic transformation process and the development of AML.

The Transcriptomic Portrait of Locally Advanced Breast Cancer and Its Prognostic Value in a Multi-Country Cohort of Latin American Patients.

Purposes: Most molecular-based published studies on breast cancer do not adequately represent the unique and diverse genetic admixture of the Latin American population. Searching for similarities and differences in molecular pathways associated with these tumors and evaluating its impact on prognosis may help to select better therapeutic approaches. Patients and Methods: We collected clinical, pathological, and transcriptomic data of a multi-country Latin American cohort of 1,071 stage II-III breast cancer patients of the Molecular Profile of Breast Cancer Study (MPBCS) cohort. The 5-year prognostic ability of intrinsic (transcriptomic-based) PAM50 and immunohistochemical classifications, both at the cancer-specific (OSC) and disease-free survival (DFS) stages, was compared. Pathway analyses (GSEA, GSVA and MetaCore) were performed to explore differences among intrinsic subtypes. Results: PAM50 classification of the MPBCS cohort defined 42·6% of tumors as LumA, 21·3% as LumB, 13·3% as HER2E and 16·6% as Basal. Both OSC and DFS for LumA tumors were significantly better than for other subtypes, while Basal tumors had the worst prognosis. While the prognostic power of traditional subtypes calculated with hormone receptors (HR), HER2 and Ki67 determinations showed an acceptable performance, PAM50-derived risk of recurrence best discriminated low, intermediate and high-risk groups. Transcriptomic pathway analysis showed high proliferation (i.e. cell cycle control and DNA damage repair) associated with LumB, HER2E and Basal tumors, and a strong dependency on the estrogen pathway for LumA. Terms related to both innate and adaptive immune responses were seen predominantly upregulated in Basal tumors, and, to a lesser extent, in HER2E, with respect to LumA and B tumors. Conclusions: This is the first study that assesses molecular features at the transcriptomic level in a multicountry Latin American breast cancer patient cohort. Hormone-related and proliferation pathways that predominate in PAM50 and other breast cancer molecular classifications are also the main tumor-driving mechanisms in this cohort and have prognostic power. The immune-related features seen in the most aggressive subtypes may pave the way for therapeutic approaches not yet disseminated in Latin America. Clinical Trial Registration: ClinicalTrials.gov (Identifier: NCT02326857).

SUZ12 is a candidate target of the non-canonical WNT pathway in the progression of chronic myeloid leukemia.

Polycomb proteins form multiprotein complexes that repress target genes by chromatin remodeling. In this work, we report that the SUZ12 polycomb gene is over-expressed in bone marrow samples of patients at the blastic phase of chronic myeloid leukemia. We also found a direct interaction between polycomb group genes and the WNT signaling pathway in chronic myeloid leukemia transformation. Electrophoretic mobility shift assay (EMSA), Chromatin immunoprecipitation assay (ChIP), and mass spectrometry assays identified noncanonical WNT pathway members, such as WNT5A and WNT11, bound to the SUZ12 promoter. Immunohistochemistry and immunofluorescence with WNT5A and WNT11 antibodies confirmed nuclear localization. Knockdown of WNTs 1, 5A, and 11 with RNAi approaches showed that WNT members are capable of activating SUZ12 transcription with varying promoter affinities. Finally, we suggest that SUZ12 is blocking cellular differentiation, as SUZ12 knockdown release differentiation programs in chronic myeloid blastic phase (CML-BP) transformed cell line.

miRNome Profiling Reveals Shared Features in Breast Cancer Subtypes and Highlights miRNAs That Potentially Regulate MYB and EZH2 Expression.

Breast cancer (BC) has been extensively studied, as it is one of the more commonly diagnosed cancer types worldwide. The study of miRNAs has increased what is known about the complexity of pathways and signaling and has identified potential biomarkers and therapeutic targets. Thus, miRNome profiling could provide important information regarding the molecular mechanisms involved in BC. On average, more than 430 miRNAs were identified as differentially expressed between BC cell lines and normal breast HMEC cells. From these, 110 miRNAs were common to BC subtypes. The miRNome enrichment analysis and interaction maps highlighted epigenetic-related pathways shared by all BC cell lines and revealed potential miRNA targets. Quantitative evaluation of BC patient samples and GETx/TCGA-BRCA datasets confirmed MYB and EZH2 as potential targets from BC miRNome. Moreover, overall survival was impacted by EZH2 expression. The expression of 15 miRNAs, selected according to aggressiveness of BC subtypes, was confirmed in TCGA-BRCA dataset. Of these miRNAs, miRNA-mRNA interaction prediction revealed 7 novel or underexplored miRNAs in BC: miR-1271-5p, miR-130a-5p, and miR-134 as MYB regulators and miR-138-5p, miR-455-3p, miR-487a, and miR-487b as EZH2 regulators. Herein, we report a novel molecular miRNA signature for BC and identify potential miRNA/mRNAs involved in disease subtypes.