Identification of mRNA and non-coding RNA hubs using network analysis in organ tropism regulated triple negative breast cancer metastasis

Abstract

Triple negative breast cancer (TNBC) is aggressive in nature, resistant to conventional therapy and often ends in organ specific metastasis. In this study, publicly available datasets were used to identify miRNA, mRNA and lncRNA hubs. Using validated mRNA-miRNA, mRNA-mRNA and lncRNA-miRNA interaction information obtained from various databases, RNA interaction networks for TNBC and its subtype specific as well as organ tropism regulated metastasis were generated. Further, miRNA-mRNA-lncRNA triad classification was performed using social network analysis from subnetworks and visualized using Cytoscape. Survival analysis of the RNA hubs, oncoprint analysis for mRNAs and pathway analysis of the lncRNAs were also performed. Results indicated that two lncRNAs (NEAT1 and CASC7) and four miRNAs (hsa-miR-106b-5p, hsa-miR-148a-3p, hsa-miR-25-3p and hsa-let-7i-5p) were common between hubs identified in TNBC and TNBC associated metastasis. The exclusive hubs for TNBC associated metastasis were hsa-miR-200b-3p, SP1, HSPA4 and RAB1B. HMGA1 was the top ranked hub in mesenchymal subtype associated lung metastasis, while hsa-miR-27a-3p was identified as the top ranked hub mRNA in luminal androgen receptor subtype associated bone metastasis. When lncRNA associated pathway analysis was performed, Hs Cytoplasmic Ribosomal Protein pathway was found to be the most significant and among the selected hubs, CTNND1, SON and hsa-miR-29c emerged as TNBC survival markers. TP53, FOXA1, MTDH and HDGF were found as the top ranked mRNAs in oncoprint analysis. The pipeline proposed for the first time in this study with validated RNA interaction data integration and graph-based learning for miRNA-mRNA-lncRNA triad classification from RNA hubs may aid experimental cost reduction and its successful execution will allow it to be extended to other diseases too.