Detection of miRNA regulatory effect on triple negative breast cancer transcriptome
- Loredana Martignetti†1, 2, 3Email author,
- Bruno Tesson†1, 2, 3, 4,
- Anna Almeida1, 4,
- Andrei Zinovyev1, 2, 3,
- Gordon C Tucker6,
- Thierry Dubois1, 4, 5 and
- Emmanuel Barillot1, 2, 3
© Martignetti et al.; licensee BioMed Central Ltd. 2015
Published: 1 June 2015
Identifying key microRNAs (miRNAs) contributing to the genesis and development of a particular disease is a focus of many recent studies. We introduce here a rank-based algorithm to detect miRNA regulatory activity in cancer-derived tissue samples which combines measurements of gene and miRNA expression levels and sequence-based target predictions. The method is designed to detect modest but coordinated changes in the expression of sequence-based predicted target genes. We applied our algorithm to a cohort of 129 tumour and healthy breast tissues and showed its effectiveness in identifying functional miRNAs possibly involved in the disease. These observations have been validated using an independent publicly available breast cancer dataset from The Cancer Genome Atlas. We focused on the triple negative breast cancer subtype to highlight potentially relevant miRNAs in this tumour subtype. For those miRNAs identified as potential regulators, we characterize the function of affected target genes by enrichment analysis. In the two independent datasets, the affected targets are not necessarily the same, but display similar enriched categories, including breast cancer related processes like cell substrate adherens junction, regulation of cell migration, nuclear pore complex and integrin pathway. The R script implementing our method together with the datasets used in the study can be downloaded here (http://bioinfo-out.curie.fr/projects/targetrunningsum).
MicroRNAs (miRNAs) are endogenous ~22 nucleotide RNA molecules that act as fundamental repressors of gene expression in many biological systems. In animals, they target mRNAs by recognizing and directly binding to multiple partially complementary sites preferentially located in the 3' untranslated regions (UTRs) of transcripts. Watson-Crick base-pairing to the 5' end of miRNAs, especially to the so-called 'seed' region that comprises nucleotides 2-7, is considered crucial for targeting [1, 2], even if recently developed techniques for ligation and sequencing of miRNA-target RNA duplexes highlight widespread noncanonical seed interactions, containing bulged or mismatched nucleotides . Although molecular mechanisms of miRNA action remain intensely debated , multiple studies revealed that mammalian miRNAs repress genes predominantly by destabilization of target mRNAs [5, 6]. By computational and experimental approaches it was established that thousands of human protein-coding genes are regulated by miRNAs [7, 8]. Given the wide scope of their targeting, miRNAs are considered as an additional layer of regulatory circuitry in the cell. Experimental observations suggest that miRNAs are regulators of development and cellular homeostasis through their control of diverse biological processes, from differentiation and proliferation to apoptosis . Their role in regulating fundamental cell mechanisms suggests that they could be involved in cancer and indeed their expression is strongly deregulated in almost all human malignancies. Functional characterization of these aberrantly expressed microRNAs indicates that they might function as oncogenes and tumor suppressors [10, 11].
Identifying key miRNAs contributing to the genesis and development of a particular disease is a focus of many recent studies. Statistical methodology for this task is not fully established due to the mild effect of miRNAs on the expression of their targets. A major source of information to infer the actual regulatory activity of miRNAs derives from high-throughput experimental data such as transcriptome profiles. The underlying assumption is that regulatory activity by miRNAs could be reflected by the expression changes of their target transcripts. For miRNAs that promote mRNA decay, there would be a negative correlation between miRNA and mRNA expression. Existing tools based on this assumption mainly rely on case-control mRNA profile experiments involving strong perturbations such as the knockout/knockdown/overexpression of one or few miRNAs [12–14]. In data originating from less controlled conditions, such as mRNA profiles of pathological tissue collected from patients, detecting miRNA-mediated target destabilization is more challenging due to presence of multiple cell types in samples, the activity of additional regulatory factors and complex RNA cross-regulation such as the miRNA sponge effect [15–17]. The recent study of a large cohort of breast tumour samples  suggests that miRNAs exert their effect by modulating mRNAs rather than acting as on-off switches. Other studies inferring miRNA regulation on tumour sample transcriptome exploit additional molecular information such as AGO2-PAR-CLIP binding-site data  or DNA copy number and promoter methylation at the mRNA gene locus .
We applied our method to elucidate the regulatory effect of miRNAs on the breast cancer transcriptome. Breast cancer is classified into various subtypes mainly based on the immunohistochemical staining of estrogen (ER), progesterone (PR) and HER-2 (ERBB2) receptors. The complex nature and heterogeneity of this disease, particularly with regard to gene expression profiles, make it difficult to detect the shaping effect of miRNAs on the transcriptome. We applied our algorithm to a breast carcinoma dataset including gene and miRNA expression from normal and breast tumour samples (which we refer to as Maire dataset [22, 23]) and we show that it is able to identify miRNAs with statistically significant enrichment score. These results are then compared with those obtained on an independent dataset of normal and breast tumour samples from The Cancer Genome Atlas (TCGA) . We further focus on the triple negative breast cancer (TNBC) to highlight miRNAs potentially relevant in this particular tumour subgroup. This subtype is intensively studied due to the lack of effective targeted therapies. We ran our algorithm including only samples characterized as triple-negative breast tumours and identified a set of miRNAs showing statistically significant signal in both Maire and TCGA datasets.
Finally, we investigated miRNAs identified as potential regulators to characterize the function of their targets. In the proposed algorithm, hundreds of genes account for the enrichment signal of a single miRNA and we expect a subset of them to participate in common cellular functions. We use multiple annotation databases to infer biological processes affected by the identified miRNAs. For those identified in both datasets, the corresponding sets of target genes were analyzed separately. We observe that even if the specific genes accounting for the enrichment of biological categories among miRNA target genes detected in the Maire dataset and in the TCGA datasets are not necessarily the same, they are associated to common cancer-related pathways
Functional miRNAs in breast cancer
The analyzed dataset includes 129 tumour and healthy breast tissues for which both miRNA and mRNA expression profiles are available (see Material and Methods). In this study, we used sequence-based target sets obtained from TargetScan version 6.2 [8, 25], a widely used algorithm which takes into account sequences that match the seed region of each miRNA and evaluate their conservation in several species. The confidence of target predictions is calculated as described in Methods. We include in our study 394 miRNAs for which both expression data and TargetScan predictions are available.
We identified 136 miRNAs as potential regulators with FDR < 0.1 (results are reported in Additional File 1, Table S1). Among the top significant miRNAs, we found several ones that are known to function as oncomirs, such as the members of miR 17-92 cluster and its paralog cluster miR-106b-25, miR-15 and miR-16 .
Common predicted miRNAs in Maire cohort and TCGA for the analysis including all samples.
PV Adj MAIRE
PV Adj TCGA
Triple negative breast cancer specific study
Common predicted miRNAs in Maire cohort and TCGA for the TNBC specific analysis.
PV Adj MAIRE
PV Adj TCGA
A group of 20 miRNAs is identified specifically in the TNBC study. Remarkably, some of these are already associated with aggressive breast cancer: miR-29a, miR-29c and miR-148a have been shown to be downregulated and associated to aberrant hypermethylation in basal-like cell line , miR-27a is involved in endothelial differentiation of breast cancer in a basal-like cell line  and in the MDA-MB-231 basal-like cell line  and miR-139-5p is described as a regulator of breast cancer cell motility and invasion .
Identification of biological processes targeted by miRNAs
To assess the biological relevance of miRNAs identified as potential regulators, we investigated whether the genes that account for the enrichment signal of a given miRNA participate in the same cellular process or signalling pathway. The assumption that some miRNAs downregulate a group of genes participating in the same pathway is supported by multiple experimental studies [36–38]. Based on this hypothesis, for each significant miRNA we tested the leading-edge subset of genes for functional enrichment using curated annotation databases, such as Gene Ontology , KEGG , BioCarta , Reactome  and ACSN .
For miRNAs identified in both Maire and TCGA datasets, the corresponding subsets of leading-edge targets were analyzed separately. Interestingly, these subsets of genes display highly significant overlap and similar enriched categories, supporting the relevance of miRNA regulatory role in breast cancer. We report in Additional File 3, Table S2 the complete list of enriched categories with Bonferroni corrected p-values below 10−2. The most enriched categories include breast cancer related processes like cell substrate adherens junction, regulation of cell migration, nuclear pore complex and integrin pathway.
High-throughput mRNA and miRNA expression data from large cohorts of normal and pathological tissue samples can be exploited to detect miRNA regulatory activity. The rank-based algorithm introduced here is able to detect miRNA-mediated target destabilization from normal and breast cancer expression profiles. Reproducible results were obtained in two independent datasets, providing a list of miRNAs potentially relevant in breast cancer. Moreover, the association of these miRNAs to cancer related processes is supported by functional enrichment of affected target genes.
The fact that a better overlap was obtained between Maire cohort and TCGA when restricting the analysis to the TNBC subtype compared to when using all samples may at first seem counterintuitive since using more samples should allow better power to detect correlations. However, this observation may be explained by the fact that when considering all breast cancer subtypes and healthy samples together, a larger part of the variation in the transcriptome data will arise from factors that are not directly linked to miRNA activity. For example, when putting together data from the luminal and non luminal subtypes, much of the variation will be associated with the status of the estrogen receptor pathway. Such variation can induce an important correlation structure in the data that may confuse the detection of the much subtler variation associated with miRNA regulation. The proposed algorithm can be considered a suitable tool to elucidate the regulatory role of miRNAs in physiological conditions.
All genes are ranked according to the correlation between gene expression and the expression of miRNA m
The enrichment score defined in Equation (1) is computed for the sequence-based target set S m associated to the miRNA m
The significance of the enrichment score ES(S m ) is evaluated by a permutation procedure estimating an empirical p-value PV(ES(S m ))
The ranking scheme and the enrichment score definition are described as follows.
where r mj is the correlation between the expression of gene j and the expression of miRNA m, < s m · |r m | α > is the average of prediction confidence-weighted correlations for the set S m . The parameter α controls the contribution of the correlation r mj such that correlation values can contribute linearly (α = 1) or non-linearly in the running sum. Similar to the original GSEA algorithm, in this study, we set α = 1. According to this equation, the running sum is incremented by value (s mj · |r mj | α − < s m · |r m | α >) when encountering a gene in S m and decreased by < s m · |r m | α > when not.
The statistical significance of the enrichment score ES(S m ) is assessed by a permutation based p-value. The enrichment score of the randomly shuffled list G m is computed for N p =1000 permutation rounds and the empirical null distribution of ES(S m ) is generated. An empirical p-value PV(ES(S m )) is estimated for the positive and the negative region of the distribution by the proportion of permutations that result in larger ES(S m ) than originally observed (or lower ES(S m ) for the negative region). Once PV(ES(S m )) for a miRNA m is below a fixed threshold, the regulatory activity of miRNA m on the transcriptome is inferred. The PV threshold is set according to False Dicovery Rate (FDR) obtained by the Benjamini-Hochberg method . In our study, the PV threshold was set according to FDR < 0.1.
Healthy samples from mammary plastic surgery and tumor samples were obtained from patients treated at the hospital of Institut Curie (Biological Resource Center, Paris, France) as described previously [22, 23]. Experiments were performed in agreement with the Bioethic Law No. 2004-800 and the Ethic Charter from the French National Institute of Cancer (INCa), and after approval of the ethics committee of our Institution. From the initial dataset, we retained only the samples for which both mRNA and miRNA data were available, numbering to 129, including 11 healthy breast tissue samples, 37 TNBC, 28 ER-/HER2+, 24 Luminal A and 29 Luminal B breast tumour samples as characterized by immunohistochemical staining.
MiRNA expression data
Samples were hybridized on the Agilent miRNA microarray kit (V3). One hundred ng of total RNAs was hybridized to the microarrays according to the manufacturer's instructions. Hybridized microarrays were scanned with a DNA microarray scanner (Agilent G2565BA) and features were extracted using the Agilent Feature Extraction image analysis tool with default protocols and settings. Data were first transformed using the reverse hyperbolic sine function and quantile normalized. The data were then corrected for a hybridization batch effect using a linear model including the hybridization series as a fixed effect. Next, probes with negative intensity values in 95% or more of the arrays were discarded, leaving 516 miRNAs for analysis. When technical replicates for a sample were present, they were subsequently averaged. Two samples displayed an outlier behavior evident from Principal Component Analysis (data not shown) and were discarded from this analysis. The dataset can be downloaded here at the following address: http://bioinfo-out.curie.fr/projects/targetrunningsum.
Gene expression data
For the protein-coding transcriptome, the data from Affymetrix U133plus2 arrays was processed as described . In summary, we used the brainarray HGU133Plus2-Hs-Entrez version 13 custom chipset definition file , data were then normalized with GC-RMA , technical batch artefacts were corrected using a linear model, and genes with noise-level expression in 95% or more arrays were filtered out leaving 11543 probesets each corresponding to a unique Entrez Identifier.
MiRNA and gene expression data from TCGA
We conducted our study on the publicly available data of tumour and healthy breast tissues from TCGA described in . We selected 500 tumors and 21 tumor-adjacent normal breast tissue samples for which both mRNA and miRNA data were available. Of the 500 tumors, 82 were assigned to the TNBC subtype according to ER, PR and HER2 negative status. In this dataset mRNA expression levels were determined by Agilent custom 244K whole genome microarrays and miRNA abundance was measured by Illumina sequencing technology. Level 3 released data contain quality controlled and processed data done by Broad Institute's TCGA workgroup with expression call for genes per samples. Gene level expression data were normalized by using Robust Multi-array Average (RMA) and expression values were gene centered. MiRNA expression was given as read counts normalized to relative read frequency in each sample. Detected miRNAs were defined as having more than 10 reads in at least 10% of the samples, leaving 332 miRNAs for analysis.
Target predictions by TargetScan
In the TargetScan algorithm, the prediction score of a miRNA binding site depends on the level of conservation and sequence context criteria such as the distance of the target from the 3'UTR end and the AU composition of the flanking area. For each miRNA, we take as prediction confidence weight of its targets the total context scores generated by the algorithm for 3'UTRs aggregating the binding site scores. By construction, the total context scores range between -1 and 0. When total context score for multiple 3'UTR isoforms are predicted, we take the total context score of the longest 3'UTR isoform.
The authors would like to thank Dr. Guillem Rigaill for suggesting the experimental design. The research leading to these results and the publication of this article have received funding from the European Union Seventh Framework Programme (FP7/2007-2013) ASSET project under grant agreement number FP7-HEALTH-2010-259348.
This article has been published as part of BMC Genomics Volume 16 Supplement 6, 2015: Proceedings of the Italian Society of Bioinformatics (BITS): Annual Meeting 2014: Genomics. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcgenomics/supplements/16/S6.
- Brennecke J, Stark A, Russell RB, Cohen SM: Principles of microRNA-target recognition. PLOS Biology. 2005, 3 (3): e85-10.1371/journal.pbio.0030085.PubMed CentralView ArticlePubMedGoogle Scholar
- Brodersen P, Voinnet O: Revisiting the principles of microRNA target recognition and mode of action. Nature Reviews Molecular Cell Biology. 2009, 10 (2): 141-148. 10.1038/nrm2619.View ArticlePubMedGoogle Scholar
- Helwak A, Kudla G, Dudnakova T, Tollervey D: Mapping the human miRNA interactome by CLASH reveals frequent noncanonical binding. Cell. 2013, 153 (3): 654-65. 10.1016/j.cell.2013.03.043.PubMed CentralView ArticlePubMedGoogle Scholar
- Morozova N, Zinovyev A, Nonne N, Pritchard LL, Gorban AN, Harel-Bellan A: Kinetic signatures of microRNA modes of action. RNA. 2012, 18 (9): 1635-55. 10.1261/rna.032284.112.PubMed CentralView ArticlePubMedGoogle Scholar
- Selbach M, Schwanh¨ausser B, Thierfelder N, Fang Z, Khanin R, Rajewsky N: Widespread changes in protein synthesis induced by microRNAs. Nature. 2008, 455 (7209): 58-63. 10.1038/nature07228.View ArticlePubMedGoogle Scholar
- Guo H, Ingolia NT, Weissman JS, Bartel DP: Mammalian microRNAs predominantly act to decrease target mRNA levels. Nature. 2010, 466 (7308): 835-840. 10.1038/nature09267.PubMed CentralView ArticlePubMedGoogle Scholar
- Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116 (2): 281-297. 10.1016/S0092-8674(04)00045-5.View ArticlePubMedGoogle Scholar
- Friedman RC, Farh KKH, Burge CB, Bartel DP: Most mammalian mRNAs are conserved targets of microRNAs. Genome Res. 2009, 19: 92-105.PubMed CentralView ArticlePubMedGoogle Scholar
- Kim Dh, Gru¨n D, van Oudenaarden A: Dampening of expression oscillations by synchronous regulation of a microRNA and its target. Nat Genet. 2013, 45 (11): 1337-44. 10.1038/ng.2763.View ArticlePubMedGoogle Scholar
- Schickel R, Boyerinas B, Park SM, Peter ME: MicroRNAs: key players in the immune system, differentiation, tumorigenesis and cell death. Oncogene. 2008, 27 (45): 5959-5974. 10.1038/onc.2008.274.View ArticlePubMedGoogle Scholar
- Esquela-Kerscher A, Slack FJ: Oncomirs - microRNAs with a role in cancer. Nat Rev Cancer. 2006, 6 (4): 259-269. 10.1038/nrc1840.View ArticlePubMedGoogle Scholar
- van Dongen S, Abreu-Goodger C, Enright AJ: Detecting microRNA binding and siRNA off-target effects from expression data. Nat Methods. 2008, 5 (12): 1023-5. 10.1038/nmeth.1267.PubMed CentralView ArticlePubMedGoogle Scholar
- Liang Z, Zhou H, He Z, Zheng H, Wu J: mirAct: a web tool for evaluating microRNA activity based on gene expression data. Nucleic Acids Res. 2011, 39: W139-44. 10.1093/nar/gkr351.PubMed CentralView ArticlePubMedGoogle Scholar
- Rasmussen SH, Jacobsen A, Krogh A: cWords - systematic microRNA regulatory motif discovery from mRNA expression data. Silence. 2013, 4 (1): 2-10.1186/1758-907X-4-2.PubMed CentralView ArticlePubMedGoogle Scholar
- Sumazin P, Yang X, Chiu HS, Chung WJ, Iyer A, Llobet-Navas D, Rajbhandari P, Bansal M, Guarnieri P, Silva J, Califano A: An extensive microRNA-mediated network of RNA-RNA interactions regulates established oncogenic pathways in glioblastoma. Cell. 2011, 147 (2): 370-81. 10.1016/j.cell.2011.09.041.PubMed CentralView ArticlePubMedGoogle Scholar
- Hansen TB, Kjems J, Damgaard CK: Circular RNA and miR-7 in cancer. Cancer Res. 2013, 73 (18): 5609-12. 10.1158/0008-5472.CAN-13-1568.View ArticlePubMedGoogle Scholar
- Tay Y, Rinn J, Pandolfi PP: The multilayered complexity of ceRNA crosstalk and competition. Nature. 2014, 505 (7483): 344-52. 10.1038/nature12986.PubMed CentralView ArticlePubMedGoogle Scholar
- Dvinge H, Git A, Graf S, Salmon-Divon M, Curtis C, Sottoriva A, Zhao Y, Hirst M, Armisen J, Miska EA, Chin SF, Provenzano E, Turashvili G, Green A, Ellis I, Aparicio S, Caldas C: The shaping and functional consequences of the microRNA landscape in breast cancer. Nature. 2013, 497 (7449): 378-82. 10.1038/nature12108.View ArticlePubMedGoogle Scholar
- Farazi TA, Ten Hoeve JJ, Brown M, Mihailovic A, Horlings HM, van de Vijver MJ, Tuschl T, Wessels LF: Identification of distinct miRNA target regulation between breast cancer molecular subtypes using AGO2-PAR-CLIP and patient datasets. Genome Biol. 2014, 7 (1): 15-R29.Google Scholar
- Jacobsen A1, Silber J, Harinath G, Huse JT, Schultz N, Sander C: Analysis of microRNA-target interactions across diverse cancer types. Nat Struct Mol Biol. 2013, 20 (11): 1325-32. 10.1038/nsmb.2678.PubMed CentralView ArticlePubMedGoogle Scholar
- Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005, 102 (43): 15545-50. 10.1073/pnas.0506580102.PubMed CentralView ArticlePubMedGoogle Scholar
- Maire V, Baldeyron C, Richardson M, Tesson B, Vincent-Salomon A, Gravier E, Marty-Prouvost B, De Koning L, Rigaill G, Dumont A, Gentien D, Barillot E, Roman-Roman S, Depil S, Cruzalegui F, Pierré A, Tucker GC, Dubois T: TTK/hMPS1 is an attractive therapeutic target for triple-negative breast cancer. PLoS One. 2013, 8 (5): e63712-10.1371/journal.pone.0063712.PubMed CentralView ArticlePubMedGoogle Scholar
- Maire V, Nemati F, Richardson M, Vincent-Salomon A, Tesson B, Rigaill G, Gravier E, Marty-Prouvost B, De Koning L, Lang G, Gentien D, Dumont A, Barillot E, Marangoni E, Decaudin D, Roman-Roman S, Pierré A, Cruzalegui F, Depil S, Tucker GC, Dubois T: Polo-like kinase 1: a potential therapeutic option in combination with conventional chemotherapy for the management of patients with triple-negative breast cancer. Cancer Res. 2013, 73 (2): 813-23. 10.1158/0008-5472.CAN-12-2633.View ArticlePubMedGoogle Scholar
- Cancer Genome Atlas Network: Comprehensive molecular portraits of human breast tumours. Nature. 2012, 490 (7418): 61-70. 10.1038/nature11412.View ArticleGoogle Scholar
- Grimson A, Farh KK, Johnston WK, Garrett-Engele P, Lim LP, Bartel DP: MicroRNA targeting specificity in mammals: determinants beyond seed pairing. Mol Cell. 2007, 27 (1): 91-105. 10.1016/j.molcel.2007.06.017.PubMed CentralView ArticlePubMedGoogle Scholar
- Cho WC: OncomiRs: the discovery and progress of microRNAs in cancers. Mol Cancer. 2007, 25 (6): 60-View ArticleGoogle Scholar
- Enerly E, Steinfeld I, Kleivi K, Leivonen SK, Aure MR, Russnes HG, Rønneberg JA, Johnsen H, Navon R, Rødland E, Makela R, Naume B, Perala M, Kallioniemi O, Kristensen VN, Yakhini Z, Børresen-Dale AL: miRNA-mRNA integrated analysis reveals roles for miRNAs in primary breast tumors. PLoS One. 2011, 6 (2): e16915-10.1371/journal.pone.0016915.PubMed CentralView ArticlePubMedGoogle Scholar
- Shalgi R, Lieber D, Oren M, Pilpel Y: Global and local architecture of the mammalian microRNA-transcription factor regulatory network. PLoS Comput Biol. 2007, 3 (7): e131-10.1371/journal.pcbi.0030131.PubMed CentralView ArticlePubMedGoogle Scholar
- Re A, Corá D, Taverna D, Caselle M: Genome-wide survey of microRNA-transcription factor feed-forward regulatory circuits in human. Mol Biosyst. 2009, 5 (8): 854-67. 10.1039/b900177h.PubMed CentralView ArticlePubMedGoogle Scholar
- Riba A, Bosia C, El Baroudi M, Ollino L, Caselle M: A combination of transcriptional and microRNA regulation improves the stability of the relative concentrations of target genes. PLoS Comput Biol. 2014, 10 (2): e1003490-10.1371/journal.pcbi.1003490.PubMed CentralView ArticlePubMedGoogle Scholar
- Toft DJ, Cryns VL: Minireview: Basal-like breast cancer: from molecular profiles to targeted therapies. Mol Endocrinol. 2011, 25 (2): 199-211. 10.1210/me.2010-0164.PubMed CentralView ArticlePubMedGoogle Scholar
- Sandhu R, Rivenbark AG, Mackler RM, Livasy CA, Coleman WB: Dysregulation of microRNA expression drives aberrant DNA hypermethylation in basal-like breast cancer. Int J Oncol. 2014, 44 (2): 563-72.PubMed CentralPubMedGoogle Scholar
- Tang W, Yu F, Yao H, Cui X, Jiao Y, Lin L, Chen J, Yin D, Song E, Liu Q: miR-27a regulates endothelial differentiation of breast cancer stem like cells. Oncogene. 2013, 10: JunGoogle Scholar
- Mertens-Talcott SU1, Chintharlapalli S, Li X, Safe S: The oncogenic microRNA-27a targets genes that regulate specificity protein transcription factors and the G2-M checkpoint in MDA-MB-231 breast cancer cells. Cancer Res. 2007, 15;67 (22): 11001-11.View ArticleGoogle Scholar
- Krishnan K, Steptoe AL, Martin HC, Pattabiraman DR, Nones K, Waddell N, Mariasegaram M, Simpson PT, Lakhani SR, Vlassov A, Grimmond SM, Cloonan N: miR-139-5p is a regulator of metastatic pathways in breast cancer. RNA. 2013, 19 (12): 1767-80. 10.1261/rna.042143.113.PubMed CentralView ArticlePubMedGoogle Scholar
- Lal A, Navarro F, Maher CA, Maliszewski LE, Yan N, O'Day E, Chowdhury D, Dykxhoorn DM, Tsai P, Hofmann O, Becker KG, Gorospe M, Hide W, Lieberman J: miR-24 inhibits cell proliferation by targeting E2F2, MYC, and other cell-cycle genes via binding to "seedless" 3'UTR microRNA recognition elements. Mol Cell. 2009, 35 (5): 610-625. 10.1016/j.molcel.2009.08.020.PubMed CentralView ArticlePubMedGoogle Scholar
- Xu N, Papagiannakopoulos T, Pan G, Thomson JA, Kosik KS: MicroRNA-145 regulates OCT4, SOX2, and KLF4 and represses pluripotency in human embryonic stem cells. Cell. 2009, 137 (4): 647-658. 10.1016/j.cell.2009.02.038.View ArticlePubMedGoogle Scholar
- Li Z, Hassan MQ, Jafferji M, Aqeilan RI, Garzon R, Croce CM, van Wijnen AJ, Stein JL, Stein GS, Lian JB: Biological functions of miR-29b contribute to positive regulation of osteoblast differentiation. J Biol Chem. 2009, 284 (23): 15676-15684. 10.1074/jbc.M809787200.PubMed CentralView ArticlePubMedGoogle Scholar
- Harris MA, Clark J, Ireland A, Lomax J, Ashburner M, Foulger R, Eilbeck K, Lewis S, Marshall B, Mungall C, Richter J, Rubin GM, Blake JA, Bult C, Dolan M, Drabkin H, Eppig JT, Hill DP, Ni L, Ringwald M, Balakrishnan R, Cherry JM, Christie KR, Costanzo MC, Dwight SS, Engel S, Fisk DG, Hirschman JE, Hong EL, Nash RS, Sethuraman A, Theesfeld CL, Botstein D, Dolinski K, Feierbach B, Berardini T, Mundodi S, Rhee SY, Apweiler R, Barrell D, Camon E, Dimmer E, Lee V, Chisholm R, Gaudet P, Kibbe W, Kishore R, Schwarz EM, Sternberg P, Gwinn M, Hannick L, Wortman J, Berriman M, Wood V, de la Cruz N, Tonellato P, Jaiswal P, Seigfried T, White R, Consortium GO: The Gene Ontology (GO) database and informatics resource. Nucleic Acids Res. 2004, 32 (Database): D258-D261.PubMedGoogle Scholar
- Kanehisa M, Goto S, Sato Y, Kawashima M, Furumichi M, Tanabe M: Data, information, knowledge and principle: back to metabolism in KEGG. Nucleic Acids Res. 2014, 42 (Database): D199-205.PubMed CentralView ArticlePubMedGoogle Scholar
- Nishimura D: BioCarta. Biotech Software Internet Report. 2001, 2 (3): 117-120. 10.1089/152791601750294344.View ArticleGoogle Scholar
- Croft D, Mundo AF, Haw R, Milacic M, Weiser J, Wu G, Caudy M, Garapati P, Gillespie M, Kamdar MR, Jassal B, Jupe S, Matthews L, May B, Palatnik S, Rothfels K, Shamovsky V, Song H, Williams M, Birney E, Hermjakob H, Stein L, D'Eustachio P: The Reactome pathway knowledgebase. Nucleic Acids Res. 2014, 42 (Database): 472-7.View ArticleGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Statist Soc B. 1995, 57 (1): 289-300.Google Scholar
- Dai M, Wang P, Boyd AD, Kostov G, Athey B, Jones EG, Bunney WE, Myers RM, Speed TP, Akil H, Watson SJ, Meng F: Evolving gene/transcript definitions significantly alter the interpretation of GeneChip data. Nucleic Acids Res. 2005, 33 (20): e175-10.1093/nar/gni179.PubMed CentralView ArticlePubMedGoogle Scholar
- Wu ZI RA, Gentleman R, Martinez-Murillo F, Spencer F: Model-based background adjustment for oligonucleotide expression arrays. Journal of the American Statistical Association. 2004, 99: 909-917. 10.1198/016214504000000683.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.