Global Analysis of miRNA–mRNA Interaction Network in Breast Cancer with Brain Metastasis

Background: MicroRNAs (miRNAs) have been linked to a number of cancer types including breast cancer. The rate of brain metastases is 10-30% in patients with advanced breast cancer which is associated with poor prognosis. The potential application of miRNAs in the diagnostics and therapeutics of breast cancer with brain metastasis is an area of intense interest. In an initial effort to systematically address the differential expression of miRNAs and mRNAs in primary breast cancer which may provide clues for early detection of brain metastasis, we analyzed the consequent changes in global patterns of gene expression in Gene Expression Omnibus (GEO) data set obtained by microarray from patients with in situ carcinoma and patients with brain metastasis. Materials and Methods: The miRNA-pathway regulatory network and miRNA–mRNA regulatory network were investigated in breast cancer specimens from patients with brain metastasis to screen for significantly dysregulated miRNAs followed by prediction of their target genes and pathways by Gene Ontology (GO) analysis. Results: Functional coordination of the changes of gene expression can be modulated by individual miRNAs. Two miRNAs, hsa-miR-17-5p and hsa-miR-16-5p, were identified as having the highest associations with targeted mRNAs [such as B-cell lymphoma 2 (BCL2), small body size/mothers against decapentaplegic 3 (SMAD3) and suppressor of cytokine signaling 1 (SOCS1)] and pathways associated with epithelial–mesenchymal transitions and other processes linked with cancer metastasis (including cell cycle, adherence junctions and extracellular matrix–receptor interaction). mRNAs for two genes [HECT, UBA and WWE domain containing 1 (HUWE1) and BCL2] were found to have the highest associations with miRNAs, which were down-regulated in brain metastasis specimens of breast cancer. The change of 11 selected miRNAs was verified in The Cancer Genome Atlas (TCGA) breast cancer dataset. Up-regulation of hsa-miR-17-5p was detected in triple-negative breast cancer tissues in TCGA. Furthermore, a negative correlation of hsa-miR-17-5p with overall survival and phosphatase and tensin homolog (PTEN) and BCL2 target genes was found in TCGA breast cancer specimens. Conclusion: Our findings provide a functionally coordinated expression pattern of different families of miRNAs that may have potential to provide clinicians with a strategy to treat breast cancer with brain metastasis from a systems-rather than a single-gene perspective.