MiR-206 Target Prediction in Breast Cancer Subtypes by Bioinformatics Tools

authors:

avatar Mahnaz Seifi-Alan 1 , avatar Roshanak Shamsi 1 , avatar Ali Behmanesh 2 , avatar Reza Mirfakhraie ORCID 1 , avatar Mir Davood Omrani 1 , avatar Soudeh Ghafouri-Fard ORCID 1 , *

Department of Medical Genetics, Shahid Beheshti University of Medical Sciences, Tehran, Iran
Iran University of Medical Sciences, Tehran, Iran

how to cite: Seifi-Alan M, Shamsi R, Behmanesh A, Mirfakhraie R, Omrani M D, et al. MiR-206 Target Prediction in Breast Cancer Subtypes by Bioinformatics Tools. Int J Cancer Manag. 2018;11(7):e69554. https://doi.org/10.5812/ijcm.69554.

Abstract

Background:

MicroRNAs (miRNAs) are small endogenous non-coding RNAs with fundamental roles in the regulation of protein expression that is involved in the pathogenesis of many cancers including breast cancer. Among them is miR-206, whose role as a tumor suppressor gene has been demonstrated in breast cancer. Consequently, the identification of its putative target in breast cancer is of practical value.

Methods:

In the present study, we have suggested a new approach for the identification of miR-206 target genes with possible role in breast cancer pathogenesis. We used 15 online tools for the prediction of miR-206 target genes as well as gene expression data produced by DNA microarray technology.

Results:

By combining these two sets of data, we suggested a list of miR-206 target genes with possible involvement in breast cancer. In addition, we depicted an interaction network including miR-206 and its putative targets.

Conclusions:

Considering the complexity of miR-206 interactions with several targets, such in silico analyses would considerably lessen the work load of laboratory experiments.

1. Background

As the most frequently diagnosed type of cancer, breast cancer among women can be classified based on expression patterns to luminal A, luminal B, Her2+, and triple negative (TN) subtypes, which are correlated with patients’ survival and prognosis (1). Several researches have focused on the exploration of gene expression patterns among these subtypes to find a biologically relevant biomarker in breast cancer (2-5). The expression pattern analysis of non-coding RNAs has also been the focus of researchers (6, 7). Pathway-based expression analysis has further been suggested as a more systematic strategy for biomarker discovery (8, 9). MicroRNAs (miRNAs), as the major regulators of gene expression, and pathway regulation are thought to participate in tumorigenesis process by changing the expression of several mRNA coding genes, non-coding RNA as well as cancer-related signaling pathways. MiR-206 is among miRNAs, whose role has been highlighted as a tumor suppressor. It has been the first miRNA detected in breast cancer. A miR-206-binding site has been found within the 3’-untranslated region (3’UTR) of estrogen receptor (ER)-α. ER negative breast cancer cells such as MDA-MB-231 cells have higher expression of this miRNA compared with ER positive cells such as MCF-7 cells. Notably, the forced over-expression of this miRNA in triple negative breast cancer cells has diminished their metastatic potential and decreased the expression levels of matrix metalloproteinase while increased expression of breast cancer metastatic suppressor (BRMS)-1. Functional studies have shown Cx43 to be a target of miR-206. Expression analysis in human breast cancer specimens has demonstrated the associations between miR-206 levels and both lymph node status and Cx43 expression. While miR-206 significantly decreased the proliferation and metastatic potential of cancer cells, it could not suppress tumor initiation in a mouse xenograft model due to the contraindicatory roles of Cx43 during tumorigenesis (10). Its down-regulation in breast cancer tissues has been associated with increased tumor size and advanced tumor stage. Up-regulation of miR-206 in MCF-7 breast cancer cells prevented cell proliferation and cell growth by the inhibition of the G1/S transition. MiR-206 has an inhibitory effect on the expression of cyclinD2 at both the transcript and protein levels. Taken together, miR-206 has been recognized as a critical tumor suppressor gene in breast tissues, a probable prognostic biomarker, or a beneficial target for the treatment of patients with breast cancer (11).

Based on the complex interaction network between miRNAs and mRNAs, which is implicated in physiologic processes and cancer, several studies aimed at predicting the miRNA targets. However, the valid prediction of miRNA targets is still problematic. Merging experimental and computational tools for decoding miRNA functions and targets has been suggested as an attractive approach in this regard (12). Consequently, in the present study, we have suggested a bioinformatics approach for the identification of miRNA-206 target genes implicated in breast cancer based on miRNA target prediction tools and microarray data analysis. This miRNA has been chosen based on its evident role in the suppression of breast cancer and its contribution in the regulation of tumor-associated pathways.

2. Methods

The workflow of data analysis is demonstrated in Figure 1.

The bioinformatics pipeline used in the current study
The bioinformatics pipeline used in the current study

2.1. Breast Cancer Expression Data

mRNA expression profiles for 285 breast cancer samples and 22 normal tissues were collected from gene expression omnibus (GEO) repository as series GSE65194 and GSE45827. The GEO database launched in 2000 by the national center for biotechnology information (NCBI) incorporates gene expression data provided from microarray technology (13). By entering “breast cancer” and “transcriptome” key words and choosing “Expression profiling by array” as the filter and the minimum sample size of 100 for datasets, 10 datasets were retrieved. GSE45827 and GSE65194 datasets with similar array platforms (GPL570) and inclusion of 4 molecular subtypes of breast cancer (luminal A, luminal B, Her2+ and TN) have been selected for further analyses. The pipeline applied for the selection of these datasets is demonstrated in Figure 2.

Flow chart of the protocol used for the search of breast cancer microarray datasets from the GEO database
Flow chart of the protocol used for the search of breast cancer microarray datasets from the GEO database

2.2. Prediction of miRNA-206 Targets from miRNA Prediction Databases

In this study, systematic miRNA-206 target search was carried out on 14 prediction databases (Mirwalk, MirMAP, Mirbase, PITA, Microrna.org conserved, RNA22, MBStar, DIANA, MirTarbase, mirDIP, MirTar2, Mirza G, PACCMIT, Targetscan, TargetSpy) (Table 1).

Table 1.

Summary of miRNA Prediction Tools Used in the Present Study

Target Prediction ToolOrganismFeatures of ToolWebsite
DIANAAnyConservation, seed match, and free energyhttp://www.microrna.gr/microT-CDS
miRWalkHuman, mouse, and ratConservation, seed match and free energyhttp://mirwalk.uni-hd.de/
Target scanHuman, mouse, Fly, Fish, and WormConservation, seed matchhttp://www.targetscan.org/
PITAHuman, mouse, Fly, and WormConservation, seed match, free energy, and site accessibilityhttp://genie.weizmann.ac.il/pubs/mir07/
ExiqonHuman, mouse, and ratConservation, seed match, and Target-Site Abundancehttps://www.exiqon.com/miRSearch
RNA22Human, Fruit Fly, Mouse, and WormSeed match and free energyhttps://cm.jefferson.edu/rna22/
TargetSpyHuman, mouse, rat, Fruit Fly, and ChickenSeed match and site accessibilityhttps://omictools.com/targetspy-tool
MirMAPHuman, Chimpanzee, Mouse, Rat, Cow, Chicken, Zebrafish, and OpossumConservation, seed match, and free energyhttp://mirmap.ezlab.org
MirbaseHuman, mouse, rat, worm, and flySeed match and Target-Site Abundancehttp://microrna.sanger.ac.uk
MBStarHumanSeed match and free energyhttp://www.isical.ac.in/~bioinfo_miu/MBStar30
Mirza-GHumanConservation, seed match, and site accessibilityhttps://omictools.com/mirza-genome-wide-tool
PACCMITAnyConservation, seed match, and site accessibilityhttp://lcpt.epfl.ch
MirTar2Human, mouse, rat, dog, and chickenConservation, seed match, free energy, and site accessibilityhttp://mirdb.org
Microrna. orgHuman, mouse, Fruit Fly, and ratConservation, seed match, free energyhttp://www.microrna.org

The fundamental criteria used in some of these prediction databases (for instance TargetScan and DIANA) are complementarity between miRNA and mRNA in seed regions, folding free energy for the miRNA-mRNA duplex, and evolutionary conservation (14). In order to increase the sensitivity of our approach, we chose the mRNAs predicted as miR-206 target by at least 2 databases out of the 15 databases. A list of candidate mRNAs was selected.

2.3. Experimentally Validated miR-206 Target Genes

Due to the drawbacks associated with the predictions of miRNA target sites with the bioinformatics tools (15), in order to find functionally relevant targets in breast cancer, we searched the literature to find functional studies for the identification of miR-206 targets in breast cancer as well as all types of cancer.

2.4. Detection of Differentially Expressed Genes in Breast Cancer Subtypes

We compared the expression profile of 4 breast cancer subtype samples with normal samples by GEO2R web tool (https://www.ncbi.nlm.nih.gov/geo/info/geo2r.html). Based on Log2-fold change between 2 experimental conditions (LogFC) and adjusted P values calculated by the software, 300 genes have been chosen in each subtype with the highest differential expression among tumoral and normal samples.

2.5. Identification of miR-206 Targets Among Differentially Expressed Genes in Breast Cancer Subtypes

The R statistical program (16) was applied to find miR-206 target genes, which are differentially expressed in each cancer subtype compared with normal samples.

2.6. Enrichment Annotation Analysis and Network Construction

FunRich analysis tool version 3 (17) was used for the network construction and the enrichment annotation step. This software depicts a schematic clusterization of the gene list with pathway annotations.

3. Results

3.1. Experimentally Validated miR-206 Target Genes

With the purpose of assessment of the sensitivity of miRNA target prediction tools, we compared the list of experimentally validated miR-206 target genes with those predicted by computational algorithms. The lists of experimentally validated miR-206 target genes in all cancer types and breast cancer as well as the number of bioinformatics tools predicted each gene are presented in Tables 2 and 3, respectively. CCND1 and CCND2 were predicted as targets of miR-206 by a high proportion of tools and were also experimentally validated to be targets of this miRNA in diverse cancer types.

Table 2.

Experimentally Validated Target Genes of miR-206 in All Cancers

Target GeneNumber of Tools Predicted the Gene as miR-206 TargetValidation TechniqueDiseaseReference
CCND210 toolsReporter assayGastric cancer(18)
MET8 toolsReporter assay, Western blot, qPCRRhabdomyosarcoma(19)
GPD28 toolsReporter assay, Western blotMany tumor cell lines(20)
G6PD8 toolsReporter assay, Western blotMany tumor cell lines(20)
CCND16 toolsReporter assay, Western blotDifferent cancer cell lines(21)
CDK45 toolsReporter assayMelanoma(22)
PGD5 toolsReporter assay, Western blotMany tumor cell lines(20)
TKT3 toolsReporter assay, Western blotMany tumor cell lines(20)
EGFR3 toolsReporter assay, Western blot, qPCRLung Squamous Cell Carcinoma(23)
Table 3.

Experimentally Validated Target Genes of miR-206 in Breast Cancer

Target GeneNumber of Tools Predicted the Gene as miR-206 TargetValidation TechniqueReference
ESR110 toolsReporter assay, Western blot, qPCR(24)
CCND210 toolsReporter assay(11)
TWF110 toolsReporter assay, qPCR and Western blot(25)
SPRED110 toolsReporter assay(26)
VEGF8 toolsReporter assay(27)
RASA18 toolsReporter assay(26)
NRP18 toolsReporter assay(28)
SMAD28 toolsReporter assay(28)
MAP3K136 toolsReporter assay, qPCR(29)
CCND16 toolsReporter assay(30)
Ncl (Nucleolin)6 toolsReporter assay, qPCR and Western blot(31)
Tbx36 toolsReporter assay(32)
PDCD44 toolsReporter assay(33)
CORO1C4 toolsReporter assay, qPCR and Western blot(34)
PFKFB33 toolsReporter assay, qPCR and Western blot(35)
Connexin43TargetScan and PicTarReporter assay and qPCR(10)

3.2. miR-206 Target Genes Implicated in Breast Cancer

Based on the proposed approach of combining microarray data analysis and bioinformatics prediction tools, we identified miR-206 target genes, which are differentially expressed in breast cancer tissues compared with normal tissues and are possibly implicated in breast cancer (Table 4). The high degrees of similarity in expression profile were observed for SH3GL3, KANK3, and CWF19L2 (overexpressed in all subtypes) and TPM3 (down-regulated in all subtypes). However, some other genes had a subtype specific expression profile. For instance, WDR48 and PTPRS were down-regulated and up-regulated only in luminal B subtype respectively.

Table 4.

MiR-206 Predicted Target Genes with Differential Expression in Breast Cancer Subtypes Based on the Log2-Transformed Fold Changea

Gene SymbolGene DescriptionHER2Triple NegativeLuminal ALuminal BNumber of Tools Predicted the Gene as miR-206 Target
SH3GL3SH3 domain containing GRB2 like 3, endophilin A3****6 tools
PNPPurine nucleoside phosphorylase###6 tools
CORO1CCoronin 1C#6 tools
WDR48WD repeat domain 48#6 tools
PTPRSProtein tyrosine phosphatase, receptor type S*6 tools
SULF1Sulfates 1###5 tools
ATP6V1AATPase H+ transporting V1 subunit A###5 tools
YWHAZTyrosine 3-monooxygenase activation protein zeta##5 tools
FN1Fibronectin 1###5 tools
KIF2AKinesin family member 2A#5 tools
RAB5AMember RAS oncogene family#5 tools
UBE2HUbiquitin conjugating enzyme E2 H###4 tools
TPM3Tropomyosin 3####4 tools
KDELR2KDEL endoplasmic reticulum protein retention receptor2##4 tools
SFRP1Secreted frizzled related protein 1*4 tools
TWF1Twinfilin actin binding protein 1#4 tools
FPGTFucose-1-phosphate guanylyltransferase#4 tools
SRSF1Serine and arginine rich splicing factor 1#4 tools
GMFBGlia maturation factor beta###3 tools
KANK3KN motif and ankyrin repeat domains 3****3 tools
CALUCalcium-binding protein***3 tools
HOXA5Homeobox A5**3 tools
VMP1Vacuole membrane protein 1#3 tools
TRIM59Tripartite motif containing 59###3 tools
SRPK2SRSF protein kinase 2##3 tools
PGM5Phosphoglucomutase 5***3 tools
MAPK1Mitogen-activated protein kinase 1###3 tools
ACTR3ARP3 actin related protein 3 homolog###3 tools
CWF19L2CWF19 like 2, cell cycle control****3 tools
NCBP1Nuclear cap binding protein subunit 1#3 tools
SRSF3Serine and arginine rich splicing factor 3##3 tools
INMTIndolethylamine N-methyltransferase**3 tools
SPTBN1Spectrin beta, non-erythrocytic 1*3 tools
ANP32EAcidic nuclear phosphoprotein 32 family member E#3 tools
SDPRCaveolae associated protein 2*3 tools
AMTAminomethyltransferase*3 tools
TBC1D9TBC1 domain family member 9#3 tools
DDX3XDEAD-box helicase 3, X-linked#3 tools
PREX1Phosphatidylinositol-3,4,5-trisphosphate dependent Rac exchange factor 1#3 tools
PCDH19Protocadherin 19*3 tools
ITIH5Inter-alpha-trypsin inhibitor heavy chain family 5**3 tools
MAPRE1Microtubule associated protein RP/EB family member 1#3 tools
PRPF4BPre-mRNA processing factor 4B###3 tools
GABRPGamma-aminobutyric acid type A receptor pi subunit*3 tools
LMOD1Leiomodin 1**3 tools
POLR3KRNA polymerase III subunit K##3 tools
UCHL5Ubiquitin C-terminal hydrolase L5#3 tools
CA12Carbonic anhydrase 12#3 tools
EDEM3ER degradation enhancing alpha-mannosidase like pro3#3 tools
WISP1WNT1 inducible signaling pathway protein 1#3 tools
RBM47RNA binding motif protein 47#3 tools
FBXO22F-box protein 22#3 tools
ARIH1Ariadne RBR E3 ubiquitin protein ligase 1#3 tools
ZNF146Zinc finger protein 146#3 tools
ARF4ADP ribosylation factor 4#3 tools
IKBKBInhibitor of nuclear factor kappa B kinase subunit beta#3 tools
RAB18RAB18, member RAS oncogene family#3 tools
NIPBLNIPBL, cohesin loading factor*3 tools
SS18nBAF chromatin remodeling complex subunit#3 tools
ZNF326Zinc finger protein 326#3 tools

3.3. Network Construction

FunRich tool provided both functional enrichment and interaction network for miR-206 and associated mRNAs. The final dataset obtained from the combinatory approach was analyzed against 3 different background databases, namely FunRich, UniProt, and Custom. Figure 3 demonstrates the interaction diagram. Interaction network provided by this tool shows that miR-206 is involved in fundamental pathways in breast cancer. The most enriched one was TGF-beta pathway.

The schematic interaction diagram of miR-206-target interactions in breast cancer
The schematic interaction diagram of miR-206-target interactions in breast cancer

4. Discussion

The expression pattern analysis of distinct subtypes of breast cancer might be helpful in the identification of specific biomarkers as well as the development of personalized treatment modalities (36, 37). The role of miRNAs and their putative targets might be different in each cancer type or even among cancer subtypes. In the present study, we proposed a computational method for the identification of miR-206 targets in distinct breast cancer subtypes. Our suggested approach might guide researchers to arrange their hypotheses of miRNA function on a certain group of microarray samples to choose the candidate mRNAs and miRNAs prior to laboratory experiment. Such approach would also help in the identification of the role of predicted targets as tumor suppressor or oncogenes in signaling pathways. In the final list provided by the proposed approach, there are some genes with expression change in all breast cancer subtypes as well as some genes with alterations in a certain subtype. Such data would help researchers to choose the more suitable cell line for the validation of miR-206 targets.

MiR-206 targets have been validated in experimental studies. For instance, in ER positive MCF-7 cell line, luciferase assays showed that miR-206 decreases cyclinD2 expression by targeting 2 binding sites in the 3’-UTR of cyclinD2 mRNA. The results have been confirmed by qRT-PCR and Western blot assays (11). Another study in the same cell line has shown miR-201 inhibitory effect on the expression of transforming growth factor (TGF)-β, neuropilin-1 (NRP1), and SMAD2, which participate in migration, invasion, and EMT in these cells (28). In TN breast cancer cell line, miR-206 up-regulation has decreased the expression of matrix metalloproteinase (MMP)-2, MMP-9 and increased the expression level of breast cancer metastatic suppressor (BRMS)-1. However, luciferase assays showed only GJA1 (Cx43) as a target of miR-206 (10). Considering the complex interaction network between certain miRNA and numerous mRNAs as well as among various mRNAs, the identification of the net target of each miRNA is of practical significance. Currently, no high-throughput and low-cost miRNA target screening method is available. Although a number of computational methods based on sequence complementarity of the miRNA and the mRNAs have been proposed, the predicted interactions using these computational techniques are unreliable according to high false positive rates. Incorporation of the expression values of miRNAs and mRNAs (and/or proteins) has been suggested as a method to improve the results of sequence-based predictions (12). Our proposed approach provides a tool for combination of miRNA target prediction with mRNA expression data. Although this method is not claimed to be perfect, the suggested list of targets produced by this method is anticipated to significantly lessen laboratory experimental load and the number of interactions to be validated. Homogeneity in experimental conditions should be maintained in our proposed approach. Considering the more significant role of other factors such as transcription factors in regulation of gene expression compared with miRNAs role, heterogeneity of information sources is regarded as a limitation in such approaches (12).

In conclusion, our analysis demonstrated novel mRNA candidates as probable miR-206 targets, which possibly participate in the pathogenesis of breast cancer. The gene list provided by our proposed approach contains some genes previously validated as miR-206 targets in breast cancer (CORO1C, TWF1). However, many of the experimentally validated miR-206 targets are not among gene list provided by our approach. The main cause of such discrepancy is that in our proposed approach, we included genes with high differential expression between tumor tissues and normal samples to find the most biologically relevant targets in breast cancer. However, in the experimental validation of miRNAs targets, the level of expression change is much lower than our defined level in differential expression analyses. For instance, in Yin et al.’s study, the enforced expression of miR-206 resulted in about 50% reduction in the expression of NRP1 and SMAD2 (28). Fu et al. (10) observed an approximately similar reduction in Cx43 levels following the over-expression of miR-206. However, based on our proposed approach, we considered 300 genes with the most differentially expression between tumoral and normal samples, which had at least 10 fold change. Such threshold definition is anticipated to cover the most biologically relevant mRNA targets. A possible limitation of our approach is that by using microarray expression data, miRNA targets can only be detected if the mRNA is degraded by the miRNA. Therefore, our approach may leave out some targets that are influenced at the translation level.

Acknowledgements

References

  • 1.

    Ghafouri-Fard S, Shamsi R, Seifi-Alan M, Javaheri M, Tabarestani S. Cancer-testis genes as candidates for immunotherapy in breast cancer. Immunotherapy. 2014;6(2):165-79. [PubMed ID: 24491090]. https://doi.org/10.2217/imt.13.165.

  • 2.

    Dianatpour M, Mehdipour P, Nayernia K, Mobasheri MB, Ghafouri-Fard S, Savad S, et al. Expression of Testis Specific Genes TSGA10, TEX101 and ODF3 in Breast Cancer. Iran Red Crescent Med J. 2012;14(11):722-6. [PubMed ID: 23396665]. [PubMed Central ID: PMC3560543]. https://doi.org/10.5812/ircmj.3611.

  • 3.

    Iranpour M, Soudyab M, Geranpayeh L, Mirfakhraie R, Azargashb E, Movafagh A, et al. Expression analysis of four long noncoding RNAs in breast cancer. Tumour Biol. 2016;37(3):2933-40. [PubMed ID: 26409453]. https://doi.org/10.1007/s13277-015-4135-2.

  • 4.

    Kazemi-Oula G, Ghafouri-Fard S, Mobasheri MB, Geranpayeh L, Modarressi MH. Upregulation of RHOXF2 and ODF4 Expression in Breast Cancer Tissues. Cell J. 2015;17(3):471-7. [PubMed ID: 26464818]. [PubMed Central ID: PMC4601867].

  • 5.

    Sarrafzadeh S, Geranpayeh L, Tasharrofi B, Soudyab M, Nikpayam E, Iranpour M, et al. Expression Study and Clinical Correlations of MYC and CCAT2 in Breast Cancer Patients. Iran Biomed J. 2017;21(5):303-11. [PubMed ID: 28480695]. [PubMed Central ID: PMC5548962].

  • 6.

    Tasharrofi B, Soudyab M, Nikpayam E, Iranpour M, Mirfakhraie R, Sarrafzadeh S, et al. Comparative expression analysis of hypoxia-inducible factor-alpha and its natural occurring antisense in breast cancer tissues and adjacent noncancerous tissues. Cell Biochem Funct. 2016;34(8):572-8. [PubMed ID: 27862063]. https://doi.org/10.1002/cbf.3230.

  • 7.

    Soudyab M, Iranpour M, Ghafouri-Fard S. The Role of Long Non-Coding RNAs in Breast Cancer. Arch Iran Med. 2016;19(7):508-17. [PubMed ID: 27362246].

  • 8.

    Taherian-Esfahani Z, Abedin-Do A, Nouri Z, Mirfakhraie R, Ghafouri-Fard S, Motevaseli E. Lactobacilli Differentially Modulate mTOR and Wnt/ beta-Catenin Pathways in Different Cancer Cell Lines. Iran J Cancer Prev. 2016;9(3). e5369. [PubMed ID: 27703648]. [PubMed Central ID: PMC5038836]. https://doi.org/10.17795/ijcp-5369.

  • 9.

    Esfandiary A, Taherian-Esfahani Z, Abedin-Do A, Mirfakhraie R, Shirzad M, Ghafouri-Fard S, et al. Lactobacilli Modulate Hypoxia-Inducible Factor (HIF)-1 Regulatory Pathway in Triple Negative Breast Cancer Cell Line. Cell J. 2016;18(2):237-44. [PubMed ID: 27540529]. [PubMed Central ID: PMC4988423].

  • 10.

    Fu Y, Shao ZM, He QZ, Jiang BQ, Wu Y, Zhuang ZG. Hsa-miR-206 represses the proliferation and invasion of breast cancer cells by targeting Cx43. Eur Rev Med Pharmacol Sci. 2015;19(11):2091-104. [PubMed ID: 26125274].

  • 11.

    Zhou J, Tian Y, Li J, Lu B, Sun M, Zou Y, et al. miR-206 is down-regulated in breast cancer and inhibits cell proliferation through the up-regulation of cyclinD2. Biochem Biophys Res Commun. 2013;433(2):207-12. [PubMed ID: 23466356]. https://doi.org/10.1016/j.bbrc.2013.02.084.

  • 12.

    Muniategui A, Pey J, Planes FJ, Rubio A. Joint analysis of miRNA and mRNA expression data. Brief Bioinform. 2013;14(3):263-78. [PubMed ID: 22692086]. https://doi.org/10.1093/bib/bbs028.

  • 13.

    Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30(1):207-10. [PubMed ID: 11752295]. [PubMed Central ID: PMC99122].

  • 14.

    Li L, Xu J, Yang D, Tan X, Wang H. Computational approaches for microRNA studies: a review. Mamm Genome. 2010;21(1-2):1-12. [PubMed ID: 20012966]. https://doi.org/10.1007/s00335-009-9241-2.

  • 15.

    Reyes-Herrera PH, Ficarra E. One decade of development and evolution of microRNA target prediction algorithms. Genomics Proteomics Bioinformatics. 2012;10(5):254-63. [PubMed ID: 23200135]. [PubMed Central ID: PMC5054202]. https://doi.org/10.1016/j.gpb.2012.10.001.

  • 16.

    Team RC. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2014.

  • 17.

    Pathan M, Keerthikumar S, Ang CS, Gangoda L, Quek CY, Williamson NA, et al. FunRich: An open access standalone functional enrichment and interaction network analysis tool. Proteomics. 2015;15(15):2597-601. [PubMed ID: 25921073]. https://doi.org/10.1002/pmic.201400515.

  • 18.

    Zhang L, Liu X, Jin H, Guo X, Xia L, Chen Z, et al. miR-206 inhibits gastric cancer proliferation in part by repressing cyclinD2. Cancer Lett. 2013;332(1):94-101. [PubMed ID: 23348698]. https://doi.org/10.1016/j.canlet.2013.01.023.

  • 19.

    Yan D, Dong Xda E, Chen X, Wang L, Lu C, Wang J, et al. MicroRNA-1/206 targets c-Met and inhibits rhabdomyosarcoma development. J Biol Chem. 2009;284(43):29596-604. [PubMed ID: 19710019]. [PubMed Central ID: PMC2785592]. https://doi.org/10.1074/jbc.M109.020511.

  • 20.

    Singh A, Happel C, Manna SK, Acquaah-Mensah G, Carrerero J, Kumar S, et al. Transcription factor NRF2 regulates miR-1 and miR-206 to drive tumorigenesis. J Clin Invest. 2013;123(7):2921-34. [PubMed ID: 23921124]. [PubMed Central ID: PMC3696551]. https://doi.org/10.1172/JCI66353.

  • 21.

    Alteri A, De Vito F, Messina G, Pompili M, Calconi A, Visca P, et al. Cyclin D1 is a major target of miR-206 in cell differentiation and transformation. Cell Cycle. 2013;12(24):3781-90. [PubMed ID: 24107628]. [PubMed Central ID: PMC3905070]. https://doi.org/10.4161/cc.26674.

  • 22.

    Georgantas R3, Streicher K, Luo X, Greenlees L, Zhu W, Liu Z, et al. MicroRNA-206 induces G1 arrest in melanoma by inhibition of CDK4 and Cyclin D. Pigment Cell Melanoma Res. 2014;27(2):275-86. [PubMed ID: 24289491]. https://doi.org/10.1111/pcmr.12200.

  • 23.

    Mataki H, Seki N, Chiyomaru T, Enokida H, Goto Y, Kumamoto T, et al. Tumor-suppressive microRNA-206 as a dual inhibitor of MET and EGFR oncogenic signaling in lung squamous cell carcinoma. Int J Oncol. 2015;46(3):1039-50. [PubMed ID: 25522678]. https://doi.org/10.3892/ijo.2014.2802.

  • 24.

    Adams BD, Furneaux H, White BA. The micro-ribonucleic acid (miRNA) miR-206 targets the human estrogen receptor-alpha (ERalpha) and represses ERalpha messenger RNA and protein expression in breast cancer cell lines. Mol Endocrinol. 2007;21(5):1132-47. [PubMed ID: 17312270]. https://doi.org/10.1210/me.2007-0022.

  • 25.

    Samaeekia R, Adorno-Cruz V, Bockhorn J, Chang YF, Huang S, Prat A, et al. miR-206 Inhibits Stemness and Metastasis of Breast Cancer by Targeting MKL1/IL11 Pathway. Clin Cancer Res. 2017;23(4):1091-103. [PubMed ID: 27435395]. [PubMed Central ID: PMC5247402]. https://doi.org/10.1158/1078-0432.CCR-16-0943.

  • 26.

    Sharma SB, Lin CC, Farrugia MK, McLaughlin SL, Ellis EJ, Brundage KM, et al. MicroRNAs 206 and 21 cooperate to promote RAS-extracellular signal-regulated kinase signaling by suppressing the translation of RASA1 and SPRED1. Mol Cell Biol. 2014;34(22):4143-64. [PubMed ID: 25202123]. [PubMed Central ID: PMC4248710]. https://doi.org/10.1128/MCB.00480-14.

  • 27.

    Liang Z, Bian X, Shim H. Downregulation of microRNA-206 promotes invasion and angiogenesis of triple negative breast cancer. Biochem Biophys Res Commun. 2016;477(3):461-6. [PubMed ID: 27318091]. [PubMed Central ID: PMC4955785]. https://doi.org/10.1016/j.bbrc.2016.06.076.

  • 28.

    Yin K, Yin W, Wang Y, Zhou L, Liu Y, Yang G, et al. MiR-206 suppresses epithelial mesenchymal transition by targeting TGF-beta signaling in estrogen receptor positive breast cancer cells. Oncotarget. 2016;7(17):24537-48. [PubMed ID: 27014911]. [PubMed Central ID: PMC5029720]. https://doi.org/10.18632/oncotarget.8233.

  • 29.

    Han H, Chen Y, Cheng L, Prochownik EV, Li Y. microRNA-206 impairs c-Myc-driven cancer in a synthetic lethal manner by directly inhibiting MAP3K13. Oncotarget. 2016;7(13):16409-19. [PubMed ID: 26918941]. [PubMed Central ID: PMC4941324]. https://doi.org/10.18632/oncotarget.7653.

  • 30.

    Elliman SJ, Howley BV, Mehta DS, Fearnhead HO, Kemp DM, Barkley LR. Selective repression of the oncogene cyclin D1 by the tumor suppressor miR-206 in cancers. Oncogenesis. 2014;3. e113. [PubMed ID: 25111862]. [PubMed Central ID: PMC5189965]. https://doi.org/10.1038/oncsis.2014.26.

  • 31.

    Bose S, Tholanikunnel TE, Reuben A, Tholanikunnel BG, Spicer EK. Regulation of nucleolin expression by miR-194, miR-206, and HuR. Mol Cell Biochem. 2016;417(1-2):141-53. [PubMed ID: 27221739]. https://doi.org/10.1007/s11010-016-2721-2.

  • 32.

    Amir S, Simion C, Umeh-Garcia M, Krig S, Moss T, Carraway K3, et al. Regulation of the T-box transcription factor Tbx3 by the tumour suppressor microRNA-206 in breast cancer. Br J Cancer. 2016;114(10):1125-34. [PubMed ID: 27100732]. [PubMed Central ID: PMC4865973]. https://doi.org/10.1038/bjc.2016.73.

  • 33.

    Lin CC, Sharma SB, Farrugia MK, McLaughlin SL, Ice RJ, Loskutov YV, et al. Kruppel-like factor 4 signals through microRNA-206 to promote tumor initiation and cell survival. Oncogenesis. 2015;4. e155. [PubMed ID: 26053033]. [PubMed Central ID: PMC4753526]. https://doi.org/10.1038/oncsis.2015.8.

  • 34.

    Wang J, Tsouko E, Jonsson P, Bergh J, Hartman J, Aydogdu E, et al. miR-206 inhibits cell migration through direct targeting of the actin-binding protein coronin 1C in triple-negative breast cancer. Mol Oncol. 2014;8(8):1690-702. [PubMed ID: 25074552]. [PubMed Central ID: PMC5528580]. https://doi.org/10.1016/j.molonc.2014.07.006.

  • 35.

    Ge X, Lyu P, Cao Z, Li J, Guo G, Xia W, et al. Overexpression of miR-206 suppresses glycolysis, proliferation and migration in breast cancer cells via PFKFB3 targeting. Biochem Biophys Res Commun. 2015;463(4):1115-21. [PubMed ID: 26093295]. https://doi.org/10.1016/j.bbrc.2015.06.068.

  • 36.

    Esfandiary A, Ghafouri-Fard S. MAGE-A3: an immunogenic target used in clinical practice. Immunotherapy. 2015;7(6):683-704. [PubMed ID: 26100270]. https://doi.org/10.2217/imt.15.29.

  • 37.

    Esfandiary A, Ghafouri-Fard S. New York esophageal squamous cell carcinoma-1 and cancer immunotherapy. Immunotherapy. 2015;7(4):411-39. [PubMed ID: 25917631]. https://doi.org/10.2217/imt.15.3.